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

Molecular analysis of the env, LTR, and pX regions of bovine leukemia virus in dairy cattle of Türkiye

Bovine leukemia virus is a retrovirus that causes enzootic bovine leukosis and is associated with global economic losses in the livestock industry. The aim of this study was to investigate the genotype determination of BLVs from cattle housed in 6 different farms in Türkiye and the characterization of their LTR and pX (tax, rex, R3, and G4 gene) regions. For this purpose, blood samples from 48 cattle infected with BLV were used. The phylogenetic analysis based on the env gene sequences revealed that all BLVs were clustered in genotype 1 (G1), and the sequences of the LTR (n = 48) and the pX region (n = 33) of BLVs were obtained. Also, analysis of these nucleic acid and amino acid sequences allowed assessments similar to those reported in earlier studies to be relevant to transactivation and pathogenesis. This study reports the molecular analysis of the LTR and pX region of BLVs in Türkiye for the first time.

Detection and genetic characterization of the novel torque teno virus group 6 in Taiwanese general population

Torque teno virus (TTV) is one of the most common human viruses and can infect an individual with multiple genotypes chronically and persistently. TTV group 6 is a recently discovered phylogenetic group first isolated from eastern Taiwan indigenes, but whether the TTV group 6 was also prevalent in the general population still unknown. One hundred and three randomly collected blood samples from general population and 66 TTV positive DNA samples extracted from Taiwan indigenes were included. A group-6-specific PCR was developed for re-screen over TTV positive samples. Two TTV group 6 positive samples from general population were cloned and sequenced for identifying mix-infected TTVs and confirming their classification by maximum-likelihood and Bayesian inference phylogeny. TTV group 6 can be detected in 4.5% (4/89) and 7.6% (5/66) of TTV positive samples from Taiwanese general population and eastern Taiwan indigenes, respectively. Sample VC09 was mix-infected with TTV groups 3 and 6. Sample VC99 was mix-infected with TTV groups 3, 4 and 6. A highly diverse triple overlapping region was observed, which may represent a unique phenomenon of TTV. The group-6-specific PCR can successfully detect TTV group 6. TTV group 6 may be prevalent worldwide regardless of the geographic region and/or ethnic groups.

Genetic analysis of the pX region of bovine leukemia virus genotype 1 in Holstein Friesian cattle with different stages of infection

The pX genetic region of bovine leukemia virus (BLV) includes four genes with overlapping reading frames that code for the Tax, Rex, R3, and G4 proteins. These proteins are involved in the regulation of transcriptional and post-transcriptional viral expression, as well as having oncogenic potential. Our goal was to investigate the pathogenicity of the pX region of BLV genotype 1 in terms of lymphocytosis, lymphomas, and proviral DNA load. We screened 724 serological samples from mixed-age Holstein Friesian cattle from six states in Mexico. Peripheral blood leukocytes (PBLs) were isolated from whole blood with anticoagulant, and genomic DNA was extracted from the PBLs using a commercial kit. Then, a set of primers that hybridize in conserved regions of the BLV pX region were used, which allowed for PCR standardization to detect proviral DNA in infected cells. Positive amplicons were sequenced using the Sanger method, resulting in 1156-nucleotide-long final sequences that included the four pX region genes. The experimental group consisted of 30 animals. Twelve of these had lymphocytosis, six had lymphoma, and 12 were apparently healthy cattle without any signs of lymphocytosis or lymphoma. The presence of lymphoma was detected in six bovine tumor tissues using histopathology, and the presence of BLV was detected by in situ hybridization. Phylogenetic analysis demonstrated that the 30 sequences were associated with genotype 1, and the genetic distance between the sequences ranged from 0.2% to 2.09%. We identified two sequences in the G4 gene: one with a three-nucleotide deletion resulting in the loss of a leucine (AGU_7488L, in a cow with lymphocytosis), and one with a nine-nucleotide deletion resulting in the loss of leucine, proline, and leucine (AGU_18A, in a cow without lymphocytosis). Analysis of the PX region indicated that positive selection had occurred in the G4, rex, and R3 genes, and we found no difference in proviral DNA load between the studied groups. We were unable to establish an association between variations in the pX region and the development of lymphocytosis, lymphoma, asymptomatic status, or proviral DNA load in BLV-infected cattle.

