DOI: 10.5281/zenodo.14825318
title: Supplemental data for "The genome of a sea spider corroborates a shared Hox cluster motif in arthropods with reduced posterior tagma"

description: # List of files on Zenodo

 

Data accompanying the manuscript have been uploaded on Zenodo (10.5281/zenodo.14185694). Here we

present a brief description of each file and put them in meaningful groups.

 

## referenced supplement

 

Supplementary material from the manuscript.

 

- suppl-fig1-alignment.pdf

- suppl-fig2-genomescope.png

- suppl-fig3-abdA_araneae.png

- suppl-fig4-genomic_context.png





- suppl-fig5-oxford.png

- suppl-fig6-correlation.png

- suppl-fig7-heatmap.png

 

 



- suppl-table1-data_overview.tsv

- suppl-table2-Plit_COI-alignment_distances.pdf

- suppl-table3-isoseq.tsv

- suppl-table4-genome_progress.tsv

- suppl-table5-abdA.tsv

- suppl-table6-r2_g3735-Alignment-HitTable.csv

- suppl-table7-chelicerate_genomes.tsv

- suppl-table8-arthropod_repeat_content.tsv

- suppl-table9-chelicerate_repeat_content.tsv

- suppl-table10-Pycnognonum_microRNAs.xlsx

 

- suppl-file1-araneae_abdA.fasta

- suppl-file2-arthropod_genomes.tsv





- suppl-file3-hox_tree.pdf

- suppl-file4-hro_tree.pdf

- suppl-file5-irx_tree.pdf

- suppl-file6-sine_tree.pdf

- suppl-file7-nk_tree.pdf



 

- unref-table-mapping_rates.tsv

 

## figures

 

- figs.zip: archive of raw and processed figures for the manuscript in full resolution

 

## analysis

 

### genomic context

 

