Immune-induced evolutionary selection focused on a single reading frame in overlapping hepatitis B virus proteins

Viruses employ various means to evade immune detection. Reduction of CD8(+) T cell epitopes is one of the common strategies used for this purpose. Hepatitis B virus (HBV), a member of the Hepadnaviridae family, has four open reading frames, with about 50% overlap between the genes they encode. We computed the CD8(+) T cell epitope density within HBV proteins and the mutations within the epitopes. Our results suggest that HBV accumulates escape mutations that reduce the number of epitopes. These mutations are not equally distributed among genes and reading frames. While the highly expressed core and X proteins are selected to have low epitope density, polymerase, which is expressed at low levels, does not undergo the same selection. In overlapping regions, mutations in one protein-coding sequence also affect the other protein-coding sequence. We show that mutations lead to the removal of epitopes in X and surface proteins even at the expense of the addition of epitopes in polymerase. The total escape mutation rate for overlapping regions is lower than that for nonoverlapping regions. The lower epitope replacement rate for overlapping regions slows the evolutionary escape rate of these regions but leads to the accumulation of mutations more robust in the transfer between hosts, such as mutations preventing proteasomal cleavage into epitopes.     

Annotation of selection strengths in viral genomes

MOTIVATION: 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. RESULTS: 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.     

Overlapping of MINK and CHRNE gene loci in the course of mammalian evolution

Overlapping of genes, especially in an anti-parallel fashion, is quite rare in eukaryotic genomes. We have found a rare instance of exon overlapping involving CHRNE and MINK gene loci on chromosome 17 in humans. CHRNE codes for the subunit of the nicotinic acetylcholine receptor (AChR) whereas MINK encodes a serine/threonine kinase belonging to the GCK family. To elucidate the evolutionary trail of this gene overlapping event, we examined the genomes of a number of primates and found that mutations in the polyadenylation signal of the CHRNE gene in early hominoids led to the overlap. Upon extending this analysis to genomes of other orders of placental mammals, we observed that the overlapping occurred at least three times independently during the course of mammalian evolution. Because CHRNE and MINK are differentially expressed, the potentially hazardous mutations responsible for the exon overlap seem to have escaped evolutionary pressures by differential temporo-spatial expression of the two genes.     

Genetic Architecture Promotes the Evolution and Maintenance of Cooperation

When cooperation has a direct cost and an indirect benefit, a selfish behavior is more likely to be selected for than an altruistic one. Kin and group selection do provide evolutionary explanations for the stability of cooperation in nature, but we still lack the full understanding of the genomic mechanisms that can prevent cheater invasion. In our study we used Aevol, an agent-based, in silico genomic platform to evolve populations of digital organisms that compete, reproduce, and cooperate by secreting a public good for tens of thousands of generations. We found that cooperating individuals may share a phenotype, defined as the amount of public good produced, but have very different abilities to resist cheater invasion. To understand the underlying genetic differences between cooperator types, we performed bio-inspired genomics analyses of our digital organisms by recording and comparing the locations of metabolic and secretion genes, as well as the relevant promoters and terminators. Association between metabolic and secretion genes (promoter sharing, overlap via frame shift or sense-antisense encoding) was characteristic for populations with robust cooperation and was more likely to evolve when secretion was costly. In mutational analysis experiments, we demonstrated the potential evolutionary consequences of the genetic association by performing a large number of mutations and measuring their phenotypic and fitness effects. The non-cooperating mutants arising from the individuals with genetic association were more likely to have metabolic deleterious mutations that eventually lead to selection eliminating such mutants from the population due to the accompanying fitness decrease. Effectively, cooperation evolved to be protected and robust to mutations through entangled genetic architecture. Our results confirm the importance of second-order selection on evolutionary outcomes, uncover an important genetic mechanism for the evolution and maintenance of cooperation, and suggest promising methods for preventing gene loss in synthetically engineered organisms.     

Purifying and directional selection in overlapping prokaryotic genes

In overlapping genes, the same DNA sequence codes for two proteins using different reading frames. Analysis of overlapping genes can help in understanding the mode of evolution of a coding region from noncoding DNA. We identified 71 pairs of convergent genes, with overlapping 3' ends longer than 15 nucleotides, that are conserved in at least two prokaryotic genomes. Among the overlap regions, we observed a statistically significant bias towards the 123:132 phase (i.e. the second codon base in one gene facing the degenerate third position in the second gene). This phase ensures the least mutual constraint on nonconservative amino acid replacements in both overlapping coding sequences. The excess of this phase is compatible with directional (positive) selection acting on the overlapping coding regions. This could be a general evolutionary mode for genes emerging from noncoding sequences, in which the protein sequence has not been subject to selection.     

