ORFan Genes Research Articles

Distinct genome trichotomy in members of Hapalosiphonaceae is guided by habitat adaptation with Mastigocladus laminosus UU774 as a case study

The family Hapalosiphonaceae is grouped under subsection IV complying with International Code of Nomenclature for algae, fungi, and plants (ICN). This is a complex taxonomic group with often overlapping morphological traits. The members of this group have heterogeneous habitats and need urgent attention for de-entangling the taxon. The genome comparisons among 44 high quality drafts belonging to Hapalosiphonaceae suggests the presence of three distinct groups, one largely dominated by thermophiles (Group 1) and the other comprising of members of diverse habitats (Group 2). The third group is small with only 2 representative members (Group 3) are considered as outliers. The group containing members of diverse habitats are characterized by the presence of higher number of ORFan genes, mobilome and secondary metabolites genes. The thermophiles belonging to Group 1 are characterized by compact genome whereas the Group 2 genomes having diverse ecological habitat show genome expansion possibly due to recent acquisition. Comparison of synonymous and non-synonymous substitution of conserved genes suggest, members of Group 2 are undergoing relaxed purifying selection. There is a significant increase in the number of core genes (2075 for n = 35) when Group 1 organisms are considered alone compared to only 46 core genes when all three groups are taken together (n = 44). The presence of a huge pan genome (47,756) and a relatively tiny core genome (46) suggests that the taxon is an open and expanding. Members of Group 2 are distinct in coding for commercial metabolite shinorine that was visibly absent from other members of the Hapalosiphonaceae. Our RNAseq investigations on Mastigocladus laminosus UU774, a candidate of Group 2, showed that metabolic clusters overexpressed during stress. UU774 metabolic extracts shown potential antibacterial and anticancer cell line activity.

CURE on yeast genes of unknown function increases students' bioinformatics proficiency and research confidence.

Course-based undergraduate research experiences (CUREs) can reduce barriers to research opportunities while increasing student knowledge and confidence. However, the number of widely adopted, easily transferable CUREs is relatively small. Here, we describe a CURE aimed at determining the function of poorly characterized Saccharomyces cerevisiae genes. More than 20 years after sequencing of the yeast genome, nearly 10% of open reading frames (ORFs) still have at least one uncharacterized Gene Ontology (GO) term. We refer to these genes as "ORFans" and formed a consortium aimed at assigning functions to them. Specifically, over 70 faculty members attended summer workshops to learn the bioinformatics workflow and basic laboratory techniques described herein. Ultimately, this CURE was adapted for implementation at 34 institutions, resulting in over 1,300 students conducting course-based research on ORFans. Pre-/post-tests confirmed that students gained both (i) an understanding of gene ontology and (ii) knowledge regarding the use of bioinformatics to assign gene function. After using these data to craft their own hypotheses, then testing their predictions by constructing and phenotyping deletion strains, students self-reported significant gains in several areas, including computer modeling and exposure to a project where no one knows the outcome. Interestingly, most net gains self-reported by ORFan Gene Project participants were greater than published findings for CUREs assessed with the same survey instrument. The surprisingly strong impact of this CURE may be due to the incoming lack of experience of ORFan Project participants and/or the independent thought required to develop testable hypotheses from complex data sets.

Targeted Genome Editing of Virulent Pseudomonas Phages Using CRISPR-Cas3.

The vast number of unknown phage-encoded ORFan genes and limited insights into the genome organization of phages illustrate the need for efficient genome engineering tools to study bacteriophage genes in their natural context. In addition, there is an application-driven desire to alter phage properties, which is hampered by time constraints for phage genome engineering in the bacterial host. We here describe an optimized CRISPR-Cas3 system in Pseudomonas for straightforward editing of the genome of virulent bacteriophages. The two-vector system combines a broad host range CRISPR-Cas3 targeting plasmid with a SEVA plasmid for homologous directed repair, which enables the creation of clean deletions, insertions, or substitutions in the phage genome within a week. After creating the two plasmids separately, a co-transformation to P. aeruginosa cells is performed. A subsequent infection with the targeted phage allows the CRISPR-Cas3 system to cut the DNA specifically and facilitate or select for homologous recombination. This system has also been successfully applied for P. aeruginosa and Pseudomonas putida genome engineering. The method is straightforward, efficient, and universal, enabling to extrapolate the system to other phage-host pairs.

