NBDC Research ID: hum0214.v9

 

研究内容の概要

目的: 免疫細胞における遺伝子発現制御と免疫疾患に対する寄与を解明する

方法:

JGAS000220(JGAD000309、JGAD000310):全身性強皮症21症例、antineutrophil cytoplasmic antibody(ANCA)関連血管炎26症例、対照健常者28名の末梢血免疫細胞画分を分取し各画分から抽出したtotal RNAを用いたRNA-seqを実施した。

E-GEAD-397 / E-GEAD-398 / E-GEAD-420:全身性エリテマトーデス、筋炎、全身性強皮症、混合性結合組織病、シェーグレン症候群、関節リウマチ、ベーチェット病、成人スチル病、ANCA関連血管炎、高安動脈炎、対照健常者の全血及び末梢血免疫細胞28画分を分取し、全血より全ゲノムシーケンス、各画分サンプルを用いてRNA-seqを実施した。サンプルのフィルタリングと遺伝子発現量の標準化の上で、各免疫細胞画分でeQTL解析を実施した。

JGAS000296:全身性強皮症50症例、対照健常者48名の末梢血免疫細胞24画分を分取し、各画分サンプルを用いてRNA-seqを実施した。遺伝子発現量を定量し、サンプルのフィルタリングを行った。

JGAS000220(JGAD000371、JGAD000372、JGAD000373):全身性エリテマトーデス、対照健常者の末梢血免疫細胞19画分を分取し、各画分サンプルを用いてRNA-seqを実施した。末梢血免疫細胞15画分でATAC-seqを実施した。

JGAS000486:全身性エリテマトーデス、対照健常者の末梢血免疫細胞27画分を分取し、各画分サンプルを用いてRNA-seqを実施した。

JGAS000598:関節リウマチ50症例、対照健常者39名の末梢血免疫細胞18画分を分取し、各画分サンプルを用いてRNA-seqを実施した。関節リウマチ15症例については、経時的な解析を実施した。

JGAS000485:全身性エリテマトーデス、全身性強皮症、筋炎、混合性結合組織病、シェーグレン症候群、関節リウマチ、ベーチェット病、成人スチル病、ANCA関連血管炎、高安動脈炎、対照健常者の末梢血B細胞5画分を分取した。各画分サンプルを用いてRNA-seqを実施し、B細胞受容体配列のマッピングを実施した。

JGAS000626:関節リウマチ19症例、全身性エリテマトーデス62症例、特発性炎症性筋疾患40症例、対照健常者64名の末梢血免疫細胞9画分を分取し、各画分サンプルを用いてRNA-seqを実施した。遺伝子発現量を定量し、サンプルのフィルタリングを行った。

JGAS000627:全身性エリテマトーデス136症例、対照健常者89名の末梢血免疫細胞27画分を分取し、各画分サンプルを用いてRNA-seqを実施した。遺伝子発現量を定量し、サンプルのフィルタリングを行った。

JGAS000648:炎症性筋疾患6症例、対照健常者1名から末梢血・筋組織・気管支洗浄液を採取した。末梢血についてはCD4T細胞を分取し、scRNA-seqを実施し遺伝子発現量を定量した。筋組織・気管支洗浄液についてはCD45陽性細胞を分取し、scRNA-seqを実施し遺伝子発現量を定量した上でCD4Tクラスターに属する細胞のみを抽出した。

対象: 全身性強皮症、全身性エリテマトーデス、筋炎、混合性結合組織病、シェーグレン症候群、関節リウマチ、ベーチェット病、成人スチル病、ANCA関連血管炎、高安動脈炎、対照健常者

URL: http://ryumachi.umin.jp/

 

