Simcluster: clustering enumeration gene expression data on the simplex space

Abstract Background Transcript enumeration methods such as SAGE, MPSS, and sequencing-by-synthesis EST "digital northern", are important high-throughput techniques for digital gene expression measurement. As other counting or voting processes, these measurements constitute compositional data exhibiting properties particular to the simplex space where the summation of the components is constrained. These properties are not present on regular Euclidean spaces, on which hybridization-based microarray data is often modeled. Therefore, pattern recognition methods commonly used for microarray data analysis may be non-informative for the data generated by transcript enumeration techniques since they ignore certain fundamental properties of this space. Results Here we present a software tool, Simcluster, designed to perform clustering analysis for data on the simplex space. We present Simcluster as a stand-alone command-line C package and as a user-friendly on-line tool. Both versions are available at: http://xerad.systemsbiology.net/simcluster. Conclusion Simcluster is designed in accordance with a well-established mathematical framework for compositional data analysis, which provides principled procedures for dealing with the simplex space, and is thus applicable in a number of contexts, including enumeration-based gene expression data.

Comparative performance evaluation of hepatitis C virus genotyping based on the 5' untranslated region versus partial sequencing of the NS5B region of brazilian patients with chronic hepatitis C

Background Genotyping of hepatitis C virus (HCV) has become an essential tool for prognosis and prediction of treatment duration. The aim of this study was to compare two HCV genotyping methods: reverse hybridization line probe assay (LiPA v.1) and partial sequencing of the NS5B region. Methods Plasma of 171 patients with chronic hepatitis C were screened using both a commercial method (LiPA HCV Versant, Siemens, Tarrytown, NY, USA) and different primers targeting the NS5B region for PCR amplification and sequencing analysis. Results Comparison of the HCV genotyping methods showed no difference in the classification at the genotype level. However, a total of 82/171 samples (47.9%) including misclassification, non-subtypable, discrepant and inconclusive results were not classified by LiPA at the subtype level but could be discriminated by NS5B sequencing. Of these samples, 34 samples of genotype 1a and 6 samples of genotype 1b were classified at the subtype level using sequencing of NS5B. Conclusions Sequence analysis of NS5B for genotyping HCV provides precise genotype and subtype identification and an accurate epidemiological representation of circulating viral strains.

A note on the use of the generalized odds ratio in meta-analysis of association studies involving bi- and tri-allelic polymorphisms

Abstract Background The generalized odds ratio (GOR) was recently suggested as a genetic model-free measure for association studies. However, its properties were not extensively investigated. We used Monte Carlo simulations to investigate type-I error rates, power and bias in both effect size and between-study variance estimates of meta-analyses using the GOR as a summary effect, and compared these results to those obtained by usual approaches of model specification. We further applied the GOR in a real meta-analysis of three genome-wide association studies in Alzheimer's disease. Findings For bi-allelic polymorphisms, the GOR performs virtually identical to a standard multiplicative model of analysis (e.g. per-allele odds ratio) for variants acting multiplicatively, but augments slightly the power to detect variants with a dominant mode of action, while reducing the probability to detect recessive variants. Although there were differences among the GOR and usual approaches in terms of bias and type-I error rates, both simulation- and real data-based results provided little indication that these differences will be substantial in practice for meta-analyses involving bi-allelic polymorphisms. However, the use of the GOR may be slightly more powerful for the synthesis of data from tri-allelic variants, particularly when susceptibility alleles are less common in the populations (≤10%). This gain in power may depend on knowledge of the direction of the effects. Conclusions For the synthesis of data from bi-allelic variants, the GOR may be regarded as a multiplicative-like model of analysis. The use of the GOR may be slightly more powerful in the tri-allelic case, particularly when susceptibility alleles are less common in the populations.

On the morphological differentiation between Libinia spinosa and L. ferreirae (Crustacea: Brachyura: Majoidea: Epialtidae)

Libinia spinosa H. Milne Edwards in Guérin, 1832 and L. ferreirae Brito Capello, 1871, inhabit very similar environments, and their geographic and bathymetric distributions overlap for about 3000 km along the southwestern Atlantic. Both species are commonly caught in the same haul and differentiating between them can often be difficult. Traditionally, morphological differentiation between L. spinosa and L. ferreirae has been based exclusively on the number of spines along the median, longitudinal line of the carapace and the development of a process at the anterolateral angle of the basal segment of the antenna. Because Libinia spinosa and L. ferreirae share similar numbers of median spines (7 and 6, respectively), and the number of median spines of the carapace and the process at the anterolateral angle of the basal antennal segment are variable, they are of little value in separating these species. It is shown herein that unequivocal identification can be easily achieved based on features of the male and female thoracic sternum, pereiopod dactyli, and infraorbital notch. A lectotype is designated for L. spinosa and its authorship and date are corrected. Libinia gibbosa A. Milne-Edwards, 1878, is demonstrated to be a junior synonym of L. ferreirae. The holotype of L. gibbosa is figured for the first time.

Bifurcation effects on the plasma edge of tokamaks with ergodic limiter.

In the present paper, we solve a twist symplectic map for the action of an ergodic magnetic limiter in a large aspect-ratio tokamak. In this model, we study the bifurcation scenarios that occur in the remnants regular islands that co-exist with chaotic magnetic surfaces. The onset of atypical local bifurcations created by secondary shearless tori are identified through numerical profiles of internal rotation number and we observe that their rupture can reduce the usual magnetic field line escape at the tokamak plasma edge.

On the three-finger protein domain fold and CD59-like proteins in Schistosoma mansoni

Background: It is believed that schistosomes evade complement-mediated killing by expressing regulatory proteins on their surface. Recently, six homologues of human CD59, an important inhibitor of the complement system membrane attack complex, were identified in the schistosome genome. Therefore, it is important to investigate whether these molecules could act as CD59-like complement inhibitors in schistosomes as part of an immune evasion strategy. Methodology/Principal Findings: Herein, we describe the molecular characterization of seven putative SmCD59-like genes and attempt to address the putative biological function of two isoforms. Superimposition analysis of the 3D structure of hCD59 and schistosome sequences revealed that they contain the three-fingered protein domain (TFPD). However, the conserved amino acid residues involved in complement recognition in mammals could not be identified. Real-time RT-PCR and Western blot analysis determined that most of these genes are up-regulated in the transition from free-living cercaria to adult worm stage. Immunolocalization experiments and tegument preparations confirm that at least some of the SmCD59-like proteins are surface-localized; however, significant expression was also detected in internal tissues of adult worms. Finally, the involvement of two SmCD59 proteins in complement inhibition was evaluated by three different approaches: (i) a hemolytic assay using recombinant soluble forms expressed in Pichia pastoris and E. coli; (ii) complement-resistance of CHO cells expressing the respective membrane-anchored proteins; and (iii) the complement killing of schistosomula after gene suppression by RNAi. Our data indicated that these proteins are not involved in the regulation of complement activation. Conclusions: Our results suggest that this group of proteins belongs to the TFPD superfamily. Their expression is associated to intra-host stages, present in the tegument surface, and also in intra-parasite tissues. Three distinct approaches using SmCD59 proteins to inhibit complement strongly suggested that these proteins are not complement inhibitors and their function in schistosomes remains to be determined.