ORFans in Mitochondrial Genomes of Marine Polychaete Polydora.

Most characterized metazoan mitochondrial genomes are compact and encode a small set of proteins that are essential for oxidative phosphorylation, as well as rRNA and tRNA for their expression. However, in rare cases, invertebrate taxa have additional open reading frames (ORFs) in their mtDNA sequences. Here, we sequenced and analyzed the mitochondrial genome of a polychaete worm, Polydora cf. ciliata, part of whose life cycle takes place in low-oxygen conditions. In the mitogenome, we found three "ORFan" regions (544, 1,060, and 427 bp) that have no resemblance to any standard metazoan mtDNA gene but lack stop codons in one of the reading frames. Similar regions are found in the mitochondrial genomes of three other Polydora species and Bocardiella hamata. All five species share the same gene order in their mitogenomes, which differ from that of other known Spionidae mitogenomes. By analyzing the ORFan sequences, we found that they are under purifying selection pressure and contain conservative regions. The codon adaptation indices (CAIs) of the ORFan genes were in the same range of values as the CAI of conventional protein-coding genes in corresponding mitochondrial genomes. The analysis of the P. cf. ciliata mitochondrial transcriptome showed that ORFan-544, ORFan-427, and a portion of the ORFan-1060 are transcribed. Together, this suggests that ORFan-544 and ORFan-427 encode functional proteins. It is likely that the ORFans originated when the Polydora/Bocardiella species complex separated from the rest of the Spionidae, and this event coincided with massive gene rearrangements in their mitochondrial genomes and tRNA-Met duplication.

Using unique ORFan genes as strain-specific identifiers for Escherichia coli

BackgroundBacterial identification at the strain level is a much-needed, but arduous and challenging task. This study aimed to develop a method for identifying and differentiating individual strains among multiple strains of the same bacterial species. The set used for testing the method consisted of 17 Escherichia coli strains picked from a collection of strains isolated in Germany, Spain, the United Kingdom and Vietnam from humans, cattle, swine, wild boars, and chickens. We targeted unique or rare ORFan genes to address the problem of selective and specific strain identification. These ORFan genes, exclusive to each strain, served as templates for developing strain-specific primers.ResultsMost of the experimental strains (14 out of 17) possessed unique ORFan genes that were used to develop strain-specific primers. The remaining three strains were identified by combining a PCR for a rare gene with a selection step for isolating the experimental strains. Multiplex PCR allowed the successful identification of the strains both in vitro in spiked faecal material in addition to in vivo after experimental infections of pigs and recovery of bacteria from faecal material. In addition, primers for qPCR were also developed and quantitative readout from faecal samples after experimental infection was also possible.ConclusionsThe method described in this manuscript using strain-specific unique genes to identify single strains in a mixture of strains proved itself efficient and reliable in detecting and following individual strains both in vitro and in vivo, representing a fast and inexpensive alternative to more costly methods.

Saccharomyces cerevisiae YIL164C/NIT1 is a Putative Nitrilase

