A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Genotypic Characteristics of HIV Type 1 Based on gp120 Hypervariable Region 3 of Isolates from Southern Brazil

The aim of this study was to investigate HIV-1 molecular diversity and the epidemiological profile of HIV-1-infected patients from Ribeirao Preto, Brazil. A nested PCR followed by sequencing of a 302-base pair fragment of the env gene (C2-V3 region) was performed in samples from HIV-1-positive patients. A total of 45 sequences were aligned with final manual adjustments. The phylogenetic analyses showed a higher prevalence of HIV-1 subtype B in the studied population (97.8%) with only one sample yielding an F1 subtype. The viral genotyping prediction showed that CCR5 tropism was the most prevalent in the studied cohort. Geno2pheno analysis showed that R5 and CXCR4 prediction were 69% and 31%, respectively. There was no statistical significance, either in viral load or in CD4(+) T cell count when R5 and X4 prediction groups were compared. Moreover, the GPGR tetramer was the most common V3 loop core motif identified in the HIV-1 strains studied (34.1%) followed by GWGR, identified in 18.1% of the samples. The high level of B subtype in this Brazilian population reinforces the nature of the HIV epidemic in Brazil, and corroborates previous data obtained in the Brazilian HIV-infected population.

Genetic polymorphisms in TLR4, CR1 and Duffy genes are not associated with malaria resistance in patients from Baixo Amazonas region, Brazil

The main purpose of this research was to analyze the relation of the genetic polymorphisms frequently expressed by antigen-presenting cells, erythrocytes and malaria susceptibility/resistance with the human malaria infection cases. The sample used consisted of 23 Plasmodium vivax ( Pv)- and P. falciparum ( Pf)-infected patients, and 21 healthy individuals as a control group, from the Baixo Amazonas population in Para, Brazil. The Asp299Gly polymorphisms in the Toll-like receptor 4 ( TLR4), and Gly42Asp, Arg89Cys, Ala100Thr, and T-33C in the Duffy gene ( FY) were analyzed by restriction fragment length polymorphism-polymerase chain reaction. The Lys1590Glu and Arg1601Gly polymorphisms in the complement receptor type 1 (CR1) were analyzed by DNA sequencing. According to the results obtained and statistical analysis considering a significance level or alpha = 0.01, we conclude that the low heterozygote frequency (2.27%) for the Asp299Gly mutation, detected in the TLR4 gene, is not related to the Pv and Pf infections in the patients analyzed. Also, the promoter region GATA-1 analysis of the FY gene in the Pv-infected patients showed that the heterozygote frequency for the T-33C mutation (11.36% of the infected patients and 20.45% of the control patients) is not related to infection resistance. Regarding the CR1 gene, the observed heterozygote frequency (9.09%) for the Arg1601Gly mutation in Pf-infected patients when compared to heterozygote frequency in the control group (18.18%) suggests that there is no correlation with infection resistance.

Concurrent sequence variation of TP53 and TP73 genes in anaplastic astrocytoma

Disruption or loss of tumor suppressor gene TP53 is implicated in the development or progression of almost all different types of human malignancies. Other members of the p53 family have been identified. One member, p73, not only shares a high degree of similarity with p53 in its primary sequence, but also has similar functions. Like p53, p73 can bind to DNA and activate transcription. Using PCR-SSCP and gene sequencing, we analyzed the TP53 and TP73 genes in a case of a grade III anaplastic astrocytoma that progressed to glioblastoma. We found a deletion of AAG at position 595-597 of TP53 (exon 6), resulting in the deletion of Glu 199 in the protein and a genomic polymorphism of TP73, identified as an A-to-G change, at position E8/+15 at intron 8 (IVS8-15A>G). The mutation found at exon 6 of the gene TP53 could be associated with the rapid tumoral progression found in this case, since the mutated p53 may inactivate the wild-type p53 and the p73 alpha protein, which was conserved here, leading to an increase in cellular instability.

First Report of the Isolation and Molecular Characterization of Rickettsia amblyommii and Rickettsia felis in Central America

