Wed. Nov 14, 2018

Geoduck Broodstock Experiment

FFAR meeting notes

  • Can’t untangle effects from high/low pH or static/fluctuating, so include additional fluctuating low pH treatment with new geoduck
    • Hollie has a code for this
Read More

Thur. Nov 1, 2018

Methylation analysis:

Geoduck alignments too slow, should be directional, and many chromosomal sequence extraction errors

Read More

Mon. Oct 29, 2018

from mox rsync –archive –progress –verbose 20181004/

Read More

Sun. Oct 28, 2018

Methylation analysis: P. generosa juvenile OA response

bismark currently running on full data set

mapping and coverage analysis with subsetted samples (10K reads each)

bismark summary
- all samples have > 50% mapping efficiency with Bismark settings –score_min L,0,-1.2 -I 60 –non_directional
- deduplication removed very few reads, per usual
- % methylation is on average 20%

methylkit summary - coverage plots don’t appear to obviously show any particular sample or condition should be excluded due to low coverage.this may become more obvious in the full analysis or downstream analysis
- the genome region > 60000 seems to have low coverage in all samples
plot of genome regions with no coverage organized by length of treatment

Read More

Tues. Oct 9, 2018

ASCA analysis with oyster proteomics temperature x time series dataset

  • Goal: To determine if proteins drive difference between temperature groups over time
    • we need to understand the model that ASCA creates from the data so we want to extract the loadings from the temperature factor ASCA and try to understand the behavior of the proteins with high loadings values
      • PC1 explains 100% of the variation; this seems strange so need to understand what is up with that
Read More

Mon. Oct 1, 2018

My goals today

  1. Get Husky card
  2. Get bus pass
  3. Geoduck meeting @ 1:30p
    • notes from meeting:
      • How frequently to sample
      • Why are animals limited to 80?
        • because feeding of live algae is limiting
          • benefits of live algae vs. paste
            • paste better because doesn’t mess with pH and we can have more animals that way
Read More

Friday. Sept 28, 2018

##My goals for today

  • work on Bismark alignment comparison
    • how does Bismark compare to Methylpy?
      • reference papers:
      • jupyter notebook
      • install methylpy
      • try bowtie vs. minimap
      • make genomes
    • Insert size alignment setting
      • in bowtie, this is the -I parameter
        • default is 0; in methylpy it’s 20.
  • Roberto’s last day lunch @12:30
Read More

Thurs. Sept 27, 2018

My goals for today

  • FISH 546 9-10
  • Epigenetics reading group @ NOAA NWFSC 12:30-1:30
    • transgenerational of envronmental info
    • thoughts: stress opens heterochromatin => transposon/pseudo gene?/protogene? expression; does the cell just go crazy under stress and hope something works? Or is there a non-random mechanism that selects for advantageous changes?
    • siRNAs silence TEs: piRNAs specificly target transposons and this is heritable
    • should we be doing ChIPseq of H3K9me (heterochromatin mark) to complement methylation data if H3K9me marks change under stress like they do in the paper?
  • TGIT 4-6
Read More

Wed. Sept 26, 2018

My goals for today

  • Work on Bismark alignments
    • Trying different alignment score settings:
      • Looked into bowtie2 alignment manual
      • alignement score is dependent on matches and mismatches
        • alignment score of 0 = 100% perfect match
        • 1 mismatch = -6 score
        • 1 ambinguous character = 1
        • if there’s a gap in the reference seq. or the read seq: 5 + 3*gap length
          • a 2bp gap in the read seq. = 2 insertions (relative to the read seq.)
          • a 2bp gap in the ref. seq. = 2 deletions (relative to the read seq.)
      • default in normal alignment mode is L, 0, -0.2 (where L = linear, the first number = the intercept, and the last number = the slope, and x = read length). The default specifies for reads that are 100bp the alignment score must be >= -20.
        • in –local alignment mode, G is used and G = natural log (i.e. f(x) = 1.0 + 5.4ln(x), where x is the length of the read). The min score default in local mode is 26 for 100 bp read.
Read More

Tues. Sept 25, 2018

My goals for today

  • Get personal items from Mukilteo field station
    • done
  • Farewell lunch at noon
    • done
  • Get everyone’s thoughts on water chemistry data table for crab metabolomics paper
    • done
Read More

Mon. Sept 24, 2018

My goals for today

  1. Lab meeting @ 10
  2. Meet with Steven to decide projects to focus on
  3. SAFS kickoff @ 2-5

Notes from lab meeting

  • Sam’s working on oyster transcriptome and genome assembly
  • Kaitlyn’s working on oyster proteomics temperature x time series analysis
    • pharmacokinetics toolkit applicable? (Brent asked)
    • RNA-seq time series analysis tool applicable?
  • Yaamini’s working on reviewer comments for JSR submission on oyster broodstock DMR comparison under OA conditions
    • what is the best way to vaildate methylation levels?
      • MSP vs. MBD or both? Need data relevant to physiology
  • Brent’s working on lots of grants and coordinating Pt. Whitney setup
    • Pt. Whitney: need to set a date, set up CO2 tanks
      • monitoring Geoduck over time by
        • measuring respiration weekly
        • hemolymph sampling (how contaminated with sea water will this be?)
          • where’s the best place to poke?
            • heart (Brent tried, not best survival rate)
            • aductor muscle (maybe too small, size of an eraser)
            • siphon
          • ELISA for vitellogenin to tell dev. stage
            • can we order an commercial ELISA kit to test on Geoduck?
            • can we test kit with old hemolymph sample where vitellogenin was expressed?
          • [SRM]
        • non-lethal gonad biopsy (Emma’s paper)
        • histology? at least on some for validation of dev. stage
  • Laura’s working on summer histology analysis on OWL (protocol on Git). Austrailia work was temp. association with spawning and using histology to quantify eggs to tell stage: brooder vs. non-brooder. Also doing OA Oly. RNAseq data analysis pipeline
  • Emma’s working on hatchery microbiome with OA conditioned Geoduck. Comparing traditional assembly to MOCAT to 6-Gill and digging into unaligned reads. Larvae at pH 8.2 had ciliates in seawater, ciliates didn’t grow at lower pH treatments.
Read More

Wed. Sept 19, 2018

Progress on first Bismark analysis:

-Found a way to do my analysis/write files to the server but keep my Jupyter notebook in my github desktop folder see Jupyter notebook
-still working on code for generating Bismark summary files

Read More