Search

Abstract The Israeli acute paralysis virus (IAPV) is a honeybee-infecting virus that was found to be associated with colony collapse disorder. The IAPV genome contains two genes encoding a structural and a nonstructural polyprotein. We applied a recently developed method for the estimation of selection in overlapping genes to detect purifying selection and, hence, functionality. We provide evolutionary evidence for the existence of a functional overlapping gene, which is translated in the +1 reading frame of the structural polyprotein gene. Conserved orthologs of this putative gene, which we provisionally call pog (predicted overlapping gene), were also found in the genomes of a monophyletic clade of dicistroviruses that includes IAPV, acute bee paralysis virus, Kashmir bee virus, and Solenopsis invicta (red imported fire ant) virus 1.

Abstract Background Same-strand overlapping genes may occur in frameshifts of one (phase 1) or two nucleotides (phase 2). In previous studies of bacterial genomes, long phase-1 overlaps were found to be more numerous than long phase-2 overlaps. This bias was explained by either genomic location or an unspecified selection advantage. Models that focused on the ability of the two genes to evolve independently did not predict this phase bias. Here, we propose that a purely compositional model explains the phase bias in a more parsimonious manner. Same-strand overlapping genes may arise through either a mutation at the termination codon of the upstream gene or a mutation at the initiation codon of the downstream gene. We hypothesized that given these two scenarios, the frequencies of initiation and termination codons in the two phases may determine the number for overlapping genes. Results We examined the frequencies of initiation- and termination-codons in the two phases, and found that termination codons do not significantly differ between the two phases, whereas initiation codons are more abundant in phase 1. We found that the primary factors explaining the phase inequality are the frequencies of amino acids whose codons may combine to form start codons in the two phases. We show that the frequencies of start codons in each of the two phases, and, hence, the potential for the creation of overlapping genes, are determined by a universal amino-acid frequency and species-specific codon usage, leading to a correlation between long phase-1 overlaps and genomic GC content. Conclusion Our model explains the phase bias in same-strand overlapping genes by compositional factors without invoking selection. Therefore, it can be used as a null model of neutral evolution to test selection hypotheses concerning the evolution of overlapping genes. Reviewers This article was reviewed by Bill Martin, Itai Yanai, and Mikhail Gelfand.

Cancers that are poorly immune infiltrated pose a substantial challenge, with current immunotherapies yielding limited clinical success. Stem-like memory T cells (TSCM) have been identified as a subgroup of T cells that possess strong proliferative capacity and that can expand and differentiate following interactions with dendritic cells (DCs). In this study, we explored the pattern of expression of a recently discovered inhibitory receptor poliovirus receptor-related immunoglobulin domain protein (PVRIG) and its ligand, poliovirus receptor-related ligand 2 (PVRL2), in the human tumor microenvironment. Using spatial and single-cell RNA transcriptomics data across diverse cancer indications, we found that among the T-cell checkpoints, PVRIG is uniquely expressed on TSCM and PVRL2 is expressed on DCs in immune aggregate niches in tumors. PVRIG blockade could therefore enhance TSCM-DC interactions and efficiently drive T-cell infiltration to tumors. Consistent with these data, following PVRIG blockade in patients with poorly infiltrated tumors, we observed immune modulation including increased tumor T-cell infiltration, T-cell receptor (TCR) clonality, and intratumoral T-cell expansion, all of which were associated with clinical benefit. These data suggest PVRIG blockade as a promising strategy to induce potent antitumor T-cell responses, providing a novel approach to overcome resistance to immunotherapy in immune-excluded tumors., (©2024 American Association for Cancer Research.)

Proteostasis collapse, the diminished ability to maintain protein homeostasis, has been established as a hallmark of nematode aging. However, whether proteostasis collapse occurs in humans has remained unclear. Here, we demonstrate that proteostasis decline is intrinsic to human senescence. Using transcriptome-wide characterization of gene expression, splicing, and translation, we found a significant deterioration in the transcriptional activation of the heat shock response in stressed senescent cells. Furthermore, phosphorylated HSF1 nuclear localization and distribution were impaired in senescence. Interestingly, alternative splicing regulation was also dampened. Surprisingly, we found a decoupling between different unfolded protein response (UPR) branches in stressed senescent cells. While young cells initiated UPR-related translational and transcriptional regulatory responses, senescent cells showed enhanced translational regulation and endoplasmic reticulum (ER) stress sensing; however, they were unable to trigger UPR-related transcriptional responses. This was accompanied by diminished ATF6 nuclear localization in stressed senescent cells. Finally, we found that proteasome function was impaired following heat stress in senescent cells, and did not recover upon return to normal temperature. Together, our data unraveled a deterioration in the ability to mount dynamic stress transcriptional programs upon human senescence with broad implications on proteostasis control and connected proteostasis decline to human aging., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

