tag:blogger.com,1999:blog-41710865483238458272024-02-19T09:00:19.462-08:00Plants, Systems & BioinformaticsAnonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.comBlogger8125tag:blogger.com,1999:blog-4171086548323845827.post-80383552869202313042016-07-09T15:53:00.001-07:002016-07-09T15:57:43.454-07:00Scientific Computing Resources For the Biologist<span style="font-family: "times" , "times new roman" , serif;">In scientific computing there are typically two types of computing referenced:</span><br />
<ol>
<li><span style="font-family: "times" , "times new roman" , serif;">High performance computing (HPC)</span></li>
<li><span style="font-family: "times" , "times new roman" , serif;">High throughput computing (HTC)</span></li>
</ol>
<span style="font-family: "times" , "times new roman" , serif;">Both are meant for working with "big data". Which of these is more important for your research depends on the how parallel the code being executed on these large datasets is. Additionally there are also cases where a hybrid is between HPC and HTC is best. The problem you will face is typically when to use HPC, HTC or a combination. Knowing which computing cluster is best for your research will allow you to choose resources to apply for.</span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<br />
<h3>
<span style="font-family: "times" , "times new roman" , serif; font-size: x-large;"><span style="font-size: x-large;">High Performance</span><span style="font-size: x-large;"> Computing (HPC)</span></span></h3>
<div>
<span style="font-family: "times" , "times new roman" , serif; font-size: x-large;"><span style="font-size: x-small;"><br /></span></span></div>
<span style="font-family: "times" , "times new roman" , serif;">High performance computers (Supercomputers) can run code at or near its optimal rate. In other words HPC is capable of running code that normally takes days to running it in matter of seconds.</span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;">An example would be using a 15 base pair sliding window to call differentially methylated regions between two samples. Because this dataset contains information (methylated or not methylated) for every cytosine, each sample datasets is around 100 MB (in this example). Also the sliding window approach does not allow us to easily split up the data and run parallel tasks (the most we can split is by Chromosome). Also I (A Biologist by training) wrote this code and know very little about optimizing my algorithms. So I can wait for my slow code to run on a small server (where I may run into memory swapping issues) or use brute force to run this code on a HPC that can handle the work load. In summary, HPC servers are best used when you have a large (memory and time) task to run that can not be optimized further.</span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;">HPC is facilitated through access to large computing power containing thousands of processors in close proximity working together (cluster) with a large shared memory capacity and large working memory distributed across many processors. Because of the large computational capacity (memory and processors) needed for HPC they are expensive and usually shared with many different researcher groups. Some examples of these may include your university's cluster or resources like <a href="https://www.xsede.org/high-performance-computing">XSEDE</a>. XSEDE is a powerful resource because it gives researchers access to a variety of clusters meant to individual computing needs. One of the advantages and disadvantage to sharing resources and using a cluster is that you are not granted full control of the system. You therefore do not have sudo rights and must locally install everything. Many clusters have pre-installed programs for commonly used programs like BLAST but they may not be the versions you are looking for. XSEDE's pre-installed software can be found <a href="https://portal.xsede.org/software#/">here</a>.</span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;">Amazon and Microsoft are dramatically cutting the cost of using large computing resources that are not shared and can easily be personalized (The advantages of these are explained more here).</span><br />
<span style="font-family: "times" , "times new roman" , serif;"><span style="font-size: x-large;"><br /></span>
</span><br />
<h3>
<span style="font-family: "times" , "times new roman" , serif;"><span style="font-size: x-large;">High Throughput</span><span style="font-size: x-large;"> Computing (HTC)</span></span></h3>
<span style="font-family: "times" , "times new roman" , serif;"><span style="font-size: x-large;"><br /></span>
High throughput computing also takes advantage of many processors but rather than all processors working together on one task they work separately on many tasks. This type of computing is best for code that that runs the same task repeatedly (can run many tasks in parallel).</span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;">An example of this is code that needs to run 1 million blast queries against the same reference genome. Each individual blast command is small and can be ran with little memory. But running all one million blasts is a huge task that would take days to run on one processor. HTC allows each task to be ran on a separate processor, breaking a large file of 1 million queries into 1 million separate fasta files that are ran simultaneously. Running each task separately not only cuts down the memory needed on each processor but also dramatically speeds the process.</span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<br />
<div>
<span style="font-family: "times" , "times new roman" , serif;">The open science grid (<a href="https://www.opensciencegrid.org/">OSG</a>) is an example of a HTC resource. OSG is a collection of computers connected throughout the US in a grid. This network is extremely large and therefore powerful (put number of computers here). The grid is composed of a variety of machines donating time to the grid. These machines range from personal computers to university clusters. Because the OSG has to deal with a large variety of software and machines they developed a program (HTCondor) for matching idle computers with tasks best suited for them. HTCondor carries out this matching process using user defined system requirements to find computers on the grid that are available, have the system requirements requested (memory, unix) and have requested software installed (blast or matlab). DAGman also used by OSG allows for optimization in submitting jobs to determining the series of functions carried out by which node and when. </span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<br />
<h4>
<span style="font-family: "times" , "times new roman" , serif; font-size: x-large;">Using both/either HPC and HTC </span></h4>
<div>
<span style="font-family: "times" , "times new roman" , serif; font-size: x-large;"><br /></span></div>
<span style="font-family: "times" , "times new roman" , serif;">Amazon and Microsoft also provide computational resources to researchers at a low cost (or through educational grants). These resources were developed for the general public and therefore are much more flexible. Here you rent a computer of the desired type and size for a desired period of time. Both XSEDE and OSG require using their HPC or HTC service for the length of the grant (1 year access to machines) where Amazon and Microsoft charge the grant by the hour of machines and number of machines used. This allows you to use a system with both HTC and HPC characteristics depending on the number of machines and size of the machines rented out at any given time. Of course renting out many large machines and/or cost much more then renting out a single large machine or multiple small machines. </span><br />
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;">An advantage to paying hourly for computational time is that resources can easily be man available to users when needed because everyone is held accountable for time spent on them. Both XSEDE and OSG allow users unlimited usage to users. This often leads to users waiting for requested resources to </span><span style="font-family: "times" , "times new roman" , serif;">become available. When resources are in high depend OSG uses a fair use policy preventing users from hogging a resources. </span><br />
<span style="font-family: "times" , "times new roman" , serif;"><span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;">I hope this article help clarified the differences to help you decide when to use one resource over the other or if you need to use both in combination.</span></span><br />
<span style="font-family: "times" , "times new roman" , serif;"><span style="font-family: "times" , "times new roman" , serif;"><br /></span>
<span style="font-family: "times" , "times new roman" , serif;"><br /></span></span></div>
<div>
<h3>
<span style="font-family: "times" , "times new roman" , serif; font-size: x-large;">Access to resources </span></h3>
<span style="font-family: "times" , "times new roman" , serif;"><span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span>
</span><br />
<div style="line-height: normal;">
</div>
<span style="font-family: "times" , "times new roman" , serif;"><span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span>
</span><br />
<ul><span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="font-family: "times" , "times new roman" , serif; line-height: 24px;">
<li><b><a href="https://portal.xsede.org/allocations-overview#writing"><span style="font-size: medium;">XSEDE (</span>Mostly HPC but now has access to the OSG)</a></b></li>
<ol>
<li><span style="font-size: medium;">Who can apply:</span></li>
<ul>
<li><span style="font-size: medium;"><span style="font-size: small;">Posdocs and Professors</span></span></li>
<li>Graduate Students with NSF-GRFP </li>
</ul>
<li><span style="font-size: medium;">How do you apply:</span></li>
<ul>
<li>First year requires an abstract and information on computational needs</li>
<li>After a year a 2-page grant is required for access to XSEDE</li>
</ul>
</ol>
</span></span></span></ul>
<ol style="line-height: 24px;"><span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="font-family: "times" , "times new roman" , serif; line-height: 24px;">
</span></span></span></ol>
<span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;"><span style="font-family: "times" , "times new roman" , serif;">
</span></span></span></span>
<br />
<ol style="line-height: 24px;"><ol><ul><span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span></ul>
<span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span></ol>
<span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span></ol>
<span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span>
<br />
<ul style="line-height: 24px;"><span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
<li><span style="font-family: "times" , "times new roman" , serif; font-size: medium;"><b><a href="https://twiki.opensciencegrid.org/bin/view/Education/OSGUserSchool2016">OSG (HTC)</a></b></span></li>
<ol>
<li><span style="font-family: "times" , "times new roman" , serif; font-size: medium;">Who can apply:</span></li>
<ul>
<li><span style="font-family: "times" , "times new roman" , serif;">Any researcher</span></li>
</ul>
<li><span style="font-family: "times" , "times new roman" , serif; font-size: medium;">How do you apply:</span></li>
<ul>
<li><span style="font-family: "times" , "times new roman" , serif;">Researchers can apply to OSG user school where they receive training plus 1-year free on the grid after that they need to join an organized group on OSG.</span></li>
<li><span style="font-family: "times" , "times new roman" , serif;">This application is a few pages on how the OSG will help facilitate your research needs.</span></li>
</ul>
</ol>
</span></span></span></ul>
<span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span>
<br />
<ul><span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
<li style="line-height: 24px;"><span style="font-family: "times" , "times new roman" , serif; font-size: medium;"><b><a href="https://aws.amazon.com/grants/">Amazon</a> & <a href="https://azure.microsoft.com/en-us/community/education/">Microsoft</a> (HPC, HTC, MTC)</b></span></li>
<ol>
<li style="line-height: 24px;"><span style="font-size: medium;"><span style="font-family: "times" , "times new roman" , serif;">Who can apply:</span></span></li>
<ul>
<li><span style="font-family: "times" , "times new roman" , serif;"><span style="font-size: medium; line-height: 24px;">Students and Postdocs at research universities</span><span style="font-size: medium; line-height: 24px;">. </span></span></li>
</ul>
<li style="line-height: 24px;"><span style="font-size: medium;"><span style="font-family: "times" , "times new roman" , serif;">How do you apply:</span></span></li>
</ol>
<ul style="line-height: 24px;"><ul>
<li><span style="font-family: "times" , "times new roman" , serif;">Submit an abstract and information on computational needs </span></li>
</ul>
</ul>
</span></span></span></ul>
<span style="color: #474747; font-family: "raleway" , "helvetica neue" , "helvetica" , sans-serif;"><span style="background-color: white;"><span style="line-height: 24px;">
</span></span></span></div>
Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com3tag:blogger.com,1999:blog-4171086548323845827.post-19495236627551324322016-03-13T18:09:00.000-07:002016-06-04T17:57:05.128-07:00Using AWS for Research<!-- SyntaxHighlighter -->
<!-- 1- Core javascript and css -->
<link href="http://alexgorbatchev.com/pub/sh/current/styles/shCore.css" rel="stylesheet" type="text/css"></link>
<link href="http://alexgorbatchev.com/pub/sh/current/styles/shCoreDefault.css" rel="stylesheet" type="text/css"></link>
<link href="http://alexgorbatchev.com/pub/sh/current/styles/shThemeDefault.css" rel="stylesheet" type="text/css"></link>
<script src="http://alexgorbatchev.com/pub/sh/current/scripts/shCore.js" type="text/javascript">
<script src='http://alexgorbatchev.com/pub/sh/current/scripts/shAutoloader.js' type='text/javascript'/>
<!-- 2- Brushes -->
<script src='http://alexgorbatchev.com/pub/sh/current/scripts/shBrushPlain.js' type='text/javascript'/>
<script src='http://alexgorbatchev.com/pub/sh/current/scripts/shBrushBash.js' type='text/javascript'/>
<!-- 3- Activate SyntaxHighlighter -->
<script language='javascript' type='text/javascript'>
SyntaxHighlighter.config.bloggerMode = true;
SyntaxHighlighter.all();
</script>
<br />
<h2>
<b><span style="font-size: x-large;">What is Amazon Web Services (AWS)? </span></b></h2>
<br />
AWS is a cloud computing service provided by Amazon. Amazon EC2 is the service we'll go over in this post. It allows users to launch their own virtual instances from a variety of operating systems.
Amazon provides these computational resources to the public and private sector at an extremely low price. I have listed below both the advantages and disadvantages of using cloud computing resources. For me the advantages definitely outweigh the disadvantages :) AWS is generally much more flexible than your universities local cluster in that you can choose the size, type and number of machines and have the rights to download and run whenever needed.<br />
<h3>
<span style="font-size: large;"> Advantages</span> </h3>
<ol>
<li>None of the maintenance of hosting your own server </li>
<li>Access to high computing machines </li>
<li>Access to high parallel computing machines (like those running hadoop) </li>
<li>The cost is extremely cheap </li>
<li>Sudo user rights (can install whatever you want without asking your sysadmin) </li>
<li>The computers available to you come in a variety shapes and sizes</li>
</ol>
<h3>
<span style="font-size: large;"> Disadvantages</span> </h3>
<ol>
<li>There might not be enough machines for everyone and their maybe a wait time (This has yet happen) </li>
<li>It cost money (but amazon provides educational grants) </li>
<li>AWS is a blank slate therefore you need to install everything and copy files over (They make this easy with thing like s3 for storing files and allow you to create images of an instance so you can relaunch one with software pre-installed) </li>
</ol>
<h3>
<span style="font-size: large;">AWS Getting Started </span></h3>
Using EC2 is easy once you've launched your instance you can ssh using your ssh keys:<br />
If your machine is ubuntu then <code>ssh -i location_to_pem_file.pem ubuntu@ec2-54-153-7-122.us-west-1.compute.amazonaws.com</code>
All inputs and outputs should be saved in your /mnt/ directory (this is where all of the system storage is).<br />
<br />
By default you are not the owner of this directory so you need to change the permissions.<br />
<br />
<code>
cd /mnt/<br />
sudo chown ubuntu:ubuntu .<br />
mkdir data<br />
</code>
<br />
<h3>
If you need to mount more storage on your EC2 machine: </h3>
<ol>
<li> Go to AWS Console </li>
<li>Create an EBS volume of the size needed </li>
<li>Attach this volume to your running EC2 machine </li>
<li>Then mount the volume to your EC2 machine using the following commands </li>
</ol>
<br />
<pre class="brush: bash">sudo mkfs -t ext4 /dev/xvdf
sudo mount /dev/xvdf /mnt/data
# /dev/xvdf might be diff depending on where volume was added
cd /mnt/data/
mount
# mount while in data dir
sudo chown ubuntu:ubuntu .
df
# df to confirm that your volume was successfully added
</pre>
<pre class="brush: bash"></pre>
<br />
<h3>
All data can be easily on saved on <a href="https://aws.amazon.com/s3/">S3</a> and instances or machines with preinstalled packages can be saved and reused as <a href="http://docs.aws.amazon.com/AWSEC2/latest/UserGuide/AMIs.html">images</a></h3>
Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com15tag:blogger.com,1999:blog-4171086548323845827.post-11900166689918275092015-01-18T12:22:00.001-08:002015-05-11T12:23:29.585-07:00Change node and edge colors in Cytoscape with PythonOften with Cytoscape we are interested in overlaying other information such as expression data. In order to easily do this without changing the layout/format of the original file it requires a bit of programming. Here I will explain how this can be easily done in python and have provided some of my own code (although it is still a bit messy).<br />
<br />
First export your Cytoscape file to an XML file. This can be done in Cytoscape under file> export> network. XML is a superset of HTML, this file contains the positional information for each node in your network along with the edges, color and shape of each edge and node and any other graphical information you can think of.<br />
<br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjl7kVna7P1es7QuNISRvrP42NBzmRvZjGxZi066TwI9FZckR5ua8xsF9albSuuphqZFR_nsxl0oidDUmajlJtswUPgWZTc2Dz3IUhysgPtLMXR0XA09HY7HQ15jcRSNif4NYu31-Ib42GE/s1600/Screen+Shot+2014-10-12+at+5.25.23+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjl7kVna7P1es7QuNISRvrP42NBzmRvZjGxZi066TwI9FZckR5ua8xsF9albSuuphqZFR_nsxl0oidDUmajlJtswUPgWZTc2Dz3IUhysgPtLMXR0XA09HY7HQ15jcRSNif4NYu31-Ib42GE/s1600/Screen+Shot+2014-10-12+at+5.25.23+PM.png" height="240" width="400" /></a></div>
<br />
<br />
<br />
With the Cytoscape file in XML format we can then parse this file in python using <a href="http://www.crummy.com/software/BeautifulSoup/bs4/doc/">Beautiful Soup</a>!<br />
Beautiful Soup is an XML parser written for python. First import Beautiful Soup and tell it the format.<br />
<br />
<pre><code style="color: black; word-wrap: normal;"> from bs4 import BeautifulSoup
soup = BeautifulSoup(xgmml)
</code></pre>
<br />
You can then select all the nodes or edges from you XML file using soup.find_all<br />
<br />
<pre><code style="color: black; word-wrap: normal;">
edges = soup.find_all("edge")
print edges[0]
</code></pre>
<br />
<pre><code style="color: black; word-wrap: normal;">
<edge cy:directed="1" id="412" label="MYB46 (PD) C4H" source="145" target="123">
<att name="interaction" type="string" value="PD"></att>
<att name="shared name" type="string" value="MYB46 (PD) C4H"></att>
<att name="selected" type="boolean" value="0"></att>
<att name="name" type="string" value="MYB46 (PD) C4H"></att>
<att name="shared interaction" type="string" value="PD"></att>
</code></pre>
<br />
And Search for nodes or edges that match using with by getting the nodes label<br />
<br />
<pre><code style="color: black; word-wrap: normal;">
edge = edge[0]
print edge.get("label")
</code></pre>
<br />
<pre><code style="color: black; word-wrap: normal;">
MYB46 (PD) C4H
</code></pre>
<br />
You can then change the color of a node for example by grabbing the graphics tag and<br />
<br />
<pre><code style="color: black; word-wrap: normal;">
graphics_tag = edges.graphics
graphics_tag["fill"] = "#eb2200"
</code></pre>
<br />
The graphics tag also contain other information inducing arrow types and others shown below:<br />
<br />
<pre><code style="color: black; word-wrap: normal;">
<graphics fill="#999999" width="2.0">
<att name="EDGE_TARGET_ARROW_SHAPE" type="string" value="NONE"></att>
<att name="EDGE_LABEL_FONT_SIZE" type="string" value="10"></att>
<att name="EDGE_LABEL_COLOR" type="string" value="#000000"></att>
<att name="EDGE_LABEL_TRANSPARENCY" type="string" value="255"></att>
<att name="EDGE_SELECTED" type="string" value="false"></att>
<att name="EDGE_TOOLTIP" type="string" value=""></att>
<att name="EDGE_SOURCE_ARROW_UNSELECTED_PAINT" type="string" value="#000000"></att>
<att name="EDGE_TARGET_ARROW_UNSELECTED_PAINT" type="string" value="#999999"></att>
<att name="EDGE_LABEL_FONT_FACE" type="string" value="Dialog,plain,10"></att>
<att name="EDGE_LINE_TYPE" type="string" value="SOLID"></att>
<att name="EDGE_VISIBLE" type="string" value="true"></att>
<att name="EDGE_BEND" type="string" value=""></att>
<att name="EDGE_SOURCE_ARROW_SELECTED_PAINT" type="string" value="#ffff00"></att>
<att name="EDGE_CURVED" type="string" value="true"></att>
<att name="EDGE_SOURCE_ARROW_SHAPE" type="string" value="NONE"></att>
<att name="EDGE_TARGET_ARROW_SELECTED_PAINT" type="string" value="#ffff00"></att>
<att name="EDGE_STROKE_SELECTED_PAINT" type="string" value="#ff0000"></att>
<att name="EDGE_TRANSPARENCY" type="string" value="255"></att>
<att name="EDGE_LABEL" type="string" value=""></att>
</graphics>
</code></pre>
<br />
I have used this code for many things and made if available on github here. I used this code to color feed forward loops and stress expression response in our Xylem Secondary Cell wall Network recently published in Nature.<br />
<br />
Here is our Xylem Secondary Cell wall Network before changes:<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi8DIouIAHqvNX5uTtXMdrhT4t9dhj9OZZz7qRJrq0D2HS_YZog_JMhsjKTNBDEFhggmbPYHXYlHq5DGZvYtASX6sgvhVckXTZ3Qq667J0bVnu7sWiECVifM7k1-SKb4FnMKLRNU9QKJkaz/s1600/Screen+Shot+2014-10-12+at+6.07.28+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEi8DIouIAHqvNX5uTtXMdrhT4t9dhj9OZZz7qRJrq0D2HS_YZog_JMhsjKTNBDEFhggmbPYHXYlHq5DGZvYtASX6sgvhVckXTZ3Qq667J0bVnu7sWiECVifM7k1-SKb4FnMKLRNU9QKJkaz/s1600/Screen+Shot+2014-10-12+at+6.07.28+PM.png" height="276" width="640" /></a></div>
<br />
<br />
The feed forward loops in our network are highlighted in Red here:<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhl76ncIAmyntw_z4RA6deVUMkgS8_EEFc6eCg9TlKlxQpRayluUTCEvp2ZKFyaxuP5jMEsZaSGYH5GBRt1urQ757rx6A0zkG3IIfFboO85tr5-DrZhz7dPYx9OcqvBkX9nfvmtbXiOWT8_/s1600/Screen+Shot+2014-10-12+at+6.07.20+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhl76ncIAmyntw_z4RA6deVUMkgS8_EEFc6eCg9TlKlxQpRayluUTCEvp2ZKFyaxuP5jMEsZaSGYH5GBRt1urQ757rx6A0zkG3IIfFboO85tr5-DrZhz7dPYx9OcqvBkX9nfvmtbXiOWT8_/s1600/Screen+Shot+2014-10-12+at+6.07.20+PM.png" height="273" width="640" /></a></div>
<br />
<br />
<br />
Here is the same network with salt stress microarray data overlaid where green shows positive correlation between transcription factors/ nodes and red negative.<br />
<br />
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjjwC4M-kAqBWy_cJplwYYpGzEQSeux3QIswLOQMYSjUvZoqCwFCCR1moBqHc80cPyD2Gr0NYrI6N_AwKuomt3jszedf9eMiDJD24GI_QOCr6TEAbZnwn2ZUjO9IS4CC3rEWPJ_ELtZHLH1/s1600/Screen+Shot+2014-10-12+at+6.10.30+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjjwC4M-kAqBWy_cJplwYYpGzEQSeux3QIswLOQMYSjUvZoqCwFCCR1moBqHc80cPyD2Gr0NYrI6N_AwKuomt3jszedf9eMiDJD24GI_QOCr6TEAbZnwn2ZUjO9IS4CC3rEWPJ_ELtZHLH1/s1600/Screen+Shot+2014-10-12+at+6.10.30+PM.png" height="312" width="640" /></a></div>
<br />
<br />
<br />
<br />
<br />
<a href="http://www.nature.com/nature/journal/vaop/ncurrent/full/nature14099.html">An Arabidopsis gene regulatory network for secondary cell wall synthesis</a>
M Taylor-Teeples, L Lin, M de Lucas, G Turco, TW Toal… - Nature, 2014
<br />
<br />
<br />Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com1tag:blogger.com,1999:blog-4171086548323845827.post-26239972709410575912014-10-09T22:54:00.002-07:002015-01-18T12:23:40.549-08:00Get 23andme health information from 23andme data<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
One of the best things about 23andme is its transparency!! In addition to its awesome <a href="https://api.23andme.com/">api</a> that allows users to build tools to analyze their genetic data, the 23andme website also contains an <a href="https://www.23andme.com/you/explorer/"><span style="color: #6aa84f;">exploratory</span></a><span style="color: #cc0000;"> </span>mode which allows users to search for specific single nucleotide polymorphisms (SNPs). SNPs associated with 23andme health reports in addition to other traits are available on <a href="http://www.snpedia.com/index.php/Main_Page"><span style="color: #6aa84f;">SNPedia</span></a>. Some of these traits include SNPs for <a href="http://www.snpedia.com/index.php/Rs1446109"><span style="color: #6aa84f;">left</span></a><span style="color: #6aa84f;"> </span>handedness or even SNPs for <a href="http://www.snpedia.com/index.php/Rs53576"><span style="color: #6aa84f;">empathy</span></a>. To get this information simply search for a SNP or trait of interest on <span style="color: #6aa84f;"><a href="http://www.snpedia.com/index.php/Main_Page"><span style="color: #6aa84f;">SNPedia</span></a> </span>such as empathy: Rs53576. Then login to your 23andme and go to <a href="https://www.23andme.com/you/explorer/"><span style="color: #6aa84f;">https://www.23andme.com/you/explorer/</span></a><span style="color: #6aa84f;">.</span></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjOhidjvacg9PO1oNLcj7ief0hAvOFp790qzJXrfcrtH5Duem9dbEs5mfMPdc2xHusFfunTLUser9uhEGydo6EzvgN-7OZRhyphenhyphen7TWHPPJdCMgx0wCcsEK2t4O-R00iVl-m2eFYGnxIEpsQ9_/s1600/Screen+Shot+2014-10-09+at+10.06.51+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjOhidjvacg9PO1oNLcj7ief0hAvOFp790qzJXrfcrtH5Duem9dbEs5mfMPdc2xHusFfunTLUser9uhEGydo6EzvgN-7OZRhyphenhyphen7TWHPPJdCMgx0wCcsEK2t4O-R00iVl-m2eFYGnxIEpsQ9_/s1600/Screen+Shot+2014-10-09+at+10.06.51+PM.png" height="400" width="386" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
Here you can search for any SNP of interest like Rs53576. My genotype at this SNP is GG.</div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhD6YQH2YO3JeSpBEebn-liy67cP294gJr6xSFqsUzu1Rq1GsfoODobQzpOLXSd8b33JHuMbla-JvZDHc871w8p4D8Isv3eGQ1SZzzs4rU9asyC15IWLb1oVX723NjNEXI-GYtvoKPgsU24/s1600/Screen+Shot+2014-10-09+at+10.08.05+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhD6YQH2YO3JeSpBEebn-liy67cP294gJr6xSFqsUzu1Rq1GsfoODobQzpOLXSd8b33JHuMbla-JvZDHc871w8p4D8Isv3eGQ1SZzzs4rU9asyC15IWLb1oVX723NjNEXI-GYtvoKPgsU24/s1600/Screen+Shot+2014-10-09+at+10.08.05+PM.png" height="235" width="640" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
If we go back to <a href="http://www.snpedia.com/index.php/Main_Page">SNPedia</a> in the left hand corner we can see that having GG at the position indicates I may have my optimistic and empathetic tendencies encoded in my genome :)</div>
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgecb47w9RYtBZ1Y69mc4p9ZibO09RKbBBRFieworPg0U973yYpVAdTE6kgdTAG5uxocoypLBVBfpLSY9FmsVM-n8EwgsiwXwJnbR1uvv05dFgJoWGKCaGfca0aHEL4VETDYILdclgwqrq0/s1600/Screen+Shot+2014-10-09+at+10.38.42+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgecb47w9RYtBZ1Y69mc4p9ZibO09RKbBBRFieworPg0U973yYpVAdTE6kgdTAG5uxocoypLBVBfpLSY9FmsVM-n8EwgsiwXwJnbR1uvv05dFgJoWGKCaGfca0aHEL4VETDYILdclgwqrq0/s1600/Screen+Shot+2014-10-09+at+10.38.42+PM.png" height="191" width="400" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: left;">
Some of these SNPs arnt as well studied like the SNP for left handedness <span class="selflink" style="background-color: white; font-family: sans-serif; font-size: 14px; line-height: 22.3999996185303px;"><b><span class="selflink" style="line-height: 22.3999996185303px;"><a href="http://www.snpedia.com/index.php/Rs1446109"><span style="color: #6aa84f;">rs1446109</span></a></span><span style="color: #38761d; line-height: 22.3999996185303px;">-</span><a href="http://www.snpedia.com/index.php/Rs1007371" style="background-attachment: initial; background-clip: initial; background-image: none; background-origin: initial; background-position: initial; background-repeat: initial; background-size: initial; line-height: 22.3999996185303px; text-decoration: none;" title="Rs1007371"><span style="color: #6aa84f;">rs1007371</span></a><span style="color: #38761d; line-height: 22.3999996185303px;">-</span><a href="http://www.snpedia.com/index.php/Rs723524" style="background-attachment: initial; background-clip: initial; background-image: none; background-origin: initial; background-position: initial; background-repeat: initial; background-size: initial; line-height: 22.3999996185303px; text-decoration: none;" title="Rs723524"><span style="color: #6aa84f;">rs723524</span></a><span style="color: #6aa84f;"> </span></b></span><span class="selflink" style="background-color: white; font-family: sans-serif; font-size: 14px; line-height: 22.3999996185303px;"><span style="color: #252525;">here if you have something like AA-TT-GG you are likely to be left handed. Also at the top of the image above, if the orientation does not say "plus" but instead "minus" you need to use the </span><span style="color: #6aa84f;"><a href="http://reverse-complement.com/"><span style="color: #6aa84f;">reverse complement</span></a> </span></span><span class="selflink" style="background-color: white; color: #252525; font-family: sans-serif; font-size: 14px; line-height: 22.3999996185303px;">of the 23andme genotype to determine your genotype at that SNP. </span></div>
<div class="separator" style="clear: both; text-align: left;">
<span class="selflink" style="background-color: white; color: #252525; font-family: sans-serif; font-size: 14px; line-height: 22.3999996185303px;"><br /></span></div>
<div style="clear: both; text-align: center;">
<br /></div>
<div style="clear: both; text-align: left;">
</div>
<div style="text-align: center;">
<span style="color: #252525; font-family: sans-serif; font-size: large; line-height: 22.3999996185303px;"><br /></span></div>
<div style="text-align: center;">
<span style="background-color: white; color: #252525; font-family: sans-serif; line-height: 22.3999996185303px;"><span style="font-size: x-large;">Happy Hacking :) </span></span></div>
<br />
<div class="separator" style="clear: both; text-align: left;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<br /></div>
<br />
<div class="separator" style="clear: both; text-align: center;">
</div>
<br />Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com1tag:blogger.com,1999:blog-4171086548323845827.post-4659558018296281912014-01-20T15:50:00.003-08:002014-01-20T15:50:39.410-08:00And so it begins As a second year PhD student in Biology at UC Davis, we are required to taking one last exam before we can move on to candidacy. I take mine sometime in the next 3-6 months!!! This test is about 3 hours in length half based on classes and half on your research project. As initiative to start prepping for my exam I am going to start blogging on topics I should know and need to know. So here it goes #BlogingForMyQE (qualifying exam)Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com0tag:blogger.com,1999:blog-4171086548323845827.post-82168513174115053682013-07-20T15:11:00.003-07:002013-07-23T11:58:56.792-07:00The Five Obstacles You Will Face in Your First Year of Grad SchoolOfficially finished my first year at UC Davis as a Biochemistry, Molecular, Cellular and Developmental Biology (BMCDB) graduate student. When I think back to the start of my first year I realize I was very anxious and confused. I had heard so many horror stories about graduate school and was really not sure what to expect. Hopefully this will help with advice on the graduate school environment.<br />
<br />
<br />
<h3>
1) Finding the one (your advisor):</h3>
<br />
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjWhIcwpe0TAalGOqh2NQ-K7WdwZNE_7zI36StxoXj82e_Ep391OsD13BHOY8uRdnxVv27wcBmkx4dL3UA8rlUspfAu93Xy_HGHmRYmgL1lJjlfXItL1xAR9nLpTaptEjGV72Wy4A2xK7J8/s1600/Screen+Shot+2013-07-20+at+2.04.43+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="298" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEjWhIcwpe0TAalGOqh2NQ-K7WdwZNE_7zI36StxoXj82e_Ep391OsD13BHOY8uRdnxVv27wcBmkx4dL3UA8rlUspfAu93Xy_HGHmRYmgL1lJjlfXItL1xAR9nLpTaptEjGV72Wy4A2xK7J8/s640/Screen+Shot+2013-07-20+at+2.04.43+PM.png" width="640" /></a></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="http://www.phdcomics.com/comics/archive.php?comicid=1296">http://www.phdcomics.com/comics/archive.php?comicid=1296</a></div>
<br />
This one is stressful and tough, after all you will be working with this person for 4+ years. Even before deciding on a school make sure there are lots of options. Many of my top rotation choices did not have funding or had just left the school or the lab just wasn't what I thought it would be. Also make sure to contact older graduate students in your program and get advice on different labs and what the work ethic is like. In hindsight rotating in the summer really made a huge difference. I not only gained an extra rotation but also could easily meet with professors before the chaos began. It is best when leaving a lab to let them know you are weighting your options but really enjoyed your time there. After rotations have finished you should discuss potential thesis projects with each lab before making your decision. Just remember choose your rotations wisely, start early, and work hard :)<br />
<br />
<h3>
2) Pass your classes:</h3>
<div>
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="http://www.phdcomics.com/comics/archive/phd062104s.gif" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="275" src="http://www.phdcomics.com/comics/archive/phd062104s.gif" width="640" /></a></div>
<br />
<br />
It is actually easier then it seems to pass as long as you complete all the assignments and actually try. This does not mean that the workload is not heavy! Juggling classes, assignments and rotations is enough to make you want to nap 24/7. It's all about focusing your energy in the right place. Keep in mind that grades do not matter, while where you will be spending your Ph.D does! In other words do not stress grades too much. Form study groups, try hard and you will be fine. Part of grad school is the realizing that there are a lot of subjects that you know very little about where others in your cohort maybe experts. Within study groups show up with your work done, so you can not only defend you answer but admit to your mistakes and learn from others.<br />
<br />
<h3>
3) Get funding:</h3>
<div>
<br /></div>
<iframe allowfullscreen="" frameborder="0" height="360" src="//www.youtube.com/embed/XTvajOQ_xak" width="480"></iframe>
<br />
<br />
Start looking, writing and applying early, as in the summer before you start. I was forced to start my NSF application the summer before and was very grateful that I had. Once school starts there are so many other things to worry about. Also apply for everything even if you don't feel competitive enough. No matter how secure your funding is, many of us learned throughout rotations just how rare a professor with funding is. The best funding advice I could give is to double/triple check that all letter of rec writers will write outstanding letters and have proposals read by multiple people from varying departments and sell yourself!<br />
<br />
<h3>
4) Make Friends:</h3>
<div>
<br /></div>
<div class="separator" style="clear: both; text-align: center;">
<a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEj-iD3dQkSFNEFTFkWceyJYWMVjMBNUbn0HaDPqZlhN6_ECtm6maXBCnlEw92foAXyx-XZWh87y7IbrXWUeA3VddBSTt77zHC_sAjUZQxLqb5j7280LMj-C55qLLHszjcc4RLbdtyhSHmIx/s1600/Screen+Shot+2013-07-20+at+2.46.38+PM.png" imageanchor="1" style="margin-left: 1em; margin-right: 1em;"><img border="0" height="411" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEj-iD3dQkSFNEFTFkWceyJYWMVjMBNUbn0HaDPqZlhN6_ECtm6maXBCnlEw92foAXyx-XZWh87y7IbrXWUeA3VddBSTt77zHC_sAjUZQxLqb5j7280LMj-C55qLLHszjcc4RLbdtyhSHmIx/s640/Screen+Shot+2013-07-20+at+2.46.38+PM.png" width="640" /></a></div>
<a href="http://whatshouldwecallgradschool.tumblr.com/post/55377281333/dancing-at-grad-student-parties">http://whatshouldwecallgradschool.tumblr.com/post/55377281333/dancing-at-grad-student-parties</a><br />
<br />
<br />
Last but not least you should be making friends in graduate school, these are your soon to be collegues and possibly people you will encounter in your future. These are people you will see everyday, who encounter the same struggles, be friends not enemies. The friendship part was very unexpected for me yet helped tremendously. You are going to be stressed and confused and have worries. While it's nice at the end to celebrate your hard work with these people they are also great support throughout. These people are smart and will be able to help you throughout your graduate career. After all you all were accepted into the same program and one thing I realized over time was that everyone in our cohort has their own thing they excel at at. Also it's amazing to finally have friends just as nerdy as you! Shout out to @annaphase<br />
<br />
<br />
<br />
<b>HAVE FUN!</b><br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com0tag:blogger.com,1999:blog-4171086548323845827.post-10728712547529374102013-05-01T18:36:00.001-07:002013-07-21T21:55:25.305-07:00RNASeq Analysis: The Basics<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"></b><br />
<div dir="ltr" style="line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"><span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"><u>Strange and Mysterious File Types (you might encounter)</u></span></b></div>
<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"><br /><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"></span></b>
<br />
<div dir="ltr" style="line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"><span style="font-family: Arial; font-size: 15px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;">Sequence Read Archive format (SRA)</span><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">: This is an NCBI specific file format used because of its ability to compress read sequence information. This is often the output of many illumina sequencing pipelines.</span></b></div>
<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;">
<br /><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"></span></b><br />
<div dir="ltr" style="line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"><span style="font-family: Arial; font-size: 15px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;">Fastq file:</span><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"> SRA files can be converted to Fastq files, these are similar to Fatsta files and contain a header, associated genomic sequence and a quality score for the sequence. This is often encoded in binary and needs to be read by quality control algorithms. Thus, Fastq files contain your raw sequence information.</span></b></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
<br /><b style="font-weight: normal;"><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"></span></b></span><br />
<div dir="ltr" style="line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; font-size: 15px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;">BAM and SAM Files: </span><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">These are your alignment files where SAM stands for sequence alignment file and BAM is the binary form of this. While BAM is unreadable by humans it is often used because it is more memory efficient and quicker for computer algorithms to read in. These are obtained by often by Bowtie or Tophat after your Fastq reads have been aligned to your reference genome</span></span></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
<br /><b style="font-weight: normal;"><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"></span></b></span><br />
<div dir="ltr" style="margin-bottom: 0pt; margin-top: 0pt;">
<div style="line-height: 1.15;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; font-size: 15px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;">GTF/GFF:</span><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"><b> </b>These are often used as reference for counting how many reads map to genomic regions. These are tab delimited files containing the , start, stop, chromosome, and strand information along with name of a genomic region (such as gene, transposon, mRNA).</span></span><br />
<br /></div>
<div style="font-weight: bold; line-height: 1.15;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"><br /></span></span></div>
<div style="font-weight: bold; line-height: 1.15;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"><br /></span></span></div>
<table align="center" cellpadding="0" cellspacing="0" class="tr-caption-container" style="margin-left: auto; margin-right: auto; text-align: center;"><tbody>
<tr><td style="text-align: center;"><a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgf5DTegkgSHBOkTS5FVwtSBbtqLYcJ7d0xlmjKPF8S3gpZdh_lgR5XehllXPn8a9_A5TxIznvrPRwnAo3N6giFPh5lF9gjxF7TYaicMXOJ_yg-ylH3dZrFUePvZzt0Jww3rvbs3MBivTXb/s1600/Screen+Shot+2013-05-04+at+11.07.53+PM.png" imageanchor="1" style="margin-left: auto; margin-right: auto;"><img border="0" height="291" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEgf5DTegkgSHBOkTS5FVwtSBbtqLYcJ7d0xlmjKPF8S3gpZdh_lgR5XehllXPn8a9_A5TxIznvrPRwnAo3N6giFPh5lF9gjxF7TYaicMXOJ_yg-ylH3dZrFUePvZzt0Jww3rvbs3MBivTXb/s640/Screen+Shot+2013-05-04+at+11.07.53+PM.png" width="640" /></a></td></tr>
<tr><td class="tr-caption" style="text-align: center;">RNASeq process in flow chart on top. Bottom indicated different file formats and types of analysis for each step.</td></tr>
</tbody></table>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; vertical-align: baseline; white-space: pre-wrap;"><span style="font-size: 15px; font-weight: normal; line-height: 1.15;"></span></span><br /></span>
<br />
<div style="display: inline !important;">
<div style="display: inline !important;">
<div style="display: inline !important;">
<div dir="ltr" style="display: inline !important; margin-bottom: 0pt; margin-top: 0pt;">
<div style="display: inline !important;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; vertical-align: baseline; white-space: pre-wrap;"><span style="font-weight: normal; line-height: 1.15;"><span style="font-weight: normal;"><span style="font-weight: normal;"><span style="font-family: Arial; vertical-align: baseline;"><span style="font-size: 19px; font-weight: bold;"><u>Pre-processing:</u></span></span></span></span></span></span></span></div>
</div>
</div>
</div>
<div style="font-weight: bold;">
</div>
</div>
</div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
</span>
<br />
<div style="font-weight: bold;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">If files are obtained from NCBI’s Gene Expression Omnibus (<a href="http://www.ncbi.nlm.nih.gov/geo/delimited">geo</a>) </span><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">then they are most likely in .SRA format and need to be converted to fastq using the fastaq-dump package.</span></b><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"></b></span></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
</span>
<br />
<div style="font-weight: bold;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"></b></span></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
<br /><u style="font-family: Arial; font-size: 19px; font-weight: bold; line-height: 1.15; white-space: pre-wrap;">Quality control:</u><br />
</span><br />
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><br /></span></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
</span>
<br />
<div dir="ltr" style="margin-bottom: 0pt; margin-top: 0pt;">
Before processing the data reads must first be trimmed for adapter sequences and reads of bad quality should be filtered out. A great tool for easily viewing the distribution of quality scores is <a href="http://www.bioinformatics.babraham.ac.uk/projects/fastqc/">FASTQC</a>. I use <a href="https://github.com/tanghaibao/trimReads">trimreads</a> to both filter low quality reads and trim adaptor sequences. For more detailed analysis others have used many components of the f<a href="http://hannonlab.cshl.edu/fastx_toolkit/">astx toolkit</a>.<br />
<br /></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"></span></span><br />
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"><u>Alignment to reference genome:</u></span></span></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
<br /><span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"></span><div dir="ltr" style="line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span style="background-color: white; font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">Your fastq file contains your RNA small sequence reads. These need to be aligned back to a reference genome. </span><span style="background-color: white; font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">Reference genome</span><span style="background-color: white; font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">: Model organism you are using, sequence containing coding sequence only, microRNA only or entire genome. BOWTIE excels both in speed and memory efficiency for aligning short sequence reads back to a long reference sequence. In order to obtain this efficiency the reference genome must first be indexed. This creates a database of keys (in this case they would be cds sequences) for each record and positions them into a memory efficient tree. An index is created simply by:</span></div>
<div dir="ltr" style="font-weight: bold; line-height: 1.425; margin-bottom: 11pt; margin-top: 11pt;">
<span style="background-color: #f8f8f8; color: #333333; font-family: Consolas; font-size: 16px; vertical-align: baseline; white-space: pre-wrap;">bowtie-build reference_genome name_for_refrence_genome</span></div>
<div dir="ltr" style="line-height: 1.425; margin-bottom: 11pt; margin-top: 11pt;">
<span style="background-color: white; font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">TopHat/Bowtie can then be used to align the fastq sequences to the index reads. TopHat is often preferred in RNAseq analysis because while it is build on top of Bowtie it is also able to identify splice junctions between exons.</span></div>
<div dir="ltr" style="font-weight: bold; line-height: 1.425; margin-bottom: 11pt; margin-top: 11pt;">
<span style="background-color: #f8f8f8; color: #333333; font-family: Consolas; font-size: 16px; vertical-align: baseline; white-space: pre-wrap;">Tophat --bowtie1 name_for_refrence_genome fastq_reads</span></div>
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"><u>Determining differentially expressed gene:</u></span></div>
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"><u><br /></u></span></div>
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-family: Times; font-size: medium; font-weight: normal; line-height: normal; white-space: normal;"></b></span></div>
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-family: Times; font-size: medium; font-weight: normal; line-height: normal; white-space: normal;"><b id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b" style="font-weight: normal; line-height: normal;"></b></b></div>
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-family: Times; font-size: medium; font-weight: normal; line-height: normal; white-space: normal;"><b id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b" style="font-weight: normal; line-height: normal;"><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">Once files have been aligned to a reference genome count data can be obtained.</span></b></b></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><br /><span style="font-weight: normal;"><span style="font-family: Times; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"></span></span><div dir="ltr" style="font-family: Times; font-size: medium; font-weight: normal; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt; white-space: normal;">
<span style="font-family: Arial; font-size: 15px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;">Count Data: </span><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;">These are counts of RNASeq reads that align to a gene (or other genomic regions indicated by the GTF file). Therefore this tab-delimited file contains gene_ids and their corresponding RNA counts from your data.</span></div>
<br /><span style="font-weight: normal;"><span style="font-family: Times; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"></span></span><div dir="ltr" style="margin-bottom: 0pt; margin-top: 0pt;">
<span style="font-family: Arial, Helvetica, sans-serif; vertical-align: baseline;"><span style="font-weight: bold;"><span style="font-weight: normal; line-height: 1.15; white-space: pre-wrap;">To </span></span><span style="line-height: 17px; white-space: pre-wrap;">analyze</span><span style="font-weight: bold;"><span style="font-weight: normal; line-height: 1.15; white-space: pre-wrap;"> for differential expression two programs are commonly used DeSeq or EdgeR. You can read up on comparisons of these two here: </span></span></span></div>
</span></span></span><br />
<ul>
<li><a href="http://www.ncbi.nlm.nih.gov/pubmed/22988256"><span style="color: blue; font-family: Arial, Helvetica, sans-serif;">A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis.</span></a></li>
<li><a href="http://www.biomedcentral.com/1471-2164/13/484/"><span style="color: blue; font-family: Arial, Helvetica, sans-serif;">Efficient experimental design and analysis strategies for the detection of differential expression using RNA-Sequencing</span></a></li>
<li><a href="http://gettinggeneticsdone.blogspot.com/2012/09/deseq-vs-edger-comparison.html"><span style="color: blue; font-family: Arial, Helvetica, sans-serif;">DESeq vs edgeR Comparison</span></a></li>
</ul>
<span style="font-family: Arial, Helvetica, sans-serif;"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><br /></b></b></span>
</span><br />
<div dir="ltr" style="line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt; white-space: normal;">
<span style="font-family: Arial, Helvetica, sans-serif;"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><span style="vertical-align: baseline; white-space: pre-wrap;">Both programs assume that the number of reads in the sample can be modeled using a negative bionomial distribution.</span></span></span></span></span><br />
<div dir="ltr" style="display: inline !important; font-weight: normal; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-family: Arial, Helvetica, sans-serif;"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b">
</span></span></span></span></span></span></span></div>
<div style="font-family: Times; font-size: medium; font-weight: normal;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><br /></b></b></span></div>
</div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b">
</b></b></span>
<br />
<div style="font-family: Times; font-size: medium; font-weight: normal; line-height: normal; white-space: normal;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><span style="font-family: Arial; font-size: 15px; vertical-align: baseline; white-space: pre-wrap;"><br /></span></b></b></span></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b">
</b></b><br />
<span style="font-family: Arial; font-size: 19px; font-weight: bold; vertical-align: baseline; white-space: pre-wrap;"></span></span><br />
<div dir="ltr" style="font-weight: bold; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-weight: normal;"></b></span></div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
</span>
<br />
<div dir="ltr" style="display: inline !important; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><b id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span style="font-family: Arial; vertical-align: baseline;"><span class="Apple-tab-span"> <span style="font-size: 19px; white-space: pre;">Additional</span> <span style="font-size: 19px; white-space: pre;">Resources</span>:</span></span></b></span><br />
<b><span style="font-family: Arial; vertical-align: baseline;"><span class="Apple-tab-span"><br /></span></span></b>
<b><span style="font-family: Arial; vertical-align: baseline;"><span class="Apple-tab-span"><br /></span></span></b></div>
<div style="font-weight: bold;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"></span><br />
<div dir="ltr" style="display: inline !important; font-weight: normal; line-height: 1.15; margin-bottom: 0pt; margin-top: 0pt;">
<div style="margin: 0px;">
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de" style="font-family: Arial, Helvetica, sans-serif;"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"></span></span></span></span></span></span></span></span></div>
</div>
<div style="-webkit-text-stroke-width: 0px; color: black; font-family: Times; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 18px; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px;">
</div>
<br />
<div style="-webkit-text-stroke-width: 0px; color: black; font-family: Times; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 18px; margin: 0px; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px;">
<span style="font-family: Arial, Helvetica, sans-serif;"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de"><span id="docs-internal-guid-345044b6-62d3-cae1-6da5-f8ad6de0ac8b"><span style="vertical-align: baseline; white-space: pre-wrap;"><a href="http://sfg.stanford.edu/quality.html">Amazing Webpage for RNASeq help</a></span></span></span></span></span></div>
<div style="-webkit-text-stroke-width: 0px; color: black; font-family: Times; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 18px; margin: 0px; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px;">
<span style="font-family: Arial, Helvetica, sans-serif;"><a href="https://github.com/BradyLab/RNASeq">My Github page</a></span></div>
<div style="-webkit-text-stroke-width: 0px; color: black; font-family: Times; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 18px; margin: 0px; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px;">
<span style="font-family: Arial, Helvetica, sans-serif;"><span style="vertical-align: baseline; white-space: pre-wrap;"><br /></span></span></div>
<div style="-webkit-text-stroke-width: 0px; color: black; font-family: Times; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: 18px; margin: 0px; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px;">
<br /></div>
</div>
<span id="docs-internal-guid-1a253e31-62d1-4a1a-1d36-06cc738ed9de">
</span>Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com5tag:blogger.com,1999:blog-4171086548323845827.post-66323145720077005212013-03-26T00:25:00.001-07:002013-03-26T00:25:02.514-07:005 Apps All Graduate Students Should Have<h4>
<span style="color: #444444; font-size: x-large;">1) <a href="http://www.6wunderkinder.com/wunderlist">Wunderlist</a></span></h4>
<strong><br /></strong>
The one thing I have repeatedly learned in graduate school is that you will always be adding to your TODO list. This list is never ending and often bogs me down. Wunderlist works well for me because it allows me to set deadlines and has a friendly interface. I keep everything from school work to twitter post on it =)<br />
<br />
<strong><br /></strong>
<iframe frameborder="0" height="315" src="http://www.youtube.com/embed/PhYruJX4ek8" width="560"></iframe><br />
<strong><br /></strong>
<br />
<h4>
<strong><span style="color: #444444; font-size: x-large;">2)<a href="http://scholar.google.com/"> Google Scholar </a></span></strong></h4>
<strong><br /></strong>
If you have published a paper you should make a Google Scholar account. If you haven't published yet you should make a Google Scholar too! Google Scholar is like a Facebook for scientist, it allows others to easily find papers published by you and who you have collaborated with. It also recommends papers for you based on papers you have read and what you have published. <br />
<br />
<br />
<br />
<strong><img alt="image" src="http://media.tumblr.com/2b9b621d189048aa860db5d308cad123/tumblr_inline_mk994e2tTW1qz4rgp.png" /></strong><br />
<strong><br /></strong>
<br />
<h4>
<strong><span style="color: #444444; font-size: x-large;">3) <a href="http://brettterpstra.com/projects/nvalt/">nvALT</a></span></strong></h4>
<strong><br /></strong>
You are not always going to have time to take notes, sometimes you just want to jot something down quickly. This app allows users to take quick notes and find them easily in the future. Also it allows you to write (and preview) your notes in HTML and Markdown so that information can be easily be transferred to a website or maybe even a blog.<br />
<br />
<br />
<strong><img alt="image" src="http://media.tumblr.com/ed7deb6138a822351dc1e8a46a22b1fb/tumblr_inline_mk99ms5IgL1qz4rgp.png" /></strong><br />
<strong><br /></strong>
<br />
<h4>
<strong><span style="color: #444444; font-size: x-large;">4) <a href="https://www.dropbox.com/">Dropbox</a>/Google Drive</span></strong></h4>
<br />
If you have not started using Dropbox or Google Docs yet, get on it! This allows you to access your information from any computer anywhere any time. That research paper you were putting off reading, put off no more. Another bonus to these services is that information is stored in the cloud, free backup. Although you do not ever think your computer will crash and 4 hours of work will be destroyed, it will happen, TRUST ME!<br />
<strong><img alt="image" src="http://media.tumblr.com/a9a59d1e98a10af8e3a40ed917498f7f/tumblr_inline_mk99t0HVLK1qz4rgp.jpg" /></strong><br />
<strong><br /></strong>
<br />
<h4>
<strong><br /></strong></h4>
<h4>
<span style="color: #444444; font-size: x-large;"><strong>5) <a href="https://chrome.google.com/webstore/detail/stayfocusd/laankejkbhbdhmipfmgcngdelahlfoji?hl=en)">StayFocused Chrome Extension</a></strong> </span></h4>
Also available for firefox and safari browsers. This app allows you to spend a limited amount of time on distracting sites (such as Facebook) before telling you to get back to work and blocking you from that site for the rest of the night.<br />
<br />
<br />
<strong> <img alt="image" src="http://media.tumblr.com/250e0229ac235779aab53b491455a26c/tumblr_inline_mk996deoZ01qz4rgp.png" /></strong><br />
<strong><br /></strong>
<strong><br /></strong>
<strong><br /></strong>
<strong><br /></strong>
<strong><br /></strong>Anonymoushttp://www.blogger.com/profile/01673256848690551590noreply@blogger.com0