Recombination detection and analysis using RDP3

Immunopathological investigation and genetic evolution of Avian leukosis virus Subgroup-J associated with myelocytomatosis in broiler flocks in Egypt

BACKGROUND

Avian leukosis virus Subgroup-J (ALV-J) is a rapidly oncogenic evolving retrovirus infecting a variety of avian species; causing severe economic losses to the local poultry industry.

METHODS

To investigate ALV-J, a total of 117 blood samples and 57 tissue specimens of different organs were collected for virological, and pathological identification, serological examinations, molecular characterization, and sequencing analysis. To the best of our knowledge, this is the first detailed report recorded in broiler flocks in Egypt. The present study targets the prevalence of a viral tumor disease circulating in broiler flocks in the El-Sharqia, El-Dakahliya, and Al-Qalyubiyya Egyptian governorates from 2021 to 2023 using different diagnostic techniques besides ALV-J gp85 genetic diversity determination.

RESULT

We first isolated ALV-J on chicken embryo rough cell culture; showing aggregation, rounding, and degeneration. Concerning egg inoculation, embryonic death, stunting, and curling were observed. Only 79 serum samples were positive for ALV-J (67.52%) based on the ELISA test. Histopathological investigation showed tumors consist of uniform masses, usually well-differentiated myelocytes, lymphoid cells, or both in the liver, spleen, and kidneys. Immunohistochemical examination showed that the myelocytomatosis-positive signals were in the spleen, liver, and kidney. The PCR assay of ALV-J gp85 confirmed 545 base pairs with only 43 positive samples (75.4%). Two positive samples were sequenced and submitted to the Genbank with accession numbers (OR509852-OR509853). Phylogenetic analysis based on the gp85 gene showed that the ALV-J Dakahlia-2 isolate is genetically related to ALV-EGY/YA 2021.3, ALV-EGY/YA 2021.4, ALV-EGY/YA 2021.14, and ALV-EGY/YA 2021.9 with amino acid identity percentage 96%, 97%; 96%, 96%; respectively. Furthermore, ALV-J Sharqia-1 isolate is highly genetically correlated to ALV-EGY/YA 2021.14, and ALV-EGY/YA 2021.9, ALV-J isolate QL1, ALV-J isolate QL4, ALV-J isolate QL3, ALV-EGY/YA 2021.4 with amino acid identity percentage 97%, 97%; 98%, 97%, 97%, 95%; respectively.

CONCLUSIONS

This study confirmed that ALV-J infection had still been prevalent in broilers in Egypt, and the genetic characteristics of the isolates are diverse.

Antigenic Characterization of Infectious Bronchitis Virus in the South China during 2021-2022

Avian infectious bronchitis is a serious and highly contagious disease that is caused by the infectious bronchitis virus (IBV). From January 2021 to June 2022, 1008 chicken tissue samples were collected from various regions of southern China, and 15 strains of the IBV were isolated. Phylogenetic analysis revealed that the strains mainly comprised the QX type, belonging to the same genotype as the currently prevalent LX4 type, and identified four recombination events in the S1 gene, among which lineages GI-13 and GI-19 were most frequently involved in recombination. Further study of seven selected isolates revealed that they caused respiratory symptoms, including coughing, sneezing, nasal discharge, and tracheal sounds, accompanied by depression. Inoculation of chicken embryos with the seven isolates resulted in symptoms such as curling, weakness, and bleeding. Immunization of specific pathogen-free (SPF) chickens with inactivated isolates produced high antibody levels that neutralized the corresponding strains; however, antibodies produced by vaccine strains were not effective in neutralizing the isolates. No unambiguous association was found between IBV genotypes and serotypes. In summary, a new trend in IBV prevalence has emerged in southern China, and currently available vaccines do not provide protection against the prevalent IBV strains in this region, facilitating the continued spread of IBV.

Insights into the Epidemiology, Phylodynamics, and Evolutionary Changes of Lineage GI-7 Infectious Bronchitis Virus

Infectious bronchitis virus (IBV) is distributed worldwide and causes significant losses in the poultry industry. In recent decades, lineages GI-19 and GI-7 have become the most prevalent IBV strains in China. However, the molecular evolution and phylodynamics of the lineage GI-7 IBV strains remain largely unknown. In this study, we identified 19 IBV strains from clinical samples from January 2021 to June 2022 in China, including 12 strains of GI-19, 3 strains of GI-7, and 1 strain each of GI-1, GI-9, GI-13, and GI-28. These results indicated that lineages GI-19 and GI-7 IBVs are still the most prevalent IBVs in China. Here, we investigated the evolution and transmission dynamics of lineage GI-7 IBVs. Our results revealed that the Taiwan province might be the origin of lineage GI-7 IBVs and that South China plays an important role in the spread of IBV. Furthermore, we found low codon usage bias of the S1 gene in lineage GI-7 IBVs. This allowed IBV to replicate in the host during evolution as a result of reduced competition, mainly driven by natural selection and mutational pressure, where the role of natural selection is more prominent. Collectively, our results reveal the genetic diversity and evolutionary dynamics of lineage GI-7 IBVs, which could assist in the prevention and control of viral infection.

Phylogenetic and Spatiotemporal Analyses of the Complete Genome Sequences of Avian Coronavirus Infectious Bronchitis Virus in China During 1985-2020: Revealing Coexistence of Multiple Transmission Chains and the Origin of LX4-Type Virus

Infectious bronchitis (IB) virus (IBV) causes considerable economic losses to poultry production. The data on transmission dynamics of IBV in China are limited. The complete genome sequences of 212 IBV isolates in China during 1985–2020 were analyzed as well as the characteristics of the phylogenetic tree, recombination events, dN/dS ratios, temporal dynamics, and phylogeographic relationships. The LX4 type (GI-19) was found to have the highest dN/dS ratios and has been the most dominant genotype since 1999, and the Taiwan-I type (GI-7) and New type (GVI-1) showed an increasing trend. A total of 59 recombinants were identified, multiple recombination events between the field and vaccine strains were found in 24 isolates, and the 4/91-type (GI-13) isolates were found to be more prone to being involved in the recombination. Bayesian phylogeographic analyses indicated that the Chinese IBVs originated from Liaoning province in the early 1900s. The LX4-type viruses were traced back to Liaoning province in the late 1950s and had multiple transmission routes in China and two major transmission routes in the world. Viral phylogeography identified three spread regions for IBVs (including LX4 type) in China: Northeastern China (Heilongjiang, Liaoning, and Jilin), north and central China (Beijing, Hebei, Shanxi, Shandong, and Jiangsu), and Southern China (Guangxi and Guangdong). Shandong has been the epidemiological center of IBVs (including LX4 type) in China. Overall, our study highlighted the reasons why the LX4-type viruses had become the dominant genotype and its origin and transmission routes, providing more targeted strategies for the prevention and control of IB in China.

Experimental evolution of cowpea mild mottle virus reveals recombination-driven reduction in virulence accompanied by increases in diversity and viral fitness

Major themes in pathogen evolution are emergence, evolution of virulence, host adaptation and the processes that underlie them. RNA viruses are of particular interest due to their rapid evolution. The in vivo molecular evolution of an RNA plant virus was demonstrated here using a necrotic isolate of cowpea mild mottle virus (CPMMV) and a susceptible soybean genotype submitted to serial inoculations. We show that the virus lost the capacity to cause necrosis after six passages through the host plant. When a severe bottleneck was imposed, virulence reduction occurred in the second passage. The change to milder symptoms had fitness benefits for the virus (higher RNA accumulation) and for its vector, the whitefly Bemisia tabaci. Genetic polymorphisms were highest in ORF1 (viral replicase) and were independent of the symptom pattern. Recombination was a major contributor to this diversity - even with the strong genetic bottleneck, recombination events and hot spots were detected within ORF1. Virulence reduction was associated with different sites in ORF1 associated to recombination events in both experiments. Overall, the results demonstrate that the reduction in virulence was a consequence of the emergence of new variants, driven by recombination. Besides providing details of the evolutionary mechanisms behind a reduction in virulence and its effect under viral and vector fitness, we propose that this recombination-driven switch in virulence allows the pathogen to rapidly adapt to a new host and, potentially, switch back.