Is There a Twelfth Protein-Coding Gene in the Genome of Influenza A? A Selection-Based Approach to the Detection of Overlapping Genes in Closely Related Sequences

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.     

Evolution of X-Degenerate Y Chromosome Genes in Greater Apes: Conservation of Gene Content in Human and Gorilla,But Not Chimpanzee

Compared with the X chromosome, the mammalian Y chromosome is considerably diminished in size and has lost most of its ancestral genes during evolution. Interestingly, for the X-degenerate region on the Y chromosome, human has retained all 16 genes, while chimpanzee has lost 4 of the 16 genes since the divergence of the two species. To uncover the evolutionary forces governing ape Y chromosome degeneration, we determined the complete sequences of the coding exons and splice sites for 16 gorilla Y chromosome genes of the X-degenerate region. We discovered that all studied reading frames and splice sites were intact, and thus, this genomic region experienced no gene loss in the gorilla lineage. Higher nucleotide divergence was observed in the chimpanzee than the human lineage, particularly for genes with disruptive mutations, suggesting a lack of functional constraints for these genes in chimpanzee. Surprisingly, our results indicate that the human and gorilla orthologues of the genes disrupted in chimpanzee evolve under relaxed functional constraints and might not be essential. Taking mating patterns and effective population sizes of ape species into account, we conclude that genetic hitchhiking associated with positive selection due to sperm competition might explain the rapid decline in the Y chromosome gene number in chimpanzee. As we found no evidence of positive selection acting on the X-degenerate genes, such selection likely targets other genes on the chimpanzee Y chromosome. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A Method for the Simultaneous Estimation of Selection Intensities in Overlapping Genes

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.     

The role of population size, pleiotropy and fitness effects of mutations in the evolution of overlapping gene functions

Sheltered from deleterious mutations, genes with overlapping or partially redundant functions may be important sources of novel gene functions. While most partially redundant genes originated in gene duplications, it is much less clear why genes with overlapping functions have been retained, in some cases for hundreds of millions of years. A case in point is the many partially redundant genes in vertebrates, the result of ancient gene duplications in primitive chordates. Their persistence and ubiquity become surprising when it is considered that duplicate and original genes often diversify very rapidly, especially if the action of natural selection is involved. Are overlapping gene functions perhaps maintained because of their protective role against otherwise deleterious mutations? There are two principal objections against this hypothesis, which are the main subject of this article. First, because overlapping gene functions are maintained in populations by a slow process of "second order" selection, population sizes need to be very high for this process to be effective. It is shown that even in small populations, pleiotropic mutations that affect more than one of a gene''s functions simultaneously can slow the mutational decay of functional overlap after a gene duplication by orders of magnitude. Furthermore, brief and transient increases in population size may be sufficient to maintain functional overlap. The second objection regards the fact that most naturally occurring mutations may have much weaker fitness effects than the rather drastic "knock-out" mutations that lead to detection of partially redundant functions. Given weak fitness effects of most mutations, is selection for the buffering effect of functional overlap strong enough to compensate for the diversifying force exerted by mutations? It is shown that the extent of functional overlap maintained in a population is not only independent of the mutation rate, but also independent of the average fitness effects of mutation. These results are discussed with respect to experimental evidence on redundant genes in organismal development.     

Stability and evolution of overlapping genes

Abstract.— When the same sequence of nucleotides codes for regions of more than one functional polypeptide, this sequence contains overlapping genes. Overlap is most common in rapidly evolving genomes with high mutation rates such as viruses, bacteria, and mitochondria. Overlap is thought to be important as: (1) a means of compressing a maximum amount of information into short sequences of structural genes; and (2) as a mechanism for regulating gene expression through translational coupling of functionally related polypeptides. The stability of overlapping codes is examined in relation to the information cost of overlap and the mutation rate of the genome. The degree of overlap in a given population will tend to become monomorphic. Evolution toward partial overlap of genes is shown to depend on a convex cost function of overlap. Overlap does not evolve when expression of overlapping genes is mutually exclusive and produced by rare mutations to the wild-type genome. Assuming overlap increases coupling between functionally related genes, the conditions favoring overlap are explored in relation to the kinetics of gene activation and decay. Coupling is most effective for genes in which the gene overlapping at its 5'end (leading gene) decays rapidly, while the gene overlapping at the 3'end (induced gene) decays slowly. If gene expression can feedback on itself (autocatalysis), then high rates of activation favor overlap.     

