Mastering seeds for genomic size nucleotide BLAST searches

One of the most common activities in bioinformatics is the search for similar sequences. These searches are usually carried out with the help of programs from the NCBI BLAST family. As the majority of searches are routinely performed with default parameters, a question that should be addressed is how reliable the results obtained using the default parameter values are, i.e. what fraction of potential matches have been retrieved by these searches. Our primary focus is on the initial hit parameter, also known as the seed or word, used by the NCBI BLASTn, MegaBLAST and other similar programs in searches for similar nucleotide sequences. We show that the use of default values for the initial hit parameter can have a big negative impact on the proportion of potentially similar sequences that are retrieved. We also show how the hit probability of different seeds varies with the minimum length and similarity of sequences desired to be retrieved and describe methods that help in determining appropriate seeds. The experimental results described in this paper illustrate situations in which these methods are most applicable and also show the relationship between the various BLAST parameters.     

Diversity of Helicobacter pylori cagA and vacA genes in Costa Rica: its relationship with atrophic gastritis and gastric cancer

BACKGROUND: Associations between Helicobacter pylori gene diversity and gastric cancer have not been reported on in Costa Rica, despite its being one of the countries with the highest gastric cancer incidence and mortality rates in the world. The aim of this study was to determine the prevalence of H. pylori cagA and vacA genes and investigate whether it could be correlated with atrophic gastritis (AG) and gastric cancer (GC) in Costa Rica. MATERIALS AND METHODS: Genomic DNAs from isolates of 104 patients classified into two groups: non-atrophic gastritis group (n = 68) and atrophic gastritis group (n = 36), were subjected to PCR-based genotyping of cagA and vacA genes and their correlation with clinical outcome was investigated. Total DNA extractions from gastric tissues of 25 H. pylori-infected gastric cancer patients were utilized for comparative purposes. RESULTS: The presence of cagA (75.3%), vacA s1b (75.3%), and vacA m1 (74.2%) was detected, and colonization by strains with different vacA genotypes in the same stomach was found in 9.7% of the patients. Age- and sex-adjusted vacA s1b and vacA m1 were associated with GC while only vacA m1 was significantly associated with AG. A tendency for association between cagA and vacA s1b, and AG was reported. CONCLUSIONS: The prevalence status of the cagA and vacA (s1/m1) genes in Costa Rica seems to fall between that found in European/North American and East Asian countries, and both cagA and vacA seem to have clinical relevance in this country.     

Single-cell analysis of yeast, mammalian cells, and fungal spores with a microfluidic pressure-driven chip-based system.

BACKGROUND: Cytomics aims at understanding the function of cellular systems by analysis of single cells. Recently, there has been a growing interest in single cell measurements being performed in microfluidic systems. These systems promise to integrate staining, measurement, and analysis in a single system. One important aspect is the limitation of allowable cell sizes due to microfluidic channel dimensions. Here we want to demonstrate the broad applicability of microfluidic chip technology for the analysis of many different cell types. METHODS: We have developed a microfluidic chip and measurement system that allows flow cytometric analysis of fluorescently stained cells from different organisms. In this setup, the cells are moved by pressure-driven flow inside a network of microfluidic channels and are analyzed individually by fluorescence detection. RESULTS: We have successfully applied the system to develop a methodology to detect viable and dead cells in yeast cell populations. Also, we have measured short interfering RNA (siRNA) mediated silencing of protein expression in mammalian cells. In addition, we have characterized the infection state of Magnaportae grisea fungal spores. CONCLUSIONS: Results obtained with the microfluidic system demonstrate a broad applicability of microfluidic flow cytometry to measurement of various cell types.     

Expression of Hbs,Kirre, and Rst during Drosophila ovarian development

The Irre cell‐recognition module (IRM) is a group of evolutionarily conserved and structurally related transmembrane glycoproteins of the immunoglobulin superfamily. In Drosophila melanogaster, it comprises the products of the genes roughest (rst; also known as irreC‐rst), kin‐of‐irre (kirre; also known as duf), sticks‐and‐stones (sns), and hibris (hbs). In this model organism, the behavior of this group of proteins as a partly redundant functional unit mediating selective cell recognition was demonstrated in a variety of developmental contexts, but their possible involvement in ovarian development and oogenesis has not been investigated, notwithstanding the fact that some rst mutant alleles are also female sterile. Here, we show that IRM genes are dynamically and, to some extent, coordinately transcribed in both pupal and adult ovaries. Additionally, the spatial distribution of Hbs, Kirre, and Rst proteins indicates that they perform cooperative, although largely nonredundant, functions. Finally, phenotypical characterization of three different female sterile rst alleles uncovered two temporally separated and functionally distinct requirements for this locus in ovarian development: one in pupa, essential for the organization of peritoneal and epithelial sheaths that maintain the structural integrity of the adult organ and another, in mature ovarioles, needed for the progression of oogenesis beyond stage 10.     

Rho GTPase expression in tumourigenesis: evidence for a significant link

Rho proteins belong to the small GTPases superfamily. They function as molecular switches that, in response to diverse stimuli, control key signaling and structural aspects of the cell. Although early studies proposed a role for Rho GTPases in cellular transformation, this effect was underestimated due to the fact that no genetic mutations affecting Rho-encoding genes were found in tumors. Recently, it has become evident that Rho GTPases participate in the carcinogenic process by either overexpression of some of the members of the family with oncogenic activity, downmodulation of other members with suggested tumor suppressor activity, or by alteration of upstream modulators or downstream effectors. Thus, alteration of the levels of expression of different members of the family of Rho GTPases has been detected in many types of human tumors leading to a great interest in the cellular effects elicited by these oncoproteins. This essay reviews the current evidence of dysregulation of Rho signaling by overexpression in human tumors.     

Comparative genome-wide DNA methylation analysis of colorectal tumor and matched normal tissues

Aberrant DNA methylation often occurs in colorectal cancer (CRC). In our study we applied a genome-wide DNA methylation analysis approach, MethylCap-seq, to map the differentially methylated regions (DMRs) in 24 tumors and matched normal colon samples. In total, 2687 frequently hypermethylated and 468 frequently hypomethylated regions were identified, which include potential biomarkers for CRC diagnosis. Hypermethylation in the tumor samples was enriched at CpG islands and gene promoters, while hypomethylation was distributed throughout the genome. Using epigenetic data from human embryonic stem cells, we show that frequently hypermethylated regions coincide with bivalent loci in human embryonic stem cells. DNA methylation is commonly thought to lead to gene silencing; however, integration of publically available gene expression data indicates that 75% of the frequently hypermethylated genes were most likely already lowly or not expressed in normal tissue. Collectively, our study provides genome-wide DNA methylation maps of CRC, comprehensive lists of DMRs, and gives insights into the role of aberrant DNA methylation in CRC formation.