Analysis of codon usage bias in potato virus Y non-recombinant strains

Potato virus Y (PVY) is a member of the genus Potyvirus, family Potyviridae, is considered one of the most devastating pest affecting economically important crops, such as potato, tobacco, tomato and pepper, representing a serious threat due to high incidence and worldwide distribution. Its economic significance as well as it biological and molecular complexities have aroused great attention, thus several studies have explore it genetic characteristics. However, little is known about PVY codon usage. To shed light on the relation of codon usage among viruses and their hosts is extremely important to understand virus survival, fitness and evolution. In this study, we performed a comprehensive analysis of codon usage and composition of PVY non-recombinant strains (PVY^N-NA, PVY^Eu-N, PVY^O, PVY^O5, PVY^C) based on 130 complete open reading frame sequences extracted from public databases. Furthermore, similarities between the synonymous codon usage of PVY and its main hosts were investigated. The results obtained in the current study suggest that the overall codon usage among PVY genotypes is similar and slightly biased. PVY codon usage is strongly influenced by mutational bias, but also by G + C compositional constraint and dinucleotide composition. Furthermore, similarities among codon usage preferences between PVY strains and analyzed hosts were observed.

Fowl Adenovirus (FAdV) Recombination with Intertypic Crossovers in Genomes of FAdV-D and FAdV-E, Displaying Hybrid Serological Phenotypes

After analyzing 27 new genomes from fowl adenovirus (FAdV) field isolates and so-far unsequenced prototypes, we report the first evidence for recombination in FAdVs. Recombination was confined to species FAdV-D and FAdV-E, accommodating the largest number of, and the intraspecies-wise most differentiated, types. The majority of detected events occurred in FAdV-E, involving segments with parental origin of all constitutive types. Together with the diversity of breakpoints, this suggests widespread recombination in this species. With possible constraints through species-specific genes and diversification patterns, the recombinogenic potential of FAdVs attains particular interest for inclusion body hepatitis (IBH), an important disease in chickens, caused by types from the recombination-prone species. Autonomously evolving, recombinant segments were associated with major sites under positive selection, among them the capsid protein hexon and fiber genes, the right-terminal ORFs 19, 25, and the ORF20/20A family. The observed mosaicism in genes indicated as targets of adaptive pressures points toward an immune evasion strategy. Intertypic hexon/fiber-recombinants demonstrated hybrid neutralization profiles, retrospectively explaining reported controversies on reference strains B3-A, T8-A, and X11-A. Furthermore, cross-neutralization supported sequence-based evidence for interdomain recombination in fiber and contributed to a tentatively new type. Overall, our findings challenge the purported uniformity of types responsible for IBH, urging more complete identification strategies for FAdVs. Finally, important consequences arise for in vivo studies investigating cross-protection against IBH.

The evolutionary dynamics of DENV 4 genotype I over a 60-year period

Dengue virus serotype 4 (DENV 4) has had a relatively low prevalence worldwide for decades; however, likely due to data paucity, no study has investigated the epidemiology and evolutionary dynamics of DENV 4 genotype I (DENV 4-I). This study aims to understand the diversity, epidemiology and dynamics of DENV 4-I. We collected 404 full length DENV4-1 envelope (E) gene sequences from 14 countries using two sources: Yunnan Province in China (15 strains during 2013–2016) and GenBank (489 strains up to 2018-01-11). Conducting phylogenetic and phylogeographical analyses, we estimated the virus spread, population dynamics, and selection pressures using different statistical analysis methods (substitution saturation, likelihood mapping, Bayesian coalescent inference, and maximum likelihood estimation). Our results show that during the last 60 years (1956–2016), DENV 4-I was present in mainland and maritime Southeast Asia, the Indian subcontinent, the southern provinces of China, parts of Brazil and Australia. The recent spread of DENV 4-I likely originated in the Philippines and later spread to Thailand. From Thailand, it spread to adjacent countries and eventually the Indian subcontinent. Apparently diverging around years 1957, 1963, 1976 and 1990, the different Clades (Clade I-V) were defined. The mean overall evolution rate of DENV 4-I was 9.74 (95% HPD: 8.68–10.82) × 10⁻⁴ nucleotide substitutions/site/year. The most recent common ancestor for DENV 4-I traces back to 1956. While the demographic history of DENV 4-I fluctuated, peaks appeared around 1982 and 2006. While purifying selection dominated the majority of E-gene evolution of DENV 4-I, positive selection characterized Clade III (Vietnam). DENV 4-I evolved in situ in Southeast Asia and the Indian subcontinent. Thailand and Indian acted as the main and secondary virus distribution hubs globally and regionally. Our phylogenetic analysis highlights the need for strengthened regional cooperation on surveillance and sharing of sample sequences to improve global dengue control and cross-border transmission prevention efforts.

Dengue virus (DENV) can be classified into four serotypes, DENV 1, 2, 3 and 4. Although DENV 4 is the first dengue serotype to diverge in phylogenetic analyses of the genus Flavivirus, this serotype occurs at a low prevalence worldwide and spreads the least rapidly. Similar to other serotypes, DENV 4 can also cause severe dengue (SD) disease manifestations, such as dengue haemorrhagic fever (DHF) and dengue shock syndrome (DSS). To date, no study has investigated the epidemiology and dynamics of DENV 4 genotype I comprehensively. In this study, we seek to address this gap. Our study shows that the distribution of DENV 4-I is mainly restricted to Southeast Asia and the Indian subcontinent. The most recent spread of DENV 4-I likely originated from Southeast Asia–initially circulating in the Philippines, then Thailand and later on the Indian subcontinent. Viruses evolved in situ in Southeast Asia and the Indian subcontinent, respectively. Although DENV 4-I occasionally spread elsewhere, this genotype did not become widely established. The overall evolution rate of DENV 4-I was comparable with that of DENV 2–4. The nucleotide sequences indicates that the demographic history of DENV 4-I fluctuated with peaks apparent during parts of the 1980s and 2000s. Although a weak positive selection existed in Clade III -predominately in Vietnam, purifying selection dominated the E-gene evolution of DENV 4-I.

Molecular characteristics of avian leukosis viruses isolated from indigenous chicken breeds in China

To assess the status of avian leukosis virus (ALV) infection in indigenous chicken breeds in China, 121 plasma samples collected from various indigenous chicken breeds were tested for the presence of ALV from 2015 to 2016. A total of 14 ALV strains were isolated and identified, including two ALV-A strains, one ALV-B strain, eight ALV-J strains, and three ALV-K strains. To study the genome structure, biological characteristics, and the evolutionary relationships of the ALV-K strains with other known subgroup strains from infected chickens, we determined the complete genome sequence of the three ALV-K strains and performed comparative analysis using the whole genome sequence or selected sequence elements. The replication rates of the three ALV-K strains were markedly lower than the rates of other ALVs, and they shared a common mutation in the pol gene, which had not been previously observed. In addition, nine putative recombinant events were detected in the genomes of the three newly isolated ALV-K strains, with high statistical support. This was the first report of an ALV-K reorganization event, which has contributed to its genetic evolution. In summary, we established a robust classification system for ALV, especially for ALV-K, and revealed additional genomic diversity for the ALV strains in indigenous chicken breeds. Therefore additional works are warranted to explore ALV genomics and epidemiology.

Molecular characterization of chicken infectious anemia virus from contaminated live-virus vaccines