Conservation of mRNA secondary structures may filter out mutations in Escherichia coli evolution

Recent reports indicate that mutations in viral genomes tend to preserve RNA secondary structure, and those mutations that disrupt secondary structural elements may reduce gene expression levels, thereby serving as a functional knockout. In this article, we explore the conservation of secondary structures of mRNA coding regions, a previously unknown factor in bacterial evolution, by comparing the structural consequences of mutations in essential and nonessential Escherichia coli genes accumulated over 40 000 generations in the course of the ‘long-term evolution experiment’. We monitored the extent to which mutations influence minimum free energy (MFE) values, assuming that a substantial change in MFE is indicative of structural perturbation. Our principal finding is that purifying selection tends to eliminate those mutations in essential genes that lead to greater changes of MFE values and, therefore, may be more disruptive for the corresponding mRNA secondary structures. This effect implies that synonymous mutations disrupting mRNA secondary structures may directly affect the fitness of the organism. These results demonstrate that the need to maintain intact mRNA structures imposes additional evolutionary constraints on bacterial genomes, which go beyond preservation of structure and function of the encoded proteins.     

The cognitive and speech genes are jointly shaped by both positive and relaxed selection in the human lineage

The emergence of a coordinated network of cognitive and speech genes in the human lineage performing overlapping functions is a great evolutionary puzzle. Prior studies on the speech gene FOXP2 are inconclusive on the nature of selection operating on this gene in the human lineage. Here, I show that the evolution of FOXP2 is accelerated in the human lineage due to relaxation of purifying selection (relaxed selection). Five potential genes associated with human-specific intelligence and speech genes have evolved under the impact of positive selection and three genes including FOXP2 have undergone relaxation of purifying selection in the human lineage. Overall, three evolutionary processes namely positive selection, relaxation of purifying selection and neutral evolution have contributed for the genomic evolution of extraordinary cognitive ability and speech in the hominin lineage. The cognitive and speech genes subjected to natural selection in the human lineage have demonstrated a coevolutionary trend.     

Stepwise Evolution and Exceptional Conservation of ORF1a/b Overlap in Coronaviruses

The programmed frameshift element (PFE) rerouting translation from ORF1a to ORF1b is essential for the propagation of coronaviruses. The combination of genomic features that make up PFE—the overlap between the two reading frames, a slippery sequence, as well as an ensemble of complex secondary structure elements—places severe constraints on this region as most possible nucleotide substitution may disrupt one or more of these elements. The vast amount of SARS-CoV-2 sequencing data generated within the past year provides an opportunity to assess the evolutionary dynamics of PFE in great detail. Here, we performed a comparative analysis of all available coronaviral genomic data available to date. We show that the overlap between ORF1a and ORF1b evolved as a set of discrete 7, 16, 22, 25, and 31 nucleotide stretches with a well-defined phylogenetic specificity. We further examined sequencing data from over 1,500,000 complete genomes and 55,000 raw read data sets to demonstrate exceptional conservation and detect signatures of selection within the PFE region.     

Potential evolutionary influences on overlapping reading frames in the bovine leukemia virus pXBL region

Bovine leukemia virus contains a pXBL region encoding the 3' parts of four regulatory proteins (Tax, Rex, G4, R3) in overlapping reading frames. Here we report the pXBL polymorphisms of 30 isolates from four countries. Rates of overall and synonymous substitutions were consistently lower, and nucleotide/amino acid composition bias and codon bias higher, in more-overlapped than in less-overlapped regions. Ratios of nonsynonymous/synonymous substitutions were lowest in the tax gene and its subregions. The 5' parts of the four genes showed selection patterns corresponding to their genomic context outside of the pXBL region. Longer G4 variants due to a natural stop codon mutation had additional triple overlap with reduced sequence variability. These data support the concept that a higher level of overlapping in coding regions correlates with greater evolutionary constraint. Tax, the most conserved among the four regulatory proteins, showed purifying selection consistent with its importance in the viral life cycle.     

Mito-nuclear interactions as drivers of gene movement on and off the X-chromosome

Background