データID内容制限公開日
JGAS000220 NGS(RNA-seq:全身性強皮症) 制限公開(Type I) 2020/10/09
JGAS000220 NGS(RNA-seq:ANCA関連血管炎) 制限公開(Type I) 2021/03/05
E-GEAD-397 RNA-seqのリードカウント 非制限公開 2021/04/28
E-GEAD-398 Conditional eQTLサマリーデータ(有意データのみ) 非制限公開 2021/04/28
E-GEAD-420 Nominal eQTLデータ(全データ) 非制限公開 2021/04/28
JGAS000296 NGS(RNA-seq) 制限公開(Type I) 2022/01/21
JGAS000220 NGS(RNA-seq)全身性エリテマトーデス 制限公開(Type I) 2022/03/09
JGAS000220 NGS(ATAC-seq) 制限公開(Type I) 2022/03/09
JGAS000486 NGS(RNA-seq)解析によるリードカウントデータ 制限公開(Type I) 2022/08/25
JGAS000598 NGS(RNA-seq)解析によるリードカウントデータ 制限公開(Type I) 2023/03/28
JGAS000485 RNA-seqによるB細胞免疫細胞サブセットのB細胞受容体レパトアデータ 制限公開(Type I) 2023/07/06
JGAS000626 NGS(RNA-seq)解析によるリードカウントデータ 制限公開(Type I) 2024/02/08
JGAS000627 NGS(RNA-seq)解析によるリードカウントデータ 制限公開(Type I) 2024/02/08
JGAS000648 NGS(scRNA-seq)解析によるカウントデータ 制限公開(Type I) 2024/02/08

※リリース情報はこちら

※制限公開データの利用にあたっては、利用申請が必要です。申請方法はこちら

※ 論文等でデータベースからダウンロードしたデータを含む結果を公表する際には、NBDCヒトデータベースに登録されたデータを利用した旨について謝辞(Acknowledgement)に記載して下さい。記載例はこちら

 

分子データ

JGAS000220(RNA-seq:全身性強皮症)

対象

全身性強皮症(ICD10:M340):21症例

対照健常者:13名

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞19画分(Naive_B, SM_B, USM_B, DN_B, Plasmablast, Th1, Th2, Th17, Tfh, Naive_CD4, Mem_CD4, Fr._II_eTreg, Naive_CD8, Mem_CD8, mDC, pDC, CD16p_Mono, CD16n_Mono, NK)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp
マッピング方法 STAR(hg38)
フィルタリング(QC)方法 アダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。50塩基未満のリード、及び低クオリティーの塩基を多く含むリード(20%を超える塩基で、Phred quality score < 20)を除去した。ユニークマップ率が80%未満、もしくはユニークマップリード数が5.00 x 106未満のサンプルは除去した。同一サブセット内の遺伝子発現の相関係数(Di)を計算した。Diが0.9未満のサンプルは除去した。
遺伝子数 26353
Japanese Genotype-phenotype Archive Dataset ID JGAD000309
総データ量 125 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000220(RNA-seq:ANCA関連血管炎)

対象

ANCA関連血管炎(ICD10:M318):26症例

対照健常者:28名(上記13名を含む)

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞20画分(Naive_B, SM_B, USM_B, DN_B, Plasmablast, Th1, Th2, Th17, Tfh, Naive_CD4, Mem_CD4, Fr._II_eTreg, Naive_CD8, Mem_CD8, mDC, pDC, CD16p_Mono, CD16n_Mono, NK, Neu)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp
マッピング方法 STAR(hg38)
フィルタリング(QC)方法 アダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。50塩基未満のリード、及び低クオリティーの塩基を多く含むリード(20%を超える塩基で、Phred quality score < 20)を除去した。ユニークマップ率が80%未満、もしくはユニークマップリード数が5.00 x 106未満のサンプルは除去した。同一サブセット内の遺伝子発現の相関係数(Di)を計算した。Diが0.9未満のサンプルは除去した。
遺伝子数 26353
Japanese Genotype-phenotype Archive Dataset ID JGAD000310
総データ量 125 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

E-GEAD-397(RNA-seq)

対象

全身性エリテマトーデス(ICD10:M329):62症例