The aim of this study was to investigate possible causes of the pervasiveness of chicken infectious anemia virus (CIAV) infection in chickens in recent years in China. A total of 14 batches of live-virus vaccines were examined using PCR to detect CIAV contamination, of which only 2 samples (a Newcastle disease vaccine and a fowl pox vaccine) tested positive for CIAV. These Newcastle and fowl pox vaccines were then inoculated into 1-day-old specific-pathogen-free chickens. Serum samples were collected from chickens infected with the PCR-positive vaccines, and these tested positive for CIAV-specific antibodies as tested using ELISA. In addition, DNA samples isolated from the serum samples also tested positive by PCR. The results indicated that the samples were contaminated with CIAV and identified 2 exogenous CIAV strains, designated CIAV-N22 and CIAV-F10, in the respective samples. The full genome sequences of these novel CIAV strains were sequenced and analyzed. Phylogenetic tree analysis indicated that the CIAV-F10 strain might represent a recombinant viral strain arising from the parental CIAV strains JQ690762 and KJ728816. Overall, the results suggested that vaccination with CIAV-contaminated vaccines contributed to the prevalence and spread of CIAV infection in chickens. Furthermore, the CIAV contaminant was likely subsequently transmitted to commercial chickens through congenital transmission. Our findings therefore highlight the need for more extensive screening of live-virus vaccines for poultry in China to reduce the threat of contamination with exogenous viruses.

GET_PHYLOMARKERS, a Software Package to Select Optimal Orthologous Clusters for Phylogenomics and Inferring Pan-Genome Phylogenies, Used for a Critical Geno-Taxonomic Revision of the Genus Stenotrophomonas

The massive accumulation of genome-sequences in public databases promoted the proliferation of genome-level phylogenetic analyses in many areas of biological research. However, due to diverse evolutionary and genetic processes, many loci have undesirable properties for phylogenetic reconstruction. These, if undetected, can result in erroneous or biased estimates, particularly when estimating species trees from concatenated datasets. To deal with these problems, we developed GET_PHYLOMARKERS, a pipeline designed to identify high-quality markers to estimate robust genome phylogenies from the orthologous clusters, or the pan-genome matrix (PGM), computed by GET_HOMOLOGUES. In the first context, a set of sequential filters are applied to exclude recombinant alignments and those producing anomalous or poorly resolved trees. Multiple sequence alignments and maximum likelihood (ML) phylogenies are computed in parallel on multi-core computers. A ML species tree is estimated from the concatenated set of top-ranking alignments at the DNA or protein levels, using either FastTree or IQ-TREE (IQT). The latter is used by default due to its superior performance revealed in an extensive benchmark analysis. In addition, parsimony and ML phylogenies can be estimated from the PGM. We demonstrate the practical utility of the software by analyzing 170 Stenotrophomonas genome sequences available in RefSeq and 10 new complete genomes of Mexican environmental S. maltophilia complex (Smc) isolates reported herein. A combination of core-genome and PGM analyses was used to revise the molecular systematics of the genus. An unsupervised learning approach that uses a goodness of clustering statistic identified 20 groups within the Smc at a core-genome average nucleotide identity (cgANIb) of 95.9% that are perfectly consistent with strongly supported clades on the core- and pan-genome trees. In addition, we identified 16 misclassified RefSeq genome sequences, 14 of them labeled as S. maltophilia, demonstrating the broad utility of the software for phylogenomics and geno-taxonomic studies. The code, a detailed manual and tutorials are freely available for Linux/UNIX servers under the GNU GPLv3 license at https://github.com/vinuesa/get_phylomarkers. A docker image bundling GET_PHYLOMARKERS with GET_HOMOLOGUES is available at https://hub.docker.com/r/csicunam/get_homologues/, which can be easily run on any platform.

Characterization of tomato leaf curl purple vein virus, a new monopartite New World begomovirus infecting tomato in Northeast Brazil

A new begomovirus species was identified from tomato plants with upward leaf curling and purple vein symptoms, which was first identified in the Piaui state of Northeast (NE) Brazil in 2014. Tomato leaf samples were collected in 2014 and 2016, and PCR with degenerate primers revealed begomovirus infection. Rolling circle amplification and restriction enzyme digestion indicated a single genomic DNA of ~ 2.6 kb. Cloning and sequencing revealed a genome organization similar to DNA-A components of New World (NW) bipartite begomoviruses, with no DNA-B. The complete nucleotide sequence had the highest identity (80%) with the DNA-A of Macroptilium yellow spot virus (MacYSV), and phylogenetic analyses showed it is a NW begomovirus that clusters with MacYSV and Blainvillea yellow spot virus, also from NE Brazil. Tomato plants agroinoculated with a dimeric clone of this genomic DNA developed upward leaf curling and purple vein symptoms, indistinguishable from those observed in the field. Based on agroinoculation, this virus has a narrow host range, mainly within the family Solanaceae. Co-inoculation experiments with tomato severe rugose virus and tomato mottle leaf curl virus, the two predominant begomoviruses infecting tomato in Brazil, revealed a synergistic interaction among these begomoviruses. The name Tomato leaf curl purple vein virus (ToLCPVV) is proposed for this new begomovirus.

Genomic Characterization of Recent Chicken Anemia Virus Isolates in China

Chicken anemia virus (CAV) causes diseases in young chickens, which include increased pathogenicity of secondary infectious agents, generalized lymphoid depletion, and immunodepression. In the present study, we have identified 22 CAV strains isolated from several commercial chicken farms in Northern China during 2014–2015. In addition, two CAVs were also isolated from stray mouse and dog feces, respectively. To our knowledge, this is the first report of identification of CAV from mouse and dog feces. Phylogenetic analysis of 121 full-length CAV genome sequences showed that all available CAV could be classified into eight lineages, supported by phylogenetic trees estimated using different methods. Furthermore, the 24 novel CAV sequences scattered across different branches, lack of clear spatio-temporal distribution characterization. Analysis of the 450 amino acids of VP1 protein identified 33 amino acid substitutions that were specific for CAVs from northern China. Putative gene recombination events were also detected in the genomes of newly isolated CAVs. In particular, a putative recombinant event was detected in the CAV-Dog genome with high statistical support. In summary, we established a robust classification system for CAV, revealed additional genomic diversity of CAV, and therefore, warranted additional efforts to explore CAV genomics and epidemiology.

[Discovery of DNA viruses in wildlife in Zambia and Indonesia]

Zoonoses originate from pathogens harbored in domestic and wild animals and therefore it is likely impossible to completely eradicate zoonotic diseases. For pre-emptive measures to attempt to predict the emergence of zoonosis outbreaks and the prevention of future epidemics and pandemics, it is imperative to identify natural host animals carrying potential pathogens and elucidate the routes of pathogen transmission into the human population. Our research team is conducting epidemiological research studies in Zambia and Indonesia for the control of viral zoonotic diseases. In this review, we present the research findings, including the discovery of orthopoxviruses and polyomaviruses in wildlife in Zambia and the identification of herpesviruses in bats in Indonesia among our activities.

Molecular evolution of emerging Banna virus

Banna virus (BAV) is an emerging pathogen that causes human viral encephalitis and has been isolated from types of blood-sucking insects and mammals in Asia. However, there are no reported systematic studies that describe the origin and evolution of BAV. Here, a phylogenetic analysis of BAVs isolated from a variety of potential vectors and vertebrate hosts worldwide revealed that BAVs emerged in the beginning of the 20th century and do not exhibit a species barrier. The mean substitution rate of BAVs was 2.467×10^-2substitution/site/year (95% HPD, 1.093×10^-3 to 5.628×10^-2). The lineage is mainly composed of BAVs from high-latitude regions, which are the most recently emerged viruses with significantly higher substitution rates compared with the lineage comprised of the isolates from middle or low-latitude regions. The genetic differences between BAV strains are positively correlated with the geographic distribution. Strains from the same latitude regions are almost 100% identical, whereas the differences between strains from long distance regions with different latitudes could be >60%. Our results demonstrate that BAV is an emerging virus at a stage that involves rapid evolution and has great potential for introduction into non-endemic areas. Thus, enhanced surveillance of BAV is highly recommended worldwide.

Characterization of VP1 sequence of Coxsackievirus A16 isolates by Bayesian evolutionary method

BACKGROUND

