Determination and Molecular Analysis of the Complete Genome Sequence of Two Wild-Type Rabies Viruses Isolated from a Haematophagous Bat and a Frugivorous Bat in Brazil

The complete genome sequences of two Brazilian wild-type rabies viruses (RABV), a BR-DR1 isolate from a haematophagous bat (Desmodus rotundus) and a BR-AL1 isolate from a frugivorous bat (Artibeus lituratus), were determined. The genomes of the BR-DR1 and RR-AL1 had 11,923 and 11,922 nt, respectively, and both encoded the five standard genes of rhabdoviruses. The complete nucleotide sequence identity between the BR-DR1 and BR-AL1 isolates was 97%. The BR-DR1 and BR-AL1 isolates had some conserved functional sites revealed by the fixed isolates, whereas both isolates had unique amino acid substitutions in the antigenic region IV of the nucleocapsid gene. Therefore, it is speculated that both isolates were nearly identical in virologic character. According to our phylogenetic analysis based on the complete genomes, both isolates belonged to genotype 1, and to the previously defined ""vampire bat-related RABV lineage"" which consisted of mainly D. rotundus- and A. lituratus- isolates; however, a branch pattern with high bootstrap values suggested that BR-DR1 was more closely related to the 9001FRA isolate, which was collected from a dog bitten by a bat in French Guiana, than to BR-AL1. This result suggests that the vampire bat-related RABV lineage includes Brazilian vampire bat and Brazilian frugivorous bat RABV and is further divided into Brazilian vampire bat and Brazilian frugivorous bat RABV sub-lineages. The phylogenetic analysis based on the complete genomes was valuable in discriminating among very closely related isolates.

Complete Genome Sequence of Mycoplasma suis and Insights into Its Biology and Adaption to an Erythrocyte Niche

Mycoplasma suis, the causative agent of porcine infectious anemia, has never been cultured in vitro and mechanisms by which it causes disease are poorly understood. Thus, the objective herein was to use whole genome sequencing and analysis of M. suis to define pathogenicity mechanisms and biochemical pathways. M. suis was harvested from the blood of an experimentally infected pig. Following DNA extraction and construction of a paired end library, whole-genome sequencing was performed using GS-FLX (454) and Titanium chemistry. Reads on paired-end constructs were assembled using GS De Novo Assembler and gaps closed by primer walking; assembly was validated by PFGE. Glimmer and Manatee Annotation Engine were used to predict and annotate protein-coding sequences (CDS). The M. suis genome consists of a single, 742,431 bp chromosome with low G+C content of 31.1%. A total of 844 CDS, 3 single copies, unlinked rRNA genes and 32 tRNAs were identified. Gene homologies and GC skew graph show that M. suis has a typical Mollicutes oriC. The predicted metabolic pathway is concise, showing evidence of adaptation to blood environment. M. suis is a glycolytic species, obtaining energy through sugars fermentation and ATP-synthase. The pentose-phosphate pathway, metabolism of cofactors and vitamins, pyruvate dehydrogenase and NAD(+) kinase are missing. Thus, ribose, NADH, NADPH and coenzyme A are possibly essential for its growth. M. suis can generate purines from hypoxanthine, which is secreted by RBCs, and cytidine nucleotides from uracil. Toxins orthologs were not identified. We suggest that M. suis may cause disease by scavenging and competing for host nutrients, leading to decreased life-span of RBCs. In summary, genome analysis shows that M. suis is dependent on host cell metabolism and this characteristic is likely to be linked to its pathogenicity. The prediction of essential nutrients will aid the development of in vitro cultivation systems.

Complete Genome Viral Phylogenies Suggests the Concerted Evolution of Regulatory Cores and Accessory Satellites

