忍者ブログ

[PR]

×

[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。

書けん日

科研費が書けん日

という、しょうもないけど消耗しすぎて必ず頭をよぎるギャグがあるのですが。

もちろん今日も科研費が書けん日なんですが。
科研費が書ける立場であることが実はとても恵まれているというのも事実なので、もうひと踏ん張りしたいと思います。

PR

PSMCの前処理

毎度おなじみ、僕です。
このブログが生きていることに驚いています。

1個体のゲノムから個体群動態を推定するPSMCという方法。
https://github.com/lh3/psmc

今ではいろんな発展系もありますが、やっぱり理論的にも使い方的にもシンプルでわかりやすい(相対的に)原作を使いたくなるときもあります。

日本語解説はこちらが詳しいです。
https://sites.google.com/site/hiromimatsumae/過去の情報/集団遺伝学


が、新しいsamtools系だとこの通りにはいかず、古いバージョンを併用する必要が出てきて面倒です。

特に公式のこの部分。
"""
samtools mpileup -C50 -uf ref.fa aln.bam | bcftools view -c - | vcfutils.pl vcf2fq -d 10 -D 100 | gzip > diploid.fq.gz
"""

これの新しいバージョンに対応したコマンドは以下です。
"""
bcftools mpileup -C50 -Ou -f ref.fa aln.bam | bcftools call -c - | vcfutils.pl vcf2fq -d 10 -D 100 | gzip > diploid.fq.gz
"""
bcftoolsで代用できました。

お試しあれ。

フローサイトメトリー用語解説

フローサイトメトリーによるゲノムサイズや倍数性の推定に関する専門用語
Estimation of nuclear DNA content in plants using flow cytometry
Jaroslav Dolezel, Johann Greilhuber & Jan Suda
Nature protocol (2007)

haplophase 単相世代(染色体数 n) 蘚苔類によくある。
diplophase 複相世代(染色体数 2n)維管束植物によくある。

haplophaseの個体の複製してない核のDNA含有量は"1C"。染色体数nのハプロイドゲノムのDNA量にあたる。
DNA複製後haplophaseの核のDNA含有量は"2C"。diplophase個体の複製していない核と同じ。
endopolyploidy核内倍数性は普通Cのレベルで表す。

基本染色体数がxのmonoploidのゲノムサイズはCx-valueで表す。
6倍体植物パンコムギの体細胞は2n=6x=42で、2C-valueは6Cx-valueと等しい。

haploid, diploidなど染色体数についてのワードは、DNA含有量を示すのに使うべきではない。
また、フローサイトメトリーで得られる倍数性データは染色体カウントで得られるものとは区別すべきだ。
核DNA含有量から倍数性/異数性に言及するなら、"DNA 倍数性"、"DNA 異数性"とすべきだ。

菌類の分類

今日は菌の分類について、以下の論文のアブストだけ読んだのでメモ。

Fungal evolution: diversity, taxonomy and phylogeny of the Fungi
Biol. Rev. (2019)

菌類のゲノムを使った系統解析から、以下の9つの門に再定義しました。日本語はwikipediaを参照しています。

Opisthosporidia
Chytridiomycota ツボカビ門
Neocallimastigomycota ネオカリマスティクス門
Blastocladiomycota コウマクノウキン門
Zoopagomycota
Mucoromycota ケカビ亜門
Glomeromycota グロムス門
Ascomycota 子嚢菌門
Basidiomycota 担子菌門

日本語のないものが新しく分けられた分類群でしょうか。
論文には各門の特徴や代表的な種について書いてあるので気が向いたら追記します。

転職します

転職?異動?します。

クソ環境からは抜け出さなければならない。
続けてあげる義理もない。

奴隷にしたいんだったら給料の桁がたりんのよ。

一桁多かったら奴隷になるか?それもどうだろう。ならないかもな。

カレンダー

11 2024/12 01
S M T W T F S
1 2 3 4 5 6 7
8 9 10 11 12 13 14
15 16 17 18 19 20 21
22 23 24 25 26 27 28
29 30 31

フリーエリア

最新コメント

プロフィール

HN:
No Name Ninja
性別:
非公開

バーコード

ブログ内検索

P R