Data Science to Analyze Big Genomic Data

Posted on Jul 28, 2018

Finding the neural stem cell populations in mouse brain



The main objective of this project is to build an interactive app to process and analyze genomic data with the objective of identifying new stem cell populations in mouse brain. Characterizing the unique gene expression signature of these cells could be the starting point for finding and defining cancer stem cells (CSC) in human tumors, particularly glioblastoma,   the most common malignant brain tumor diagnosed in adults. Despite numerous advances in cancer therapy, these tumors remain incurable, suggesting that current treatments fail to target the cells responsible for tumor growth. One possible explanation for this would be that CSCs are resistant to current therapies and responsible for tumor remission.

In order to develop new therapies that target CSC, reliable marker genes (a marker gene is a gene that is uniquely expressed in a given cell type) that define the “CSC signature” (a molecular signature is a group of marker genes of a given cell type), are needed. The first step is to study mouse brain and define their Neural Stem Cell (NSC) signature. Based on the mouse NSC Signature, we would define the Cancer Stem Cell signature in human tumors. And finally we would use the Cancer Stem Cell marker genes to screen for new drugs that could lead to novel Glioblastoma treatments.

The dataset that I used for this project is the result of dissecting a specific part of a mouse’s brain (dentate gyrus, responsible for learning and memory), dissociating the cells into a single cell suspension and performing Next Generation Sequencing on individual cells. In particular, single cell RNA sequencing. This technique allows the quantification of all the RNA transcript of each individual cell. We quantify RNA instead of DNA because whereas the DNA content is the same in every cell of our body, the RNA expression profile uniquely defines a cell’s identity and function. It is what distinguish a skin cell from a neuron.

After filtering low quality cells, the final data object had 18572 rows and 17922 columns. Each row corresponds to a gene, and each column corresponds to a cell. The values in the database are a measure of the RNA expression of each gene in each cell. Each cell expresses an average of 2048 genes. The size of next generation sequencing genomic data -- in this case the processed database --consumed 3.2Gb of memory, makes Data science and machine learning algorithms indispensable tools for handling this kind of  big data.



Who is who

Next, we are going to determine the identity of each cluster based on known marker genes. As this is unpublished data, the real names of the clusters are not shown. After filtering the data, I performed a clustering analysis based on K-nearest neighbor KNN algorithm. The resulting plots are called t-SNE (T-distributed Stochastic Neighbor Embedding) plots. In these plots, each cell is represented by a dot (hence each dot summarizes the information of all the rows and columns referred to the expression levels for one single cell). The spatial proximity between the dots is related to the molecular similarities. Therefore, cells that cluster together belong to the same cell type.

In order to define “who is who” and where  our stem cell populations are, we have to interrogate the gene expression of the different clusters. For example the genes Gfap and Id4 are known to be expressed in astrocytes and stem cells; therefore, we know that the highlighted groups are astrocytes (here “cell type b”) and stem cells (here “stem cell 1” and “stem cell 2”). Other genes   such as Rbfox3 and Syt1, are known to be expressed in neurons. Consequently, by selecting those genes in our app, we can assign the highlighted clusters to neuronal populations.

Going through different gene expression profiles, with the help of this interactive tSNE plot and dot plot, we figure out which cluster belongs to which cell type and, most importantly, where our stem cell clusters are. In particular, I was able to find a new SC cluster, here called “stem cell 1” that is related to the previously identified, here called “stem cell 2”.



Developmental time

After clustering, I used machine learning algorithms to further determine the relative location of each cluster in the neurogenic differentiation trajectory. I used a semi-supervised ordering mode to order the cells along the trajectory based on known marker genes that define the biological progress through the neuronal differentiation process and other (1000) genes that covary with them. Once the algorithm has learned the specific trajectory of expression changes, it places each cell at its proper position in the trajectory all the way from the most primitive stem cell state to a more differentiated neuroblast, going through the intermediate NPC (Neural Progenitor Cell, here called “cell type b”) state.

The left panel shows the stem cell to neuron differentiation trajectory, and the right panel shows the expression of specific genes along developmental Pseudotime. As seen in the plot below, stem cell genes like Gfap and Id4 start high and end low.  When selecting neurogenic genes in the app, such as Dcx and Calb2, we can confirm that they start low and end up high as expected.

In addition to finding the “developmental” order of our clusters, I could see that the newly discovered stem cell cluster (“stem cell 1”) starts at the beginning of the trajectory, supporting that it indeed corresponds to a stem cell cell population.




So far, we have identified the different cell types, including the stem cell populations, we have reconstructed developmental time and finally it is time to find new genes, new NSC markers.