The endoplasmic reticulum (ER) stress response, also known as the unfolded protein response (UPR), is a complex cellular response to ER protein misfolding that involves transcriptional regulatory branches and a PERK-mediated translational regulatory branch. Here we revealed that amino acid biosynthesis regulation is coupled to protein synthesis demands during ER stress. Specifically, we demonstrated that the UPR leads to PERK-dependent induction in the biosynthesis of specific amino acids, and to upregulation of their corresponding tRNA synthetases. Furthermore, we found that sequences of UPR-upregulated proteins are significantly enriched with these UPR-induced amino acids. Interestingly, whereas the UPR leads to repression of ER target proteins, we showed that secreted proteins tended to escape this repression and were highly enriched for the UPR-induced amino acids. Our results unravel coordination between amino acid supply, namely, biosynthesis and tRNA loading, and demand from UPR-induced proteins under ER stress, thus revealing an additional regulatory layer of protein synthesis., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Systemic RNAi, initiated by double-stranded RNA (dsRNA) ingestion, has been reported in diverse invertebrates, including honey bees, demonstrating environmental RNA uptake that undermines homologous gene expression. However, the question why any organism would take up RNA from the environment has remained largely unanswered. Here, we report on horizontal RNA flow among honey bees mediated by secretion and ingestion of worker and royal jelly diets. We demonstrate that transmission of jelly-secreted dsRNA to larvae is biologically active and triggers gene knockdown that lasts into adulthood. Worker and royal jellies harbor differential naturally occurring RNA populations. Jelly RNAs corresponded to honey bee protein-coding genes, transposable elements, and non-coding RNA, as well as bacteria, fungi, and viruses. These results reveal an inherent property of honey bees to share RNA among individuals and generations. Our findings suggest a transmissible RNA pathway, playing a role in social immunity and signaling between members of the hive., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

The UPR (Unfolded Protein Response) is a well-orchestrated response to ER protein folding and processing overload, integrating both transcriptional and translational outputs. Its three arms in mammalian cells, the PERK translational response arm, together with the ATF6 and IRE1-XBP1-mediated transcriptional arms, have been thoroughly investigated. Using ribosome footprint profiling, we performed a deep characterization of gene expression programs involved in the early and late ER stress responses, within WT or PERK -/- Mouse Embryonic Fibroblasts (MEFs). We found that both repression and activation gene expression programs, affecting hundreds of genes, are significantly hampered in the absence of PERK. Specifically, PERK -/- cells do not show global translational inhibition, nor do they specifically activate early gene expression programs upon short exposure to ER stress. Furthermore, while PERK -/- cells do activate/repress late ER-stress response genes, the response is substantially weaker. Importantly, we highlight a widespread PERK-dependent repression program, consisting of ER targeted proteins, including transmembrane proteins, glycoproteins, and proteins with disulfide bonds. This phenomenon occurs in various different cell types, and has a major translational regulatory component. Moreover, we revealed a novel interplay between PERK and the XBP1-ATF6 arms of the UPR, whereby PERK attenuates the expression of a specific subset of XBP1-ATF6 targets, further illuminating the complexity of the integrated ER stress response.

Deciphering the global distribution of polyploid plants is fundamental for understanding plant evolution and ecology. Many factors have been hypothesized to affect the uneven distribution of polyploid plants across the globe. Nevertheless, the lack of large comparative datasets has restricted such studies to local floras and to narrow taxonomical scopes, limiting our understanding of the underlying drivers of polyploid plant distribution. We present a map portraying the worldwide polyploid frequencies, based on extensive spatial data coupled with phylogeny-based polyploidy inference for tens of thousands of species. This allowed us to assess the potential global drivers affecting polyploid distribution. Our data reveal a clear latitudinal trend, with polyploid frequency increasing away from the equator. Climate, especially temperature, appears to be the most influential predictor of polyploid distribution. However, we find this effect to be mostly indirect, mediated predominantly by variation in plant lifeforms and, to a lesser extent, by taxonomical composition and species richness. Thus, our study presents an emerging view of polyploid distribution that highlights attributes that facilitate the establishment of new polyploid lineages by providing polyploids with sufficient time (that is, perenniality) and space (low species richness) to compete with pre-adapted diploid relatives.