Mito-nuclear gene interactions regulate energy conversion, and are fundamental to eukaryotes. Generally, mito-nuclear coadaptation would be most efficient if the interacting nuclear genes were X-linked, because this maximizes the probability of favorable mito-nuclear allelic combinations co-transmitting across generations. Thus, under a coadaptation (CA) hypothesis, nuclear genes essential for mitochondrial function might be under selection to relocate to the X-chromosome. However, maternal inheritance predisposes the mitochondrial DNA (mtDNA) to accumulate variation that, while male-harming, is benign to females. Numerous nuclear genes were recently reported in Drosophila melanogaster, which exhibit male-specific patterns of differential expression when placed alongside different mtDNA haplotypes, suggesting that nuclear genes are sensitive to an underlying male-specific mitochondrial mutation load. These genes are thus candidates for involvement in mito-nuclear interactions driven by sexual conflict (SC), and selection might have moved them off the X-chromosome to facilitate an optimal evolutionary counter-response, through males, to the presence of male-harming mtDNA mutations. Furthermore, the presence of male-harming mtDNA mutations could exert selection for modifiers on the Y-chromosome, thus placing these mito-sensitive nuclear genes at the center of an evolutionary tug-of-war between mitochondrion and Y-chromosome.We test these hypotheses by examining the chromosomal distributions of three distinct sets of mitochondrial-interacting nuclear genes in D. melanogaster; the first is a list of genes with mitochondrial annotations by Gene Ontologies, the second is a list comprising the core evolutionary-conserved mitochondrial proteome, and the third is a list of genes involved in male-specific responses to maternally-inherited mitochondrial variation and which might be putative targets of Y-chromosomal regulation.

Results

Genes with mitochondrial annotations and genes representing the mitochondrial proteome do not exhibit statistically-significant biases in chromosomal representation. However, genes exhibiting sex-specific sensitivity to mtDNA are under-represented on the X-chromosome, over-represented among genes known to be sensitive to Y-chromosomal variation, and among genes previously associated with male fitness, but under-represented among genes associated with direct sexual antagonism.

Conclusions

Our results are consistent with the SC hypothesis, suggesting that mitochondrial mutational pressure selects for gene movement off-the-X, hence enabling mito-nuclear coadaptation to proceed along trajectories that result in optimized fitness in both sexes.     

Trade-Offs Predicted by Metabolic Network Structure Give Rise to Evolutionary Specialization and Phenotypic Diversification

Mitigating trade-offs between different resource-utilization functions is key to an organism’s ecological and evolutionary success. These trade-offs often reflect metabolic constraints with a complex molecular underpinning; therefore, their consequences for evolutionary processes have remained elusive. Here, we investigate how metabolic architecture induces resource-utilization constraints and how these constraints, in turn, elicit evolutionary specialization and diversification. Guided by the metabolic network structure of the bacterium Lactococcus cremoris, we selected two carbon sources (fructose and galactose) with predicted coutilization constraints. By evolving L. cremoris on either fructose, galactose, or a mix of both sugars, we imposed selection favoring divergent metabolic specializations or coutilization of both resources, respectively. Phenotypic characterization revealed the evolution of either fructose or galactose specialists in the single-sugar treatments. In the mixed-sugar regime, we observed adaptive diversification: both specialists coexisted, and no generalist evolved. Divergence from the ancestral phenotype occurred at key pathway junctions in the central carbon metabolism. Fructose specialists evolved mutations in the fbp and pfk genes that appear to balance anabolic and catabolic carbon fluxes. Galactose specialists evolved increased expression of pgmA (the primary metabolic bottleneck of galactose metabolism) and silencing of ptnABCD (the main glucose transporter) and ldh (regulator/enzyme of downstream carbon metabolism). Overall, our study shows how metabolic network architecture and historical contingency serve to predict targets of selection and inform the functional interpretation of evolved mutations. The elucidation of the relationship between molecular constraints and phenotypic trade-offs contributes to an integrative understanding of evolutionary specialization and diversification.     

Environmental Dependence of Genetic Constraint

The epistatic interactions that underlie evolutionary constraint have mainly been studied for constant external conditions. However, environmental changes may modulate epistasis and hence affect genetic constraints. Here we investigate genetic constraints in the adaptive evolution of a novel regulatory function in variable environments, using the lac repressor, LacI, as a model system. We have systematically reconstructed mutational trajectories from wild type LacI to three different variants that each exhibit an inverse response to the inducing ligand IPTG, and analyzed the higher-order interactions between genetic and environmental changes. We find epistasis to depend strongly on the environment. As a result, mutational steps essential to inversion but inaccessible by positive selection in one environment, become accessible in another. We present a graphical method to analyze the observed complex higher-order interactions between multiple mutations and environmental change, and show how the interactions can be explained by a combination of mutational effects on allostery and thermodynamic stability. This dependency of genetic constraint on the environment should fundamentally affect evolutionary dynamics and affects the interpretation of phylogenetic data.