While the yeast Saccharomyces cerevisiae is a well‐studied organism with its genome fully sequenced two decades ago, there remain several hundred uncharacterized open reading frames, ORFans, with no significant sequence homology to other genes. As ORFan roles cannot be determined solely through sequence homology, their functions remain a mystery until both in vitroand in vivobiochemical and cell biology experiments are performed. In collaboration with the Yeast ORFan Gene Project, we used bioinformatics modules to initially characterize one such ORFan, YIL165C.We found that together with the gene immediately preceding it, YIL164C (NIT1), the two genes have significant homology to the nitrilase superfamily of proteins but have low sequence homology to the other two S. cerevisiae nitrilases NIT2 and NIT3. Interestingly, the YIL164C protein alone contains both the conserved catalytic triad and a region that shares homology with prokaryotic nitrilases that specifically target aromatic, rather than aliphatic, nitrile substrates. We additionally experimentally verified that lab strains of S. cerevisiae, such as BY4741, SEY6210, and GA74, have a bonafide stop codon breaking up the two proteins in clones from both genomic DNA and mRNA, while the brewing yeast strain WLP001 (White Labs) and other related yeast strains have a single nucleotide polymorphism (SNP; A to G transition) allowing these two proteins to be expressed as a single, continuous polypeptide chain. To understand the function of both the full‐length YIL164C/YIL165C and the truncated YIL164C protein, we expressed and purified the proteins in bacteria and began initial studies to understand their substrate preference. We hope to identify distinctions in activity, or the lack thereof, between the longer and the truncated proteins. Understanding these enzymes in vitro will provide clues to their in vivo functions and the evolutionary relationships among genes in this superfamily.

Transcriptional Organization of the Salmonella Typhimurium Phage P22 pid ORFan Locus.

Many phage genes lack sequence similarity to any other open reading frame (ORF) in current databases. These enigmatic ORFan genes can have a tremendous impact on phage propagation and host interactions but often remain experimentally unexplored. We previously revealed a novel interaction between phage P22 and its Salmonella Typhimurium host, instigated by the ORFan gene pid (for phage P22 encoded instigator of dgo expression) and resulting in derepression of the host dgoRKAT operon. The pid gene is highly expressed in phage carrier cells that harbor a polarly located P22 episome that segregates asymmetrically among daughter cells. Here, we discovered that the pid locus is fitted with a weak promoter, has an exceptionally long 5′ untranslated region that is instructive for a secondary pid mRNA species, and has a 3′ Rho-independent termination loop that is responsible for stability of the pid transcript.

Variability of mitochondrial ORFans hints at possible differences in the system of doubly uniparental inheritance of mitochondria among families of freshwater mussels (Bivalvia: Unionida)

BackgroundSupernumerary ORFan genes (i.e., open reading frames without obvious homology to other genes) are present in the mitochondrial genomes of gonochoric freshwater mussels (Bivalvia: Unionida) showing doubly uniparental inheritance (DUI) of mitochondria. DUI is a system in which distinct female-transmitted and male-transmitted mitotypes coexist in a single species. In families Unionidae and Margaritiferidae, the transition from dioecy to hermaphroditism and the loss of DUI appear to be linked, and this event seems to affect the integrity of the ORFan genes. These observations led to the hypothesis that the ORFans have a role in DUI and/or sex determination. Complete mitochondrial genome sequences are however scarce for most families of freshwater mussels, therefore hindering a clear localization of DUI in the various lineages and a comprehensive understanding of the influence of the ORFans on DUI and sexual systems. Therefore, we sequenced and characterized eleven new mitogenomes from poorly sampled freshwater mussel families to gather information on the evolution and variability of the ORFan genes and their protein products.ResultsWe obtained ten complete plus one almost complete mitogenome sequence from ten representative species (gonochoric and hermaphroditic) of families Margaritiferidae, Hyriidae, Mulleriidae, and Iridinidae. ORFan genes are present only in DUI species from Margaritiferidae and Hyriidae, while non-DUI species from Hyriidae, Iridinidae, and Mulleriidae lack them completely, independently of their sexual system. Comparisons among the proteins translated from the newly characterized ORFans and already known ones provide evidence of conserved structures, as well as family-specific features.ConclusionsThe ORFan proteins show a comparable organization of secondary structures among different families of freshwater mussels, which supports a conserved physiological role, but also have distinctive family-specific features. Given this latter observation and the fact that the ORFans can be either highly mutated or completely absent in species that secondarily lost DUI depending on their respective family, we hypothesize that some aspects of the connection among ORFans, sexual systems, and DUI may differ in the various lineages of unionids.

