BACKGROUND
teleosts, which roughly constitute 96% of all living fishes and half of the extant vertebrate, are the most phenotypically diversified and species-rich group of all the vertebrate species  <cit> .
the vast morphological and species diversity of teleosts have received intense attention worldwide because of their importance in both scientific research and aquaculture.
however, the genomic basis of the complex phenotype of teleosts during evolution remains obscure, despite the increasing amount of genome and transcriptome sequence data available.
one important mechanism for the phenotypic diversity of species is the duplication of genes and entire genomes  <cit> .
evidence has recently been accumulated to allow a consensus to be reached that all teleosts experienced an additional whole genome duplication , which occurred after the basal ray-finned fishes separate from the actinopterygian stem lineage but before the teleosts began radiation  <cit> .
combining the absolute dates and phylogenetic timing of the 3r duplication, some groups have thought that the fsgd might be causally related to an increase in the number of species as well as their biological diversity  <cit> .
however, if the fsgd  was responsible for the evolutionary success and astounding biological diversification of teleosts, it must have occurred prior to the radiation of teleosts.
with the fossil record, paleontological evidence have suggested that the first appearance of most of the extant teleosts was only about  <dig> million years ago  <cit> , which is shorter than the fsgd that occurred at least  <dig> million years ago.
thus, it would be inappropriate to think that the fsgd was a major driving force behind the rapid radiation of teleosts  <cit> .
furthermore, considering the large amount of fossil data for basal ray-finned fishes, a consequential fsgd would not provide sufficient explanation for the morphological complexity of teleosts  <cit> .
besides the fish-specific genome duplication, alternative explanations for the increasing morphological complexity of teleosts include their experience with a higher rate of chromosomal rearrangements  <cit>  and a faster evolution of protein sequences  <cit>  and conserved noncoding elements   <cit>  compared to cartilaginous fishes and mammals.
their implications for the evolution and diversity of teleosts have been intensively discussed  <cit> .
however, conserved teleost lineage-specific genes have been poorly characterized.
lineage-specific genes, also referred to as taxonomically restricted genes   <cit>  are defined as genes found in one particular taxonomic group but share no sequence similarity with genes from other lineages  <cit> .
with the advent of large-scale genome sequencing projects for a wide range of species, lineage-specific genes have been extensively studied in mammals  <cit> , insects  <cit> , plants  <cit> , and microbial species  <cit> .
lineage-specific genes are a significantly abundant component of all genomes sequenced to-date  <cit> , which defies an early hypothesis that an increasing database size would eventually reduce the number of lineage-specific genes  <cit> .
orphan genes were first discussed when analyzing the yeast genome; approximately one-third of the identified genes fell into this category  <cit> .
likewise in drosophila melanogaster, the most accurate and complete genome analyzed, lineage-specific genes were found to make up nearly  <dig> % of the total genes  <cit> .
apart from being abundant, lineage-specific genes have also been thought to be important for lineage specific traits and adaptations  <cit> .
in hydra, for example, interspecific differences in tentacle formation are closely related to the changes in the expression of taxonomically restricted genes  <cit> .
and in drosophila, the flightin gene is specifically important for increasing the frequency of the flight muscle to deliver the maximum power to the wing, which is a rather specific adaptation for the dipterans  <cit> .
although abundant in quantity and important in functionality, the evolutionary origin of lineage-specific genes is still enigmatic.
several hypotheses about the origin of lineage-specific genes, including gene duplication followed by rapid sequence divergence, lateral gene transfer, accelerated evolutionary rate, artifacts from genome annotation, as well as de novo evolution from noncoding sequences have been proposed  <cit> .
despite the fact that the origin and evolution of lineage-specific genes is still poorly understood, the identification, characterization, function, and expression analysis of lineage-specific genes may provide a better understand for lineage-specific adaptation, such as the successful diversification of teleosts.
in this study, we identified conserved teleost-specific genes  and orphan genes in zebrafish using comparative genomics.
we then characterized each set of these genes by diverse features, including gene size, protein size, exon number, gc content, transcript support, and chromosomal locations.
as a large portion of the ctsgs and orphan genes have no known function, ab initio predictions using protfun were performed to infer possible biological functions.
we then explored the evolutionary origin of lineage-specific genes and performed a comprehensive analysis of their tissues and developmental stages specific expression patterns using the wealth of available expression data, including est and rna-seq data, which in turn provided important complementary datasets that may be used to uncover their functions in the future.
collectively, identification of lineage-specific genes as well as orphan genes and future studies of their function by means of target gene knockdown  <cit>  or knockout  <cit>  will no doubt help increase our understanding of the molecular basis of the successful diversification of teleosts.
RESULTS
identification of ctsgs and orphan genes
the procedure used to identify orphan genes and ctsgs was modified from previous studies  <cit> .
we used the blast-basic local alignment search tool  <cit>  with default parameters for all assignments.
firstly, we used blastp to search for homologs of all  <dig>  proteins annotated of zebrafish in each of the other  <dig> species-levels outside the teleostei protein sets with an e-value cutoff of 10- <dig>  those proteins for which we could not find any homolog in any alternative-spliced forms were used for the next set of searches, which was done using tblastn to search for homologs against the genome of the same species.
in this step of the search, we only performed the tblastn search when there were no annotated proteins of the genome of the species, such as the genomes of elephant shark and salmo salar.
a total of  <dig>  danio rerio proteins with significant sequence similarity  to either a protein or genomic sequence from a species outside the group of teleosts were defined as the evolutionarily conserved proteins .
secondly, the remaining  <dig>  d. rerio proteins with no significant similarity to any sequence  outside the group of teleosts were further searched against the proteins and the genome of the group of teleosts, including oryzias latipes, gasterosteus aculeatus, takifugu rubripes, tetraodon nigroviridis, oreochromis niloticus, gadus morhua, salmo salar, using blastp and tblastn with e-values the same as the first step.
this step produced two datasets: 1)  <dig> proteins with significant sequence similarity only to the sequences from the group of teleosts and 2)  <dig> proteins that had no significant sequence similarity to any sequences within the species we used in this study .
in order to further eliminate false positives due to the incompleteness of the annotated protein sets and genomes, these two sets of proteins were then searched against the uniprot knowledgebase  and the current non redundant protein database in ncbi using blastp.
after manual inspection of the alignments, several genes initially assigned to the taxonomic categories were moved to the correct classification.
finally, the protein id of evolutionarily conserved proteins were transformed to their coding genes’ id and the final sets of ctsgs, orphan genes, and evolutionarily conserved genes  contained  <dig>   <dig>  and  <dig>  d. rerio genes, respectively .
characterization of ctsgs and orphan genes
several general tendencies of lineage-specific genes and orphan genes  in different taxa have been observed in previous studies.
first, younger genes tend to have shorter gene size, shorter proteins, and fewer exon numbers than older ones  <cit> .
second, younger genes show lower expression on average  <cit> .
to determine whether this was the case for the two sets of lineage-specific genes within zebrafish, we characterized and compared the genic properties of the ctsgs and orphan genes with those of the evolutionarily conserved gene sets .
just like the orphan genes in other species, such as primates, insects, and plants, the average gene size, average protein size, and average exon numbers per gene of the orphan genes were significantly lower than that of the ctsgs and ecs gene sets .
however, the gc content of the orphan genes was significantly higher than that of the ctsgs and the ecs gene sets , with the ctsgs having the lowest gc content.
the observation that the ctsgs had the lowest gc content was consistent with several previous studies, specifically those about lower gc content of taxonomically-restricted genes in drosophila  <cit> , honey bee  <cit> , and arabidopsis thaliana <cit> .
meanwhile, the observation that the orphan genes within zebrafish had the highest gc content was consistent with the elevated gc content within the coding sequences of the lineage-specific genes in rice  <cit> .
overall, the genic features of the two sets of lineage-specific genes within zebrafish were distinct from the evolutionarily conserved gene sets and were comparable to their counterparts in other species.
the physical locations of the two sets of lineage-specific genes within zebrafish were assigned to the  <dig> zebrafish chromosomes according to the information from ensembl version  <dig>  <cit> .
the number as well as the percentage of lineage-specific genes on each chromosome was also counted.
neither the ctsgs nor the orphan genes were distributed evenly within the different chromosomes of zebrafish .
for the orphan genes, four chromosomes  harbored a total of  <dig> orphan genes, each with seven to twelve orphan genes, which is more than fifty percent of the total number of orphan genes in zebrafish.
however, there were seven chromosomes  which harbored no orphan genes.
the chromosomal distribution of the ctsgs was basically similar to that of the orphan genes.
approximately fifty percent of ctsgs  were located within five chromosomes , whereas five chromosomes  each housed only one or no ctsgs.
in order to test whether the percentage of lineage-specific genes on each of chromosome correlated with the length of the chromosomes in zebrafish, spearman’s test was employed.
spearman’s test showed that the number of ctsgs on each chromosome significantly correlated with the length of their respective chromosome , whereas the number of orphan genes on each chromosome did not correlate with the length of their respective chromosome .
our results showed a trend similar to previous studies that found that younger genes exhibited lower expression on average  <cit> .
to determine whether a gene had evidence of expression from est or full-length cdna , we first downloaded the unigene data using biomart.
if a gene model of zebrafish was annotated with an est or fl-cdna in ensembl, then we considered that gene to have transcript support.
combining these two sets of results, the percentage of genes with transcript support was determined .
the transcript support for orphan genes  and for the ctsgs  was significantly lower than that for ec genes , with the orphan genes having the lowest transcript support .
functional inference using protfun
it has been shown that few of the highly taxonomically restricted genes  have been the focus of experimental work or could be characterized by go categories  <cit> .
in order to determine whether this is the case for the two sets of lineage-specific genes within zebrafish, the function annotations of ctsgs and orphan genes available at ensembl  <cit>  were explored using biomart  <cit> .
as expected, a significantly large percentage  of orphan genes were annotated as uncharacterized proteins  or had no description information  and there was only one gene annotated with go term accession.
likewise, about  <dig> % of ctsgs were annotated as uncharacterized proteins and hypothetical proteins or even have no description information and approximately  <dig> % of ctsgs were without go term annotations.
these observations suggest that the functions of most of the two sets of lineage-specific genes are unknown and some of the annotated genes may be the result of incorrect annotations, considering that there were seven orphan genes encoding proteins less than ten amino acids.
with respect to the remaining orphan genes with functional annotations in ensembl, there was no bias seen in specific functions; however, there were about eight genes of ctsgs involved in immunity, indicating that immune response may be very important in the radiation of teleosts.
to further explore the functions of the two sets of lineage-specific genes within zebrafish, their cellular roles and go categories were predicted using the protfun  <dig>  server .
as with the orphan genes in mice, orphan genes and ctsgs within zebrafish that were functionally associated with translation and energy metabolism were the most commonly presented, followed by those associated with transport, binding, and regulatory functions  <cit> .
furthermore, amino acid biosynthesis, replication, transcription, and the cell envelope were well represented compared with other cellular categories such as the biosynthesis of cofactors, cellular processes, purines, and pyrimidines.
the percentages of orphan genes and ctsgs with the same cellular roles were generally similar except one; specifically, about 33% of ctsgs were associated with the cell envelope, which was significantly higher than that of the orphan genes .
with respect to the go categories predicted by protfun, about 70% and 80% of orphan genes and ctsgs, respectively, were able to be assigned to a go category.
the growth factor category  was a relatively abundant category for both orphan genes and ctsgs.
strikingly and interestingly, the proteins of ctsgs involved in immune response were the most represented, exhibiting a percentage of  <dig> %, which is approximately two times that of the same go category for the orphan genes, implying a significantly larger expansion of these genes in teleosts.
in contrast, the proteins of orphan genes involved in transcription regulation were the second most represented , exhibiting a percentage about three times more than that of the same category of proteins in ctsgs.
although there was a remarkable difference in the percentage, the numbers of both orphan genes and ctsgs of this go category were very similar .
in addition, the proteins associated with signal transducer, receptor, hormone, structural protein, and stress response were well represented compared with the go category of ion channel proteins, which included associations to the ion channel, voltage-gated ion channel, cation channel, and metal ion transport.
these functional data suggest that lineage-specific genes may indeed be excellent candidates for genes involved in important housekeeping functions in a specific clade or species  <cit> .
high percentage of lineage-specific genes generated via gene duplications
gene duplication followed by rapid sequence divergence by a paralog, which is beyond the threshold of similarity searches, has long been thought to be the major mechanism that provided raw materials for the emergence of new genes since the publication of the famous monograph authored by susumu ohno  <cit> , although there are several other hypotheses regarding the origin of the lineage-specific genes, such as horizontal gene transfer  <cit> , an accelerated evolutionary rate  <cit> , de novo emergence from non-genic sequences  <cit>  as well as artifacts from genome annotation  <cit> .
in order to determine the proportion of contribution by gene duplication to the lineage-specific genes within zebrafish, we sought to identify such lineage-specific genes generated from the duplication-divergence mechanism using a simple method; that is, we determined whether any paralogs of the lineage-specific genes were widely evolutionarily conserved  <cit> .
to ascertain a minimum percentage of genes that may be generated via gene doubling, we first downloaded the information regarding the paralogs of ctsgs and orphan genes annotated in ensembl version  <dig> using biomart.
secondly, paralogs that were also found to be evolutionarily conserved were considered to be associated with lineage-specific genes created via gene duplication followed by rapid sequence divergence.
through this analysis, we found about  <dig> %  of orphan genes had paralogs, of which only five  were evolutionarily conserved.
that is, only about  <dig> % of orphan genes were generated via gene doubling.
compared with the relatively low percentage of orphan genes formed by gene duplication, the contribution of gene duplication to the generation of ctsgs was very remarkable.
as much as  <dig> % of ctsgs could have emerged as a result of the duplication-divergence mechanism, which was much higher than previous reports for insects  <cit>  and primate  <cit> , but lower than reports for drosophila  <cit> .
in addition, if a lineage-specific gene is generated by gene duplication followed by rapid sequence divergence, then the similarity between the lineage-specific gene and its evolutionarily conserved paralogs should be lower than that between the lineage-specific gene and its evolutionarily not conserved paralogs.
to test this hypothesis, we further explored the similarity between the lineage-specific genes and their paralogs, both evolutionarily conserved and not evolutionarily conserved.
consistent with the hypothesis, the average similarity of ctsgs with their evolutionarily conserved paralogs was significantly lower than their average similarity with the other group of paralogs .
however, there was no significant difference between orphan genes and their paralogs .
these data further confirm that gene duplication followed by rapid sequence divergence was a major mechanism generating new genes  <cit> .
lists of all orphan genes and ctsgs with putative parent genes in the zebrafish genome that are widely conserved as well as not conserved are provided .
orphan genes are preferentially expressed in the reproductive system
the est database is a collection of millions of ests gathered from thousands of rna libraries covering dozens of zebrafish organs or tissues at different developmental stages  <cit> .
in order to elucidate the expression patterns of the lineage-specific genes in zebrafish, we first analyzed this comprehensive dataset to detect whether the lineage-specific genes were preferentially expressed in certain tissues or organs.
one key step in achieving this goal was the reliable mapping of est and full-length cdna to genomic sequences, so we followed a relatively stringent pipeline  <cit>  to retain high quality mappings .we counted a gene as expressed in a tissue or organ as long as it was supported by one est.
through this stringent pipeline, we found that orphan genes were preferentially expressed in reproductive organs and tissues compared with other organs at a 5% significance level .
the same test for the brain showed no significant enrichment of orphan genes, which was the same observation when analyzing primate-specific genes  <cit> .
nevertheless, there was no significant enrichment in gonad expression for the ctsgs.
considering that the est database covers a large number of tissues and organs, these observations suggest that the orphan genes within zebrafish were preferentially enriched for expression of the reproductive system compared with the ctsgs.
our result is in accordance with the association of recent-origin genes with both reproductive expression  <cit>  and reproduction-related behavior  <cit> .
temporal and spatial expression profiles of lineage-specific genes
knowing the expression patterns of lineage-specific genes at different developmental stages and in different tissues or organs is essential to illuminate whether the lineage-specific genes have corresponding biological function.
therefore, we exploited the rna-seq data to assess a far more precise measurement of the temporal and spatial expression profiles of these lineage-specific genes.
rna-seq is a recently developed method to transcriptome profiling that uses high throughput technologies and is a powerful method to quantify the expression of genes  <cit> .
a time-series of rna-seq data from  <dig> time-points during early zebrafish organogenesis that mark important developmental stages was obtained from previous studies  <cit> .
rna-seq data from  <dig> different zebrafish tissues were downloaded from ncbi  <cit> .
the rpkm, defined here as the number of unique mapped reads to the coding regions divided by one thousandth of the total length of all the exons of the gene, subsequently normalized by dividing by one millionth of the total number of valid reads, was calculated .
evidence for the expression of the lineage-specific genes in zebrafish, here defined as the mapped reads to their coding regions, was found in the transcriptome data for  <dig> of the  <dig> orphan genes  and  <dig> of the  <dig> ctsgs .
these percentages were significantly higher than those observed using est/fl-cdna data.
for the remaining  <dig> orphan genes not represented by rna-seq data, one had evidence of expression in est data and one had evidence of expression in protein data.
thus, the failure to find evidence for expression of these two genes with rna-seq data may suggest that the two genes were expressed in other developmental stages and tissues or at a very low level in the analyzed developmental stages and tissues, or the expression evidence based on the est and protein data may be incorrect .
as for the other  <dig> orphan genes not represented by rna-seq data, their gene models have probably been incorrectly annotated, considering that each of them had a transcript less than 60 bp.
on the other hand, expression evidence for the  <dig> ctsgs not represented by rna-seq data had all been found from est data.
therefore, the failure to find evidence for expression using rna-seq data may suggest that these genes were expressed in other developmental stages or tissues.
for example, gene ensdarg <dig> was expressed in olfactory tissue at about 3–4 months old and gene ensdarg <dig> was expressed in kidney tissue, which was not included in our rna-seq data.
in addition to the high tissue preferences observed in previous studies, we also discovered two intriguing temporal expression patterns of the lineage-specific genes within zebrafish.
first, lineage-specific genes were more likely to be enriched in later-stage embryos and early larval stages with a higher fraction of orphan genes.
the mean normalized level of gene expression, as defined here by the sum of all the rpkm for all developmental stages divided by the total number of lineage-specific genes, was generally higher in later-stage embryos and early larval stages than in early-stage embryos .
on the other hand, a larger proportion of the lineage-specific genes were expressed during the later-stage embryos and early larval stages than during the early-stage embryos .
interestingly, upon sorting the normalized expression levels of the lineage-specific genes for the  <dig> time-points from highest to lowest, we found that the proportion of genes having the highest expression level in the developmental stage, which was defined as the numbers of genes having highest expression level in the developmental stage divided by total gene number, was also higher in later-stage embryos and early larval stages than in early-stage embryos .
a similar pattern was also observed for the proportion of genes exhibiting the second, third, and fourth highest levels of expression .
finally, the percentage of high expression level in later-stage embryos and early larval stages was higher for orphan genes compared to the ctsgs .
collectively, these results strongly suggest that lineage-specific genes were more preferentially to be enriched in later-stage embryos and early larval stages.
second, we calculated the shannon entropy-based specificity score per genes as a measure of expression level divergence during zebrafish embryogenesis  to test whether there were more pronounced changes in gene expression level between two consecutive developmental stages for lineage-specific genes than for evolutionarily conserved genes.
we found that the entire two sets of lineage-specific genes showed an increased temporal specificity compared with the evolutionarily conserved genes .
this phenomenon suggests that lineage-specific genes had a more restricted temporal expression than the ec genes.
in addition, rna-seq data showed that orphan genes were preferentially expressed in the ovary and ctsgs were enriched for brain expression.
the tissue that had the largest proportion of the genes expressed and the proportion of genes exhibiting the highest expression level in the tissue basically support the above conclusion, which is also consistent with the evidence from the est data.
several individual lineage-specific genes were also found to have intriguing expression patterns .
for example, orphan gene ensdarg <dig> exhibited significantly high expression in the ovary; thus, we speculated that that gene may play a role in reproduction.
orphan gene ensdarg <dig> was highly, and nearly specifically, expressed in the post-mbt stages, suggesting that this gene may contribute to the development of later-stage embryos and early larval stages.
one of ctsgs, ensdarg <dig>  was highly expressed in both female and male brains, which indicates that this gene may have a role in zebrafish brain development.
another ctsg, ensdarg <dig>  exhibited significantly high expression in the early-stage embryos, elucidating its important contribution to the development of early embryos development.
a list of the number of unique mapped reads and the rpkm of each lineage-specific gene is provided .
in order to validate the expression pattern of these lineage-specific genes, reverse transcription polymerase chain reaction  assay was used.
primers for  <dig> lineage-specific genes  were designed and all of these genes were amplified.
the information of the primers and the results of rt-pcr were provided .
discussion
enormous lineage-specific genes identified in other taxa with potentially important functions  <cit>  motivated our genomewide search for lineage-specific genes within zebrafish.
here, we adopted blast  <cit> , the preferred method for detecting homologs, and phylostratigraphy  <cit>  to identify two sets of lineage-specific genes within zebrafish.
then we characterized these genes, predicted their functions ab initio, inferred their evolutionary origin, and analyzed their expression patterns, making this the most comprehensive study of lineage-specific genes within teleosts and zebrafish to date.
the  <dig> ctsgs and  <dig> orphan genes obtained in this study are attractive targets for future experimental discovery, owing to their lineage specificity and to the fact that the majority encode proteins whose functions are yet to be determined .
compared with the lineage-specific genes identified in plants  <cit> , the number of lineage-specific genes within zebrafish is significantly lower, which may reflect the basic difference between animals and plants, considering the likely small number of lineage-specific genes identified in primate  <cit>  and insects  <cit> .
although yang et al.
<cit>  identified a relatively small number of lineage-specific genes in arabidopsis, oryza, and populus, whose number is close to that in animals, their criteria used to define sequence conservation was too relaxed, making the validity of their results questionionable.
for example, they restricted their analysis to only the genes with expression evidence support and employed a very relaxed criterion to define sequence conservation  that has not been used in other studies.
taken together, the dramatic difference in number of lineage-specific genes observed between the genomes of animals and plants should not be the result of the method we used and may suggest that there is a remarkable genetic difference in terms of lineage-specific genes between the genomes of animals and plants.
in addition, this difference may suggest that genome doubling followed by sequence divergence occurred in plants at a higher frequency  <cit> , which may explain to some extent why there are many more lineage-specific genes in plant genomes.
both the ctsgs and orphan genes had shorter gene length compared with the ec genes, probably owing to fewer numbers of exons per gene and higher percentage of intronless lineage-specific genes.
for example, nearly 28% of orphan genes contained only one exon, while the percentage of single exon ec genes was only 6%.
one reason for such a difference may be that intronless genes can arise via retroposition, which has been confirmed to create a large amount of new genes in the zebrafish genome  <cit> .
alternatively, this difference may be a result of the “introns late” hypothesis, which assumes intron accretion into the protein-coding genes is continuous throughout the evolutionary time of eukaryotes  <cit> .
thus, the younger the genes are, the fewer exons they have.
additionally, since orphan genes are species specific, these genes may have arisen in relatively recent years.
collectively, these reasons may partly explain why young orphan genes contain a single exon and why lineage-specific genes are shorter than older evolutionarily conserved genes.
generally speaking, lineage-specific genes are thought to play significant roles in the evolution of lineage specific phenotypes and adaptive innovation  <cit> .
although there are a large number of lineage-specific genes whose functions have not been characterized and only one orphan gene and 31% ctsgs have go term accession, we were still able to find five orphan genes and eight ctsgs whose functions are closely related to immunity.
the significant enrichment of immunity proteins in the lineage-specific genes within zebrafish indicates that defense against pathogens may be an important goal in terms of the successful diversification of fishes.
fishes are an extremely diverse group of aquatic vertebrate animals that also exhibit enormous diversity in the habitats they occupy.
fishes live in almost every conceivable type of aquatic habitat, from an elevation up to  <dig>  meters in tibet to  <dig>  meters below the surface of the ocean and some species even make short excursions onto land.
some fishes can also live in almost pure freshwater, while others reside in very salty lakes.
they can tolerate temperatures ranging from as high as  <dig> °c to −2°c under the antarctic ice sheet  <cit> .
thus, fishes should be confronted with much more diverse pathogen invasion.
therefore, lineage-specific genes involved with immunity should help fishes better adapt to various pathogens and successfully survive within their diverse habitats.
in addition, the prediction of gene function is based on homology to proteins with known function in other species.
some lineage-specific genes lack homologs in other lineages, so we predicted their function ab initio.
interestingly, the proteins of ctsgs involved in immune response were the most represented, with a percentage of  <dig> % of ctsgs, probably implying a significantly larger expansion of these genes in teleosts.
therefore, function assignment both based on homology and prediction ab initio showed a significant enrichment in proteins related to immune response, suggesting that the successful adaptation of teleosts may be explained by their conserved lineage-specific genes.
variation of gene number within different organisms suggests a general process of new gene origination  <cit> .
one basic question in biology is the molecular mechanisms involved in the creation of new genes.
there have already been several hypotheses regarding the origin of lineage-specific genes.
however, determining the exact mechanisms regarding the origin of lineage-specific genes depends on the comparative genome analysis of taxonomically closely related species.
it is extremely difficult to achieve the aforementioned goal for research on fish so far.
gene duplication followed by rapid sequence divergence in one of the paralogs is a well explored source of lineage-specific genes  <cit> .
a simple method for determining such genes is to determine whether any of the paralogs of lineage-specific genes are widely evolutionarily conserved.
through this analysis, we found that there were a significantly large number of lineage-specific genes generated by gene duplication followed by rapid sequence divergence of one of the paralogs.
it was also confirmed by observing that the similarity between the genes and their evolutionarily conserved paralogs was lower than the similarity between the genes and their paralogs not evolutionarily conserved.
as for other mechanisms forming lineage-specific genes, we will explore these questions in the future when the genome sequence of the silver carp , a relatively close species to zebrafish, is released.
previous studies have shown that young new genes generated by various mechanisms seem to have been preferentially endowed with testis-specific or testis-biased expression patterns  <cit> .
in accordance with this observation, there are a significantly large number of new genes within zebrafish expressed in the reproductive system reflecting the expectation that emergence of new, lineage-specific genes may accompany speciation or reproduction.
this suggests that this expression pattern is a general phenomenon not only in mammals and drosophila, but also in teleosts.
there are several hypotheses which can explain this propensity.
first, sex- and reproduction-related genes are generally recognized as a class of rapidly evolving genes and undergo adaptive evolution after speciation events involved in male reproduction  <cit> .
furthermore, the testis is the most rapidly evolving organ owing to the strong selective pressures to which it is subjected because of its important roles in sperm competition, sexual conflict, reproductive isolation, germline pathogens, and mutations causing segregation distortion in the male germline  <cit> .
second, the “hypertranscription” state  <cit>  caused by chromatin remodelling and rna polymerase ii complexes in the meiotic and postmeiotic spermatogenic cells would favor the initial, unprovoked transcription of newly arisen genes  <cit> .
as for the ctsgs, however, no significant reproductive expression was enriched, which further confirmed that only the young new genes were specifically expressed in the testis, since the ctsgs were relatively older than the orphan genes.
expression analyses of lineage-specific genes using est or microarrays have elucidated the fact that lineage-specific novel genes are preferentially expressed in specific tissues or organs, such as the testis or brain  <cit> .
although est data covers a large number of samples, which could be used to compare the expression between different samples, the coverage of individual genes is too low to quantify the expression level of genes.
microarrays also have some limitations, such as cross-hybridization and saturation of signals  <cit> .
therefore, we used the rna-seq data from various developmental stages and tissues to quantify the lineage-specific genes and highlight two novel properties of these genes.
first, in addition to being highly tissue-specific, lineage-specific gene expression were highly temporally restricted.
second, lineage-specific genes were preferentially expressed in later-stage embryos and early larval stages compared with early-embryos.
the higher expression level of lineage-specific genes after the mbt suggests that lineage-specific genes are important components for the zygotic transcription.
maternally deposited mrnas direct early development before the initiation of zygotic transcription during mid-blastula transition  <cit> .
however, zygotic transcription plays a more important role in the regulation of development after mbt, since a high percentage of maternally stored mrna has been degraded during the post-mbt stages.
in addition, it has been shown that all vertebrate embryos must converge towards a narrow point, called phylotypic stage at which all vertebrate show high morphogenetic resemblance, to acquire the basic scheme on which subsequent differences will emerge  <cit> .
the phenomenon that more lineage-specific genes are expressed after the phylotypic stages may probably be linked to the acquisition of species-specific morphological traits.
all vertebrate resemble each other at the phylotypic stage, so the crucial steps to form the morphological differences between species resulting from the expression product after the phylotypic stages.
therefore, lineage-specific genes within zebrafish should be crucial for the significantly morphological diversity of teleosts.
on the other hand, lineage-specific genes showed relatively higher expression levels during early larval stages, making them candidates for functions in specific tissues and organs during organogenesis.
expression analysis using rna-seq from different tissues and organs supported the observations from the est data and further showed that orphan genes are preferentially expressed in reproductive tissues, which also confirmed the potential roles of lineage-specific genes during organogenesis.
CONCLUSIONS
in the study, we have identified two sets of lineage-specific genes, ctsgs and orphan genes, which are specific to teleosts and zebrafish, respectively.
the conserved teleost-specific genes were found to be especially enriched in proteins with immunity functions, implying that defense against invasion by diverse pathogens was critical to the successful diversification of teleosts.
we also revealed that, in addition to being highly tissue-preferred expression, lineage-specific genes are also highly temporally restricted and are preferentially expressed in later-stage embryos and early larval stages compared with early-embryos.
this study provides valuable information for further analysis of the functions of these genes during zebrafish embryogenesis and will be helpful in improving the understanding of the successful diversification of teleosts.

 teleosts, which roughly constitute 96% of all living fishes and half of the extant vertebrate, are the most phenotypically diversified and species-rich group of all the vertebrate species. the vast morphological and species diversity of teleosts have received intense attention worldwide because of their importance in both scientific research and aquaculture. however, the genomic basis of the complex phenotype of teleosts during evolution remains obscure, despite the increasing amount of genome and transcriptome sequence data available. one important mechanism for the phenotypic diversity of species is the duplication of genes and entire genomes. evidence has recently been accumulated to allow a consensus to be reached that all teleosts experienced an additional whole genome duplication, which occurred after the basal ray-finned fishes separate from the actinopterygian stem lineage but before the teleosts began radiation. combining the absolute dates and phylogenetic timing of the 3r duplication, some groups have thought that the fsgd might be causally related to an increase in the number of species as well as their biological diversity. however, if the fsgd was responsible for the evolutionary success and astounding biological diversification of teleosts, it must have occurred prior to the radiation of teleosts. with the fossil record, paleontological evidence have suggested that the first appearance of most of the extant teleosts was only about million years ago, which is shorter than the fsgd that occurred at least million years ago. thus, it would be inappropriate to think that the fsgd was a major driving force behind the rapid radiation of teleosts. furthermore, considering the large amount of fossil data for basal ray-finned fishes, a consequential fsgd would not provide sufficient explanation for the morphological complexity of teleosts. besides the fish-specific genome duplication, alternative explanations for the increasing morphological complexity of teleosts include their experience with a higher rate of chromosomal rearrangements and a faster evolution of protein sequences and conserved noncoding elements compared to cartilaginous fishes and mammals. their implications for the evolution and diversity of teleosts have been intensively discussed. however, conserved teleost lineage-specific genes have been poorly characterized. lineage-specific genes, also referred to as taxonomically restricted genes are defined as genes found in one particular taxonomic group but share no sequence similarity with genes from other lineages. with the advent of large-scale genome sequencing projects for a wide range of species, lineage-specific genes have been extensively studied in mammals, insects, plants, and microbial species. lineage-specific genes are a significantly abundant component of all genomes sequenced to-date, which defies an early hypothesis that an increasing database size would eventually reduce the number of lineage-specific genes. orphan genes were first discussed when analyzing the yeast genome; approximately one-third of the identified genes fell into this category. likewise in drosophila melanogaster, the most accurate and complete genome analyzed, lineage-specific genes were found to make up nearly % of the total genes. apart from being abundant, lineage-specific genes have also been thought to be important for lineage specific traits and adaptations. in hydra, for example, interspecific differences in tentacle formation are closely related to the changes in the expression of taxonomically restricted genes. and in drosophila, the flightin gene is specifically important for increasing the frequency of the flight muscle to deliver the maximum power to the wing, which is a rather specific adaptation for the dipterans. although abundant in quantity and important in functionality, the evolutionary origin of lineage-specific genes is still enigmatic. several hypotheses about the origin of lineage-specific genes, including gene duplication followed by rapid sequence divergence, lateral gene transfer, accelerated evolutionary rate, artifacts from genome annotation, as well as de novo evolution from noncoding sequences have been proposed. despite the fact that the origin and evolution of lineage-specific genes is still poorly understood, the identification, characterization, function, and expression analysis of lineage-specific genes may provide a better understand for lineage-specific adaptation, such as the successful diversification of teleosts. in this study, we identified conserved teleost-specific genes and orphan genes in zebrafish using comparative genomics. we then characterized each set of these genes by diverse features, including gene size, protein size, exon number, gc content, transcript support, and chromosomal locations. as a large portion of the ctsgs and orphan genes have no known function, ab initio predictions using protfun were performed to infer possible biological functions. we then explored the evolutionary origin of lineage-specific genes and performed a comprehensive analysis of their tissues and developmental stages specific expression patterns using the wealth of available expression data, including est and rna-seq data, which in turn provided important complementary datasets that may be used to uncover their functions in the future. collectively, identification of lineage-specific genes as well as orphan genes and future studies of their function by means of target gene knockdown or knockout will no doubt help increase our understanding of the molecular basis of the successful diversification of teleosts.  identification of ctsgs and orphan genes the procedure used to identify orphan genes and ctsgs was modified from previous studies. we used the blast-basic local alignment search tool with default parameters for all assignments. firstly, we used blastp to search for homologs of all proteins annotated of zebrafish in each of the other species-levels outside the teleostei protein sets with an e-value cutoff of 10- those proteins for which we could not find any homolog in any alternative-spliced forms were used for the next set of searches, which was done using tblastn to search for homologs against the genome of the same species. in this step of the search, we only performed the tblastn search when there were no annotated proteins of the genome of the species, such as the genomes of elephant shark and salmo salar. a total of danio rerio proteins with significant sequence similarity to either a protein or genomic sequence from a species outside the group of teleosts were defined as the evolutionarily conserved proteins. secondly, the remaining d. rerio proteins with no significant similarity to any sequence outside the group of teleosts were further searched against the proteins and the genome of the group of teleosts, including oryzias latipes, gasterosteus aculeatus, takifugu rubripes, tetraodon nigroviridis, oreochromis niloticus, gadus morhua, salmo salar, using blastp and tblastn with e-values the same as the first step. this step produced two datasets: 1) proteins with significant sequence similarity only to the sequences from the group of teleosts and 2) proteins that had no significant sequence similarity to any sequences within the species we used in this study. in order to further eliminate false positives due to the incompleteness of the annotated protein sets and genomes, these two sets of proteins were then searched against the uniprot knowledgebase and the current non redundant protein database in ncbi using blastp. after manual inspection of the alignments, several genes initially assigned to the taxonomic categories were moved to the correct classification. finally, the protein id of evolutionarily conserved proteins were transformed to their coding genes’ id and the final sets of ctsgs, orphan genes, and evolutionarily conserved genes contained and d. rerio genes, respectively. characterization of ctsgs and orphan genes several general tendencies of lineage-specific genes and orphan genes in different taxa have been observed in previous studies. first, younger genes tend to have shorter gene size, shorter proteins, and fewer exon numbers than older ones. second, younger genes show lower expression on average. to determine whether this was the case for the two sets of lineage-specific genes within zebrafish, we characterized and compared the genic properties of the ctsgs and orphan genes with those of the evolutionarily conserved gene sets. just like the orphan genes in other species, such as primates, insects, and plants, the average gene size, average protein size, and average exon numbers per gene of the orphan genes were significantly lower than that of the ctsgs and ecs gene sets. however, the gc content of the orphan genes was significantly higher than that of the ctsgs and the ecs gene sets, with the ctsgs having the lowest gc content. the observation that the ctsgs had the lowest gc content was consistent with several previous studies, specifically those about lower gc content of taxonomically-restricted genes in drosophila, honey bee, and arabidopsis thaliana. meanwhile, the observation that the orphan genes within zebrafish had the highest gc content was consistent with the elevated gc content within the coding sequences of the lineage-specific genes in rice. overall, the genic features of the two sets of lineage-specific genes within zebrafish were distinct from the evolutionarily conserved gene sets and were comparable to their counterparts in other species. the physical locations of the two sets of lineage-specific genes within zebrafish were assigned to the zebrafish chromosomes according to the information from ensembl version. the number as well as the percentage of lineage-specific genes on each chromosome was also counted. neither the ctsgs nor the orphan genes were distributed evenly within the different chromosomes of zebrafish. for the orphan genes, four chromosomes harbored a total of orphan genes, each with seven to twelve orphan genes, which is more than fifty percent of the total number of orphan genes in zebrafish. however, there were seven chromosomes which harbored no orphan genes. the chromosomal distribution of the ctsgs was basically similar to that of the orphan genes. approximately fifty percent of ctsgs were located within five chromosomes, whereas five chromosomes each housed only one or no ctsgs. in order to test whether the percentage of lineage-specific genes on each of chromosome correlated with the length of the chromosomes in zebrafish, spearman’s test was employed. spearman’s test showed that the number of ctsgs on each chromosome significantly correlated with the length of their respective chromosome, whereas the number of orphan genes on each chromosome did not correlate with the length of their respective chromosome. our results showed a trend similar to previous studies that found that younger genes exhibited lower expression on average. to determine whether a gene had evidence of expression from est or full-length cdna, we first downloaded the unigene data using biomart. if a gene model of zebrafish was annotated with an est or fl-cdna in ensembl, then we considered that gene to have transcript support. combining these two sets of results, the percentage of genes with transcript support was determined. the transcript support for orphan genes and for the ctsgs was significantly lower than that for ec genes, with the orphan genes having the lowest transcript support. functional inference using protfun it has been shown that few of the highly taxonomically restricted genes have been the focus of experimental work or could be characterized by go categories. in order to determine whether this is the case for the two sets of lineage-specific genes within zebrafish, the function annotations of ctsgs and orphan genes available at ensembl were explored using biomart. as expected, a significantly large percentage of orphan genes were annotated as uncharacterized proteins or had no description information and there was only one gene annotated with go term accession. likewise, about % of ctsgs were annotated as uncharacterized proteins and hypothetical proteins or even have no description information and approximately % of ctsgs were without go term annotations. these observations suggest that the functions of most of the two sets of lineage-specific genes are unknown and some of the annotated genes may be the result of incorrect annotations, considering that there were seven orphan genes encoding proteins less than ten amino acids. with respect to the remaining orphan genes with functional annotations in ensembl, there was no bias seen in specific functions; however, there were about eight genes of ctsgs involved in immunity, indicating that immune response may be very important in the radiation of teleosts. to further explore the functions of the two sets of lineage-specific genes within zebrafish, their cellular roles and go categories were predicted using the protfun server. as with the orphan genes in mice, orphan genes and ctsgs within zebrafish that were functionally associated with translation and energy metabolism were the most commonly presented, followed by those associated with transport, binding, and regulatory functions. furthermore, amino acid biosynthesis, replication, transcription, and the cell envelope were well represented compared with other cellular categories such as the biosynthesis of cofactors, cellular processes, purines, and pyrimidines. the percentages of orphan genes and ctsgs with the same cellular roles were generally similar except one; specifically, about 33% of ctsgs were associated with the cell envelope, which was significantly higher than that of the orphan genes. with respect to the go categories predicted by protfun, about 70% and 80% of orphan genes and ctsgs, respectively, were able to be assigned to a go category. the growth factor category was a relatively abundant category for both orphan genes and ctsgs. strikingly and interestingly, the proteins of ctsgs involved in immune response were the most represented, exhibiting a percentage of %, which is approximately two times that of the same go category for the orphan genes, implying a significantly larger expansion of these genes in teleosts. in contrast, the proteins of orphan genes involved in transcription regulation were the second most represented, exhibiting a percentage about three times more than that of the same category of proteins in ctsgs. although there was a remarkable difference in the percentage, the numbers of both orphan genes and ctsgs of this go category were very similar. in addition, the proteins associated with signal transducer, receptor, hormone, structural protein, and stress response were well represented compared with the go category of ion channel proteins, which included associations to the ion channel, voltage-gated ion channel, cation channel, and metal ion transport. these functional data suggest that lineage-specific genes may indeed be excellent candidates for genes involved in important housekeeping functions in a specific clade or species. high percentage of lineage-specific genes generated via gene duplications gene duplication followed by rapid sequence divergence by a paralog, which is beyond the threshold of similarity searches, has long been thought to be the major mechanism that provided raw materials for the emergence of new genes since the publication of the famous monograph authored by susumu ohno, although there are several other hypotheses regarding the origin of the lineage-specific genes, such as horizontal gene transfer, an accelerated evolutionary rate, de novo emergence from non-genic sequences as well as artifacts from genome annotation. in order to determine the proportion of contribution by gene duplication to the lineage-specific genes within zebrafish, we sought to identify such lineage-specific genes generated from the duplication-divergence mechanism using a simple method; that is, we determined whether any paralogs of the lineage-specific genes were widely evolutionarily conserved. to ascertain a minimum percentage of genes that may be generated via gene doubling, we first downloaded the information regarding the paralogs of ctsgs and orphan genes annotated in ensembl version using biomart. secondly, paralogs that were also found to be evolutionarily conserved were considered to be associated with lineage-specific genes created via gene duplication followed by rapid sequence divergence. through this analysis, we found about % of orphan genes had paralogs, of which only five were evolutionarily conserved. that is, only about % of orphan genes were generated via gene doubling. compared with the relatively low percentage of orphan genes formed by gene duplication, the contribution of gene duplication to the generation of ctsgs was very remarkable. as much as % of ctsgs could have emerged as a result of the duplication-divergence mechanism, which was much higher than previous reports for insects and primate, but lower than reports for drosophila. in addition, if a lineage-specific gene is generated by gene duplication followed by rapid sequence divergence, then the similarity between the lineage-specific gene and its evolutionarily conserved paralogs should be lower than that between the lineage-specific gene and its evolutionarily not conserved paralogs. to test this hypothesis, we further explored the similarity between the lineage-specific genes and their paralogs, both evolutionarily conserved and not evolutionarily conserved. consistent with the hypothesis, the average similarity of ctsgs with their evolutionarily conserved paralogs was significantly lower than their average similarity with the other group of paralogs. however, there was no significant difference between orphan genes and their paralogs. these data further confirm that gene duplication followed by rapid sequence divergence was a major mechanism generating new genes. lists of all orphan genes and ctsgs with putative parent genes in the zebrafish genome that are widely conserved as well as not conserved are provided. orphan genes are preferentially expressed in the reproductive system the est database is a collection of millions of ests gathered from thousands of rna libraries covering dozens of zebrafish organs or tissues at different developmental stages. in order to elucidate the expression patterns of the lineage-specific genes in zebrafish, we first analyzed this comprehensive dataset to detect whether the lineage-specific genes were preferentially expressed in certain tissues or organs. one key step in achieving this goal was the reliable mapping of est and full-length cdna to genomic sequences, so we followed a relatively stringent pipeline to retain high quality mappings.we counted a gene as expressed in a tissue or organ as long as it was supported by one est. through this stringent pipeline, we found that orphan genes were preferentially expressed in reproductive organs and tissues compared with other organs at a 5% significance level. the same test for the brain showed no significant enrichment of orphan genes, which was the same observation when analyzing primate-specific genes. nevertheless, there was no significant enrichment in gonad expression for the ctsgs. considering that the est database covers a large number of tissues and organs, these observations suggest that the orphan genes within zebrafish were preferentially enriched for expression of the reproductive system compared with the ctsgs. our result is in accordance with the association of recent-origin genes with both reproductive expression and reproduction-related behavior. temporal and spatial expression profiles of lineage-specific genes knowing the expression patterns of lineage-specific genes at different developmental stages and in different tissues or organs is essential to illuminate whether the lineage-specific genes have corresponding biological function. therefore, we exploited the rna-seq data to assess a far more precise measurement of the temporal and spatial expression profiles of these lineage-specific genes. rna-seq is a recently developed method to transcriptome profiling that uses high throughput technologies and is a powerful method to quantify the expression of genes. a time-series of rna-seq data from time-points during early zebrafish organogenesis that mark important developmental stages was obtained from previous studies. rna-seq data from different zebrafish tissues were downloaded from ncbi. the rpkm, defined here as the number of unique mapped reads to the coding regions divided by one thousandth of the total length of all the exons of the gene, subsequently normalized by dividing by one millionth of the total number of valid reads, was calculated. evidence for the expression of the lineage-specific genes in zebrafish, here defined as the mapped reads to their coding regions, was found in the transcriptome data for of the orphan genes and of the ctsgs. these percentages were significantly higher than those observed using est/fl-cdna data. for the remaining orphan genes not represented by rna-seq data, one had evidence of expression in est data and one had evidence of expression in protein data. thus, the failure to find evidence for expression of these two genes with rna-seq data may suggest that the two genes were expressed in other developmental stages and tissues or at a very low level in the analyzed developmental stages and tissues, or the expression evidence based on the est and protein data may be incorrect. as for the other orphan genes not represented by rna-seq data, their gene models have probably been incorrectly annotated, considering that each of them had a transcript less than 60 bp. on the other hand, expression evidence for the ctsgs not represented by rna-seq data had all been found from est data. therefore, the failure to find evidence for expression using rna-seq data may suggest that these genes were expressed in other developmental stages or tissues. for example, gene ensdarg was expressed in olfactory tissue at about 3–4 months old and gene ensdarg was expressed in kidney tissue, which was not included in our rna-seq data. in addition to the high tissue preferences observed in previous studies, we also discovered two intriguing temporal expression patterns of the lineage-specific genes within zebrafish. first, lineage-specific genes were more likely to be enriched in later-stage embryos and early larval stages with a higher fraction of orphan genes. the mean normalized level of gene expression, as defined here by the sum of all the rpkm for all developmental stages divided by the total number of lineage-specific genes, was generally higher in later-stage embryos and early larval stages than in early-stage embryos. on the other hand, a larger proportion of the lineage-specific genes were expressed during the later-stage embryos and early larval stages than during the early-stage embryos. interestingly, upon sorting the normalized expression levels of the lineage-specific genes for the time-points from highest to lowest, we found that the proportion of genes having the highest expression level in the developmental stage, which was defined as the numbers of genes having highest expression level in the developmental stage divided by total gene number, was also higher in later-stage embryos and early larval stages than in early-stage embryos. a similar pattern was also observed for the proportion of genes exhibiting the second, third, and fourth highest levels of expression. finally, the percentage of high expression level in later-stage embryos and early larval stages was higher for orphan genes compared to the ctsgs. collectively, these results strongly suggest that lineage-specific genes were more preferentially to be enriched in later-stage embryos and early larval stages. second, we calculated the shannon entropy-based specificity score per genes as a measure of expression level divergence during zebrafish embryogenesis to test whether there were more pronounced changes in gene expression level between two consecutive developmental stages for lineage-specific genes than for evolutionarily conserved genes. we found that the entire two sets of lineage-specific genes showed an increased temporal specificity compared with the evolutionarily conserved genes. this phenomenon suggests that lineage-specific genes had a more restricted temporal expression than the ec genes. in addition, rna-seq data showed that orphan genes were preferentially expressed in the ovary and ctsgs were enriched for brain expression. the tissue that had the largest proportion of the genes expressed and the proportion of genes exhibiting the highest expression level in the tissue basically support the above conclusion, which is also consistent with the evidence from the est data. several individual lineage-specific genes were also found to have intriguing expression patterns. for example, orphan gene ensdarg exhibited significantly high expression in the ovary; thus, we speculated that that gene may play a role in reproduction. orphan gene ensdarg was highly, and nearly specifically, expressed in the post-mbt stages, suggesting that this gene may contribute to the development of later-stage embryos and early larval stages. one of ctsgs, ensdarg was highly expressed in both female and male brains, which indicates that this gene may have a role in zebrafish brain development. another ctsg, ensdarg exhibited significantly high expression in the early-stage embryos, elucidating its important contribution to the development of early embryos development. a list of the number of unique mapped reads and the rpkm of each lineage-specific gene is provided. in order to validate the expression pattern of these lineage-specific genes, reverse transcription polymerase chain reaction assay was used. primers for lineage-specific genes were designed and all of these genes were amplified. the information of the primers and the results of rt-pcr were provided. discussion enormous lineage-specific genes identified in other taxa with potentially important functions motivated our genomewide search for lineage-specific genes within zebrafish. here, we adopted blast, the preferred method for detecting homologs, and phylostratigraphy to identify two sets of lineage-specific genes within zebrafish. then we characterized these genes, predicted their functions ab initio, inferred their evolutionary origin, and analyzed their expression patterns, making this the most comprehensive study of lineage-specific genes within teleosts and zebrafish to date. the ctsgs and orphan genes obtained in this study are attractive targets for future experimental discovery, owing to their lineage specificity and to the fact that the majority encode proteins whose functions are yet to be determined. compared with the lineage-specific genes identified in plants, the number of lineage-specific genes within zebrafish is significantly lower, which may reflect the basic difference between animals and plants, considering the likely small number of lineage-specific genes identified in primate and insects. although yang et al. identified a relatively small number of lineage-specific genes in arabidopsis, oryza, and populus, whose number is close to that in animals, their criteria used to define sequence conservation was too relaxed, making the validity of their results questionionable. for example, they restricted their analysis to only the genes with expression evidence support and employed a very relaxed criterion to define sequence conservation that has not been used in other studies. taken together, the dramatic difference in number of lineage-specific genes observed between the genomes of animals and plants should not be the result of the method we used and may suggest that there is a remarkable genetic difference in terms of lineage-specific genes between the genomes of animals and plants. in addition, this difference may suggest that genome doubling followed by sequence divergence occurred in plants at a higher frequency, which may explain to some extent why there are many more lineage-specific genes in plant genomes. both the ctsgs and orphan genes had shorter gene length compared with the ec genes, probably owing to fewer numbers of exons per gene and higher percentage of intronless lineage-specific genes. for example, nearly 28% of orphan genes contained only one exon, while the percentage of single exon ec genes was only 6%. one reason for such a difference may be that intronless genes can arise via retroposition, which has been confirmed to create a large amount of new genes in the zebrafish genome. alternatively, this difference may be a result of the “introns late” hypothesis, which assumes intron accretion into the protein-coding genes is continuous throughout the evolutionary time of eukaryotes. thus, the younger the genes are, the fewer exons they have. additionally, since orphan genes are species specific, these genes may have arisen in relatively recent years. collectively, these reasons may partly explain why young orphan genes contain a single exon and why lineage-specific genes are shorter than older evolutionarily conserved genes. generally speaking, lineage-specific genes are thought to play significant roles in the evolution of lineage specific phenotypes and adaptive innovation. although there are a large number of lineage-specific genes whose functions have not been characterized and only one orphan gene and 31% ctsgs have go term accession, we were still able to find five orphan genes and eight ctsgs whose functions are closely related to immunity. the significant enrichment of immunity proteins in the lineage-specific genes within zebrafish indicates that defense against pathogens may be an important goal in terms of the successful diversification of fishes. fishes are an extremely diverse group of aquatic vertebrate animals that also exhibit enormous diversity in the habitats they occupy. fishes live in almost every conceivable type of aquatic habitat, from an elevation up to meters in tibet to meters below the surface of the ocean and some species even make short excursions onto land. some fishes can also live in almost pure freshwater, while others reside in very salty lakes. they can tolerate temperatures ranging from as high as °c to −2°c under the antarctic ice sheet. thus, fishes should be confronted with much more diverse pathogen invasion. therefore, lineage-specific genes involved with immunity should help fishes better adapt to various pathogens and successfully survive within their diverse habitats. in addition, the prediction of gene function is based on homology to proteins with known function in other species. some lineage-specific genes lack homologs in other lineages, so we predicted their function ab initio. interestingly, the proteins of ctsgs involved in immune response were the most represented, with a percentage of % of ctsgs, probably implying a significantly larger expansion of these genes in teleosts. therefore, function assignment both based on homology and prediction ab initio showed a significant enrichment in proteins related to immune response, suggesting that the successful adaptation of teleosts may be explained by their conserved lineage-specific genes. variation of gene number within different organisms suggests a general process of new gene origination. one basic question in biology is the molecular mechanisms involved in the creation of new genes. there have already been several hypotheses regarding the origin of lineage-specific genes. however, determining the exact mechanisms regarding the origin of lineage-specific genes depends on the comparative genome analysis of taxonomically closely related species. it is extremely difficult to achieve the aforementioned goal for research on fish so far. gene duplication followed by rapid sequence divergence in one of the paralogs is a well explored source of lineage-specific genes. a simple method for determining such genes is to determine whether any of the paralogs of lineage-specific genes are widely evolutionarily conserved. through this analysis, we found that there were a significantly large number of lineage-specific genes generated by gene duplication followed by rapid sequence divergence of one of the paralogs. it was also confirmed by observing that the similarity between the genes and their evolutionarily conserved paralogs was lower than the similarity between the genes and their paralogs not evolutionarily conserved. as for other mechanisms forming lineage-specific genes, we will explore these questions in the future when the genome sequence of the silver carp, a relatively close species to zebrafish, is released. previous studies have shown that young new genes generated by various mechanisms seem to have been preferentially endowed with testis-specific or testis-biased expression patterns. in accordance with this observation, there are a significantly large number of new genes within zebrafish expressed in the reproductive system reflecting the expectation that emergence of new, lineage-specific genes may accompany speciation or reproduction. this suggests that this expression pattern is a general phenomenon not only in mammals and drosophila, but also in teleosts. there are several hypotheses which can explain this propensity. first, sex- and reproduction-related genes are generally recognized as a class of rapidly evolving genes and undergo adaptive evolution after speciation events involved in male reproduction. furthermore, the testis is the most rapidly evolving organ owing to the strong selective pressures to which it is subjected because of its important roles in sperm competition, sexual conflict, reproductive isolation, germline pathogens, and mutations causing segregation distortion in the male germline. second, the “hypertranscription” state caused by chromatin remodelling and rna polymerase ii complexes in the meiotic and postmeiotic spermatogenic cells would favor the initial, unprovoked transcription of newly arisen genes. as for the ctsgs, however, no significant reproductive expression was enriched, which further confirmed that only the young new genes were specifically expressed in the testis, since the ctsgs were relatively older than the orphan genes. expression analyses of lineage-specific genes using est or microarrays have elucidated the fact that lineage-specific novel genes are preferentially expressed in specific tissues or organs, such as the testis or brain. although est data covers a large number of samples, which could be used to compare the expression between different samples, the coverage of individual genes is too low to quantify the expression level of genes. microarrays also have some limitations, such as cross-hybridization and saturation of signals. therefore, we used the rna-seq data from various developmental stages and tissues to quantify the lineage-specific genes and highlight two novel properties of these genes. first, in addition to being highly tissue-specific, lineage-specific gene expression were highly temporally restricted. second, lineage-specific genes were preferentially expressed in later-stage embryos and early larval stages compared with early-embryos. the higher expression level of lineage-specific genes after the mbt suggests that lineage-specific genes are important components for the zygotic transcription. maternally deposited mrnas direct early development before the initiation of zygotic transcription during mid-blastula transition. however, zygotic transcription plays a more important role in the regulation of development after mbt, since a high percentage of maternally stored mrna has been degraded during the post-mbt stages. in addition, it has been shown that all vertebrate embryos must converge towards a narrow point, called phylotypic stage at which all vertebrate show high morphogenetic resemblance, to acquire the basic scheme on which subsequent differences will emerge. the phenomenon that more lineage-specific genes are expressed after the phylotypic stages may probably be linked to the acquisition of species-specific morphological traits. all vertebrate resemble each other at the phylotypic stage, so the crucial steps to form the morphological differences between species resulting from the expression product after the phylotypic stages. therefore, lineage-specific genes within zebrafish should be crucial for the significantly morphological diversity of teleosts. on the other hand, lineage-specific genes showed relatively higher expression levels during early larval stages, making them candidates for functions in specific tissues and organs during organogenesis. expression analysis using rna-seq from different tissues and organs supported the observations from the est data and further showed that orphan genes are preferentially expressed in reproductive tissues, which also confirmed the potential roles of lineage-specific genes during organogenesis.  in the study, we have identified two sets of lineage-specific genes, ctsgs and orphan genes, which are specific to teleosts and zebrafish, respectively. the conserved teleost-specific genes were found to be especially enriched in proteins with immunity functions, implying that defense against invasion by diverse pathogens was critical to the successful diversification of teleosts. we also revealed that, in addition to being highly tissue-preferred expression, lineage-specific genes are also highly temporally restricted and are preferentially expressed in later-stage embryos and early larval stages compared with early-embryos. this study provides valuable information for further analysis of the functions of these genes during zebrafish embryogenesis and will be helpful in improving the understanding of the successful diversification of teleosts. 