Detection and molecular characterization of bovine leukemia virus in various regions of Iran

Purpose. Bovine leukemia virus (BLV) infects cattle worldwide, imposing an economic impact on the dairy cattle industry. The purpose of this study was to evaluate the molecular epidemiology of BLV in Iran.Methodology. Blood samples taken from 280 cows aged over 2 years old from 13 provinces of Iran were used for leukocyte count and blocking ELISA. Genomic DNA was extracted from the peripheral blood leukocytes of BLV-infected samples and fetal lamb kidney cells to perform PCR of partial env, rex and tax genes and long-terminal-repeat region. The PCR products were sequenced, the phylogenetic tree of each gene was constructed, and nucleotide and amino acid sequence pair distances were calculated.Results. The frequency of BLV infection was 32.8 % among animals and was 80 % among provinces. In BLV seropositive animals, the rate of persistent lymphocytosis was 36.9 %. The constructed phylogenetic trees showed the presence of two BLV genotypes (1 and 4) in Iranian strains. As previous studies, our results showed that the env gene was more variable than previously thought, the Rex protein could withstand more amino acid changes compared to the Tax protein, and no significant differences were observed in average changes of the nucleotide of these genes between clinical stages.Conclusions. Our data indicates an increase in the frequency of this infection in Iran. This is the first study report of the presence of BLV genotype 4 in Iranian farms. These findings may have an important role in the control and prevention of BLV infection in Iran and other countries.

Bovine leukemia virus discovered in human blood

BACKGROUND

Bovine leukemia virus (BLV) infection is widespread in cattle globally and is present in marketed beef and dairy products. Human infection with BLV has been reported in breast and lung cancer tissues and was significantly associated with breast cancer in 3 case-control studies. The purpose of this current research was to determine if BLV is present in human blood cells and if antibodies to BLV are related to blood cell infection.

METHODS

Standard liquid PCR and Sanger DNA sequencing were used to test for BLV in buffy coat cells (leukocytes and platelets) of blood specimens from 95 self-selected female subjects. Enzyme-linked immunosorbent assay (ELISA) for IgG, IgM, and IgA was used to detect antibodies to BLV in the plasma of the corresponding blood samples.

RESULTS

BLV DNA was detected in the buffy coat cells of blood in 33/95 (38%) of the subjects by PCR and DNA sequencing. IgG antibodies were detected in 30/95(32%), IgM in 55/95(58%), and IgA in 30/95(32%) of the subjects. There was no significant correlation between presence of the antibodies and presence of BLV DNA.

CONCLUSIONS

This first report of BLV in human blood raises the question of whether infection of leukocytes could conceivably lead to leukemia as it does in infected cattle. Also, system wide circulation of infected blood cells could facilitate BLV transit to various internal tissues/organs with potential for their infection and subsequent development of cancer. The most likely route of BLV transmission to humans would be zoonotic, as a foodborne infection. Although eradicated from cattle in some countries, BLV still has a high rate of infection in the Americas, the Middle East, and parts of Europe and Asia. This report of BLV in the blood layer containing human leukocytes/platelets adds important information which could be useful to elucidate possible routes of transmission of BLV to humans and to prevent further human infection.

Bovine leukemia virus long terminal repeat variability: identification of single nucleotide polymorphisms in regulatory sequences

Background

Limited data are available on the incidence of variations in nucleotide sequences of long terminal repeat (LTR) regions of Bovine Leukemia Virus (BLV). Consequently, the possible impact of SNPs on BLV LTR function are poorly elucidated. Thus, a detailed and representative study of full-length LTR sequences obtained from sixty-four BLV isolates from different geographical regions of Poland, Moldova, Croatia, Ukraine and Russia were analyzed for their genetic variability.