We consider the concerted evolution of viral genomes in four families of DNA viruses. Given the high rate of horizontal gene transfer among viruses and their hosts, it is an open question as to how representative particular genes are of the evolutionary history of the complete genome. To address the concerted evolution of viral genes, we compared genomic evolution across four distinct, extant viral families. For all four viral families we constructed DNA-dependent DNA polymerase-based (DdDp) phylogenies and in addition, whole genome sequence, as quantitative descriptions of inter-genome relationships. We found that the history of the polymerase gene was highly predictive of the history of the genome as a whole, which we explain in terms of repeated, co-divergence events of the core DdDp gene accompanied by a number of satellite, accessory genetic loci. We also found that the rate of gene gain in baculovirus and poxviruses proceeds significantly more quickly than the rate of gene loss and that there is convergent acquisition of satellite functions promoting contextual adaptation when distinct viral families infect related hosts. The congruence of the genome and polymerase trees suggests that a large set of viral genes, including polymerase, derive from a phylogenetically conserved core of genes of host origin, secondarily reinforced by gene acquisition from common hosts or co-infecting viruses within the host. A single viral genome can be thought of as a mutualistic network, with the core genes acting as an effective host and the satellite genes as effective symbionts. Larger virus genomes show a greater departure from linkage equilibrium between core and satellites functions.

Phylogenetic Analysis of Hepatitis B Virus Genotype F Complete Genome Sequences From Chilean Patients With Chronic Infection

Molecular epidemiological data concerning the hepatitis B virus (HBV) in Chile are not known completely. Since the HBV genotype F is the most prevalent in the country, the goal of this study was to obtain full HBV genome sequences from patients infected chronically in order to determine their subgenotypes and the occurrence of resistance-associated mutations. Twenty-one serum samples from antiviral drug-naive patients with chronic hepatitis B were subjected to full-length PCR amplification, and both strands of the whole genomes were fully sequenced. Phylogenetic analyses were performed along with reference sequences available from GenBank (n = 290). The sequences were aligned using Clustal X and edited in the SE-AL software. Bayesian phylogenetic analyses were conducted by Markov Chain Monte Carlo simulations (MCMC) for 10 million generations in order to obtain the substitution tree using BEAST. The sequences were also analyzed for the presence of primary drug resistance mutations using CodonCode Aligner Software. The phylogenetic analyses indicated that all sequences were found to be the HBV subgenotype F1b, clustered into four different groups, suggesting that diverse lineages of this subgenotype may be circulating within this population of Chilean patients. J. Med. Virol. 83: 1530-1536, 2011. (C) 2011 Wiley-Liss, Inc.

Complete genome analysis of a rabies virus isolate from Brazilian wild fox

The complete genome sequence of wild-type rabies virus (RABV) isolated from a wild Brazilian hoary fox (Dusicyon sp.), the BR-Pfx1 isolate, was determined and compared with fixed RABV strains. The genome structure and organization of the BR-Pfx1 isolate were composed of 11,924 nt and included the five standard genes of rhabdoviruses. Sequences of mRNA start and stop signals for transcription were highly conserved among all structural protein genes of the BR-Pfx1 isolate. All amino acid residues in the glycoprotein (G) gene associated with pathogenicity were retained in the BR-Pfx1 isolate, while unique amino acid substitutions were found in antigenic region I of the nucleoprotein gene and III of G. These results suggest that although the standard genome structure and organization of the RABV isolate are common between the BR-Pfx1 isolate and fixed RABV strains, the unique amino acid substitutions in functional sites of the BR-Pfx1 isolate may result in different biological characteristics from fixed RABV strains.

Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production

Bioethanol is a biofuel produced mainly from the fermentation of carbohydrates derived from agricultural feedstocks by the yeast Saccharomyces cerevisiae. One of the most widely adopted strains is PE-2, a heterothallic diploid naturally adapted to the sugar cane fermentation process used in Brazil. Here we report the molecular genetic analysis of a PE-2 derived diploid (JAY270), and the complete genome sequence of a haploid derivative (JAY291). The JAY270 genome is highly heterozygous (similar to 2 SNPs/kb) and has several structural polymorphisms between homologous chromosomes. These chromosomal rearrangements are confined to the peripheral regions of the chromosomes, with breakpoints within repetitive DNA sequences. Despite its complex karyotype, this diploid, when sporulated, had a high frequency of viable spores. Hybrid diploids formed by outcrossing with the laboratory strain S288c also displayed good spore viability. Thus, the rearrangements that exist near the ends of chromosomes do not impair meiosis, as they do not span regions that contain essential genes. This observation is consistent with a model in which the peripheral regions of chromosomes represent plastic domains of the genome that are free to recombine ectopically and experiment with alternative structures. We also explored features of the JAY270 and JAY291 genomes that help explain their high adaptation to industrial environments, exhibiting desirable phenotypes such as high ethanol and cell mass production and high temperature and oxidative stress tolerance. The genomic manipulation of such strains could enable the creation of a new generation of industrial organisms, ideally suited for use as delivery vehicles for future bioenergy technologies.

