Overlapping Reading Frames Research Articles

Genomic Analysis of Hepatitis B Virus Reveals Antigen State and Genotype as Sources of Evolutionary Rate Variation

Hepatitis B virus (HBV) genomes are small, semi-double-stranded DNA circular genomes that contain alternating overlapping reading frames and replicate through an RNA intermediary phase. This complex biology has presented a challenge to estimating an evolutionary rate for HBV, leading to difficulties resolving the evolutionary and epidemiological history of the virus. Here, we re-examine rates of HBV evolution using a novel data set of 112 within-host, transmission history (pedigree) and among-host genomes isolated over 20 years from the indigenous peoples of the South Pacific, combined with 313 previously published HBV genomes. We employ Bayesian phylogenetic approaches to examine several potential causes and consequences of evolutionary rate variation in HBV. Our results reveal rate variation both between genotypes and across the genome, as well as strikingly slower rates when genomes are sampled in the Hepatitis B e antigen positive state, compared to the e antigen negative state. This Hepatitis B e antigen rate variation was found to be largely attributable to changes during the course of infection in the preCore and Core genes and their regulatory elements.

Molecular biology of the hepatitis B virus and role of the X gene

The hepatitis B virus (HBV) is a widespread human pathogen and a major health problem in many countries. Molecular cloning and sequencing of the viral DNA genome has demonstrated a small and compact structure organized into four overlapping reading frames that encode the viral proteins. Besides structural proteins of the core and the envelope, HBV encodes a DNA polymerase with reverse transcriptase activity, a secreted antigen of unknown function, and a transcriptional activator that is essential for viral replication. Major steps of the viral life cycle have been unraveled, including transcription of all viral RNAs from nuclear covalently closed circular DNA (cccDNA), followed by encapsidation of pregenomic RNA, a more-than-genome length transcript, and reverse transcription of pregenomic RNA leading to asymmetric synthesis of the DNA strands. Although HBV has been recognized as a human tumor virus, no direct transforming activity could be evidenced in different cellular and animal models. However, the transcriptional regulatory protein HBx encoded by the X gene is endowed with weak oncogenic activity. HBx harbors pleiotropic activities and plays a major role in HBV pathogenesis and in liver carcinogenesis.

CD8 T cell response and evolutionary pressure to HIV-1 cryptic epitopes derived from antisense transcription

Retroviruses pack multiple genes into relatively small genomes by encoding several genes in the same genomic region with overlapping reading frames. Both sense and antisense HIV-1 transcripts contain open reading frames for known functional proteins as well as numerous alternative reading frames (ARFs). At least some ARFs have the potential to encode proteins of unknown function, and their antigenic properties can be considered as cryptic epitopes (CEs). To examine the extent of active immune response to virally encoded CEs, we analyzed human leukocyte antigen class I–associated polymorphisms in HIV-1 gag, pol, and nef genes from a large cohort of South Africans with chronic infection. In all, 391 CEs and 168 conventional epitopes were predicted, with the majority (307; 79%) of CEs derived from antisense transcripts. In further evaluation of CD8 T cell responses to a subset of the predicted CEs in patients with primary or chronic infection, both sense- and antisense-encoded CEs were immunogenic at both stages of infection. In addition, CEs often mutated during the first year of infection, which was consistent with immune selection for escape variants. These findings indicate that the HIV-1 genome might encode and deploy a large potential repertoire of unconventional epitopes to enhance vaccine-induced antiviral immunity.

Virology and clinical sequelae of drug-resistant HBV in HIV-HBV-coinfected patients on highly active antiretroviral therapy

Several of the nucleoside/nucleotide analogues used to treat HIV also inhibit HBV replication, with lamivudine being the oldest of this group. Thus, prior to licensing of tenofovir, many HIV-HBV-coinfected individuals received lamivudine as the only drug active against HBV as part of an anti-HIV regimen, which set the stage for the emergence of drug-resistant HBV. In coinfected persons, lamivudine-resistant HBV develops more rapidly than in HBV-monoinfected persons, but it is not known if this is true for the newer agents. Owing to overlapping reading frames of the HBV polymerase and surface antigens, drug-resistant changes in HBV Pol can lead to mutations in the envelope. This review will discuss studies of drug-resistant HBV in HIV-infected persons including drug-resistant mutations that have been identified and clinical sequelae of these mutations.

