An overview of dinucleotide and codon usage in all viruses.

Viruses are, by far, the most abundant biological entities on earth. They are found in all known ecological niches and are the causative agents of many important diseases in plants and animals. From an evolutionary point of view, since viruses do not share any orthologous genes, there is a general consensus that they are polyphyletic; that is, they do not have a common ancestor. This means that they appeared several times during the course of evolution. For their life cycle, they are always obligate parasites of a free cellular life form, which can be bacteria, archaea, or eukaryotes. More complexity is added to these entities by the fact that their genetic material can be DNA or RNA (double- or single-stranded) or retrotranscribed. Given these features, we wondered if some general rules can be inferred when studying two basic genomic signatures-dinucleotides and codon usage-analyzing all available complete and non-redundant viral sequences. In spite of the obviously biased sample of sequences available, some general features appear to emerge.

Resistance-associated substitutions and response to treatment in a chronic hepatitis C virus infected-patient: an unusual virological response case report.

BACKGROUND: Direct-Acting agents (DAAs) target and inhibit essential viral replication proteins. They have revolutionized the treatment of Hepatitis C virus (HCV) infection reaching high levels of sustained virologic response. However, the detection of basal resistance-associated substitutions (RASs) to DAAs in naïve patients could be important in predicting the treatment outcome in some patients exhibiting failures to DAA-based therapies. Therefore, the aim of this work was to evaluate the presence of RASs as minority variants within intra-host viral populations, and assess their relationship to response to therapy on a multiple times relapser patient infected chronically with HCV. CASE PRESENTATION: A male HCV infected-patient with a genotype 1a strain was evaluated. He had previously not responded to dual therapy (pegylated interferon-α plus ribavirin) and was going to start a direct-acting agent-based therapy (DAAs). He showed no significant liver fibrosis (F0). Viral RNA was extracted from serum samples taken prior and after therapy with DAAs (sofosbubir/ledipasvir/ribavirin). NS5A and NS5B genomic regions were PCR-amplified and the amplicons were sequenced using Sanger and next-generation sequencing (NGS) approaches. RASs were searched in in-silico translated sequences for all DAAs available and their frequencies were determined for those detected by NGS technology. Sanger sequencing did not reveal the presence of RASs in the consensus sequence neither before nor after the DAA treatment. However, several RASs were found at low frequencies, both before as well as after DAA treatment. RASs found as minority variants (particularly substitutions in position 93 within NS5A region) seem to have increased their frequency after DAA pressure. Nevertheless, these RASs did not become dominant and the patient still relapsed, despite perfect adherence to treatment and having no other complications beyond the infection (no significant fibrosis, no drug abuse). CONCLUSIONS: This report shows that some patients might relapse after a DAA-based therapy even when RASs (pre- and post-treatment) are detected in very low frequencies (< 1%) within intra-host viral populations. Increased awareness of this association may improve detection and guide towards a personalized HCV treatment, directly improving the outcome in hard-to-treat patients.

Evidence of increasing diversification of emerging Severe Acute Respiratory Syndrome Coronavirus 2 strains.

On 30th January 2020, an outbreak of atypical pneumonia caused by a novel betacoronavirus, named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was declared a public health emergency of international concern by the World Health Organization. For this reason, a detailed evolutionary analysis of SARS-CoV-2 strains currently circulating in different geographic regions of the world was performed. A compositional analysis as well as a Bayesian coalescent analysis of complete genome sequences of SARS-CoV-2 strains recently isolated in Europe, North America, South America, and Asia was performed. The results of these studies revealed a diversification of SARS-CoV-2 strains in three different genetic clades. Co-circulation of different clades in different countries, as well as different genetic lineages within different clades were observed. The time of the most recent common ancestor was established to be around 1st November 2019. A mean rate of evolution of 6.57 × 10-4 substitutions per site per year was found. A significant migration rate per genetic lineage per year from Europe to South America was also observed. The results of these studies revealed an increasing diversification of SARS-CoV-2 strains. High evolutionary rates and fast population growth characterizes the population dynamics of SARS-CoV-2 strains.

An evolutionary insight into emerging Ebolavirus strains isolated in Africa.