The broader arthropod/chelicerate context for the _P. litorale_ genome assembly. Files used/generated by the [genomic_context](https://github.com/galicae/plit-genome/blob/main/07-analysis/genomic_context.ipynb) notebook.

 

- araneae.tsv: repeat makeup of published chelicerate assemblies

- arthropoda.tsv: genome assembly statistics for arthropod genomes. From NCBI Genomes.

- modern_taxids.txt: list of taxonomic IDs for species; made to be submitted to NCBI Taxonomy.

- tax_report.txt: the full taxonomic report for each query species. Contains tax IDs for the entire lineage.

- total_repeats.tsv: total repeat content of published chelicerate genomes.

 

### Hox genes

 

Files concerning the Hox gene cluster analysis.

 

- Plit_HoxGenes.fasta: putative Hox gene sequences for P. litorale

- putative_hox-Alignment-HitTable.csv:

- araneae_abdA.fasta: arachnid abdA sequences from NCBI

- araneae_abdA.m8: alignment of arachnid abdA sequences against the draft genome

- plit.m8: MMseqs2 alignment results for the putative Hox sequences against the draft genome.

- r2_g3735-Alignment-HitTable.csv: NCBI BLAST results (nr) for gene r2_g3735

 

## processed (intermediate) data

 

### 00-kmer-jellyfish.zip

 

k-mer spectra analysis with GenomeScope and GenomeScope2.0

 

### 00-seq-qc.zip

 

quality control output for raw sequencing data (e.g. FastQC output)

 

### 01-assembly

 

- assembly_graph.gfa: Flye output

- assembly_graph.gv: Flye output

- assembly_info.txt: Flye output

- assembly.fasta: Flye output

- backmap.hifi.sort.bam.cov-hist.pdf: coverage histogram of the back-mapped PacBio data

- backmap.ont.sort.bam.cov-hist.pdf: coverage histogram of the back-mapped ONT data

- BUSCO.arthropoda_odb10.txt: BUSCO completeness report (arthropoda_odb10)

- BUSCO.metazoa_odb10.txt: BUSCO completeness report (metazoa_odb10)

- flye.log: Flye assembler log

- quast_report.pdf: assembly QC

 

### 02-scaffold

 

`yahs` output files:

 

- asm_hic.sorted.bam

- flye-yahs.fa

- yahs.out_scaffolds_final.agp

- yahs.out_scaffolds_final.fa.hic

- yahs.out_scaffolds_final.fa.assembly

 

Juicebox (manual curating) results:

 

- yahs.out_scaffolds_final.fa.review.assembly

- 02-flye-yahs-juicebox.fa

 

GAP `sort_scaffolds` pipeline outputs

 

- 03-flye-yahs-juicebox-merge.fasta

- plit_q_0_50000_0.5FracBest_unseen_scaffolds.txt

- plit_q_0_50000_0.5FracBest_insertion_stats.tsv

- plit_q_0_50000_0.5FracBest_appended_scaffolds.tsv

- plit_q_0_50000_0.5FracBest_inserted_scaffolds.tsv

 

### 03-contamination

 

Refer to the [contamination analysis](https://github.com/galicae/plit-genome/blob/main/04-contam/README.md) for details.

 

Decontaminating the draft genome from non-metazoan scaffolds:

 

- plit_q_0_50000_0.5FracBest_output_filtered.fasta: input draft genome

- contam_tax.m8: Alignment results of UniRef90 against draft genome (MMseqs2)

- scaffolds_taxonomic_distribution.tsv: summary of contam_tax.m8; number of genes from each taxonomic level per scaffold.

- scaffolds_taxonomic_distribution_collapsed_vir.tsv: scaffolds with predominantly viral hits

- scaffolds_taxonomic_distribution_suspect.tsv: scaffolds whose genes are <90% metazoan

 

Checking for widespread _Metridium_ contamination:

 

- primary_mq30.txt: list of high-quality mapping reads (presumptive "metridial")

- metridium_scaffolds.txt_summary: no. of presumptive _Metridium_ reads per draft scaffold

- metridium_scaffolds.txt: filtered SAM file with all high-quality "_Metridium_" hits on draft scaffolds

- metridium_contigs.sam_summary: no. of presumptive _Metridium_ reads per Flye contig

 

### 04-annotation

 

Repeat analysis with RepeatModeler/RepeatMasker:

 

- draft.fasta.tbl: output of RepeatModeler in tabular form

- pb.sam.flagstats: summary of mapping the repeat families to the PacBio data.

- pycno-families.fa: output of RepeatModeler - the sequences of the _P. litorale_ repeat families

- draft.fasta.out.gff: output of RepeatModeler - repeat locations on the draft genome

 

Protein coding gene annotation:

 

- annot-01-isoseq.gff: GFF file with the gene models proposed using Iso-seq isoforms

- annot-01-braker.gff: GFF file with the gene models proposed from round 1 of BRAKER3 using developmental transcriptomes

- annot-02-braker.gff3: GFF file with the gene models proposed from round 2 of BRAKER3, using developmental transcriptome reads that weren't used in round 1

- annot-03-denovo.gff3: GFF file with the gene models proposed from the de novo transcriptomes

- deep_denovo_assemblies.zip: the de-novo assembled transcriptomes from the deeply sequenced developmental time points. Also available on ENA.

 

tRNAscan output

 

- trnascan.bed

- trnascan.out

- trnascan.fasta

- trnascan.stats

 

MirMachine output

 

- Pli_september.PRE.gff: MirMachine output with permissive threshold

- Pli_september.PRE-1.gff: MirMachine output with strict threshold

- Pli_september.PRE.fasta: predicted miRNA sequences

 

## Results

 

- draft_softmasked.fasta: draft genome with repetitive regions softmasked

- draft.fasta: draft genome fasta

- hox.gff3: the position of the Hox genes in GFF3 form.

- merged_sorted.gff3: protein-coding gene models from all rounds of annotation

- transcripts.fa: TransDecoder-extracted putative transcripts

- transcripts.fa.transdecoder.pep: TransDecoder predicted peptides

- out.emapper.annotations: EggNOG-mapper functional annotation for the predicted peptides

- out.emapper.best.annotations: filtered EggNOG-mapper annotation; best hit per gene kept

- miRNA.fasta: FASTA sequences of predicted miRNAs

- miRNA.lenient.gff: GFF of miRNA positions (permissive MirMachine cutoff)

- miRNA.strict.gff: GFF of miRNA positions (strict MirMachine cutoff)