Selection characterization on overlapping reading frame of multiple-protein-encoding P gene in Newcastle disease virus

The aim of this study was to characterize the molecular evolution of P and V protein genes of the Newcastle disease virus (NDV). The P gene sequences of 55 NDV isolates, representing different chronological and geographic origins, were obtained from GenBank. In this paper, the evolution of the specific regions of the NDV P gene, encoding the P and V proteins, was analyzed. The nucleotides from the shared P/V region encoded the co-amino terminus of the two proteins, while the P–V/V–P region was respectively encoded by the nucleotides within the P ORF or the V ORF in the common sequence (after the mRNA editing site). As well, the P-cut region exclusively encoded the P protein. Finally, the P–V and V–P regions were further broken down into P1 and P2 fragments with the corresponding V1 and V2 fragments. In the P gene, the P-cut portion corresponding to the C-terminal of the P protein was the most highly conserved, while the P–V region was the most variable. This was interpreted as a lower constraint for function in the common sequence than in the unique P sequence that is known to contain an important function. Interestingly, in the common P–V/V–P function, variability of V1 was compensated by a higher conservation of the corresponding P1, and conversely for the P2/V2, which suggested that the flexibility of one ORF with less function served the purpose of allowing positive selection in the other overlapping ORF that exhibited more function.

Tula hantavirus NSs protein accumulates in the perinuclear area in infected and transfected cells

The small RNA segment of some hantaviruses (family Bunyaviridae) encodes two proteins: the nucleocapsid protein and, in an overlapping reading frame, a non-structural (NSs) protein. The hantavirus NSs protein, like those of orthobunya- and phleboviruses, counteracts host innate immunity. Here, for the first time, the NSs protein of a hantavirus (Tula virus) has been observed in infected cells and shown to localize in the perinuclear area. Transiently expressed NSs protein showed similar localization, although the kinetics was slightly different, suggesting that to reach its proper location in the infected cell, the NSs protein does not have to cooperate with other viral proteins.

A highly attenuated measles virus vaccine strain encodes a fully functional C protein.

The P, V, and C proteins of measles virus are encoded in overlapping reading frames of the P gene, which makes it difficult to analyze the functions of the individual proteins in the context of virus infection. We established a system to analyze the C protein independently from the P and V proteins by placing its gene in an additional transcription unit between the H and L genes. Analyses with this system indicated that a highly attenuated Edmonston lineage vaccine strain encodes a fully functional C protein, and the P and/or V protein is involved in the attenuated phenotype.

Stem cells, ribozymes and HIV