Methods

Overlap extension PCR, sequencing and Bayesian phylogenetic reconstruction of LTR sequences were performed. These analyses were followed by detailed sequence comparison, estimation of genetic heterogeneity and identification of transcription factor binding site (TFBS) modifications.

Results

Phylogenetic analysis of curated LTR sequences and those available in the GenBank database reflected the acknowledged env gene classification of BLV into 10 genotypes, and further clustered analysed sequences into three genotypes - G4, G7 and G8. Additional molecular studies revealed the presence of 97 point mutations distributed at 89 positions throughout all 64 LTR sequences. The highest rate of variability was noted in U3 and U5 subregions. However, the variability in regulatory sequences (V_R) was assessed as lower than the variability within non-regulatory sequences (V_NR) for both, U3 and U5 subregions. In contrast, V_R value for R subregion, as well as for the total LTR, was higher than the V_NR suggesting the existence of positive selection. Twelve unique SNPs for these LTR sequences localized in regulatory and non-regulatory elements were identified. The presence of different types of substitutions lead to the abrogation of present or to the creation of additional TFBS.

Conclusion

This study represents the largest study of LTR genetic variability of BLV field isolates from Eastern part of Europe. Phylogenetic analysis of LTRs supports the clustering BLV variants based on their geographic origin. The SNP screening showed variations modifying LTR regulatory sequences, as well as altering TFBS. These features warrant further exploration as they could be related to proviral load and distinctive regulation of BLV transcription and replication.

Electronic supplementary material

The online version of this article (10.1186/s12985-018-1062-z) contains supplementary material, which is available to authorized users.

Overlapping genes and the proteins they encode differ significantly in their sequence composition from non-overlapping genes

Overlapping genes represent a fascinating evolutionary puzzle, since they encode two functionally unrelated proteins from the same DNA sequence. They originate by a mechanism of overprinting, in which point mutations in an existing frame allow the expression (the "birth") of a completely new protein from a second frame. In viruses, in which overlapping genes are abundant, these new proteins often play a critical role in infection, yet they are frequently overlooked during genome annotation. This results in erroneous interpretation of mutational studies and in a significant waste of resources. Therefore, overlapping genes need to be correctly detected, especially since they are now thought to be abundant also in eukaryotes. Developing better detection methods and conducting systematic evolutionary studies require a large, reliable benchmark dataset of known cases. We thus assembled a high-quality dataset of 80 viral overlapping genes whose expression is experimentally proven. Many of them were not present in databases. We found that overall, overlapping genes differ significantly from non-overlapping genes in their nucleotide and amino acid composition. In particular, the proteins they encode are enriched in high-degeneracy amino acids and depleted in low-degeneracy ones, which may alleviate the evolutionary constraints acting on overlapping genes. Principal component analysis revealed that the vast majority of overlapping genes follow a similar composition bias, despite their heterogeneity in length and function. Six proven mammalian overlapping genes also followed this bias. We propose that this apparently near-universal composition bias may either favour the birth of overlapping genes, or/and result from selection pressure acting on them.

Computational analysis of envelope glycoproteins from diverse geographical isolates of bovine leukemia virus identifies highly conserved peptide motifs

BACKGROUND

Bovine leukemia virus (BLV) is a deltaretrovirus infecting bovine B cells and causing enzootic bovine leucosis. The SU or surface subunit, gp51, of its envelope glycoprotein is involved in receptor recognition and virion attachment. It contains the major neutralizing and CD4+ and CD8+ T cell epitopes found in naturally infected animals. In this study, we aimed to determine global variation and conservation within gp51 in the context of developing an effective global BLV vaccine.

RESULTS

