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.

Highlights of the XXI annual meeting of the Brazilian Society of Protozoology, the XXXII annual meeting on Basic Research in Chagas' disease & an international symposium on vesicle trafficking in parasitic Protozoa – 7 to 9 November 2005, Caxambu, Minas Gerais, Brazil

This report focuses on the 2005 Annual meeting held in Caxambu, Minas Gerais, Brazil that was convened and organized by the Brazilian Society of Protozoology http://www.sbpz.org.br/. This is an annual event and details of these meetings can be found on the Society's website. Within the space available it has been impossible to cover all the important and fascinating contributions and what is presented are our personal views of the meetings scientific highlights and new developments. The contents undoubtedly reflect each author's scientific interests and expertise. Fuller details of the round tables, seminars and posters can be consulted on line at http://www.sbpz.org.br/livroderesumos2005.php.

Independent of ErbB1 gene copy number, EGF stimulates migration but is not associated with cell proliferation in non-small cell lung cancer

Abstract Background Lung cancer often exhibits molecular changes, such as the overexpression of the ErbB1 gene. ErbB1 encodes epidermal growth factor receptor (EGFR), a tyrosine kinase receptor, involved mainly in cell proliferation and survival. EGFR overexpression has been associated with more aggressive disease, poor prognosis, low survival rate and low response to therapy. ErbB1 amplification and mutation are associated with tumor development and are implicated in ineffective treatment. The aim of the present study was to investigate whether the ErbB1 copy number affects EGFR expression, cell proliferation or cell migration by comparing two different cell lines. Methods The copies of ErbB1 gene was evaluated by FISH. Immunofluorescence and Western blotting were performed to determine location and expression of proteins mentioned in the present study. Proliferation was studied by flow cytometry and cell migration by wound healing assay and time lapse. Results We investigated the activation and function of EGFR in the A549 and HK2 lung cancer cell lines, which contain 3 and 6 copies of ErbB1, respectively. The expression of EGFR was lower in the HK2 cell line. EGFR was activated after stimulation with EGF in both cell lines, but this activation did not promote differences in cellular proliferation when compared to control cells. Inhibiting EGFR with AG1478 did not modify cellular proliferation, confirming previous data. However, we observed morphological alterations, changes in microfilament organization and increased cell migration upon EGF stimulation. However, these effects did not seem to be consequence of an epithelial-mesenchymal transition. Conclusion EGFR expression did not appear to be associated to the ErbB1 gene copy number, and neither of these aspects appeared to affect cell proliferation. However, EGFR activation by EGF resulted in cell migration stimulation in both cell lines.

Copy number variation in Williams-Beuren syndrome: suitable diagnostic strategy for developing countries

Abstract Background Williams-Beuren syndrome (WBS; OMIM 194050) is caused by a hemizygous contiguous gene microdeletion at 7q11.23. Supravalvular aortic stenosis (SVAS), mental retardation, and overfriendliness comprise typical symptoms of WBS. Although fluorescence in situ hybridization (FISH) is considered the gold standard technique, the microsatellite DNA markers and multiplex ligation-dependent probe amplification (MLPA) could be used for to confirm the diagnosis of WBS. Results We have evaluated a total cohort of 88 patients with a suspicion clinical diagnosis of WBS using a collection of five markers (D7S1870, D7S489, D7S613, D7S2476, and D7S489_A) and a commercial MLPA kit (P029). The microdeletion was present in 64 (72.7%) patients and absent in 24 (27.3%) patients. The parental origin of deletion was maternal in 36 of 64 patients (56.3%) paternal in 28 of 64 patients (43.7%). The deletion size was 1.55 Mb in 57 of 64 patients (89.1%) and 1.84 Mb in 7 of 64 patients (10.9%). The results were concordant using both techniques, except for four patients whose microsatellite markers were uninformative. There were no clinical differences in relation to either the size or parental origin of the deletion. Conclusion MLPA was considered a faster and more economical method in a single assay, whereas the microsatellite markers could determine both the size and parental origin of the deletion in WBS. The microsatellite marker and MLPA techniques are effective in deletion detection in WBS, and both methods provide a useful diagnostic strategy mainly for developing countries.