'Drc', a structurally novel ssDNA-binding transcription regulator of N4-related bacterial viruses.

Bacterial viruses encode a vast number of ORFan genes that lack similarity to any other known proteins. Here, we present a 2.20 Å crystal structure of N4-related Pseudomonas virus LUZ7 ORFan gp14, and elucidate its function. We demonstrate that gp14, termed here as Drc (ssDNA-binding RNA Polymerase Cofactor), preferentially binds single-stranded DNA, yet contains a structural fold distinct from other ssDNA-binding proteins (SSBs). By comparison with other SSB folds and creation of truncation and amino acid substitution mutants, we provide the first evidence for the binding mechanism of this unique fold. From a biological perspective, Drc interacts with the phage-encoded RNA Polymerase complex (RNAPII), implying a functional role as an SSB required for the transition from early to middle gene transcription during phage infection. Similar to the coliphage N4 gp2 protein, Drc likely binds locally unwound middle promoters and recruits the phage RNA polymerase. However, unlike gp2, Drc does not seem to need an additional cofactor for promoter melting. A comparison among N4-related phage genera highlights the evolutionary diversity of SSB proteins in an otherwise conserved transcription regulation mechanism.

The Complete Maternally and Paternally Inherited Mitochondrial Genomes of a Freshwater Mussel Potamilus alatus (Bivalvia: Unionidae).

Doubly uniparental inheritance (DUI) of mitochondrial DNA, found only in some bivalve families and characterized by the existence of gender-associated mtDNA lineages that are inherited through males (M-type) or females (F-type), is one of the very few exceptions to the general rule of strict maternal mtDNA inheritance in animals. M-type sequences are often undetected and hence still underrepresented in the GenBank, which hinders the progress of the understanding of the DUI phenomenon. We have sequenced and analyzed the complete M and F mitogenomes of a freshwater mussel, Potamilus alatus. The M-type was 493 bp longer (M = 16 560, F = 16 067 bp). Gene contents, order and the distribution of genes between L and H strands were typical for unionid mussels. Candidates for the two ORFan genes (forf and morf) were found in respective mitogenomes. Both mitogenomes had a very similar A+T bias: F = 61% and M = 62.2%. The M mitogenome-specific cox2 extension (144 bp) is much shorter than in other sequenced unionid mitogenomes (531–576 bp), which might be characteristic for the Potamilus genus. The overall topology of the phylogenetic tree is in very good agreement with the currently accepted phylogenetic relationships within the Unionidae: both studied sequences were placed within the Ambleminae subfamily clusters in the corresponding M and F clades.

Recombinant expression of Mimivirus L725 ORFan gene product.

Mimivirus was the first discovered amoebal giant virus. The Mimivirus virions are covered by a dense layer of approximately 130 nm-long fibers, the length and shape of which diverge from those of other viruses. Here, we aimed at expressing the L725 protein to further confirm and study its role as a fiber-associated protein. We report Escherichia coli expression of the L725 protein, which is encoded by a Mimivirus ORFan, was previously identified by proteomics in purified viral fibers and demonstrated to be a fiber-associated protein by RNA-silencing experiments. The expressed protein was recognized by anti-Mimivirus fiber or anti-Mimivirus L725 polyclonal antibodies. This study is the only expression, to our knowledge, of a product from a Mimiviral ORFan gene.

In silico analyses of mitochondrial ORFans in freshwater mussels (Bivalvia: Unionoida) provide a framework for future studies of their origin and function.