On July 19, 2019, the World Health Organization declared the current Ebolavirus (EBOV) outbreak in Congo Democratic Republic (COD) a public health emergency of international concern. To address the potential threat of EBOV evolution outpacing antibody treatment and vaccine efforts, a detailed evolutionary analysis of EBOV strains circulating in different African countries was performed. Genome composition of EBOV strains was studied using multivariate statistical analysis. To investigate the patterns of evolution of EBOV strains, a Bayesian Markov Chain Monte Carlo approach was used. Two different genetic lineages, with a distinct genome composition gave rise to the recent EBOV outbreaks in central and western Africa. Strains isolated in COD in 2018 fall into two different genetic clusters, according to their geographical location of isolation. Different amino acid substitutions among strains from these two clusters have been found, particularly in NP, GP, and L proteins. Significant differences in codon and amino acid usage among clusters were found. Strains isolated in COD in 2018 belong to two distinct genetic clusters, with distinct codon and amino acid usage. Geographical diversity plays an important role in shaping the molecular evolution of EBOV populations.

IL28B gene polymorphism rs12979860, but not rs8099917, contributes to the occurrence of chronic HCV infection in Uruguayan patients.

BACKGROUND: Host single-nucleotide polymorphisms (SNPs) near the interleukin 28B (IL28B) locus are associated with sustained virological response to antiviral therapy and with spontaneous Hepatitis C Virus (HCV) clearance. Prevalence of these SNPs varies depending on ethnicity. The impact of IL28B SNPs in HCV-infected patients is currently unknown in Uruguay. Therefore, the aim of this study was to evaluate and compare the distribution of polymorphisms in the IL28B gene (rs12979860 and rs8099917) among HCV-infected patients and healthy individuals in Uruguay and thus assess their possible association with the establishment of HCV infection. METHODS: DNA was recovered from 92 non-infected individuals and 78 HCV-infected patients and SNPs were determined by RFLP and allelic discrimination by real-time PCR. RESULTS: The distribution of rs12979860 genotypes for the infected population was 29.5%-CC, 47.4%-CT and 23.1%-TT and for the control group 45.7%, 42.4% and 11.9%, respectively. Prevalence in both infected and uninfected individuals is similar to that reported in other countries with admixed populations. The distribution of rs8099917 genotypes for the infected population was 57.7%-TT, 27.2%-TG and 14.1%-GG and for the control group 60.9%, 33.7% and 5.4%, respectively. The comparison of rs12979860 genotype distribution between the two populations evidenced a higher prevalence of the favourable genotype (CC) in the uninfected control group (p < 0.05). Additionally, results generated using logistic regression analysis show that individuals carrying rs12979860-TT or CT genotypes have a higher likelihood of developing chronic hepatitis upon infection with HCV, when compared to CC carriers, considering rs8099917 genotype as constant. CONCLUSION: Patients with HCV infection have a statistically significant lower prevalence of the favourable rs12979860 genotype when compared to uninfected individuals; therefore we can establish that only IL28B rs12979860-CT and TT genotypes seem to contribute to the occurrence of chronic HCV infection in the cohort of Uruguayan population studied. Considering that a trend towards a higher frequency of "good" response genotypes was observed in responder patients, we believe that IL28B rs12979860 genotyping could be a useful tool for predicting different therapies outcome, including in the DAA era.

Evidence of increasing diversification of Zika virus strains isolated in the American continent.

Zika virus (ZIKV) is a member of the family Flaviviridae. ZIKV emerged in Brazil in 2015, causing an unprecedented epidemic and since then the virus has rapidly spread throughout the Americas. These facts highlight the need of detailed phylogenetic studies to understand the emergence, spread, and evolution of ZIKV populations. For these reasons, a Bayesian coalescent Markov Chain Monte Carlo analysis of complete genome sequences of ZIKV strains recently isolated in the American continent was performed. The results of these studies revealed an increasing diversification of ZIKV strains in different genetic lineages and co-circulation of distinct genetic lineages in several countries in the region. The time of the most recent common ancestor (tMRCA) was established to be around February 20, 2014 for ZIKV strains circulating in the American region. A mean rate of evolution of 1.55 × 10-3 substitutions/site/year was obtained for ZIKV strains included in this study. A Bayesian skyline plot indicate a sharp increase in population size from February 2014 to July 2015 and a decline during 2016. These results are discussed in terms of the emergence and evolution of ZIKV populations in the American continent.

Genome-wide analysis of codon usage bias in Bovine Coronavirus.