In a recent Phase 2 clinical trial, Mitsuyasu et al. demonstrated that the delivery of an anti-HIV ribozyme from a gammaretroviral vector integrated into hematopoietic stem cells (HSCs) offers a safe and promising gene therapy strategy for the treatment of this infectious disease. This gene therapy strategy relies on the idea that, after the ex vivo transduction of an anti-HIV gene into autologous CD34+ HSCs, these human immunodeficiency virus (HIV)-resistant progenitors can be re-infused into an HIV patient, where they will differentiate and expand to permanently protect against the onset of AIDS (Figure 1). Although the inhibition of HIV has been considered a feasible and attractive candidate for gene therapy applications for over 20 years, this is the first reported randomized, double-blind cell-delivered Phase 2 clinical trial for HIV-1. The Mitsuyasu et al. study monitored 74 HIV-1-infected individuals for 100 weeks after infusion of transduced or control autologous HSCs. Thirty-eight of these patients were treated with HSCs carrying an integrated copy of a therapeutic, replication-incompetent, Moloney murine leukemia virus-based gammaretroviral vector (LNL6). The therapeutic vector (OZ1) expresses a hammerhead ribozyme against overlapping reading frames of the viral genes vpr and tat in unspliced and spliced viral transcripts, respectively. Although the researchers observed no significant differences in the viral load between OZ1 and control groups at any particular time point, they reported that transduced CD34+ cell therapy did not cause any apparent adverse effects and is safe for further exploration, thereby opening the doors for further trials with other gene therapy agents. Over the course of the 100-week clinical study, Mitsuyasu et al. observed sharp declines in the percentages of patients who maintained detectable levels of integrated DNA and mRNA from the OZ1 retroviral vector. In fact, at 20 weeks post-infusion, primary blood mononuclear cells from more than half of the treated patients had no detectable levels of OZ1 DNA or mRNA. This result is unexpected under the assumption that the OZ1-transduced HSCs successfully engrafted, differentiated and expanded into the T-cell lineage, and maintained some level of anti-HIV protection over untreated CD4+ T cells. Similar Phase I trials—including reports from the same authors—have confirmed successful engraftment, differentiation and expansion of autologous HSCs transduced with anti-HIV retroviral vectors. Therefore, although the HSCs may have failed to engraft under the conditions of the trial, it is also possible that the anti-HIV ribozyme failed to protect the OZ1 cells from HIV-1 infection and replication. This could be the result of mutational escape of the virus, suboptimal function of the ribozyme in vivo or possible shortcomings associated with gammaretroviral vectors. This Phase 2 investigation marks an important step toward a safe, cell-delivered gene therapy for HIV. However, to further advance this promising gene therapy approach, it will be necessary to identify the particular limitations and modify the existing therapeutic strategy. First, the virus might have acquired mutations in the target region of the ribozyme that would abrogate the function of the anti-HIV agent. Mutations arising at the second or third positions of the three-nucleotide ribozyme cleavage site of spliced and unspliced viral transcripts might render the anti-HIV ribozyme ineffective. Close examination of the targeted region suggests

Overlapping Genes Produce Proteins with Unusual Sequence Properties and Offer Insight into De Novo Protein Creation

It is widely assumed that new proteins are created by duplication, fusion, or fission of existing coding sequences. Another mechanism of protein birth is provided by overlapping genes. They are created de novo by mutations within a coding sequence that lead to the expression of a novel protein in another reading frame, a process called "overprinting." To investigate this mechanism, we have analyzed the sequences of the protein products of manually curated overlapping genes from 43 genera of unspliced RNA viruses infecting eukaryotes. Overlapping proteins have a sequence composition globally biased toward disorder-promoting amino acids and are predicted to contain significantly more structural disorder than nonoverlapping proteins. By analyzing the phylogenetic distribution of overlapping proteins, we were able to confirm that 17 of these had been created de novo and to study them individually. Most proteins created de novo are orphans (i.e., restricted to one species or genus). Almost all are accessory proteins that play a role in viral pathogenicity or spread, rather than proteins central to viral replication or structure. Most proteins created de novo are predicted to be fully disordered and have a highly unusual sequence composition. This suggests that some viral overlapping reading frames encoding hypothetical proteins with highly biased composition, often discarded as noncoding, might in fact encode proteins. Some proteins created de novo are predicted to be ordered, however, and whenever a three-dimensional structure of such a protein has been solved, it corresponds to a fold previously unobserved, suggesting that the study of these proteins could enhance our knowledge of protein space.

Hepatitis B Viral Surface Mutations in Patients with Adefovir Resistant Chronic Hepatitis B with A181T/V Polymerase Mutations

The hepatitis B virus (HBV) polymerase gene has overlapping reading frames with surface genes, which allows to alter the amino acid codon of the surface genes. In adefovir (ADV) treated chronic hepatitis B patients carrying rtA181T/rtA181V mutations, overlap with surface gene mutations such as sW172stop/sL173F has been reported. However, the clinical consequences of such surface mutations have not been determined. The aim of this study was to determine the surface gene sequence in ADV-resistant patients carrying the A181T/V mutation and to describe the clinical significance. Of the 22 patients included in this study, 13 were ADV-resistant with rtA181T/V mutations (polymerase mutation group, Group P) and nine were antiviral treatment-naïve (control group, Group C). The Pre-S1 gene mutation, V60A, was detected in 11 patients (Group P=8, Group C=3). A start codon mutation in the Pre-S2 gene was found in five patients (Group P=3, Group C=2). An S gene mutation, sA184V, was found in nine patients, all of whom were in group P. Although sW172stop and sL173F mutations were detected, reduced HBsAg titer was not observed. Further study of these mutations and their clinical implications are needed.