筋炎(ICD10:M339、M332):65症例

全身性強皮症(ICD10:M340):67症例

混合性結合組織病(ICD10:M351):19症例

シェーグレン症候群(ICD10:M350):18症例

関節リウマチ(ICD10:M0690):25症例

ベーチェット病(ICD10:M352):23症例

成人スチル病(ICD10:M0610):18症例

ANCA関連血管炎(ICD10:M318):26症例

高安動脈炎(ICD10:M314):16症例

対照健常者:92名

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞28画分(Naive_CD4, Mem_CD4, Fr._I_nTreg, Fr._II_eTreg, Fr._III_T, Th1, Th2, Th17, Tfh, NK, Naive_CD8, Mem_CD8, EM_CD8, CM_CD8, TEMRA_CD8, Naive_B, USM_B, SM_B, DN_B, Plasmablast, CL_Mono(別名CD16n_Mono), CD16p_Mono, Int_Mono, NC_Mono, mDC, pDC, LDG, Neu)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp
マッピング方法 STAR(GRCh38)
フィルタリング(QC)方法 cutadapt (v1.16)およびfastx-toolkit(v0.0.14)を用いてアダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。低クオリティーの塩基が20%を超えるリードを除去した。STAR(v2.5.3)のtwo-pass modeを用いてGRCh38へのマッピングを行った(Gencode version 27アノテーション)。ユニークマップ率が90%未満(Plasmablastでは70%未満、他のB細胞サブセットでは<85%未満)、もしくはユニークマップリード数が6.00 x 10^6未満のサンプルは除去した。遺伝子発現はHTSeq(v 0.11.2)で定量した。解析に不適切なサンプルの除外基準としては、各細胞分画で低発現遺伝子(90%を超えるサンプルで10リード未満)を除外し、edgeRのTMMノーマライゼーションを実施、ログcount per million(CPM)変換、ComBatによるバッチ除去を行った上で、同一サブセット内の遺伝子発現のSpearman相関係数を計算し、平均値が0.9未満のサンプルを除外した。
遺伝子数 53344
Genomic Expression Archive ID E-GEAD-397
総データ量 1.3 GB(疾患名、phase情報、年齢、性別データおよび各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

E-GEAD-398 / E-GEAD-420(eQTL)

対象

全身性エリテマトーデス(ICD10:M329):62症例

筋炎(ICD10:M339、M332):65症例

全身性強皮症(ICD10:M340):67症例

混合性結合組織病(ICD10:M351):19症例

シェーグレン症候群(ICD10:M350):18症例

関節リウマチ(ICD10:M0690):24症例

ベーチェット病(ICD10:M352):23症例

成人スチル病(ICD10:M0610):18症例

ANCA関連血管炎(ICD10:M318):25症例

高安動脈炎(ICD10:M314):16症例

対照健常者:79名

規模 eQTL
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq X Ten]
ライブラリソース 全血より抽出したDNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) TruSeq DNA PCR-Free Library prep kit
断片化の方法 TruSeq DNA PCR-Free Library prep kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 151 bp
マッピング方法 BWA-MEM(GRCh38)
フィルタリング(QC)方法 全ゲノムシーケンスデータはGATK(v 4.0.6.0)best practiceに準じて解析した。検出バリアント率が99%未満のサンプルは除外した。遺伝型のインピュテーションにはBEAGLE(v 5.1)を用いた。Call rate <85%, Hardy-Weinberg平衡検定P値<1.0x10-6, マイナーアレル頻度1%未満のバリアントは除外した。
遺伝子数

E-GEAD-398:19441

E-GEAD-420:22381