A total of 256 sequences extracted from the NCBI database and collected in different parts of the world, were studied to identify conserved segments along the env gene sequences that encode the gp51 protein. Using the MEME server and the conserved DNA Region module for analysis within DnaSP, we identified six conserved segments, referred to as A-F, and five semi-conserved segments, referred to as G-K. The amino acid conservation ranged from 98.8 to 99.8% in conserved segments A to F, while segments G to K had 89.6-95.2% conserved amino acid sequence. Selection analysis of individual segments revealed that residues of conserved segments had undergone purifying selection, whereas, particular residues in the semi-conserved segments are currently undergoing positive selection, specifically at amino acid positions 48 in segment K, 74 in segment G, 82 in segment I, 133 and 142 in segment J, and residue 291 in segment H. Each of the codons for these six residues contain the most highly variable nucleotides within their respective semi-conserved segments.

CONCLUSIONS

The data described here show that the consensus amino acid sequence constitutes a strong candidate from which a global vaccine can be derived for use in countries where eradication by culling is not economically feasible. The most conserved segments overlap with amino acids in known immunodeterminants, specifically in epitopes D-D', E-E', CD8+ T-cell epitopes, neutralizing domain 1 and CD4+ T-cell epitopes. Two of the segments reported here represent unique segments that do not overlap with previously identified antigenic determinants. We propose that evidence of positive selection in some residues of the semi-conserved segments suggests that their variation is involved in viral strategy to escape immune surveillance of the host.

Overlapping Regions in HIV-1 Genome Act as Potential Sites for Host-Virus Interaction

More than a decade, overlapping genes in RNA viruses became a subject of research which has explored various effect of gene overlapping on the evolution and function of viral genomes like genome size compaction. Additionally, overlapping regions (OVRs) are also reported to encode elevated degree of protein intrinsic disorder (PID) in unspliced RNA viruses. With the aim to explore the roles of OVRs in HIV-1 pathogenesis, we have carried out an in-depth analysis on the association of gene overlapping with PID in 35 HIV1- M subtypes. Our study reveals an over representation of PID in OVR of HIV-1 genomes. These disordered residues endure several vital, structural features like short linear motifs (SLiMs) and protein phosphorylation (PP) sites which are previously shown to be involved in massive host–virus interaction. Moreover, SLiMs in OVRs are noticed to be more functionally potential as compared to that of non-overlapping region. Although, density of experimentally verified SLiMs, resided in 9 HIV-1 genes, involved in host–virus interaction do not show any bias toward clustering into OVR, tat and rev two important proteins mediates host–pathogen interaction by their experimentally verified SLiMs, which are mostly localized in OVR. Finally, our analysis suggests that the acquisition of SLiMs in OVR is mutually exclusive of the occurrence of disordered residues, while the enrichment of PPs in OVR is solely dependent on PID and not on overlapping coding frames. Thus, OVRs of HIV-1 genomes could be demarcated as potential molecular recognition sites during host–virus interaction.

Exposure to Bovine Leukemia Virus Is Associated with Breast Cancer: A Case-Control Study

BACKGROUND

Age, reproductive history, hormones, genetics, and lifestyle are known risk factors for breast cancer, but the agents that initiate cellular changes from normal to malignant are not understood. We previously detected bovine leukemia virus (BLV), a common oncogenic virus of cattle, in the breast epithelium of humans. The objective of this study was to determine whether the presence of BLV DNA in human mammary epithelium is associated with breast cancer.

METHODS

This was a case-control study of archival formalin fixed paraffin embedded breast tissues from 239 donors, received 2002-2008 from the Cooperative Human Tissue Network. Case definition as breast cancer versus normal (women with no history of breast cancer) was established through medical records and examination of tissues by an anatomical pathologist. Breast exposure to BLV was determined by in situ-PCR detection of a biomarker, BLV DNA, localized within mammary epithelium.

RESULTS

The frequency of BLV DNA in mammary epithelium from women with breast cancer (59%) was significantly higher than in normal controls (29%) (multiply- adjusted odds ratio = 3.07, confidence interval = 1.66-5.69, p = .0004, attributable risk = 37%). In women with premalignant breast changes the frequency of BLV DNA was intermediate (38%) between that of women with breast cancer and normal controls (p for trend < .001).

CONCLUSIONS