Investigating selection on viruses: a statistical alignment approach

BackgroundTwo problems complicate the study of selection in viral genomes: Firstly, the presence of genes in overlapping reading frames implies that selection in one reading frame can bias our estimates of neutral mutation rates in another reading frame. Secondly, the high mutation rates we are likely to encounter complicate the inference of a reliable alignment of genomes. To address these issues, we develop a model that explicitly models selection in overlapping reading frames. We then integrate this model into a statistical alignment framework, enabling us to estimate selection while explicitly dealing with the uncertainty of individual alignments. We show that in this way we obtain un-biased selection parameters for different genomic regions of interest, and can improve in accuracy compared to using a fixed alignment.ResultsWe run a series of simulation studies to gauge how well we do in selection estimation, especially in comparison to the use of a fixed alignment. We show that the standard practice of using a ClustalW alignment can lead to considerable biases and that estimation accuracy increases substantially when explicitly integrating over the uncertainty in inferred alignments. We even manage to compete favourably for general evolutionary distances with an alignment produced by GenAl. We subsequently run our method on HIV2 and Hepatitis B sequences.ConclusionWe propose that marginalizing over all alignments, as opposed to using a fixed one, should be considered in any parametric inference from divergent sequence data for which the alignments are not known with certainty. Moreover, we discover in HIV2 that double coding regions appear to be under less stringent selection than single coding ones. Additionally, there appears to be evidence for differential selection, where one overlapping reading frame is under positive and the other under negative selection.

Frequent p16-Independent Inactivation of p14ARF in Human Melanoma

The tumor suppressors p14(ARF) (ARF) and p16(INK4A) (p16) are encoded by overlapping reading frames at the CDKN2A/INK4A locus on chromosome 9p21. In human melanoma, the accumulated evidence has suggested that the predominant tumor suppressor at 9p21 is p16, not ARF. However, recent observations from melanoma-prone families and murine melanoma models suggest a p16-independent tumor suppressor role for ARF. We analyzed a group of melanoma metastases and cell lines to investigate directly whether somatic alterations to the ARF gene support its role as a p16-independent tumor suppressor in human melanoma, assuming that two alterations (genetic and/or epigenetic) would be required to inactivate a gene. We examined the p16/ARF locus in 60 melanoma metastases from 58 patients and in 9 human melanoma cell lines using multiplex ligation-dependent probe amplification and multiplex polymerase chain reaction (PCR) to detect deletions, methylation-specific PCR to detect promoter methylation, direct sequencing to detect mutations affecting ARF and p16, and, in a subset of 20 tumors, immunohistochemistry to determine the effect of these alterations on p16 protein expression. All statistical tests were two-sided. We observed two or more alterations to the ARF gene in 26/60 (43%) metastases. The p16 gene sustained two or more alterations in 13/60 (22%) metastases (P = .03). Inactivation of ARF in the presence of wild-type p16 was seen in 18/60 (30%) metastases. Genetic and epigenetic analyses of the human 9p21 locus indicate that modifications of ARF occur independently of p16 inactivation in human melanoma and suggest that ARF is more frequently inactivated than p16.

Continuous Crossover(s) Events of HIV-1 CRF01_AE and B Subtype Strains in Malaysia: Evidence of Rapid and Extensive HIV-1 Evolution in the Region