eQTL検出方法 各細胞サブセットで80%を超えるサンプルで5リード未満、もしくは0.5 CPM未満の低発現遺伝子を除外した。常染色体遺伝子発現データのTMMノーマライゼーション、CPM変換、inverse normal transformによるサンプル間ノーマライゼーションを実施した。また隠れた共変量の推測のためPEERを実施した。eQTL解析の共変量として、上位2つの遺伝子主成分、サンプル回収のphase、性別、PEERによる潜在的因子を含めた。Phase 1で回収したMem CD8はPhase 2ではCM CD8とEM CD8として分割しているため、大部分の細胞の含まれるEM CD8と統合してeQTL解析を実施した。各細胞サブセットでのConditional eQTL解析では、QTLtoolsのpermutation passによる10000回のパーミュテーションを実施し、各遺伝子レベルでのFDR <0.05と対応するnominal p-value thresholdを定めた。その上で、QTL toolsのconditional passを用いてforward-backward stepwise regression eQTL analysisを実施した。Nominal eQTL解析では QTL toolsのnominal passを用いて、各遺伝子の転写開始点から1M塩基対のバリアントと遺伝子発現の関連を解析した。
Genomic Expression Archive ID

E-GEAD-398

E-GEAD-420(2021/6/9: Added columns for REF/ALT allele. Slope indicates the effect size of alternative alleles.)

総データ量

E-GEAD-398: 3.9 GB(eQTLサマリーデータ [FDR<0.05]、txt)

E-GEAD-420: 38 GB(全eQTLサマリーデータ、txt)

コメント(利用にあたっての制限事項) NBDC policy

 

JGAD000296(RNA-seq)

対象

全身性強皮症(ICD10:M340):50症例

対照健常者:48名

(対象の一部はE-GEAD-397と重複)

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞24画分(Naive_CD4, Mem_CD4, Fr._I_nTreg, Fr._II_eTreg, Fr._III_T, Th1, Th2, Th17, Tfh, NK, Naive_CD8, EM_CD8, CM_CD8, TEMRA_CD8, Naive_B, USM_B, SM_B, DN_B, Plasmablast, CL_Mono, Int_Mono, NC_Mono, mDC, pDC)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp
マッピング方法 STAR(hg38)
フィルタリング(QC)方法 アダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。50塩基未満のリード、及び低クオリティーの塩基を多く含むリード(20%を超える塩基で、Phred quality score < 20)を除去した。ユニークマップ率が80%未満、もしくはユニークマップリード数が5.00 x 106未満のサンプルは除去した。同一サブセット内の遺伝子発現の相関係数(Di)を計算した。Diが0.9未満のサンプルは除去した。
遺伝子数 26353
Japanese Genotype-phenotype Archive Dataset ID JGAD000406
総データ量 172.8 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000220(RNA-seq:全身性エリテマトーデス)

対象

全身性エリテマトーデス(ICD10:M329):107症例

対照健常者:92名

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞19画分(Naive_CD4, Mem_CD4, Fr._II_eTreg, Th1, Th2, Th17, Tfh, NK, Naive_CD8, Mem_CD8, Naive_B, USM_B, SM_B, DN_B, Plasmablast, CD16n_Mono, CD16p_Mono, mDC, pDC)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) Smart-Seq v2(test cohort)、SMART-seq v4 Ultra Low Input RNA Kit(validation cohort)
断片化の方法 Smart-Seq v2(test cohort)、SMART-seq v4 Ultra Low Input RNA Kit(validation cohort)
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー)

76 bp(test cohort)

100 bp(validation cohort)

マッピング方法 STAR(GRCh38)
フィルタリング(QC)方法 FASTQファイルはSTAR(v2.5)によってGRCh38; GenBank assembly GCA_000001405.18にマッピングした。カウントテーブルの作成にはHTSeq-count(v0.6.1)を用いた。fastx-toolkit(v0.0.14)を用いてアダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。ミトコンドリア遺伝子発現は細胞ストレスを反映するため、ミトコンドリア遺伝子発現 <8%をカットオフに追加した。同一サブセット内の遺伝子発現の相関係数(r)を計算し、r2 < 0.8のサンプルは外れ値として除外した。
遺伝子数