During 2010, 15 adult ticks, identified as Amblyomma cajennense, were collected from horses in Cahuita and Turrialba districts, whereas 7 fleas, identified as Ctenocephalides felis, were collected from a dog in San Jose city, Costa Rica. In the laboratory, three A. cajennense specimens, two from Cahuita and one from Turrialba, were individually processed for rickettsial isolation in cell culture, as was a pool of seven fleas. Rickettsiae were successfully isolated and established in Vero cell culture from the three ticks and from a pool of seven fleas in C6/36 cell culture. The three tick isolates were genotypically identified as Rickettsia amblyommii, and the flea isolate was identified as Rickettsia felis through DNA sequencing of portions of the rickettsial genes gltA, ompA, and ompB of each isolate. In addition, other seven ticks were shown to contain rickettsial DNA. Polymerase chain reaction products of at least two of these ticks were sequenced and also showed to correspond to R. amblyommii. Overall, 66.7% (10/15) of the A. cajennense adult ticks were found to be infected with rickettsiae. This is the first report of a successful isolation in cell culture of R. amblyommii and R. felis from Central America.

ISOLATION OF TOXOPLASMA GONDII FROM GOATS FROM BRAZIL

Goats are economically important in many countries, and little is known of caprine toxoplasmosis in Brazil. Anti-bodies to toxoplasma gondii were assayed in the sera of 143 goats from 3 Brazilian states, using modified agglutination test (MAT titer >= 1:25); 46 (32.2%) tested positive. Samples of brain, heart, diaphragm, and masseter of seropositive animals were pooled, digested in pepsin, and bioassayed in mice. Viable T. gondii specimens were isolated from tissue homogenates of 12 goats; the isolates were designated TgGtBr1-12. Ten of the 12 isolates killed 100% of infected mice, indicating that goats can harbor mouse-virulent T. gondii and, hence, can serve as a source of infection for humans.

SEROPREVALENCE AND ISOLATION OF TOXOPLASMA GONDII FROM SHEEP FROM SAO PAULO STATE, BRAZIL

Sheep are important in the epidemiology of Toxoplasma gondii infection, but little is known of ovine toxoplasmosis in Brazil. Antibodies to T. gondii were assayed in sera of 495 sheep from 36 countries of Sao Paulo state. Brazil, using the modified agglutination test (MAT titer >= 1:25); 120 (24.2%) sheep tested positive. Samples of brain, heart, and diaphragm of 85 seropositive sheep were pooled, digested in pepsin, and bioassayed in mice. Toxoplasma gondii was isolated from tissue homogenated of 16 sheep, and the isolated were designated TgShBr 1-16. Six of the 16 T. gondii isolated killed 100% of infected mice. Results indicate that asymptomatic sheep can harbour mouse-virulent T. gondii; hence, they can serve as a source of infection for humans.

Epigenetic Silencing of CRABP2 and MX1 in Head and Neck Tumors

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous disease affecting the epithelium of the oral cavity, pharynx and larynx. Conditions of most patients are diagnosed at late stages of the disease, and no sensitive and specific predictors of aggressive behavior have been identified yet. Therefore, early detection and prognostic biomarkers are highly desirable for a more rational management of the disease. Hypermethylation of CpG islands is one of the most important epigenetic mechanisms that leads to gene silencing in tumors and has been extensively used for the identification of biomarkers. In this study, we combined rapid subtractive hybridization and microarray analysis in a hierarchical manner to select genes that are putatively reactivated by the demethylating agent 5-aza-2'-deoxycytidine (5Aza-dC) in HNSCC cell lines (FaDu, UM-SCC-14A, UM-SCC-17A, UM-SCC-38A). This combined analysis identified 78 genes, 35 of which were reactivated in at least 2 cell lines and harbored a CpG island at their 5' region. Reactivation of 3 of these 35 genes (CRABP2, MX1, and SLC15A3) was confirmed by quantitative real-time polymerase chain reaction (PCR; fold change, >= 3). Bisulfite sequencing of their CpG islands revealed that they are indeed differentially methylated in the HNSCC cell lines. Using methylation-specific PCR, we detected a higher frequency of CRABP2 (58.1% for region 1) and MX1 (46.3%) hypermethylation in primary HNSCC when compared with lymphocytes from healthy individuals. Finally, absence of the CRABP2 protein was associated with decreased disease-free survival rates, supporting a potential use of CRABP2 expression as a prognostic biomarker for HNSCC patients.

