Decreasing the number of false positives in sequence classification

Abstract Background A large number of probabilistic models used in sequence analysis assign non-zero probability values to most input sequences. To decide when a given probability is sufficient the most common way is bayesian binary classification, where the probability of the model characterizing the sequence family of interest is compared to that of an alternative probability model. We can use as alternative model a null model. This is the scoring technique used by sequence analysis tools such as HMMER, SAM and INFERNAL. The most prevalent null models are position-independent residue distributions that include: the uniform distribution, genomic distribution, family-specific distribution and the target sequence distribution. This paper presents a study to evaluate the impact of the choice of a null model in the final result of classifications. In particular, we are interested in minimizing the number of false predictions in a classification. This is a crucial issue to reduce costs of biological validation. Results For all the tests, the target null model presented the lowest number of false positives, when using random sequences as a test. The study was performed in DNA sequences using GC content as the measure of content bias, but the results should be valid also for protein sequences. To broaden the application of the results, the study was performed using randomly generated sequences. Previous studies were performed on aminoacid sequences, using only one probabilistic model (HMM) and on a specific benchmark, and lack more general conclusions about the performance of null models. Finally, a benchmark test with P. falciparum confirmed these results. Conclusions Of the evaluated models the best suited for classification are the uniform model and the target model. However, the use of the uniform model presents a GC bias that can cause more false positives for candidate sequences with extreme compositional bias, a characteristic not described in previous studies. In these cases the target model is more dependable for biological validation due to its higher specificity.

Influence of Genomic Ancestry on the Distribution of SLCO1B1, SLCO1B3 and ABCB1 Gene Polymorphisms among Brazilians

The frequency distribution of SNPs and haplotypes in the ABCB1, SLCO1B1 and SLCO1B3 genes varies largely among continental populations. This variation can lead to biases in pharmacogenetic studies conducted in admixed populations such as those from Brazil and other Latin American countries. The aim of this study was to evaluate the influence of self-reported colour, geographical origin and genomic ancestry on distributions of the ABCB1, SLCO1B1 and SLCO1B3 polymorphisms and derived haplotypes in admixed Brazilian populations. A total of 1039 healthy adults from the north, north-east, south-east and south of Brazil were recruited for this investigation. The c.388A>G (rs2306283), c.463C>A (rs11045819) and c.521T>C (rs4149056) SNPs in the SLCO1B1 gene and c.334T>G (rs4149117) and c.699G>A (rs7311358) SNPs in the SLCO1B3 gene were determined by Taqman 5'-nuclease assays. The ABCB1 c.1236C>T (rs1128503), c.2677G>T/A (rs2032582) and c.3435C>T (rs1045642) polymorphisms were genotyped using a previously described single-base extension/termination method. The results showed that genotype and haplotype distributions are highly variable among populations of the same self-reported colour and geographical region. However, genomic ancestry showed that these associations are better explained by a continuous variable. The influence of ancestry on the distribution of alleles and haplotype frequencies was more evident in variants with large differences in allele frequencies between European and African populations. Design and interpretation of pharmacogenetic studies using these transporter genes should include genomic controls to avoid spurious conclusions based on improper matching of study cohorts from Brazilian populations and other highly admixed populations.

Distribution of human immunodeficiency virus type 1 subtypes in the State of Amazonas, Brazil, and subtype C identification

Few studies have reported the molecular epidemiological characterization of HIV-1 in the Northern region of Brazil. The present study reports the molecular and epidemiological characterization of 31 HIV-1 isolates from blood donors from the State of Amazonas who donated blood between April 2006 and March 2007. Serum/plasma samples from all donors were screened for HIV antibodies by ELISA and the results confirmed by Western blot analysis. Genomic DNA was extracted from the buffy coat using the Super Quik-Gene-DNA Isolation kit. Nested PCR was performed on the env, gag, and pol regions of HIV-1 using the Gene Amp PCR System 9700. Sequencing reactions were performed using the inner PCR primers and the DYEnamic (TM) ET Dye Terminator Kit, and phylogenetic analysis was performed using the gag, pol, and env gene sequences. We collected samples from 31 blood donors who tested positive for HIV-1 in confirmatory experiments. The male: female ratio of blood donors was 3.4:1, and the mean age was 32.4 years (range: 19 to 61 years). Phylogenetic analysis showed that subtype B is the most prevalent among Northern Brazilian HIV-1-seropositive blood donors. One HIV-1 subtype C and one circulating recombinant form (CRF_BF) of HIV-1 were identified in the State of Amazonas. This is the first study showing the occurrence of a possible "homogenous" subtype C in this region of Brazil. This finding could contribute to a better characterization of the HIV-1 strains that circulate in the country.

