Characterization of Dengue Virus Type 2: New Insights on the 2010 Brazilian Epidemic

Dengue viruses (DENV) serotypes 1, 2, and 3 have been causing yearly outbreaks in Brazil. In this study, we report the reintroduction of DENV2 in the coast of Sao Paulo State. Partial envelope viral genes were sequenced from eighteen patients with dengue fever during the 2010 epidemic. Phylogenetic analysis showed this strain belongs to the American/Asian genotype and was closely related to the virus that circulated in Rio de Janeiro in 2007 and 2008. The phylogeny also showed no clustering by clinical presentation, suggesting that the disease severity could not be explained by distinct variants or genotypes. The time of the most recent common ancestor of American/Asian genotype and the Sao Paulo and Rio de Janeiro (SP/RJ) monophyletic cluster was estimated to be around 40 and 10 years, respectively. Since this virus was first identified in Brazil in 2007, we suggest that it was already circulating in the country before causing the first documented outbreak. This is the first description of the 2010 outbreak in the State of Sao Paulo, Brazil, and should contribute to efforts to control and monitor the spread of DENVs in endemic areas.

Spatio-Temporal Tracking and Phylodynamics of an Urban Dengue 3 Outbreak in Sao Paulo, Brazil

The dengue virus has a single-stranded positive-sense RNA genome of similar to 10.700 nucleotides with a single open reading frame that encodes three structural (C, prM, and E) and seven nonstructural (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) proteins. It possesses four antigenically distinct serotypes (DENV 1-4). Many phylogenetic studies address particularities of the different serotypes using convenience samples that are not conducive to a spatio-temporal analysis in a single urban setting. We describe the pattern of spread of distinct lineages of DENV-3 circulating in Sao Jose do Rio Preto, Brazil, during 2006. Blood samples from patients presenting dengue-like symptoms were collected for DENV testing. We performed M-N-PCR using primers based on NS5 for virus detection and identification. The fragments were purified from PCR mixtures and sequenced. The positive dengue cases were geo-coded. To type the sequenced samples, 52 reference sequences were aligned. The dataset generated was used for iterative phylogenetic reconstruction with the maximum likelihood criterion. The best demographic model, the rate of growth, rate of evolutionary change, and Time to Most Recent Common Ancestor (TMRCA) were estimated. The basic reproductive rate during the epidemics was estimated. We obtained sequences from 82 patients among 174 blood samples. We were able to geo-code 46 sequences. The alignment generated a 399-nucleotide-long dataset with 134 taxa. The phylogenetic analysis indicated that all samples were of DENV-3 and related to strains circulating on the isle of Martinique in 2000-2001. Sixty DENV-3 from Sao Jose do Rio Preto formed a monophyletic group (lineage 1), closely related to the remaining 22 isolates (lineage 2). We assumed that these lineages appeared before 2006 in different occasions. By transforming the inferred exponential growth rates into the basic reproductive rate, we obtained values for lineage 1 of R(0) = 1.53 and values for lineage 2 of R(0) = 1.13. Under the exponential model, TMRCA of lineage 1 dated 1 year and lineage 2 dated 3.4 years before the last sampling. The possibility of inferring the spatio-temporal dynamics from genetic data has been generally little explored, and it may shed light on DENV circulation. The use of both geographic and temporally structured phylogenetic data provided a detailed view on the spread of at least two dengue viral strains in a populated urban area.

Syphilis at the Crossroad of Phylogenetics and Paleopathology

The origin of syphilis is still controversial. Different research avenues explore its fascinating history. Here we employed a new integrative approach, where paleopathology and molecular analyses are combined. As an exercise to test the validity of this approach we examined different hypotheses on the origin of syphilis and other human diseases caused by treponemes (treponematoses). Initially, we constructed a worldwide map containing all accessible reports on palaeopathological evidences of treponematoses before Columbus's return to Europe. Then, we selected the oldest ones to calibrate the time of the most recent common ancestor of Treponema pallidum subsp. pallidum, T. pallidum subsp. endemicum and T. pallidum subsp. pertenue in phylogenetic analyses with 21 genetic regions of different T. pallidum strains previously reported. Finally, we estimated the treponemes' evolutionary rate to test three scenarios: A) if treponematoses accompanied human evolution since Homo erectus; B) if venereal syphilis arose very recently from less virulent strains caught in the New World about 500 years ago, and C) if it emerged in the Americas between 16,500 and 5,000 years ago. Two of the resulting evolutionary rates were unlikely and do not explain the existent osseous evidence. Thus, treponematoses, as we know them today, did not emerge with H. erectus, nor did venereal syphilis appear only five centuries ago. However, considering 16,500 years before present (yBP) as the time of the first colonization of the Americas, and approximately 5,000 yBP as the oldest probable evidence of venereal syphilis in the world, we could not entirely reject hypothesis C. We confirm that syphilis seems to have emerged in this time span, since the resulting evolutionary rate is compatible with those observed in other bacteria. In contrast, if the claims of precolumbian venereal syphilis outside the Americas are taken into account, the place of origin remains unsolved. Finally, the endeavor of joining paleopathology and phylogenetics proved to be a fruitful and promising approach for the study of infectious diseases.