BACKGROUND: Bovine coronavirus (BCoV) belong to the genus Betacoronavirus of the family Coronaviridae. BCoV are widespread around the world and cause enteric or respiratory infections among cattle, leading to important economic losses to the beef and dairy industry worldwide. To study the relation of codon usage among viruses and their hosts is essential to understand host-pathogen interaction, evasion from host's immune system and evolution. METHODS: We performed a comprehensive analysis of codon usage and composition of BCoV. RESULTS: The global codon usage among BCoV strains is similar. Significant differences of codon preferences in BCoV genes in relation to codon usage of Bos taurus host genes were found. Most of the highly frequent codons are U-ending. G + C compositional constraint and dinucleotide composition also plays a role in the overall pattern of BCoV codon usage. CONCLUSIONS: The results of these studies revealed that mutational bias is a leading force shaping codon usage in this virus. Additionally, relative dinucleotide frequencies, geographical distribution, and evolutionary processes also influenced the codon usage pattern.

Bayesian coalescent inference reveals high evolutionary rates and diversification of Zika virus populations.

Zika virus (ZIKV) is a member of the family Flaviviridae. In 2015, ZIKV triggered an epidemic in Brazil and spread across Latin America. By May of 2016, the World Health Organization warns over spread of ZIKV beyond this region. Detailed studies on the mode of evolution of ZIKV strains are extremely important for our understanding of the emergence and spread of ZIKV populations. In order to gain insight into these matters, a Bayesian coalescent Markov Chain Monte Carlo analysis of complete genome sequences of recently isolated ZIKV strains was performed. The results of these studies revealed a mean rate of evolution of 1.20 × 10(-3) nucleotide substitutions per site per year (s/s/y) for ZIKV strains enrolled in this study. Several variants isolated in China are grouped together with all strains isolated in Latin America. Another genetic group composed exclusively by Chinese strains were also observed, suggesting the co-circulation of different genetic lineages in China. These findings indicate a high level of diversification of ZIKV populations. Strains isolated from microcephaly cases do not share amino acid substitutions, suggesting that other factors besides viral genetic differences may play a role for the proposed pathogenesis caused by ZIKV infection. J. Med. Virol. 88:1672-1676, 2016. © 2016 Wiley Periodicals, Inc.

Molecular epidemiology of group A rotavirus among children admitted to hospital in Salto, Uruguay, 2011-2012: first detection of the emerging genotype G12.

Group A rotavirus (RVA) is the most important etiologic agent of infant acute gastroenteritis (AGE) worldwide. Detection and molecular characterization of RVA in Salto department, Northwestern region of Uruguay, was conducted on 175 clinical samples, being 153 stool and 22 vomit samples, collected from hospitalized children with AGE, between 0-15 years old, from two hospitals of Salto city during 2011 and 2012. RVA was detected and genotyped by seminested multiplex RT-PCR in order to determine G- and P-genotypes. Positive samples were sequenced and phylogenetic analyses were carried out in order to determine lineages and sub-lineages. RVA were detected in 64 (37%) of the samples and the G and P genotypes observed were: 6% G1P[8], 23% G2P[4]/G2P[X]/GXP[4], 23% G3P[8]/G3P[X], 14% G12P[8]/G12P[X], 16% GXP[8], 1,5% G12P[9], 3% G2P[4]/[8], and 16% non-typeable. VP7 and VP4 genotypes related to DS-1 like gene constellation were prevalent during 2011 and those VP7 and VP4 genotypes related to Wa-like constellation were prevalent during 2012 (mainly represented by G3P[8]). Interestingly, RVA was detected in vomit samples in a high prevalence (41%). RVA was observed mainly in the age group between 1 and 5 years old (75% of the cases), and seasonality with a high detection rate in winter season was observed for the two consecutive years of surveillance. To our knowledge, this study represents the first detection and molecular characterization of RVA in Salto department, Northwestern region of Uruguay; and the first identification of the emerging genotype G12 in the country.

An evolutionary insight into Newcastle disease viruses isolated in Antarctica.

The disease caused by Newcastle disease virus (NDV) is a severe threat to the poultry industry worldwide. Recently, NDV has been isolated in the Antarctic region. Detailed studies on the mode of evolution of NDV strains isolated worldwide are relevant for our understanding of the evolutionary history of NDV. For this reason, we have performed Bayesian coalescent analysis of NDV strains isolated in Antarctica to study evolutionary rates, population dynamics, and patterns of evolution. Analysis of F protein cleavage-site sequences of NDV isolates from Antarctica suggested that these strains are lentogenic. Strains isolated in Antarctica and genotype I reference strain Ulster/67 diverged from ancestors that existed around 1958. The time of the most recent common ancestor (MRCA) was established to be around 1883 for all class II viruses. A mean rate of evolution of 1.78 × 10(-3) substitutions per site per year (s/s/y) was obtained for the F gene sequences of NDV strains examined in this study. A Bayesian skyline plot indicated a decline in NDV population size in the last 25 years. The results are discussed in terms of the possible role of Antarctica in emerging or re-emerging viruses and the evolution of NDV populations worldwide.