A Genome-Wide Association Study of Upper Aerodigestive Tract Cancers Conducted within the INHANCE Consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p <= 5 x 10(-7)). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1 x 10(-8)) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p = 2 x 10(-8)) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5 x 10(-8); rs1229984-ADH1B, p = 7 x 10(-9); and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.

Characterization of new IS elements and studies of their dispersion in two subspecies of Leifsonia xyli

Background: Leifsonia xyli is a xylem-inhabiting bacterial species comprised of two subspecies: L. xyli subsp. xyli (Lxx) and L. xyli subsp. cynodontis (Lxc). Lxx is the causal agent of ratoon stunting disease in sugarcane commercial fields and Lxc colonizes the xylem of several grasses causing either mild or no symptoms of disease. The completely sequenced genome of Lxx provided insights into its biology and pathogenicity. Since IS elements are largely reported as an important source of bacterial genome diversification and nothing is known about their role in chromosome architecture of L. xyli, a comparative analysis of Lxc and Lxx elements was performed. Results: Sample sequencing of Lxc genome and comparative analysis with Lxx complete DNA sequence revealed a variable number of IS transposable elements acting upon genomic diversity. A detailed characterization of Lxc IS elements and a comparative review with IS elements of Lxx are presented. Each genome showed a unique set of elements although related to same IS families when considering features such as similarity among transposases, inverted and direct repeats, and element size. Most of the Lxc and Lxx IS families assigned were reported to maintain transposition at low levels using translation regulatory mechanisms, consistent with our in silico analysis. Some of the IS elements were found associated with rearrangements and specific regions of each genome. Differences were also found in the effect of IS elements upon insertion, although none of the elements were preferentially associated with gene disruption. A survey of transposases among genomes of Actinobacteria showed no correlation between phylogenetic relatedness and distribution of IS families. By using Southern hybridization, we suggested that diversification of Lxc isolates is also mediated by insertion sequences in probably recent events. Conclusion: Collectively our data indicate that transposable elements are involved in genome diversification of Lxc and Lxx. The IS elements were probably acquired after the divergence of the two subspecies and are associated with genome organization and gene contents. In addition to enhancing understanding of IS element dynamics in general, these data will contribute to our ongoing comparative analyses aimed at understanding the biological differences of the Lxc and Lxx.

Phylogenetic nomenclature and evolution of mannose-binding lectin (MBL2) haplotypes

Background: Polymorphisms of the mannose-binding lectin gene (MBL2) affect the concentration and functional efficiency of the protein. We recently used haplotype-specific sequencing to identify 23 MBL2 haplotypes, associated with enhanced susceptibility to several diseases. Results: In this work, we applied the same method in 288 and 470 chromosomes from Gabonese and European adults, respectively, and found three new haplotypes in the last group. We propose a phylogenetic nomenclature to standardize MBL2 studies and found two major phylogenetic branches due to six strongly linked polymorphisms associated with high MBL production. They presented high Fst values and were imbedded in regions with high nucleotide diversity and significant Tajima's D values. Compared to others using small sample sizes and unphased genotypic data, we found differences in haplotyping, frequency estimation, Fu and Li's D* and Fst results. Conclusion: Using extensive testing for selective neutrality, we confirmed that stochastic evolutionary factors have had a major role in shaping this polymorphic gene worldwide.

Endophytic and entomopathogenic strains of Beauveria sp to control the bovine tick Rhipicephalus (Boophilus) microplus

