The Planemo toolkit for developing, deploying, and executing scientific data analyses in Galaxy and beyond.

There are thousands of well-maintained high-quality open-source software utilities for all aspects of scientific data analysis. For more than a decade, the Galaxy Project has been providing computational infrastructure and a unified user interface for these tools to make them accessible to a wide range of researchers. To streamline the process of integrating tools and constructing workflows as much as possible, we have developed Planemo, a software development kit for tool and workflow developers and Galaxy power users. Here we outline Planemo's implementation and describe its broad range of functionality for designing, testing, and executing Galaxy tools, workflows, and training material. In addition, we discuss the philosophy underlying Galaxy tool and workflow development, and how Planemo encourages the use of development best practices, such as test-driven development, by its users, including those who are not professional software developers.

The 2015 Bioinformatics Open Source Conference (BOSC 2015).

The Bioinformatics Open Source Conference (BOSC) is organized by the Open Bioinformatics Foundation (OBF), a nonprofit group dedicated to promoting the practice and philosophy of open source software development and open science within the biological research community. Since its inception in 2000, BOSC has provided bioinformatics developers with a forum for communicating the results of their latest efforts to the wider research community. BOSC offers a focused environment for developers and users to interact and share ideas about standards; software development practices; practical techniques for solving bioinformatics problems; and approaches that promote open science and sharing of data, results, and software. BOSC is run as a two-day special interest group (SIG) before the annual Intelligent Systems in Molecular Biology (ISMB) conference. BOSC 2015 took place in Dublin, Ireland, and was attended by over 125 people, about half of whom were first-time attendees. Session topics included "Data Science;" "Standards and Interoperability;" "Open Science and Reproducibility;" "Translational Bioinformatics;" "Visualization;" and "Bioinformatics Open Source Project Updates". In addition to two keynote talks and dozens of shorter talks chosen from submitted abstracts, BOSC 2015 included a panel, titled "Open Source, Open Door: Increasing Diversity in the Bioinformatics Open Source Community," that provided an opportunity for open discussion about ways to increase the diversity of participants in BOSC in particular, and in open source bioinformatics in general. The complete program of BOSC 2015 is available online at http://www.open-bio.org/wiki/BOSC_2015_Schedule.

Comparative transcriptomics and proteomics of three different aphid species identifies core and diverse effector sets.

BACKGROUND: Aphids are phloem-feeding insects that cause significant economic losses to agriculture worldwide. While feeding and probing these insects deliver molecules, called effectors, inside their host to enable infestation. The identification and characterization of these effectors from different species that vary in their host range is an important step in understanding the infestation success of aphids and aphid host range variation. This study employs a multi-disciplinary approach based on transcriptome sequencing and proteomics to identify and compare effector candidates from the broad host range aphid Myzus persicae (green peach aphid) (genotypes O, J and F), and narrow host range aphids Myzus cerasi (black cherry aphid) and Rhopalosiphum padi (bird-cherry oat aphid). RESULTS: Using a combination of aphid transcriptome sequencing on libraries derived from head versus body tissues as well as saliva proteomics we were able to predict candidate effectors repertoires from the different aphid species and genotypes. Among the identified conserved or core effector sets, we identified a significant number of previously identified aphid candidate effectors indicating these proteins may be involved in general infestation strategies. Moreover, we identified aphid candidate effector sequences that were specific to one species, which are interesting candidates for further validation and characterization with regards to species-specific functions during infestation. We assessed our candidate effector repertoires for evidence of positive selection, and identified 49 candidates with DN/DS ratios >1. We noted higher rates of DN/DS ratios in predicted aphid effectors than non-effectors. Whether this reflects positive selection due to co-evolution with host plants, or increased neofunctionalization upon gene duplication remains to be investigated. CONCLUSION: Our work provides a comprehensive overview of the candidate effector repertoires from three different aphid species with varying host ranges. Comparative analyses revealed candidate effectors that are most likely are involved in general aspects of infestation, whereas others, that are highly divergent, may be involved in specific processes important for certain aphid species. Insights into the overlap and differences in aphid effector repertoires are important in understanding how different species successfully infest different ranges of plant species.

The genome and genetics of a high oxidative stress tolerant Serratia sp. LCN16 isolated from the plant parasitic nematode Bursaphelenchus xylophilus.