Background: Recent studies have described a widespread induction of transcriptional readthrough as a consequence of various stress conditions in mammalian cells. This novel phenomenon, initially identified from analysis of RNA-seq data, suggests intriguing new levels of gene expression regulation. However, the mechanism underlying naturally occurring transcriptional readthrough, as well as its regulatory consequences, still remain elusive. Furthermore, the readthrough response to stress has thus far not been investigated outside of mammalian species, and the occurrence of readthrough in many physiological and disease conditions remains to be explored., Results: To facilitate a wider investigation into transcriptional readthrough, we created the DoGFinder software package, for the streamlined identification and quantification of readthrough transcripts, also known as DoGs (Downstream of Gene-containing transcripts), from any RNA-seq dataset. Using DoGFinder, we explore the dependence of DoG discovery potential on RNA-seq library depth, and show that stress-induced readthrough induction discovery is robust to sequencing depth, and input parameter settings. We further demonstrate the use of the DoGFinder software package on a new publically available RNA-seq dataset, and discover DoG induction in human PME cells following hypoxia - a previously unknown readthrough inducing stress type., Conclusions: DoGFinder will enable users to explore, in a few simple steps, the readthrough phenomenon in any condition and organism. DoGFinder is freely available at https://github.com/shalgilab/DoGFinder .

Transcription is a highly regulated process, and stress-induced changes in gene transcription have been shown to play a major role in stress responses and adaptation. Genome-wide studies reveal prevalent transcription beyond known protein-coding gene loci, generating a variety of RNA classes, most of unknown function. One such class, termed downstream of gene-containing transcripts (DoGs), was reported to result from transcriptional readthrough upon osmotic stress in human cells. However, how widespread the readthrough phenomenon is, and what its causes and consequences are, remain elusive. Here we present a genome-wide mapping of transcriptional readthrough, using nuclear RNA-Seq, comparing heat shock, osmotic stress, and oxidative stress in NIH 3T3 mouse fibroblast cells. We observe massive induction of transcriptional readthrough, both in levels and length, under all stress conditions, with significant, yet not complete, overlap of readthrough-induced loci between different conditions. Importantly, our analyses suggest that stress-induced transcriptional readthrough is not a random failure process, but is rather differentially induced across different conditions. We explore potential regulators and find a role for HSF1 in the induction of a subset of heat shock-induced readthrough transcripts. Analysis of public datasets detected increases in polymerase II occupancy in DoG regions after heat shock, supporting our findings. Interestingly, DoGs tend to be produced in the vicinity of neighboring genes, leading to a marked increase in their antisense-generating potential. Finally, we examine genomic features of readthrough transcription and observe a unique chromatin signature typical of DoG-producing regions, suggesting that readthrough transcription is associated with the maintenance of an open chromatin state., Competing Interests: The authors declare no conflict of interest.

Sexual system is a key determinant of genetic variation and reproductive success, affecting evolution within populations and within clades. Much research in plants has focused on evolutionary transitions away from the most common state of hermaphroditism and toward the rare state of dioecy (separate sexes). Rather than transitions predominantly toward greater sexual differentiation, however, evolution may proceed in the direction of lesser sexual differentiation. We analyzed the macroevolutionary dynamics of sexual system in angiosperm genera that contain both dioecious and nondioecious species. Our phylogenetic analyses encompass a total of 2145 species from 40 genera. Overall, we found little evidence that rates of sexual system transitions are greater in any direction. Counting the number of inferred state changes revealed a mild prevalence of transitions away from hermaphroditism and away from dioecy, toward states of intermediate sexual differentiation. We identify genera in which future studies of sexual system evolution might be especially productive, and we discuss how integrating genetic or population-level studies of sexual system could improve the power of phylogenetic comparative analyses. Our work adds to the evidence that different selective pressures and constraints act in different groups, helping maintain the variety of sexual systems observed among plants., (© 2017 The Author(s). Evolution © 2017 The Society for the Study of Evolution.)

