The transcript catalogue of the short-lived fish Nothobranchius furzeriprovides insights into age-dependent changes of mRNA levels
© Petzold et al.; licensee BioMed Central Ltd. 2013
Received: 19 July 2012
Accepted: 5 March 2013
Published: 16 March 2013
The African annual fish Nothobranchius furzeri has over recent years been established as a model species for ageing-related studies. This is mainly based on its exceptionally short lifespan and the presence of typical characteristics of vertebrate ageing. To substantiate its role as an alternative vertebrate ageing model, a transcript catalogue is needed, which can serve e.g. as basis for identifying ageing-related genes.
To build the N. furzeri transcript catalogue, thirteen cDNA libraries were sequenced using Sanger, 454/Roche and Solexa/Illumina technologies yielding about 39 Gb. In total, 19,875 protein-coding genes were identified and annotated. Of these, 71% are represented by at least one transcript contig with a complete coding sequence. Further, transcript levels of young and old fish of the strains GRZ and MZM-0403, which differ in lifespan by twofold, were studied by RNA-seq. In skin and brain, 85 differentially expressed genes were detected; these have a role in cell cycle control and proliferation, inflammation and tissue maintenance. An RNA-seq experiment for zebrafish skin confirmed the ageing-related relevance of the findings in N. furzeri. Notably, analyses of transcript levels between zebrafish and N. furzeri but also between N. furzeri strains differed largely, suggesting that ageing is accelerated in the short-lived N. furzeri strain GRZ compared to the longer-lived strain MZM-0403.
We provide a comprehensive, annotated N. furzeri transcript catalogue and a first transcriptome-wide insight into N. furzeri ageing. This data will serve as a basis for future functional studies of ageing-related genes.
Invertebrate model species like yeast, worm and fly have been successfully used to study the genetic background of lifespan determination and ageing, leading to the identification of a number of relevant genes and pathways . Although several of these pathways are conserved across taxa [2–4], reproducing the results in vertebrate species has not always been possible. A major drawback of vertebrate models like mouse or zebrafish is their relatively long lifespan of up to five years, which makes the validation of results obtained in invertebrates time-consuming and expensive. Alternative vertebrate model organisms are needed that allow validation of known and discovery of potentially new ageing-related pathways/genes within a reasonable time frame.
In 2005, the Turquoise killifish, Nothobranchius furzeri, was proposed as a new model for studies of ageing . N. furzeri inhabit temporary ponds in South-East Africa and were first collected in 1968 in the Game Reserve Gona Re Zhou (GRZ) in Zimbabwe. Due to the seasonal nature of the habitat, this species shows an annual life cycle in which hatching, growth and reproduction happen during the wet season and fertilised eggs survive the dry season in the dry mud by entering into diapause . Consequently, the lifespan of African N. furzeri populations is limited by the duration of the wet season, which lasts up to several months depending on location. Remarkably, the short-lived phenotype has been retained in captivity despite unlimited water and food supply, and N. furzeri laboratory strains originally derived from dry/humid habitats differ in lifespan . For example, the GRZ strain, which stems from a dry habitat, has a maximum lifespan (10% survivorship) of 12-16 weeks [8, 9]. The more recently isolated strain MZM-0403, originally collected in a humid habitat in Southern Mozambique (MZM), has a maximum lifespan of 25-32 weeks , which still is exceptionally short compared to other vertebrate species. Despite their short lifespan, the fish show typical ageing-related phenotypes such as physiological / cognitive decay and expression of ageing-related biomarkers .
Over past years, genetic and genomic resources for N. furzeri have been generated. A first insight into cytogenetic and genomic characteristics, including karyotype (19 chromosomes, 2n = 38), genome size (1.6–1.9 Gb) and composition (45% repeats), phylogeny (among fish with a sequenced genome available, medaka is the closest relative) were provided by our group . Further, genetic variability in GRZ (highly inbred) and MZM-0403 (resembling a wild population) was assessed and a genetic linkage map, based on microsatellites and gene-associated single nucleotide variations, was created and used to identify several quantitative trait loci affecting lifespan . On the other hand, sequence resources are still limited, i.e. a search for Nothobranchius furzeri in NCBI Entrez lists 13 genes, 279 nucleotide and 236 protein entries (NCBI Entrez, 18/02/13). Here, we set out to build a comprehensive, high-quality, annotated N. furzeri transcript catalogue. We used the catalogue in a pilot RNA-seq experiment aimed at identifying ageing-related genes, i.e. genes differentially expressed in the short- and longer-lived N. furzeri strains GRZ and MZM-0403 at young and old age.
Assembly of the N. furzeritranscript catalogue
N. furzeri transcriptome libraries and sequencing methods
Specimens used for RNA-isolation
Number of reads
Average length (bp)
Total length (Mb)
one male GRZ, 9 weeks, whole body; normaliseda
one male GRZ, 10 weeks, whole body; normaliseda
one male GRZ, 10 weeks, whole body
six male and female GRZ, 8 weeks, whole body; normalisedb
three GRZ, 1 week, whole body
Solexa/Illumina, 2 x 150 bp
five male GRZ, 5 weeks, skin
Solexa/Illumina, 2 x 101 bp
Solexa/Illumina, 76 bp
five male GRZ, 5 weeks, brain
Solexa/Illumina, 2 x 101 bp
Solexa/Illumina, 76 bp
three male GRZ, 14 weeks, skin
Solexa/Illumina, 76 bp
three male GRZ, 14 weeks, brain
Solexa/Illumina, 76 bp
four male MZM-0403, 5 weeks, skin
Solexa/Illumina, 76 bp
five male MZM-0403, 5 weeks, brain
Solexa/Illumina, 76 bp
five male MZM-0403, 31 weeks, skin
Solexa/Illumina, 76 bp
three male MZM-0403, 31 weeks, brain
Solexa/Illumina, 76 bp
Assembly metrics a
First assembly a
Second assembly b
Iterative assembly procedure
PAVE contigs, Solexa/Illumina
Number of contigs (≥ 300 bp)
Total length of contigs (≥ 300 bp)
Average length of contigs
Median length of contigs
Maximum length of contigs
Number of large contigs (≥ 1 kb)
Total length of large contigs
To quantify the impact of the three different sequencing technologies on our transcriptome assembly, all reads were mapped (Sanger and 454/Roche with Newbler; Solexa/Illumina with BWA ) onto the contigs of the final assembly. As a result, the vast majority of reads was mapped successfully (Sanger: 99%, 454/Roche: 92%, Solexa/Illumina: 97%). Sanger, 454/Roche and Solexa/Illumina reads covered about 20, 56 and 237 Mb, respectively, of the 253 Mb total contig length. Thereof, only 1, 7 and 178 Mb, respectively, were covered exclusively by either of the three technologies, showing the major impact of Solexa/Illumina on the final transcriptome assembly.
To assess the quality of the transcriptome assembly, transcript contigs were compared against a set of reference sequences obtained from the Ensembl genome annotation of medaka , the closest relative of N. furzeri. To avoid misinterpretation due to annotation errors of the medaka genome, we limited our analysis to the 1,750 medaka protein sequences with experimental support (Ensembl category “known”). All N. furzeri transcript contigs were compared with BLASTx (e-value ≤ 10-20) against the medaka protein sequences and evaluated for representation and coverage, assembly fragmentation and redundancy as well as putative chimerism. Altogether, 13,451 N. furzeri transcript contigs produced hits for 1,290 (74%) medaka proteins. Of these, 74% were covered by hits over ≥ 90% of their length. Moreover, for 67%, a hit by a single transcript contig was sufficient to achieve this coverage. To assess fragmentation, we determined the fractions of transcript contigs covering either the N-terminus (13%) or the C-terminus (16%) or neither one (64%) of the medaka protein. On the other hand, 7% were spanning the complete protein and included UTRs. In contrast to this relative low degree of fragmentation of the transcript catalogue, over 92% of the total identified protein sequence was hit by at least two transcript contigs with an average redundancy of nine contigs per aa. Finally, the BLASTx results facilitated identifying putative chimeric transcript contigs. To be marked as chimeric, a transcript contig was required to have two non-overlapping hits to different medaka proteins which meet the following criteria: e-value ≤ 10-20, coverage of both protein hits ≥ 75% and identity ≥ 50%. According to this definition, only 0.3% of the transcript contigs were identified as chimeric.
To identify transcript contigs derived from contaminating sources, tBLASTx searches (e-value ≤ 10-20) were performed against a database compiled from both fish-specific transcripts from Ensembl as well as non-fish transcripts from RefSeq . For 1.9% of the contigs (4 Mb, median length: 531 bp), the best hit represented a non-fish species and the e-value of the next fish-specific hit was at least 10 orders of magnitude higher. Most hits were found for mammals (0.8%), followed by invertebrates (0.7%) and plants (0.3%). As the latter two may originate from tissue contamination by fish’s environment, those tagged as mammalian contamination represent rather laboratory artefacts acquired e.g. by the multistep cDNA normalization procedure prior Sanger sequencing (data not shown). The contigs were excluded from further analyses. Accordingly, 209,660 N. furzeri transcript contigs (249 Mb) were analysed further.
Annotation of the N. furzeritranscript contigs
Annotation was based on BLAST searches against the following protein and nucleotide databases: (i) protein and transcript sequences of medaka, stickleback, tetraodon and zebrafish from Ensembl , (ii) human protein sequences from RefSeq , (iii) the UniProt database , (iv) the NCBI non-redundant protein database nr  and (v) fish-specific ESTs from NCBI UniGene . BLAST parameters and e-value thresholds were defined for each database separately (Additional file 1: Table S1).
BLAST based annotation of N. furzeri transcript contigs
Database entries hit
Database entries used
Ensembl fish proteins
NCBI nr proteins
Refseq human proteins
NCBI UniGene transcripts
Ensembl fish transcripts
For each transcript contig, BLAST hits were filtered and the best was selected for annotation. That way, 109,032 (52%) transcript contigs were annotated using 46,218 database entries. Of these, 21,728 (47%) were hit by more than one transcript contig, with an average of four per database entry. This reflects rather the redundancy than the fragmentation of the assembly (see above) and in respect to the protein hits provides an estimate of protein isoforms (see below). Subsequently, each transcript contig was assigned a gene symbol based on the information of the database entry used for its annotation. As a result, the N. furzeri transcript catalogue contains transcript contigs for 19,875 protein-coding genes. About half of all genes (10,392/52%) had multiple isoforms, with an average of three isoforms per gene.
Gene ontology of the N. furzeritranscript contigs
To strengthen the annotation of the N. furzeri transcript catalogue, contigs were assigned to Gene Ontology (GO, ) categories. The GO terms were obtained from the respective Ensembl and UniProt database entries used for the annotation of the transcript contigs. Of 109,032 annotated transcript contigs, 76,013 (72%) were linked to 340,452 GO counts. “Molecular Function” represented the largest class of assignments (70,305 contigs/ 166,476 GO counts), followed by “Biological Process” (49,972/108,336) and “Cellular Component” (40,551/65,640). After removing redundant counts, 6,796 unique GO terms were identified.
These GO terms were further mapped to GO-Slim terms, a set of 131 high-level functional annotations. In respect to the longest transcript contig per N. furzeri gene, 137,683 GO-Slim counts were assigned to 15,040 transcript contigs representing 124 unique GO-Slim terms (Additional file 2: Figure S1). Most frequent terms among the three GO domains were “metabolic process” (6,108 contigs, domain: biological process), “binding” (10,567 contigs, biological function) and “cell” (7,668 contigs, cellular component). All GO annotations are available through the “Nothobranchius furzeri Information Network transcriptome browser (see below).
Comparison of N. furzeritranscriptome annotation to model fish genomes
BLASTx comparison of N. furzeri transcript contigs to the protein/gene annotation of four other fish genomes
Number of proteins/ genes annotated
24,661 / 19,686
27,576 / 20,787
23,118 / 19,602
41,478 / 26,095
Contigs with hits
19,272 / 17,173
20,534 / 17,795
18,613 / 16,744
22,669 / 18,620
78% / 87%
74% / 86%
81% / 85%
55% / 71%
Average protein identity
Contigs with hits exclusively in this species
914 / 914
690 / 690
356 / 356
942 / 942
Average protein identity
Since the comparison of gene numbers relies on the level of annotation and therefore can be biased, we evaluated the number of gene families, i.e. groups of genes with a common evolutionary ancestor. Gene family information for medaka, stickleback, tetraodon and zebrafish was obtained from the Ensembl Compara gene trees . Zebrafish showed the highest number of gene families (8,241), followed by stickleback (7,663), medaka (7,501) and tetraodon (7,261). Based on the BLASTx results, we identified 7,489 gene families hit by at least one N. furzeri transcript contig. Moreover, for 85% (6,428 families), at least one of the four model fish species shared the same number of genes per family. These results indicate a completeness of the provided transcript catalogue comparable to that of model fish species for which an annotated genome is available.
Duplicated N. furzerigenes
Teleost fish genomes underwent an additional round of genome duplication [24, 25], compared to other vertebrates, and thus contain a high number of paralogs. Genes duplicated in N. furzeri were assessed in two steps. First, we identified those duplicated also in at least one of the model fish species. Towards this, for each N. furzeri gene the longest transcript contig was compared using BLASTx against Ensembl proteins of medaka, stickleback, tetraodon and zebrafish, and the Ensembl gene of the best hit was recorded. The corresponding Ensembl Compara gene tree was then used to find paralogs of this gene which are specific for teleost fish. If any of these paralogs were hit by another N. furzeri gene, then the two N. furzeri genes were considered as paralogous. By this, we identified among the 17,261 N. furzeri genes with BLASTx hits 4,498 paralogs grouped in 2,148 paralog families (2.1 paralogs per family, on average; Additional file 4). Second, among the remaining genes we searched redundant transcript contigs for indications of N. furzeri-specific paralogs using a clustering approach based on protein identities. For each gene, all transcript contigs were translated into proteins and aligned. Protein identities of these alignments were clustered based on a minimum distance of 10% between each cluster. For 1,249 genes more than one cluster was found, corresponding to a total of 2,619 (2.1 contig clusters per gene, on average; Additional file 4). Overall, this indicates that the currently annotated transcript catalogue consists of 17,525 (19,875-4,498 + 2,148) protein-coding genes if paralogs are counted only once. Of those, 19.1% (2,148 + 1,249) showed evidences of multiplication, whereas 90% (4,134 + 2,274) of the 7,117 (4,498 + 2,619) paralogs occurred as pairs. The remaining 709 (364 + 345) paralogs were grouped in 193 paralog families consisting of 3 and more paralogs (average: 4 paralogs per family, maximum: 23).
N. furzeritranscriptome browser - NFINtb
A public web site interface named the “Nothobranchius furzeri Information Network transcriptome browser” (NFINtb, https://gen100.imb-jena.de/EST2UNI/nfintb/, Additional file 5: Figure S3) was set up. It provides the annotated transcript catalogue, including detailed BLAST results, predicted protein domains, associated gene ontology terms, microsatellites and gene expression data. Further information on sequenced libraries, distribution of contig lengths and annotation numbers, respective GO terms and corresponding transcript contigs is given. A query mask allows keyword (e.g. annotation, gene symbol, GO category) and BLAST searches by sequence similarity. N. furzeri sequences can be retrieved by batch download including annotation details.
Differentially expressed genes in old vs. young N. furzeri
Summary of RNA-seq analysis
Mapped / Unique (%)
58 / 46
42 / 37
53 / 43
41 / 36
65 / 52
54 / 48
64 / 53
52 / 46
Median transcript level (rpkm a)
Range of transcript levels (rpkm)
0.005 – 13,043
0.006 – 40,801
0.005 – 39,533
0.01 – 5,587
0.006 – 10,409
0.003 – 5,288
0.005 – 3,878
0.009 – 3,946
To identify genes statistically significant differentially expressed with age, the available RNA-seq datasets were lacking biological replicates as the primary goal of these sequencing efforts was the completion of the transcript catalogue (3–5 specimens were pooled before sequencing, Table 1). To overcome this limitation, we decided to abstract from strain-specific differences and used the respective datasets of the two N. furzeri strains as replicates. In total, 43 differentially expressed genes (DEGs) were found in each of the two tissues, respectively (p ≤ 0.01, DESeq ; Additional file 6). In skin, 11 genes were up- and 32 genes downregulated with age. In brain, there were 24 up- and 19 downregulated genes. The intersection of the DEG sets contained only VDRB (vitamin D receptor b) which was upregulated in both tissues.
To identify biological functions potentially affected by the changes in transcript levels, DEGs were submitted to the Database for Annotation, Visualization and Integrated Discovery (DAVID, [27, 28]). DEGs were analysed for shared gene ontology terms, pathways and protein domains, and resulting enriched gene groups (minimum enrichment score ≥ 1) were considered. The 34 upregulated genes grouped into the following three annotation clusters: “activation of immune response”, “cell adhesion at the plasma membrane” and “positive regulation of apoptosis”. The first cluster was identified by genes upregulated in brain whereas the others contained genes from both tissues. The 51 downregulated genes seemed to act mostly as “collagens/proteoglycans” and were located in the “extracellular region”; these functional clusters were identified by skin DEGs only. Additional annotation clusters of downregulated genes from both tissues described processes like “cell cycle”, “cell division” and “cell proliferation” as well as “DNA replication” (Additional file 7).
To evaluate whether the observed N. furzeri ageing signature is specific to this species or more general for fish, we performed a similar analysis in zebrafish which serves a model organism for studying many aspects of biology, e.g. ageing. To consider in this case biological replicates from the beginning, we generated RNA-seq data from skin of 5 and 42 month old individual zebrafish (10 and 5 specimens, respectively). In total, 1,191 zebrafish DEGs were found (p ≤ 0.01). This number is much higher than in N. furzeri, which is likely due to the presence of a more appropriate selection of biological replicates. Of the 43 DEGs identified in N. furzeri skin, 19 (43%) were also present in zebrafish DEGs and showed the same direction in fold change (Additional file 9). This was significantly different from chance (p = 2.4x10-14, binomial test). Of the 24 remaining genes, additional 16 (33%) showed the same direction, but did not reach significance in zebrafish. One of the remaining eight genes (cytokine-like 1, CYTL1) was significant in both species with a different direction (N. furzeri: downregulated, zebrafish: upregulated). To summarise, we found that a significant number of DEGs was similarly regulated in ageing of N. furzeri and zebrafish.
Strain-specific signatures of ageing in N. furzeri
We built a comprehensive, annotated transcript catalogue for N. furzeri, the vertebrate species with the shortest lifespan recorded for vertebrates in captivity, and identified 85 genes showing significantly changed transcript levels over the N. furzeri lifetime.
Developing an assembly strategy
Transcriptome assembly is still a challenging task, mainly because most available tools are designed for assembling genomes and therefore assume an even coverage distribution and error rate . Transcript levels, however, are highly variable, which results in an uneven read coverage distribution and a variable error rate. Also, genome assembly usually is aimed at generating one consensus sequence, whereas in transcriptome assembly, alternative splicing can lead to multiple consensus sequences for one gene. Further challenges are introduced when both classical and next-generation sequencing technologies are employed, because these differ in sequencing coverage, read length and accuracy.
We addressed these problems by applying a two-step strategy. First, Sanger and 454/Roche reads were assembled using PAVE, a classical overlap-layout-consensus assembly tool. PAVE compensates well for sequencing errors and less perfect overlaps and yields longer contigs within a reasonable runtime. Second, Solexa/Illumina reads were iteratively assembled onto a backbone of PAVE contigs using a de Brujin graph assembly tool, which is better suited for large amounts of short reads and much faster than PAVE. Two additional steps were of importance for the iterative assembly procedure. One, input complexity of each next-generation library can be largely reduced by mapping reads onto existing contigs and discarding those reads that fully align to a contig. Two, since de Brujin graph assembly tools are more susceptible by sequencing errors and thus tend to assemble shorter contigs, a classical overlap-layout-consensus assembly tool was used to merge these contigs into larger ones.
Quality assessment of the transcriptome assembly
Quality control of assemblies is essential because the quality of the assembled contigs has a direct impact on subsequent analyses. However, while a number of established metrics and protocols exist for genome assembly, only a few have been recently proposed for transcriptome assembly [29, 30]. Moreover, they were originally defined for reference-based transcriptome assemblies and require a genome sequence to compare against; although, technically, a set of already assembled transcript sequences might also serve as reference. This is often not the case for de novo transcriptome assemblies and makes quality assessment much more difficult. In case of N. furzeri, only very few sequences were available, and their use as reference was prevented by their small number. Instead, we used a set of protein sequences with experimental support from medaka, the closest relative of N. furzeri with a genome available.
The majority of these medaka proteins (74%) were hit by at least one N. furzeri transcript contig, which suggested a good representation of these proteins in the transcript catalogue. Moreover, most proteins were completely covered by N. furzeri transcript contigs and, notably, a large number of proteins were completely represented by a single contig. Note that these results were obtained from a cross-species comparison, and the sequence divergence between the two species might interfere with the results, thus representing a conservative estimate.
The remaining proteins without hits might be missed for several reasons. First of all, the corresponding genes might not be expressed in the sequenced transcriptomes. Alternatively, the sequence divergence between the N. furzeri and medaka proteins might be too high so that BLASTx failed to identify these proteins. Finally, their genes might not be present in N. furzeri or were turned into pseudogenes without expression.
We also assessed the fragmentation and redundancy of the transcript catalogue. The fragmentation was relatively low, and only a few transcript contigs were regarded as fragments (non-overlapping). In contrast, the vast majority of the medaka protein sequences were covered by multiple contigs, indicative for alternative splicing and/or gene duplication events.
Chimeric transcript contigs may originate from gene or transcript fusions and assembly errors. The observed low putative-chimerism rate of 0.3% was regarded as an indication for a negligible misassembly rate not requiring further computational and / or experimental efforts. More insights into the falsity or validity of these putative chimeric transcript contigs will provide an N. furzeri genome sequence in the future.
Completeness of the transcript catalogue
The presented N. furzeri transcript catalogue represents almost 20,000 protein-coding genes, the number typically expected for a vertebrate genome . However, without a high-quality reference genome sequence, it is not possible to determine the true number of N. furzeri genes. Therefore, we first compared the gene number to respective numbers in fish species with a sequenced genome, i.e. medaka, stickleback, tetraodon and zebrafish. Hits to annotated genes indicated that between 71% and 87% of N. furzeri protein-coding genes were identified by our approach. But comparing gene numbers between species heavily depends on the quality of annotations. Therefore, we determined gene families since their number should be more constant among fish. The results of this analysis led to the conclusion that the presented N. furzeri transcript catalogue is in terms of completeness comparable to gene catalogues of model fish for which genomes have been sequenced.
The genes missing in our N. furzeri transcript catalogue are likely those expressed in early development of N. furzeri, because we used specimens of at least one week of age and were not able to include embryonic tissue. For example, in zebrafish, previous studies isolated several hundred genes involved in developmental control [32, 33]. Additionally, gene annotation in zebrafish is currently far more comprehensive than in the other three fish species medaka, stickleback and tetraodon, and is also largely supported by expression data. Thus, it is most likely that a number of genes remain to be identified also in the other species. In addition, the number of gene copies retained after the fish-specific genome duplication varies considerably between the four species . Most likely, this is also the case for N. furzeri. Thus, we cannot exclude that some of the missing genes are actually lost in N. furzeri.
In addition to overall gene number, the representation of CDS is important. For more than 70% of the N. furzeri protein-coding genes, a transcript contig with a complete (> 90%) CDS is present. When the sum of the non-redundant CDS identified in the transcript catalogue was compared to the total CDS of all reference entries used for annotation, completeness was estimated to be 83%. Future analyses on alternative UTRs and splice isoforms will greatly benefit from the transcript catalogue once the N. furzeri genome sequence is available.
Duplicated genes in N. furzeri
We found evidences for multiple copies of about 19% of the N. furzeri protein-coding genes, whereas the vast majority occurs in paralogous pairs reminiscent to the 3rd genome duplication in teleost fish [24, 25]. Conversely, this indicates that >80% of gene copies originated by this genome duplication were lost from the N. furzeri genome which corresponds to previous estimates for teleost fishes of 76–85% . Moreover, it has to be considered, that not all gene duplications may date back to the teleost-specific genome duplication and that by our transcriptome approach only paralogs could be detected whose expression was sufficiently high in the analysed samples. Once the N. furzeri genome sequence is available, a more comprehensive assessment of duplicated genes and their evolutionary fate will be possible.
Detection of ageing-related, differentially expressed genes in N. furzeri
Selected RNA-seq experiments initially aimed at completing the transcript catalogue were subsequently further exploited to detect changes in N. furzeri gene expression over lifetime. For mapping and quantification of the respective datasets, only the longest transcript contig of each gene was used (19,875 contigs, 53 Mb). Respectively, only 53% of reads mapped uniquely to this reference. When using all transcript contigs with gene annotations as reference sequences (85,431 contigs, 144 Mb), mapping efficiency increased to 91%, thus indicating that initial mapping results are not due to incompleteness of the transcript catalogue. On the other hand, the latter results included a high fraction of multiple mappings (~80%). This redundancy of the transcript catalogue is caused both by the assembly of contigs representing alternative splicing events and gene paralogs as well by assembly problems due to sequencing errors and repetitive motifs. Solutions to deal with multiple mappings were described previously [35–38], but most are designed for genome reference sequences and presumably would introduce additional bias into our transcriptome data and thus were not applied. However, we did compare mapping efficiency using the zebrafish RNA-seq data. Interestingly, when we used transcript instead of genome sequences as a reference for mapping the zebrafish data, mapping percentages were between 40% and 49% vs. up to 80% for the genome (data not shown). This suggested that the number of mapped reads is largely determined by the respective reference, and using genome sequences as reference is superior to using all (known) transcripts. Again, it is highly desirable to generate an N. furzeri genome reference to shed additional light on this question.
Despite these limitations, we quantified transcript levels in skin and brain of young and old GRZ and MZM-0403. In all samples, over 17,000 genes (Table 5) could be profiled regardless of tissue. This was surprising, because transcriptome profiles are generally believed to be specific for the analysed tissue or organ [39–41]; however, many genes had very low transcript levels, suggesting that RNA-seq detected some kind of baseline transcription. When analysing the transcriptome profiles for general correlation, samples clearly clustered according to tissues. Interestingly, additional weaker but tissue-specific correlations were observed. For skin, samples of the same age clustered next to each other, indicating that age-specific transcript patterns are similar in both short- and longer-lived strains. In contrast, in brain samples of the same strain clustered next, suggesting that in this tissue transcriptome profiles differ more between strains than between ages. Finally, PCA conducted on all samples showed that most of the variation between the N. furzeri transcriptome profiles can be explained by age (reflected by PC1 which explains 56% of the variation), strain (PC2, 22%) and tissue (PC3, 18%). Notably, based on correlation, samples of the same tissue are most similar and, consequently, PCA identified tissue as the component that introduces the least amount of variation.
We identified 43 DEGs in N. furzeri skin and brain, respectively. Interestingly, these tissues differ in the major direction of transcript level changes with age. In skin, the majority of DEGs were down regulated whereas in brain the opposite was found. This observation is in good agreement with both correlation analysis and PCA and again points to tissue-specific signatures of ageing. Further, DEGs functions reflect expected signs of an ageing phenotype . DEGs involved in apoptosis were upregulated whereas DEGs with a role in cell cycle control, cell division and proliferation were down regulated. This was, however, more prominent in brain than in skin, suggesting more dynamic/pronounced changes of transcript levels in brain over lifetime compared to skin. Transcript levels also showed an increase in inflammation in brain, which previously has been reported as another sign of age-related decline in brain [43, 44]. In skin, a number of downregulated DEGs were involved in the maintenance of collagenous tissue and cartilage.
Unfortunately, the RNA-seq datasets for N. furzeri did lack replicates, as they primarily were designed to complete the transcript catalogue, i.e. represent pooled samples of 3–5 specimens. On the other hand, those pooled datasets are less prone to fluctuations. The validation of selected DEGs by qRT-PCR showed that the rate of putative false positives was low (2 of 12). As in both cases the same qRT-PCR assay failed to reach significance (although the direction of expression change was confirmed), a technical PCR problem cannot be ruled out. However, the present approach using only two datasets derived from different strains as biological replicates most likely missed a considerable number of DEGs that show higher variation within or between strains and therefore failed to reach significance.
Using RNA-seq data from skin in young and old zebrafish we were able to confirm that transcript patterns of ageing are similar in both species. Thus, N. furzeri, in spite of its exceptionally short lifespan of only few months, exhibits general characteristics of vertebrate ageing which is in line with previously observed typical cytological and physiological ageing-related phenotypes .
Accelerated ageing in N. furzeristrain GRZ
Previous work based on the analysis of age-related biomarkers suggested that ageing is accelerated in the short-lived strain GRZ compared to the longer-lived strain MZM-0403 [7, 45]. Our analysis of whole transcriptomes of GRZ and MZM-0403 also supported the difference in ageing pace between these strains. When comparing time-normalised fold changes of identified DEGs, we found that age-related transcript level changes are stronger in GRZ than in MZM-0403. As the assumption of linear changes in transcript levels with age is most likely in many cases an oversimplification, the presented two-point analysis is a first approximation. Further studies should include additional time points to better dissect chronological and biological ageing in the N. furzeri strains.
Future work will be directed towards completing the transcript catalogue, e.g. including genes expressed in early development. In parallel, we are working on the assembly and annotation of a first working draft of the N. furzeri genome sequence which will be available in the near future. Transcriptome and genome sequences will be used to determine gene structures, detect splice isoforms and genes encoding non-protein-coding RNA species. Second, further RNA-seq experiments will be performed utilising more and strain-specific biological replicates as well as additional time points to better understand ageing-related mechanisms in N. furzeri and vertebrates in general.
We established a high-quality, comprehensive and annotated transcript catalogue of the short-lived fish species N. furzeri, a novel model organism for age research. We employed it as a mapping reference for RNA-seq experiments and identified genes which transcript levels were significantly changed upon ageing. Regulation was tissue-specific with the majority of genes upregulated in brain but downregulated in skin. A comparison to zebrafish suggested that transcript patterns of ageing are similar in both species. Further, ageing-related changes were accelerated in the short-lived N. furzeri strain with respect to the longer-lived.
Fish strains and maintenance
N. furzeri strains GRZ and MZM-0403 were described previously [7, 9]. Up to 12 fish per strain were kept in 40 litre tanks at 26°C under a light regime of 12:12 h light:dark and fed on red mosquito larvae (Chironomidae) ad libitum once a day. Water was filtered using air-driven foam filters; the water was changed once a week. Wild-type zebrafish of the TüAB strain were kept in groups in an open circulating standard zebrafish system. All experiments were performed according to the “Principles of laboratory animal care” as well as to the current version of the German Law on the Protection of Animals.
RNA isolation, library construction and sequencing
For N. furzeri, total RNA was isolated from skin and brain using the RNeasy Mini kit (Qiagen) and from whole body using the RNeasy Midi kit (Qiagen). Males and females were used as given in Table 1. For zebrafish, total RNA was isolated from skin of male specimens aged 5 and 42 months, respectively, using Trizol (Invitrogen).
For Sanger sequencing (library 1, Table 1), cDNA was normalised by Evrogen (Moscow, Russia) and cloned into pAL17.3 via SfiI. Recombinant plasmids were amplified in E. coli, purified, sequenced from both ends (BigDye Terminator v.2.1 Cycle Sequencing Kit, ABI, Weiterstadt, Germany) and separated on ABI 3730xl capillary sequencers.
For 454/Roche sequencing, normalised and non-normalised cDNA was prepared by Evrogen, Moscow, Russia (library 2 and 3, Table 1) and Vertis, Freising, Germany (library 4, Table 1). The Roche adaptor containing cDNA libraries were diluted and subjected to emPCR followed by sequencing according to the manufacturer’s instructions (Roche Diagnostics). Sequencing was performed on 70 × 75 picotiter plates by two and a half runs using the GS LR70 (FLX) Sequencing kit and two runs with the XLR70t (Titanium) kit, respectively.
For Solexa/Illumina sequencing, N. furzeri libraries 5 to 13 and zebrafish libraries were prepared using the mRNA-Seq 8 sample prep kit (RS-100-0801) (Illumina, San Diego, CA, USA) according to the manufacturer’s instruction. Clusters were generated using Illumina’s Single Read Cluster Generation Kit v4 (GD-103-4001) or Paired End Cluster Generation Kit v4 (PE-203-4001), respectively. Each library was loaded onto one lane of the flow cell at 7 pM concentration. Sequencing of N. furzeri and zebrafish libraries was performed on a Genome Analyzer IIx (Illumina) for either 76, 101 or 150 cycles (Table 1), using Illumina’s 36 Cycle Sequencing Kit v4 (FC-104-4002) following the manufacturer’s protocol.
Sanger reads were deposited at the NCBI dbEST database. 454/Roche and Solexa/Illumina reads were submitted to the NCBI sequence read archive (SRA) under submission accession SRA050046. A NCBI BioProject, accession PRJNA85613, provides links to all submitted sequence data. Zebrafish Solexa/Illumina reads were submitted to the NCBI SRA under submission accession SRA054207 (NCBI BioProject PRJNA169239).
Sanger and 454/Roche reads with associated quality values were pre-processed as follows: Removal of known vector sequences and poly (A) tails as well as masking of low complexity regions was done with SeqClean (v. 1.0, ). Repetitive elements were masked with RepeatMasker (v. 3.2.9, ) and a library of N. furzeri-specific repeats . For, Solexa/Illumina libraries, the String Graph Assembler (v. 0.9.12, ) was used to correct sequencing errors and, for assembly, to remove low-quality reads.
TGICL (v. 2.0), Newbler (v. 2.3) and PAVE (v. 2.0) were run with default parameters. Only contigs longer than or equal to 300 bp were further analysed. Solexa/Illumina libraries were assembled using the following iterative assembly procedure on each library: In each iteration, the Burrows-Wheeler Alignment tool (BWA v. 0.5.9, ) was used to map Solexa/Illumina reads of the respective library onto the previous assembly. Subsequently, remaining unmapped reads and contigs of the previous assembly were joined with the CLC Assembly Cell (v. 3.2.2, ), a de Bruijn graph assembly tool. To improve length and quality, the resulting contigs were further assembled with TGICL. To use large contigs (≥ 2,000 bp) as input for TGICL, they were split into 2,000 bp-fragments which overlap by 300 bp. After assembly, CD-HIT-EST (v. 4.5.5, ) was used to remove highly redundant contigs (99% identity and 90% coverage of the smaller sequence). Further, all contigs smaller than 300 bp were excluded. Each round resulted in a “new” assembly consisting of the previous one plus new sequence information introduced by the respective assembled Solexa/Illumina dataset and served as a backbone for the next assembly utilising the reads of a further library. Once this iterative assembly process was finished, N. furzeri GRZ Solexa/Illumina libraries were mapped with BWA onto previously generated contigs. The final GRZ-specific consensus sequence was called with SAMtools and BCFtools (v. 0.1.16, ) thereby improving the overall contig quality and correcting base calls introduced by MZM-0403 reads.
All transcript contigs were deposited at the NCBI transcriptome shotgun assembly (TSA) archive. The corresponding accessions are also available at the NCBI BioProject PRJNA85613 page.
Annotation was based on BLAST searches against several protein and nucleotide databases (Table 3). BLAST parameters and e-value thresholds were defined for each database separately (Additional file 1: Table S1). Hits were filtered for maximum e-value, minimal overlap and minimum identity. Hits with non-informative descriptions like ‘bac’, ‘fosmid’ or ‘shotgun’ were excluded. Additional information like gene symbols or species names were parsed from definition lines of hits or, in case of Ensembl databases, fetched from the Ensembl MySQL server. Transcript contigs were annotated based on the best hits of each database and description, gene symbol and species name were copied as additional information.
Gene Ontology files (OBO v1.2) and database associations for UniProt were downloaded from the Gene Ontology web site  and stored in a MySQL database. GO terms were mapped to GO-Slim terms using the map2slim script (go-perl package, see Gene Ontology web site). For each transcript, UniProt identifier were parsed from its annotations and used to assign GO terms. For transcript contigs annotated with hits from Ensembl databases, additional GO terms were fetched from the Ensembl MySQL server.
Identification of duplicated genes
To identify genes which are duplicated in teleost fish, the longest transcript contig per gene was compared using BLASTx (e-value ≤ 10-07) against Ensembl proteins of medaka, stickleback, tetraodon and zebrafish. Subsequently, information about paralogous genes in these fish species was obtained from the Ensembl Compara gene trees. In case, the two BLASTx hits point to Ensembl genes of two different species, the Ensembl gene of the first species is ‘replaced’ with the corresponding ortholog in the second species.
To identify genes exclusively duplicated in N. furzeri, all transcript contigs which belong to the same gene and share the same annotation were retrieved and translated into proteins. Only protein sequences with a minimum length of 100 aa were considered for further analysis. The protein sequence of the longest transcript contig was then aligned to each of the other (secondary) transcript contigs thereby generating pairwise protein alignments. Alignments were filtered for a minimum length of 100 aa and pairwise protein identities were calculated. The remaining (secondary) transcript contigs were clustered based on their protein divergence to the longest transcript contig with a minimum distance of 10% between each cluster. Two or more clusters indicate genes putatively duplicated exclusively in N. furzeri.
Gene expression profiling and qRT-PCR validation
To quantify transcript levels in eight N. furzeri samples, 76 bp-reads from libraries 6–13 were mapped with BWA onto a nonredundant set of transcript contigs representing the longest transcript contig per gene. Zebrafish datasets were mapped with Bowtie  against the zebrafish genome (danRer7) and a precompiled exon junction database.
Only unique mappings were counted. For each gene, coverage was determined as number of reads mapping onto the respective transcript contig normalised to its length in kilobases and to the total number of mapped reads in millions (rpkm). Analyses of correlation (Spearman) and principal component analysis were done with R. Identification of DEGs was done with R and the DESeq package. Raw read counts were provided as input for DESeq.
qRT-PCR was carried out in MZM-0403. RNA of the same skin and brain samples analysed in the RNA-seq experiment was used. cDNA synthesis was performed in a 20 μl volume using 500 ng total RNA, 10 pmol Oligo-dT primer and 200 U SuperScript II reverse transcriptase (Invitrogen). Real-time PCR was performed with the SYBR GreenER qPCR SuperMix (Invitrogen) and using the iCycler iQ5 detection system (Bio-Rad). Primer sequences are given in Additional file 8: Table S2. Cycle threshold values were normalised to the gene INSR (insulin receptor) that already previously showed to be very stably expressed at different ages . Fold changes were determined with the relative expression software tool .
We thank Tom Hofmann, Ivonne Görlich, Kathleen Seitz, Oliver Backhaus and Hella Ludewig for expert technical assistance, and Alessandro Cellerino for providing the N. furzeri strain MZM-0403. We are especially grateful to Karol Szafranski for support in statistical analysis. Moreover, we thank Katrin Sameith (University Medical Centre Utrecht) for useful comments on the manuscript. This work was supported by the Leibniz-Association (SAW2006) and the Jena Centre for Systems Biology of Ageing – JenAge, funded by the German Ministry for Education and Research (BMBF; support code: 0315581).
- Christensen K, Johnson TE, Vaupel JW: The quest for genetic determinants of human longevity: challenges and insights. Nat Rev Genet. 2006, 7: 436-448. 10.1038/nrg1871.PubMed CentralView ArticlePubMedGoogle Scholar
- Curran SP, Ruvkun G: Lifespan Regulation by Evolutionarily Conserved Genes Essential for Viability. PLoS Genet. 2007, 3: e56-10.1371/journal.pgen.0030056.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith ED, Tsuchiya M, Fox LA, Dang N, Hu D, Kerr EO, Johnston ED, Tchao BN, Pak DN, Welton KL, Promislow DEL, Thomas JH, Kaeberlein M, Kennedy BK: Quantitative evidence for conserved longevity pathways between divergent eukaryotic species. Genome Res. 2008, 18: 564-570. 10.1101/gr.074724.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Piper MDW, Selman C, McElwee JJ, Partridge L: Separating cause from effect: how does insulin/IGF signalling control lifespan in worms, flies and mice?. J Intern Med. 2008, 263: 179-191. 10.1111/j.1365-2796.2007.01906.x.View ArticlePubMedGoogle Scholar
- Genade T, Benedetti M, Terzibasi E, Roncaglia P, Valenzano DR, Cattaneo A, Cellerino A: Annual fishes of the genus Nothobranchius as a model system for aging research. Aging Cell. 2005, 4: 223-233. 10.1111/j.1474-9726.2005.00165.x.View ArticlePubMedGoogle Scholar
- Jubb RA: A new Nothobranchius (Pisces, Cyprinodontidae) from Southeastern Rhodesia. J Am Killifish Assoc. 1971, 8: 314-321.Google Scholar
- Terzibasi E, Valenzano DR, Benedetti M, Roncaglia P, Cattaneo A, Domenici L, Cellerino A: Large Differences in Aging Phenotype between Strains of the Short-Lived Annual Fish Nothobranchius furzeri. PLoS One. 2008, 3 (12): e3866-10.1371/journal.pone.0003866.PubMed CentralView ArticlePubMedGoogle Scholar
- Valdesalici S, Cellerino A: Extremely short lifespan in the annual fish Nothobranchius furzeri. Proc. Biol. Sci. 2003, 270 (Suppl 2): S189-S191.PubMed CentralView ArticlePubMedGoogle Scholar
- Reichwald K, Lauber C, Nanda I, Kirschner J, Hartmann N, Schories S, Gausmann U, Taudien S, Schilhabel M, Szafranski K, Glöckner G, Schmid M, Cellerino A, Schartl M, Englert C, Platzer M: High tandem repeat content in the genome of the short-lived annual fish Nothobranchius furzeri: a new vertebrate model for aging research. Genome Biol. 2009, 10: R16-10.1186/gb-2009-10-2-r16.PubMed CentralView ArticlePubMedGoogle Scholar
- Valenzano DR, Terzibasi E, Cattaneo A, Domenici L, Cellerino A: Temperature affects longevity and age-related locomotor and cognitive decay in the short-lived fish Nothobranchius furzeri. Aging Cell. 2006, 5: 275-278. 10.1111/j.1474-9726.2006.00212.x.View ArticlePubMedGoogle Scholar
- Kirschner J, Weber D, Neuschl C, Franke A, Böttger M, Zielke L, Powalsky E, Groth M, Shagin D, Petzold A, Hartmann N, Englert C, Brockmann GA, Platzer M, Cellerino A, Reichwald K: Mapping of quantitative trait loci controlling lifespan in the short-lived fish Nothobranchius furzeri- a new vertebrate model for age research. Aging Cell. 2012, 11 (2): 252-261. 10.1111/j.1474-9726.2011.00780.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Myers EW: Toward simplifying and accurately formulating fragment assembly. J Comput Biol. 1995, 2: 275-290. 10.1089/cmb.1995.2.275.View ArticlePubMedGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.View ArticlePubMedGoogle Scholar
- 454 Life Sciences, a Roche Company. http://www.454.com/,
- Soderlund C, Johnson E, Bomhoff M, Descour A: PAVE: program for assembling and viewing ESTs. BMC Genomics. 2009, 10: 400-10.1186/1471-2164-10-400.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralView ArticlePubMedGoogle Scholar
- Flicek P, Amode MR, Barrell D, Beal K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fairley S, Fitzgerald S, Gil L, Gordon L, Hendrix M, Hourlier T, Johnson N, Kahari AK, Keefe D, Keenan S, Kinsella R, Komorowska M, Koscielny G, Kulesha E, Larsson P, Longden I, McLaren W, Muffato M, Overduin B, Pignatelli M, Pritchard B, Riat HS, Ritchie GRS, Ruffier M, Schuster M, Sobral D, Tang YA, Taylor K, Trevanion S, Vandrovcova J, White S, Wilson M, Wilder SP, Aken BL, Birney E, Cunningham F, Dunham I, Durbin R, Fernandez-Suarez XM, Harrow J, Herrero J, Hubbard TJP, Parker A, Proctor G, Spudich G, Vogel J, Yates A, Zadissa A, Searle SMJ: Ensembl 2012. Nucleic Acids Res. 2011, 40: D84-D90.PubMed CentralView ArticlePubMedGoogle Scholar
- Pruitt KD, Tatusova T, Maglott DR: NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005, 33: D501-D504. 10.1093/nar/gki476.PubMed CentralView ArticlePubMedGoogle Scholar
- Boeckmann B, Bairoch A, Apweiler R, Blatter M-C, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O’Donovan C, Phan I, Pilbout S, Schneider M: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31: 365-370. 10.1093/nar/gkg095.PubMed CentralView ArticlePubMedGoogle Scholar
- National Center for Biotechnology. http://www.ncbi.nlm.nih.gov/,
- Pontius J, Wagner L, Schuler G: UniGene: a unified view of the transcriptome. The NCBI Handbook. 2003, Bethesda (MD): National Center for Biotechnology Information, Available: http://www.ncbi.nlm.nih.gov/books/NBK21101/Google Scholar
- The Gene Ontology. http://www.geneontology.org/,
- Vilella AJ, Severin J, Ureta-Vidal A, Heng L, Durbin R, Birney E: EnsemblCompara GeneTrees: Complete, Duplication-Aware Phylogenetic Trees in Vertebrates. Genome Res. 2009, 19: 327-335.PubMed CentralView ArticlePubMedGoogle Scholar
- Amores A, Force A, Yan YL, Joly L, Amemiya C, Fritz A, Ho RK, Langeland J, Prince V, Wang YL, Westerfield M, Ekker M, Postlethwait JH: Zebrafish hox clusters and vertebrate genome evolution. Science. 1998, 282: 1711-1714.View ArticlePubMedGoogle Scholar
- Taylor JS, Braasch I, Frickey T, Meyer A, Van de Peer Y: Genome duplication, a trait shared by 22000 species of ray-finned fish. Genome Res. 2003, 13: 382-390. 10.1101/gr.640303.PubMed CentralView ArticlePubMedGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMed CentralView ArticlePubMedGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticleGoogle Scholar
- Martin JA, Wang Z: Next-generation transcriptome assembly. Nat Rev Gen. 2011, 12: 671-682. 10.1038/nrg3068.View ArticleGoogle Scholar
- Martin J, Bruno VM, Fang Z, Meng X, Blow M, Zhang T, Sherlock G, Snyder M, Wang Z: Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads. BMC Genomics. 2010, 11: 663-10.1186/1471-2164-11-663.PubMed CentralView ArticlePubMedGoogle Scholar
- Volff J-N: Vertebrate genomes. Basel, New York: KargerGoogle Scholar
- Geisler R, Rauch G-J, Geiger-Rudolph S, Albrecht A, Van Bebber F, Berger A, Busch-Nentwich E, Dahm R, Dekens MPS, Dooley C, Elli AF, Gehring I, Geiger H, Geisler M, Glaser S, Holley S, Huber M, Kerr A, Kirn A, Knirsch M, Konantz M, Küchler AM, Maderspacher F, Neuhauss SC, Nicolson T, Ober EA, Praeg E, Ray R, Rentzsch B, Rick JM, Rief E, Schauerte HE, Schepp CP, Schönberger U, Schonthaler HB, Seiler C, Sidi S, Söllner C, Wehner A, Weiler C, Nüsslein-Volhard C: Large-scale mapping of mutations affecting zebrafish development. BMC Genomics. 2007, 8: 11-10.1186/1471-2164-8-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Haffter P, Granato M, Brand M, Mullins MC, Hammerschmidt M, Kane DA, Odenthal J, Van Eeden FJ, Jiang YJ, Heisenberg CP, Kelsh RN, Furutani-Seiki M, Vogelsang E, Beuchle D, Schach U, Fabian C, Nüsslein-Volhard C: The identification of genes with unique and essential functions in the development of the zebrafish, Danio rerio. Development. 1996, 123: 1-36.PubMedGoogle Scholar
- Brunet FG, Roest Crollius H, Paris M, Aury J-M, Gibert P, Jaillon O, Laudet V, Robinson-Rechavi M: Gene loss and evolutionary rates following whole-genome duplication in teleost fishes. Mol Biol Evol. 2006, 23: 1808-1816. 10.1093/molbev/msl049.View ArticlePubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Kwan G, Van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28: 511-515. 10.1038/nbt.1621.PubMed CentralView ArticlePubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, Di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29: 644-652. 10.1038/nbt.1883.PubMed CentralView ArticlePubMedGoogle Scholar
- Robertson G, Schein J, Chiu R, Corbett R, Field M, Jackman SD, Mungall K, Lee S, Okada HM, Qian JQ, Griffith M, Raymond A, Thiessen N, Cezard T, Butterfield YS, Newsome R, Chan SK, She R, Varhol R, Kamoh B, Prabhu A-L, Tam A, Zhao Y, Moore RA, Hirst M, Marra MA, Jones SJM, Hoodless PA, Birol I: De novo assembly and analysis of RNA-seq data. Nat Methods. 2010, 7: 909-912. 10.1038/nmeth.1517.View ArticlePubMedGoogle Scholar
- Li B, Dewey CN: RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011, 12: 323-10.1186/1471-2105-12-323.PubMed CentralView ArticlePubMedGoogle Scholar
- Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. PNAS. 2004, 101: 6062-6067. 10.1073/pnas.0400782101.PubMed CentralView ArticlePubMedGoogle Scholar
- Whitehead A, Crawford DL: Variation in tissue-specific gene expression among natural populations. Genome Biol. 2005, 6: R13-10.1186/gb-2005-6-2-r13.PubMed CentralView ArticlePubMedGoogle Scholar
- Enard W, Khaitovich P, Klose J, Zöllner S, Heissig F, Giavalisco P, Nieselt-Struwe K, Muchmore E, Varki A, Ravid R, Doxiadis GM, Bontrop RE, Pääbo S: Intra- and Interspecific Variation in Primate Gene Expression Patterns. Science. 2002, 296: 340-343. 10.1126/science.1068996.View ArticlePubMedGoogle Scholar
- De Magalhães JP, Curado J, Church GM: Meta-analysis of age-related gene expression profiles identifies common signatures of aging. Bioinformatics. 2009, 25: 875-881. 10.1093/bioinformatics/btp073.PubMed CentralView ArticlePubMedGoogle Scholar
- Schumacher B, Van der Pluijm I, Moorhouse MJ, Kosteas T, Robinson AR, Suh Y, Breit TM, Van Steeg H, Niedernhofer LJ, Van Ijcken W, Bartke A, Spindler SR, Hoeijmakers JHJ, Van der Horst GTJ, Garinis GA: Delayed and accelerated aging share common longevity assurance mechanisms. PLoS Genet. 2008, 4: e1000161-10.1371/journal.pgen.1000161.PubMed CentralView ArticlePubMedGoogle Scholar
- Park S-K, Kim K, Page GP, Allison DB, Weindruch R, Prolla TA: Gene expression profiling of aging in multiple mouse strains: identification of aging biomarkers and impact of dietary antioxidants. Aging Cell. 2009, 8: 484-495. 10.1111/j.1474-9726.2009.00496.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Di Cicco E, Tozzini ET, Rossi G, Cellerino A: The short-lived annual fish Nothobranchius furzeri shows a typical teleost aging process reinforced by high incidence of age-dependent neoplasias. Exp Gerontol. 2011, 46: 249-256. 10.1016/j.exger.2010.10.011.View ArticlePubMedGoogle Scholar
- TGI Software Tools - SeqClean. http://compbio.dfci.harvard.edu/tgi/software/,
- RepeatMasker Open-3.0. http://www.repeatmasker.org/,
- Simpson JT, Durbin R: Efficient construction of an assembly string graph using the FM-index. Bioinformatics. 2010, 26: i367-i373. 10.1093/bioinformatics/btq217.PubMed CentralView ArticlePubMedGoogle Scholar
- CLC bio. http://www.clcbio.com/,
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.View ArticlePubMedGoogle Scholar
- SAM - next generation alignment tools. http://samtools.sourceforge.net/,
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Hartmann N, Reichwald K, Lechel A, Graf M, Kirschner J, Dorn A, Terzibasi E, Wellner J, Platzer M, Rudolph KL, Cellerino A, Englert C: Telomeres shorten while Tert expression increases during ageing of the short-lived fish Nothobranchius furzeri. Mech Ageing Dev. 2009, 130: 290-296. 10.1016/j.mad.2009.01.003.View ArticlePubMedGoogle Scholar
- Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.PubMed CentralView ArticlePubMedGoogle Scholar
- Korf I, Yandell M, Bedell J: BLAST. Sebastopol, CA, USA: O’Reilly & Associates, IncGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.