Mixed Infections Unravel Novel HCV Inter-Genotypic Recombinant Forms within the Conserved IRES Region.

Hepatitis C virus (HCV) remains a significant global health challenge, affecting millions of people worldwide, with chronic infection a persistent threat. Despite the advent of direct-acting antivirals (DAAs), challenges in diagnosis and treatment remain, compounded by the lack of an effective vaccine. The HCV genome, characterized by high genetic variability, consists of eight distinct genotypes and over ninety subtypes, underscoring the complex dynamics of the virus within infected individuals. This study delves into the intriguing realm of HCV genetic diversity, specifically exploring the phenomenon of mixed infections and the subsequent detection of recombinant forms within the conserved internal ribosome entry site (IRES) region. Previous studies have identified recombination as a rare event in HCV. However, our findings challenge this notion by providing the first evidence of 1a/3a (and vice versa) inter-genotypic recombination within the conserved IRES region. Utilizing advanced sequencing methods, such as deep sequencing and molecular cloning, our study reveals mixed infections involving genotypes 1a and 3a. This comprehensive approach not only confirmed the presence of mixed infections, but also identified the existence of recombinant forms not previously seen in the IRES region. The recombinant sequences, although present as low-frequency variants, open new avenues for understanding HCV evolution and adaptation.

A detailed molecular analysis of complete bovine leukemia virus genomes isolated from B-cell lymphosarcomas.

It is widely accepted that the majority of cancers result from multiple cellular events leading to malignancy after a prolonged period of clinical latency, and that the immune system plays a critical role in the control of cancer progression. Bovine leukemia virus (BLV) is an oncogenic member of the Retroviridae family. Complete genomic sequences of BLV strains isolated from peripheral blood mononuclear cells (PBMC) from cattle have been previously reported. However, a detailed characterization of the complete genome of BLV strains directly isolated from bovine tumors is much needed in order to contribute to the understanding of the mechanisms of leukemogenesis induced by BLV in cattle. In this study, we performed a molecular characterization of BLV complete genomes from bovine B-cell lymphosarcoma isolates. A nucleotide substitution was found in the glucocorticoid response element (GRE) site of the 5' long terminal repeat (5'LTR) of the BLV isolates. All amino acid substitutions in Tax previously found to be related to stimulate high transcriptional activity of 5'LTR were not found in these studies. Amino acid substitutions were found in the nucleocapsid, gp51 and G4 proteins. Premature stop-codons in R3 were observed. Few mutations or amino acid substitutions may be needed to allow BLV provirus to achieve silencing. Substitutions that favor suppression of viral expression in malignant B cells might be a strategy to circumvent effective immune attack.

Molecular characterisation of hepatitis B virus in the resident Chinese population in Panama City.

Despite the effectiveness of current hepatitis B virus (HBV) vaccines, it is estimated that 350 million individuals suffer from chronic HBV infection and more than 50% of these affected individuals live on the Asian continent. Panama is a country with a great diversity of foreign groups; the Chinese community is a large example of this phenomenon. There is an urgent need to perform studies that evaluate the prevalence and the genetic diversity of HBV in this community. This study aimed to evaluate the prevalence of HBV and its genotypes and mutant variants in the Chinese population residing in Panama. In total, 320 subjects were enrolled in the study. Forty-two subjects (13.1%) were positive for HBsAg and HBV-DNA from 18 subjects revealed the presence of genotypes B2 and C1. Secondary mutations associated with drug resistance at positions rtV207L and rtN239T of the reverse transcriptase gene were identified. Additionally, the mutation pair A1762T/G1764A was found in three samples and the mutation G1896A was detected in an HBeAg-negative subject. In conclusion, to our knowledge, this is the first study to report high HBV prevalence rates in resident ethnic Chinese in Central America and the presence of genotypes B2 and C1 in this region.

Spatial, temporal and evolutionary insight into seasonal epidemic Influenza A virus strains near the equatorial line: The case of Ecuador.