The future of medical diagnostics: review paper

While histopathology of excised tissue remains the gold standard for diagnosis, several new, non-invasive diagnostic techniques are being developed. They rely on physical and biochemical changes that precede and mirror malignant change within tissue. The basic principle involves simple optical techniques of tissue interrogation. Their accuracy, expressed as sensitivity and specificity, are reported in a number of studies suggests that they have a potential for cost effective, real-time, in situ diagnosis.

Semantic integration of gene expression analysis tools and data sources using software connectors

Abstract Background The study and analysis of gene expression measurements is the primary focus of functional genomics. Once expression data is available, biologists are faced with the task of extracting (new) knowledge associated to the underlying biological phenomenon. Most often, in order to perform this task, biologists execute a number of analysis activities on the available gene expression dataset rather than a single analysis activity. The integration of heteregeneous tools and data sources to create an integrated analysis environment represents a challenging and error-prone task. Semantic integration enables the assignment of unambiguous meanings to data shared among different applications in an integrated environment, allowing the exchange of data in a semantically consistent and meaningful way. This work aims at developing an ontology-based methodology for the semantic integration of gene expression analysis tools and data sources. The proposed methodology relies on software connectors to support not only the access to heterogeneous data sources but also the definition of transformation rules on exchanged data. Results We have studied the different challenges involved in the integration of computer systems and the role software connectors play in this task. We have also studied a number of gene expression technologies, analysis tools and related ontologies in order to devise basic integration scenarios and propose a reference ontology for the gene expression domain. Then, we have defined a number of activities and associated guidelines to prescribe how the development of connectors should be carried out. Finally, we have applied the proposed methodology in the construction of three different integration scenarios involving the use of different tools for the analysis of different types of gene expression data. Conclusions The proposed methodology facilitates the development of connectors capable of semantically integrating different gene expression analysis tools and data sources. The methodology can be used in the development of connectors supporting both simple and nontrivial processing requirements, thus assuring accurate data exchange and information interpretation from exchanged data.

MYC, FBXW7 and TP53 copy number variation and expression in Gastric Cancer

Abstract Background MYC deregulation is a common event in gastric carcinogenesis, usually as a consequence of gene amplification, chromosomal translocations, or posttranslational mechanisms. FBXW7 is a p53-controlled tumor-suppressor that plays a role in the regulation of cell cycle exit and reentry via MYC degradation. Methods We evaluated MYC, FBXW7, and TP53 copy number, mRNA levels, and protein expression in gastric cancer and paired non-neoplastic specimens from 33 patients and also in gastric adenocarcinoma cell lines. We also determined the invasion potential of the gastric cancer cell lines. Results MYC amplification was observed in 51.5% of gastric tumor samples. Deletion of one copy of FBXW7 and TP53 was observed in 45.5% and 21.2% of gastric tumors, respectively. MYC mRNA expression was significantly higher in tumors than in non-neoplastic samples. FBXW7 and TP53 mRNA expression was markedly lower in tumors than in paired non-neoplastic specimens. Moreover, deregulated MYC and FBXW7 mRNA expression was associated with the presence of lymph node metastasis and tumor stage III-IV. Additionally, MYC immunostaining was more frequently observed in intestinal-type than diffuse-type gastric cancers and was associated with MYC mRNA expression. In vitro studies showed that increased MYC and reduced FBXW7 expression is associated with a more invasive phenotype in gastric cancer cell lines. This result encouraged us to investigate the activity of the gelatinases MMP-2 and MMP-9 in both cell lines. Both gelatinases are synthesized predominantly by stromal cells rather than cancer cells, and it has been proposed that both contribute to cancer progression. We observed a significant increase in MMP-9 activity in ACP02 compared with ACP03 cells. These results confirmed that ACP02 cells have greater invasion capability than ACP03 cells. Conclusion In conclusion, FBXW7 and MYC mRNA may play a role in aggressive biologic behavior of gastric cancer cells and may be a useful indicator of poor prognosis. Furthermore, MYC is a candidate target for new therapies against gastric cancer.