Characterization of Hepatitis B virus (HBV) genotypes in patients from Rondonia, Brazil

Background: Hepatitis B virus (HBV) can be classified into nine genotypes (A-I) defined by sequence divergence of more than 8% based on the complete genome. This study aims to identify the genotypic distribution of HBV in 40 HBsAg-positive patients from Rondonia, Brazil. A fragment of 1306 bp partially comprising surface and polymerase overlapping genes was amplified by PCR. Amplified DNA was purified and sequenced. Amplified DNA was purified and sequenced on an ABI PRISM (R) 377 Automatic Sequencer (Applied Biosystems, Foster City, CA, USA). The obtained sequences were aligned with reference sequences obtained from the GenBank using Clustal X software and then edited with Se-Al software. Phylogenetic analyses were conducted by the Markov Chain Monte Carlo (MCMC) approach using BEAST v.1.5.3. Results: The subgenotypes distribution was A1 (37.1%), D3 (22.8%), F2a (20.0%), D4 (17.1%) and D2 (2.8%). Conclusions: These results for the first HBV genotypic characterization in Rondonia state are consistent with other studies in Brazil, showing the presence of several HBV genotypes that reflects the mixed origin of the population, involving descendants from Native Americans, Europeans, and Africans.

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.

A Component of the Xanthomonadaceae Type IV Secretion System Combines a VirB7 Motif with a N0 Domain Found in Outer Membrane Transport Proteins

Type IV secretion systems (T4SS) are used by Gram-negative bacteria to translocate protein and DNA substrates across the cell envelope and into target cells. Translocation across the outer membrane is achieved via a ringed tetradecameric outer membrane complex made up of a small VirB7 lipoprotein (normally 30 to 45 residues in the mature form) and the C-terminal domains of the VirB9 and VirB10 subunits. Several species from the genera of Xanthomonas phytopathogens possess an uncharacterized type IV secretion system with some distinguishing features, one of which is an unusually large VirB7 subunit (118 residues in the mature form). Here, we report the NMR and 1.0 angstrom X-ray structures of the VirB7 subunit from Xanthomonas citri subsp. citri (VirB7(XAC2622)) and its interaction with VirB9. NMR solution studies show that residues 27-41 of the disordered flexible N-terminal region of VirB7(XAC2622) interact specifically with the VirB9 C-terminal domain, resulting in a significant reduction in the conformational freedom of both regions. VirB7(XAC2622) has a unique C-terminal domain whose topology is strikingly similar to that of N0 domains found in proteins from different systems involved in transport across the bacterial outer membrane. We show that VirB7(XAC2622) oligomerizes through interactions involving conserved residues in the N0 domain and residues 42-49 within the flexible N-terminal region and that these homotropic interactions can persist in the presence of heterotropic interactions with VirB9. Finally, we propose that VirB(7XAC2622) oligomerization is compatible with the core complex structure in a manner such that the N0 domains form an extra layer on the perimeter of the tetradecameric ring.

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.

European 1: A globally important clonal complex of Mycobacterium bovis

We have identified a globally important clonal complex of Mycobacterium bovis by deletion analysis of over one thousand strains from over 30 countries. We initially show that over 99% of the strains of M. bovis, the cause of bovine tuberculosis, isolated from cattle in the Republic of Ireland and the UK are closely related and are members of a single clonal complex marked by the deletion of chromosomal region RDEu1 and we named this clonal complex European 1 (Eu1). Eu1 strains were present at less than 14% of French, Portuguese and Spanish isolates of M. bovis but are rare in other mainland European countries and Iran. However, strains of the Eu1 clonal complex were found at high frequency in former trading partners of the UK (USA, South Africa, New Zealand, Australia and Canada). The Americas, with the exception of Brazil, are dominated by the Eu1 clonal complex which was at high frequency in Argentina, Chile, Ecuador and Mexico as well as North America. Eu1 was rare or absent in the African countries surveyed except South Africa. A small sample of strains from Taiwan were non-Eu1 but, surprisingly, isolates from Korea and Kazakhstan were members of the Eu1 clonal complex. The simplest explanation for much of the current distribution of the Eu1 clonal complex is that it was spread in infected cattle, such as Herefords, from the UK to former trading partners, although there is evidence of secondary dispersion since. This is the first identification of a globally dispersed clonal complex M. bovis and indicates that much of the current global distribution of this important veterinary pathogen has resulted from relatively recent International trade in cattle. Crown Copyright (C) 2011 Published by Elsevier B.V. All rights reserved.