The Asian HIV epidemic appears to be complex, characterized by the prevalence of multiple subtypes and circulating recombinant forms with gradual replacement of pure HIV-1 subtypes in several geographical regions. The main objectives of the present study are to identify and analyse the full-length viral genomes of three unique recombinant forms (URFs); the HIV-1 isolates 07MYKLD47, 07MYKLD48 and 07MYKLD49 from Malaysia. Long-range polymerase chain reaction (PCR) amplification of seven overlapping reading frames was used to derive near full-length HIV-1 genomes. Detailed phylogenetic and bootscanning analyses were performed to determine phylogenetic associations and subtypic assignments. We further confirmed the mosaic composition of these CRF01_AE/B inter-subtype recombinant forms, which are composed of B-subtype fragment(s) in the backbone of CRF01_AE. Both 07MYKLD47 and 07MYKLD48 have an insertion of B subtype (880 bp and 532 bp) in the gag-pol and gp41-env gene regions, respectively. Whereas the isolate 07MYKLD49 has three B-subtype fragments inserted in different gene region along the genome; one each in the gag-pol (1862 bp) and pol-vif (1935 bp) regions, and a short B-subtype insertion (541 bp) in the 5' LTR-gag region. This highlights the public health relevance of newly emerging second generation HIV-1 recombinant forms and their dispersal, along with their rapid and continuous evolution in the region.

Large gene overlaps in prokaryotic genomes: result of functional constraints or mispredictions?

BackgroundAcross the fully sequenced microbial genomes there are thousands of examples of overlapping genes. Many of these are only a few nucleotides long and are thought to function by permitting the coordinated regulation of gene expression. However, there should also be selective pressure against long overlaps, as the existence of overlapping reading frames increases the risk of deleterious mutations. Here we examine the longest overlaps and assess whether they are the product of special functional constraints or of erroneous annotation.ResultsWe analysed the genes that overlap by 60 bps or more among 338 fully-sequenced prokaryotic genomes. The likely functional significance of an overlap was determined by comparing each of the genes to its respective orthologs. If a gene showed a significantly different length from its orthologs it was considered unlikely to be functional and therefore the result of an error either in sequencing or gene prediction. Focusing on 715 co-directional overlaps longer than 60 bps, we classified the erroneous ones into five categories: i) 5'-end extension of the downstream gene due to either a mispredicted start codon or a frameshift at 5'-end of the gene (409 overlaps), ii) fragmentation of a gene caused by a frameshift (163), iii) 3'-end extension of the upstream gene due to either a frameshift at 3'-end of a gene or point mutation at the stop codon (68), iv) Redundant gene predictions (4), v) 5' & 3'-end extension which is a combination of i) and iii) (71). We also studied 75 divergent overlaps that could be classified as misannotations of group i). Nevertheless we found some convergent long overlaps (54) that might be true overlaps, although an important part of convergent overlaps could be classified as group iii) (124).ConclusionAmong the 968 overlaps larger than 60 bps which we analysed, we did not find a single real one among the co-directional and divergent orientations and concluded that there had been an excessive number of misannotations. Only convergent orientation seems to permit some long overlaps, although convergent overlaps are also hampered by misannotations. We propose a simple rule to flag these erroneous gene length predictions to facilitate automatic annotation.

Detection and characterization of a functional insertion sequence, ISVpa 2, in Vibrio parahaemolyticus

PCR analysis of the pandemic strain of Vibrio parahaemolyticus, KX-V237 (total genome sequenced) showed a subculture where the size of the amplicons had increased. The purpose of this study was to analyze the mechanism of this change. We found a 1,243-bp DNA sequence inserted in one of the pandemic marker genes in this strain. The inserted DNA sequence possessed the genetic structures shared by insertion sequences (ISs) of the IS 3 family. This IS had 26-bp imperfect terminal inverted repeats (IRs) and two partially overlapping reading frames, orfA and orfB. OrfA codes for a helix-turn-helix, OrfA and OrfAB produced by translational frameshifting code for leucine zipper motifs, and OrfB codes for a DDE motif. orfA and orfB were homologous to those in the IS 3 family. This IS was named ISVpa 2. Southern blot analysis showed the copy number of ISVpa 2 in our stock culture and its subculture of KX-V237 was three and four, respectively; whereas it was only one in the reported genome sequence. Analysis of the flanking sequences for seven ISVpa 2 copies showed ISVpa 2 is capable of inserting at multiple sites and ISVpa 2 causes genetic rearrangements including insertional inactivation of the target gene and adjacent deletion. ISVpa 2 created 3-base duplications upon insertion. PCR, hybridization, and nucleotide sequence analyses showed ISVpa 2 homologs were detected in all of the 62 other strains of V. parahaemolyticus examined; and in some strains of Vibrio vulnificus (98% identity), Vibrio penaeicida (86% identity), and Vibrio splendidus (87% identity); but was not in 25 other species in the genus Vibrio. The data demonstrate that ISVpa 2 is a transpositionally active IS discovered for the first time in V. parahaemolyticus and suggest that ISVpa 2 may be transferred among the species of the genus Vibrio.