A mathematical model for optimizing the indications of liver transplantation in patients with hepatocellular carcinoma

Abstract Background The criteria for organ sharing has developed a system that prioritizes liver transplantation (LT) for patients with hepatocellular carcinoma (HCC) who have the highest risk of wait-list mortality. In some countries this model allows patients only within the Milan Criteria (MC, defined by the presence of a single nodule up to 5 cm, up to three nodules none larger than 3 cm, with no evidence of extrahepatic spread or macrovascular invasion) to be evaluated for liver transplantation. This police implies that some patients with HCC slightly more advanced than those allowed by the current strict selection criteria will be excluded, even though LT for these patients might be associated with acceptable long-term outcomes. Methods We propose a mathematical approach to study the consequences of relaxing the MC for patients with HCC that do not comply with the current rules for inclusion in the transplantation candidate list. We consider overall 5-years survival rates compatible with the ones reported in the literature. We calculate the best strategy that would minimize the total mortality of the affected population, that is, the total number of people in both groups of HCC patients that die after 5 years of the implementation of the strategy, either by post-transplantation death or by death due to the basic HCC. We illustrate the above analysis with a simulation of a theoretical population of 1,500 HCC patients with tumor size exponentially. The parameter λ obtained from the literature was equal to 0.3. As the total number of patients in these real samples was 327 patients, this implied in an average size of 3.3 cm and a 95% confidence interval of [2.9; 3.7]. The total number of available livers to be grafted was assumed to be 500. Results With 1500 patients in the waiting list and 500 grafts available we simulated the total number of deaths in both transplanted and non-transplanted HCC patients after 5 years as a function of the tumor size of transplanted patients. The total number of deaths drops down monotonically with tumor size, reaching a minimum at size equals to 7 cm, increasing from thereafter. With tumor size equals to 10 cm the total mortality is equal to the 5 cm threshold of the Milan criteria. Conclusion We concluded that it is possible to include patients with tumor size up to 10 cm without increasing the total mortality of this population.

Transposable elements: structure and dynamic of the autonomous retrotransposon R2 in Arthropoda with non-canonical reproduction

Eukaryotic ribosomal DNA constitutes a multi gene family organized in a cluster called nucleolar organizer region (NOR); this region is composed usually by hundreds to thousands of tandemly repeated units. Ribosomal genes, being repeated sequences, evolve following the typical pattern of concerted evolution. The autonomous retroelement R2 inserts in the ribosomal gene 28S, leading to defective 28S rDNA genes. R2 element, being a retrotransposon, performs its activity in the genome multiplying its copy number through a “copy and paste” mechanism called target primed reverse transcription. It consists in the retrotranscription of the element’s mRNA into DNA, then the DNA is integrated in the target site. Since the retrotranscription can be interrupted, but the integration will be carried out anyway, truncated copies of the element will also be present in the genome. The study of these truncated variants is a tool to examine the activity of the element. R2 phylogeny appears, in general, not consistent with that of its hosts, except some cases (e.g. Drosophila spp. and Reticulitermes spp.); moreover R2 is absent in some species (Fugu rubripes, human, mouse, etc.), while other species have more R2 lineages in their genome (the turtle Mauremys reevesii, the Japanese beetle Popilia japonica, etc). R2 elements here presented are isolated in 4 species of notostracan branchiopods and in two species of stick insects, whose reproductive strategies range from strict gonochorism to unisexuality. From sequencing data emerges that in Triops cancriformis (Spanish gonochoric population), in Lepidurus arcticus (two putatively unisexual populations from Iceland) and in Bacillus rossius (gonochoric population from Capalbio) the R2 elements are complete and encode functional proteins, reflecting the general features of this family of transposable elements. On the other hand, R2 from Italian and Austrian populations of T. cancriformis (respectively unisexual and hermaphroditic), Lepidurus lubbocki (two elements within the same Italian population, gonochoric but with unfunctional males) and Bacillus grandii grandii (gonochoric population from Ponte Manghisi) have sequences that encode incomplete or non-functional proteins in which it is possible to recognize only part of the characteristic domains. In Lepidurus couesii (Italian gonochoric populations) different elements were found as in L. lubbocki, and the sequencing is still in progress. Two hypothesis are given to explain the inconsistency of R2/host phylogeny: vertical inheritance of the element followed by extinction/diversification or horizontal transmission. My data support previous study that state the vertical transmission as the most likely explanation; nevertheless horizontal transfer events can’t be excluded. I also studied the element’s activity in Spanish populations of T. cancriformis, in L. lubbocki, in L. arcticus and in gonochoric and parthenogenetic populations of B. rossius. In gonochoric populations of T. cancriformis and B. rossius I found that each individual has its own private set of truncated variants. The situation is the opposite for the remaining hermaphroditic/parthenogenetic species and populations, all individuals sharing – in the so far analyzed samples - the majority of variants. This situation is very interesting, because it isn’t concordant with the Muller’s ratchet theory that hypothesizes the parthenogenetic populations being either devoided of transposable elements or TEs overloaded. My data suggest a possible epigenetic mechanism that can block the retrotransposon activity, and in this way deleterious mutations don’t accumulate.