Molecular characterization of the Hepatitis B virus genotypes in Colombia: A Bayesian inference on the genotype F

Hepatitis B is a worldwide health problem affecting about 2 billion people and more than 350 million are chronic carriers of the virus. Nine HBV genotypes (A to I) have been described. The geographical distribution of HBV genotypes is not completely understood due to the limited number of samples from some parts of the world. One such example is Colombia, in which few studies have described the HBV genotypes. In this study, we characterized HBV genotypes in 143 HBsAg-positive volunteer blood donors from Colombia. A fragment of 1306 bp partially comprising HBsAg and the DNA polymerase coding regions (S/POL) was amplified and sequenced. Bayesian phylogenetic analyses were conducted using the Markov Chain Monte Carlo (MCMC) approach to obtain the maximum clade credibility (MCC) tree using BEAST v.1.5.3. Of all samples, 68 were positive and 52 were successfully sequenced. Genotype F was the most prevalent in this population (77%) - subgenotypes F3 (75%) and Fib (2%). Genotype G (7.7%) and subgenotype A2 (15.3%) were also found. Genotype G sequence analysis suggests distinct introductions of this genotype in the country. Furthermore, we estimated the time of the most recent common ancestor (TMRCA) for each HBV/F subgenotype and also for Colombian F3 sequences using two different datasets: (i) 77 sequences comprising 1306 bp of S/POL region and (ii) 283 sequences comprising 681 bp of S/POL region. We also used two other previously estimated evolutionary rates: (i) 2.60 x 10(-4) s/s/y and (ii) 1.5 x 10(-5) s/s/y. Here we report the HBV genotypes circulating in Colombia and estimated the TMRCA for the four different subgenotypes of genotype F. (C) 2010 Elsevier B.V. All rights reserved.

Molecular epidemiology and genetic diversity of hepatitis B virus genotype E in an isolated Afro-Colombian community

Hepatitis B virus (HBV) infection is a significant public health concern with 350 million chronic carriers worldwide. Eight HBV genotypes (A-H) have been described so far. Genotype E (HBV/E) is widely distributed in West Africa and has rarely been found in other continents, except for a few cases in individuals with an African background. In this study, we characterized HBV genotypes in Quibdo, Colombia, by partial S/P gene sequencing, and found, for the first time, HBV/E circulating in nine Afro-Colombian patients who had no recent contact with Africa. The presence of HBV/E in this community as a monophyletic group suggests that it was a result of a recent introduction by some Afro-descendent contact or, alternatively, that the virus came with slaves brought to Colombia. By using sequences with sampling dates, we estimated the substitution rate to be about 3.2x10(-4) substitutions per site per year, which resulted in a time to the most recent common ancestor (TMRCA) of 29 years. In parallel, we also estimated the TMRCA for HBV/E by using two previously estimated substitution rates (7.7x10(-4) and 1.5x10(-5) substitutions per site per year). The TMRCA was around 35 years under the higher rate and 1500 years under the slower rate. In sum, this work reports for the first time the presence of an exclusively African HBV genotype circulating in South America. We also discuss the time of the entry of this virus into America based on different substitution rates estimated for HBV.

Dynamics of Hepatitis D (delta) virus genotype 3 in the Amazon region of South America

Hepatitis delta virus (HDV) is widely distributed and associated with fulminant hepatitis epidemics in areas with high prevalence of HBV. Several studies performed in the 1980s showed data on HDV infection in South America, but there are no studies on the viral dynamics of this virus. The aim of this study was to conduct an evolutionary analysis of hepatitis delta genotype 3 (HDV/3) prevalent in South America: estimate its nucleotide substitution rate, determine the time of most recent ancestor (TMRCA) and characterize the epidemic history and evolutionary dynamics. Furthermore, we characterized the presence of HBV/HDV infection in seven samples collected from patients who died due to fulminant hepatitis from Amazon region in Colombia and included them in the evolutionary analysis. This is the first study reporting HBV and HDV sequences from the Amazon region of Colombia. Of the seven Colombian patients, five were positive for HBV-DNA and HDV-RNA. Of them, two samples were successfully sequenced for HBV (subgenotypes F3 and Fib) and the five samples HDV positive were classified as HDV/3. By using all HDV/3 available reference sequences with sampling dates (n = 36), we estimated the HDV/3 substitution rate in 1.07 x 10(-3) substitutions per site per year (s/s/y), which resulted in a time to the most recent common ancestor (TMRCA) of 85 years. Also, it was determined that HDV/3 spread exponentially from early 1950s to the 1970s in South America. This work discusses for the first time the viral dynamics for the HDV/3 circulating in South America. We suggest that the measures implemented to control HBV transmission resulted in the control of HDV/3 spreading in South America, especially after the important raise in this infection associated with a huge mortality during the 1950s up to the 1970s. The differences found among HDV/3 and the other HDV genotypes concerning its diversity raises the hypothesis of a different origin and/or a different transmission route. (C) 2011 Elsevier B.V. All rights reserved.