Pathogenicity of strains of the entomopathogenic fungus Beauveria bassiana and endophytic strains of Beauveria sp against the bovine tick Rhipicephalus (Boophilus) microplus was tested in laboratory bioassays and under field conditions. Suspensions containing 10(5), 10(7) and 10(9) conidia/mL were prepared of each fungal strain for laboratory bioassays. The ticks were maintained at 28 degrees C, 90 +/- 5% relative humidity, and the following variables were evaluated: initial female weight, egg weight, hatching percentage, reproductive efficiency, and percentage control. For tests under field conditions, a Beauveria suspension containing 10(6) conidia/mL was sprayed on tick-infested cows. After 72 h, the ticks were collected to estimate mortality under field conditions. Laboratory bioassays showed a mortality of 20 to 50% of the ticks seven days after inoculation with 10(7) Beauveria conidia/mL. Under field conditions 10(6) Beauveria conidia/mL induced 18-32% mortality. All Beauveria strains were effective in biological control of R. (Boophilus) microplus under laboratory and field test conditions. This is the first demonstration that endophytic fungi can be used for biological control of the cattle tick; this could help reduce environmental contamination by diminishing the need for chemical acaricides. Two endophytic strains were isolated from maize leaves and characterized by molecular sequencing of 5.8S rDNA ITS1 and ITS2 and morphological analyses of conidia. We found that these two endophytic Beauveria isolates, designated B95 and B157, are close to Beauveria amorpha.

Detection of polymorphisms of the mtDNA control region of Caretta caretta (Testudines: Cheloniidae) by PCR-SSCP

Marine turtles are increasingly being threatened worldwide by anthropogenic activities. Better understanding of their life cycle, behavior and population structure is imperative for the design of adequate conservation strategies. The mtDNA control region is a fast-evolving matrilineal marker that has been employed in the study of marine turtle populations. We developed and tested a simple molecular tracing system for Caretta caretta mtDNA haplotypes by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Using this technique, we were able to distinguish the SSCP patterns of 18 individuals of the haplotypes CC-A4, CC-A24 and CCxLO, which are commonly found in turtles sampled on the Brazilian coast. When we analyzed 15 turtles with previously unknown sequences, we detected two other haplotypes, in addition to the other four. Based on DNA sequencing, they were identified as the CC-A17 and CC-A1 haplotypes. Further analyses were made with the sea turtles, Chelonia mydas (N = 8), Lepidochelys olivacea (N = 3) and Eretmochelys imbricata (N = 1), demonstrating that the PCR-SSCP technique is able to distinguish intra-and interspecific variation in the family Cheloniidae. We found that this technique can be useful for identifying sea turtle mtDNA haplotypes, reducing the need for sequencing.

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.

High Prevalence of bla(CTX-M) Extended Spectrum Beta-Lactamase Genes in Klebsiella pneumoniae Isolates from a Tertiary Care Hospital: First report of bla(SHV-12), bla(SHV-31), bla(SHV-38), and bla(CTX-M-15) in Brazil

The aim of this study was to investigate the presence and prevalence of bla(TEM), bla(SHV), and bla(CTX-M) and bla(GES)-like genes, responsible for extended spectrum beta-lactamases (ESBLs) production in clinical isolates of Klebsiella pneumoniae collected from a Brazilian tertiary care hospital. Sixty-five ESBL producing K. pneumoniae isolates, collected between 2005 and 2007, were screened by polymerase chain reaction (PCR). Identification of bla genes was achieved by sequencing. Genotyping of ESBL producing K. pneumoniae was performed by the enterobacterial repetitive intergenic consensus-PCR with cluster analysis by the Dice coefficient. The presence of genes encoding ESBLs was confirmed in 59/65 (90.8%) isolates, comprising 20 bla(CTX-M-2), 14 bla(CTX-M-59), 12 bla(CTX-M-15), 9 bla(SHV-12), 1 bla(SHV-2), 1 bla(SHV-2a), 1 bla(SHV-5), and 1 bla(SHV-31) genes. The ESBL genes bla(SHV-12), bla(SHV-31), and bla(CTX-M-15), and the chromosome-encoded SHV-type beta-lactamase capable of hydrolyzing imipenem were detected in Brazil for the first time. The analysis of the enterobacterial repetitive intergenic consensus-PCR band patterns revealed a high rate of multiclonal bla(CTX-M) carrying K. pneumoniae isolates (70.8%), suggesting that dissemination of encoding plasmids is likely to be the major cause of the high prevalence of these genes among the K. pneumoniae isolates considered in this study.