Mammalian Gene PEG10 Expresses Two Reading Frames by High Efficiency –1 Frameshifting in Embryonic-associated Tissues

Paternally expressed gene 10 (PEG10) is a mammalian gene that is essential for embryonic development in mice. The gene contains two overlapping open reading frames (ORF1 and ORF2) and is derived from a retroelement that acquired a cellular function. It is not known if both reading frames are required for PEG10 function. Synthesis of ORF2 would be possible only if programmed -1 frameshifting occurred during ORF1 translation. In this study the frameshifting activity of PEG10 was analyzed in vivo, and a potential role for ORF2 was investigated. Phylogenetic analysis demonstrated that PEG10 is highly conserved in therian mammals, with all species retaining the elements necessary for frameshifting as well as functional motifs in each ORF. The frameshift site of PEG10 was highly active in cultured cells and produced the ORF1-2 protein. In mice, endogenous ORF1 and an ORF1-2 frameshift protein were detected in the developing placenta and amniotic membrane from 9.5 days post-coitus through to term with a very high frameshift efficiency (>60%). Mutagenesis of the active site motif of a putative protease within ORF2 showed that this enzyme is active and participates in post-translational processing of PEG10 ORF1-2. Both PEG10 proteins were also detected in first trimester human placenta. By contrast, neither protein expression nor frameshifting was detected in adult mouse tissues. These studies imply that the ORF1-2 protein, synthesized utilizing the most efficient -1 frameshift mechanism yet documented in vivo, will have an essential function that is intrinsic to the importance of PEG10 in mammals.

A Nested Gene in Streptomyces Bacteria Encodes a Protein Involved in Quaternary Complex Formation

The gene pgaM is involved in the biosynthesis of an angucycline-type polyketide antibiotic in Streptomyces sp. PGA64. It encodes a two-domain polypeptide consisting of an N-terminal flavoprotein oxygenase and a C-terminal short-chain alcohol dehydrogenase/reductase, which are fused together at the translational level as a result of end codon deletion. Here we show that translation also initiates at an internal start codon that enables independent expression of a separate reductase subunit, PgaMred. This confirms that the gene exhibits a rare viral-like arrangement of two overlapping reading frames that allows simultaneous expression of two alternative forms of the protein. Together, these two proteins associate to form a stable non-covalent complex, the native form of PgaM. The reductase subunit PgaMred is shown to provide enzyme stability and to affect the redox state of the oxygenase domain FAD. Finally, a model for the quaternary structure of the complex that explains the necessity for a nested gene system and the unusual behaviour of the protein subunits in vitro is presented.

A Bifunctional Locus (BIO3-BIO1) Required for Biotin Biosynthesis in Arabidopsis