Among the specimens in this study, the presence of amplified BLV DNA was significantly associated with breast cancer. The odds ratio magnitude was comparable to those of well-established breast cancer risk factors related to reproductive history, hormones, and lifestyle and was exceeded only by risk factors related to genetics (familial breast cancer), high dose ionizing radiation, and age. These findings have the potential for primary and secondary prevention of breast cancer.

Viral proteins originated de novo by overprinting can be identified by codon usage: application to the "gene nursery" of Deltaretroviruses

A well-known mechanism through which new protein-coding genes originate is by modification of pre-existing genes, e.g. by duplication or horizontal transfer. In contrast, many viruses generate protein-coding genes de novo, via the overprinting of a new reading frame onto an existing (“ancestral”) frame. This mechanism is thought to play an important role in viral pathogenicity, but has been poorly explored, perhaps because identifying the de novo frames is very challenging. Therefore, a new approach to detect them was needed. We assembled a reference set of overlapping genes for which we could reliably determine the ancestral frames, and found that their codon usage was significantly closer to that of the rest of the viral genome than the codon usage of de novo frames. Based on this observation, we designed a method that allowed the identification of de novo frames based on their codon usage with a very good specificity, but intermediate sensitivity. Using our method, we predicted that the Rex gene of deltaretroviruses has originated de novo by overprinting the Tax gene. Intriguingly, several genes in the same genomic region have also originated de novo and encode proteins that regulate the functions of Tax. Such “gene nurseries” may be common in viral genomes. Finally, our results confirm that the genomic GC content is not the only determinant of codon usage in viruses and suggest that a constraint linked to translation must influence codon usage.

How does novelty originate in nature? It is commonly thought that new genes are generated mainly by modifications of existing genes (the “tinkering” model). In contrast, we have shown recently that in viruses, numerous genes are generated entirely de novo (“from scratch”). The role of these genes remains underexplored, however, because they are difficult to identify. We have therefore developed a new method to detect genes originated de novo in viral genomes, based on the observation that each viral genome has a unique “signature”, which genes originated de novo do not share. We applied this method to analyze the genes of Human T-Lymphotropic Virus 1 (HTLV1), a relative of the HIV virus and also a major human pathogen that infects about twenty million people worldwide. The life cycle of HTLV1 is finely regulated – it can stay dormant for long periods and can provoke blood cancers (leukemias) after a very long incubation. We discovered that several of the genes of HTLV1 have originated de novo. These novel genes play a key role in regulating the life cycle of HTLV1, and presumably its pathogenicity. Our investigations suggest that such “gene nurseries” may be common in viruses.

The effect of gene overlapping on the rate of RNA virus evolution

Gene overlapping is widely employed by RNA viruses to generate genetic novelty while retaining a small genome size. However, gene overlapping also increases the deleterious effect of mutations as they affect more than one gene, thereby reducing the evolutionary rate of RNA viruses and hence their adaptive capacity. Although there is general agreement on the benefits of gene overlapping as a mechanism of genomic compression for rapidly evolving organisms, its effect on the pace of RNA virus evolution remains a source of debate. To address this issue, we collected sequence data from 117 instances of gene overlapping across 19 families, 30 genera, and 55 species of RNA viruses. On these data, we analyzed how genetic distances, selective pressures, and the distribution of RNA secondary structures and conserved protein functional domains vary between overlapping (OV) and nonoverlapping (NOV) regions. We show that gene overlapping generally results in a decrease in the rate of RNA virus evolution through a reduction in the frequency of synonymous mutations. However, this effect is less pronounced in genes with a terminal rather than an internal gene overlap, which might result from a greater proportion of protein functional conserved domains in NOV than in OV regions, in turn reducing the number of nonsynonymous mutations in the former. Overall, our analyses clarify the role of gene overlapping as a modulator of the evolutionary rates exhibited by RNA viruses and shed light on the factors that shape the genetic diversity of this important group of pathogens.

The molecular characterization of bovine leukaemia virus isolates from Eastern Europe and Siberia and its impact on phylogeny

