# Siberian sturgeon multi-tissue reference transcriptome database

Christophe C. Kropp, Cédric Cabau, Gonzalo Greif, André Lasalle, Santiago Di Landro, Denise Vizziano-Cantonnet

► To cite this version:

Christophe C. Kropp, Cédric Cabau, Gonzalo Greif, André Lasalle, Santiago Di Landro, et al. Siberian sturgeon multi-tissue reference transcriptome database. Database - The journal of Biological Databases and Curation, 2020, 2020, pp.1-8. 10.1093/database/baaa082. hal-03157672

HAL Id: hal-03157672

https://hal.inrae.fr/hal-03157672

Submitted on 5 Mar 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Original article

# Siberian sturgeon multi-tissue reference transcriptome database

Christophe Klopp¹, Cédric Cabau², Gonzalo Greif³, André Lasalle⁴,
Santiago Di Landro⁴ and Denise Vizziano-Cantonnet⁴,†

¹SIGENAE, Genotoul Bioinfo, MIAT UR875, INRAe, Chemin de Borde-Rouge - Auzeville BP 52627, 31326 CASTANET-TOLOSAN CEDEX, France, ²SIGENAE, GenPhySE, Université de Toulouse, INRAe, ENVT, Chemin de Borde-Rouge - Auzeville BP 52627, 31326 CASTANET-TOLOSAN CEDEX, France, ³Laboratorio de Interacción Hospedero-Patógeno/Unidad de Biología Molecular, Instituto Pasteur de Montevideo, Mataojo 2020, Montevideo 11400, Uruguay and ⁴Laboratorio de Fisiología de la Reproducción y Ecología de Peces, Instituto de Biología, Facultad de Ciencias, Universidad de la República Oriental del Uruguay, Iguá 4225, Montevideo 11400, Uruguay

† Corresponding author: Tel: +598 2 525 86 18 to 21: Fax: +598 2 525 86 17; Ext. 7150; Email: vizziano@gmail.com

Citation details: Klopp, C., Cabau, C., Greif, G., et al. Siberian sturgeon multi-tissue reference transcriptome database

Database (2020) Vol. 2016: article ID baaa082; doi:10.1093/database/baaa082

Received 22 June 2020; Revised 25 August 2020; Accepted 1 September 2020

## Abstract

**Motivation:** Siberian sturgeon is a long lived and late maturing fish farmed for caviar production in 50 countries. Functional genomics enable to find genes of interest for fish farming. In the absence of a reference genome, a reference transcriptome is very useful for sequencing based functional studies.

**Results:** We present here a high-quality transcriptome assembly database built using RNA-seq reads coming from brain, pituitary, gonadal, liver, stomach, kidney, anterior kidney, heart, embryonic and pre-larval tissues. It will facilitate crucial research on topics such as puberty, reproduction, growth, food intake and immunology. This database represents a major contribution to the publicly available sturgeon transcriptome reference datasets.

**Availability:** The database is publicly available at http://siberiansturgeontissuedb.sigenae.org

**Supplementary information:** Supplementary data are available at *Database* online.

## Introduction

The Siberian sturgeon, *Acipenser baerii*, is a non-teleost ray-finned fish (Actinopterygii) species of the order Acipenseriformes, which includes sturgeons and paddle-fishes (1) and faces critical conservation problems (2, 3). Siberian sturgeons from the easternmost of the three great Siberian rivers have been cultured in Europe since the early 1980s and this culture is now performed in 50 countries in the Northern and Southern Hemispheres (4). Its production

remains on the craft scale with worldwide production of only 27 500 tons per year (4), as compared to industrial levels of salmonid and tilapia production at 3 281 100 and 5 977 000 tons per year, respectively (5).

Knowledge of sturgeon physiology and genetics is less advanced than for other aquaculture species. For example, salmonids and African ciclids benefit from over 100 years of biological research, with genome assemblies published for Atlantic salmon (6), trout (7) and tilapia (8). In contrast
to the economic value of caviar, only one single sturgeon species, the sterlet, Acipenser ruthenus (9) received enough attention to have is genome sequenced. The sterlet has one of the smallest sturgeon genome sizes of 1.9 Gb compared to 2.3, 3.3 and 4.2 Gb of *Acipenser stellatus*, *Acipenser oxyrhynchus* and *Acipenser baerii*, respectively. The large variation between sturgeon genome sizes comes from recent whole genome duplications. The authors of (9) have shown that slow genome nucleotide content evolution renders genome assembly more difficult because of the large sections which are not separated in the assembly process and have to be split afterwards. Assembly complexity and genome size are causes for the lack of a Siberian sturgeon reference genome.

