jpHMM: Improving the reliability of recombination prediction in HIV-1

Previously, we developed jumping profile hidden Markov model (jpHMM), a new method to detect recombinations in HIV-1 genomes. The jpHMM predicts recombination breakpoints in a query sequence and assigns to each position of the sequence one of the major HIV-1 subtypes. Since incorrect subtype assignment or recombination prediction may lead to wrong conclusions in epidemiological or vaccine research, information about the reliability of the predicted parental subtypes and breakpoint positions is valuable. For this reason, we extended the output of jpHMM to include such information in terms of ‘uncertainty’ regions in the recombination prediction and an interval estimate of the breakpoint. Both types of information are computed based on the posterior probabilities of the subtypes at each query sequence position. Our results show that this extension strongly improves the reliability of the jpHMM recombination prediction. The jpHMM is available online at http://jphmm.gobics.de/

Detection of novel astroviruses in urban brown rats and previously known astroviruses in humans

Several novel astroviruses have been recently discovered in humans and in other animals. Here, we report results from our surveillance of astroviruses in human and rodent faecal samples in Hong Kong. Classical human astroviruses (n=9) and a human MLB1 astrovirus were detected in human faecal samples (n=622). Novel astroviruses were detected from 1.6 % of the faecal samples of urban brown rat (Rattus norvegicus) (n=441), indicating the prevalence of astrovirus infection in rats might be much lower than that recently observed in bats. These rat astroviruses were phylogenetically related to recently discovered human astroviruses MLB1 and MLB2, suggesting that the MLB viruses and these novel rat astroviruses may share a common ancestor

Turnip mosaic potyvirus probably first spread to Eurasian brassica crops from wild orchids about 1000 years ago

Turnip mosaic potyvirus (TuMV) is probably the most widespread and damaging virus that infects cultivated brassicas worldwide. Previous work has indicated that the virus originated in western Eurasia, with all of its closest relatives being viruses of monocotyledonous plants. Here we report that we have identified a sister lineage of TuMV-like potyviruses (TuMV-OM) from European orchids. The isolates of TuMV-OM form a monophyletic sister lineage to the brassica-infecting TuMVs (TuMV-BIs), and are nested within a clade of monocotyledon-infecting viruses. Extensive host-range tests showed that all of the TuMV-OMs are biologically similar to, but distinct from, TuMV-BIs and do not readily infect brassicas. We conclude that it is more likely that TuMV evolved from a TuMV-OM-like ancestor than the reverse. We did Bayesian coalescent analyses using a combination of novel and published sequence data from four TuMV genes [helper component-proteinase protein (HC-Pro), protein 3(P3), nuclear inclusion b protein (NIb), and coat protein (CP)]. Three genes (HC-Pro, P3, and NIb), but not the CP gene, gave results indicating that the TuMV-BI viruses diverged from TuMV-OMs around 1000 years ago. Only 150 years later, the four lineages of the present global population of TuMV-BIs diverged from one another. These dates are congruent with historical records of the spread of agriculture in Western Europe. From about 1200 years ago, there was a warming of the climate, and agriculture and the human population of the region greatly increased. Farming replaced woodlands, fostering viruses and aphid vectors that could invade the crops, which included several brassica cultivars and weeds. Later, starting 500 years ago, inter-continental maritime trade probably spread the TuMV-BIs to the remainder of the world

Accurate Detection of Recombinant Breakpoints in Whole-Genome Alignments

We propose a novel method for detecting sites of molecular recombination in multiple alignments. Our approach is a compromise between previous extremes of computationally prohibitive but mathematically rigorous methods and imprecise heuristic methods. Using a combined algorithm for estimating tree structure and hidden Markov model parameters, our program detects changes in phylogenetic tree topology over a multiple sequence alignment. We evaluate our method on benchmark datasets from previous studies on two recombinant pathogens, Neisseria and HIV-1, as well as simulated data. We show that we are not only able to detect recombinant regions of vastly different sizes but also the location of breakpoints with great accuracy. We show that our method does well inferring recombination breakpoints while at the same time maintaining practicality for larger datasets. In all cases, we confirm the breakpoint predictions of previous studies, and in many cases we offer novel predictions

HIV-1 Epidemic in the Caribbean Is Dominated by Subtype B