[This corrects the article DOI: 10.1002/ece3.2277.].

Premise of the Study: Flowering plants display a variety of sexual systems, ranging from complete cosexuality (hermaphroditism) to separate-sexed individuals (dioecy). While dioecy is relatively rare, it has evolved many times and is present in many plant families. Transitions in sexual systems are hypothesized to be affected by large genomic events such as whole-genome duplication, or polyploidy, and several models have been proposed to explain the observed patterns of association., Methods: In this study, we assessed the association between ploidy and sexual system (separate or combined sexes). To this end, we assembled a database of ploidy levels and sexual systems for ∼1000 species, spanning 18 genera and 15 families. We applied several phylogenetic comparative approaches, including Pagel's coevolutionary framework and sister clade analyses, for detecting correlations between ploidy level and sexual system., Key Results: Our results indicate a broad association between polyploidy and sexual system dimorphism, with low evolutionary stability of the diploid-dioecious condition observed in several clades. A detailed examination of the clades exhibiting this correlation reveals that it is underlain by various patterns of transition rate asymmetry., Conclusions: We conclude that the long-hypothesized connection between ploidy and sexual system holds in some clades, although it may well be affected by factors that differ from clade to clade. Our results further demonstrate that to better understand the evolutionary processes involved, more sophisticated methods and extensive and detailed data sets are required for both broad and focused inquiry., (© 2016 Botanical Society of America.)

Vertebrate sex-determining mechanisms (SDMs) are triggered by the genotype (GSD), by temperature (TSD), or occasionally, by both. The causes and consequences of SDM diversity remain enigmatic. Theory predicts SDM effects on species diversification, and life-span effects on SDM evolutionary turnover. Yet, evidence is conflicting in clades with labile SDMs, such as reptiles. Here, we investigate whether SDM is associated with diversification in turtles and lizards, and whether alterative factors, such as lifespan's effect on transition rates, could explain the relative prevalence of SDMs in turtles and lizards (including and excluding snakes). We assembled a comprehensive dataset of SDM states for squamates and turtles and leveraged large phylogenies for these two groups. We found no evidence that SDMs affect turtle, squamate, or lizard diversification. However, SDM transition rates differ between groups. In lizards TSD-to-GSD surpass GSD-to-TSD transitions, explaining the predominance of GSD lizards in nature. SDM transitions are fewer in turtles and the rates are similar to each other (TSD-to-GSD equals GSD-to-TSD), which, coupled with TSD ancestry, could explain TSD's predominance in turtles. These contrasting patterns can be explained by differences in life history. Namely, our data support the notion that in general, shorter lizard lifespan renders TSD detrimental favoring GSD evolution in squamates, whereas turtle longevity permits TSD retention. Thus, based on the macro-evolutionary evidence we uncovered, we hypothesize that turtles and lizards followed different evolutionary trajectories with respect to SDM, likely mediated by differences in lifespan. Combined, our findings revealed a complex evolutionary interplay between SDMs and life histories that warrants further research that should make use of expanded datasets on unexamined taxa to enable more conclusive analyses.

Dioecy, the sexual system in which male and female organs are found in separate individuals, allows greater specialization for sex-specific functions and can be advantageous under various ecological and environmental conditions. However, dioecy is rare among flowering plants. Previous studies identified contradictory trends regarding the relative diversification rates of dioecious lineages vs their nondioecious counterparts, depending on the methods and data used. We gathered detailed species-level data for dozens of genera that contain both dioecious and nondioecious species. We then applied a probabilistic approach that accounts for differential speciation, extinction, and transition rates between states to examine whether there is an association between dioecy and lineage diversification. We found a bimodal distribution, whereby dioecious lineages exhibited higher diversification in certain genera but lower diversification in others. Additional analyses did not uncover an ecological or life history trait that could explain a context-dependent effect of dioecy on diversification. Furthermore, in-depth simulations of neutral characters demonstrated that such bimodality is also found when simulating neutral characters across the observed trees. Our analyses suggest that - at least for these genera with the currently available data - dioecy neither consistently places a strong brake on diversification nor is a strong driver., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