Novel insights into the genomic basis of citrus canker based on the genome sequences of two strains of Xanthomonas fuscans subsp. aurantifolii

Abstract Background Citrus canker is a disease that has severe economic impact on the citrus industry worldwide. There are three types of canker, called A, B, and C. The three types have different phenotypes and affect different citrus species. The causative agent for type A is Xanthomonas citri subsp. citri, whose genome sequence was made available in 2002. Xanthomonas fuscans subsp. aurantifolii strain B causes canker B and Xanthomonas fuscans subsp. aurantifolii strain C causes canker C. Results We have sequenced the genomes of strains B and C to draft status. We have compared their genomic content to X. citri subsp. citri and to other Xanthomonas genomes, with special emphasis on type III secreted effector repertoires. In addition to pthA, already known to be present in all three citrus canker strains, two additional effector genes, xopE3 and xopAI, are also present in all three strains and are both located on the same putative genomic island. These two effector genes, along with one other effector-like gene in the same region, are thus good candidates for being pathogenicity factors on citrus. Numerous gene content differences also exist between the three cankers strains, which can be correlated with their different virulence and host range. Particular attention was placed on the analysis of genes involved in biofilm formation and quorum sensing, type IV secretion, flagellum synthesis and motility, lipopolysacharide synthesis, and on the gene xacPNP, which codes for a natriuretic protein. Conclusion We have uncovered numerous commonalities and differences in gene content between the genomes of the pathogenic agents causing citrus canker A, B, and C and other Xanthomonas genomes. Molecular genetics can now be employed to determine the role of these genes in plant-microbe interactions. The gained knowledge will be instrumental for improving citrus canker control.

Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range

Abstract Background Citrus bacterial canker is a disease that has severe economic impact on citrus industries worldwide and is caused by a few species and pathotypes of Xanthomonas. X. citri subsp. citri strain 306 (XccA306) is a type A (Asiatic) strain with a wide host range, whereas its variant X. citri subsp. citri strain Aw12879 (Xcaw12879, Wellington strain) is restricted to Mexican lime. Results To characterize the mechanism for the differences in host range of XccA and Xcaw, the genome of Xcaw12879 that was completed recently was compared with XccA306 genome. Effectors xopAF and avrGf1 are present in Xcaw12879, but were absent in XccA306. AvrGf1 was shown previously for Xcaw to cause hypersensitive response in Duncan grapefruit. Mutation analysis of xopAF indicates that the gene contributes to Xcaw growth in Mexican lime but does not contribute to the limited host range of Xcaw. RNA-Seq analysis was conducted to compare the expression profiles of Xcaw12879 and XccA306 in Nutrient Broth (NB) medium and XVM2 medium, which induces hrp gene expression. Two hundred ninety two and 281 genes showed differential expression in XVM2 compared to in NB for XccA306 and Xcaw12879, respectively. Twenty-five type 3 secretion system genes were up-regulated in XVM2 for both XccA and Xcaw. Among the 4,370 common genes of Xcaw12879 compared to XccA306, 603 genes in NB and 450 genes in XVM2 conditions were differentially regulated. Xcaw12879 showed higher protease activity than XccA306 whereas Xcaw12879 showed lower pectate lyase activity in comparison to XccA306. Conclusions Comparative genomic analysis of XccA306 and Xcaw12879 identified strain specific genes. Our study indicated that AvrGf1 contributes to the host range limitation of Xcaw12879 whereas XopAF contributes to virulence. Transcriptome analyses of XccA306 and Xcaw12879 presented insights into the expression of the two closely related strains of X. citri subsp. citri. Virulence genes including genes encoding T3SS components and effectors are induced in XVM2 medium. Numerous genes with differential expression in Xcaw12879 and XccA306 were identified. This study provided the foundation to further characterize the mechanisms for virulence and host range of pathotypes of X. citri subsp. citri.

