Aug 15 - Sept 9, 2019 Salmon-sea lice BS seq analysis

Sea lice bismark alignments:

Read More

Wed. Jul. 10, 2019 Salmon-Sea lice Zymo Pico Methyl prep + RRBS digest repeat

Zymo Pico Methyl kit prep

  • set up programs on PCR machine in 209 (called ZYM1,ZYM2, ZYM4, and ZYM5 under STRIGG folder)
  • Into a 48-well plate, I aliquoted out 20-50ng DNA and up to 20uL of nanopure H2O according to this sheet
  • I added 50ng of sea lice DNA:
    • Sea lice Female 1 = 67.2ng/uL (qubit HS, 7/10)
      • added 0.75uL + 19.25uL nanopure H2O
    • Sea lice Female 2 = 17ng/ul (qubit HS, 7/10)
      • added 3uL + 17uL nanopure H2O
Read More

May 8, 2019 Geoduck broodstock hemolymph DMR calling

Running methylkit:

I had planned to run MethylKit on Emu’s RStudio web interface to find DMRs in the WGBS hemolymph. In reviewing a previous MethylKit R script, I realized I need sorted bam files, which my initial bismark mox script didn’t include. So I submitted a new mox job to sort my bams. Once the job finishes (after Mox maintenance), I will copy the sorted bams over to gannet and then on to Emu so I can run the Rstudio web interface (a little painful, but getting there!).

Read More

Thurs. Mar 14, 2019, Oyster Seed Proteomics

Mapping to SR lab GO slim terms

  • I was unable to get simantic similarity for these terms (which is needed in order to relate proteins to one another through their terms) because they don’t map to terms in the goslim_generic.obo file or to the go data that ontologyX parses. So I’m not going to use these.

Cleaning up analysis for poster figs

Read More

Wed. Feb 20, 2019, Geoduck genome versions, histology scoring, pt. whitney supplies

Geoduck gonad histology scoring:

  • worked with Kaitlyn to figure out scoring for geoduck gonad histology
    • we looked at Robert Marshall’s dissertation(he was in the Pearce lab at UBC), Grace’s geoduck gonad project, and Molly’s manuscript (from Brent) for some criteria on scoring .
    • so far, Kaitlyn has recorded the data here
    • for pictures and to find a scope with a ruler for measuring oocyte diameter, I contacted Lindsay in Jackie’s lab to see if we can use theirs.

Collected supplies to bring to Pt.Whitney tomorrow

  • heath stacks
  • liquid nitrogen in the dry shipper
  • pipette tips
  • DI water
  • Mercuric chloride waste container
  • microfuge tube racks (3)
  • 1 box large gloves
Read More

Fri. Jan 11, 2019, Oyster Seed Proteomics

In trying to run NMDS analysis on technical replicate ADJNSAF data, I found discrepencies between the ADJNSAF values in Steven’s ABACUS_output021417NSAF.tsv and Sean’s Abacus_output.tsv. I compared Steven’s ABACUS_output021417.tsv file (from which he made ABACUS_output021417NSAF.tsv, see his jupyter notebook with Sean’s Abacus_output.tsv and found no difference:

Read More

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