In general, community similarity is thought to decay with distance; however, this view may be complicated by the relative roles of different ecological processes at different geographical scales, and by the compositional perspective (e.g. species, functional group and phylogenetic lineage) used. Coastal salt marshes are widely distributed worldwide, but no studies have explicitly examined variation in salt marsh plant community composition across geographical scales, and from species, functional and phylogenetic perspectives. Based on studies in other ecosystems, we hypothesized that, in coastal salt marshes, community turnover would be more rapid at local versus larger geographical scales; and that community turnover patterns would diverge among compositional perspectives, with a greater distance decay at the species level than at the functional or phylogenetic levels. We tested these hypotheses in salt marshes of two regions: The southern Atlantic and Gulf Coasts of the United States. We examined the characteristics of plant community composition at each salt marsh site, how community similarity decayed with distance within individual salt marshes versus among sites in each region, and how community similarity differed among regions, using species, functional and phylogenetic perspectives. We found that results from the three compositional perspectives generally showed similar patterns: there was strong variation in community composition within individual salt marsh sites across elevation; in contrast, community similarity decayed with distance four to five orders of magnitude more slowly across sites within each region. Overall, community dissimilarity of salt marshes was lowest on the southern Atlantic Coast, intermediate on the Gulf Coast, and highest between the two regions. Our results indicated that local gradients are relatively more important than regional processes in structuring coastal salt marsh communities. Our results also suggested that in ecosystems with low species diversity, functional and phylogenetic approaches may not provide additional insight over a species-based approach.

Prokaryotic genomes are small and compact. Either this feature is caused by neutral evolution or by natural selection favoring small genomes-genome streamlining. Three separate prior lines of evidence argue against streamlining for most prokaryotes. We find that the same three lines of evidence argue for streamlining in the genomes of thermophile bacteria. Specifically, with increasing habitat temperature and decreasing genome size, the proportion of genomic DNA in intergenic regions decreases. Furthermore, with increasing habitat temperature, generation time decreases. Genome-wide selective constraints do not decrease as in the reduced genomes of host-associated species. Reduced habitat variability is not a likely explanation for the smaller genomes of thermophiles. Genome size may be an indirect target of selection due to its association with cell volume. We use metabolic modeling to demonstrate that known changes in cell structure and physiology at high temperature can provide a selective advantage to reduce cell volume at high temperatures.

New protein-coding genes can originate either through modification of existing genes or de novo. Recently, the importance of de novo origination has been recognized in eukaryotes, although eukaryotic genes originated de novo are relatively rare and difficult to identify. In contrast, viruses contain many de novo genes, namely those in which an existing gene has been "overprinted" by a new open reading frame, a process that generates a new protein-coding gene overlapping the ancestral gene. We analyzed the evolution of 12 experimentally validated viral genes that originated de novo and estimated their relative ages. We found that young de novo genes have a different codon usage from the rest of the genome. They evolve rapidly and are under positive or weak purifying selection. Thus, young de novo genes might have strain-specific functions, or no function, and would be difficult to detect using current genome annotation methods that rely on the sequence signature of purifying selection. In contrast to young de novo genes, older de novo genes have a codon usage that is similar to the rest of the genome. They evolve slowly and are under stronger purifying selection. Some of the oldest de novo genes evolve under stronger selection pressure than the ancestral gene they overlap, suggesting an evolutionary tug of war between the ancestral and the de novo gene.

An overlapping open reading frame (ORF) with a potential to encode a functional protein has been identified within the 3'-proximal ORF of Solenopsis invicta virus 1 (SINV-1) and three bee viruses. This ORF has been referred to as predicted overlapping gene (pog). Protein motif searches of POG revealed weak relationships precluding assignment of a potential function. Neither a transcript nor a protein encoded by the pog ORF has been detected. However, recently, a protein encoded by the corresponding +1 overlapping ORF (termed ORFx) in the Israeli acute paralysis virus (IAPV) was demonstrated by recombinant means as well as in IAPV-infected honey bees. The objective of our study was to attempt to provide empirical evidence for the presence of a pog-derived protein from SINV-1-infected fire ants. A number of different laboratory and field SINV-1-infected Solenopsis invicta preparations were examined by western blotting for the presence of a POG protein sequence. In every case, these preparations failed to yield any detectable bands when probed with a polyclonal antibody preparation raised to a portion of the pog predicted protein sequence. Although impossible to prove a negative result, proper controls used in these studies suggested that the pog ORF is not translated into a functional protein in SINV-1.

