Wednesday, 04/15/15

9:00 am – 10:00 pm


Code development

  • discuss Jeff’s new approach to Library construction
  • Fix indexing in Assemble Oligo Array probes.
    • complete. see Library11CompileOligoArrayResults_150414.m
  • BLAST combos of primers and new Fast Secondaries
    • completed for 5′ fusions. killed half of ~500 remaining secondaries.
    • now doing for 5′ and 3′ rev-complement primer fusions.
  • library 12 passed checks finally! Sent to Jeff and team
  • lib12 sent to CustomArray, requested PO

Ph project

  • revising text
  • writing model description
  • writing figure caption
  • adjusting figure
  • sent to NF and AW
Posted in Summaries | Comments Off on Wednesday, 04/15/15

Tuesday 04/14/15

10:00 am – 8:20 pm, 9:20pm – 12:00am


  • reply to Quake lab with raw data
  • working on fixing lib12 build
  • found bug in my matlab BLAST — allhits field not being populated if no varargin’s passed. fixed.
  • updating MERFISH team on discussions with collaborators, 12:00pm-12:40pm
  • adjusting Library building script to allow buffering of total probe numbers 1:30pm – 3:00pm
  • need to shuffle pri-secondary combos and secondaries when changing to new secondaries. (primary encoding sequences will be the same).
  • running script to assemble library 12, 3:00pm – 4:00pm

adding remove probes so as to balance bits 4:00pm – 7:00 pm

  • better to do based on bit combos
  • potential bug in probe assembly: probe names are non-unique. Not certain targeting sequences selected at random are being selected without replacement… Need to investigate this.
  • Investigated: these are from different 1kb sections, that happen by chance to have the same base end. Hence our sequence IDs are non-unique. should fix this for ID purposes, but makes no difference for probes.

Running probe building (again), 7:00pm – 12:00am

  • serious problems with FAST libraries losing whole bits or large number of bits when fused to primers in the final BLAST.
  • wrote loop to run through primers if we don’t keep at least 75% of the least abundant bit combo (since everyone will be pruned down to the least abundant).
  • FAST secondaries + lib5 primers are getting destroyed by BLAST to abundant genes and rRNA/tRNA.
  • memory is getting destroyed by BLAST
  • launching on Morgan, 10:15 pm
  • Writing to TSTORMdata2 is at least 10x slower than writing straight to disk, even from Morgan. Nope. 10x was too optimistic. Monet lapped Morgan in printing libraries. TSTORMdata2 is probably closer to 100x slower.
  • On Monet, added breakpoint to start writing code to ID primers and secondaries that are highly repeated and killed by BLAST to abundant genes.
  • BLAST is extremely slow for this purpose. (30 min per attempted primer and secondary set). Investigating BLAT alternative.
  • probe names are non-unique.
    • should have recorded which pt number of the sequence the probe was from and this would not have been a problem
    • tried adding Tm. – no luck, still non-unique
    • adding instead random 12-number integer id per probe.
    • everything is unique. BLAST still refuses to build a database. darn it.

checking probe numbers:

  • variable length 4490 > 1, 5374 > .1, Tm 66-76
  • 35-55% GC content gives 4057 of 8750 with FPKM > 1

Multiplex chromatin

  • discussion with Bogdan on strategy 11:30 am 12:00 pm

Discussion with XZ

  • discussed setup distribution

Ph project

  • reran simulations on cluster size in Ph-ML to match expression range. Finished much faster this time.
  • zoomed in on sub-panels
Posted in Summaries | Comments Off on Tuesday 04/14/15

SDB 2015 abstract

Super-resolution Imaging of Chromatin Nanostructure
Links Epigenetic State and 3D Genome Packaging

Metazoan genomes are packaged at multiple scales, and regulation of this spatial organization may play an important role in the developmental and cell fate selection. Unfortunately, current methods provide little information about 3D organization of the genome at the length scale of genes (kilobases) and regulatory domains (hundreds of kilobases) in single cells. I will present a new super-resolution imaging approach to study the structural organization of the genome at the kilobase to megabase scale in individual cells at 20 nm resolution. These domains largely occupy diffraction limited volumes and thus their structures cannot be resolved with conventional imaging approaches. From thousands of images of dozens of epigenetic domains from across the Drosophila genome, we have discovered, within a single cell type, a substantial diversity of structural patterns: compact and diffuse domains, branched and linear domains, domains that are highly entangled with one-another and domains which are strictly segregated. These different structural features are closely correlated to certain differences in the epigenetic state of the chromatin. I will focus on the organization of Polycomb bound domains, which exhibit a surprising, entangled structure and length-dependent compaction. Computational models suggest this organization could contribute to the repressive nature of the domains. Preliminary comparisons between developing tissues suggests differential regulation of chromatin structure associated with developmental fate selection. This work suggests further super-resolution imaging studies of chromatin structure at this scale may greatly aid our understanding of the role of genome structural regulation in development.

Posted in Summaries | Comments Off on SDB 2015 abstract

Monday 04/13/15

9:00 am – 6:20 pm, 8:00 pm- 10:30 pm


  • design new libraries
  • discussed design with team (9:30-10:30am)
  • working on library assembly (8:00pm-10:00pm)

New Libraries (140 gene).

  1. 20mer secondaries 30mer targeting regions 1-16 (shorter than before)
  2. 20mer secondaries 30mer targeting regions 17-32
  3. FAST 20mer secondaries 30mer targeting regions 1-16
  4. FAST 20mer secondaries 30mer targeting regions 1-16
  5. 20mer secondaries 6 on-bits 30mer targeting 1-16. (2 bits per target)
  6. 20mer secondaries 30mer targeting regions 1-16 (4 bits per target)