BACKGROUND: Pine wilt disease (PWD) is a worldwide threat to pine forests, and is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus. Bacteria are known to be associated with PWN and may have an important role in PWD. Serratia sp. LCN16 is a PWN-associated bacterium, highly resistant to oxidative stress in vitro, and which beneficially contributes to the PWN survival under these conditions. Oxidative stress is generated as a part of the basal defense mechanism used by plants to combat pathogenic invasion. Here, we studied the biology of Serratia sp. LCN16 through genome analyses, and further investigated, using reverse genetics, the role of two genes directly involved in the neutralization of H2O2, namely the H2O2 transcriptional factor oxyR; and the H2O2-targeting enzyme, catalase katA. RESULTS: Serratia sp. LCN16 is phylogenetically most closely related to the phytosphere group of Serratia, which includes S. proteamaculans, S. grimessi and S. liquefaciens. Likewise, Serratia sp. LCN16 shares many features with endophytes (plant-associated bacteria), such as genes coding for plant polymer degrading enzymes, iron uptake/transport, siderophore and phytohormone synthesis, aromatic compound degradation and detoxification enzymes. OxyR and KatA are directly involved in the high tolerance to H2O2 of Serratia sp. LCN16. Under oxidative stress, Serratia sp. LCN16 expresses katA independently of OxyR in contrast with katG which is under positive regulation of OxyR. Serratia sp. LCN16 mutants for oxyR (oxyR::int(614)) and katA (katA::int(808)) were sensitive to H2O2 in relation with wild-type, and both failed to protect the PWN from H2O2-stress exposure. Moreover, both mutants showed different phenotypes in terms of biofilm production and swimming/swarming behaviors. CONCLUSIONS: This study provides new insights into the biology of PWN-associated bacteria Serratia sp. LCN16 and its extreme resistance to oxidative stress conditions, encouraging further research on the potential role of this bacterium in interaction with PWN in planta environment.

Evidence for an Opportunistic and Endophytic Lifestyle of the Bursaphelenchus xylophilus-Associated Bacteria Serratia marcescens PWN146 Isolated from Wilting Pinus pinaster.

Pine wilt disease (PWD) results from the interaction of three elements: the pathogenic nematode, Bursaphelenchus xylophilus; the insect-vector, Monochamus sp.; and the host tree, mostly Pinus species. Bacteria isolated from B. xylophilus may be a fourth element in this complex disease. However, the precise role of bacteria in this interaction is unclear as both plant-beneficial and as plant-pathogenic bacteria may be associated with PWD. Using whole genome sequencing and phenotypic characterization, we were able to investigate in more detail the genetic repertoire of Serratia marcescens PWN146, a bacterium associated with B. xylophilus. We show clear evidence that S. marcescens PWN146 is able to withstand and colonize the plant environment, without having any deleterious effects towards a susceptible host (Pinus thunbergii), B. xylophilus nor to the nematode model C. elegans. This bacterium is able to tolerate growth in presence of xenobiotic/organic compounds, and use phenylacetic acid as carbon source. Furthermore, we present a detailed list of S. marcescens PWN146 potentials to interfere with plant metabolism via hormonal pathways and/or nutritional acquisition, and to be competitive against other bacteria and/or fungi in terms of resource acquisition or production of antimicrobial compounds. Further investigation is required to understand the role of bacteria in PWD. We have now reinforced the theory that B. xylophilus-associated bacteria may have a plant origin.

Newly identified RNAs of raspberry leaf blotch virus encoding a related group of proteins.

Members of the genus Emaravirus, including Raspberry leaf blotch virus (RLBV), are enveloped plant viruses with segmented genomes of negative-strand RNA, although the complete genome complement for any of these viruses is not yet clear. Currently, wheat mosaic virus has the largest emaravirus genome comprising eight RNAs. Previously, we identified five genomic RNAs for RLBV; here, we identify a further three RNAs (RNA6-8). RNA6-8 encode proteins that have clear homologies to one another, but not to any other emaravirus proteins. The proteins self-interacted in yeast two-hybrid and bimolecular fluorescence complementation (BiFC) experiments, and the P8 protein interacted with the virus nucleocapsid protein (P3) using BiFC. Expression of two of the proteins (P6 and P7) using potato virus X led to an increase in virus titre and symptom severity, suggesting that these proteins may play a role in RLBV pathogenicity; however, using two different tests, RNA silencing suppression activity was not detected for any of the RLBV proteins encoded by RNA2-8.

