Repository navigation
Expand file tree
/
Copy pathdata-availability.qmd
More file actions
134 lines (85 loc) · 7.93 KB
/
Copy pathdata-availability.qmd
File metadata and controls
134 lines (85 loc) · 7.93 KB
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
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
---
title: "Data Availability"
description: |
On this page, you can find information on obtaining sequence data, data products, and processing scripts. Code is embedded in indivudual workflows on the website or can be accessed through the project GitHub repo.
---
## Archived Sequence Data
Raw **16S rRNA** data for each sample (before removing primers) is available on figshare at [doi:10.25573/data.14686665](https://doi.org/10.25573/data.14686665). Raw **ITS** data for each sample (before removing primers) is available on figshare at [doi:10.25573/data.14686755](https://doi.org/10.25573/data.14686755).
All trimmed sequence data (primers removed) is deposited at the European Nucleotide Archive (ENA) under the study accession number [PRJEB45074 (ERP129199)](https://www.ebi.ac.uk/ena/browser/view/PRJEB45074). The trimmed **16S rRNA** data (primers removed) are deposited under sample accession numbers [ERS6485270-ERS6485284](https://www.ebi.ac.uk/ena/browser/view/ERS6485270-ERS6485284). The trimmed **ITS** data (primers removed) are deposited under sample accession numbers [ERS6485285-ERS6485299](https://www.ebi.ac.uk/ena/browser/view/ERS6485285-ERS6485299).
## Pipeline Data
Data for each individual pipeline are available through the Smithsonian figshare under a single collection at [doi:10.25573/data.c.5667571](https://doi.org/10.25573/data.c.5667571). In addition, data from each pipeline are available for download from figshare using the links at the bottom of each workflow page (where applicable).
</br>
:::{.callout-note collapse="true" icon=false}
## **Show** descriptions and DOI's for each workflow
Scott, Jarrod (2022): (high-temp/16s) No 0. Raw 16S rRNA fastq data files for the high-temp SWELTR study. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.14686665.v1](https://doi.org/10.25573/data.14686665.v1)
Scott, Jarrod (2022): (high-temp/ITS) No 0. Raw ITS fastq data files for the high-temp SWELTR study. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.14686755.v1](https://doi.org/10.25573/data.14686755.v1)
Scott, Jarrod (2022): (high-temp) No 1. [DADA2 Workflow (16S rRNA/ITS)](https://sweltr.github.io/high-temp/dada2.html) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.14687184.v1](https://doi.org/10.25573/data.14687184.v1)
Scott, Jarrod (2022): (high-temp) No 2. [Data Preparation](https://sweltr.github.io/high-temp/data-prep.html) (16S rRNA/ITS) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.14690739.v1](https://doi.org/10.25573/data.14690739.v1)
Scott, Jarrod (2022): (high-temp) No 3. [Filtering](https://sweltr.github.io/high-temp/filtering.html): (16S rRNA/ITS) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.14701440.v1](https://doi.org/10.25573/data.14701440.v1)
Scott, Jarrod (2022): (high-temp) No 4. [Taxonomic diversity](https://sweltr.github.io/high-temp/taxa.html): (16S rRNA/ITS) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.16826737.v1](https://doi.org/10.25573/data.16826737.v1)
Scott, Jarrod (2022): (high-temp) No 5. [Aplha diversity](https://sweltr.github.io/high-temp/alpha.html) (16S rRNA/ITS) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.16826779.v1](https://doi.org/10.25573/data.16826779.v1)
Scott, Jarrod (2022): (high-temp) No 6. [Beta diversity](https://sweltr.github.io/high-temp/beta.html) (16S rRNA/ITS) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.16828063.v1](https://doi.org/10.25573/data.16828063.v1)
Scott, Jarrod (2022): (high-temp) No 7. [Differentially Abundant ASVs & Taxa](https://sweltr.github.io/high-temp/da.html) (16S rRNA/ITS) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.16828249.v1](https://doi.org/10.25573/data.16828249.v1)
Scott, Jarrod (2022): (high-temp) No 8. [Metadata Analysis](https://sweltr.github.io/high-temp/metadata.html) (16S rRNA/ITS) Output. Smithsonian Tropical Research Institute. Dataset. [https://doi.org/10.25573/data.16828294.v1](https://doi.org/10.25573/data.16828294.v1)
:::
## Submitting Sequence Data
We submitted out data to the [European Nucleotide Archive (ENA)](https://www.ebi.ac.uk/ena). The ENA does not like RAW data and prefers to have primers removed. So we submitted the trimmed Fastq files to the ENA. You can find these data under the study accession number **PRJEB45074 (ERP129199)**. The RAW files on our figshare site (see above).
To submit to the ENA you need two data tables (plus your sequence data). One file describes the samples and the other file describes the sequencing data.
You can download our submission data tables here:
* [Sample submission form](include/data-availability/ENA_Sample_SUBMISSION_SWELTR-high-temp.tsv)
* [FASTq submission form](include/data-availability/ENA_PAIRED_FASTQ_SUBMISSION_SWELTR-high-temp.tsv)
Note, these forms are study specific, so please use these as guides only.
### Instructions for Submitting to the ENA
Even though I have done this dozens of times, the process of submitting sequence data to read archives still baffles me.
> Note: I submit data with primers, barcodes, etc **removed**
#### Register Project & Upload Sample Data
1) go to https://www.ebi.ac.uk/ena/submit and select **Submit to ENA Interactively**.
2) Login or Register.
3) Go to **New Submission** tab and, if this is a new project, select **Register study (project)**.
4) Hit Next
5) Enter details and hit **Submit**.
6) Next, [Select Checklist](https://www.ebi.ac.uk/ena/submit/checklists). This will be specific to the type of samples you have and basically will create a template so you can add your sample metadata. For this study I chose **GSC MIxS soil**, checklist accession number [ERC000022](https://www.ebi.ac.uk/ena/data/view/ERC000022)
7) Next
8) Now go through and select/deselect fields as needed. Note, some fields are mandatory.
9) Once finished, hit **Next** to fill in any details that apply to *All* samples and the download the template. Alternatively, you can download the template and fill in the data by hand.
10) Upload the sample sheet.
11) Once everything looks good and uploaded, click the **New Submission** tab.
#### Upload Sequence Data
12) Hit **Skip** and then select **Two Fastq files (Paired)**, Download the template and fill in the details.
13) Next, make sure all the trimmed fastq files are gzipped **.gz** (these are what you submit).
14) Navigate to the directory with the trimmed, compressed fastq files and run:
```
md5sum *.gz
```
16) Add the checksums and file names to the fastq submission form. You can read more about [Preparing A File For Upload here](https://ena-docs.readthedocs.io/en/latest/submit/fileprep/preparation.html).
**STOP**
Before uploading the sheet, you first must upload the fastq files. Documentation for Uploading Files To ENA can be found [here](https://ena-docs.readthedocs.io/en/latest/submit/fileprep/upload.html).
17) In the same directory, run:
```
lftp webin2.ebi.ac.uk -u Webin-XXXXX
```
where `Webin-XXXXX` is your user name. Enter your password
at the prompt run:
```
mput *.gz
# when finished
bye
```
The files should begin uploading. Depending on internet speed and/or file sizes/numbers, this could take a while.
18) Once the upload is finished, upload and submit the fastq submission form. If the `sample alias` field dows not autofill you may need to upload the sample form again.
```{r}
#| message: false
#| results: hide
#| eval: true
#| echo: false
remove(list = ls())
### COmmon formatting scripts
### NOTE: captioner.R must be read BEFORE captions_XYZ.R
source(file.path("assets", "functions.R"))
```
#### Source Code {.appendix}
The source code for this page can be accessed on GitHub `r fa(name = "github")` by [clicking this link](`r source_code()`).
#### Last updated on {.appendix}
```{r,echo=FALSE, eval=TRUE}
Sys.time()
```