Partial reconstruction of the trajectories of a discretely observed branching diffusion with immigration and an application to inference

In this treatise we consider finite systems of branching particles where the particles move independently of each other according to d-dimensional diffusions. Particles are killed at a position dependent rate, leaving at their death position a random number of descendants according to a position dependent reproduction law. In addition particles immigrate at constant rate (one immigrant per immigration time). A process with above properties is called a branching diffusion withimmigration (BDI). In the first part we present the model in detail and discuss the properties of the BDI under our basic assumptions. In the second part we consider the problem of reconstruction of the trajectory of a BDI from discrete observations. We observe positions of the particles at discrete times; in particular we assume that we have no information about the pedigree of the particles. A natural question arises if we want to apply statistical procedures on the discrete observations: How can we find couples of particle positions which belong to the same particle? We give an easy to implement 'reconstruction scheme' which allows us to redraw or 'reconstruct' parts of the trajectory of the BDI with high accuracy. Moreover asymptotically the whole path can be reconstructed. Further we present simulations which show that our partial reconstruction rule is tractable in practice. In the third part we study how the partial reconstruction rule fits into statistical applications. As an extensive example we present a nonparametric estimator for the diffusion coefficient of a BDI where the particles move according to one-dimensional diffusions. This estimator is based on the Nadaraya-Watson estimator for the diffusion coefficient of one-dimensional diffusions and it uses the partial reconstruction rule developed in the second part above. We are able to prove a rate of convergence of this estimator and finally we present simulations which show that the estimator works well even if we leave our set of assumptions.

Ergodicity and regularity of invariant measure for branching Markov processes with immigration