Recent studies have shown that bovine leukemia virus (BLV) sequences can be classified into seven distinct genotypes based on full gp51 sequence. This classification was based on available sequence data that mainly represented the BLV population that is circulating in cattle from the US and South America. In order to aid with a global perspective inclusion of data from Eastern Europe is required. In this study we examined 44 BLV isolates from different geographical regions of Poland, Belarus, Ukraine, and Russia. Phylogenetic analysis based on a 444bp fragment of env gene revealed that most of isolates belonged to genotypes 4 and 7. Furthermore, we confirmed the existence of a new genotype, genotype 8, which was highly supported by phylogenetic computations. A significant number of amino acid substitutions were found in the sequences of the studied Eastern European isolates, of which 71% have not been described previously. The substitutions encompassed mainly the C-part of the CD4+ epitope, zinc binding peptide region, CD8+ T cell epitope, and overlapping linear epitope E. These observations highlight the use of sequence data to both elucidate phylogenetic relationships and the potential effect on serological detection of geographically diverse isolates.

Probable presence of an ubiquitous cryptic mitochondrial gene on the antisense strand of the cytochrome oxidase I gene

BACKGROUND

Mitochondria mediate most of the energy production that occurs in the majority of eukaryotic organisms. These subcellular organelles contain a genome that differs from the nuclear genome and is referred to as mitochondrial DNA (mtDNA). Despite a disparity in gene content, all mtDNAs encode at least two components of the mitochondrial electron transport chain, including cytochrome c oxidase I (Cox1).

PRESENTATION OF THE HYPOTHESIS

A positionally conserved ORF has been found on the complementary strand of the cox1 genes of both eukaryotic mitochondria (protist, plant, fungal and animal) and alpha-proteobacteria. This putative gene has been named gau for gene antisense ubiquitous in mtDNAs. The length of the deduced protein is approximately 100 amino acids. In vertebrates, several stop codons have been found in the mt gau region, and potentially functional gau regions have been found in nuclear genomes. However, a recent bioinformatics study showed that several hypothetical overlapping mt genes could be predicted, including gau; this involves the possible import of the cytosolic AGR tRNA into the mitochondria and/or the expression of mt antisense tRNAs with anticodons recognizing AGR codons according to an alternative genetic code that is induced by the presence of suppressor tRNAs. Despite an evolutionary distance of at least 1.5 to 2.0 billion years, the deduced Gau proteins share some conserved amino acid signatures and structure, which suggests a possible conserved function. Moreover, BLAST analysis identified rare, sense-oriented ESTs with poly(A) tails that include the entire gau region. Immunohistochemical analyses using an anti-Gau monoclonal antibody revealed strict co-localization of Gau proteins and a mitochondrial marker.

TESTING THE HYPOTHESIS

This hypothesis could be tested by purifying the gau gene product and determining its sequence. Cell biological experiments are needed to determine the physiological role of this protein.

IMPLICATIONS OF THE HYPOTHESIS

Studies of the gau ORF will shed light on the origin of novel genes and their functions in organelles and could also have medical implications for human diseases that are caused by mitochondrial dysfunction. Moreover, this strengthens evidence for mitochondrial genes coded according to an overlapping genetic code.

Long-term evolution of the Luteoviridae: time scale and mode of virus speciation

Despite their importance as agents of emerging disease, the time scale and evolutionary processes that shape the appearance of new viral species are largely unknown. To address these issues, we analyzed intra- and interspecific evolutionary processes in the Luteoviridae family of plant RNA viruses. Using the coat protein gene of 12 members of the family, we determined their phylogenetic relationships, rates of nucleotide substitution, times to common ancestry, and patterns of speciation. An associated multigene analysis enabled us to infer the nature of selection pressures and the genomic distribution of recombination events. Although rates of evolutionary change and selection pressures varied among genes and species and were lower in some overlapping gene regions, all fell within the range of those seen in animal RNA viruses. Recombination breakpoints were commonly observed at gene boundaries but less so within genes. Our molecular clock analysis suggested that the origin of the currently circulating Luteoviridae species occurred within the last 4 millennia, with intraspecific genetic diversity arising within the last few hundred years. Speciation within the Luteoviridae may therefore be associated with the expansion of agricultural systems. Finally, our phylogenetic analysis suggested that viral speciation events tended to occur within the same plant host species and country of origin, as expected if speciation is largely sympatric, rather than allopatric, in nature.