BackgroundMany species of bivalves exhibit a unique system of mtDNA transmission named Doubly Uniparental Inheritance (DUI). Under this system, species have two distinct, sex-linked mitochondrial genomes: the M-type mtDNA, which is transmitted by males to male offspring and found in spermatozoa, and the F-type mtDNA, which is transmitted by females to all offspring, and found in all tissues of females and in somatic tissues of males. Bivalves with DUI also have sex-specific mitochondrial ORFan genes, (M-orf in the M mtDNA, F-orf in the F mtDNA), which are open reading frames having no detectable homology and no known function. DUI ORFan proteins have previously been characterized in silico in a taxonomically broad array of bivalves including four mytiloid, one veneroid and one unionoid species. However, the large evolutionary distance among these taxa prevented a meaningful comparison of ORFan properties among these divergent lineages. The present in silico study focuses on a suite of more closely-related Unionoid freshwater mussel species to provide more reliably interpretable information on patterns of conservation and properties of DUI ORFans. Unionoid species typically have separate sexes, but hermaphroditism also occurs, and hermaphroditic species lack the M-type mtDNA and possess a highly mutated version of the F-orf in their maternally transmitted mtDNA (named H-orf in these taxa). In this study, H-orfs and their respective proteins are analysed for the first time.ResultsDespite a rapid rate of evolution, strong structural and functional similarities were found for M-ORF proteins compared among species, and among the F-ORF and H-ORF proteins across the studied species. In silico analyses suggest that M-ORFs have a role in transport and cellular processes such as signalling, cell cycle and division, and cytoskeleton organisation, and that F-ORFs may be involved in cellular traffic and transport, and in immune response. H-ORFs appear to be structural glycoproteins, which may be involved in signalling, transport and transcription. Our results also support either a viral or a mitochondrial origin for the ORFans.ConclusionsOur findings reveal striking structural and functional similarities among proteins encoded by mitochondrial ORFans in freshwater mussels, and strongly support a role for these genes in the DUI mechanism. Our analyses also support the possibility of DUI systems with elements of different sources/origins and different mechanisms of action in the distantly-related DUI taxa. Parallel situations to the novel mitochondrially-encoded functions of freshwater mussel ORFans present in some other eukaryotes are also discussed.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2986-6) contains supplementary material, which is available to authorized users.

Identification of giant Mimivirus protein functions using RNA interference.

Genomic analysis of giant viruses, such as Mimivirus, has revealed that more than half of the putative genes have no known functions (ORFans). We knocked down Mimivirus genes using short interfering RNA as a proof of concept to determine the functions of giant virus ORFans. As fibers are easy to observe, we targeted a gene encoding a protein absent in a Mimivirus mutant devoid of fibers as well as three genes encoding products identified in a protein concentrate of fibers, including one ORFan and one gene of unknown function. We found that knocking down these four genes was associated with depletion or modification of the fibers. Our strategy of silencing ORFan genes in giant viruses opens a way to identify its complete gene repertoire and may clarify the role of these genes, differentiating between junk DNA and truly used genes. Using this strategy, we were able to annotate four proteins in Mimivirus and 30 homologous proteins in other giant viruses. In addition, we were able to annotate >500 proteins from cellular organisms and 100 from metagenomic databases.

55.2, a phage T4 ORFan gene, encodes an inhibitor of Escherichia coli topoisomerase I and increases phage fitness.

Topoisomerases are enzymes that alter the topological properties of DNA. Phage T4 encodes its own topoisomerase but it can also utilize host-encoded topoisomerases. Here we characterized 55.2, a phage T4 predicted ORF of unknown function. High levels of expression of the cloned 55.2 gene are toxic in E. coli. This toxicity is suppressed either by increased topoisomerase I expression or by partial inactivation of the ATPase subunit of the DNA gyrase. Interestingly, very low-level expression of 55.2, which is non-lethal to wild type E. coli, prevents the growth of a deletion mutant of the topoisomerase I (topA) gene. In vitro, gp55.2 binds DNA and blocks specifically the relaxation of negatively supercoiled DNA by topoisomerase I. In vivo, expression of gp55.2 at low non-toxic levels alters the steady state DNA supercoiling of a reporter plasmid. Although 55.2 is not an essential gene, competition experiments indicate that it is required for optimal phage growth. We propose that the role of gp55.2 is to subtly modulate host topoisomerase I activity during infection to insure optimal T4 phage yield.