To study the spatial and temporal patterns of Influenza A virus (IAV) is essential for an efficient control of the disease caused by IAV and efficient vaccination programs. However, spatiotemporal patterns of spread as well as genetic lineage circulation of IAV on a countrywide scale have not been clearly determined for many tropical regions of the world. In order to gain insight into these matters, the spatial and temporal patterns of IAV in six different geographic regions of Ecuador, from 2011 to 2021, were determined and the timing and magnitude of IAV outbreaks in these localities investigated. The results of these studies revealed that although Ecuador is a South American country situated in the Equator line, its IAV epidemiology resembles that of temperate Northern Hemisphere countries. Phylogenetic analysis of H1N1pdm09 and H3N2 IAV strains isolated in five different localities of Ecuador revealed that provinces in the south of this country have the largest effective population size by comparison with provinces in the north, suggesting that the southern provinces may be acting as a source of IAV. Co-circulation of different H1N1pdm09 and H3N2 genetic lineages was observed in different geographic regions of Ecuador.

Phylogenetic analysis reveals that the H5N1 avian influenza A outbreak in poultry in Ecuador in November 2022 is associated with the highly pathogenic clade 2.3.4.4b.

The ongoing H5N1 outbreak in the Americas caused by clade 2.3.4.4 is causing unprecedented impact in poultry and wild birds. In November 2022, a highly pathogenic avian influenza A outbreak was declared in poultry in Ecuador, affecting more than 1.1 million heads of poultry in two farms by February 2023. Phylogenetic analysis shows that the virus clade is 2.3.4.4b, and to the best of our knowledge, this is the first scientific publication reporting this clade in South America.

Pandemic influenza A virus codon usage revisited: biases, adaptation and implications for vaccine strain development.

BACKGROUND: Influenza A virus (IAV) is a member of the family Orthomyxoviridae and contains eight segments of a single-stranded RNA genome with negative polarity. The first influenza pandemic of this century was declared in April of 2009, with the emergence of a novel H1N1 IAV strain (H1N1pdm) in Mexico and USA. Understanding the extent and causes of biases in codon usage is essential to the understanding of viral evolution. A comprehensive study to investigate the effect of selection pressure imposed by the human host on the codon usage of an emerging, pandemic IAV strain and the trends in viral codon usage involved over the pandemic time period is much needed. RESULTS: We performed a comprehensive codon usage analysis of 310 IAV strains from the pandemic of 2009. Highly biased codon usage for Ala, Arg, Pro, Thr and Ser were found. Codon usage is strongly influenced by underlying biases in base composition. When correspondence analysis (COA) on relative synonymous codon usage (RSCU) is applied, the distribution of IAV ORFs in the plane defined by the first two major dimensional factors showed that different strains are located at different places, suggesting that IAV codon usage also reflects an evolutionary process. CONCLUSIONS: A general association between codon usage bias, base composition and poor adaptation of the virus to the respective host tRNA pool, suggests that mutational pressure is the main force shaping H1N1 pdm IAV codon usage. A dynamic process is observed in the variation of codon usage of the strains enrolled in these studies. These results suggest a balance of mutational bias and natural selection, which allow the virus to explore and re-adapt its codon usage to different environments. Recoding of IAV taking into account codon bias, base composition and adaptation to host tRNA may provide important clues to develop new and appropriate vaccines.

Evolutionary history and spatiotemporal dynamic of GIII norovirus: From emergence to classification in four genotypes.

Noroviruses belong to a genetically diverse group of viruses infecting a wide range of mammalian host species, and those detected in cattle and sheep are classified within genogroup III (GIII). The current classification of norovirus in genogroups and genotypes is based on phylogenetic clustering and average distances within and between these phylogenetic clusters; however, the classification studies have been focused mainly on human norovirus, being GIII norovirus relegated. Due to the increasing number of studies on GIII norovirus, the need of an updated and extensive classification is evident. The aim of this study was to update the classification of norovirus within GIII, to describe the emergence of a circulating recombinant strain, and to reconstruct the evolutionary history of this genogroup. Two P-types (GIII.P1-2) and four genotypes (GIII.1-4) were described. For the genogroup GIII, the evolutionary rate estimated was 2.78E-3 s/s/y (95%HPD, 1.79E-3 s/s/y-3.78E-3 s/s/y), and the tMRCA was estimated around 1500 (95%HPD, 1247-1688). Despite the long history of this genogroup, the genotypes detected at present emerged in the last 100 years. Interestingly, most of the recombinant GIII.2P[1] strains detected worldwide were originated from a single recombination event and this recombinant strain was later dispersed through the world. Finally, our results indicate that a scenario of genotypes replacement through the time is highly probable.