In this thesis we consider systems of finitely many particles moving on paths given by a strong Markov process and undergoing branching and reproduction at random times. The branching rate of a particle, its number of offspring and their spatial distribution are allowed to depend on the particle's position and possibly on the configuration of coexisting particles. In addition there is immigration of new particles, with the rate of immigration and the distribution of immigrants possibly depending on the configuration of pre-existing particles as well. In the first two chapters of this work, we concentrate on the case that the joint motion of particles is governed by a diffusion with interacting components. The resulting process of particle configurations was studied by E. Löcherbach (2002, 2004) and is known as a branching diffusion with immigration (BDI). Chapter 1 contains a detailed introduction of the basic model assumptions, in particular an assumption of ergodicity which guarantees that the BDI process is positive Harris recurrent with finite invariant measure on the configuration space. This object and a closely related quantity, namely the invariant occupation measure on the single-particle space, are investigated in Chapter 2 where we study the problem of the existence of Lebesgue-densities with nice regularity properties. For example, it turns out that the existence of a continuous density for the invariant measure depends on the mechanism by which newborn particles are distributed in space, namely whether branching particles reproduce at their death position or their offspring are distributed according to an absolutely continuous transition kernel. In Chapter 3, we assume that the quantities defining the model depend only on the spatial position but not on the configuration of coexisting particles. In this framework (which was considered by Höpfner and Löcherbach (2005) in the special case that branching particles reproduce at their death position), the particle motions are independent, and we can allow for more general Markov processes instead of diffusions. The resulting configuration process is a branching Markov process in the sense introduced by Ikeda, Nagasawa and Watanabe (1968), complemented by an immigration mechanism. Generalizing results obtained by Höpfner and Löcherbach (2005), we give sufficient conditions for ergodicity in the sense of positive recurrence of the configuration process and finiteness of the invariant occupation measure in the case of general particle motions and offspring distributions.

Reporting Experiments in Homeopathic Basic Research—Description of the Checklist Development

The objective of this study was to develop a criteria catalogue serving as a guideline for authors to improve quality of reporting experiments in basic research in homeopathy. A Delphi Process was initiated including three rounds of adjusting and phrasing plus two consensus conferences. European researchers who published experimental work within the last 5 years were involved. A checklist for authors provide a catalogue with 23 criteria. The “Introduction” should focus on underlying hypotheses, the homeopathic principle investigated and state if experiments are exploratory or confirmatory. “Materials and methods” should comprise information on object of investigation, experimental setup, parameters, intervention and statistical methods. A more detailed description on the homeopathic substances, for example, manufacture, dilution method, starting point of dilution is required. A further result of the Delphi process is to raise scientists' awareness of reporting blinding, allocation, replication, quality control and system performance controls. The part “Results” should provide the exact number of treated units per setting which were included in each analysis and state missing samples and drop outs. Results presented in tables and figures are as important as appropriate measures of effect size, uncertainty and probability. “Discussion” in a report should depict more than a general interpretation of results in the context of current evidence but also limitations and an appraisal of aptitude for the chosen experimental model. Authors of homeopathic basic research publications are encouraged to apply our checklist when preparing their manuscripts. Feedback is encouraged on applicability, strength and limitations of the list to enable future revisions.

Management, feeding, production, reproduction and udder health on organic and conventional Swiss dairy farms

Organic dairy farms (OP; n=60) and conventional dairy farms (integrated production, IP; n=60), matched in size, location, and agricultural zone (altitude), were studied for possible differences in management, feeding, production, reproduction and udder health. OP and IP farms were similar in size (17.7 and 16.9 ha), milk quota (65900 and 70,000 kg/year), cow number (14 and 15), cow age (5.3 and 5.2 years), housing of cows of the Simmental x Red Holstein or Holstein breeds (87 and 75%; 45 and 60%), but differed significantly with respect to loose housing systems (18 and 7%), outside paddocks (98 and 75%), energy-corrected 305-d milk yield (5,695 and 6,059 kg), milk protein content (31.8 and 32.7 g/kg), use of bucket milking systems (73 and 33%), observance of regular (12-h) milking intervals (47 and 68%), routine application of the California-Mastitis-Test (10 and 28%), teat dipping after milking (25 and 43%) and blanket dry cow treatments (0 and 45%). Milk somatic cell counts on OP and IP farms (119 000 and 117,000/mL) and reproduction data were similar and there were no significant differences between OP and IP farms as concerns available feeds, planning and management of feeding. Alternative veterinary treatments were used more often on OP than IP farms (55 and 17%). Main causes for cow replacements on OP and IP farms were fertility disorders (both 45%), age (40 and 42%), sale (30 and 37%) and udder health (35 and 13%).Between OP and IP Swiss dairy farms thus relatively few larger differences were found.