<!DOCTYPE html> <html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta charset="utf-8" /> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <meta name="generator" content="pandoc" /> <title>Home</title> <script src="site_libs/jquery-1.11.3/jquery.min.js"></script> <meta name="viewport" content="width=device-width, initial-scale=1" /> <link href="site_libs/bootstrap-3.3.5/css/cosmo.min.css" rel="stylesheet" /> <script src="site_libs/bootstrap-3.3.5/js/bootstrap.min.js"></script> <script src="site_libs/bootstrap-3.3.5/shim/html5shiv.min.js"></script> <script src="site_libs/bootstrap-3.3.5/shim/respond.min.js"></script> <script src="site_libs/navigation-1.1/tabsets.js"></script> <link href="site_libs/highlightjs-9.12.0/textmate.css" rel="stylesheet" /> <script src="site_libs/highlightjs-9.12.0/highlight.js"></script> <link href="site_libs/font-awesome-4.5.0/css/font-awesome.min.css" rel="stylesheet" /> <style type="text/css">code{white-space: pre;}</style> <style type="text/css"> pre:not([class]) { background-color: white; } </style> <script type="text/javascript"> if (window.hljs) { hljs.configure({languages: []}); hljs.initHighlightingOnLoad(); if (document.readyState && document.readyState === "complete") { window.setTimeout(function() { hljs.initHighlighting(); }, 0); } } </script> <style type="text/css"> h1 { font-size: 34px; } h1.title { font-size: 38px; } h2 { font-size: 30px; } h3 { font-size: 24px; } h4 { font-size: 18px; } h5 { font-size: 16px; } h6 { font-size: 12px; } .table th:not([align]) { text-align: left; } </style> </head> <body> <style type = "text/css"> .main-container { max-width: 940px; margin-left: auto; margin-right: auto; } code { color: inherit; background-color: rgba(0, 0, 0, 0.04); } img { max-width:100%; height: auto; } .tabbed-pane { padding-top: 12px; } button.code-folding-btn:focus { outline: none; } </style> <style type="text/css"> /* padding for bootstrap navbar */ body { padding-top: 51px; padding-bottom: 40px; } /* offset scroll position for anchor links (for fixed navbar) */ .section h1 { padding-top: 56px; margin-top: -56px; } .section h2 { padding-top: 56px; margin-top: -56px; } .section h3 { padding-top: 56px; margin-top: -56px; } .section h4 { padding-top: 56px; margin-top: -56px; } .section h5 { padding-top: 56px; margin-top: -56px; } .section h6 { padding-top: 56px; margin-top: -56px; } </style> <script> // manage active state of menu based on current page $(document).ready(function () { // active menu anchor href = window.location.pathname href = href.substr(href.lastIndexOf('/') + 1) if (href === "") href = "index.html"; var menuAnchor = $('a[href="' + href + '"]'); // mark it active menuAnchor.parent().addClass('active'); // if it's got a parent navbar menu mark it active as well menuAnchor.closest('li.dropdown').addClass('active'); }); </script> <div class="container-fluid main-container"> <!-- tabsets --> <script> $(document).ready(function () { window.buildTabsets("TOC"); }); </script> <!-- code folding --> <div class="navbar navbar-default navbar-fixed-top" role="navigation"> <div class="container"> <div class="navbar-header"> <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar"> <span class="icon-bar"></span> <span class="icon-bar"></span> <span class="icon-bar"></span> </button> <a class="navbar-brand" href="index.html">fucci-seq</a> </div> <div id="navbar" class="navbar-collapse collapse"> <ul class="nav navbar-nav"> <li> <a href="index.html">Home</a> </li> <li> <a href="about.html">About</a> </li> <li> <a href="license.html">License</a> </li> </ul> <ul class="nav navbar-nav navbar-right"> <li> <a href="https://github.com/jdblischak/fucci-seq"> <span class="fa fa-github"></span> </a> </li> </ul> </div><!--/.nav-collapse --> </div><!--/.container --> </div><!--/.navbar --> <!-- Add a small amount of space between sections. --> <style type="text/css"> div.section { padding-top: 12px; } </style> <div class="fluid-row" id="header"> <h1 class="title toc-ignore">Home</h1> </div> <div id="TOC"> <ul> <li><a href="#overview">Overview</a></li> <li><a href="#microscopy-image-analysis">Microscopy image analysis</a></li> <li><a href="#rna-seq-data-preprcessing">RNA-seq data preprcessing</a></li> <li><a href="#cell-cycle-signal-in-gene-expression-data">Cell cycle signal in gene expression data</a></li> <li><a href="#model-fitting">Model fitting</a></li> <li><a href="#one-time-investigations">One-time investigations</a></li> </ul> </div> <hr /> <div id="overview" class="section level2"> <h2>Overview</h2> <ul> <li><a href="data-overview.html">Data description</a></li> </ul> <hr /> </div> <div id="microscopy-image-analysis" class="section level2"> <h2>Microscopy image analysis</h2> <ul> <li><a href="images-process.html">Processing images - from images to intensities</a></li> </ul> <p>We evaluated and pre-processed the results of image analysis as follows:</p> <ol style="list-style-type: decimal"> <li>We visually inspect images deteced to have none or more than one nucleus. For cases that are inconsistent with visual inspection, we correct the number of nuclei detected. <ul> <li><a href="images-multiple-nuclei.html">Inspect images with multiple nuclei</a></li> <li><a href="images-zero-nuclei.html">Inspect images with no nucleus</a></li> </ul></li> <li>We applied background correction to the intensity measurements of GFP, RFP and DAPI based on the following analyses. <ul> <li><a href="confess-prelim.html">CONFESS results</a></li> <li><a href="images-qc.html">QC analysis including no. nuclei detected, DAPI, and intensity variation</a></li> <li><a href="images-qc-followup.html">Explore using log10 sum pixel intensity for signal metrics</a></li> <li><a href="images-metrics.html">Compare correction approaches using median versus mean background</a></li> <li><a href="images-metrics-cell-shape.html">Explore associations between nucleus shape metrics vs intensities</a></li> </ul></li> <li>We analyzed intensity variation across individuals and batches and determined on an approach that removes batch effect in the data. <ul> <li><a href="images-qc-labels.html">Visualize signal variation by plate and individual identity</a></li> <li><a href="images-qc-variation.html">Visualize the structure of signal variation by individual identity</a></li> <li><a href="images-normalize-quantile.html">Quantile normalization for GFP, RFP and DAPI</a></li> <li><a href="images-normalize-anova.html">Estimate variance explained in IBD and correct for batch effects in intensities</a></li> </ul></li> </ol> <hr /> </div> <div id="rna-seq-data-preprcessing" class="section level2"> <h2>RNA-seq data preprcessing</h2> <ol style="list-style-type: decimal"> <li>The first step in preprocessing RNA-seq data consists of QC and filtering. <ul> <li>Sample QC and filtering <ul> <li><a href="sampleqc.html">Sample QC criteria</a><br /> </li> <li><a href="totals.html">Sequencing depth</a><br /> </li> <li><a href="reads-v-molecules.html">Reads versus molecules</a></li> </ul></li> <li>Gene QC and filtering <ul> <li><a href="gene-filtering.html">gene filtering</a></li> <li><a href="pca-tf.html">PCA with technical fators</a></li> </ul></li> </ul></li> <li>We then analyzed and corrected for batch effect due to C1 plate in the sequencing data <ul> <li><a href="seqdata-batch-correction.html">Estimate variance explained in IBD and correct for batch effects</a></li> </ul></li> </ol> <hr /> </div> <div id="cell-cycle-signal-in-gene-expression-data" class="section level2"> <h2>Cell cycle signal in gene expression data</h2> <ol style="list-style-type: decimal"> <li>We investigated cell cycle signals in the sequencing data alone.</li> </ol> <ul> <li>Consider <a href="images-transgene.html">transgene count in sequencing data</a></li> <li>Compute linear correlation between gene expression levels <a href="images-seq-correlation.html">with DAPI and FUCCI intensities</a></li> </ul> <ol start="2" style="list-style-type: decimal"> <li>We then assign categorical labels of cell cycle and explored the expresson profiles of these categories.</li> </ol> <ul> <li>Cluster samples by <a href="images-pam.html">Partition around medoids(PAM)</a></li> <li>Cluster samples by <a href="images-mclust.html">Guassian mixture modeling</a></li> <li>Select a subset of samples that are closet to the cluster centers (cluster representives) <a href="images-subset-silhouette.html">using silhouette index</a></li> <li>Examine gene pression scores defined by the Macosko paper in the selected cluster representatives <a href="images-classify-fucci.html">before confounding correction</a> and <a href="images-classify-fucci-adjusted.html">after confounding correction</a></li> <li>Examine gene expression scores defined by the Macosk paer <a href="images-classify-leng.html">in the sorted cells of the Leng et al. 2015 paper</a></li> </ul> <ol start="3" style="list-style-type: decimal"> <li>We ordered cells on a circle using FUCCI intensities alone.</li> </ol> <ul> <li>First, I used GFP and RFP intensities to estimate a <a href="images-circle-ordering.html">least-square fit of unit circle which approximate the relative ordering of cells on cell cycle</a></li> <li>I computed the PCs of GFP and RFP and used these to infer the relative ordering on a circle (i.e., polar coordinates), and <a href="images-circle-ordering-eval.html">to evaluate the circle fit, I computed circular-circular and circular-linear correlation with DAPI and gene expression</a></li> <li>Here I put together lists of genes identified as signfiicantly correlatd with the PC-based fit <a href="images-circle-ordering-sigcorgenes.html">by linear correlation and circular-linear correlation, also cyclical genes by smash</a></li> </ul> <hr /> </div> <div id="model-fitting" class="section level2"> <h2>Model fitting</h2> <ul> <li>CellcycleR 0.1.6 <ul> <li><a href="images-cellcycleR-convergence.html">Model convergence assessment</a></li> <li><a href="images-cellcycleR.html">Fitting on intensities across plates and individuals</a></li> <li><a href="cellcycler-seqdata-leng.html">Fitting on Leng data)</a></li> <li><a href="cellcycler-seqdata-fucci.html">Fitting on fucci-seq RNA-seq data)</a></li> </ul></li> <li><p>Correlation for circular response variable: some simulations to check <a href="circ-correlation-simulation.html">its property and also relations with Fisher’s z transformation</a></p></li> <li>Nonparametric circular regression <ul> <li>First, evaluate smash for smoothing gene expression data and compare <a href="images-circle-ordering-npreg.html">results of smash with NPCirc</a></li> </ul></li> </ul> <hr /> </div> <div id="one-time-investigations" class="section level2"> <h2>One-time investigations</h2> <ul> <li><p>Why some gene symbols (genes) correspond to <a href="ensembl.html">multiple Ensembl IDs?</a></p></li> <li><p>I selected a set of cell cycle <a href="seqdata-select-cellcyclegenes.html">genes that belong to GO term Cell Cycle and looked at the overlap with the detected genes in our data.</a></p></li> </ul> </div> <!-- Adjust MathJax settings so that all math formulae are shown using TeX fonts only; see http://docs.mathjax.org/en/latest/configuration.html. This will make the presentation more consistent at the cost of the webpage sometimes taking slightly longer to load. Note that this only works because the footer is added to webpages before the MathJax javascript. --> <script type="text/x-mathjax-config"> MathJax.Hub.Config({ "HTML-CSS": { availableFonts: ["TeX"] } }); </script> <hr> <p> This <a href="http://rmarkdown.rstudio.com">R Markdown</a> site was created with <a href="https://github.com/jdblischak/workflowr">workflowr</a> </p> <hr> <!-- To enable disqus, uncomment the section below and provide your disqus_shortname --> <!-- disqus <div id="disqus_thread"></div> <script type="text/javascript"> /* * * CONFIGURATION VARIABLES: EDIT BEFORE PASTING INTO YOUR WEBPAGE * * */ var disqus_shortname = 'rmarkdown'; // required: replace example with your forum shortname /* * * DON'T EDIT BELOW THIS LINE * * */ (function() { var dsq = document.createElement('script'); dsq.type = 'text/javascript'; dsq.async = true; dsq.src = '//' + disqus_shortname + '.disqus.com/embed.js'; (document.getElementsByTagName('head')[0] || document.getElementsByTagName('body')[0]).appendChild(dsq); })(); </script> <noscript>Please enable JavaScript to view the <a href="http://disqus.com/?ref_noscript">comments powered by Disqus.</a></noscript> <a href="http://disqus.com" class="dsq-brlink">comments powered by <span class="logo-disqus">Disqus</span></a> --> </div> <script> // add bootstrap table styles to pandoc tables function bootstrapStylePandocTables() { $('tr.header').parent('thead').parent('table').addClass('table table-condensed'); } $(document).ready(function () { bootstrapStylePandocTables(); }); </script> <!-- dynamically load mathjax for compatibility with self-contained --> <script> (function () { var script = document.createElement("script"); script.type = "text/javascript"; script.src = "https://mathjax.rstudio.com/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"; document.getElementsByTagName("head")[0].appendChild(script); })(); </script> </body> </html>