An evolutionary insight into Severe Acute Respiratory Syndrome Coronavirus 2 Omicron variant of concern.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel virus that belongs to the family Coronaviridae. This virus produces a respiratory illness known as coronavirus disease 2019 (COVID-19) and is to blame for the pandemic of COVID-19. Due to its massive circulation around the world and the capacity of mutation of this virus, genomic studies are much needed in to order to reveal new variants of concern (VOCs). On November 26th, 2021, the WHO announced that a new SARS-CoV-2 VOC, named Omicron, had emerged. In order to get insight into the emergence, spread and evolution of Omicron SARS-CoV-2 variants, a comprehensive phylogenetic study was performed. The results of these studies revealed significant differences in codon usage among the S genes of SARS-CoV-2 VOCs Alfa, Beta, Gamma, Delta and Omicron, which can be linked to SARS-CoV-2 genotypes. Omicron variant did not evolve out of one of the early VOCs, but instead it belongs to a complete different genetic lineage from previous ones. Strains classified as Omicron variants evolved from ancestors that existed around May 15th, 2020, suggesting that this VOC may have been circulating undetected for a period of time until its emergence was observed in South Africa. A rate of evolution of 5.61 × 10-4 substitutions/site/year was found for Omicron strains enrolled in these analyses. The results of these studies demonstrate that S genes have suitable genetic information for clear assignment of emerging VOCs to its specific genotypes.

Phylogenetic analysis of SARS-CoV-2 viruses circulating in the South American region: Genetic relations and vaccine strain match.

The pandemic of coronavirus disease 2019 (COVID-19) is caused by a novel member of the family Coronaviridae, now known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Recent studies revealed the emergence of virus variants with substitutions in the spike and/or nucleocapsid and RNA-dependent RNA polymerase proteins that are partly responsible for enhanced transmission and reduced or escaped anti-SARS-CoV-2 antibodies that may reduce the efficacy of antibodies and vaccines against the first identified SARS-CoV-2 strains. In order to gain insight into the emergence and evolution of SARS-CoV-2 variants circulating in the South American region, a comprehensive phylogenetic study of SARS-CoV-2 variants circulating in this region was performed. The results of these studies revealed sharp increase in virus effective population size from March to April of 2020. At least 62 different genotypes were found to circulate in this region. Variants of concern (VOCs) Alpha, Beta, Gamma and Delta co-circulate in the region, together with variants of interest (VOIs) Lambda, Mu and Zeta. Most of SARS-CoV-2 variants circulating in the South American region belongs to B.1 genotypes and have substitutions in the spike and/or nucleocapsid and polymerase proteins that confer high transmissibility and/or immune resistance. 148 amino acid positions of the spike protein and 70 positions of the nucleocapsid were found to have substitutions in different variants isolated in the region by comparison with reference strain Wuhan-Hu-1. Significant differences in codon usage among spike genes of SARS-CoV-2 strains circulating in South America was found, which can be linked to SARS-CoV-2 genotypes.

Rotavirus A genotype P[4]G2: genetic diversity and reassortment events among strains circulating in Brazil between 2005 and 2009.

Group A rotaviruses (RV-A) are the leading cause of severe gastroenteritis in infants and young children worldwide. Due to the epidemiologic complexity of RV-A, especially in developing countries, it is important to determine which genotypes are circulating, principally after the introduction in March 2006 of the monovalent (P[8]G1) Rotarix® vaccine in Brazil by the National Immunization Program. In Phase III trials with Rotarix®, the prevalence of genotype P[4]G2 was extremely low, and therefore, evaluation of heterotypic immunization against this genotype was performed by meta-analysis statistics tests. Different studies have shown the re-emergence of genotype P[4]G2 in Brazil, since 2005, as well as in other countries, suggesting that it could be a continental phenomenon related to the temporal variability in the genotype's naturally occurring distribution. It is important to note that genotype P[4]G2 does not share VP4 or VP7 antigens with the vaccine strain. Therefore, we performed a phylogenetic analysis based on VP4 (VP8), VP7, VP6, and NSP4 genes of RV-A genotype P[4]G2 samples isolated from the five regions of Brazil between 2005 and 2009. This study revealed that different genetic variants of RV-A genotype P[4]G2 circulated in Brazil between 2005 and 2009, and that this variability is determined mainly by: occurrence of point mutations; reassortment events; and widespread global gene flow. The results obtained in this study are important to our understanding of the epidemiology and evolution of RV-A genotype P[4]G2 and demonstrate the importance of continuous monitoring and molecular characterization of RV-A strains circulating in human and animal populations.