For later

  1. 20mer secondaries variable targeting

Chromatin tracing

  • Prep for chromatin tracing meeting 10:40 – 2:00pm
  • discussion with Bogdan on strategy 11am-12pm
  • see slides
  • tracing discussion with TW and XZ and team 2:00pm – 6:15 pm

Ph project

  • new draft of figures and text from NF

Abstract Submission

  • submitted abstract for SDB meeting, ~2 hours before deadline (I think). 10:00-10:30pm
Posted in Summaries | Comments Off on Monday 04/13/15

Saturday 04/11/15

Ph Project

  • working on simulations

Assume certain amount of open sites. All will be fill with no Ph.
Add PhM. The same fraction of PhM will be bound as the fraction of Ph bound, since the kinetics are the same.

start with 500 unbound and 1500 bound with no PhM, all sites occupied, eq constant 7:1 bound
as we add PhM, it fills and kills sites, only need 1 per site to kill it (if capped early). Sites tend to either fill to max or get capped as created.
Strong transition point at PhM = PhE/2.


  • improved rendering: scatter is proportional to cluster size
  • should zoom in more — difficult because matlab makes chain look disconnected where it goes outside of frame.
  • should look up how to set zoom by command line (manual scroll does not work smoothly as it should).


  • See email to Nicole and Ajaz describing today’s results.
Posted in Summaries | Comments Off on Saturday 04/11/15

Friday 04/10/15

10:00 am – 10:15 pm

group meeting 10am-1pm

MERFISH library design

  • would like a new human transcriptome database, build off of multiple cell lines.
  • Steven has ENCODE lines downloaded and cufflinks run. Just need to write matlab script to parse data

Ideas for increasing probe count:

  • allow flexible probe length

Ph simulations

  • wt can not join a cluster with equal Ph.
  • fixed bug in modification (empty sites still need to allow Ph to join)
  • started working on prettier model rendering for paper figures
Posted in Summaries | Comments Off on Friday 04/10/15

MERFISH published! (in Science Magazine)

Posted in Genomics | Comments Off on MERFISH published! (in Science Magazine)

Wednesday 04/08/15

10:00 am – 8:00pm

10 am – 11 am: Ph simulations

  • migrating Ph simulations onto Morgan.
  • running Ph simulations
  • maybe much shorter equilibration time and just many more rounds of simulation would be better
  • running lots of short simulations across a broader range of Ph. This round wont finish before tomorrow.

11:00 am – 12:45 am

  • revising MERFISH proofs

MERFISH meeting and probe design 1pm – 5pm

  • 20mers 66C – 76C not enough

OligoArray batch launch issues

  • read-write speed of fasta files seems to be limiting new launch — 30 jobs quickly launched but none complete
  • troubleshooting OligoArray speed
  • at crosshybeT of 99, still get lots of ‘nonspecific oligo ignored at this time. Later get
  • launched a 40 to 50% GC content version of scan (this is going to get less than 1000 genes) [no other stringency parameters]
  • launched a 15 to 75% GC content (this isn’t much of a cut).
  • stopped previous 20 mer building at 35869 of 37442 probes. (lib11 prep folder: 15-04-03_Probes_20mers_t66Cx99Cs99C).

Probe numbers

  • Previous 30mers, FPKM > 0.1 and >100 probes = 2197
  • New 20mers, FPKM > 0.1 and >100 probes = 3929
  • Previous 30mers, FPKM > 1 and >100 probes = 1856
  • New 20mers, FPKM > 1 and >100 probes = 3026
  • New 20mers, FPKM > 1 and >50 probes = 6920 (4 2ndaries per probe we still have 50 per hybe).
  • New 20mers cuts: Tm 66 to 76.
  • Of 37608 genes in our database, 12,016 are longer than 2 kb (i.e. 100 probes if we just uniformly tile with 20mers). Not sure why I am losing so many…

5pm – 8pm: more Ph simulations

  • need to explain limited binding sites (K27 model) and its effects on cluster size distribution.
  • I think actually binding is too strong, data suggests there are fewer
Posted in Summaries | Comments Off on Wednesday 04/08/15

Tuesday 04/07/15

9:00 am – 11:45 pm


  • revise Text for chromatin paper
  • Extended data figures for chromatin entanglement (done)
  • revise Fig 4
  • Histone images for Doory (sent)
  • Revised model figure for Ajaz
  • write abstract for ASBMB on MERFISH (see post)
  • write abstract for Soc Dev Bio on chromatin

9am – 2pm

  • revising chromatin text
  • polishing abstract with Bogdan, 10:30am – 11:30am

5pm – 9pm

  • entanglement supp figures
  • revising figure 4
    • update colors, labels, graph sizes
    • add images of subdomains
    • cleaned up model fig code
  • rewrote caption
  • updated into text, now version 4-4.

9pm – 11:45 pm

  • waiting for BB to read MS
  • should we change to all volume no Rg?
  • generated semilog-y plots for internal scaling data and linear plots. All in dropbox\Figures\Version7\Figure2
  • sent to XZ

Data Analysis

  • running E04 E05 data through Chromatin Cropper
  • errors in automatic load after splitdax. needs troubleshooting

Potential conferences

    • Ren, Zhang, Wysocka
  • Keystone symposia (but nothing looks relevant. This would have been good:
Posted in Summaries | Comments Off on Tuesday 04/07/15

Protected: draft ASBMB 2015 abstract

This content is password protected. To view it please enter your password below:

Posted in Conference Notes | Comments Off on Protected: draft ASBMB 2015 abstract