Coxsackievirus A16 (CV-A16), a major etiopathologic cause of pediatric hand, foot, and mouth disease (HFMD) worldwide, has been reported to have caused several fatalities. Revealing the evolutionary and epidemiologic dynamics of CV-A16 across time and space is central to understanding its outbreak potential.

METHODS

In this study, we isolated six CV-A16 strains in China's Jilin province and construct a maximum clade credibility (MCC) tree for CV-A16 VP1 gene by the Bayesian Markov Chain Monte Carlo method using 708 strains from GenBank with epidemiological information. The evolution characteristics of CV-A16 VP1 gene was also analysed dynamicly through Bayesian skyline plot.

RESULTS

All CV-A16 strains identified could be classified into five major genogroups, denoted by GI-GV. GIV and GV have co-circulated in China since 2007, and the CV-A16 epidemic strain isolated in the Jilin province, China, can be classified as GIV-3. The CV-A16 genogroups circulating recently in China have the same ancestor since 2007. The genetic diversity of the CV-A16 VP1 gene shows a continuous increase since the mid-1990s, with sharp increases in genetic diversity in 1997 and 2007 and reached peak in 2007. Very low genetic diversity existed after 2010. The CV-A16 VP1 gene evolutionary rate was 6.656E-3 substitutions per site per year.

CONCLUSIONS

We predicted the dynamic phylogenetic trends, which indicate outbreak trends of CV-A16, and provide theoretical foundations for clinical prevention and treatment of HFMD which caused by a CV-A16.

Complete genome analysis of a novel recombinant isolate of pepper veinal mottle virus from mainland China

BACKGROUND

Pepper veinal mottle virus (PVMV) was well established in Africa, and also reported infecting pepper (Capsicum annuum L) in Taiwan and India. However, there is not available of PVMV in mainland China. Here, the first complete genome sequence of PVMV isolated from pepper in mainland China was reported.

FINDING

The complete genomic sequence of isolate PVMV-HN isolated from pepper in mainland China is reported in this study. The genome of PVMV-HN is 9793 nucleotides (nt) excluding the poly (A) tail, shares 98-99 % nucleotide sequence identity with those two PVMV isolates from Ghana and Taiwan. Recombinant analysis showed that PVMV-HN probably represents a novel recombinant of PVMV. The phylogenetic relationship of PVMV-HN isolate to other PVMV isolates and other potyviruses based on genome or polyprotein sequence level and CP gene level, was also analyzed in this study.

CONCLUSION

The current study will help to understand phylogenetic relationship of isolate PVMV-HN.

Enzootic transmission of yellow fever virus, Venezuela

Phylogenetic analysis of yellow fever virus (YFV) strains isolated from Venezuela strongly supports YFV maintenance in situ in Venezuela, with evidence of regionally independent evolution within the country. However, there is considerable YFV movement from Brazil to Venezuela and between Trinidad and Venezuela.

Evidence of recombinant strains of porcine epidemic diarrhea virus, United States, 2013

To investigate the evolutionary process by which porcine epidemic diarrhea virus (PEDV) in the United States hypothetically descended from strains in China, we analyzed PEDV-positive samples collected in China during January 2012–July 2013. Recombination in 2 strain sublineages was likely associated with identification of PEDV in the United States in 2013.

Dengue is still an imported disease in China: a case study in Guangzhou

Dengue virus and its four serotypes (DENV 1-4) infect approximately 390 million people worldwide each year, with most cases in tropical and subtropical regions. Because of repeated introduction of DENV from epidemic regions and suitable weather conditions, many regions have shifted from hypo-endemicity to hyper-endemicity over recent decades. Since the first dengue outbreak in 1978, it is crucial to understand the current situation in China over nearly 40 years. The purpose of the study was to examine whether dengue in China was endemic or not, which is essential for relevant dengue control and prevention strategy implementation in China. The study, combining epidemiological characteristics of dengue from the disease notification system, phylogenetic and phylogeographic analyses, showed that all four serotypes had been detected in Guangzhou, China, which was dominated by DENV 1-2. The Maximum Likelihood tree analytic results showed that the virus detected in Guangzhou localized in different clades, except of virus of 2002 and 2003 clustered together. There existed the mutual introductions between Guangzhou and Southeast Asia. Most of the viruses were imported from Southeast Asia and the sources of outbreaks in Guangzhou mainly originated from Thailand, Indonesia, and the Philippines. The study indicates that dengue in China still remains as an imported disease, with the possibility of localization.

Comparison of evolutionary patterns of different HBeAg status in HBV genotype B and C

The aim of this work was to investigate the different evolutionary patterns between HBeAg(+) and HBeAg(-) status of genotype B and genotype C of Hepatitis B virus (HBV). Complete HBV genome sequences were obtained from GenBank, and were classified as HBeAg(+) of genotype B, HBeAg(-) of genotype B, HBeAg(+) of genotype C, and HBeAg(-) of genotype C. The results showed that genotype B had lower diversity and a longer phylogenetic tree than genotype C. This means that genotype B was more conserved than genotype C. The genetic distance of HBeAg(+) was shorter than that of HBeAg(-) in genotype B. dN, dS, dN/dS of HBeAg(+) were lower than that of HBeAg(-). In addition, the Core and S genes had lower diversity, and were under lower selective pressure than other genes. We propose that by having different HBeAg HBV receives disparate immune selective pressures, however these different genotypes also receive disparate natural selective pressures. The evolutionary pattern of HBV could be one of the important factors that may affect viral replication in HBV infected patients and could therefore play an important role in the diseases progression.

Molecular and Biological Characterization of Distinct Strains of Jatropha mosaic virus from the Dominican Republic Reveal a Potential to Infect Crop Plants

In the Dominican Republic (DO), jatropha plants with yellow mosaic symptoms are commonly observed in and around fields of various crop plants. Complete nucleotide sequences of DNA-A and DNA-B components of four bipartite begomovirus isolates associated with symptomatic jatropha plants collected from three geographical locations in the DO were determined. Sequence comparisons revealed highest identities (91 to 92%) with the DNA-A component of an isolate of Jatropha mosaic virus (JMV) from Jamaica, indicating that the bipartite begomovirus isolates from the DO are strains of JMV. When introduced into jatropha seedlings by particle bombardment, the cloned components of the JMV strains from the DO induced stunting and yellow mosaic, indistinguishable from symptoms observed in the field, thereby fulfilling Koch's postulates for the disease. The JMV strains also induced disease symptoms in Nicotiana benthamiana, tobacco, and several cultivars of common bean from the Andean gene pool, including one locally grown in the DO. Asymmetry in the infectivity and symptomatology of pseudorecombinants provided further support for the strain designation of the JMV isolates from the DO. Thus, JMV in the DO is a complex of genetically distinct strains that have undergone local evolution and have the potential to cause disease in crop plants.

HGT turbulence: Confounding phylogenetic influence of duplicative horizontal transfer and differential gene conversion

Horizontal gene transfer (HGT) often leads to phylogenetic incongruence. When "duplicative HGT" introduces a second copy of a pre-existing gene, the two copies may then engage in gene conversion, leading to phylogenetically mosiac genes. When duplicative HGT is followed by differential gene conversion among descendant lineages, as under the DH-DC model, phylogenetic analysis is further complicated. To explore the effects of DH-DC on phylogeny reconstruction, we analyzed two sets of sequences: (1) an augmented set of plant mitochondrial atp1 sequences for which we recently published evidence of DH-DC; and (2) a set of simulated sequences for which we varied the extent of chimerism, the number of chimeric genes and nucleotide substitution rates. We show that the phylogenetic behavior of evolutionarily chimeric genes is highly volatile and depends on both the degree of chimerism and the number of differentially chimeric genes present in the analysis. Furthermore, we show that the presence of chimeric genes in gene trees can spuriously affect the phylogenetic position of purely native sequences, especially by attracting these sequences toward basal positions in trees. We propose the term "HGT turbulence" to describe these complex effects of evolutionarily chimeric genes on phylogenetic results.