26354(test cohort)

26485(validation cohort)

Japanese Genotype-phenotype Archive Dataset ID

JGAD000371

JGAD000372

総データ量 250 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000220(ATAC-seq)

対象

全身性エリテマトーデス(ICD10:M329):8症例

対照健常者:8名

規模 ATAC-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞15画分(Naive_B, SM_B, USM_B, DN_B, Plasmablast, Th1, Th2, Th17, Tfh, Naive_CD4, Mem_CD4, Naive_CD8, CD16p_Mono, CD16n_Mono, NK)より抽出したDNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) Fast-ATAC-seq protocol
断片化の方法 Fast-ATAC-seq protocol
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 102 bp
Japanese Genotype-phenotype Archive Dataset ID JGAD000373
総データ量 27 GB(tdf、bed)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000486(RNA-seq:全身性エリテマトーデス)

対象

全身性エリテマトーデス(ICD10:M329):136症例

(22症例については、経時的な解析を実施、計159検体)

対照健常者:89名

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500、NovaSeq 6000]
ライブラリソース 末梢血免疫細胞27画分(Naive_CD4, Mem_CD4, Th1, Th2, Th17, Tfh, Fr._I_nTreg, Fr._II_eTreg, Fr._III_T, Naive_CD8, EM_CD8, CM_CD8, TEMRA_CD8, NK, Naive_B, USM_B, SM_B, DN_B, Plasmablast, CL_Mono (or CD16n_Mono), CD16p_Mono, Int_Mono, NC_Mono, mDC, pDC, Neu, LDG)より抽出したRNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp (HiSeq 2500)、150 bp (NovaSeq 6000)
マッピング方法 STAR(GRCh38)
フィルタリング(QC)方法 アダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。低クオリティーの塩基が20%を超えるリードを除去した。ユニークマップ率が80%未満、もしくはユニークマップリード数が5.00 x 106未満のサンプルは除去した。同一サブセット内の遺伝子発現の相関係数(Di)が0.9未満のサンプルは除去した。同一個人のwhole genome sequencingのhetero locusに対するgenotypeの一致率が0.9未満のサンプルは除去した。
遺伝子数 26353
Japanese Genotype-phenotype Archive Dataset ID JGAD000603
総データ量 478.3 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000598(RNA-seq:関節リウマチ)

対象

関節リウマチ(ICD10:M0690):50症例

(15症例については、経時的な解析を実施)

対照健常者:39名

(対象の一部はE-GEAD-397と重複)

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞18画分(CD16p_Mono, CL_Mono, DN_B, Fr_II_eTreg, mDC, Mem_CD4, Naive_B, Naive_CD4, Neu, NK, pDC, Plasmablast, SM_B, Tfh, Th1, Th17, Th2, USM_B)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp
マッピング方法 STAR(hg38)
フィルタリング(QC)方法 アダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。50塩基未満のリード、及び低クオリティーの塩基を多く含むリード(20%を超える塩基で、Phred quality score < 20)を除去した。ユニークマップ率が80%未満、もしくはユニークマップリード数が5.00 x 106未満のサンプルは除去した。同一サブセット内の遺伝子発現の相関係数(Di)を計算した。Diが平均-2標準偏差未満のサンプルは除去した。
遺伝子数 26353
Japanese Genotype-phenotype Archive Dataset ID JGAD000727
総データ量 113.9 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000485(B細胞受容体レパトア解析)

対象

全身性エリテマトーデス(ICD10:M329):136症例(159検体、23検体は経時サンプル)

全身性強皮症(ICD10:M340):90症例(90検体)

筋炎(ICD10:M339、M332):85症例(85検体)

混合性結合組織病(ICD10:M351):21症例(21検体)

シェーグレン症候群(ICD10:M350):18症例(18検体)

関節リウマチ(ICD10:M0690):27症例(27検体)

ベーチェット病(ICD10:M352):23症例(23検体)