Protein-coding genes often contain long overlapping open-reading frames (ORFs), which may or may not be functional. Current methods that utilize the signature of purifying selection to detect functional overlapping genes are limited to the analysis of sequences from divergent species, thus rendering them inapplicable to genes found only in closely related sequences. Here, we present a method for the detection of selection signatures on overlapping reading frames by using closely related sequences, and apply the method to several known overlapping genes, and to an overlapping ORF on the negative strand of segment 8 of influenza A virus (NEG8), for which the suggestion has been made that it is functional. We find no evidence that NEG8 is under selection, suggesting that the intact reading frame might be non-functional, although we cannot fully exclude the possibility that the method is not sensitive enough to detect the signature of selection acting on this gene. We present the limitations of the method using known overlapping genes and suggest several approaches to improve it in future studies. Finally, we examine alternative explanations for the sequence conservation of NEG8 in the absence of selection. We show that overlap type and genomic context affect the conservation of intact overlapping ORFs and should therefore be considered in any attempt of estimating the signature of selection in overlapping genes.

As far as protein-coding genes are concerned, there is a non-zero probability that at least one of the five possible overlapping sequences of any gene will contain an open-reading frame (ORF) of a length that may be suitable for coding a functional protein. It is, however, very difficult to determine whether or not such an ORF is functional. Recently, we proposed a method that predicts functionality of an overlapping ORF if it can be shown that it has been subject to purifying selection during its evolution. Here, we use simulation to test this method under several conditions and compare it with the method of Firth and Brown. We found that under most conditions, our method detects functional overlapping genes with higher sensitivity than Firth and Brown's method, while maintaining high specificity. Further, we tested the hypothesis that the two aminoacyl tRNA synthetase classes have originated from a pair of overlapping genes. A central piece of evidence ostensibly supporting this hypothesis is the assertion that an overlapping ORF of a heat-shock protein-70 gene, which exhibits some similarity to class 2 aminoacyl tRNA synthetases, is functional. We found signature of purifying selection only in highly divergent sequences, suggesting that the method yields false-positives in high sequence divergence and that the overlapping ORF is not a functional gene. Finally, we examined three cases of overlap in the human genome. We find varying signatures of purifying selection acting on these overlaps, raising the possibility that two of the overlapping genes may not be functional.

Published estimates of the proportion of positively selected genes (PSGs) in human vary over three orders of magnitude. In mammals, estimates of the proportion of PSGs cover an even wider range of values. We used 2,980 orthologous protein-coding genes from human, chimpanzee, macaque, dog, cow, rat, and mouse as well as an established phylogenetic topology to infer the fraction of PSGs in all seven terminal branches. The inferred fraction of PSGs ranged from 0.9% in human through 17.5% in macaque to 23.3% in dog. We found three factors that influence the fraction of genes that exhibit telltale signs of positive selection: the quality of the sequence, the degree of misannotation, and ambiguities in the multiple sequence alignment. The inferred fraction of PSGs in sequences that are deficient in all three criteria of coverage, annotation, and alignment is 7.2 times higher than that in genes with high trace sequencing coverage, "known" annotation status, and perfect alignment scores. We conclude that some estimates on the prevalence of positive Darwinian selection in the literature may be inflated and should be treated with caution.

Inferring the intensity of positive selection in protein-coding genes is important since it is used to shed light on the process of adaptation. Recently, it has been reported that overlapping genes, which are ubiquitous in all domains of life, seem to exhibit inordinate degrees of positive selection. Here, we present a new method for the simultaneous estimation of selection intensities in overlapping genes. We show that the appearance of positive selection is caused by assuming that selection operates independently on each gene in an overlapping pair, thereby ignoring the unique evolutionary constraints on overlapping coding regions. Our method uses an exact evolutionary model, thereby voiding the need for approximation or intensive computation. We test the method by simulating the evolution of overlapping genes of different types as well as under diverse evolutionary scenarios. Our results indicate that the independent estimation approach leads to the false appearance of positive selection even though the gene is in reality subject to negative selection. Finally, we use our method to estimate selection in two influenza A genes for which positive selection was previously inferred. We find no evidence for positive selection in both cases.