Genomic heterogeneity and ecological speciation within one subspecies of Bacillus subtilis

Closely related bacterial genomes usually differ in gene content, suggesting that nearly every strain in nature may be ecologically unique. We have tested this hypothesis by sequencing the genomes of extremely close relatives within a recognized taxon and analyzing the genomes for evidence of ecological distinctness. We compared the genomes of four Death Valley isolates plus the laboratory strain W23, all previously classified as Bacillus subtilis subsp. spizizenii and hypothesized through multilocus analysis to be members of the same ecotype (an ecologically homogeneous population), named putative ecotype 15 (PE15). These strains showed a history of positive selection on amino acid sequences in 38 genes. Each of the strains was under a different regimen of positive selection, suggesting that each strain is ecologically unique and represents a distinct ecological speciation event. The rate of speciation appears to be much faster than can be resolved with multilocus sequencing. Each PE15 strain contained unique genes known to confer a function for bacteria. Remarkably, no unique gene conferred a metabolic system or subsystem function that was not already present in all the PE15 strains sampled. Thus, the origin of ecotypes within this clade shows no evidence of qualitative divergence in the set of resources utilized. Ecotype formation within this clade is consistent with the nanoniche model of bacterial speciation, in which ecotypes use the same set of resources but in different proportions, and genetic cohesion extends beyond a single ecotype to the set of ecotypes utilizing the same resources.

Phylogenetic and recombination analysis of the homing protein domain of grapevine fanleaf virus (GFLV) isolates associated with 'yellow mosaic' and 'infectious malformation' syndromes in grapevine

The RNA2 of seven grapevine fanleaf virus (GFLV) isolates from vines with yellow mosaic (YM) symptoms from different origin were sequenced. These sequences showed a high variability in the homing protein (2A(HP)) and, in five of them, a putative recombination with arabis mosaic virus (ArMV) was detected. To investigate recombination frequency, the partial sequences of the 2A(HP) of 28 additional GFLV isolates from nine different countries, showing either YM or infectious malformations (MF) symptoms, were obtained and compared with those of GFLV isolates from GenBank. The analysis confirmed the high level of sequence variability (up to 41 % at the nucleotide level) among isolates. In phylogenetic trees constructed using different approaches, the sequenced isolates always clustered in four conserved groups, three of which comprised YM strains (groups 1, 2 and 3), and one (group 4) the MF strains. Potential interspecific recombination sites between GFLV and ArMV were predicted in the 2A(HP) gene of several isolates, all of which were associated with YM symptoms.

Identification of a chicken anemia virus variant-related gyrovirus in stray cats in china, 2012

The chicken anemia virus (CAV), is a known member of the genus Gyrovirus and was first isolated from chickens in Japan in 1979. Some reports have also demonstrated that CAV can be identified in human stool specimens. In this study, a variant of CAV was detected using PCR with CAV-based primers in fecal samples of stray cats. The genome of CAV variant was sequenced and the results suggest that it could be a recombinant viral strain from parental CAV strains JQ690762 and AF311900. Recombination is an important evolutionary mechanism that contributes to genetic diversification. These findings indicate that CAV variant might have originated from CAV-infected chickens. The epidemiology and pathogenesis of this novel virus remains to be elucidated. This study underscores the importance of CAV surveillance and it presents the first evidence suggesting the possibility of CAV homologous recombination in cat.

Phylogenetic and molecular characterization of chicken anemia virus in southern China from 2011 to 2012

Chicken anemia virus (CAV) is an important pathogen that causes severe immunosuppression in young chickens. We have characterized 13 CAVs isolated from different commercial farms in southern China between 2011 and 2012. We discovered 92 variable residues compared to 37 other CAV complete genome sequences from other parts of the world listed in GenBank; these residues have not been previously observed. All of the Chinese CAV genomes that were characterized in this study had a glutamine at position 394, a hallmark of highly pathogenic CAVs. We also discovered that intra-group genetic recombination plays a role in generating genetic diversity in natural populations of CAV. The GD-J-12 isolate was a possible recombinant between GD-C-12 and GD-M-12 in the genomic region that encompassed both the coding and non-coding regions.

New avian paramyxoviruses type I strains identified in Africa provide new outcomes for phylogeny reconstruction and genotype classification

Newcastle disease (ND) is one of the most lethal diseases of poultry worldwide. It is caused by an avian paramyxovirus 1 that has high genomic diversity. In the framework of an international surveillance program launched in 2007, several thousand samples from domestic and wild birds in Africa were collected and analyzed. ND viruses (NDV) were detected and isolated in apparently healthy fowls and wild birds. However, two thirds of the isolates collected in this study were classified as virulent strains of NDV based on the molecular analysis of the fusion protein and experimental in vivo challenges with two representative isolates. Phylogenetic analysis based on the F and HN genes showed that isolates recovered from poultry in Mali and Ethiopia form new groups, herein proposed as genotypes XIV and sub-genotype VIf with reference to the new nomenclature described by Diel's group. In Madagascar, the circulation of NDV strains of genotype XI, originally reported elsewhere, is also confirmed. Full genome sequencing of five African isolates was generated and an extensive phylogeny reconstruction was carried out based on the nucleotide sequences. The evolutionary distances between groups and the specific amino acid signatures of each cluster allowed us to refine the genotype nomenclature.

Southernmost Asia is the source of Japanese encephalitis virus (genotype 1) diversity from which the viruses disperse and evolve throughout Asia

BACKGROUND

Although a previous study predicted that Japanese encephalitis virus (JEV) originated in the Malaysia/Indonesia region, the virus is known to circulate mainly on the Asian continent. However, there are no reported systematic studies that adequately define how JEV then dispersed throughout Asia.

METHODOLOGY/PRINCIPAL FINDINGS

In order to understand the mode of JEV dispersal throughout the entire Asian continent and the factors that determine the dispersal characteristics of JEV, a phylogenetic analysis using Bayesian Markov chain Monte Carlo simulations was conducted on all available JEV E gene sequences in GenBank, plus strains recently isolated in China. Here we demonstrate for the first time that JEV lineages can be divided into four endemic cycles, comprising southern Asia, eastern coastal Asia, western Asia, and central Asia. The isolation places of the viruses in each endemic cycle were geographically independent regardless of years, vectors, and hosts of isolation. Following further analysis, we propose that the southernmost region (Thailand, Vietnam, and Yunnan Province, China) was the source of JEV transmission to the Asian continent following its emergence. Three independent transmission routes from the south to north appear to define subsequent dispersal of JEV. Analysis of JEV population dynamics further supports these concepts.

CONCLUSIONS/SIGNIFICANCE

These results and their interpretation provide new insights into our understanding of JEV evolution and dispersal and highlight its potential for introduction into non-endemic areas.

Conservation of meningococcal antigens in the genus Neisseria

Neisseria meningitidis, one of the major causes of bacterial meningitis and sepsis, is a member of the genus Neisseria, which includes species that colonize the mucosae of many animals. Three meningococcal proteins, factor H-binding protein (fHbp), neisserial heparin-binding antigen (NHBA), and N. meningitidis adhesin A (NadA), have been described as antigens protective against N. meningitidis of serogroup B, and they have been employed as vaccine components in preclinical and clinical studies. In the vaccine formulation, fHbp and NHBA were fused to the GNA2091 and GNA1030 proteins, respectively, to enhance protein stability and immunogenicity. To determine the possible impact of vaccination on commensal neisseriae, we determined the presence, distribution, and conservation of these antigens in the available genome sequences of the genus Neisseria, finding that fHbp, NHBA, and NadA were conserved only in species colonizing humans, while GNA1030 and GNA2091 were conserved in many human and nonhuman neisseriae. Sequence analysis showed that homologous recombination contributed to shape the evolution and distribution of both NHBA and fHbp, three major variants of which have been defined. fHbp variant 3 was probably the ancestral form of meningococcal fHbp, while fHbp variant 1 from N. cinerea was introduced into N. meningitidis by a recombination event. fHbp variant 2 was the result of a recombination event inserting a stretch of 483 bp from variant 1 into the variant 3 background. These data indicate that a high rate of exchange of genetic material between neisseriae that colonize the human upper respiratory tract exists.