Strong seed-bank effects in bacterial evolution

Bacterial genomes are mosaics with fragments showing distinct phylogenetic origins or even being unrelated to any other genetic information (ORFan genes). Thus the analysis of bacterial population genetics is in large part a collection of explanations for anomalies in relation to classical population genetic models such as the Wright-Fisher model and the Kingman coalescent that do not adequately describe bacterial population genetics, genomics or evolution. The concept of “species” as an evolutionary coherent biological group that is genetically isolated and shares genetic information through recombination among its members cannot be applied to any bacterial group. Recently, a simple probabilistic model considering the role of strong seed-bank effects in population genetics has been proposed by Blath et al. This model suggests the existence of a genetic pool with high diversity that is not subject to classical selection and extinction. We reason that certain bacterial population genetics anomalies could be explained by the prevalence of strong seed-bank effects among bacteria. To address this possibility we analyzed the genome of the bacterium Azotobacter vinelandii and show that genes that code for functions that are essential for the bacterium biology do not have a relation of ancestry with closely related bacteria, or are ORFan genes. The existence of essential genes that are not inherited from the most recent ancestor cannot be explained by classical population genetics models and is irreconcilable with the current view of genes acquired by horizontal transfer as being accessory or adaptive.

Evidence for Moonlighting Functions of the Subunit of Escherichia coli DNA Polymerase III

The holE gene is an enterobacterial ORFan gene (open reading frame [ORF] with no detectable homology to other ORFs in a database). It encodes the θ subunit of the DNA polymerase III core complex. The precise function of the θ subunit within this complex is not well established, and loss of holE does not result in a noticeable phenotype. Paralogs of holE are also present on many conjugative plasmids and on phage P1 (hot gene). In this study, we provide evidence indicating that θ (HolE) exhibits structural and functional similarities to a family of nucleoid-associated regulatory proteins, the Hha/YdgT-like proteins that are also encoded by enterobacterial ORFan genes. Microarray studies comparing the transcriptional profiles of Escherichia coli holE, hha, and ydgT mutants revealed highly similar expression patterns for strains harboring holE and ydgT alleles. Among the genes differentially regulated in both mutants were genes of the tryptophanase (tna) operon. The tna operon consists of a transcribed leader region, tnaL, and two structural genes, tnaA and tnaB. Further experiments with transcriptional lacZ fusions (tnaL::lacZ and tnaA::lacZ) indicate that HolE and YdgT downregulate expression of the tna operon by possibly increasing the level of Rho-dependent transcription termination at the tna operon's leader region. Thus, for the first time, a regulatory function can be attributed to HolE, in addition to its role as structural component of the DNA polymerase III complex.

Expression of a novel P22 ORFan gene reveals the phage carrier state in Salmonella typhimurium.

We discovered a novel interaction between phage P22 and its host Salmonella Typhimurium LT2 that is characterized by a phage mediated and targeted derepression of the host dgo operon. Upon further investigation, this interaction was found to be instigated by an ORFan gene (designated pid for phage P22 encoded instigator of dgo expression) located on a previously unannotated moron locus in the late region of the P22 genome, and encoding an 86 amino acid protein of 9.3 kDa. Surprisingly, the Pid/dgo interaction was not observed during strict lytic or lysogenic proliferation of P22, and expression of pid was instead found to arise in cells that upon infection stably maintained an unintegrated phage chromosome that segregated asymmetrically upon subsequent cell divisions. Interestingly, among the emerging siblings, the feature of pid expression remained tightly linked to the cell inheriting this phage carrier state and became quenched in the other. As such, this study is the first to reveal molecular and genetic markers authenticating pseudolysogenic development, thereby exposing a novel mechanism, timing, and populational distribution in the realm of phage–host interactions.