成人スチル病(ICD10:M0610):18症例(18検体)

ANCA関連血管炎(ICD10:M318):24症例(24検体)

高安動脈炎(ICD10:M314):16症例(16検体)

対照健常者:13名(13検体)

規模 RNA-seq / B細胞受容体レパトア解析
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血B細胞5画分(Naive_B, USM_B, SM_B, DN_B, Plasmablast)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-Seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-Seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp
マッピング方法 MiXCR
フィルタリング(QC)方法 アダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。50塩基未満のリード、及び低クオリティーの塩基を多く含むリード(20%を超える塩基で、Phred quality score < 20)を除去した。 ユニークなCDR-H3配列を 500より多く有するサンプルのみ解析に使用した。
B細胞受容体レパトア解析方法(ソフトウェア) MiXCR
Japanese Genotype-phenotype Archive Dataset ID JGAD000602
総データ量 12.9 GB(B細胞受容体レパトアデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000626(RNA-seq)

対象

関節リウマチ(ICD10:M0690):19症例

全身性エリテマトーデス(ICD10:M329):62症例

特発性炎症性筋疾患(ICD10:M339、M332):40症例

対照健常者:64名

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500]
ライブラリソース 末梢血免疫細胞9画分(Naive_CD4, Th1, Th2, Th17, Tfh, Fr._I_nTreg, Fr._II_eTreg, Fr._III_T, ThA)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp
マッピング方法 STAR(GRCh38)
フィルタリング(QC)方法 cutadapt (v1.16)およびfastx-toolkit(v0.0.14)を用いてアダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。低クオリティーの塩基が20%を超えるリードを除去した。STAR(v2.5.3)のtwo-pass modeを用いてGRCh38へのマッピングを行った(Gencode version 27アノテーション)。ユニークマップ率が90%未満(Plasmablastでは70%未満、他のB細胞サブセットでは<85%未満)、もしくはユニークマップリード数が6.00 x 10^6未満のサンプルは除去した。遺伝子発現はHTSeq(v 0.11.2)で定量した。各細胞分画で低発現遺伝子(90%を超えるサンプルで10リード未満)を除外し、edgeRのTMMノーマライゼーションを実施、ログcount per million(CPM)変換、ComBatによるバッチ除去を行った上で、同一サブセット内の遺伝子発現のSpearman相関係数を計算し、平均値が0.9未満のサンプルを除外した。
遺伝子数 26353
Japanese Genotype-phenotype Archive Dataset ID JGAD000755
総データ量 36.3 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000627(RNA-seq)

対象

全身性エリテマトーデス(ICD10:M329):136症例

対照健常者:89名

規模 RNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [HiSeq 2500、NovaSeq 6000]
ライブラリソース 末梢血免疫細胞27画分(Naive_CD4, Th1, Th2, Th17, Tfh, Fr._I_nTreg, Fr._II_eTreg, Fr._III_T, ThA, Naive_CD8, Naive_B, SM_B, USM_B, DN_B, Plasmablast, NK, CD16p_Mono, CL_Mono, Neu, mDC, pDC, TEMRA_CD8, CM_CD8, EM_CD8, NC_Mono, Int_Mono, LDG)より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) SMART-seq v4 Ultra Low Input RNA Kit
断片化の方法 SMART-seq v4 Ultra Low Input RNA Kit
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 100 bp (HiSeq 2500)、150 bp (NovaSeq 6000)
マッピング方法 STAR(GRCh38)
フィルタリング(QC)方法 アダプター配列及び3’の低クオリティー塩基をトリムした(Phred quality score < 20)。低クオリティーの塩基が20%を超えるリードを除去した。ユニークマップ率が80%未満、もしくはユニークマップリード数が5.00 x 106未満のサンプルは除去した。同一サブセット内の遺伝子発現の相関係数(Di)が0.9未満のサンプルは除去した。同一個人のwhole genome sequencingのhetero locusに対するgenotypeの一致率が0.9未満のサンプルは除去した。
遺伝子数 26353
Japanese Genotype-phenotype Archive Dataset ID JGAD000756
総データ量 129.4 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