The molecular epidemiology of HIV-1 in the Caribbean has been described using partial genome sequencing; subtype B is the most common subtype in multiple countries. To expand our knowledge of this, nearly full genome amplification, sequencing and analysis was conducted.Virion RNA from sera collected in Haiti, Dominican Republic, Jamaica and Trinidad and Tobago were reverse transcribed, PCR amplified, sequenced and phylogenetically analyzed. Nearly full genomes were completed for 15 strains; partial pol was done for 67 strains. All but one of the 67 strains analyzed in pol were subtype B; the exception was a unique recombinant of subtypes B and C collected in the Dominican Republic. Of the nearly full genomes of 14 strains that were subtype B in pol, all were subtype B from one end of the genome to the other and not inter-subtype recombinants. Surprisingly, the Caribbean subtype B strains clustered significantly with each other and separate from subtype B from other parts of the pandemic.The more complete analysis of HIV-1 from 4 Caribbean countries confirms previous research using partial genome analysis that the predominant subtype in circulation was subtype B. The Caribbean strains are phylogenetically distinct from other subtype B strains although the biological meaning of this finding is unclear

Classification of HIV-1 Sequences Using Profile Hidden Markov Models

Accurate classification of HIV-1 subtypes is essential for studying the dynamic spatial distribution pattern of HIV-1 subtypes and also for developing effective methods of treatment that can be targeted to attack specific subtypes. We propose a classification method based on profile Hidden Markov Model that can accurately identify an unknown strain. We show that a standard method that relies on the construction of a positive training set only, to capture unique features associated with a particular subtype, can accurately classify sequences belonging to all subtypes except B and D. We point out the drawbacks of the standard method; namely, an arbitrary choice of threshold to distinguish between true positives and true negatives, and the inability to discriminate between closely related subtypes. We then propose an improved classification method based on construction of a positive as well as a negative training set to improve discriminating ability between closely related subtypes like B and D. Finally, we show how the improved method can be used to accurately determine the subtype composition of Common Recombinant Forms of the virus that are made up of two or more subtypes. Our method provides a simple and highly accurate alternative to other classification methods and will be useful in accurately annotating newly sequenced HIV-1 strains

HIV-1 recombinants with multiple parental strains in low-prevalence, remote regions of Cameroon: Evolutionary relics?

<p>Abstract</p> <p>Background</p> <p>The HIV pandemic disseminated globally from Central West Africa, beginning in the second half of the twentieth century. To elucidate the virologic origins of the pandemic, a cross-sectional study was conducted of the genetic diversity of HIV-1 strains in villagers in 14 remote locations in Cameroon and in hospitalized and STI patients. DNA extracted from PBMC was PCR amplified from HIV(+) subjects. Partial <it>pol </it>amplicons (N = 164) and nearly full virus genomes (N = 78) were sequenced. Among the 3956 rural villagers studied, the prevalence of HIV infection was 4.9%; among the hospitalized and clinic patients, it was 8.6%.</p> <p>Results</p> <p>Virus genotypes fell into two distinctive groups. A majority of the genotyped strains (109/164) were the circulating recombinant form (CRF) known to be endemic in West Africa and Central West Africa, CRF02_AG. The second most common genetic form (9/164) was the recently described CRF22_01A1, and the rest were a collection of 4 different subtypes (A2, D, F2, G) and 6 different CRFs (-01, -11, -13, -18, -25, -37). Remarkably, 10.4% of HIV-1 genomes detected (17/164) were heretofore undescribed unique recombinant forms (URF) present in only a single person. Nearly full genome sequencing was completed for 78 of the viruses of interest. HIV genetic diversity was commonplace in rural villages: 12 villages each had at least one newly detected URF, and 9 villages had two or more.</p> <p>Conclusions</p> <p>These results show that while CRF02_AG dominated the HIV strains in the rural villages, the remainder of the viruses had tremendous genetic diversity. Between the trans-species transmission of SIV_cpzand the dispersal of pandemic HIV-1, there was a time when we hypothesize that nascent HIV-1 was spreading, but only to a limited extent, recombining with other local HIV-1, creating a large variety of recombinants. When one of those recombinants began to spread widely (i.e. became epidemic), it was recognized as a subtype. We hypothesize that the viruses in these remote Cameroon villages may represent that pre-epidemic stage of viral evolution.</p

Genetic characterization of 2008 reassortant influenza A virus (H5N1), Thailand