The upper respiratory tract of healthy individuals is a complex ecosystem colonized by many bacterial species. Among these, there are representatives of the genus Neisseria, including Neisseria meningitidis, a major cause of bacterial meningitis and sepsis. Given the close relationship between commensal and pathogenic species, a protein-based vaccine against N. meningitidis has the potential to impact the other commensal species of Neisseria. For this reason, we have studied the distribution and evolutionary history of the antigen components of a recombinant vaccine, 4CMenB, that recently received approval in Europe under the commercial name of Bexsero®. We found that fHbp, NHBA, and NadA can be found in some of the human commensal species and that the evolution of these antigens has been essentially shaped by the high rate of genetic exchange that occurs between strains of neisseriae that cocolonize the same environment.

Characterization of a new world monopartite begomovirus causing leaf curl disease of tomato in Ecuador and Peru reveals a new direction in geminivirus evolution

All characterized whitefly-transmitted geminiviruses (begomoviruses) with origins in the New World (NW) have bipartite genomes composed of a DNA-A and DNA-B component. Recently, an NW begomovirus lacking a DNA-B component was associated with tomato leaf curl disease (ToLCD) in Peru, and it was named Tomato leaf deformation virus (ToLDeV). Here, we show that isolates of ToLDeV associated with ToLCD in Ecuador and Peru have a single, genetically diverse genomic DNA that is most closely related to DNA-A components of NW bipartite begomoviruses. Agroinoculation of multimeric clones of the genomic DNA of three ToLDeV genotypes (two variants and a strain) resulted in the development of tomato leaf curl symptoms indistinguishable from those of ToLCD in Ecuador and Peru. Biological properties of these ToLDeV genotypes were similar to those of Old World (OW) monopartite tomato-infecting begomoviruses, including lack of sap transmissibility, phloem limitation, a resistance phenotype in tomato germplasm with the Ty-1 gene, and functional properties of the V1 (capsid protein) and C4 genes. Differences in symptom phenotypes induced by the ToLDeV genotypes in tomato and Nicotiana benthamiana plants were associated with a highly divergent left intergenic region and C4 gene. Together, these results establish that ToLDeV is an emergent NW monopartite begomovirus that is causing ToLCD in Ecuador and Peru. This is the first report of an indigenous NW monopartite begomovirus, and evidence is presented that it emerged from the DNA-A component of a NW bipartite progenitor via convergent evolution and recombination.

Highly divergent virulent isolates of Newcastle disease virus from the Dominican Republic are members of a new genotype that may have evolved unnoticed for over 2 decades

A Newcastle disease virus (NDV) outbreak in chickens was reported in the Dominican Republic in 2008. The complete genome of this isolate, chicken/DominicanRepublic(JuanLopez)/499-31/2008 (NDV-DR499-31/08), and the fusion proteins of three other related viruses from the Dominican Republic and Mexico were sequenced and phylogenetically analyzed. Genetically, these four isolates were highly distinct from all other currently known isolates of NDV, and together, they fulfill the newly established criteria for inclusion as a novel genotype of NDV (genotype XVI). The lack of any reported isolation of viruses related to this group since 1986 suggests that virulent viruses of this genotype may have evolved unnoticed for 22 years. The NDV-DR499-31/08 isolate had an intracerebral pathogenicity index (ICPI) score of 1.88, and sequencing of the fusion cleavage site identified multiple basic amino acids and a phenylalanine at position 117, indicating this isolate to be virulent. These results were further confirmed by a clinicopathological assessment in vivo. In 4-week-old chickens, NDV-DR499-31/08 behaved as a velogenic viscerotropic strain with systemic virus distribution and severe necrohemorrhagic lesions targeting mainly the intestine and the lymphoid organs. The clear phylogenetic relationship between the 2008, 1986, and 1947 ancestral viruses suggests that virulent NDV strains may have evolved in unknown reservoirs in the Caribbean and surrounding regions and underlines the importance of continued and improved epidemiological surveillance strategies to detect NDV in wild-bird species and commercial poultry.

Recombination and natural selection in hepatitis E virus genotypes

To gain new insights into the evolutionary processes that created the genetic diversity of the hepatitis E virus (HEV), the Recombination Detection Program (RDP) and SimPlot program were employed to detect recombination events in the genome, then the fixed-effects likelihood (FEL) method was used to detect natural selection effects on viral proteins. Recombination analysis provided strong evidence for both intergenotype and intragenotype recombination events in the sequences analyzed. Recombination events were found to be distributed non-randomly, with the highest frequency in the X domain and the helicase. Strain DQ450072 was identified as intergenotype-recombinant. Natural selection analysis revealed that codons under both negative selection and positive selection were distributed non-randomly. ORF1 and ORF2 have experienced strong purifying selection across genotypes. Furthermore, potentially important sites were also found under positive selection in the N-terminal end of ORF2 and the C-terminal end of ORF3. No significant difference was found among the selective pressures on different genotypes.

Molecular characterization of the measles virus genotypes in JiLin Province, China

Measles remains a severe global health threat, and nearly 30 million new cases are reported annually. Although many studies have analyzed measles viruses (MV) at the epidemiologic and phylogenetic levels, no study has yet to integrate these two types of data. To this end, we isolated 16 wild-type MV strains China's Jilin province. The MV genotype H1 was the most prevalent strain. After sequencing the nucleoprotein (N) genes of these strains, a maximum clade credibility tree was constructed by the Bayesian Markov Chain Monte Carlo method using 450 MV strains from GenBank with epidemiological information. The MV N gene evolution rate was 1.127E-3. Analysis of the time of the most recent common ancestor (TMRCA) for genotypes A/B/C/G/H revealed that genotypes D and B had the largest and smallest TMRCA (45.86 and 26.63, respectively). The highest level of genetic diversity for the MV N gene occurred around the year 2000. Here in this study, we uncovered the MV genotypes circulating in China's Jilin Province and estimated the epidemiologic and phylogenetic relationship for the six different genotypes of MV.

Complete genome sequence of arracacha mottle virus

Arracacha mottle virus (AMoV) is the only potyvirus reported to infect arracacha (Arracacia xanthorrhiza) in Brazil. Here, the complete genome sequence of an isolate of AMoV was determined to be 9,630 nucleotides in length, excluding the 3' poly-A tail, and encoding a polyprotein of 3,135 amino acids and a putative P3N-PIPO protein. Its genomic organization is typical of a member of the genus Potyvirus, containing all conserved motifs. Its full genome sequence shared 56.2 % nucleotide identity with sunflower chlorotic mottle virus and verbena virus Y, the most closely related viruses.

A combined multi-virulence-locus sequence typing and Staphylococcal Cassette Chromosome mec typing scheme possesses enhanced discriminatory power for genotyping MRSA

Methicillin-resistant Staphylococcus aureus (MRSA) remains a major threat to human populations worldwide. Knowing the extent of MRSA genetic diversity within a healthcare facility may provide important insights into the epidemiology of this important pathogen. MRSA isolates recovered from nasal swabs of patients entering the Intensive Care Unit of the Penn State Milton S. Hershey Medical Center, USA, from 2008 to 2009 were genotyped using Staphylococcal Cassette Chromosome mec (SCCmec), multilocus sequence typing (MLST) and a newly developed multi-virulence-locus sequence typing (MVLST) scheme. Sequence data for seven housekeeping genes (arcC, aroE, glpF, gmk, pta, tpi and yqiL) and six virulence genes (alt, essC, geh, hlgA, htrA and sdrC) were used for MLST and MVLST analyses, respectively. MLST identified 12 sequence types (STs) within the hospital isolates. One ST designated ST5 was the most common subtype (38.8%) followed by ST105 (22.4%) and ST8 (16.4%). In contrast, MVLST identified 29 STs (Virulence Types, VTs) from the same set of isolates, with VT6 (32.8%) being the predominant subtype followed by VT9 (8.9%) and VT2 (8.9%). Subsequent analysis of 25 MRSA isolates associated with an outbreak at a Pennsylvania state prison revealed all isolates were VT2 and SCCmec type IVa. These results suggest that a combination of MVLST and SCCmec typing may clarify the epidemiology of MRSA. Additional research with a more diverse set of strains and correlation with conventional epidemiologic data are needed to validate this new subtyping strategy.