Using Tablet for visual exploration of second-generation sequencing data.

The advent of second-generation sequencing (2GS) has provided a range of significant new challenges for the visualization of sequence assemblies. These include the large volume of data being generated, short-read lengths and different data types and data formats associated with the diversity of new sequencing technologies. This article illustrates how Tablet-a high-performance graphical viewer for visualization of 2GS assemblies and read mappings-plays an important role in the analysis of these data. We present Tablet, and through a selection of use cases, demonstrate its value in quality assurance and scientific discovery, through features such as whole-reference coverage overviews, variant highlighting, paired-end read mark-up, GFF3-based feature tracks and protein translations. We discuss the computing and visualization techniques utilized to provide a rich and responsive graphical environment that enables users to view a range of file formats with ease. Tablet installers can be freely downloaded from http://bioinf.hutton.ac.uk/tablet in 32 or 64-bit versions for Windows, OS X, Linux or Solaris. For further details on the Tablet, contact tablet@hutton.ac.uk.

Distinct circular single-stranded DNA viruses exist in different soil types.

The potential dependence of virus populations on soil types was examined by electron microscopy, and the total abundance of virus particles in four soil types was similar to that previously observed in soil samples. The four soil types examined differed in the relative abundances of four morphological groups of viruses. Machair, a unique type of coastal soil in western Scotland and Ireland, differed from the others tested in having a higher proportion of tailed bacteriophages. The other soils examined contained predominantly spherical and thin filamentous virus particles, but the Machair soil had a more even distribution of the virus types. As the first step in looking at differences in populations in detail, virus sequences from Machair and brown earth (agricultural pasture) soils were examined by metagenomic sequencing after enriching for circular Rep-encoding single-stranded DNA (ssDNA) (CRESS-DNA) virus genomes. Sequences from the family Microviridae (icosahedral viruses mainly infecting bacteria) of CRESS-DNA viruses were predominant in both soils. Phylogenetic analysis of Microviridae major coat protein sequences from the Machair viruses showed that they spanned most of the diversity of the subfamily Gokushovirinae, whose members mainly infect obligate intracellular parasites. The brown earth soil had a higher proportion of sequences that matched the morphologically similar family Circoviridae in BLAST searches. However, analysis of putative replicase proteins that were similar to those of viruses in the Circoviridae showed that they are a novel clade of Circoviridae-related CRESS-DNA viruses distinct from known Circoviridae genera. Different soils have substantially different taxonomic biodiversities even within ssDNA viruses, which may be driven by physicochemical factors.

A metagenetic approach to determine the diversity and distribution of cyst nematodes at the level of the country, the field and the individual.

Distinct populations of the potato cyst nematode (PCN) Globodera pallida exist in the UK that differ in their ability to overcome various sources of resistance. An efficient method for distinguishing between populations would enable pathogen-informed cultivar choice in the field. Science and Advice for Scottish Agriculture (SASA) annually undertake national DNA diagnostic tests to determine the presence of PCN in potato seed and ware land by extracting DNA from soil floats. These DNA samples provide a unique resource for monitoring the distribution of PCN and further interrogation of the diversity within species. We identify a region of mitochondrial DNA descriptive of three main groups of G. pallida present in the UK and adopt a metagenetic approach to the sequencing and analysis of all SASA samples simultaneously. Using this approach, we describe the distribution of G. pallida mitotypes across Scotland with field-scale resolution. Most fields contain a single mitotype, one-fifth contain a mix of mitotypes, and less than 3% contain all three mitotypes. Within mixed fields, we were able to quantify the relative abundance of each mitotype across an order of magnitude. Local areas within mixed fields are dominated by certain mitotypes and indicate towards a complex underlying 'pathoscape'. Finally, we assess mitotype distribution at the level of the individual cyst and provide evidence of 'hybrids'. This study provides a method for accurate, quantitative and high-throughput typing of up to one thousand fields simultaneously, while revealing novel insights into the national genetic variability of an economically important plant parasite.

Community-driven development for computational biology at Sprints, Hackathons and Codefests.