Duck hepatitis B virus expresses a regulatory HBx-like protein from a hidden open reading frame

Duck hepatitis B viruses (DHBV), unlike mammalian hepadnaviruses, are thought to lack X genes, which encode transcription-regulatory proteins believed to contribute to the development of hepatocellular carcinoma. A lack of association of chronic DHBV infection with hepatocellular carcinoma development supports this belief. Here, we demonstrate that DHBV genomes have a hidden open reading frame from which a transcription-regulatory protein, designated DHBx, is expressed both in vitro and in vivo. We show that DHBx enhances neither viral protein expression, intracellular DNA synthesis, nor virion production when assayed in the full-length genome context in LMH cells. However, similar to mammalian hepadnavirus X proteins, DHBx activates cellular and viral promoters via the Raf-mitogen-activated protein kinase signaling pathway and localizes primarily in the cytoplasm. The functional similarities as,well as the weak sequence homologies of DHBx and the X proteins of mammalian hepadnaviruses strongly suggest a common ancestry of ortho- and avihepadnavirus X genes. In addition, our data disclose similar intracellular localization and transcription regulatory functions of the corresponding proteins, raise new questions as to their presumed role in hepatocarcinogenesis, and imply unique opportunities for deciphering of their still-enigmatic in vivo functions.

Inferring genome trees by using a filter to eliminate phylogenetically discordant sequences and a distance matrix based on mean normalized BLASTP scores

Darwin's paradigm holds that the diversity of present-day organisms has arisen via a process of genetic descent with modification, as on a bifurcating tree. Evidence is accumulating that genes are sometimes transferred not along lineages but rather across lineages. To the extent that this is so, Darwin's paradigm can apply only imperfectly to genomes, potentially complicating or perhaps undermining attempts to reconstruct historical relationships among genomes (i.e., a genome tree). Whether most genes in a genome have arisen via treelike (vertical) descent or by lateral transfer across lineages can be tested if enough complete genome sequences are used. We define a phylogenetically discordant sequence (PDS) as an open reading frame (ORF) that exhibits patterns of similarity relationships statistically distinguishable from those of most other ORFs in the same genome. PDSs represent between 6.0 and 16.8% (mean, 10.8%) of the analyzable ORFs in the genomes of 28 bacteria, eight archaea, and one eukaryote (Saccharomyces cerevisiae). In this study we developed and assessed a distance-based approach, based on mean pairwise sequence similarity, for generating genome trees. Exclusion of PDSs improved bootstrap support for basal nodes but altered few topological features, indicating that there is little systematic bias among PDSs. Many but not all features of the genome tree from which PDSs were excluded are consistent with the 16S rRNA tree.

Distributional profiles of homologous open reading frames among bacterial phyla: implications for vertical and lateral transmission

If open reading frames (ORFs) have been transmitted primarily by vertical descent, the distributional profile of orthologues of each ORF should be congruent with the organismal tree or a subtree thereof. Distributional patterns not reconciled parsimoniously with tree-like descent and loss are prima facie evidence of lateral gene transfer. Herein, a rigorous criterion for recognizing ORF distributions is described and implemented; it does not require the inference of phylogenetic trees, nor does it assume any specific tree. Because lineage-specific differences in rates of sequence change can also generate unexpected distributional patterns, rate artefacts, were controlled for by requiring pairwise matches between ORFs to exceed a rigorous inclusion threshold, but absence of a match was assessed against a more-permissive exclusion threshold. Applying this dual-threshold criterion to cross-domain and cross-phylum distributional patterns for ORFs in 23 bacterial genomes, a relative abundance of ORFs was observed that find a match in exactly seven other bacterial phyla; 94-99% of these ORFs also find matches among the Archaea and/or Eukarya. In the larger (and some smaller) bacterial genomes, ORFs that find matches in exactly one other bacterial phylum are also relatively abundant, but fewer of these have non-bacterial homologues; most of their matches within the Bacteria are to the Proteobacteria and/or Firmicutes, which cannot be sister lineages to all bacteria. ORFs that are neither distributed universally among the Bacteria, nor necessarily shared with topologically adjacent lineages, are preferentially enriched in large bacterial genomes.

