Name: Towards AI Legal Name: Towards AI, Inc. Description: Towards AI is the world's leading artificial intelligence (AI) and technology publication. Read by thought-leaders and decision-makers around the world. Phone Number: +1-650-246-9381 Email: [email protected]
228 Park Avenue South New York, NY 10003 United States
Website: Publisher: https://towardsai.net/#publisher Diversity Policy: https://towardsai.net/about Ethics Policy: https://towardsai.net/about Masthead: https://towardsai.net/about
Name: Towards AI Legal Name: Towards AI, Inc. Description: Towards AI is the world's leading artificial intelligence (AI) and technology publication. Founders: Roberto Iriondo, , Job Title: Co-founder and Advisor Works for: Towards AI, Inc. Follow Roberto: X, LinkedIn, GitHub, Google Scholar, Towards AI Profile, Medium, ML@CMU, FreeCodeCamp, Crunchbase, Bloomberg, Roberto Iriondo, Generative AI Lab, Generative AI Lab Denis Piffaretti, Job Title: Co-founder Works for: Towards AI, Inc. Louie Peters, Job Title: Co-founder Works for: Towards AI, Inc. Louis-François Bouchard, Job Title: Co-founder Works for: Towards AI, Inc. Cover:
Towards AI Cover
Logo:
Towards AI Logo
Areas Served: Worldwide Alternate Name: Towards AI, Inc. Alternate Name: Towards AI Co. Alternate Name: towards ai Alternate Name: towardsai Alternate Name: towards.ai Alternate Name: tai Alternate Name: toward ai Alternate Name: toward.ai Alternate Name: Towards AI, Inc. Alternate Name: towardsai.net Alternate Name: pub.towardsai.net
5 stars – based on 497 reviews

Frequently Used, Contextual References

TODO: Remember to copy unique IDs whenever it needs used. i.e., URL: 304b2e42315e

Resources

Unlock the full potential of AI with Building LLMs for Productionβ€”our 470+ page guide to mastering LLMs with practical projects and expert insights!

Publication

Introduction To Genomic Data Science
Latest

Introduction To Genomic Data Science

Last Updated on November 12, 2021 by Editorial Team

Author(s): Γ–mer Γ–zgΓΌr

Introduction To Genomic DataΒ Science

From https://www.bezelyedergi.net/post/biyoinformatik-sekanslama-teknolojileri

Every living organism has its genome at its heart. If the cell is a computer, the genome sequence is the software it executes. If we think of DNA as software run by the cell, we can use our computers to analyze it with similar logic. DNA is not just a storehouse of information.

It is a physical structure that can behave in a complex way. Genomes are incredibly complex machines made up of thousands of parts. Although we know how some genes work today, it is not understood how many genes work together.

  • Genetics sees DNA only as information, looking for patterns in data, investigating the relationships between genes and physical appearance.
  • On the other hand, Genomics sees the genome as a machine and tries to understand how its parts work together.

Genomic DataΒ Science

Genomic Data Science applies methods such as statistics and machine learning found in data science to genomic problems.

Thanks to next-generation sequencing, we can sequence genomes faster, cheaper, and more successfully. For example, the Human Genome Project was carried out for $2.7 billion, and today you can sequence your genome for $1000. For this reason, petabytes of genomic data can be found on the internet.

A living organism is more complex than any machine that humans have ever produced. Solving this complexity is beyond classical algorithms or human understanding. Most genetic diseases and cancers are caused by the interaction of more than one gene and contain genetic variation.

One of the most potent weapons for analyzing incredibly complex and data point data is Deep Learning, as classical methods assume linear relationships between genes, which is not theΒ case.

With Genomic Data Science, personalized drug studies can be accelerated, or if you’re bored, you can create a 3D version of their faces from their genomes to detect people who have thrown cigarettes or gum on theΒ streets.

Life andΒ DNA

From bacteria to whales, life works on similar principles. All living things are linked together in the tree of life and most likely descended from a single common ancestor.

DNA is a long polymer consisting of 4 basics (A, T, G, C) bases. Almost all information about how to build the organism is stored here.
This information itself and how it is processed (epigenetics) changes over timeΒ evolves.

DNA can carry information in near-infinite combinations. The code for people who were never born or for new creatures that never came into existence is waiting to be discovered in thisΒ space.

Central Dogma

If DNA is software, proteins are the most critical hardware. Proteins are tiny machines that do most of the work inside the cell. A molecule called mRNA is needed when converting the information in DNA into proteins. The mRNA goes to the Ribosome, and, according to the information it contains, amino acids are combined, and protein is synthesized. Ribosomes are organic 3D printers.

The Breakdown OfΒ Dogma