BACKGROUND: Computational biology comprises a wide range of technologies and approaches. Multiple technologies can be combined to create more powerful workflows if the individuals contributing the data or providing tools for its interpretation can find mutual understanding and consensus. Much conversation and joint investigation are required in order to identify and implement the best approaches. Traditionally, scientific conferences feature talks presenting novel technologies or insights, followed up by informal discussions during coffee breaks. In multi-institution collaborations, in order to reach agreement on implementation details or to transfer deeper insights in a technology and practical skills, a representative of one group typically visits the other. However, this does not scale well when the number of technologies or research groups is large. Conferences have responded to this issue by introducing Birds-of-a-Feather (BoF) sessions, which offer an opportunity for individuals with common interests to intensify their interaction. However, parallel BoF sessions often make it hard for participants to join multiple BoFs and find common ground between the different technologies, and BoFs are generally too short to allow time for participants to program together. RESULTS: This report summarises our experience with computational biology Codefests, Hackathons and Sprints, which are interactive developer meetings. They are structured to reduce the limitations of traditional scientific meetings described above by strengthening the interaction among peers and letting the participants determine the schedule and topics. These meetings are commonly run as loosely scheduled "unconferences" (self-organized identification of participants and topics for meetings) over at least two days, with early introductory talks to welcome and organize contributors, followed by intensive collaborative coding sessions. We summarise some prominent achievements of those meetings and describe differences in how these are organised, how their audience is addressed, and their outreach to their respective communities. CONCLUSIONS: Hackathons, Codefests and Sprints share a stimulating atmosphere that encourages participants to jointly brainstorm and tackle problems of shared interest in a self-driven proactive environment, as well as providing an opportunity for new participants to get involved in collaborative projects.

Genomic insights into the origin of parasitism in the emerging plant pathogen Bursaphelenchus xylophilus.

Bursaphelenchus xylophilus is the nematode responsible for a devastating epidemic of pine wilt disease in Asia and Europe, and represents a recent, independent origin of plant parasitism in nematodes, ecologically and taxonomically distinct from other nematodes for which genomic data is available. As well as being an important pathogen, the B. xylophilus genome thus provides a unique opportunity to study the evolution and mechanism of plant parasitism. Here, we present a high-quality draft genome sequence from an inbred line of B. xylophilus, and use this to investigate the biological basis of its complex ecology which combines fungal feeding, plant parasitic and insect-associated stages. We focus particularly on putative parasitism genes as well as those linked to other key biological processes and demonstrate that B. xylophilus is well endowed with RNA interference effectors, peptidergic neurotransmitters (including the first description of ins genes in a parasite) stress response and developmental genes and has a contracted set of chemosensory receptors. B. xylophilus has the largest number of digestive proteases known for any nematode and displays expanded families of lysosome pathway genes, ABC transporters and cytochrome P450 pathway genes. This expansion in digestive and detoxification proteins may reflect the unusual diversity in foods it exploits and environments it encounters during its life cycle. In addition, B. xylophilus possesses a unique complement of plant cell wall modifying proteins acquired by horizontal gene transfer, underscoring the impact of this process on the evolution of plant parasitism by nematodes. Together with the lack of proteins homologous to effectors from other plant parasitic nematodes, this confirms the distinctive molecular basis of plant parasitism in the Bursaphelenchus lineage. The genome sequence of B. xylophilus adds to the diversity of genomic data for nematodes, and will be an important resource in understanding the biology of this unusual parasite.

THAPBI PICT-a fast, cautious, and accurate metabarcoding analysis pipeline.

THAPBI PICT is an open source software pipeline for metabarcoding analysis of Illumina paired-end reads, including cases of multiplexing where more than one amplicon is amplified per DNA sample. Initially a Phytophthora ITS1 Classification Tool (PICT), we demonstrate using worked examples with our own and public data sets how, with appropriate primer settings and a custom database, it can be applied to other amplicons and organisms, and used for reanalysis of existing datasets. The core dataflow of the implementation is (i) data reduction to unique marker sequences, often called amplicon sequence variants (ASVs), (ii) dynamic thresholds for discarding low abundance sequences to remove noise and artifacts (rather than error correction by default), before (iii) classification using a curated reference database. The default classifier assigns a label to each query sequence based on a database match that is either perfect, or a single base pair edit away (substitution, deletion or insertion). Abundance thresholds for inclusion can be set by the user or automatically using per-batch negative or synthetic control samples. Output is designed for practical interpretation by non-specialists and includes a read report (ASVs with classification and counts per sample), sample report (samples with counts per species classification), and a topological graph of ASVs as nodes with short edit distances as edges. Source code available from https://github.com/peterjc/thapbi-pict/ with documentation including installation instructions.

