Tuesday, April 26, 2016

QUBEHUB software for simulations

Including RStudio, XPPAUT, PPLane etc.


PLAS power law analysis and simulation

Voit book, software support for ODE analysis


Atlanta QBIO book orders

Books to order:

A first course in systems biology, By Eberhard Voit
ISBN: 9780815344674, Garland

Physical models of living systems,  by Philip Nelson
ISBN-10: 1-4641-4029-4; ISBN-13: 978-1-4641-4029-7;

Mathematics for the life sciences,  by Erin N. Bodine, Suzanne Lenhart, Louis Gross
ISBN-13: 9780691150727, Publisher: Princeton University Press
A Biologist's Guide to Mathematical Modeling in Ecology and Evolution, by Sarah P. Otto Troy Day , ISBN: 781400840915,  Publisher: Princeton University Press

Mathematical Modeling in Systems Biology: An Introduction, by Brian P Ingalls,
Publisher: The MIT Press; 1 edition (July 5, 2013), ISBN-10: 0262018888, ISBN-13: 978-0262018883

Data Wise, Revised and Expanded Edition: A Step-by-Step Guide to Using Assessment Results to Improve Teaching and Learning


Discipline-Based Education Research: A Guide for Scientists Paperback – July 16, 2015



-Schneider B, Carnoy M, Kilpatrick J, Schmidt WH, & Shavelson R. (2007). Estimating causal effects: Using experimental and observational designs. American Educational Research Association: Washington DC.
-Weimer, M. (2006). Enhancing scholarly work on teaching and learning. Jossey-Bass: San Francisco.

(to read) Kaya, Ma, sanger strain lifespan study

 Defining Molecular Basis for Longevity Traits in Natural Yeast Isolates
Alaattin Kaya1,#, Siming Ma1,#, Brian Wasko2, Mitchell Lee2, Matt Kaeberlein2, and Vadim N. Gladyshev1,*
1Division of Genetics, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA, 02115, USA

RLS were provided

todo: download data, and Sanger resource data

Growth rates were determined using a Bioscreen C MBR machine by analysis of optical density in the OD420-580 range as previously described in combination with the YODA Software package18. The data on transcripts, peptides (proteins), metabolites and morphology were downloaded from Yeast Resource Center http://www.yeastrc.org/g2p/download.do. Values corresponding to the 22 strains were extracted; metabolite data were not available for 378604X. Metabolites with missing values in more than one strain (other than 378604X) were discarded; the remaining missing values (6 out of 107 metabolites) were imputed based on 10 nearest neighbors, using “knnImputation” function of R package “DMwR”. For comparison across the phenotypic data, the values were standardized across the strain by setting mean = 0 and standard deviation = 1. In addition, for genes represented by multiple peptides, we calculated the mean standardized values to perform the regression.

(to read) Yang 2015 PNAS yeast GFP screen, asymmetric partition

GFP library screen, Cy5 labeling of cell wall, daughters have few cy5 because their walls are newly synthesized.

Most proteins are symmetrically distributed.

All aging factors are 'non-essential'.  This means that their 'interactions' to essential genes should be of importance in my network aging model.

Thursday, April 21, 2016

yeast DNA repair genes

Gene Ontology Term: DNA repair

Biological Process


protein  NP_010242.1



desktop flow cytometer survey

Main fluophores: 
  Propidum Iodide (PI), dihydroethydium (DHE), and dihydrorhodamine (DHR)

My understanding of the basic configuration are (based on CALIBUR):
 blue laser 488nm + red laser 635nm (optional)
Detector channels: FSC, SSC, FL1 green 530nm, FL2 yellow 585, FL3 red 670nm, with an optional  FL4 661nm. 


Vendors and models

BioRad S3e Cell Sorter




ACEA NovoCyte 1000 $40K ?
Novote Cyte model 2000, 3000 have more optins.

BeckmanFlow.com CytoFlex

Accuri C6PLUS  $49K

EmdMillipore Muse: $15K
Con: Closed system. Cannot be modified ?

Attune NxT Flow Cytomer

Stratedigm FLow Cytometer


Sales rep from BD on Gallios and CytoFlex
As per my voice message, I understand you are looking for a quote for 6 color Gallios Flow Cytometer, but we now have a new, much smaller flow cytometer with greater capability than the Gallios, and less expensive.   It is called the CytoFLEX Flow Cytometer.   The CytoFLEX is available in 21 different configurations from a 1 laser/4 color system to a 3 laser/13 color system in a small box (16in x 16in x 13in).   It is completely upgradable to any configuration at any time.   A 96 well plate loader option is available as are 561 nm and 375 nm laser options.   Independent testing has proven this instrument to be one of the most sensitive instruments in the market today.
 If the 561 or 375 nm lasers interest you, please let me know and I will forward additional information on our CytoFLEX “S” System.

Wednesday, April 20, 2016

bio125 20160420 flow cytometer data analysis in R,

Section 3
video recording

9:10-10:10am, Flow cytometer data analysis in R using my own laptop. Helped some students with RStudio package installations.

Problems: In Windows 10, Rstudio has to be run as administrator to install packages.

10:15am, course evaluation
10:25am post-computing survey
10:35am lab 11.1_group

section 4

1-2pm. R Rstudio on flow data anlysis
Problems: R3.2.3. installation.  a/s/n warning.   Jumping lines. Setting working direcories.

2:11pm. Post survey

by 3pm.  review bioinfor_1

Tuesday, April 19, 2016

toread, yeast RLS screen paper