Endosymbiosis in trypanosomatids: the genomic cooperation between bacterium and host in the synthesis of essential amino acids is heavily influenced by multiple horizontal gene transfers

Background Trypanosomatids of the genera Angomonas and Strigomonas live in a mutualistic association characterized by extensive metabolic cooperation with obligate endosymbiotic Betaproteobacteria. However, the role played by the symbiont has been more guessed by indirect means than evidenced. Symbiont-harboring trypanosomatids, in contrast to their counterparts lacking symbionts, exhibit lower nutritional requirements and are autotrophic for essential amino acids. To evidence the symbiont’s contributions to this autotrophy, entire genomes of symbionts and trypanosomatids with and without symbionts were sequenced here. Results Analyses of the essential amino acid pathways revealed that most biosynthetic routes are in the symbiont genome. By contrast, the host trypanosomatid genome contains fewer genes, about half of which originated from different bacterial groups, perhaps only one of which (ornithine cyclodeaminase, EC:4.3.1.12) derived from the symbiont. Nutritional, enzymatic, and genomic data were jointly analyzed to construct an integrated view of essential amino acid metabolism in symbiont-harboring trypanosomatids. This comprehensive analysis showed perfect concordance among all these data, and revealed that the symbiont contains genes for enzymes that complete essential biosynthetic routes for the host amino acid production, thus explaining the low requirement for these elements in symbiont-harboring trypanosomatids. Phylogenetic analyses show that the cooperation between symbionts and their hosts is complemented by multiple horizontal gene transfers, from bacterial lineages to trypanosomatids, that occurred several times in the course of their evolution. Transfers occur preferentially in parts of the pathways that are missing from other eukaryotes. Conclusion We have herein uncovered the genetic and evolutionary bases of essential amino acid biosynthesis in several trypanosomatids with and without endosymbionts, explaining and complementing decades of experimental results. We uncovered the remarkable plasticity in essential amino acid biosynthesis pathway evolution in these protozoans, demonstrating heavy influence of horizontal gene transfer events, from Bacteria to trypanosomatid nuclei, in the evolution of these pathways.

Distribution of human immunodeficiency virus type 1 subtypes in the state of Amazonas, Brazil, and subtype C identification

Few studies have reported the molecular epidemiological characterization of HIV-1 in the Northern region of Brazil. The present study reports the molecular and epidemiological characterization of 31 HIV-1 isolates from blood donors from the State of Amazonas who donated blood between April 2006 and March 2007. Serum/plasma samples from all donors were screened for HIV antibodies by ELISA and the results confirmed by Western blot analysis. Genomic DNA was extracted from the buffy coat using the Super Quik-Gene-DNA Isolation kit. Nested PCR was performed on the env, gag, and pol regions of HIV-1 using the Gene Amp PCR System 9700. Sequencing reactions were performed using the inner PCR primers and the DYEnamic™ ET Dye Terminator Kit, and phylogenetic analysis was performed using the gag, pol, and env gene sequences. We collected samples from 31 blood donors who tested positive for HIV-1 in confirmatory experiments. The male:female ratio of blood donors was 3.4:1, and the mean age was 32.4 years (range: 19 to 61 years). Phylogenetic analysis showed that subtype B is the most prevalent among Northern Brazilian HIV-1-seropositive blood donors. One HIV-1 subtype C and one circulating recombinant form (CRF_BF) of HIV-1 were identified in the State of Amazonas. This is the first study showing the occurrence of a possible "homogenous" subtype C in this region of Brazil. This finding could contribute to a better characterization of the HIV-1 strains that circulate in the country.