Bio.Phylo: a unified toolkit for processing, analyzing and visualizing phylogenetic trees in Biopython.

BACKGROUND: Ongoing innovation in phylogenetics and evolutionary biology has been accompanied by a proliferation of software tools, data formats, analytical techniques and web servers. This brings with it the challenge of integrating phylogenetic and other related biological data found in a wide variety of formats, and underlines the need for reusable software that can read, manipulate and transform this information into the various forms required to build computational pipelines. RESULTS: We built a Python software library for working with phylogenetic data that is tightly integrated with Biopython, a broad-ranging toolkit for computational biology. Our library, Bio.Phylo, is highly interoperable with existing libraries, tools and standards, and is capable of parsing common file formats for phylogenetic trees, performing basic transformations and manipulations, attaching rich annotations, and visualizing trees. We unified the modules for working with the standard file formats Newick, NEXUS and phyloXML behind a consistent and simple API, providing a common set of functionality independent of the data source. CONCLUSIONS: Bio.Phylo meets a growing need in bioinformatics for working with heterogeneous types of phylogenetic data. By supporting interoperability with multiple file formats and leveraging existing Biopython features, this library simplifies the construction of phylogenetic workflows. We also provide examples of the benefits of building a community around a shared open-source project. Bio.Phylo is included with Biopython, available through the Biopython website, http://biopython.org.

Genome-wide identification of genes regulated by the Rcs phosphorelay system in Erwinia amylovora.

The exopolysaccharide amylovoran is one of the major pathogenicity factors in Erwinia amylovora, the causal agent of fire blight of apples and pears. We have previously demonstrated that the RcsBCD phosphorelay system is essential for virulence by controlling amylovoran biosynthesis. We have also found that the hybrid sensor kinase RcsC differentially regulates amylovoran production in vitro and in vivo. To further understand how the Rcs system regulates E. amylovora virulence gene expression, we conducted genome-wide microarray analyses to determine the regulons of RcsB and RcsC in liquid medium and on immature pear fruit. Array analyses identified a total of 648 genes differentially regulated by RcsCB in vitro and in vivo. Consistent with our previous findings, RcsB acts as a positive regulator in both conditions, while RcsC positively controls expression of amylovoran biosynthetic genes in vivo but negatively controls expression in vitro. Besides amylovoran biosynthesis and regulatory genes, cell-wall and cell-envelope (membrane) as well as regulatory genes were identified as the major components of the RcsBC regulon, including many novel genes. We have also demonstrated that transcripts of rcsA, rcsC, and rcsD genes but not the rcsB gene were up-regulated when bacterial cells were grown in minimal medium or following infection of pear fruits compared with those grown in Luria Bertani medium. Furthermore, using the genome of E. amylovora ATCC 49946, a hidden Markov model predicted 60 genes with a candidate RcsB binding site in the intergenic region, 28 of which were identified in the microarray assay. Based on these findings as well as previous reported data, a working model has been proposed to illustrate how the Rcs phosphorelay system regulates virulence gene expression in E. amylovora.

Identification and localisation of the NB-LRR gene family within the potato genome.

BACKGROUND: The potato genome sequence derived from the Solanum tuberosum Group Phureja clone DM1-3 516 R44 provides unparalleled insight into the genome composition and organisation of this important crop. A key class of genes that comprises the vast majority of plant resistance (R) genes contains a nucleotide-binding and leucine-rich repeat domain, and is collectively known as NB-LRRs. RESULTS: As part of an effort to accelerate the process of functional R gene isolation, we performed an amino acid motif based search of the annotated potato genome and identified 438 NB-LRR type genes among the ~39,000 potato gene models. Of the predicted genes, 77 contain an N-terminal toll/interleukin 1 receptor (TIR)-like domain, and 107 of the remaining 361 non-TIR genes contain an N-terminal coiled-coil (CC) domain. Physical map positions were established for 370 predicted NB-LRR genes across all 12 potato chromosomes. The majority of NB-LRRs are physically organised within 63 identified clusters, of which 50 are homogeneous in that they contain NB-LRRs derived from a recent common ancestor. CONCLUSIONS: By establishing the phylogenetic and positional relationship of potato NB-LRRs, our analysis offers significant insight into the evolution of potato R genes. Furthermore, the data provide a blueprint for future efforts to identify and more rapidly clone functional NB-LRR genes from Solanum species.