Now let’s examine how genomesΒ work.

  • Eukaryotic DNA is wrapped around proteins called histones to fit inside the cell, tightly packed parts are unreadable, and methylated regions are difficult to read. The mechanisms that regulate when DNA should be opened are not fully understood.
  • There is no unidirectional flow of information from DNA to protein. Proteins can also act as regulators by binding toΒ DNA.
  • mRNA carries information about which protein is being synthesized, but it does not know when it should be synthesized. Transcription factors come into play here. Transcription factors bind to specific points in DNA and regulate the expression of the genes locatedΒ nearby.
  • miRNA, siRNA, Riboswitches can also participate in the editingΒ task.

Note: Not all editing mechanisms in the genome are understood, and the process is more complex than described here.

Predicting Transcription FactorΒ Binding

We found that a cell is complex and challenging to work with classical methods so that we will practice deep learning techniques. We will use the HepG2 cell line and JunD Transcription factor asΒ data.

HepG2 is an immortal cell line derived from the liver tissue of a 15-year-old African-American child. JunD Transcription factor encoded by the JUND gene. It can downregulate or activate otherΒ genes.

We chose the 22nd Chromosome so that the genomic data could be processed. This Chromosome contains about 50 million base pairs. There are about 3 billion base pairs in humans. The largest Chromosome is the 1st, and the smallest is theΒ 22nd.

Genomic data is stored in FASTA or FASTQ format. Data will look like bookΒ rows.

>chr22
NNNNNNNNNNNNNNNNNNNNNNNNNNN
TCCCAAATTGTGGAAGGAATGTACATTTGAC

Here we see the base sequences contained in a single strand of DNA. We can ask what the β€œN” symbol is doing in DNA. This means that it cannot be decided which base is read in the sequencing process.

Next, we need to convert the genomic information into a mathematical format that we can use. As we can see, our data consists of parts A, T, G, C, and N. We can express these parts with One-hot encoding. For example A:[1,0,0,0] or G:[0, 0,1,0] and N:[0.25,0.25,0.25,0.25] can be represented.

[[1., 0., 0., 0.],
[1., 0., 0., 0.],
[0., 1., 0., 0.]
[0., 0., 1., 0.],
[0., 0., 0., 1.]] == ["AATGC"]

After DNA one-hot is encoded, we can start training. Our training data consists of 101 bases and whether the transcription factor is bound to this sequence. For example, if there are 2000 bases, we will have 20 rows. After this point, the problem turns into supervised learning.

Our x[0] == [[1., 0., 0., 0.],
[1., 0., 0., 0.],
[0., 1., 0., 0.]
[0., 0., 1., 0.],
[0., 0., 0., 1.]]..
our y[0] == 1 

We can use many models to find the relationship between the x’s and y’s. For this example we will useΒ CNN.

From https://www.researchgate.net/figure/Causal-convolution-operation-in-a-1D-convolutional-layer-with-k-3-kernel-size-Input_fig1_337703712

Since DNA is a string, it is one-dimensional. We can use different architectures. One of the most useful architectures for finding patterns in 1D data is Convolutional Neural Networks. They are generally used on images (2D information) and 1D CNNs can be used to processΒ texts.

2D CNNs extract features from the image, learning about edges, vertices, color changes, and more complex patterns. Similarly, we can learn filters that can extract features from text data. Once our model is trained, we can calculate the probability that the JunD protein will bind to DNA fragments 101 basesΒ long.

Generating More Realistic Results

There are many factors in JunD protein binding to DNA: accessibility, methylation, shape, presence of other molecules.

From https://www.scinexx.de/dossierartikel/auch-die-verpackung-machts/

Chromatin accessibility: Defines how accessible the DNA is to molecules from outside. When DNA is tightly wrapped around histones, Transcription Factors or other molecules are inaccessible. Similarly, genes to be silenced can be packaged.

The accessibility of a region is not constant. It changes over time. If we add the chromatin accessibility information to our dataset, we can make more realistic predictions.

Conclusion and ExampleΒ Code

The complexity of biology is beyond human understanding. This is why we need machine learning to understand and improve ourselves.

When the genome was sequenced, scientists thought they had solved all the secrets, but they had only scratched the surface. Genome editing tools like CRISPR don’t solve our problem either. We cannot edit without predicting the consequences of genetic changes. Unraveling the secret of the genome, if not in an inexplicable way, would provide an opportunity to design new living things to eradicate many diseases.

Data and sample code can be found in my GitHubΒ account.

GitHub – OmerOzgur271/Genomic-Datascience


Data Science was originally published in Towards AI on Medium, where people are continuing the conversation by highlighting and responding to this story.

Published via Towards AI

Feedback ↓