toread, Yeast longevity promoted by reversing aging-associated decline in heavy isotope content

Yeast longevity promoted by reversing aging-associated decline in heavy isotope content

stochastic modeling of histone modification in yeast

mating phenotype maintenance?

CR SIR2 histone modification modeling

J Xing's similar work on histone modification

CR, rapamycin effect on LOH in CLS and H2O2 treatment, NIH R15?

CR, rapamycin effect on LOH in CLS and H2O2 treatment

Human essential genes?

Mining  shalem 14 data set

References: NEST in in shirley liu's lab.

Project Achilles

Monday, April 18, 2016

*** human essential gene

NEST, sherley liu lab,

CRISP screen paper from Feng Zhang lab:

Try to load the xlsx file into Rstudio. I waited for more than 50 minutes on Byte (4 G RAM), I then have to kill it. I then converted the xlsx to csv, and it worked in less than 1 minute!!!

> length(unique(tb$sgRNA.sequence))
[1] 64751
There are 64.7K CRISP shots.

Only 18.7K genes are tagged. So, the missing one contain essential genes. Some other criteria are needed. 
> length(unique(tb$Gene.name))
[1] 18736

.gitignore file



Resazurin Cell Viability assay


Elastic net regression

Elastic net regression


***, compound profiling, initial thoughts

These compounds are 'newly designed and synthesized' and their targets are needed to be identified and verified.

PCA, then, Euclean distance,  MCL
PCA, then K-means

For comparison: Elastic net

Score criteria: jaccard index with original annotation.


Friday, April 8, 2016

bio125 cancer and ucsc genome broswer

Section 3:
2015 bio125 assignments are based on GRC37/hg19 assembled in 2009.
4 exercises in class
10:30am, class is over.

Wednesday, April 6, 2016

bio125 HU treatment of yeast cells

Get the HU experiment start as quickly as possible.
During the 1.5 hour incubaion time,
check FOA plate, take pictures
go over group exercises. 

Section 3:
Projector is set at a wrong station. A technician came and figure it out. 

During the incubation, flip tube occasionally to make sure cells are suspended. 

HU made cells arrested at G1 phase with round shapes. 
Section 4:

Problems: Because we did not sonicate the cells, so often cells in G1 phases are often chained together.  It is easy for students to count big circle with tiny dots as S phases, and everthing else as unbudded.
https://youtu.be/yxA1-MvwZfQ  2014 HU video

Tuesday, April 5, 2016

bio325 guest lecture, yeast oxidative stress analysis

Try to lead students work on Rmd file for yeast gene expression data analysis

It turns out that the bioconductor packages are not supported by the recent R3.2.4. I had to lead students to use R3.2.3.

Among the students, there are 7 Apple laptops and 2 windows laptops.

I made a mistake when going through probes on Affymetrix Yeast 2.0 chip. It contains 2 MSH2 genes, on for S cerevisiae and one for S pombe.

Saturday, April 2, 2016

23andMe in courses

23andMe in courses

Udacity course

MSH2, DNA mismatch repair, and sporadic colon cancer

BIO125 used MSH2
A Howard U and Johns Hopkins study on sporadic colon cancer and DNA mismatch repair

Primer design to amplify your DNA signature, identify RE to distinguish the polymorphic sites.

open human, data source


Personal Genome Project,

Personal Genome Project


Why the data vary from 5M to 200M?

spelman biology curricula

Competence-based core curriculum. reading, writing, talking and doing science.

SAWOK, science as a way of knowing


Gene and genealogy resources

23andme Data Format, technology

The file 23andme provides has four columns: 
rs ID, chromosome, position, and genotype

Beadchip for genotyping

Bioconductor 23andme


Friday, April 1, 2016

bio233 order list fall 2016

Wish list for bio233 orders

** marker pens
** 10 tally counters
** 95% enthanol large order for 70% enthanol bottle
**latex-free gloves
** lens cleaners

swim noodles
colored bacteria for edutainment

cover slips, fix

slides for gram stains: spirilum +/-, bacillus +/-, cocci +/-

x brightline hemacytometer chambers

scotch tapes
play doughs

nutrient plate agars

small bottles for ethanol. Idodine. Crystal violet.

Yeast selective media, drop-out media

bibulous papers
filter papers

BPB pH color indicator for fermentation lab. 

BD Calibrite beads

15ml falcon tubes

2016 research day preparation

Research posters:
Maya Jones
Faith Lyons
Maya Bryant
Jessica Corley

BIO125 course posters
Gaina-Yvvan Pierre, Elsie, Unique Hayes
Jaliyah Peterson, Lady Nwdike

Oral presentations:
 Kierra Parker
 Imani-Michelle White

bio125, 20160401Fri FOA

Materials: FOA plates, cells, serial water, Eppendorf tubes.

Section 3:
Ask the purpose of URA3 and FOA,
9:20, draw a diagram on the board. 
Go over previous student mistakes, not mix cells well, label the cover, adding cells to the cover. 
9:30, students started the lab. Students should have 9 tubes, and can be balanced triangularly in centrifuge. 
9:47. Students cannot see the pellet for OD=0.01. I remind them to orient their tubes the same way, so they can know where the pellet positions, even they cannot see the pellet.
One group only labeled OD on the tube but not strain names, and are confused after centrifuging
by 10:30, four groups are still working on the lab, because they have to redo the dilutions.
Forgot to ask students to label their section, and group on plates.

Section 4:
Emphasize that tubes during spin should be arranged 

One group was not sure that tips need to be changed each time for a different dilution.