In January and November 2008, outbreaks of avian influenza have been reported in 4 provinces of Thailand. Eight Influenza A H5N1 viruses were recovered from these 2008 AI outbreaks and comprehensively characterized and analyzed for nucleotide identity, genetic relatedness, virulence determinants, and possible sites of reassortment. The results show that the 2008 H5N1 viruses displayed genetic drift characteristics (less than 3% genetic differences), as commonly found in influenza A viruses. Based on phylogenetic analysis, clade 1 viruses in Thailand were divided into 3 distinct branches (subclades 1, 1.1 and 1.2). Six out of 8 H5N1 isolates have been identified as reassorted H5N1 viruses, while other isolates belong to an original H5N1 clade. These viruses have undergone inter-lineage reassortment between subclades 1.1 and 1.2 and thus represent new reassorted 2008 H5N1 viruses. The reassorted viruses have acquired gene segments from H5N1, subclade 1.1 (PA, HA, NP and M) and subclade 1.2 (PB2, PB1, NA and NS) in Thailand. Bootscan analysis of concatenated whole genome sequences of the 2008 H5N1 viruses supported the reassortment sites between subclade 1.1 and 1.2 viruses. Based on estimating of the time of the most recent common ancestors of the 2008 H5N1 viruses, the potential point of genetic reassortment of the viruses could be traced back to 2006. Genetic analysis of the 2008 H5N1 viruses has shown that most virulence determinants in all 8 genes of the viruses have remained unchanged. In summary, two predominant H5N1 lineages were circulating in 2008. The original CUK2-like lineage mainly circulated in central Thailand and the reassorted lineage (subclades 1.1 and 1.2) predominantly circulated in lower-north Thailand. To prevent new reassortment, emphasis should be put on prevention of H5N1 viruses circulating in high risk areas. In addition, surveillance and whole genome sequencing of H5N1 viruses should be routinely performed for monitoring the genetic drift of the virus and new reassorted strains, especially in light of potential reassortment between avian and mammalian H5N1 viruses

Mixed emotions to near-miss outcomes: a psychophysiological study with facial electromyography

Near-misses occur across many forms of gambling and are rated as unpleasant while simultaneously increasing the motivation to continue playing. On slot machines, the icon position relative to the payline moderates the effects of near-misses, with near-misses before the payline increasing motivation, and near-misses after the payline being rated as aversive. Near-misses are also known to increase physiological arousal compared to full-misses, but physiological measures to date have not been able to dissociate positive and negative emotional responses. The present study measured facial electromyography at the corrugator (brow) and zygomaticus (cheek) sites, as well as electrodermal activity (EDA), following gambling outcomes on a two-reel slot machine simulation in 77 novice gamblers. Behavioral data was collected using trial-by-trial ratings of motivation and valence. Wins were rated as more pleasant and increased motivation to continue playing, compared to non-win outcomes. Wins were also accompanied by increased EDA and zygomaticus activity. Near-misses after the payline were rated as more aversive than other non-wins, and this was accompanied by increased EDA and zygomaticus activity. Near-misses before the payline increased motivation to continue playing, and were accompanied by increased EDA. Thus, both subjective and physiological responses to near-misses differ for events falling either side of the payline. The ‘near-miss effect’ is not a unitary phenomenon. Facial EMG has differential sensitivity to positive and negative valence and may be a useful measure for future studies of gambling behavior

Geographical and Ethnic Distribution of the HBV C/D Recombinant on the Qinghai-Tibet Plateau

Two forms of hepatitis B virus (HBV) C/D recombinant have been identified in western China, but little is known about their geographical and ethnic distributions, and particularly the clinical significance and specific mutations in the pre-core region. To address these questions, a total of 624 chronic HBV carriers from four ethnic populations representing five provinces in western China were enrolled in this study. Genotypes were firstly determined by restriction fragment length polymorphism, and then confirmed by full or partial genome nucleotide sequencing. The distribution of HBV genotypes was as follows: HBV/B: 40 (6.4%); HBV/C: 221 (35.4%); HBV/D: 39 (6.3%); HBV/CD: 324 (51.9%). In the 324 HBV C/D recombinant infections, 244 (75.3%) were infected with the “CD1” and 80 (24.7%) were infected with the “CD2.” The distribution of HBV genotypes exhibited distinct patterns in different regions and ethnic populations. Geographically, the C/D recombinant was the most prevalent HBV strain on the Qinghai-Tibet Plateau. Ethnically, the C/D recombinant had a higher prevalence in Tibetan patients than in other populations. Clinically, patients with HBV/CD1 showed significantly lower levels of serum total bilirubin than patients with HBV/C2. The prevalence of HBeAg was comparable between patients with HBV/CD1 and HBV/C2 (63.3% vs 50.0%, P = 0.118) whether patients were taken together or stratified by age into three groups (65.6% vs 58.8% in <30 years, P = 0.758; 61.9% vs 48.0% in 30–50 years, P = 0.244; 64.3% vs 33.3%, P = 0.336). Virologically HBV/CD1 had a significantly lower frequency of G1896A than HBV/C2. In conclusion, the HBV C/D recombinant is restricted to the Qinghai-Tibet Plateau in western China and is found predominantly in Tibetans. The predominance of the premature pre-core stop mutation G1896A in patients with the HBV C/D recombinant may account for the higher prevalence of HBeAg in these patients