Raspberry leaf blotch virus, a putative new member of the genus Emaravirus, encodes a novel genomic RNA.

A new, segmented, negative-strand RNA virus with morphological and sequence similarities to other viruses in the genus Emaravirus was discovered in raspberry plants exhibiting symptoms of leaf blotch disorder, a disease previously attributed to the eriophyid raspberry leaf and bud mite (Phyllocoptes gracilis). The virus, tentatively named raspberry leaf blotch virus (RLBV), has five RNAs that each potentially encode a single protein on the complementary strand. RNAs 1, 2 and 3 encode, respectively, a putative RNA-dependent RNA polymerase, a glycoprotein precursor and the nucleocapsid. RNA4 encodes a protein with sequence similarity to proteins of unknown function that are encoded by the genomes of other emaraviruses. When expressed transiently in plants fused to green or red fluorescent protein, the RLBV P4 protein localized to the peripheral cell membrane and to punctate spots in the cell wall. These spots co-localized with GFP-tagged tobacco mosaic virus 30K cell-to-cell movement protein, which is itself known to associate with plasmodesmata. These results suggest that the P4 protein may be a movement protein for RLBV. The fifth RLBV RNA, encoding the P5 protein, is unique among the sequenced emaraviruses. The amino acid sequence of the P5 protein does not suggest any potential function; however, when expressed as a GFP fusion, it localized as small aggregates in the cytoplasm near to the periphery of the cell.

The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants.

FASTQ has emerged as a common file format for sharing sequencing read data combining both the sequence and an associated per base quality score, despite lacking any formal definition to date, and existing in at least three incompatible variants. This article defines the FASTQ format, covering the original Sanger standard, the Solexa/Illumina variants and conversion between them, based on publicly available information such as the MAQ documentation and conventions recently agreed by the Open Bioinformatics Foundation projects Biopython, BioPerl, BioRuby, BioJava and EMBOSS. Being an open access publication, it is hoped that this description, with the example files provided as Supplementary Data, will serve in future as a reference for this important file format.

DNA Metabarcoding and Isolation by Baiting Complement Each Other in Revealing Phytophthora Diversity in Anthropized and Natural Ecosystems.

Isolation techniques supplemented by sequencing of DNA from axenic cultures have provided a robust methodology for the study of Phytophthora communities in agricultural and natural ecosystems. Recently, metabarcoding approaches have emerged as new paradigms for the detection of Phytophthora species in environmental samples. In this study, Illumina DNA metabarcoding and a conventional leaf baiting isolation technique were compared to unravel the variability of Phytophthora communities in different environments. Overall, 39 rhizosphere soil samples from a natural, a semi-natural and a horticultural small-scale ecosystem, respectively, were processed by both baiting and metabarcoding. Using both detection techniques, 28 out of 39 samples tested positive for Phytophthora. Overall, 1,406,613 Phytophthora internal transcribed spacer 1 (ITS1) sequences and 155 Phytophthora isolates were obtained, which grouped into 21 taxa, five retrieved exclusively by baiting (P. bilorbang; P. cryptogea; P. gonapodyides; P. parvispora and P. pseudocryptogea), 12 exclusively by metabarcoding (P. asparagi; P. occultans; P. psycrophila; P. syringae; P. aleatoria/P. cactorum; P. castanetorum/P. quercina; P. iranica-like; P. unknown sp. 1; P. unknown sp. 2; P. unknown sp. 3; P. unknown sp. 4; P. unknown sp. 5) and four with both techniques (P. citrophthora, P. multivora, P. nicotianae and P. plurivora). Both techniques complemented each other in describing the variability of Phytophthora communities from natural and managed ecosystems and revealing the presence of rare or undescribed Phytophthora taxa.

Evolution of relative reading frame bias in unidirectional prokaryotic gene overlaps.

Pairs of unidirectional (same strand) genes can overlap in one of two phases (relative reading frames). There is a striking bias in the relative abundance of prokaryotic gene overlaps in the two possible phases. A simple model is presented based on unidirectional gene overlaps evolving from nonoverlapping gene pairs, through the adoption of alternative start codons by the downstream genes. Potential alternative start codons within upstream gene sequences were found to occur at greater frequencies in one phase, corresponding to the most prevalent phase of gene overlaps. We therefore suggest that the phase bias of overlapping genes is primarily a consequence of the N-terminal extension of downstream genes through adoption of new start codons.