alpha-Hydroxynitrile lyase protein from Xylella fastidiosa: Cloning, expression, and characterization

Xylella fastidiosa is a xylem-restricted plant pathogen that causes a range of diseases in several and important crops. Through comparative genomic sequence analysis many genes were identified and, among them, several potentially involved in plant-pathogen interaction. The experimental determination of the primary sequence of some markedly expressed proteins for X fastidiosa and the comparison with the nucleic acids sequence of genome identified one of them as being SCJ21.16 (XFa0032) gene product. The comparative analysis of this protein against SWISSPROT database, in special, resulted in similarity with a-hydroxynitrile lyase enzyme (HNL) from Arabidopsis thaliana, causing interest for being one of the most abundant proteins both in the whole cell extract as well as in the extracellular protein fraction. It is known that HNL enzyme are involved in a process termed ""cyanogenesis"", which catalyzes the dissociation of alpha-hydroxinitrile into carbonyle and HCN when plant tissue is damaged. Although the complete genome sequences of X.fastidiosa are available and the cyanogenesis process is well known, the biological role of this protein in this organism is not yet functionally characterized. In this study we presented the cloning, expression, characterization of recombinant HNL from X fastidiosa, and its probable function in the cellular metabolism. The successful cloning and heterologous expression in Escherichia coli resulted in a satisfactory amount of the recombinant HNL expressed in a soluble, and active form giving convenient access to pure enzyme for biochemical and structural studies. Finally, our results confirmed that the product of the gene XFa0032 can be positively assigned as FAD-independent HNLs. (C) 2009 Elsevier Ltd. All rights reserved.

An oligonucleotide primer set for PCR amplification of the complete honey bee mitochondrial genome

Mitochondrial DNA markers have been widely used to address population and evolutionary questions in the honey bee Apis mellifera. Most of the polymorphic markers are restricted to few mitochondrial regions. Here we describe a set of 24 oligonucleotides that allow PCR amplification of the entire mitochondrial genome of the honey bee A. mellifera in 12 amplicons. These fragments have important applications for the study of mitochondrial genes in different subspecies of A. mellifera and as heterospecific probes to characterize mitochondrial genomes in other bee species.

The mitochondrial genome of the stingless bee Melipona bicolor (Hymenoptera, Apidae, Meliponini): sequence, gene organization and a unique tRNA translocation event conserved across the tribe Meliponini

At present a complete mtDNA sequence has been reported for only two hymenopterans, the Old World honey bee, Apis mellifera and the sawfly Perga condei. Among the bee group, the tribe Meliponini (stingless bees) has some distinction due to its Pantropical distribution, great number of species and large importance as main pollinators in several ecosystems, including the Brazilian rain forest. However few molecular studies have been conducted on this group of bees and few sequence data from mitochondrial genomes have been described. In this project, we PCR amplified and sequenced 78% of the mitochondrial genome of the stingless bee Melipona bicolor (Apidae, Meliponini). The sequenced region contains all of the 13 mitochondrial protein-coding genes, 18 of 22 tRNA genes, and both rRNA genes (one of them was partially sequenced). We also report the genome organization (gene content and order), gene translation, genetic code, and other molecular features, such as base frequencies, codon usage, gene initiation and termination. We compare these characteristics of M. bicolor to those of the mitochondrial genome of A. mellifera and other insects. A highly biased A+T content is a typical characteristic of the A. mellifera mitochondrial genome and it was even more extreme in that of M. bicolor. Length and compositional differences between M. bicolor and A. mellifera genes were detected and the gene order was compared. Eleven tRNA gene translocations were observed between these two species. This latter finding was surprising, considering the taxonomic proximity of these two bee tribes. The tRNA Lys gene translocation was investigated within Meliponini and showed high conservation across the Pantropical range of the tribe.

An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.