JGAS000648(scRNA-seq)

対象

炎症性筋疾患(ICD10:M339、M332):6症例

     筋生検検体:3症例3検体、気管支洗浄液検体:3症例3検体

対照健常者:1名

     末梢血:1検体

規模 scRNA-seq
対象領域(Target Captureの場合) -
Platform Illumina [NovaSeq 6000]
ライブラリソース 筋組織・気管支洗浄液より分取したCD45陽性細胞および末梢血より分取したCD4T細胞より抽出したtotal RNA
検体情報(購入の場合) -
ライブラリ作製方法(キット名) Chromium Next GEM Single Cell 5’ Reagent Kit v2 Dual Index
断片化の方法 Chromium Next GEM Single Cell 5’ Reagent Kit v2 Dual Index
ライブラリ構築方法 Paired-end
リード長(除:バーコード、アダプタ、プライマー、リンカー) 150 bp
マッピング方法 Cell Ranger (GRCh38)
フィルタリング(QC)方法 検出遺伝子数が、末梢血検体は500~3,500、筋生検検体は1,500~5,500、気管支洗浄液検体は1,000~5,500の細胞、およびミトコンドリア遺伝子のリードはいずれも5%未満である細胞を解析に使用した。
遺伝子数

末梢血: 17,441 genes

筋生検: 18,184 genes

気管支洗浄液: 20,314 genes

Japanese Genotype-phenotype Archive Dataset ID JGAD000778
総データ量 40.2 MB(各遺伝子にマップされたリードのカウントデータ、txt)
コメント(利用にあたっての制限事項) NBDC policy

 

提供者情報

研究代表者: 藤尾 圭志

所 属 機 関: 東京大学大学院 医学系研究科 アレルギー・リウマチ学

プロジェクト/研究グループ名: 免疫疾患の免疫細胞マルチオミクス解析

URL: http://ryumachi.umin.jp/

科研費/助成金(Research Project Number):

科研費・助成金名タイトル研究課題番号
中外製薬株式会社との共同研究費 - -
日本医療研究開発機構(AMED) 難治性疾患実用化研究事業 ゲノムおよび遺伝子発現情報の統合的解析に基づく全身性エリテマトーデスの治療標的の同定とその制御法の開発 JP17ek0109103
日本医療研究開発機構(AMED) ゲノム医療実現バイオバンク利活用プログラム 免疫担当細胞eQTLデータを用いた免疫介在性疾患ゲノム情報からの層別化および予後予測モデルの構築 JP21tm0424221
日本医療研究開発機構(AMED) 健康・医療分野におけるムーンショット型研究開発等事業 病気につながる血管周囲の微小炎症を標的とする量子技術、ニューロモデュレーション医療による未病時治療法の開発 JP21zf0127004
日本医療研究開発機構(AMED)免疫アレルギー疾患実用化研究事業 自己免疫性疾患の臓器病変局所におけるシングルセルRNAシークエンスを用いたマルチオミックス解析による病態解明基盤の構築 JP22ek0410074
国立研究開発法人日本医療研究開発機構(AMED) 革新的先端研究開発支援事業 自己免疫疾患における細胞傷害応答と適応免疫応答双方の免疫記憶に関連するT細胞サブセットの解明 JP23gm1810005
科学研究費助成事業 基盤研究(B) マルチオームシングルセル解析を用いた自己免疫疾患における加齢関連T細胞の機能解析 JP22H03110

 

関連論文

タイトルDOIデータID
1 Integrated bulk and single-cell RNA-sequencing identified disease-relevant monocytes and a gene network module underlying systemic sclerosis doi: 10.1016/j.jaut.2020.102547

JGAD000309

E-GEAD-344