Grapevine deformation virus: completion of the sequence and evidence on its origin from recombination events between Grapevine fanleaf virus and Arabis mosaic virus

The complete nucleotide (nt) sequence of Grapevine deformation virus (GDefV) RNA-1 has been determined. It consists of 7386 nt, excluding the poly(A) tail, and contains a single open reading frame (ORF) encoding a polyprotein (p1) of 252 kDa. P1 comprises the 1A(Pro-cof) proteinase cofactor, the 1B(Hel) NTP-binding protein, the 1C(VPg) viral protein genome-linked, the 1D(Prot) proteinase and the 1E(Pol) RNA-dependent RNA polymerase, all of which are conserved domains in polyproteins of different members of the order Picornavirales. The amino acid (aa) sequence of GDefV RNA1 p1 has the highest identity with the homologous products of Grapevine fanleaf virus (GFLV, 86-88%) and Arabis mosaic virus (ArMV, 73-74%), two nepoviruses of subgroup A. Four cleavage sites for proteins processing were predicted (C/A, C/S, G/E and R/G) and found similar to those of GFLV RNA1. Phylogenetic trees constructed with the complete aa sequences of protein p1 and the RNA2-encoded protein p2 of GDeFV, GFLV and ArMV, showed an incongruent allocation of GDefV in these trees. Pairwise alignment and prediction of recombination sites of both RNA segments showed that GDefV RNA2 has a mosaic structure resulting from recombination events between GFLV and ArMV at the level of the 2A(HP) (homing protein), 2B(MP) (movement protein), 2C(CP) (capsid protein) and the 3'NCR (non coding region). This strongly suggests that GDefV originated from the interspecific recombination between isolates of GFLV and ArMV.

Phylogenomics of T4 cyanophages: lateral gene transfer in the 'core' and origins of host genes

The last two decades have revealed that phages (viruses that infect bacteria) are abundant and play fundamental roles in the Earth System, with the T4-like myoviruses (herein T4-like phages) emerging as a dominant 'signal' in wild populations. Here we examine 27 T4-like phage genomes, with a focus on 17 that infect ocean picocyanobacteria (cyanophages), to evaluate lateral gene transfer (LGT) in this group. First, we establish a reference tree by evaluating concatenated core gene supertrees and whole genome gene content trees. Next, we evaluate what fraction of these 'core genes' shared by all 17 cyanophages appear prone to LGT. Most (47 out of 57 core genes) were vertically transferred as inferred from tree tests and genomic synteny. Of those 10 core genes that failed the tree tests, the bulk (8 of 10) remain syntenic in the genomes with only a few (3 of the 10) having identifiable signatures of mobile elements. Notably, only one of these 10 is shared not only by the 17 cyanophages, but also by all 27 T4-like phages (thymidylate synthase); its evolutionary history suggests cyanophages may be the origin of these genes to Prochlorococcus. Next, we examined intragenic recombination among the core genes and found that it did occur, even among these core genes, but that the rate was significantly higher between closely related phages, perhaps reducing any detectable LGT signal and leading to taxon cohesion. Finally, among 18 auxiliary metabolic genes (AMGs, a.k.a. 'host' genes), we found that half originated from their immediate hosts, in some cases multiple times (e.g. psbA, psbD, pstS), while the remaining have less clear evolutionary origins ranging from cyanobacteria (4 genes) or microbes (5 genes), with particular diversity among viral TalC and Hsp20 sequences. Together, these findings highlight the patterns and limits of vertical evolution, as well as the ecological and evolutionary roles of LGT in shaping T4-like phage genomes.

Ultrafast evolution and loss of CRISPRs following a host shift in a novel wildlife pathogen, Mycoplasma gallisepticum

Measureable rates of genome evolution are well documented in human pathogens but are less well understood in bacterial pathogens in the wild, particularly during and after host switches. Mycoplasma gallisepticum (MG) is a pathogenic bacterium that has evolved predominantly in poultry and recently jumped to wild house finches (Carpodacus mexicanus), a common North American songbird. For the first time we characterize the genome and measure rates of genome evolution in House Finch isolates of MG, as well as in poultry outgroups. Using whole-genome sequences of 12 House Finch isolates across a 13-year serial sample and an additional four newly sequenced poultry strains, we estimate a nucleotide diversity in House Finch isolates of only ∼2% of ancestral poultry strains and a nucleotide substitution rate of 0.8-1.2×10(-5) per site per year both in poultry and in House Finches, an exceptionally fast rate rivaling some of the highest estimates reported thus far for bacteria. We also found high diversity and complete turnover of CRISPR arrays in poultry MG strains prior to the switch to the House Finch host, but after the invasion of House Finches there is progressive loss of CRISPR repeat diversity, and recruitment of novel CRISPR repeats ceases. Recent (2007) House Finch MG strains retain only ∼50% of the CRISPR repertoire founding (1994-95) strains and have lost the CRISPR-associated genes required for CRISPR function. Our results suggest that genome evolution in bacterial pathogens of wild birds can be extremely rapid and in this case is accompanied by apparent functional loss of CRISPRs.

Complete genome and clinicopathological characterization of a virulent Newcastle disease virus isolate from South America

Newcastle disease (ND) is one of the most important diseases of poultry, negatively affecting poultry production worldwide. The disease is caused by Newcastle disease virus (NDV) or avian paramyxovirus type 1 (APMV-1), a negative-sense single-stranded RNA virus of the genus Avulavirus, family Paramyxoviridae. Although all NDV isolates characterized to date belong to a single serotype of APMV-1, significant genetic diversity has been described between different NDV isolates. Here we present the complete genome sequence and the clinicopathological characterization of a virulent Newcastle disease virus isolate (NDV-Peru/08) obtained from poultry during an outbreak of ND in Peru in 2008. Phylogenetic reconstruction and analysis of the evolutionary distances between NDV-Peru/08 and other isolates representing established NDV genotypes revealed the existence of large genomic and amino differences that clearly distinguish this isolate from viruses of typical NDV genotypes. Although NDV-Peru/08 is a genetically distinct virus, pathogenesis studies conducted with chickens revealed that NDV-Peru/08 infection results in clinical signs characteristic of velogenic viscerotropic NDV strains. Additionally, vaccination studies have shown that an inactivated NDV-LaSota/46 vaccine conferred full protection from NDV-Peru/08-induced clinical disease and mortality. This represents the first complete characterization of a virulent NDV isolate from South America.

Recombination in avian gamma-coronavirus infectious bronchitis virus

Recombination in the family Coronaviridae has been well documented and is thought to be a contributing factor in the emergence and evolution of different coronaviral genotypes as well as different species of coronavirus. However, there are limited data available on the frequency and extent of recombination in coronaviruses in nature and particularly for the avian gamma-coronaviruses where only recently the emergence of a turkey coronavirus has been attributed solely to recombination. In this study, the full-length genomes of eight avian gamma-coronavirus infectious bronchitis virus (IBV) isolates were sequenced and along with other full-length IBV genomes available from GenBank were analyzed for recombination. Evidence of recombination was found in every sequence analyzed and was distributed throughout the entire genome. Areas that have the highest occurrence of recombination are located in regions of the genome that code for nonstructural proteins 2, 3 and 16, and the structural spike glycoprotein. The extent of the recombination observed, suggests that this may be one of the principal mechanisms for generating genetic and antigenic diversity within IBV. These data indicate that reticulate evolutionary change due to recombination in IBV, likely plays a major role in the origin and adaptation of the virus leading to new genetic types and strains of the virus.