A Bacteriophage-Encoded J-Domain Protein Interacts with the DnaK/Hsp70 Chaperone and Stabilizes the Heat-Shock Factor σ32 of Escherichia coli

The universally conserved J-domain proteins (JDPs) are obligate cochaperone partners of the Hsp70 (DnaK) chaperone. They stimulate Hsp70's ATPase activity, facilitate substrate delivery, and confer specific cellular localization to Hsp70. In this work, we have identified and characterized the first functional JDP protein encoded by a bacteriophage. Specifically, we show that the ORFan gene 057w of the T4-related enterobacteriophage RB43 encodes a bona fide JDP protein, named Rki, which specifically interacts with the Escherichia coli host multifunctional DnaK chaperone. However, in sharp contrast with the three known host JDP cochaperones of DnaK encoded by E. coli, Rki does not act as a generic cochaperone in vivo or in vitro. Expression of Rki alone is highly toxic for wild-type E. coli, but toxicity is abolished in the absence of endogenous DnaK or when the conserved J-domain of Rki is mutated. Further in vivo analyses revealed that Rki is expressed early after infection by RB43 and that deletion of the rki gene significantly impairs RB43 proliferation. Furthermore, we show that mutations in the host dnaK gene efficiently suppress the growth phenotype of the RB43 rki deletion mutant, thus indicating that Rki specifically interferes with DnaK cellular function. Finally, we show that the interaction of Rki with the host DnaK chaperone rapidly results in the stabilization of the heat-shock factor σ32, which is normally targeted for degradation by DnaK. The mechanism by which the Rki-dependent stabilization of σ32 facilitates RB43 bacteriophage proliferation is discussed.

Population Diversity of ORFan Genes in Escherichia coli

The origin and evolution of “ORFans” (suspected genes without known relatives) remain unclear. Here, we take advantage of a unique opportunity to examine the population diversity of thousands of ORFans, based on a collection of 35 complete genomes of isolates of Escherichia coli and Shigella (which is included phylogenetically within E. coli). As expected from previous studies, ORFans are shorter and AT-richer in sequence than non-ORFans. We find that ORFans often are very narrowly distributed: the most common pattern is for an ORFan to be found in only one genome. We compared within-species population diversity of ORFan genes with those of two control groups of non-ORFan genes. Patterns of population variation suggest that most ORFans are not artifacts, but encode real genes whose protein-coding capacity is conserved, reflecting selection against nonsynonymous mutations. Nevertheless, nonsynonymous nucleotide diversity is higher than for non-ORFans, whereas synonymous diversity is roughly the same. In particular, there is a several-fold excess of ORFans in the highest decile of diversity relative to controls, which might be due to weaker purifying selection, positive selection, or a subclass of ORFans that are decaying.

Composition bias and the origin of ORFan genes

Motivation: Intriguingly, sequence analysis of genomes reveals that a large number of genes are unique to each organism. The origin of these genes, termed ORFans, is not known. Here, we explore the origin of ORFan genes by defining a simple measure called ‘composition bias’, based on the deviation of the amino acid composition of a given sequence from the average composition of all proteins of a given genome.Results: For a set of 47 prokaryotic genomes, we show that the amino acid composition bias of real proteins, random ‘proteins’ (created by using the nucleotide frequencies of each genome) and ‘proteins’ translated from intergenic regions are distinct. For ORFans, we observed a correlation between their composition bias and their relative evolutionary age. Recent ORFan proteins have compositions more similar to those of random ‘proteins’, while the compositions of more ancient ORFan proteins are more similar to those of the set of all proteins of the organism. This observation is consistent with an evolutionary scenario wherein ORFan genes emerged and underwent a large number of random mutations and selection, eventually adapting to the composition preference of their organism over time.Contact: ron@biocoml.ls.biu.ac.ilSupplementary information: Supplementary data are available at Bioinformatics online.