Mixed infections of Pepino mosaic virus strains modulate the evolutionary dynamics of this emergent virus

Pepino mosaic virus (PepMV) is an emerging pathogen that causes severe economic losses in tomato crops (Solanum lycopersicum L.) in the Northern hemisphere, despite persistent attempts of control. In fact, it is considered one of the most significant viral diseases for tomato production worldwide, and it may constitute a good model for the analysis of virus emergence in crops. We have combined a population genetics approach with an analysis of in planta properties of virus strains to explain an observed epidemiological pattern. Hybridization analysis showed that PepMV populations are composed of isolates of two types (PepMV-CH2 and PepMV-EU) that cocirculate. The CH2 type isolates are predominant; however, EU isolates have not been displaced but persist mainly in mixed infections. Two molecularly cloned isolates belonging to each type have been used to examine the dynamics of in planta single infections and coinfection, revealing that the CH2 type has a higher fitness than the EU type. Coinfections expand the range of susceptible hosts, and coinfected plants remain symptomless several weeks after infection, so a potentially important problem for disease prevention and management. These results provide an explanation of the observed epidemiological pattern in terms of genetic and ecological interactions among the different viral strains. Thus, mixed infections appear to be contributing to shaping the genetic structure and dynamics of PepMV populations.

Natural genetic variations in bovine leukemia virus envelope gene: Possible effects of selection and escape

Bovine leukemia virus (BLV) is an oncogenic virus widespread in cattle. It belongs to the genus Deltaretrovirus of the family Retroviridae along with human and simian T-lymphotropic viruses. Here we report the addition of 28 new sequences to the current literature of 16 full-length BLV envelope gene sequences. The phylogenetic clustering, genotyping, and geographic distribution of BLV env variations corresponded in most cases. Most natural variations are mapped to the surface of the proposed conformational models of BLV gp51 N-terminus and gp30 external domain, overlapping with or adjacent to immunogenic epitopes. Analyses for evidence of possible selection pressures suggest the BLV env is under stringent negative selection overall, while strong positive selection is indicated for immunogenic epitope G. Natural env deletions bounded by similar flanking sequences were observed in multiple isolates and would result in truncated signal peptides, missing gp51, and aberrant coding frames for other proteins.

Pattern of nucleotide substitution in the overlapping nonstructural genes of influenza A virus and implication for the genetic diversity of the H5N1 subtype

In viruses under strong pressure to minimize genome size, overlapping genes represent a fine strategy to condense a maximum amount of information into short nucleotide sequences. Here, we investigated the evolution of the genes encoding the nonstructural proteins NS1 and NS2 of influenza A virus (IAV), which are one of the best characterized cases of gene overlap. By a detailed analysis of about four hundred sequences grouped into 11 IAV subtypes, we found that the overlapping coding region of the NS1 gene shows a significant increase of the rate of nonsynonymous change, with respect to its nonoverlapping counterpart. The same feature was observed in the overlapping coding region of the NS2 gene. Such a variation pattern, which implies the occurrence of several amino acid substitutions in the protein regions encoded by overlapping frames, is different from the pattern of constrained evolution typical of other viral overlapping-gene systems. Amino acid sequence analysis of the NS1 and NS2 proteins revealed that some nonsynonymous substitutions, located in the region of gene overlap, play a critical role in shaping the genetic diversity of the highly pathogenic subtype H5N1. Since both proteins contribute to disease pathogenesis by affecting many virus and host-cell processes, information provided by this study should be useful to highlight the impact of nonstructural gene variation on the pathogenicity of H5N1 viruses.