We identify here the Arabidopsis (Arabidopsis thaliana) gene encoding the third enzyme in the biotin biosynthetic pathway, dethiobiotin synthetase (BIO3; At5g57600). This gene is positioned immediately upstream of BIO1, which is known to be associated with the second reaction in the pathway. Reverse genetic analysis demonstrates that bio3 insertion mutants have a similar phenotype to the bio1 and bio2 auxotrophs identified using forward genetic screens for arrested embryos rescued on enriched nutrient medium. Unexpectedly, bio3 and bio1 mutants define a single genetic complementation group. Reverse transcription-polymerase chain reaction analysis demonstrates that separate BIO3 and BIO1 transcripts and two different types of chimeric BIO3-BIO1 transcripts are produced. Consistent with genetic data, one of the fused transcripts is monocistronic and encodes a bifunctional fusion protein. A splice variant is bicistronic, with distinct but overlapping reading frames. The dual functionality of the monocistronic transcript was confirmed by complementing the orthologous auxotrophs of Escherichia coli (bioD and bioA). BIO3-BIO1 transcripts from other plants provide further evidence for differential splicing, existence of a fusion protein, and localization of both enzymatic reactions to mitochondria. In contrast to most biosynthetic enzymes in eukaryotes, which are encoded by genes dispersed throughout the genome, biotin biosynthesis in Arabidopsis provides an intriguing example of a bifunctional locus that catalyzes two sequential reactions in the same metabolic pathway. This complex locus exhibits several unusual features that distinguish it from biotin operons in bacteria and from other genes known to encode bifunctional enzymes in plants.

Annotation of selection strengths in viral genomes

Viral genomes tend to code in overlapping reading frames to maximize informational content. This may result in atypical codon bias and particular evolutionary constraints. Due to the fast mutation rate of viruses, there is additional strong evidence for varying selection between intra- and intergenomic regions. The presence of multiple coding regions complicates the concept of K(a)/K(s) ratio, and thus begs for an alternative approach when investigating selection strengths. Building on the paper by McCauley and Hein, we develop a method for annotating a viral genome coding in overlapping reading frames. We introduce an evolutionary model capable of accounting for varying levels of selection along the genome, and incorporate it into our prior single sequence HMM methodology, extending it now to a phylogenetic HMM. Given an alignment of several homologous viruses to a reference sequence, we may thus achieve an annotation both of coding regions as well as selection strengths, allowing us to investigate different selection patterns and hypotheses. We illustrate our method by applying it to a multiple alignment of four HIV2 sequences, as well as of three Hepatitis B sequences. We obtain an annotation of the coding regions, as well as a posterior probability for each site of the strength of selection acting on it. From this we may deduce the average posterior selection acting on the different genes. Whilst we are encouraged to see in HIV2, that the known to be conserved genes gag and pol are indeed annotated as such, we also discover several sites of less stringent negative selection within the env gene. To the best of our knowledge, we are the first to subsequently provide a full selection annotation of the Hepatitis B genome by explicitly modelling the evolution within overlapping reading frames, and not relying on simple K(a)/K(s) ratios.

Eliminating the Requirement of an Essential Gene Product in an Already Very Small Virus: Scaffolding Protein B-free øX174, B-free

Unlike most viral assembly systems, two scaffolding proteins, B and D, mediate bacteriophage øX174 morphogenesis. The external scaffolding protein D is highly ordered in the atomic structure and proper function is very sensitive to mutation. In contrast, the internal scaffolding protein B is relatively unordered and extensive alterations do not eliminate function. Despite this genetic laxity, protein B is absolutely required for virus assembly. Thus, this system, with its complex arrangements of overlapping reading frames, can be regarded as an example of “irreducible complexity.” To address the biochemical functions of a dual scaffolding protein system and the evolution of complexity, progressive and targeted genetic selections were employed to lessen and finally eliminate B protein-dependence. The biochemical and genetic bases of adaptation were characterized throughout the analysis that led to the sextuple mutant with a B-independent phenotype, as evaluated by plaque formation in wild-type cells. The primary adaptation appears to be the over-expression of a mutant external scaffolding protein. Progeny production was followed in lysis-resistant cells. The ability to produce infectious virions does not require all six mutations. However, the lag phase before progeny production is shortened as mutations accumulate. The results suggest that the primary function of the internal scaffolding protein may be to lower the critical concentration of the external scaffolding protein needed to nucleate procapsid formation. Moreover, they demonstrate a novel mechanism by which a stringently required gene product can be bypassed, even in a system encoding only eight strictly essential proteins.