In order to find new genes expressed in our clusters, we use heatmaps. The red color indicates highly expressed genes, and the blue color indicates low expressed genes. The graph shows the top 15 genes (highly expressed) on the right hand side of the graph per cluster. In the heatmap we look for genes that are highly expressed in the stem cell populations (“stem cell 1” and “stem cell 2”) but not in the other cell types. Genes that are uniquely expressed in the stem cell clusters are ideal candidates for new stem cell marker genes. Ultimately, NSC marker genes could be the first step towards the discovery and validation of CSC marker genes in human tumors. Lastly, CSC marker genes would be used to screen for new drugs that could lead to the development of new therapies for glioblastoma.




  • All the cell types in this particular region of mouse brain were identified based on expression of known marker genes
  • A new neural stem cell population (“stem cell 1”) was identified and characterized
  • The newly defined stem cell group is located at the beginning of the differentiation trajectory
  • Uniquely expressed genes were identified in a new stem cell population that are ideal candidates for new stem cell markers

See my GitHub code

Link to the app


Disclaimers and Acknowledgements

This blog includes unpublished data and is not intended to be used for submissions or citations.

The data presented in this blog was generated in Luis Parada’s laboratory at Memorial Sloan Kettering Cancer Center. Single cell RNA sequencing (10X Genomics Chromium) was performed at the Epigenomics Core at Weill Cornell Medicine.

Seurat package for R was used for the pre-processing and clustering and Monocle package for R was used for pseudotime analysis.

Featured image is a courtesy of National Human Genome Research Institute (NHGRI).

About Author

Elsa Vera

Elsa Vera has a PhD in Biology, Biochemistry and Biomedicine. She is currently working as a Data Scientist at JPMorgan Chase. Before joining the bootcamp she was a scientist in the field of genomics, cancer stem cells and...
View all posts by Elsa Vera >

Related Articles

Leave a Comment

No comments found.

View Posts by Categories

Our Recent Popular Posts

View Posts by Tags

#python #trainwithnycdsa 2019 2020 Revenue 3-points agriculture air quality airbnb airline alcohol Alex Baransky algorithm alumni Alumni Interview Alumni Reviews Alumni Spotlight alumni story Alumnus ames dataset ames housing dataset apartment rent API Application artist aws bank loans beautiful soup Best Bootcamp Best Data Science 2019 Best Data Science Bootcamp Best Data Science Bootcamp 2020 Best Ranked Big Data Book Launch Book-Signing bootcamp Bootcamp Alumni Bootcamp Prep boston safety Bundles cake recipe California Cancer Research capstone car price Career Career Day citibike classic cars classpass clustering Coding Course Demo Course Report covid 19 credit credit card crime frequency crops D3.js data data analysis Data Analyst data analytics data for tripadvisor reviews data science Data Science Academy Data Science Bootcamp Data science jobs Data Science Reviews Data Scientist Data Scientist Jobs data visualization database Deep Learning Demo Day Discount disney dplyr drug data e-commerce economy employee employee burnout employer networking environment feature engineering Finance Financial Data Science fitness studio Flask flight delay gbm Get Hired ggplot2 googleVis H20 Hadoop hallmark holiday movie happiness healthcare frauds higgs boson Hiring hiring partner events Hiring Partners hotels housing housing data housing predictions housing price hy-vee Income Industry Experts Injuries Instructor Blog Instructor Interview insurance italki Job Job Placement Jobs Jon Krohn JP Morgan Chase Kaggle Kickstarter las vegas airport lasso regression Lead Data Scienctist Lead Data Scientist leaflet league linear regression Logistic Regression machine learning Maps market matplotlib Medical Research Meet the team meetup methal health miami beach movie music Napoli NBA netflix Networking neural network Neural networks New Courses NHL nlp NYC NYC Data Science nyc data science academy NYC Open Data nyc property NYCDSA NYCDSA Alumni Online Online Bootcamp Online Training Open Data painter pandas Part-time performance phoenix pollutants Portfolio Development precision measurement prediction Prework Programming public safety PwC python Python Data Analysis python machine learning python scrapy python web scraping python webscraping Python Workshop R R Data Analysis R language R Programming R Shiny r studio R Visualization R Workshop R-bloggers random forest Ranking recommendation recommendation system regression Remote remote data science bootcamp Scrapy scrapy visualization seaborn seafood type Selenium sentiment analysis sentiment classification Shiny Shiny Dashboard Spark Special Special Summer Sports statistics streaming Student Interview Student Showcase SVM Switchup Tableau teachers team team performance TensorFlow Testimonial tf-idf Top Data Science Bootcamp Top manufacturing companies Transfers tweets twitter videos visualization wallstreet wallstreetbets web scraping Weekend Course What to expect whiskey whiskeyadvocate wildfire word cloud word2vec XGBoost yelp youtube trending ZORI