To account for N/A Cq values, I substituted in a Cq value of 45 wherever the Cq value was nonexistent. This allows for some initial analysis of the data to see which genes might warrant further work. Below are the adjusted qPCR plots. I also generated the COX1 Cq plot (not normalized to actin) to examine relative mitochondrial abundance between stressed and non-stressed diploids and triploids.
Daily Archives: January 31
Shelly’s Notebook: Wed. Jan 30, 2019
Geoduck Broodstock Experiment
Gonad slices sent to HCS, Inc. for histology
List of samples sent: https://drive.google.com/open?id=1Fgb5abuMR8TyiR9hbIjUW2mRg9eeoAd3pD3CdlPlqlo
Order form: https://drive.google.com/open?id=1UXb3Vc2scDMOwDsE5kWIyt8ozCF-buhj
Shipped through FedEx as hazardous goods of expected quantities, each of the 5 cassette jars contained 30 mL of 70% Ethanol (Stabilizer from PAXgene kit).
Low pH conditioning plans
Steven and I discussed stopping the treatment on Friday and combining the low pH group crates (Tank 1 and 2) into one ambient tote and Tank 3 and 4 crates into one ambient tote. This will make room for more broodstock and will hopefully rescue the low pH conditioned animals so we can get them to spawn/strip spawn. See slack discussion.
Oyster Seed Proteomics
Map the proteins to 2019 Uniprot accessions
For re-mapping the fasta to 2019 Uniprot accessions like Steven did in the past, I started the jupyter notebook linked below on Ostrich and it’s currently still running:
Github desktop on Ostrich
Also got github desktop to work on Ostrich by downloading this older version that I found here. With the newest version, the graphics don’t work when you remote desktop in. But this older version works fine, probably because Ostrich is still running El Cap OS.
from shellytrigg http://bit.ly/2RuO5sE
via IFTTT
Shelly’s Notebook: Tues. Jan 30, 2019
Concentrating Geoduck Broodstock Hemocyte DNA for WGBS
Because the minimum concentration for WGBS with Genewiz is 20ng/ul, I pooled and precipitated the DNA samples from yesterday as follows:
- I combined tank 3 samples 15 and 16, total volume ~200uL
- I combined tank 2 samples 025 and 026, total volume ~200uL
- I added 140uL isopropanol and 20uL 3M Sodium Acetate pH 5.2 to each ~200uL pooled sample
- I vortexed to mix
- I spun at 12000 rpm for 30 min at 4C
- I removed the supernatent and added 500uL of 70% ethanol to each
- I spun tubes for 10 min at 12000 rpm at 4C
- I removed the supernatent and resuspended in 20uL EB (elution buffer from the EZNA kit, which I think is 10mM Tris, pH 8.5)
Qubit concentrations: Standard 2 = 100ng/ul
Animal ID | Tank | Treatment | Conc. (ng/ul) | Total DNA (ng) |
---|---|---|---|---|
15_16 | 3 | amb | 284 | 1230 |
025_026 | 2 | low (pH 6.8) | 53 | 830 |
We are going to send all of Tank 2 ‘025_026’ DNA and half of Tank 3 ‘15_16’ DNA, which Sam now has and planning to ship on Monday. The remaining Tank 3 ‘15_16’ DNA is in my box labeled ‘Shelly Trigg start date 10/11/2018’ in the -20C in rm 213.
Oyster seed proteomics time x temperature
Creating differential network visualizations
To create differential network visualizations, we need:
- protein relationships
- this could be GO annotations, KEGG pathways, or even protein interaction information (although I’m not sure how conserved this is between species)
- a metric to quantify abundance changes
- nodes could be colored by fold change and/or effect size to show up or down
- node size could be relative to p-value, or even PCA loadings values
Calculating log fold change and p-values
Emma suggested calculating the ratio of ‘NUMSPECSTOT’ of each protein to the sum of NUMSPECSTOT for the sample. Then calculating log fold change. And using a proportions test for calculating the p-value.
I think all fold changes should be relative to day zero so that we can visualize what protein networks are active at each stage of development and compare differences between the two temperatures (sort of like what these guys did). For a figure, what I’m picturing is a 2 x 6 grid for the two different temperatures and 6 different days. And a network in each box of the grid showing proteins colored differently based on their fold change and sizes based on p-value significance.
Here is the R markdown file I made for calculating these. I used a Chi square proportions test for the p-values. So far I have only done one comparison (23C_day 3 to day 0).
Need to:
- make a loop to do the rest of the comparisons
- map the proteins to Uniprot accessions using an updated uniprot DB.
- figure out what proteins are connected either by GO or physical interaction information
from shellytrigg http://bit.ly/2Ww7Zra
via IFTTT