Analysis of a Local HIV-1 Epidemic in Portugal Highlights Established Transmission of non-B and -G Subtypes

Existing data supports Portugal as the Western Europe country with highest HIV-1 subtype diversity. However, detailed phylogenetic studies of Portuguese HIV-1 epidemics are still scarce. Thus, our main goal was to analyze the phylodynamics of a local HIV-1 infection in the Portuguese region of Minho. Molecular epidemiological analysis was applied to data from 289 HIV-1 infected individuals followed in the reference Hospital of the province of Minho, Portugal, in which isolated viruses had been sequenced between 2000 and 2012. Viruses of the G (29.1%) and B (27.0%) subtypes were the most frequent, followed by recombinant forms (17.6%), C (14.5%), F1 (7.3%) and A1 (4.2%) subtypes. Multinomial logistic regression revealed that the odds of being infected with A1 and F1 subtype increased over the years when compared with B, G, C or recombinant viruses. As expected, polyphyletic patterns suggesting multiple and old introductions of subtypes B and G were found. However, transmission clusters of non-B and -G viruses among native individuals were also found with the dates of the most recent common ancestor estimated to the early 2000s. Our study supports that the HIV-1 subtype diversity in the Portuguese region of Minho is high and has been increasing in a manner that is apparently driven by factors other than immigration and international travel. Infections with A1 and F1 viruses in the region of Minho are becoming established and were mainly found in sexually transmitted clusters, reinforcing the need for more efficacious control measures targeting this infection route.

A mega-analysis of genome-wide association studies for major depressive disorder.

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 × 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 × 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Búsqueda de ancentros comunes entre genomas de diferentes especies

La cerca de similituds entre regions de diferents genomes ofereix molta informació sobre les relaciones entre les especies d’aquest genomes. Es molt útil per a l’estudi de la conservació de gens d’una especia a un altre, de com les propietats d’un gen son assignades a un altre gen o de com es creen variacions en genomes diferents durant l’evolució d’aquestes especies. La finalitat d’aquest projecte es la creació d’una eina per a la cerca d’ancestres comuns de diferents especies basada en la comparació de la conservació entre regions dels genomes d’aquestes especies. Per a una comparació entre genomes mes eficaç una part important del projecte es destinarà a la creació d’una nova unitat de comparació. Aquestes noves unitats seran superestructures basades en agrupació dels MUMs existent per la mateixa comparació que anomenarem superMUMs. La aplicació final estarà disponible al servidor: http://revolutionresearch.uab.es

Molecular evidence that human immunodeficiency virus type 1 dissemination in a small Brazilian city was already taking place in the early 1990s

We recently performed a molecular epidemiology survey of human immunodeficiency virus type 1 (HIV-1) infection in Miracema, a small city in Southeast Brazil, and found multiple monophyletic clusters, consistent with independent introductions and spread of different viral lineages in the city. Here we apply Bayesian coalescent-based methods to the two largest subtype B clusters and estimate that the most recent common ancestors that gave rise to these two transmission chains were in circulation around 1991-1992. The finding that HIV-1 spread in this Brazilian small city was already taking place at a time Aids was considered a problem restricted to large urban centers may have important public health implications.

Comparative genomics of the odorant-binding and chemosensory protein gene families across the arthropoda: origin and evolutionary history of the chemosensory system

Chemoreception is a biological process essential for the survival of animals, as it allows the recognition of important volatile cues for the detection of food, egg-laying substrates, mates or predators, among other purposes. Furthermore, its role in pheromone detection may contribute to evolutionary processes such as reproductive isolation and speciation. This key role in several vital biological processes makes chemoreception a particularly interesting system for studying the role of natural selection in molecular adaptation. Two major gene families are involved in the perireceptor events of the chemosensory system: the odorant-binding protein (OBP) and chemosensory protein (CSP) families. Here, we have conducted an exhaustive comparative genomic analysis of these gene families in twenty Arthropoda species. We show that the evolution of the OBP and CSP gene families is highly dynamic, with a high number of gains and losses of genes, pseudogenes and independent origins of subfamilies. Taken together, our data clearly support the birth-and-death model for the evolution of these gene families with an overall high gene-turnover rate. Moreover, we show that the genome organization of the two families is significantly more clustered than expected by chance and, more important, that this pattern appears to be actively maintained across the Drosophila phylogeny. Finally, we suggest the homologous nature of the OBP and CSP gene families, dating back their MRCA (most recent common ancestor) to 380¿420 Mya, and we propose a scenario for the origin and diversification of these families.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardi/pulcher species complex endemic to Lake Tanganyika.

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past.