<title>Background</title> <p>Understanding functional innovation is one of the most interesting questions in biology. One important mechanism of functional innovation involves changes in gene expression [##REF##10449761##1##] caused by cis-regulatory mutations [##REF##17304246##2##]. While structural mutations within existing genes are an alternative mechanism to generate new functions [##REF##17492956##3##], another possibility is the emergence of new genes. Several possible mechanisms are known to be involved in creating novel genes [##REF##14634634##4##]. The best described origins of novel genes are gene duplication [##UREF##0##5##] and exon shuffling [##UREF##1##6##,##REF##10570989##7##]. Recently it has been shown that novel genes could also originate <italic>de novo </italic>from non-coding regions [##REF##16777968##8##]. Comparative genome analyses permit the identification of previously uncharacterized genes through sequence conservation, but the identification of rapidly evolving genes or genes of very recent origin is frequently restricted to <italic>in silico </italic>predictions. As novel genes are typically short [##REF##16777968##8##,##REF##17435230##9##], these may be easily missed. Alternatively, gene expression could serve as a good indicator for the presence of a gene. Hence, either Expressed Sequence Tag (EST) databases or reverse SAGE [##REF##16456128##10##,##REF##12213963##11##] could be used to identify novel transcripts.</p> <p>Drosophila served as model for the identification of novel genes since the 1990s. One of first novel genes in this genus was <italic>jingwei </italic>in <italic>D. melanogaster </italic>[##REF##7682012##12##], which is a fusion of two genes, a retroposed copy of the <italic>alcohol dehydrogenase </italic>(<italic>Adh</italic>) gene and a duplicated copy of the <italic>yellow emperor </italic>(<italic>ymp</italic>) gene [##REF##10958846##13##]. Since then several studies applied phylogenetic methods to the growing databases aiming for the identification of novel genes. The majority of the novel genes have a sex-biased gene expression and some reports suggested that sex-biased genes change their expression pattern more rapidly than unbiased genes [##REF##12805547##14##,##REF##16612543##15##]. Furthermore, male-biased genes were shown to have a higher rate of protein evolution than unbiased genes [##REF##12907700##16##, ####REF##15282334##17##, ##REF##18039869##18####18039869##18##]. In a recent report comparing the pattern of gene expression in <italic>D. melanogaster </italic>and <italic>D. pseudoobscura </italic>we failed to find evidence for an unconditionally faster rate of sequence evolution of male-biased genes. Rather, only genes with a male-biased gene expression in <italic>D. melanogaster </italic>were found to evolve faster. Genes with a male-biased gene expression in <italic>D. pseudoobscura </italic>only were evolving at a similar rate as unbiased genes [##REF##16816428##19##]. As a large proportion of the sex-biased tags could not be mapped to the corresponding genes in <italic>D. pseudoobscura</italic>, the analysis of these tags should shed further light onto the pattern of protein evolution of sex-biased genes in <italic>D. pseudoobscura</italic>.</p> <p>In this study, we identified eight novel genes with sex-biased gene expression in <italic>D. pseudoobscura </italic>using GLGI (Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification). Consistent, with previous results [##REF##16777968##8##,##REF##17435230##9##], we observed significantly more novel genes with a male bias than with a female bias in gene expression. Interestingly, we found no significant excess of X-linked novel genes, as has been reported in the previous studies [##REF##16777968##8##,##REF##17435230##9##].</p>
<title>Methods</title> <title>GLGI analysis</title> <p><italic>D. pseudoobscura </italic>SAGE tags were obtained from Metta et al. [##REF##16816428##19##]. The GLGI procedure for the 3' extension of the SAGE tags was performed according to Chen et al. [##REF##11807982##52##]. In brief, total RNA was extracted separately from 15 male and 15 female virgin flies using Trizol (Invitrogen, Carlsbad, CA) and treated with DNAseI (MBI Fermentas) to digest the genomic DNA residues. PolyA mRNA was isolated using a 5' biotinylated and anchored oligo dT primer (5' biotin-ATCTAGAGCGGCCGC(T)<sub>16</sub>V) and streptavidin beads (Dynal). Double strand cDNA was synthesized using the "Double strand cDNA synthesis kit" (Invitrogen, Carlsbad, CA) and digested with the NlaIII isoschizomer Hin1II (MBI Fermentas). The 3' end fragments of the digested cDNA was recovered with Dynal beads. PCR amplification was performed using SAGE tag as 5' primer and anchoring sequence to the polyA was used as 3' primer. The PCR product was then cloned into TA vector (Invitrogen, Carlsbad, CA) and sequenced.</p> <title>Modified SAGE-GLGI method for detecting novel genes</title> <p>To verify the consistency of the abundance of novel genes, we also performed an independent analysis, which is a modification of the SAGE-GLGI approach. We used mRNA isolated from virgin females. Double stranded cDNA was digested with NlaIII and an adapter (5'-GCCTCCCTCGCGCCATCAGCATG-3' and 5'-CTGATGGCGCGAGGGAGGC-3') was ligated to the 5' end of the digested cDNA fragments. The resulting template was amplified using an adapter specific primer and the primer anchoring to oligo dT. The products were cloned into TA vector and sequenced using M13 primers.</p> <title>Bioinformatics</title> <p>BLASTN search was performed using the <italic>D. pseudoobscura </italic>database to identify the location of the GLGI (and modified SAGE-GLGI) fragments in the genome. As the <italic>D. pseudoobscura </italic>genome is only coarsely annotated we re-annotated the <italic>D. pseudoobscura </italic>genome for the regions of interest. We noted that a GeneWise analysis using the <italic>D. melanogaster </italic>gene model often results in a more complete protein prediction than the one available in the annotation database [##UREF##3##22##]. Hence we first performed GeneWise using those <italic>D. melanogaster </italic>genes that mapped in the proximity of the GLGI fragments to test if the <italic>D. pseudoobscura </italic>gene prediction may have been incomplete and the GLGI fragment matches to the extended protein coding region. As neither our GeneWise annotations nor the available annotation of <italic>D. pseudoobscura </italic>contained UTR sequences for most genes, we also used <italic>ab initio </italic>gene predictions. Specifically, we submitted the genomic region spanning the 5' part of the closest gene in the correct orientation and the genomic region matching to the GLGI fragment to a Genescan prediction. The coding potential of the transcripts was assessed using QRNA program [##REF##11801179##26##].</p> <title>Structure of novel genes</title> <p>To identify the 5' boundaries of the novel GLGI sequences, we performed 5' RACE (rapid amplification of cDNA ends) using GeneRacer RACE ready cDNA kit (Invitrogen, Carlsbad, CA). The 5' RACE analysis was performed on an independent extraction of total RNA.</p> <title>Evolution of novel genes</title> <p>Eight individuals of <italic>D. pseudoobscura </italic>from Mesa-Verde (Colarado, USA) were sequenced for the five novel genes putatively coding for a protein (see Additional File ##SUPPL##0##1##, Table 2 for the primer sequences). Orthologous regions in <italic>D. persimilis </italic>were obtained from flybase [##UREF##3##22##]. One individual of <italic>D. miranda</italic>, a close relative of <italic>D. pseudoobscura</italic>, was obtained from the Tucson stock center (stock number: 14011-0101.08) and sequenced for these loci at both genomic as well as transcriptomic level. A typical cycling profile consisted of 3 min denaturation at 94°C followed by 35 cycles of 94°C for 40 sec, 50°C for 50 sec, and extension at 72°C for 1 min. PCR products were sequenced for both strands using DYEnamic ET Terminator Sequencing Kit (GE Health care Bio-Sciences AB, Sweden) according to manufacturer's instructions. The extension products were purified with Sephadex G-50 fine (GE Health care Bio-Sciences AB, Sweden) and separated on a MegaBACE 500 automated capillary sequencer. Forward and reverse strands were assembled using CodonCode Aligner version 2.0.1 [##UREF##5##53##]. The sequences were aligned using clustalW [##REF##7984417##54##] and standard neutrality tests and McDonald-Kreitman test were performed using DnaSP version 4.10 [##REF##14668244##55##]. The pair wise synonymous and non synonymous substitutions between <italic>D. pseudoobscura </italic>and <italic>D. miranda </italic>were obtained using yn00 program of PAML [##REF##17483113##56##]. Likelihood ratio test was performed using codeml program of PAML using three species tree to test if the dN/dS ratio is significantly lower than 1. Multi locus HKA test [##REF##3110004##30##] was performed using HKA program written by Jody Hey [##UREF##6##57##]. All the sequences were submitted to GenBank (accession numbers <ext-link ext-link-type="gen" xlink:href="EU379026">EU379026</ext-link>–<ext-link ext-link-type="gen" xlink:href="EU379076">EU379076</ext-link>).</p>
<title>Results</title> <title>GLGI analysis</title> <p>We used recently published SAGE data to identify sex-biased tags in <italic>D. pseudoobscura </italic>[##REF##16816428##19##]. Previous analysis showed a substantially higher efficiency of tag to gene mapping for male-biased tags than for female-biased tags [##REF##16816428##19##]. As the <italic>D. pseudoobscura </italic>genome annotation is heavily based on <italic>D. melanogaster </italic>gene models, this may be due to a higher proportion of novel genes among the genes with a female-biased gene expression. To test this, we selected 20 male-biased and 64 female-biased tags that were previously not mapped, relatively highly expressed and showed significant difference in expression between the sexes (<italic>p </italic>&lt; 0.001), for further analysis. Using the GLGI method, we successfully generated longer 3'cDNA fragments for 44 SAGE tags. This success rate is in agreement with a previous GLGI analysis [##REF##15247327##20##]. The GLGI fragments include 11 male-biased and 33 female-biased tags (Table ##TAB##0##1##). Thirty female-biased (91%) tags were mapped close to putative orthologs of <italic>D. melanogaster </italic>in the 3' end while only six male-biased (55%) tags were in this category (Additional file ##SUPPL##0##1##, Tables 1a and 1b). The remaining eight tags are falling into genomic regions that show no sequence conservation between <italic>D. melanogaster </italic>and <italic>D. pseudoobscura </italic>and no genes are annotated in <italic>D. melanogaster</italic>. Furthermore, we could not identify these regions by either nucleotide or protein BLAST search in any other species whose genome is completely sequenced except for <italic>D. persimilis</italic>. This suggests that these tags potentially represent novel genes that arose in the obscura group and are absent in other Drosophila species. Alternatively, these genes may have diverged from their orthologs in other Drosophila species to such an extent that they could not be identified by BLAST algorithm. We identified five male-biased and three female-biased transcripts that are putatively novel. Hence, 46% of the male-biased and 9% of the female-biased tags are putatively novel in the genome of <italic>D. pseudoobscura</italic>. The difference between male- and female-biased tags is significant (<italic>p </italic>= 0.02), suggesting that the majority of the novel genes show male biased expression [##REF##16777968##8##,##REF##17435230##9##].</p> <title>Modified SAGE-GLGI approach</title> <p>As the GLGI analysis of sex-biased genes indicated a high number of novel genes, we were interested if this applies only to highly sex-biased genes or is a more general phenomenon. We sequenced 126 clones using a cDNA library that mimics the generation of GLGI fragments from SAGE tags. Of the 126 clones sequenced, 89 were unique and the remaining clones were redundant. From these 89 clones, 46 clones could be matched to the SAGE expression data [##REF##16816428##19##]. Out of 89, 35 clones could be mapped to the annotated genes of <italic>D. pseudoobscura</italic>. Thirty-two clones mapped close to predicted genes in the 3' end. Six clones fell unambiguously into genomic regions for which genes are predicted in <italic>D. melanogaster </italic>but not in <italic>D. pseudoobscura</italic>. For these clones we noted that the incomplete annotation in <italic>D. pseudoobscura </italic>is probably due to duplications of the genes in these regions. One clone was mapped to a mitochondrial rRNA gene that is not included in the annotation of <italic>D. pseudoobscura</italic>. Three clones were missing in the <italic>D. pseudoobscura </italic>genome assembly but exist in the trace archive database [##UREF##2##21##] and have homologous sequences in <italic>D. melanogaster</italic>. Four clones were too short to be mapped unambiguously and two clones were anti-sense transcripts. For four clones of the six remaining ones, no prediction exists in either of the species despite sequence conservation. Hence, these genes are not novel to <italic>D. pseudoobscura</italic>. The remaining two clones are putatively novel genes as these sequences are located in genomic regions with no gene predicted in <italic>D. pseudoobscura </italic>and no orthologous region could be identified in <italic>D. melanogaster </italic>by a BLAST search. Both of these clones map to <italic>D. pseudoobscura </italic>chromosome arm XR. One of them matched a SAGE tag that was present in the male and female SAGE library [##REF##16816428##19##] and showed female-biased expression. This sequence is located 72-nucleotides away from a predicted protein coding gene GA23511 in <italic>D. pseudoobscura</italic>. However, TBLASTN search of this protein against other completely sequenced Drosophila genomes did not reveal any hit except in the sibling species <italic>D. persimilis</italic>. The other clone did not match to any of the SAGE tags and was located in a region of the genome where no gene is predicted.</p> <title>Chromosomal distribution of novel sex-biased genes</title> <p>Six of the eight GLGI sequences representing putatively novel genes are located on autosomes, while only two were located on the X-chromosome. Based on the number of predicted genes on X- and autosomes in the <italic>D. pseudoobscura </italic>annotation release 2.0, we tested for an overrepresentation of novel genes on the X-chromosome and did not find support for an enrichment of new genes on one chromosome (<italic>p </italic>= 0.72, Fisher's exact test). Even after accounting for the two novel genes identified by the modified GLGI approach, we still found no evidence for an enrichment of novel genes on the X-chromosome (<italic>p </italic>= 0.75, Fisher's exact test).</p> <title>Structure of the novel genes</title> <p>In order to obtain full-length transcripts of the novel genes identified by GLGI, we performed 5' RACE analysis. All the five male biased tags were successfully amplified with 5' RACE. Out of three female-biased clones, only one could be cloned and sequenced. Based on the 5' RACE analysis, we could show that the transcripts of these novel genes are very short, ranging from 138 base pairs to 442 base pairs, yielding very short coding sequences ranging from 29 amino acids to 95 amino acids (Figure ##FIG##0##1##). These transcripts have start and stop codons at the same position in the orthologous regions of the closely relatives <italic>D. persimilis </italic>and <italic>D. miranda </italic>and the complete ORF of these transcripts is conserved in the three species. Interestingly, the gene containing the shortest ORF was found to possess an intron. Nevertheless, further experiments are required to determine whether or not these short ORFs are translated into peptides. For one of the male-biased transcripts (SAGE_M_22_3), no potential coding sequence could be obtained and hence we assumed that it could be a non protein coding transcript.</p> <title>Cross-species conservation of novel genes</title> <p>We took advantage of the 12 sequenced Drosophila genomes and performed TBLASTN searches with default parameters provided by flybase [##UREF##3##22##] to identify orthologs of the novel genes. Apart from <italic>D. persimilis</italic>, a very close relative of <italic>D. pseudoobscura</italic>, we did not detect any sequence similarity in the other species. Using PCR primers conserved between <italic>D. pseudoobscura </italic>and <italic>D. persimilis </italic>(See Additional file ##SUPPL##0##1##, Table 2), we successfully amplified the novel genes in <italic>D. miranda</italic>. Successful amplification of orthologs of the new, putatively protein coding genes from cDNA indicated their conservation in <italic>D. miranda</italic>. This implies that the novel genes arose before the split of <italic>D. pseudoobscura </italic>and <italic>D. miranda</italic>, which occurred about 2–5 MYA [##REF##8913754##23##, ####REF##9383055##24##, ##REF##17555574##25####17555574##25##]. We assessed the protein coding potential of these short transcripts using the QRNA software, which distinguishes between coding alignments, conserved non-coding alignments and RNA secondary structure [##REF##11801179##26##]. For none of the six transcripts did we find support for a protein coding potential (Table ##TAB##1##2##). However, as the sequences were short and the divergence between <italic>D. miranda </italic>and <italic>D. pseudoobscura </italic>is low, it is not clear how powerful this method is. Furthermore, it has been noted previously that failure of QRNA to identify a coding potential does not exclude the possibility of small ORFs [##REF##11553332##27##]. Except for one (SAGE_M_79), all putatively protein coding genes had a ratio of synonymous to non-synonymous substitutions smaller than one, indicating purifying selection on a protein sequence. SAGE_M_79 contained only non-synonymous substitutions and no synonymous substitutions (Table ##TAB##1##2##). Supporting this hypothesis, the ratio of polymorphic non-synonymous substitutions to synonymous substitutions is also smaller than one for the polymorphic loci (Table ##TAB##2##3##). On the other hand, the likelihood ratio test (LRT), which tests if the dN/dS ratio is significantly lower than 1, was not significant for any of the loci except SAGE_M_71 (Table ##TAB##1##2##). However, it should be noted that the power of LRT is affected by sequence length, divergence and number of species [##REF##11470850##28##] and may lead to a lower power in our case.</p> <title>Evolution of novel genes</title> <p>Recently evolved novel genes were shown to be under adaptive evolution in <italic>D. melanogaster </italic>lineage species [##REF##16777968##8##,##REF##17435230##9##]. To understand the pattern of variation and to test the deviations from neutral expectations in the novel genes of <italic>D. pseudoobscura</italic>, we sequenced the protein coding genes in one <italic>D. pseudoobscura </italic>population from Mesa-Verde (Colorado, USA, n = 8). Neutrality tests like Tajima's D (Table ##TAB##2##3##), Fu &amp; Li's D, Fu &amp; Li's F and Fu's Fs did not show a significant deviation from neutral expectations. We also performed a McDonald-Kreitman test [##REF##1904993##29##] using <italic>D. miranda</italic>, which basically assumes that under neutrality, the ratio of fixed differences between species to polymorphic differences within species should be similar for both the synonymous and non-synonymous sites. None of the novel loci showed significant deviations from neutral expectations (Table ##TAB##3##4##). It should be noted, however, that due to the short sequence of the loci analyzed the statistical power to detect selection is limited. A multi locus HKA test [##REF##3110004##30##], which tests for the deviations between the polymorphism within species versus divergence between species across different loci was also not significant, suggesting that none of the genes deviates significantly from the remaining ones.</p>
<title>Discussion</title> <p>In this report we showed that the extension of SAGE tags by GLGI is a powerful approach for the identification of novel genes.</p> <title>Origin of the novel genes</title> <p>Novel genes can be generated through different evolutionary trajectories: retrotransposition [##REF##11961099##31##], gene/genome duplication [##UREF##0##5##,##REF##15064762##32##] and exon shuffling [##REF##10570989##7##]. All these processes have in common that novel genes are built from existing genes or exons. Hence, it is expected that these building blocks should be detectable in the genome. We performed BLAST search of the novel genes against the <italic>D. pseudoobscura </italic>genome and detected only for one gene (SAGE_M_310) two BLAST hits. These two hits with complete sequence conservation were separated by 4.2-kilobase and the entire duplicated region spans approximately 1.1-kilobase with 92% identity. While we cannot rule out an assembly error, this observation suggests a recent duplication of the entire region encompassing the novel gene. Nevertheless, it is also apparent that both copies qualify as novel genes by our criteria.</p> <p>Given that we lack support for the origin of the described novel genes from already existing ones, we favour the hypothesis that they are derived from previously non-coding regions. Recently, several cases of such novel genes derived from non-coding sequences were described in other Drosophila species [##REF##16777968##8##,##REF##17435230##9##]. Like in our study, the functional evidence came from gene expression in the focal species. In order to obtain additional support for functionality of these transcripts, we also tested for gene expression in a related species and found all of the novel genes also to be transcribed in <italic>D. miranda</italic>. But, we cannot exclude that these transcripts are non-coding or the regions are spuriously expressed in the closely related species. However, such small ORFs were previously reported in Drosophila and other species and are known to be functional [##REF##17486114##33##, ####REF##16901788##34##, ##REF##17439302##35##, ##REF##15809421##36####15809421##36##]. The generation of novel genes by mutation appears a very unlikely event and the chance of generating a functional gene by random mutations decreases with the length of a gene. Hence, it is particularly interesting that all the novel genes identified in this study and previous ones [##REF##16777968##8##,##REF##17435230##9##] tend to be very short. Once more novel genes become available in combination with information about their date of origin, it will be possible to infer the probability to generate novel genes by random mutations. In particular, we will gain more insight if some sequences are more prone to develop into functional genes.</p> <title>Novel genes are male-biased</title> <p>In this report we searched for novel genes among SAGE tags showing either a strong male or female expression bias. Despite that we screened a considerably larger number of female-biased SAGE tags, the majority of the novel genes had a male-biased gene expression pattern. This pattern is consistent with previous results [##REF##16777968##8##,##REF##7682012##12##,##REF##9859991##37##, ####REF##12871908##38##, ##REF##11279525##39####11279525##39##]. It appears to be a general trend that many novel genes possess functional significance in males relative to females, but the reason for this trend is still not completely clear. Possible reasons include a higher transcription rate in male germline, greater functional pleiotropy of genes expressed in females and/or sexual competition [##REF##17994089##40##]. Alternatively, there could be more testis specific promoter sequences in the genome.</p> <title>Non-preferential X chromosomal location of novel genes</title> <p>There has been considerable controversy about the genomic location of male-biased genes. Earlier theoretical predictions [##UREF##4##41##] as well as expression studies in mice [##REF##11279525##39##] suggested that recessive genes conferring an advantage to males should be located on the X-chromosome. Expression studies in Drosophila [##REF##12511656##42##] and <italic>C. elegans </italic>[##REF##11030340##43##] showed, however, that male biased genes are preferentially located on the autosomes. This discrepancy could be explained by X-chromosome inactivation during spermatogenesis [##REF##4621547##44##, ####REF##11807039##45##, ##REF##12711214##46####12711214##46##]. The preferential X-linkage of novel genes, which show a male expression pattern [##REF##16777968##8##] requires further explanation. It has been suggested that genes conferring an advantage to males first originate on the X-chromosome and move to an autosomal location later in evolution [##REF##12511656##42##,##REF##12560539##47##]. Hence, it is conceivable that these novel male-biased genes do not serve functions that are essential during male spermatogenesis. Once such novel genes become essential, they could move to the autosomes. It was also suggested that mutations generating <italic>de novo </italic>genes might occur more often on the X chromosome or fix more readily [##REF##16777968##8##]. The novel genes we identified in this analysis can be dated from a minimum of 2–5 MYA (<italic>D. pseudoobscura </italic>and <italic>D. miranda </italic>divergence time [##REF##8913754##23##, ####REF##9383055##24##, ##REF##17555574##25####17555574##25##]) to a maximum of 13–15 MYA (divergence of obscura group [##REF##17555574##25##,##REF##1630303##48##]). It is possible that the novel genes identified in this study may be older compared to the novel genes reported by Levine et al. [##REF##16777968##8##], which was 2.5 MYA and if this is true, then these genes might have already moved to autosomes [##REF##12511656##42##]. More data are required to see if the discrepancy between our study and the results of Levine et al. [##REF##16777968##8##] are related to the difference in age or could be simply attributed to sampling effects due to the small number of genes in both studies.</p> <title>High incidence of novel genes in <italic>D. pseudoobscura</italic></title> <p>The discovery of previously uncharacterized transcripts is an observation common to many SAGE and EST sequencing experiments [##REF##12213963##11##,##REF##11279525##39##,##REF##12593804##49##,##REF##12466851##50##]. Nevertheless, the identification of evolutionary novelties requires also the absence of sequence conservation in related species. This approach was pioneered by Schmid and Tautz [##REF##9275195##51##], who studied the conservation of cDNAs across various species using hybridization to genomic Southern blots. More recently, Levine et al. [##REF##16777968##8##] specifically searched for lineage specific genes by BLAST search of all cDNA sequences available in <italic>D. melanogaster </italic>against <italic>D. yakuba</italic>, <italic>D. erecta </italic>and <italic>D. ananassae</italic>. A similar approach was also pursued for <italic>D. yakuba </italic>by BLAST search of ESTs against the genomic sequences of <italic>D. yakuba</italic>, <italic>D. melanogaster</italic>, <italic>D. erecta </italic>and <italic>D. ananassae </italic>[##REF##17435230##9##]. In this study, we focused on novel genes in <italic>D. pseudoobscura</italic>. Our previous SAGE analysis indicated that about 12% of the genes surveyed showed a significant sex bias. Of these, 73% could not be mapped to the <italic>D. pseudoobscura </italic>transcriptome and the majority of them possessed female-biased gene expression. This discrepancy could arise from a high number of novel female-biased genes in <italic>D. pseudoobscura</italic>. Contrary to our expectations and in agreement with the other studies we observed a high number of male-biased novel genes compared to female-biased novel genes. This result could be the outcome of different rates of evolution of male and female-biased genes [##REF##16612543##15##,##REF##18039869##18##] in <italic>D. pseudoobscura</italic>.</p> <p>In this study, we showed that 18% of the unmapped SAGE tags originate from putatively novel genes that arose in the <italic>obscura </italic>group. Hence, among the sex-biased genes approximately 2% are novel. In a cDNA library from <italic>D. pseudoobscura </italic>females we observed about 4% novel genes. This high incidence of novel genes in <italic>D. pseudoobscura </italic>underlines the need to supplement available genomic sequences with a thorough characterization of their transcriptome. In wake of the recent advances in the sequencing technology, we anticipate that this is already in close reach.</p>
<title>Conclusion</title> <p>We identified eight novel genes with sex-biased gene expression in <italic>D. pseudoobscura </italic>from unmapped SAGE tags using the GLGI method. In agreement with previous results from the other Drosophila species, a majority of these novel genes are male-biased in gene expression. Interestingly, we found no significant excess of X-linked novel genes. Overall, our data show that a considerable number of novel sex-biased genes exists in <italic>D. pseudoobscura </italic>that could not be predicted based on <italic>D. melanogaster </italic>gene models, which underlines the need to supplement available genomic sequences with a thorough characterization of their transcriptome.</p>
<title>Background</title> <p>The origin of functional innovation is among the key questions in biology. Recently, it has been shown that new genes could arise from non-coding DNA and that such novel genes are often involved in male reproduction.</p> <title>Results</title> <p>With the aim of identifying novel genes, we used the technique "generation of longer cDNA fragments from serial analysis of gene expression (SAGE) tags for gene identification (GLGI)" to extend 84 sex-biased 3'end SAGE tags that previously could not be mapped to the <italic>D. pseudoobscura </italic>transcriptome. Eleven male-biased and 33 female-biased GLGI fragments were obtained, of which 5 male-biased and 3 female-biased tags corresponded to putatively novel genes. This excess of novel genes with a male-biased gene expression pattern is consistent with previous results, which found novel genes to be primarily expressed in male reproductive tissues. 5' RACE analysis indicated that these novel transcripts are very short in length and could contain introns. Interspecies comparisons revealed that most novel transcripts show evidence for purifying selection.</p> <title>Conclusion</title> <p>Overall, our data indicate that among sex-biased genes a considerable number of novel genes (~2–4%) exist in <italic>D. pseudoobscura</italic>, which could not be predicted based on <italic>D. melanogaster </italic>gene models.</p>
<title>Abbreviations</title> <p>SAGE: serial analysis of gene expression; GLGI: generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification; RACE: rapid amplification of cDNA ends; UTR: untranslated region; ORF: open reading frame; MYA: million years ago.</p> <title>Authors' contributions</title> <p>MM carried out the experiments, analyzed the data and drafted the preliminary manuscript. CS conceived of the study, and participated in its design and coordination and drafted the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank V. Nolte for assistance with the GLGI experiments. We are thankful to the CS lab members for continuous support and discussion. Special thanks to S. Schaeffer for sharing the <italic>D. pseudoobscura </italic>population. This work has been supported by Fonds zur Förderung der wissenschaftlichen Forschung (FWF) grants (No. P17005, P17373) to CS. MM is supported by a Doktoratsstipendium of the Veterinärmedizinischen Universität Wien.</p>
<fig id="F1" position="float"><label>Figure 1</label><caption><p><bold>Novel genes in <italic>D. pseudoobscura</italic></bold>. List of novel genes obtained by analyzing SAGE, GLGI and 5'RACE with their protein length and chromosomal location. Black boxes represent coding regions, grey boxes represent non-coding regions and UTRs and white boxes represent introns. For the gene SAGE_M_71 (*), the chromosome location is determined based on its synteny in <italic>D. persimilis</italic>.</p></caption></fig>
<table-wrap id="T1" position="float"><label>Table 1</label><caption><p>Overall statistics of tags used in this analysis</p></caption><table frame="hsides" rules="groups"><thead><tr><th/><th align="center">Male biased</th><th align="center">Female biased</th></tr></thead><tbody><tr><td align="left">Total tags for GLGI analysis</td><td align="center">20</td><td align="center">64</td></tr><tr><td align="left">Successfully amplified tags</td><td align="center">11</td><td align="center">33</td></tr><tr><td align="left">Putative novel genes</td><td align="center">5</td><td align="center">3</td></tr><tr><td align="left">Successfully amplified novel genes using 5'RACE</td><td align="center">5</td><td align="center">1</td></tr></tbody></table></table-wrap> <table-wrap id="T2" position="float"><label>Table 2</label><caption><p>Synonymous (dS) and non-synonymous (dN) ratios between <italic>D. pseudoobscura </italic>and <italic>D. miranda </italic>and QRNA predictions of the novel genes</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Gene</th><th align="center">dN</th><th align="center">dS</th><th align="center">dN/dS</th><th align="center"><bold>LRT</bold><sup>#</sup></th><th align="center">QRNA*</th></tr></thead><tbody><tr><td align="left">SAGE_M_99</td><td align="center">0.029</td><td align="center">0.063</td><td align="center">0.462</td><td align="center">0.570</td><td align="center">O287</td></tr><tr><td align="left">SAGE_M_79</td><td align="center">0.045</td><td align="center">0.000</td><td align="center">-</td><td align="center">0.208</td><td align="center">O328</td></tr><tr><td align="left">SAGE_M_310</td><td align="center">0.014</td><td align="center">0.036</td><td align="center">0.394</td><td align="center">0.270</td><td align="center">O202</td></tr><tr><td align="left">SAGE_M_71</td><td align="center">0.013</td><td align="center">0.134</td><td align="center">0.097</td><td align="center">0.049</td><td align="center">O154</td></tr><tr><td align="left">SAGE_M_22_3</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">O416</td></tr><tr><td align="left">SAGE_F_108</td><td align="center">0.024</td><td align="center">0.056</td><td align="center">0.430</td><td align="center">0.102</td><td align="center">R460</td></tr></tbody></table></table-wrap> <table-wrap id="T3" position="float"><label>Table 3</label><caption><p>Neutrality tests for the novel genes in Mesa-Verde population (n = 8)</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Gene</th><th align="center">H</th><th align="center">π</th><th align="center">Tajima's D*</th><th align="center">πN</th><th align="center">πS</th><th align="center">πN/πS</th></tr></thead><tbody><tr><td align="left">SAGE_M_71</td><td align="center">6</td><td align="center">0.012</td><td align="center">0.687</td><td align="center">0.0000</td><td align="center">0.0725</td><td align="center">0.000</td></tr><tr><td align="left">SAGE_M_79</td><td align="center">1</td><td align="center">0.000</td><td align="center">0.000</td><td align="center">0.0000</td><td align="center">0.0000</td><td align="center">0.000</td></tr><tr><td align="left">SAGE_M_99</td><td align="center">5</td><td align="center">0.006</td><td align="center">-0.503</td><td align="center">0.0061</td><td align="center">0.0153</td><td align="center">0.393</td></tr><tr><td align="left">SAGE_M_310</td><td align="center">1</td><td align="center">0.000</td><td align="center">0.000</td><td align="center">0.0000</td><td align="center">0.0000</td><td align="center">0.000</td></tr><tr><td align="left">SAGE_F_108</td><td align="center">7</td><td align="center">0.007</td><td align="center">-0.414</td><td align="center">0.0025</td><td align="center">0.0155</td><td align="center">0.163</td></tr></tbody></table></table-wrap> <table-wrap id="T4" position="float"><label>Table 4</label><caption><p>McDonald-Kreitman test contingency table for the novel genes</p></caption><table frame="hsides" rules="groups"><thead><tr><th align="left">Gene</th><th align="center" colspan="2">Fixed differences</th><th align="center" colspan="2">Polymorphic differences</th></tr><tr><th/><th colspan="2"><hr/></th><th colspan="2"><hr/></th></tr><tr><th/><th align="center">Synonymous</th><th align="center">Non-synonymous</th><th align="center">Synonymous</th><th align="center">Non-synonymous</th></tr></thead><tbody><tr><td align="left">SAGE_M_71</td><td align="center">2</td><td align="center">2</td><td align="center">0</td><td align="center">4</td></tr><tr><td align="left">SAGE_M_79</td><td align="center">0</td><td align="center">4</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">SAGE_M_99</td><td align="center">0</td><td align="center">2</td><td align="center">1</td><td align="center">1</td></tr><tr><td align="left">SAGE_M_310</td><td align="center">1</td><td align="center">1</td><td align="center">0</td><td align="center">0</td></tr><tr><td align="left">SAGE_F_108</td><td align="center">5</td><td align="center">4</td><td align="center">3</td><td align="center">1</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>1a – List of female-biased tags with successful GLGI amplification. 1b – List of male-biased tags with successful GLGI amplification</bold>. <bold>†</bold>The gene is not predicted in the database of <italic>D. pseudoobscura</italic>. "Within the gene" indicates that the tag is located within the transcript but as it belongs to a mitochondrial gene or an unpredicted gene, it was not mapped previously. "Splice form" indicates that the SAGE tag is falling within an intron of a gene and so it is likely to be a splice variant. "Anti-sense" indicates that the GLGI sequence is matched in the anti-sense direction of the predicted gene. <bold>*</bold>This tag did not give any hit in <italic>D. pseudoobscura </italic>genome assembly but exists in the trace archive. The chromosome location is determined based on the sequence synteny in <italic>D. persimilis</italic>. "Sequencing artefact" indicates that the tag is falling within the gene but the region is having no sequence coverage (represented by 'n' in the database). <bold>2 – Primers used for cross-species amplification as well as population analysis of <italic>D. pseudoobscura</italic>.</bold></p></caption></supplementary-material>
<table-wrap-foot><p><sup>#</sup>Likelihood ratio test was performed using <italic>D. pseudoobscura</italic>, <italic>D. persimilis </italic>and <italic>D. miranda</italic>.</p><p>* O – other conserved non-coding sequences; R – RNA secondary structure. The number after the letter indicates the number of nucleotides assigned by QRNA to this category, which is in all cases equivalent to length of the full transcript.</p></table-wrap-foot> <table-wrap-foot><p>H = Number of Haplotypes.</p><p>π = Nucleotide diversity.</p><p>* <italic>p </italic>&gt; 0.05</p></table-wrap-foot>
<graphic xlink:href="1471-2148-8-182-1"/>
<media xlink:href="1471-2148-8-182-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Ohno"], "given-names": ["S"], "source": ["Evolution by gene duplication."], "year": ["1970"], "publisher-name": ["New York, Springer-Verlag"]}, {"surname": ["Patthy"], "given-names": ["L"], "source": ["Protein Evolution by Exon-shuffling"], "year": ["1995"], "publisher-name": ["New York, Springer-Verlag"]}, {"article-title": ["The NCBI database"], "ext-link": ["http://www.ncbi.nlm.nih.gov"]}, {"article-title": ["Flybase"], "ext-link": ["http://www.flybase.org"]}, {"surname": ["Rice"], "given-names": ["WR"], "article-title": ["Sex-Chromosomes and the Evolution of Sexual Dimorphism"], "source": ["Evolution"], "year": ["1984"], "volume": ["38"], "fpage": ["735"], "lpage": ["742"], "pub-id": ["10.2307/2408385"]}, {"article-title": ["CodonCode Aligner"], "ext-link": ["http://www.codoncode.com"]}, {"article-title": ["HKA Program"], "ext-link": ["http://lifesci.rutgers.edu/~heylab/index.html"]}]
BMC Evol Biol. 2008 Jun 24; 8:182