2 Identifying the most influential gene expression profile in distinguishing ANCA-associated vasculitis from healthy controls doi: 10.1016/j.jaut.2021.102617 JGAD000310
3 Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases doi: 10.1016/j.cell.2021.03.056

E-GEAD-397

E-GEAD-398

E-GEAD-420

4 Dysregulation of the gene signature of effector regulatory T cells in the early phase of systemic sclerosis doi: 10.1093/rheumatology/keac031 JGAD000406
5 Immune cell multiomics analysis reveals contribution of oxidative phosphorylation to B-cell functions and organ damage of lupus doi: 10.1136/annrheumdis-2021-221464

JGAD000371

JGAD000372

JGAD000373

6 Distinct transcriptome architectures underlying lupus establishment and exacerbation doi: 10.1016/j.cell.2022.07.021 JGAD000603
7 Immunomics analysis of rheumatoid arthritis identified precursor dendritic cells as a key cell subset of treatment resistance doi: 10.1136/ard-2022-223645

JGAD000727

JGAD000731

8 Multimodal repertoire analysis unveils B cell biology in immune-mediated diseases doi: 10.1136/ard-2023-224421 JGAD000602
9 Age-associated CD4+ T cells with B cell-promoting functions are regulated by ZEB2 in autoimmunity doi: 10.1126/sciimmunol.adk1643

JGAD000755

JGAD000756

JGAD000778

 

制限公開データの利用者一覧

研究代表者所属機関国・州名研究題目利用データID利用期間
Jiucun Wang Department of anthropology and human genetics, Fudan University 中華人民共和国 The role and mechanism study of glycosyltransferase-B3GNT2 in regulating macrophages of Ankylosing Spondylitis JGAD000309,
JGAD000603
2023/03/07-2024/03/01
Jacob D Jaffe Odyssey Therapeutics アメリカ合衆国 Investigation of immune cell transition states in autoimmune disease JGAD000309,
JGAD000371,
JGAD000372,
JGAD000603
2023/03/19-2023/10/02
武田 吉人 大阪大学大学院医学系研究科 呼吸器・免疫内科学 日本 呼吸器系・免疫疾患における細胞外小胞に着目した病態解明および診断法の開発 JGAD000310,
JGAD000371,
JGAD000372,
JGAD000373
2023/04/10-2025/03/31
Ana Rita Grosso Universidade Nova de Lisboa - NOVA School of Science and Technology ポルトガル Assessing transcriptional dyregulation of repetitive elements and monoallelic-expressed genes in lupus JGAD000603 2023/04/25-2026/11/30
Yong-Fei WANG The Chinese University of Hong Kong 香港 Investigating the Molecular Mechanisms of Systemic Lupus Erythematosus Using Functional Genomics Data JGAD000371,
JGAD000372,
JGAD000373,
JGAD000603
2023/06/21-2025/05/22
Timothy Vyse King's College London
Molecular and Medical genetics
イギリス Sequencing Based Genetic Analysis of Systemic Lupus JGAD000603 2023/08/25-2025/07/28
Ana Rita Grosso Universidade Nova de Lisboa - NOVA School of Science and Technology ポルトガル Assessing transcriptional dyregulation of repetitive elements and monoallelic-expressed genes in lupus JGAD000603 2023/04/25-2026/11/30
久保田 晋平 北海道大学 遺伝子病制御研究所 日本 関節リウマチおよび全身性エリテマトーデスにおける遺伝子発現および炎症増幅回路活性化因子の解析 JGAD000309,
JGAD000310,
JGAD000371,
JGAD000372,
JGAD000373,
JGAD000406,
JGAD000602,
JGAD000603,
JGAD000727
2024/02/26-2025/03/31
鈴木 宏尚 科研製薬株式会社 新薬創生センター 薬理部 日本 SLE患者検体の遺伝子発現データを用いたバイオインフォマティクス解析 JGAD000603 2024/05/09-2025/03/31