in this study, we performed comparative genomic analysis to discover the conserved teleost-specific genes and orphan genes within zebrafish and found that these two sets of lineage-specific genes may have played important roles during zebrafish embryogenesis. lineage-specific genes within zebrafish share many of the characteristics of their counterparts in other species: shorter length, fewer exon numbers, higher gc content, and fewer of them have transcript support. the evolutionary origin of the lineage-specific genes was determined and a very high percentage of lineage-specific genes were generated via gene duplications. the temporal and spatial expression profile of lineage-specific genes obtained by expressed sequence tags and rna-seq data revealed two novel properties: in addition to being highly tissue-preferred expression, lineage-specific genes are also highly temporally restricted, namely they are expressed in narrower time windows than evolutionarily conserved genes and are specifically enriched in later-stage embryos and early larval stages. our study provides the first systematic identification of two different sets of lineage-specific genes within zebrafish and provides valuable information leading towards a better understanding of the molecular mechanisms of the genomic basis of teleost phenotypic complexity for future studies. fish-specific genome duplication cannot provide sufficient explanation for the morphological complexity of teleosts, considering the relatively large number of extinct basal ray-finned fishes. chromosomal location analysis indicated that neither the ctsgs nor the orphan genes were distributed evenly in the chromosomes of zebrafish. the genomic basis of teleost phenotypic complexity remains obscure, despite increasing availability of genome and transcriptome sequence data. the significant enrichment of immunity proteins in ctsgs annotated by gene ontology or predicted ab initio may imply that defense against pathogens may be an important reason for the diversification of teleosts.
 the genomic basis of teleost phenotypic complexity remains obscure, despite increasing availability of genome and transcriptome sequence data. fish-specific genome duplication cannot provide sufficient explanation for the morphological complexity of teleosts, considering the relatively large number of extinct basal ray-finned fishes.  in this study, we performed comparative genomic analysis to discover the conserved teleost-specific genes and orphan genes within zebrafish and found that these two sets of lineage-specific genes may have played important roles during zebrafish embryogenesis. lineage-specific genes within zebrafish share many of the characteristics of their counterparts in other species: shorter length, fewer exon numbers, higher gc content, and fewer of them have transcript support. chromosomal location analysis indicated that neither the ctsgs nor the orphan genes were distributed evenly in the chromosomes of zebrafish. the significant enrichment of immunity proteins in ctsgs annotated by gene ontology or predicted ab initio may imply that defense against pathogens may be an important reason for the diversification of teleosts. the evolutionary origin of the lineage-specific genes was determined and a very high percentage of lineage-specific genes were generated via gene duplications. the temporal and spatial expression profile of lineage-specific genes obtained by expressed sequence tags and rna-seq data revealed two novel properties: in addition to being highly tissue-preferred expression, lineage-specific genes are also highly temporally restricted, namely they are expressed in narrower time windows than evolutionarily conserved genes and are specifically enriched in later-stage embryos and early larval stages.  our study provides the first systematic identification of two different sets of lineage-specific genes within zebrafish and provides valuable information leading towards a better understanding of the molecular mechanisms of the genomic basis of teleost phenotypic complexity for future studies. teleostlineage-specific genetranscriptomezebrafish embryogenesis 