Reference transcriptomes have been used in many different species as for example the PhyloFish database (10), for fishes. Gene expression being tissue specific, reference transcriptomes are built with a mix of tissues. This mix can directly be performed on the RNA samples or *in silico* by combining the reads or the contigs produced for each tissue. Several software packages enable transcriptome assembly. The most used software package these days is Trinity (11) but others such as rnaSPAdes (12) and Oases (13) are also available. Once the reference is built its quality can be checked by verifying its compaction, read realignment rate and protein content. The read alignment rate should be close or over 90% for all of the good quality read sets. Protein content is usually checked with BUSCO (14). BUSCO searches the contigs for a set of nearly universal single-copy ortholog proteins, this for a given branch of the tree of life. For reference transcriptomes, the awaited BUSCO complete score is 90% or more. Assembly compaction is checked using contig count and length metrics as well as TransRate score (15).

Even if trinity is the most used RNA-Seq de novo assem-
bler these days, the transcriptome assembly problem is still
an open question. Several strategies can be used with differ-
ent combinations of read sets, software packages, filtering
procedure. The aim is to have the smallest possible number
of contigs, each one harboring a unique gene or transcript,
with an high read realignment rate. Several assemblies
should be performed to find the one.

In order to be able to understand sturgeon molecular
biological processes, in the absence of a reference genome,
several teams have produced transcriptomic assemblies
(16–36). Eight different sturgeon species already have
one publication including transcriptomic data. *Acipenser
schrenckii*, *Acipenser sinensis* and *Acipenser dabryanus*
have the highest count of this type of publications with
six, four and three, respectively. *Acipenser baerii* has only
two. Nine out of 21 of these publications are related to
gonad transcriptome, liver being the second most analyzed

tissue. The studies are usually organ specific because they
are focused on a given biological question. Only two stud-
ies include four tissues or more. Half of the studies have
published their raw data in an archive and only one-third
have published their contigs. No multi-tissue reference
transcriptome is available today for Siberian sturgeon.

Our aim is to present a high-quality Siberian stur-
geon reference transcriptome assembly database (SSTdb)
built using RNA-seq of different tissues to facilitate crucial
research on topics such as puberty, reproduction, growth,
food intake and immunology to improve management of
both wild and aquaculture populations.

Methods

Ethics statement, experimental procedures, rearing proce-
dures, RNA extraction, cDNA library construction and
Illumina sequencing have been described in (37). The
raw read data files have been published in SRA in the
PRJNA589958 BioProject and have the following identi-
fiers: SRX7158196 to SRX7158205. The contigs have
been published in TSA (transcriptome shotgun assembly
sequence) with prefix GICB01.

Transcriptome assembly and annotation

The RNA-seq datasets used to build the Siberian sturgeon multi-tissue database include brain, pituitary, immature testis, immature ovaries, liver, stomach, kidney, anterior kidney, heart, embryonic, pre-larval tissues and a pool of tissues. To remove read count variability between samples, assemblies were performed on a 20 million read-pairs subset per-sample for samples with read-pair counts above this value or all the reads for samples with counts below this value. Four assemblies were performed using the de novo RNA-Seq assembly pipeline (DRAP) 1.9 (38), using a combination of two assembly programs: trinity or Oases, and two assembly strategies: global assembly of the merged read sets or tissue-specific assemblies plus reconciliations of the contig sets. These strategies were performed using runDrap for the assemblies and runMeta for the reconciliations, both modules of the DRAP software package. The corresponding contig sets were named All_trinity, All_Oases, Meta_trinity and Meta_Oases.

