Package: dnadiffusion
DNA DiffusionGenerative modeling of regulatory DNA sequences with diffusion probabilistic models.Documentation:https://pinellolab.github.io/DNA-DiffusionSource Code:https://github.com/pinellolab/DNA-DiffusionAbstractThe Human Genome Project has laid bare the DNA sequence of the entire human genome, revealing the blueprint for tens of thousands of genes involved in a plethora of biological process and pathways.
In addition to this (coding) part of the human genome, DNA contains millions of non-coding elements involved in the regulation of said genes.Such regulatory elements control the expression levels of genes, in a way that is, at least in part, encoded in their primary genomic sequence.
Many human diseases and disorders are the result of genes being misregulated.
As such, being able to control the behavior of such elements, and thus their effect on gene expression, offers the tantalizing opportunity of correcting disease-related misregulation.Although such cellular programming should in principle be possible through changing the sequence of regulatory elements, the rules for doing so are largely unknown.
A number of experimental efforts have been guided by preconceived notions and assumptions about what constitutes a regulatory element, essentialy resulting in a "trial and error" approach.Here, we instead propose to use a large-scale data-driven approach to learn and apply the rules underlying regulatory element sequences, applying the latest generative modelling techniques.Introduction and Prior WorkThe goal of this project is to investigate the application and adaptation of recent diffusion models (seehttps://lilianweng.github.io/posts/2021-07-11-diffusion-models/for a nice intro and references) to genomics data. Diffusion models are powerful models that have been used for image generation (e.g. stable diffusion, DALL-E), music generation (recent version of the magenta project) with outstanding results.
A particular model formulation called "guided" diffusion allows to bias the generative process toward a particular direction if during training a text or continuous/discrete labels are provided. This allows the creation of "AI artists" that, based on a text prompt, can create beautiful and complex images (a lot of examples here:https://www.reddit.com/r/StableDiffusion/).Some groups have reported the possibility of generating synthetic DNA regulatory elements in a context-dependent system, for example, cell-specific enhancers.
(https://elifesciences.org/articles/41279,https://www.biorxiv.org/content/10.1101/2022.07.26.501466v1)Step 1: generative modelWe propose to develop models that can generate cell type specific or context specific DNA-sequences with certain regulatory properties based on an input text prompt.
For example:"A sequence that will correspond to open (or closed) chromatin in cell type X""A sequence that will activate a gene to its maximum expression level in cell type X""A sequence active in cell type X that contains binding site(s) for the transcription factor Y""A sequence that activates a gene in liver and heart, but not in brain"Step 2: extensions and improvementsBeyond individual regulatory elements, so called "Locus Control Regions" are known to harbour multiple regulatory elements in specific configurations, working in concert to result in more complex regulatory rulesets. Having parallels with "collaging" approaches, in which multiple stable diffusion steps are combined into one final (graphical) output, we want to apply this notion to DNA sequences with the goal of designing larger regulatory loci. This is a particularly exciting and, to our knowledge, hitherto unexplored direction.Besides synthetic DNA creations, a diffusion model can help understand and interpret regulatory sequence element components and for instance be a valuable tool for studying single nucleotide variations (https://www.biorxiv.org/content/10.1101/2022.08.22.504706v1) and evolution.
(https://genomebiology.biomedcentral.com/articles/10.1186/s13059-018-1502-5)Taken together, we believe our work can accelerate our understanding of the intrinsic properties of DNA-regulatory sequence in normal development and different diseases.Proposed frameworkFor this work we propose to build a Bit Diffusion model based on the formulation proposed by Chen, Zhang and Hintonhttps://arxiv.org/abs/2208.04202. This model is a generic approach for generating discrete data with continuous diffusion models. An implementation of this approach already exists, and this is a potential code base to build upon:https://github.com/lucidrains/bit-diffusionTasks and potential roadmap:Collecting genomic datasetsImplementing the guided diffusion based on the code baseThinking about the best encoding of biological information for the guided diffusion (e.g. cell type: K562, very strong activating sequence for chromatin, or cell type: GM12878, very open chromatin)Plans for validation based on existing datasets or how to perform new biological experiments (we need to think about potential active learning strategies).DeliverablesDataset:compile and provide a complete database of cell-specific regulatory regions (DNAse assay) to allow scientists to train and generate different diffusion models based on the regulatory sequences.Models:Provide a model that can generate regulatory sequences given a specific cell type and genomic context.API:Provide an API to make it possible to manipulate DNA regulatory models and a visual playground to generate synthetic contextual sequences.DatasetsDHS Index:Chromatin (DNA + associated proteins) that is actively used for the regulation of genes (i.e. "regulatory elements") is typically accessible to DNA-binding proteins such as transcription factors (review,relevant paper).
Through the use of a technique calledDNase-seq, we've measured which parts of the genome are accessible across 733 human biosamples encompassing 438 cell and tissue types and states, resulting in more than 3.5 million DNase Hypersensitive Sites (DHSs).
Using Non-Negative Matrix Factorization, we've summarized these data into 16components, each corresponding to a different cellular context (e.g. 'cardiac', 'neural', 'lymphoid').For the efforts described in this proposal, and as part of an earlierongoing projectin the research group of Wouter Meuleman,
we've put together smaller subsets of these data that can be used to train models to generate synthetic sequences for each NMF component.Please find these data, along with a data dictionary,here.Other potential datasets:DNA-sequences data corresponding to annotated regulatory sequences such as gene promoters or distal regulatory sequences such as enhancers annotated (based on chromatin marks or accessibility) for hundreds of cells by the NHGRI funded projects like ENCODE or Roadmap Epigenomics.Data from MPRA assays that test the regulatory potential of hundred of DNA sequences in parallel (https://elifesciences.org/articles/69479.pdf,https://www.nature.com/articles/s41588-021-01009-4, ... )MIAA assays that test the ability of open chromatin within a given cell type.ModelsInput modality:A) Cell type + regulatory element ex: Brain tumor cell weak Enhancer
B) Cell type + regulatory elements + TF combination (presence or absence) Ex: Prostate cell, enhancer , AR(present), TAFP2a (present) and ER (absent),
C) Cell type + TF combination + TF positions Ex: Blood Stem cell GATA2(presence) and ER(absent) + GATA1 (100-108)
D) Sequencing having a GENETIC VARIANT -> low number diffusion steps = nucleotide importance predictionOutput:DNA-sequenceModel size:The number of enhancers and biological sequences isn’t bigger than the number of available images on the Lion dataset. The dimensionality of our generated DNA outputs should not be longer than 4 bases [A,C,T,G] X ~1kb. The final models should be bigger than ~2 GB.Models:Different models can be created based on the total sequence length.APIsTBD depending on interestPaperCan the project be turned into a paper? What does the evaluation process for such a paper look like? What conferences are we targeting? Can we release a blog post as well as the paper?Yes, We intend to have a mix of our in silico generations and experimental validations to study our models' performance on classic regulatory systems ( ex: Sickle cell and Cancer).
Our group and collaborators present a substantial reputation in the academic community and different publications in high-impact journals, such as Nature and Cell.Resources RequirementsWhat kinds of resources (e.g. GPU hours, RAM, storage) are needed to complete the project?Our initial model can be trained with small datasets (~1k sequences) in about 3 hours ( ~500 epochs) on a colab PRO (24GB ram ) single GPU Tesla K80. Based on this we expect that to train this or similar models on the large dataset mentioned above ( ~3 million sequences (4x200) we will need several high-performant GPUs for about 3 months. ( Optimization suggestions are welcome!)TimelineWhat is a (rough) timeline for this project?6 months to 1 year.Broader ImpactHow is the project expected to positively impact biological research at large?We believe this project will help to better understand genomic regulatory sequences: their composition and the potential regulators acting on them in different biological contexts and with the potential to create therapeutics based on this knowledge.ReproducibilityWe will use best practices to make sure our code is reproducible and with versioning. We will release data processing scripts and conda environments/docker to make sure other researchers can easily run it.We have several assays and technologies to test the synthetic sequences generated by these models at scale based on CRISPR genome editing or massively parallel reporter assays (MPRA).Failure CaseRegardless of the performance of the final models, we believe it is important to test diffusion models on novel domains and other groups can build on top of our investigations.Preliminary FindingsUsing the Bit Diffusion model we were able to reconstruct 200 bp sequences that presented very similar motif composition to those trained sequences. The plan is to add the cell conditional variables to the model to check how different regulatory regions depend on the cell-specific context.Next StepsExpand the model length to generate complete regulatory regions (enhancers + Gene promoter pairs)
Use our synthetic enhancers on in-vivo models and check how they can regulate the transcriptional dynamics in biological scenarios (Besides the MPRA arrays).How to contributeIf this project sounds exciting to you,please join us!
Join the OpenBioML discord:https://discord.gg/Y9CN2dUzQJ, we are discussing this project in thedna-diffusionchannel and we will provide instructions on how to get involved.Known contributorsYou can access the contributor listhere.DevelopmentSetup environmentWe usehatchto manage the development environment and production build. It is often convenient to install hatch withpipx.Run unit testsYou can run all the tests with:hatchruntestFormat the codeExecute the following command to apply linting and check typing:hatchrunlintPublish a new versionYou can check the current version with:hatchversionYou can bump the version with commands such ashatch version devorpatch,minorormajor. Or edit thesrc/dnadiffusion/__about__.pyfile. After changing the version, when you push to github, the Test Release workflow will automatically publish it on Test-PyPI and a github release will be created as a draft.Serve the documentationYou can serve the mkdocs documentation with:hatchrundocs-serveThis will automatically watch for changes in your code.