ウイルスゲノム（SARS-CoV-2）をTerraで解析する

この記事は[バイオインフォマティクス Advent Calendar 2020](https://qiita.com/advent-calendar/2020/bioinfo)の19日目の記事です

今年はCOVID19のパンデミックによりSARS-CoV-2の解析が注目されましたが、その解析はどのようなパイプラインになっているのでしょうか。Terraという解析プラットフォームに[COVID19のWorkSpace](https://support.terra.bio/hc/en-us/articles/360041068771--COVID-19-workspaces-data-and-tools-in-Terra)があったので、Terraを使ってみながら解析の流れを紹介したいと思います。

# SARS-CoV-2ゲノムの構造

SARS-CoV-2ゲノムは３０ｋ塩基長の一本鎖RNAです。
![](https://upload.wikimedia.org/wikipedia/commons/a/a8/SARS-CoV-2_genome.svg)

すべてのβコロナウイルスゲノムに共通するオープンリーディングフレーム（ORF）と、スパイク（S）、エンベロープ（E）、マトリクス（M）、ヌクレオカプシド（N）タンパク質をコードしています。

# Terra

[Terra](https://terra.bio/)はGATKでおなじみの[Broad Institute](https://www.broadinstitute.org/)とAlphabet傘下の[Verily](https://verily.com/)が共同で開発したバイオメディカル研究のためのデータ解析プラットフォームです。

![](https://terra.bio/wp-content/uploads/2020/12/Ecosystem_Graphic.png)

GCPのアカウントを用意すれば（従量課金ですが）、ワークフローを用いたデータ解析をGUIで実行することができ、データや結果の共有を簡単に行うことができます。
ワークフローはWDLで記述され、[Dockstore](https://dockstore.org/)と[Broad Method Repository](https://support.terra.bio/hc/en-us/articles/360025674392-Finding-the-workflow-method-you-need-and-its-JSON-in-the-Methods-Repository)から追加することができます。
また、1000 GenomesやENCODE Project、その他[様々なデータ](https://app.terra.bio/#library/datasets)を参照することができるようになっています。

# 解析ワークフロー

ウイルスゲノムを解析する場合にはDe novoアセンブリを用いたワークフローもありますが、ここではSARS-CoV-2のみを対象しているのでより効率的なリファレンス配列を用いた解析になっています。
![](http://drive.google.com/uc?export=view&id=1L5YgnnksOAFkdccZZFTIok1dbXYJB6YT)

1. SRAまたはFASTQファイルからuBAMファイルに変換する。
DATAタブであらかじめサンプルテーブルを作成しておく必要があります。SRAの場合はIDを、ファイルの場合はGoogle Storageへのファイルパスを指定します。

2. uBAMファイルからリファレンスゲノムを用いたアセンブルを行う。
1で作成されたuBAMファイルを指定して、リファレンスゲノムへのアライメントを行います。

3. 系統樹を作成する。
[Augur](https://docs.nextstrain.org/projects/augur/en/stable/)というツールでアセンブルされたFASTAから系統樹を作成します。

各ワークフローは以下のようなWDL形式で記述されていますが、Terra上では実行したいワークフローを選択するだけです。

```wdl
version 1.0

import "../tasks/tasks_ncbi_tools.wdl" as ncbi_tools

workflow multi_Fetch_SRA_to_BAM {
    input {
        Array[String]  SRR_accessions
    }
    scatter(SRA_ID in SRR_accessions) {
        call ncbi_tools.Fetch_SRA_to_BAM {
            input:
                SRA_ID = SRA_ID
        }
    }
    output {
        Array[File]   reads_ubam = Fetch_SRA_to_BAM.reads_ubam
        Array[String] sequencing_center = Fetch_SRA_to_BAM.sequencing_center
        Array[String] sequencing_platform = Fetch_SRA_to_BAM.sequencing_platform
        Array[String] sequencing_platform_model = Fetch_SRA_to_BAM.sequencing_platform_model
        Array[String] biosample_accession = Fetch_SRA_to_BAM.biosample_accession
        Array[String] library_id = Fetch_SRA_to_BAM.library_id
        Array[String] run_date = Fetch_SRA_to_BAM.run_date
        Array[String] sample_collection_date = Fetch_SRA_to_BAM.sample_collection_date
        Array[String] sample_collected_by = Fetch_SRA_to_BAM.sample_collected_by
        Array[String] sample_strain = Fetch_SRA_to_BAM.sample_strain
        Array[String] sample_geo_loc = Fetch_SRA_to_BAM.sample_geo_loc
        Array[File]   sra_metadata = Fetch_SRA_to_BAM.sra_metadata
    }
}
```



選択したWorkflowで入力ファイルを指定して`RUN ANALYSIS`を押すと解析が実行されます。

**fetch_sra_to_bamの例**
![run_analysis.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/99878/6dacf41d-ef41-81ac-0c76-3e088edcf974.png)


## 系統樹の描画

系統樹を描画するには[Nextstrain](https://docs.nextstrain.org/en/latest/tutorials/quickstart.html)を独自に立ち上げてもいいですが、[auspice.us](http://auspice-us.herokuapp.com/)というサイトが公開されているのでそこに3で作成された系統樹のJSONファイルをDrag&Dropすることで以下のように描画することができます。
![Screen Shot 2020-12-20 at 23.55.00.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/99878/e35004ea-7134-cbbb-7717-7809e98f08b6.png)

**参考**

- [pathogen-genomic-surveillance/COVID-19](https://app.terra.bio/#workspaces/pathogen-genomic-surveillance/COVID-19)