The assembly metrics of the four resulting contig sets
were compared using runAssessment, the third DRAP
module, to choose the set with the best balance between
contig count, total size, read mapping rate, number of
matching sterlet proteins having 80% identity and 80%
coverage, BUSCO metrics and TransRate assembly met-
rics. The contig set selected had a low number of con-
tigs, a high read-mapping rate, a large number of aligned
proteins and high BUSCO and TransRate scores. The
actinopterygii BUSCO protein reference which comprises 4584 proteins has been used in this analysis. The chosen contig set was aligned on UniProt, Swiss-Prot, RefSeq and Ensembl Lepisosteus databases using BLASTX (39) for annotation and processed with InterProScan (40) to collect structural and functional annotations. The read sets were realigned on the contigs with BWA-MEM version 7.12-r1039 (standard parameters) (41). The SAM alignment files were compressed, sorted and indexed with SAMtools version 1.3 (42) standard parameters. The contig read counts were generated using the BAM files with SAMtools idxstats version 1.3 (standard parameters) and merged into a unique expression file with Unix Bash commands. The BAM file was processed with GATK version 3.0-0-g6bad1c6 (standard RNA-Seq parameters) (43) in order to find variants. All annotations, variations and expression measures were uploaded to RNAbrowse (44). The database is accessible through the following URL [http://siberiansturgeonontissuedb.sigenae.org/].

## Reference transcriptome validation

The reference transcriptome validation was performed by first realigning the reads of three different publicly available Siberian sturgeon datasets found in SRA and verifying the alignment rates, second by comparing contig quality metrics to other public sturgeon transcriptome assemblies found in TSA and finally by manually checking the gene content of the contig set. The read alignments were performed as described previously using BWA and SAMtools. The contig sets were compared with DRAP runAssessment module with the sterlet (Acipenser runthenus) NCBI reference protein file GCF_902713435.1_fAciRut3.1_paternal_haplotype_protein.faa and BUSCO. The manual check was performed by realigning the contigs on publicly available NCBI databases and verifying their completeness.

## Results

### Transcriptome assembly

We performed four assemblies combining two assemblers (trinity and Oases) and two strategies (merged reads and tissue assemblies reconciliation). The resulting contig sets were compared using different metrics produced by DRAP runAssessment module and are presented in Table 1. The number of contigs ranges from 57 996 to 105 556. The sum of contigs length is comprised between 104 981 273 and 146 660 170 bases pairs. These figures show the high assembly compaction variability even when using a pipeline, such as DRAP which reduces contig redundancy. The All_Oases (merged read files assembled by Oases) produced the lowest number of contigs and the largest sum of

contig lengths. The read mapping rates are much less variable and range from 81.58% to 85.94%. The Meta_trinity assembly performs best on this criterion but is also the assembly with the highest number of contigs. Its performance is just slightly over the Meta_Oases one, 85.94% and 85.64%, respectively. The mapping rates are lower than the expected 90% presented in the methods section because runAssessment uses bwa aln which is much more stringent than BWA-MEM. These assemblies are also very close on the properly paired alignment criteria, 79.45% versus 79.17%. In both cases, they are better than the all assemblies which have between two and four percents lower alignment rates. The Meta_oases assembly shows also the lowest number of sequence pairs having one end on one contig and the other end of the same pair on another contig. The Meta assemblies have the highest number of proteins alignment 8597 for Oases and 8486 for trinity significantly higher than the 8007 and 7212 of All_Oases and All_trinity, respectively. We have the same pattern for the BUSCO metrics with 89.5% of reconstructed genes for Meta_Oases and 89.3% for Meta_trinity, which is much higher than the 84.4% of All_Oases and the 80.3% of All_trinity. It is noticeable that the repartition of single copy and duplicated genes is highly in favor of the meta-assembly strategy. The All_Oases and in a lesser manner the All_trinity have a very high percentage of genes falling in BUSCO duplicated category, indicating that these genes are in multi-copy in the contig sets. Meta_Oases outperforms the others contig sets on the BUSCO criteria. Meta_Oases has also the best TransRate score (3143). Meta_Oases is first on many metrics and in second place on others and was therefore chosen as the reference contig set for the Siberian sturgeon multi-tissues database presented in this work. These results show a case in which the commonly found read merging trinity assembly strategy is out performed by others.

### Acipenser baerii public read sets realignment rates

Twenty-six Siberian sturgeon read sets have been retrieved from NCBI SRA. They come from three different projects PRJNA357627, PRJNA274436 and PRJNA589957 which have 1, 6 and 18 sample data file pairs, respectively. They include immature and mature male and female gonads, developing jaw, bone plate and five samples labeled as stemming from different developmental stages without further tissue information. Supplementary table ST1 presents the datasets with their project name, tissue number of reads, number of reads mapped on the Siberian sturgeon reference, number of properly paired reads as well as the corresponding mapping and pairing rates. The mapping
<table><caption>**Table 1.** Comparing the four contig sets build with two assemblers (trinity and Oases) and two strategies (one assembly for all the reads, one assembly per sample plus contig reconciliation)</caption><thead><tr><th>Metrics</th><th>Sample</th><th>All_Oases</th><th>All_trinity</th><th>Meta_Oases</th><th>Meta_trinity</th></tr></thead><tbody><tr><td rowspan="5">Assembly</td><td>N seq</td><td>57 996</td><td>75 514</td><td>71 263</td><td>105 556</td></tr><tr><td>N50</td><td>3431</td><td>1957</td><td>2501</td><td>1844</td></tr><tr><td>L50</td><td>14 316</td><td>16 772</td><td>14 784</td><td>21 757</td></tr><tr><td>Length sum</td><td>146 660 170</td><td>104 981 273</td><td>118 836 837</td><td>134 182 723</td></tr><tr><td>Length mean</td><td>2529</td><td>1390</td><td>1668</td><td>1271</td></tr><tr><td rowspan="2">Chimera</td><td>N chimeric contigs</td><td>140</td><td>140</td><td>228</td><td>123</td></tr><tr><td>N chimeric nt</td><td>35 448</td><td>40 559</td><td>72 474</td><td>35 620</td></tr><tr><td rowspan="3">Read alignment</td><td>Mapped</td><td>84.57%</td><td>81.58%</td><td>85.64%</td><td>85.94%</td></tr><tr><td>Properly paired</td><td>77.16%</td><td>75.00%</td><td>79.17%</td><td>79.45%</td></tr><tr><td>Mate mapped to different chromosome</td><td>8 062 249</td><td>11 643 098</td><td>6 720 793</td><td>7 454 996</td></tr><tr><td rowspan="5">BUSCO vertebrata_odb9</td><td>Complete</td><td>3872</td><td>3682</td><td>4103</td><td>4093</td></tr><tr><td>Complete single-copy</td><td>1675</td><td>2653</td><td>3415</td><td>3382</td></tr><tr><td>Complete duplicated</td><td>2197</td><td>1029</td><td>688</td><td>711</td></tr><tr><td>Fragmented</td><td>132</td><td>293</td><td>115</td><td>117</td></tr><tr><td>Missing</td><td>580</td><td>609</td><td>366</td><td>374</td></tr><tr><td>Proteins</td><td>N protein(s) aligned on contig</td><td>48 925</td><td>41 836</td><td>54 172</td><td>53 009</td></tr><tr><td rowspan="2">TransRate</td><td>Score</td><td>3115</td><td>2294</td><td>3143</td><td>2546</td></tr><tr><td>Optimal score</td><td>4356</td><td>366</td><td>4607</td><td>4395</td></tr></tbody></table>

rates range from 90.02% to 99.23% with an average of 94.28% (2.20%). The pairing rates range from 80.52% to 93.96% with an average of 88.52% (3.34%).

## Comparisons with other public sturgeon transcriptome assemblies

To assess the global quality of the assembly, we compared it with eight other publicly available sturgeon assemblies found in TSA: GEUL01, GGQL01, GGWJ01, GGWK01, GGYF01, GGZT01, GGZX01 and GICD01. The comparisons were based on two elements: assembly metrics including BUSCO and sterlet protein representation. Table 2 presents the assembly metrics and protein representation metrics for the eight public assemblies as well as for the Siberian sturgeon reference transcriptome named SSTdb. The assemblies have been performed using one or several

tissues, different assemblers, different cleaning or filtering strategies which explains the large variability of the figures found in the table. The number of contigs ranges from 53 634 to 641 485. The N50 varies five fold between 604 and 2891 base pairs. The total contig length starts at 34 807 151 base pairs and reaches 342 285 171 base pairs. SSTdb shows the best N50, contig mean length metrics as well as highest protein count. Figure 1 presents the BUSCO scores of these references. SSTdb has the best BUSCO scores with the fewest missing proteins and the lowest duplication rate.

## Manual validation

To validate the multi-tissue transcriptome assembly, we searched for key genes known to be expressed in the sequenced organs and tissues controlling functions such as reproduction, growth, food intake and immunity (see

<table><caption>**Table 2.** Public sturgeon transcriptome assembly comparison</caption><thead><tr><th>Assembly</th><th>Nb seq</th><th>N50</th><th>L50</th><th>Lg sum</th><th>Lg mean</th><th>Lg max</th><th>proteins</th></tr></thead><tbody><tr><td>GEUL01.1</td><td>179 564</td><td>1946</td><td>25 772</td><td>166 715 666</td><td>928</td><td>54 403</td><td>46 841</td></tr><tr><td>GGQL01.1</td><td>77 634</td><td>2523</td><td>18 243</td><td>135 647 056</td><td>1747</td><td>16 131</td><td>41 703</td></tr><tr><td>GGWJ01.1</td><td>53 624</td><td>1086</td><td>8322</td><td>34 807 151</td><td>649</td><td>15 639</td><td>17 213</td></tr><tr><td>GGWK01.1</td><td>369 441</td><td>763</td><td>69 627</td><td>208 011 161</td><td>563</td><td>15 644</td><td>29 402</td></tr><tr><td>GGYF01.1</td><td>121 398</td><td>2211</td><td>23 671</td><td>168 641 140</td><td>1389</td><td>2596</td><td>49 218</td></tr><tr><td>GGZT01.1</td><td>203 131</td><td>1874</td><td>40 438</td><td>254 007 803</td><td>125</td><td>34 023</td><td>30 551</td></tr><tr><td>GGZX01.1</td><td>641 485</td><td>604</td><td>132 361</td><td>342 285 171</td><td>533</td><td>1664</td><td>20 915</td></tr><tr><td>GICD01.1</td><td>91 579</td><td>1011</td><td>18 288</td><td>136 604 581</td><td>1492</td><td>20 419</td><td>53 729</td></tr><tr><td>SSTdb</td><td>79 217</td><td>2891</td><td>16 359</td><td>150 824 770</td><td>1903</td><td>45 872</td><td>54 172</td></tr></tbody></table>

The public assembly names with their TSA prefix. The Siberian sturgeon reference transcriptome database is named SSTdb.
BUSCO Assessment Results

<figure><img src="image_2.png" /><figcaption>Figure 1. BUSCO scores of the eight sturgeon transcriptome assemblies found in TSA plus the sturgeon reference database named SSTbd.</figcaption></figure>

<table><caption>Table 3. Reference transcriptome manual validation table</caption><thead><tr><td>Type</td><td>Number of genes searched</td><td>Number of contigs found<sup>a</sup></td><td>Confirmed genes (%)</td></tr></thead><tbody><tr><td>Hypophysiotropic peptides</td><td>17</td><td>24</td><td>88</td></tr><tr><td>Hypophysiotropic peptide receptors</td><td>17</td><td>18</td><td>88</td></tr><tr><td>Pituitary hormones</td><td>8</td><td>8</td><td>100</td></tr><tr><td>Gonad related</td><td>11</td><td>11</td><td>100</td></tr><tr><td>Liver</td><td>7</td><td>10</td><td>100</td></tr><tr><td>Gastrointestinal hormones genes</td><td>10</td><td>11</td><td>100</td></tr><tr><td>Kidney and anterior kidney</td><td>7</td><td>8</td><td>86</td></tr><tr><td>Immunologically-relevant genes</td><td>8</td><td>13</td><td>100</td></tr></tbody></table>

Supplementary table). We searched 85 genes. Table 3 presents types and number of manually verified genes. For five types out of eight, 100% of the searched genes were confirmed in the Tissue database. For the three other types, the confirmation rates were above 86%.

Among hypophysiotropic peptides related to reproductive control and produced in the brain (Supplementary table, Section 2), we identified novel sequences for sturgeons for the gonadotropin-inhibitory hormone (gnih), a relevant factor inhibiting puberty in birds (Tsutsui et al. 2000), together with two receptors of GnIH: neuropeptide FF receptor 1-like, neuropeptide FF receptor 2-like. We are also reporting for the first time in sturgeons the GnRH1 receptor, and the GPR54 or kiss-receptor 1 (Supplementary

table, Section 2) that will help to complete some pathways to understand the process of reproduction in sturgeons at molecular level. We found all the sequences searched for the pituitary hormones (Supplementary table, Section 3), and the main receptors for pituitary hormones (Supplementary table, Section 4). In addition, sequences of steroid-related genes and vitellogenin receptor have been found in gonads (Supplementary table, Section 5) together with sequences of estrogen receptors and of the vitellogenin in the liver (Supplementary table, Section 6).

Concerning the central growth hormone (GH)-releasing and inhibiting factors, we found GH-releasing factor sequences for pituitary adenylate cyclase-activating polypeptide, thyrotropin-releasing hormone receptor (thr),

<footnote>ªFor some genes more than one contig have been found in the assembly.</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Home page of the Siberian sturgeon transcriptome database website. The main menu enables to visualize and query contigs and variants as well as download the raw and result files.</figcaption></figure>

corticotropin-releasing hormone (crh) and cholecystokinin (cck) (Supplementary table, Section 1, 7). The new receptor sequences reported in sturgeons are: the pituitary adenylate cyclase-activating polypeptide type I receptor, thyrotropin-releasing hormone receptor, corticotropin-releasing factor receptor 1 and somatostatin receptor type 5-like (Supplementary table, Section 2). We also confirmed several factors affecting the peripheral control of GH release as insulin-like growth factor 1, ghrelin, enzymes involved in cortisol and estradiol production (Supplementary table, Section 5, 6, 7, 8).

The present study identified sequences for all of the major food intake-stimulating neuropeptides (neuropeptide Y, pro-melanin-concentrating hormone), and intake-inhibiting peptides (proopiomelanocortin A precursor, cocaine- and amphetamine-regulated transcript protein, corticotropin-releasing factor receptor 1).

Among the immune system genes known to be expressed in the liver, this work describes three genes from the serum amyloid A protein family for the first time in sturgeons (serum amyloid A protein precursor, acute phase serum amyloid A, serum amyloid A-5 protein, amyloid A) and two contigs of the serum amyloid P-component-like. We also identified and lysozyme-C like, and lysozyme sequences,

(Supplementary table, Section 9) previously reported in GenBank (MF135 537.1, MF280 234.1). In terms of stress response, we identified the mineralocorticoid receptor and glucocorticoid receptor sequences (Supplementary table, Section 8).

## Discussion

Transcriptome assembly is widely used to study molecular mechanisms of species lacking a reference genome. Trinity has become one of the most used software packages to assemble de novo transcriptomes. Our assembly test comparing a combination of two assemblers and two assembly strategies shows that, in this case, the most used solution is not the one providing the best metrics. It is important for a reference assembly to show a low redundancy and at the same time to have the widest possible gene representation. The Meta_Oases assembly we chose as representative assembly for Sturgeon multi-tissue transcriptome database may not be the best for all metrics but it shows the best compromise between compactness and completeness. DRAP runMeta procedure proved to be able to select the best contigs coming from different tissues to make a good quality assembly reconciliation.
A reference transcriptome to be used in numerous projects should show high alignment metrics for new-read sets produced in conditions or tissues not used in its building process. Displaying over 90% of read alignment for the 26 publicly available read sets used in the validation, our transcriptome complies with this expectation. This verification can be performed for novel projects in order to confirm its qualification.

Compared to the other public sturgeon transcriptome references available in TSA, SSTdb has the best N50 and mean length metrics as well as the best aligned sterlet protein count, this with the third lowest contig count. The second and third assemblies regarding the contig length metrics are GGQL01.1 and GGYF01.1, respectively. SSTdb had also the best BUSCO profile with the highest complete gene count and the low duplication rate. Only GICD01.1 has a similar profile but not as good. GICD01.1 has been assembled from multiple male, female differentiated and undifferentiated gonad using the same strategy as SSTdb.

The manual biological validation shows that SSTdb contains many well-structured expected genes which will facilitate the study of different functions to improve sturgeon farming, i.e. puberty, reproduction, growth and food intake and immune system.

The public website enables to query and visualize different contig features including different annotations, expression along the contig for all included conditions and variation locations. Text queries can be performed through the biomart interface and sequence queries using blastx or tblastn depending if the query bait is a protein or nucleic sequence. The raw result files including contigs in fasta, reads in fastq, alignments in bam, variants in VCF format can also be downloaded.

The Siberian sturgeon multi-tissue assembly represents a major contribution to the publicly available sturgeon transcriptome references. This is the first work that releases a high-quality multi-tissue database in order to facilitate basic gene studies to sustain Siberian sturgeon aquaculture. Read sets as well as contigs will also be a valuable input for the Siberian sturgeon genome annotation.

## Acknowledgement

Many thanks are due to the chief executive officer of Estuario del Plata Facundo Márquez, to Dr Vet A.R., Mr E.O., and all of the Estuario del Plata (Uruguay) staff for their ample support and kindness during the development of this work.

## Author contribution

Christophe Klopp (CK) and Denise Vizziano Cantonnet (DVC) conceived and designed the experiments; DVC directed the research project; André Lasalle, Gonzalo Greif

and Santiago Di Landro performed the experiments; CK, Cédric Cabau (CC) and DVC analyzed the data and organized the datasets. The manuscript was written by CK, CC and DVC.

## Funding

This work was supported by the Comisión Sectorial de Investigación Científica (CSIC), Universidad de la República Oriental del Uruguay, (Grant number C225-348), the Agencia Nacional de Investigación e Innovación-Dirección Nacional de Recursos Acuáticos ANII-DINARA [Grant number FPA 9975] and Agencia Nacional de Investigación e Innovación [Grant number FMV_1_2017_1_135 908].

## Competing interests

The authors have no potential conflicts of interest to declare.

## References

1. Nelson, J. (2006). *Fishes of the World*, John Wiley, Hoboken, N.J.

2. Birstein, V.J. (1993) Sturgeons and paddlefishes: threatened fishes in need of conservation. *Conservation Biology*, **7**, 773–787.

3. Sokolov, L.I. et al.. (1989) Acipenser baeri brandt, 1869. the freshwater fishes of Europe: general introduction to fishes acipenseriformes.. *Book*, Vol. 1, pp. 263–284.

4. Chebanov, M. and Williot, P.. (2018) *An Assessment of the Characteristics of World Production of Siberian Sturgeon Des-tined to Human Consumption*, Springer International Publishing, Cham, pp. 217–286.

5. Tveteras, R., (2016) Global fish production data & analysis. *Presented at the Global Outlook for Aquaculture Leadership Conference in Guangzhou*, China

6. Lien, S. et al.. (2016) The Atlantic salmon genome provides insights into rediploidization. *Nature*, **533**, 200–5.

7. Berthelot, C. et al.. (2014) The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates. *Nat Commun*, **5**, 3657.

8. Brawand, D. et al.. (2014) The genomic substrate for adaptive radiation in African cichlid fish. *Nature*, **513**, 375–381.

9. Du, K. et al. (2020). The sterlet sturgeon genome sequence and the mechanisms of segmental rediploidization. *Nat Ecol Evol*, **4**, 841–852.

10. Pasquier, J. et al.. (2016) Gene evolution and gene expression after whole genome duplication in fish: the PhyloFish database.. *BMC Genomics*, **17**, 368.

11. Grabherr, M.G. et al.. (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome.. *Nat. Biotechnol.*, **29**, 644–52.

12. Bushmanova, E. et al. (2019). rnaSPAdes: a de novo transcriptome assembler and its application to RNA-Seq data. *Gigascience*, **8**(9), September 2019, giz100.

13. Schulz, M.H. et al.. (2012) Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels.. *Bioinformatics*, **28**, 1086–92.
14. Seppey, M. et al. (2019) BUSCO: Assessing genome assembly and annotation completeness. *Methods Mol. Biol.*, **1962**, 227-245.

15. Smith-Unna, R. et al. (2016) TransRate: reference-free quality assessment of de novo transcriptome assemblies. *Genome Res.*, **26**, 1134-44.

16. Chen, Y. et al. (2018) Gonadal transcriptome sequencing of the critically endangered Acipenser dabryanus to discover candidate sex-related genes. *PeerJ*, **6**, e5389.

17. Chen, H. et al. (2019) Transcriptome analysis and the effects of polyunsaturated fatty acids on the immune responses of the critically endangered angtze sturgeon (Acipenser dabryanus). *Fish Shellfish Immunol.*, **94**, 199-210.

18. Degani, G. et al. (2019) Sex-related gonadal gene expression differences in the Russian sturgeon (Acipenser gueldenstaedtii) grown in stable aquaculture conditions. *Anim. Reprod. Sci.*, **200**, 75-85.

19. Du, H. et al. (2019) Hypothalamus-pituitary-gonad axis transcriptome profiling for sex differentiation in Acipenser sinensis. *Sci Data*, **6**, 87.

20. Hale, M.C. et al. (2009) Next-generation pyrosequencing of gonad transcriptomes in the polyploid lake sturgeon (Acipenser fulvescens): the relative merits of normalization and rarefaction in gene discovery. *BMC Genomics*, **10**, 203.

21. Jiang, N. et al. (2018) Transcriptome analysis of Aeromonas hydrophila infected hybrid sturgeon (Huso dauricus x Acipenser schrenckii). *Sci Rep*, **8**, 17925.

22. Jin, S. B. et al. (2015). Comparative transcriptome analysis of testes and ovaries for the discovery of novel genes from Amur sturgeon (Acipenser schrenckii). *Genet. Mol. Res.*, **14**, 18913-27.

23. Kaitetzidou, E. et al. (2017) Expression Patterns of Atlantic Sturgeon (Acipenser oxyrinchus) During Embryonic Development. *G3 (Bethesda)*, **7**, 533-542.

24. Luo, K. et al. (2018). Transcriptome analysis of the critically endangered Dabry's sturgeon (Acipenser dabryanus) head kidney response to Aeromonas hydrophila. *Fish Shellfish Immunol.*, **83**, 249-261.

25. Mugue, N. et al. (2019). Transcriptome sequencing of hybrid bester sturgeon: Responses to poly (I:C) in the context of comparative immunogenomics. *Fish Shellfish Immunol.*, **93**, 888-894.

26. Song, W. et al. (2016). RNA-sequencing of the sturgeon Acipenser baeri provides insights into expression dynamics of morphogenic differentiation and developmental regulatory genes in early versus late developmental stages. *BMC Genomics*, **17**, 564.

27. Vidotto, M. et al. (2013) Transcriptome sequencing and de novo annotation of the critically endangered Adriatic sturgeon. *BMC Genomics*, **14**, 407.

28. Vizziano-Cantonnet, D. et al. (2018) De novo transcriptome analysis to search for sex-differentiation genes in the Siberian sturgeon. *Gen. Comp. Endocrinol.*, **268**, 96-109.

29. Yang, J. et al. (2018) Transcriptome profiling reveals candidate cleft palate-related genes in cultured Chinese sturgeons (Acipenser sinensis). *Gene*, **666**, 1-8.

30. Yebra-Pimentel, E.S. et al. (2019) Deep transcriptome analysis of the heat shock response in an Atlantic sturgeon (Acipenser oxyrinchus) cell line. *Fish Shellfish Immunol.*, **88**, 508-517.

31. Yue, H. et al. (2015) Sequencing and De Novo Assembly of the Gonadal Transcriptome of the Endangered Chinese Sturgeon (Acipenser sinensis). *PLoS ONE*, **10**, e0127332.

32. Zhang, Q. et al. (2019a) De novo Assembly and Analysis of Amur Sturgeon (Acipenser schrenckii) Transcriptome in Response to Mycobacterium Marinum Infection to Identify Putative Genes Involved in Immunity. *J. Microbiol. Biotechnol.*, **29**, 1324-1334.

33. Zhang, X. et al. (2019b) Transcriptome analysis provides insights into differentially expressed genes and long noncoding RNAs involved in sex-related differences in Amur sturgeon (Acipenser schrenckii). *Mol. Reprod. Dev.*, **86**, 132-144.

34. Zhang, X. et al. (2020) Full-length transcriptome sequencing and comparative transcriptomic analysis to uncover genes involved in early gametogenesis in the gonads of Amur sturgeon (Acipenser schrenckii). *Front. Zool.*, **17**, 11.

35. Zhang, X. J. et al. (2016). Transcriptome analysis and de novo annotation of the critically endangered Amur sturgeon (Acipenser schrenckii). *Genet. Mol. Res.*, **15**(2), 15027999.

36. Zhu, R. et al. (2016) De novo annotation of the immune-enriched transcriptome provides insights into immune system genes of Chinese sturgeon (Acipenser sinensis). *Fish Shellfish Immunol.*, **55**, 699-716.

37. Klopp, C. et al. (2020) Multi-tissue siberian sturgeon rna sequencing data. *Data in Brief*, **31**, 105820.

38. Cabau, C. et al. (2017) Compacting and correcting Trinity and Oases RNA-Seq de novo assemblies. *PeerJ*, **5**, e2988.

39. Altschul, S.F. et al. (1990) Basic local alignment search tool. *J. Mol. Biol.*, **215**, 403-10.

40. Hunter, S. et al. (2012) InterPro in 2011: new developments in the family and domain prediction database. *Nucleic Acids Res.*, **40**, (Database issue) D306-12.

41. Li, H. and Durbin, R. (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics*, **25**, 1754-60.

42. Li, H. et al. (2009) The Sequence Alignment/Map format and SAMtools. *Bioinformatics*, **25**, 2078-9.

43. McKenna, A. et al. (2010) The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. *Genome Res.*, **20**, 1297-303.

44. Mariette, J. et al. (2014) RNAbrowse: RNA-Seq de novo assembly results browser. *PLoS ONE*, **9**, e96821.