New insights into the evolution of Wolbachia infections in filarial nematodes inferred from a large range of screened species

BACKGROUND

Wolbachia are intriguing symbiotic endobacteria with a peculiar host range that includes arthropods and a single nematode family, the Onchocercidae encompassing agents of filariases. This raises the question of the origin of infection in filariae. Wolbachia infect the female germline and the hypodermis. Some evidences lead to the theory that Wolbachia act as mutualist and coevolved with filariae from one infection event: their removal sterilizes female filariae; all the specimens of a positive species are infected; Wolbachia are vertically inherited; a few species lost the symbiont. However, most data on Wolbachia and filaria relationships derive from studies on few species of Onchocercinae and Dirofilariinae, from mammals.

METHODOLOGY/PRINCIPAL FINDINGS

We investigated the Wolbachia distribution testing 35 filarial species, including 28 species and 7 genera and/or subgenera newly screened, using PCR, immunohistochemical staining, whole mount fluorescent analysis, and cocladogenesis analysis. (i) Among the newly screened Onchocercinae from mammals eight species harbour Wolbachia but for some of them, bacteria are absent in the hypodermis, or in variable density. (ii) Wolbachia are not detected in the pathological model Monanema martini and in 8, upon 9, species of Cercopithifilaria. (iii) Supergroup F Wolbachia is identified in two newly screened Mansonella species and in Cercopithifilaria japonica. (iv) Type F Wolbachia infect the intestinal cells and somatic female genital tract. (v) Among Oswaldofilariinae, Waltonellinae and Splendidofilariinae, from saurian, anuran and bird respectively, Wolbachia are not detected.

CONCLUSIONS/SIGNIFICANCE

The absence of Wolbachia in 63% of onchocercids, notably in the ancestral Oswaldofilariinae estimated 140 mya old, the diverse tissues or specimens distribution, and a recent lateral transfer in supergroup F Wolbachia, modify the current view on the role and evolution of the endosymbiont and their hosts. Further genomic analyses on some of the newly sampled species are welcomed to decipher the open questions.

Detection of known and novel adenoviruses in cattle wastes via broad-spectrum primers

The critical assessment of bovine adenoviruses (BAdV) as indicators of environmental fecal contamination requires improved knowledge of their prevalence, shedding dynamics, and genetic diversity. We examined DNA extracted from bovine and other animal waste samples collected in Wisconsin for atadenoviruses and mastadenoviruses using novel, broad-spectrum PCR primer sets. BAdV were detected in 13% of cattle fecal samples, 90% of cattle urine samples, and 100% of cattle manure samples; 44 percent of BAdV-positive samples contained both Atadenovirus and Mastadenovirus DNA. Additionally, BAdV were detected in soil, runoff water from a cattle feedlot, and residential well water. Overall, we detected 8 of 11 prototype BAdV, plus bovine, rabbit, and porcine mastadenoviruses that diverged significantly from previously reported genotypes. The prevalence of BAdV shedding by cattle supports targeting AdV broadly as indicators of the presence of fecal contamination in aqueous environments. Conversely, several factors complicate the use of AdV for fecal source attribution. Animal AdV infecting a given livestock host were not monophyletic, recombination among livestock mastadenoviruses was detected, and the genetic diversity of animal AdV is still underreported. These caveats highlight the need for continuing genetic surveillance for animal AdV and for supporting data when BAdV detection is invoked for fecal source attribution in environmental samples. To our knowledge, this is the first study to report natural BAdV excretion in urine, BAdV detection in groundwater, and recombination in AdV of livestock origin.

Novel molecular markers of Chlamydia pecorum genetic diversity in the koala (Phascolarctos cinereus)

Background

Chlamydia pecorum is an obligate intracellular bacterium and the causative agent of reproductive and ocular disease in several animal hosts including koalas, sheep, cattle and goats. C. pecorum strains detected in koalas are genetically diverse, raising interesting questions about the origin and transmission of this species within koala hosts. While the ompA gene remains the most widely-used target in C. pecorum typing studies, it is generally recognised that surface protein encoding genes are not suited for phylogenetic analysis and it is becoming increasingly apparent that the ompA gene locus is not congruent with the phylogeny of the C. pecorum genome. Using the recently sequenced C. pecorum genome sequence (E58), we analysed 10 genes, including ompA, to evaluate the use of ompA as a molecular marker in the study of koala C. pecorum genetic diversity.

Results

Three genes (incA, ORF663, tarP) were found to contain sufficient nucleotide diversity and discriminatory power for detailed analysis and were used, with ompA, to genotype 24 C. pecorum PCR-positive koala samples from four populations. The most robust representation of the phylogeny of these samples was achieved through concatenation of all four gene sequences, enabling the recreation of a "true" phylogenetic signal. OmpA and incA were of limited value as fine-detailed genetic markers as they were unable to confer accurate phylogenetic distinctions between samples. On the other hand, the tarP and ORF663 genes were identified as useful "neutral" and "contingency" markers respectively, to represent the broad evolutionary history and intra-species genetic diversity of koala C. pecorum. Furthermore, the concatenation of ompA, incA and ORF663 sequences highlighted the monophyletic nature of koala C. pecorum infections by demonstrating a single evolutionary trajectory for koala hosts that is distinct from that seen in non-koala hosts.

Conclusions

While the continued use of ompA as a fine-detailed molecular marker for epidemiological analysis appears justified, the tarP and ORF663 genes also appear to be valuable markers of phylogenetic or biogeographic divisions at the C. pecorum intra-species level. This research has significant implications for future typing studies to understand the phylogeny, genetic diversity, and epidemiology of C. pecorum infections in the koala and other animal species.

New perspectives on the evolutionary history of hepatitis B virus genotype F

Hepatitis B virus (HBV) is a globally distributed human pathogen. The aim of this work was to analyze the evolutionary history of HBV genotype F, emphasizing on the study of subgenotypes prevalent in the Southern area of South America. Complete genomes of HBV genotype F from 36 samples from Argentina and Chile were sequenced and analyzed by phylogenetic and Bayesian coalescent methods along with sequences obtained from GenBank database. The phylogeography separated not only Central American from South American isolates but also revealed that different subgenotypes are distributed in constrained although not exclusive areas of the continent. The result obtained with time-stamped complete genomes failed to explain the wide geographical distribution and the clustering observed in this genotype. Conversely, the use of Bayesian coalescent analyses with substitution rates as priors, instead of the co-estimation of tMRCA and substitution rate, allowed us to propose a far origin for the HBV genotype F based on the phylogeographical and epidemiological data.

The evolutionary history of Beet necrotic yellow vein virus deduced from genetic variation, geographical origin and spread, and the breaking of host resistance

Beet necrotic yellow vein virus (BNYVV) is an economically important pathogen of sugar beet and has been found worldwide, probably as the result of recent worldwide spread. The BNYVV genome consists of four or five RNA components. Here, we report analysis of sequence variation in the RNA3-p25, RNA4-p31, RNA2-CP, and RNA5-p26 genes of 73 worldwide isolates. The RNA3-p25 gene encodes virulence and avirulence factors. These four sets of gene sequences each fell into two to four groups, of which the three groups of p25 formed eight subgroups with different geographical distributions. Each of these subgroup isolates (strains) could have arisen from four original BNYVV population and their mixed infections. The genetic diversity for BNYVV was relatively small. Selection pressure varied greatly depending on the BNYVV gene and geographical location. Isolates of the Italy strain, in which p25 was subject to the strongest positive selection, were able to overcome the Rz1-host resistance gene to differing degrees, whereas other geographically limited strains could not. Resistance-breaking variants were generated by p25 amino acid changes at positions 67 and 68. Our studies suggest that BNYVV originally evolved in East Asia and has recently become a pathogen of cultivated sugar beet followed by the emergence of new resistance-breaking variants.