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
- discussion with Bogdan on strategy 11:30 am 12:00 pm
Discussion with XZ
- discussed setup distribution
- reran simulations on cluster size in Ph-ML to match expression range. Finished much faster this time.
- zoomed in on sub-panels
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.
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).
- 20mer secondaries 30mer targeting regions 1-16 (shorter than before)
- 20mer secondaries 30mer targeting regions 17-32
- FAST 20mer secondaries 30mer targeting regions 1-16
- FAST 20mer secondaries 30mer targeting regions 1-16
- 20mer secondaries 6 on-bits 30mer targeting 1-16. (2 bits per target)
- 20mer secondaries 30mer targeting regions 1-16 (4 bits per target)
- 20mer secondaries variable targeting
- 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
- new draft of figures and text from NF
- submitted abstract for SDB meeting, ~2 hours before deadline (I think). 10:00-10:30pm
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.