GEDI: a user-friendly toolbox for analysis of large-scale gene expression data

Abstract Background Several mathematical and statistical methods have been proposed in the last few years to analyze microarray data. Most of those methods involve complicated formulas, and software implementations that require advanced computer programming skills. Researchers from other areas may experience difficulties when they attempting to use those methods in their research. Here we present an user-friendly toolbox which allows large-scale gene expression analysis to be carried out by biomedical researchers with limited programming skills. Results Here, we introduce an user-friendly toolbox called GEDI (Gene Expression Data Interpreter), an extensible, open-source, and freely-available tool that we believe will be useful to a wide range of laboratories, and to researchers with no background in Mathematics and Computer Science, allowing them to analyze their own data by applying both classical and advanced approaches developed and recently published by Fujita et al. Conclusion GEDI is an integrated user-friendly viewer that combines the state of the art SVR, DVAR and SVAR algorithms, previously developed by us. It facilitates the application of SVR, DVAR and SVAR, further than the mathematical formulas present in the corresponding publications, and allows one to better understand the results by means of available visualizations. Both running the statistical methods and visualizing the results are carried out within the graphical user interface, rendering these algorithms accessible to the broad community of researchers in Molecular Biology.

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.

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.

RNA interference-mediated knockdown of CD49e (α5 integrin chain) in human thymic epithelial cells modulates the expression of multiple genes and decreases thymocyte adhesion

Abstract Background The thymus is a central lymphoid organ, in which bone marrow-derived T cell precursors undergo a complex process of maturation. Developing thymocytes interact with thymic microenvironment in a defined spatial order. A component of thymic microenvironment, the thymic epithelial cells, is crucial for the maturation of T-lymphocytes through cell-cell contact, cell matrix interactions and secretory of cytokines/chemokines. There is evidence that extracellular matrix molecules play a fundamental role in guiding differentiating thymocytes in both cortical and medullary regions of the thymic lobules. The interaction between the integrin α5β1 (CD49e/CD29; VLA-5) and fibronectin is relevant for thymocyte adhesion and migration within the thymic tissue. Our previous results have shown that adhesion of thymocytes to cultured TEC line is enhanced in the presence of fibronectin, and can be blocked with anti-VLA-5 antibody. Results Herein, we studied the role of CD49e expressed by the human thymic epithelium. For this purpose we knocked down the CD49e by means of RNA interference. This procedure resulted in the modulation of more than 100 genes, some of them coding for other proteins also involved in adhesion of thymocytes; others related to signaling pathways triggered after integrin activation, or even involved in the control of F-actin stress fiber formation. Functionally, we demonstrated that disruption of VLA-5 in human TEC by CD49e-siRNA-induced gene knockdown decreased the ability of TEC to promote thymocyte adhesion. Such a decrease comprised all CD4/CD8-defined thymocyte subsets. Conclusion Conceptually, our findings unravel the complexity of gene regulation, as regards key genes involved in the heterocellular cell adhesion between developing thymocytes and the major component of the thymic microenvironment, an interaction that is a mandatory event for proper intrathymic T cell differentiation.

Cell organisation, sulphur metabolism and ion transport-related genes are differentially expressed in Paracoccidioides brasiliensis mycelium and yeast cells

Abstract Background Mycelium-to-yeast transition in the human host is essential for pathogenicity by the fungus Paracoccidioides brasiliensis and both cell types are therefore critical to the establishment of paracoccidioidomycosis (PCM), a systemic mycosis endemic to Latin America. The infected population is of about 10 million individuals, 2% of whom will eventually develop the disease. Previously, transcriptome analysis of mycelium and yeast cells resulted in the assembly of 6,022 sequence groups. Gene expression analysis, using both in silico EST subtraction and cDNA microarray, revealed genes that were differential to yeast or mycelium, and we discussed those involved in sugar metabolism. To advance our understanding of molecular mechanisms of dimorphic transition, we performed an extended analysis of gene expression profiles using the methods mentioned above. Results In this work, continuous data mining revealed 66 new differentially expressed sequences that were MIPS(Munich Information Center for Protein Sequences)-categorised according to the cellular process in which they are presumably involved. Two well represented classes were chosen for further analysis: (i) control of cell organisation – cell wall, membrane and cytoskeleton, whose representatives were hex (encoding for a hexagonal peroxisome protein), bgl (encoding for a 1,3-β-glucosidase) in mycelium cells; and ags (an α-1,3-glucan synthase), cda (a chitin deacetylase) and vrp (a verprolin) in yeast cells; (ii) ion metabolism and transport – two genes putatively implicated in ion transport were confirmed to be highly expressed in mycelium cells – isc and ktp, respectively an iron-sulphur cluster-like protein and a cation transporter; and a putative P-type cation pump (pct) in yeast. Also, several enzymes from the cysteine de novo biosynthesis pathway were shown to be up regulated in the yeast form, including ATP sulphurylase, APS kinase and also PAPS reductase. Conclusion Taken together, these data show that several genes involved in cell organisation and ion metabolism/transport are expressed differentially along dimorphic transition. Hyper expression in yeast of the enzymes of sulphur metabolism reinforced that this metabolic pathway could be important for this process. Understanding these changes by functional analysis of such genes may lead to a better understanding of the infective process, thus providing new targets and strategies to control PCM.

A quantitative view of the transcriptome of Schistosoma mansoni adult-worms using SAGE

Abstract Background Five species of the genus Schistosoma, a parasitic trematode flatworm, are causative agents of Schistosomiasis, a disease that is endemic in a large number of developing countries, affecting millions of patients around the world. By using SAGE (Serial Analysis of Gene Expression) we describe here the first large-scale quantitative analysis of the Schistosoma mansoni transcriptome, one of the most epidemiologically relevant species of this genus. Results After extracting mRNA from pooled male and female adult-worms, a SAGE library was constructed and sequenced, generating 68,238 tags that covered more than 6,000 genes expressed in this developmental stage. An analysis of the ordered tag-list shows the genes of F10 eggshell protein, pol-polyprotein, HSP86, 14-3-3 and a transcript yet to be identified to be the five top most abundant genes in pooled adult worms. Whereas only 8% of the 100 most abundant tags found in adult worms of S. mansoni could not be assigned to transcripts of this parasite, 46.9% of the total ditags could not be mapped, demonstrating that the 3 sequence of most of the rarest transcripts are still to be identified. Mapping of our SAGE tags to S. mansoni genes suggested the occurrence of alternative-polyadenylation in at least 13 gene transcripts. Most of these events seem to shorten the 3 UTR of the mRNAs, which may have consequences over their stability and regulation. Conclusion SAGE revealed the frequency of expression of the majority of the S. mansoni genes. Transcriptome data suggests that alternative polyadenylation is likely to be used in the control of mRNA stability in this organism. When transcriptome was compared with the proteomic data available, we observed a correlation of about 50%, suggesting that both transcriptional and post-transcriptional regulation are important for determining protein abundance in S. mansoni. The generation of SAGE tags from other life-cycle stages should contribute to reveal the dynamics of gene expression in this important parasite.

Analysis of Schistosoma mansoni genes shared with Deuterostomia and with possible roles in host interactions

Abstract Background: Schistosoma mansoni is a blood helminth parasite that causes schistosomiasis, a disease that affects 200 million people in the world. Many orthologs of known mammalian genes have been discovered in this parasite and evidence is accumulating that some of these genes encode proteins linked to signaling pathways in the parasite that appear to be involved with growth or development, suggesting a complex co-evolutionary process. Results: In this work we found 427 genes conserved in the Deuterostomia group that have orthologs in S. mansoni and no members in any nematodes and insects so far sequenced. Among these genes we have identified Insulin Induced Gene (INSIG), Interferon Regulatory Factor (IRF) and vasohibin orthologs, known to be involved in mammals in mevalonate metabolism, immune response and angiogenesis control, respectively. We have chosen these three genes for a more detailed characterization, which included extension of their cloned messages to obtain full-length sequences. Interestingly, SmINSIG showed a 10-fold higher expression in adult females as opposed to males, in accordance with its possible role in regulating egg production. SmIRF has a DNA binding domain, a tryptophan-rich N-terminal region and several predicted phosphorylation sites, usually important for IRF activity. Fourteen different alternatively spliced forms of the S. mansoni vasohibin (SmVASL) gene were detected that encode seven different protein isoforms including one with a complete C-terminal end, and other isoforms with shorter C-terminal portions. Using S. mansoni homologs, we have employed a parsimonious rationale to compute the total gene losses/gains in nematodes, arthropods and deuterostomes under either the Coelomata or the Ecdysozoa evolutionary hypotheses; our results show a lower losses/gains number under the latter hypothesis. Conclusion: The genes discussed which are conserved between S. mansoni and deuterostomes, probably have an ancient origin and were lost in Ecdysozoa, being still present in Lophotrochozoa. Given their known functions in Deuterostomia, it is possible that some of them have been co-opted to perform functions related (directly or indirectly) to host adaptation or interaction with host signaling processes.

Sequence features responsible for intron retention in human

Abstract Background One of the least common types of alternative splicing is the complete retention of an intron in a mature transcript. Intron retention (IR) is believed to be the result of intron, rather than exon, definition associated with failure of the recognition of weak splice sites flanking short introns. Although studies on individual retained introns have been published, few systematic surveys of large amounts of data have been conducted on the mechanisms that lead to IR. Results TTo understand how sequence features are associated with or control IR, and to produce a generalized model that could reveal previously unknown signals that regulate this type of alternative splicing, we partitioned intron retention events observed in human cDNAs into two groups based on the relative abundance of both isoforms and compared relevant features. We found that a higher frequency of IR in human is associated with individual introns that have weaker splice sites, genes with shorter intron lengths, higher expression levels and lower density of both a set of exon splicing silencers (ESSs) and the intronic splicing enhancer GGG. Both groups of retained introns presented events conserved in mouse, in which the retained introns were also short and presented weaker splice sites. Conclusion Although our results confirmed that weaker splice sites are associated with IR, they showed that this feature alone cannot explain a non-negligible fraction of events. Our analysis suggests that cis-regulatory elements are likely to play a crucial role in regulating IR and also reveals previously unknown features that seem to influence its occurrence. These results highlight the importance of considering the interplay among these features in the regulation of the relative frequency of IR.

No evidence for an association between the -36A>C phospholamban gene polymorphism and a worse prognosis in heart failure

Abstract Background In Brazil, heart failure leads to approximately 25,000 deaths per year. Abnormal calcium handling is a hallmark of heart failure and changes in genes encoding for proteins involved in the re-uptake of calcium might harbor mutations leading to inherited cardiomyopathies. Phospholamban (PLN) plays a prime role in cardiac contractility and relaxation and mutations in the gene encoding PLN have been associated with dilated cardiomyopathy. In this study, our objective was to determine the presence of the -36A>C alteration in PLN gene in a Brazilian population of individuals with HF and to test whether this alteration is associated with heart failure or with a worse prognosis of patients with HF. Methods We genotyped a cohort of 881 patients with HF and 1259 individuals from a cohort of individuals from the general population for the alteration -36A>C in the PLN gene. Allele and genotype frequencies were compared between groups (patients and control). In addition, frequencies or mean values of different phenotypes associated with cardiovascular disease were compared between genotypic groups. Finally, patients were prospectively followed-up for death incidence and genotypes for the -36A>C were compared regarding mortality incidence in HF patients. Results No significant association was found between the study polymorphism and HF in our population. In addition, no association between PLN -36A>C polymorphism and demographic, clinical and functional characteristics and mortality incidence in this sample of HF patients was observed. Conclusion Our data do not support a role for the PLN -36A>C alteration in modulating the heart failure phenotype, including its clinical course, in humans.

Lactobacillus acidophilus CRL 1014 improved “gut health” in the SHIME® reactor

Abstract Background How to maintain “gut health” is a goal for scientists throughout the world. Therefore, microbiota management models for testing probiotics, prebiotics, and synbiotics have been developed. Methods The SHIME® model was used to study the effect of Lactobacillus acidophilus 1014 on the fermentation pattern of the colon microbiota. Initially, an inoculum prepared from human feces was introduced into the reactor vessels and stabilized over 2-wk using a culture medium. This stabilization period was followed by a 2-wk control period during which the microbiota was monitored. The microbiota was then subjected to a 4-wk treatment period by adding 5 mL of sterile peptone water with L. acidophilus CRL1014 at the concentration of 108 CFU/mL to vessel one (the stomach compartment). Plate counts, Denaturing Gradient Gel Electrophoresis (DGGE), short-chain fatty acid (SCFA) and ammonium analyses were carried out for monitoring of the microbial community from the colon compartments. Results A significant increase (p < 0.01) in the Lactobacillus spp. and Bifidobacterium spp. populations was observed during the treatment period. The DGGE obtained showed changes in the lactobacilli community from the colon compartments of the SHIME® reactor. The (SCFA) concentration increased (p < 0.01) during the treatment period, due mainly to significant increased levels of acetic, butyric, and propionic acids. However, ammonium concentrations decreased during the same period (p < 0.01). Conclusions This study showed the beneficial influence of L. acidophilus CRL 1014 on microbial metabolism and lactobacilli community composition for improving human health.

Chronic pain associated with the Chikungunya Fever: long lasting burden of an acute illness

Abstract Background Chikungunya virus (CHIKV) is responsible for major epidemics worldwide. Autochthonous cases were recently reported in several European countries. Acute infection is thought to be monophasic. However reports on chronic pain related to CHIKV infection have been made. In particular, the fact that many of these patients do not respond well to usual analgesics suggests that the nature of chronic pain may be not only nociceptive but also neuropathic. Neuropathic pain syndromes require specific treatment and the identification of neuropathic characteristics (NC) in a pain syndrome is a major step towards pain control. Methods We carried out a cross-sectional study at the end of the major two-wave outbreak lasting 17 months in Réunion Island. We assessed pain in 106 patients seeking general practitioners with confirmed infection with the CHIK virus, and evaluated its impact on quality of life (QoL). Results The mean intensity of pain on the visual-analogical scale (VAS) was 5.8 ± 2.1, and its mean duration was 89 ± 2 days. Fifty-six patients fulfilled the definition of chronic pain. Pain had NC in 18.9% according to the DN4 questionnaire. Conversely, about two thirds (65%) of patients with NC had chronic pain. The average pain intensity was similar between patients with or without NC (6.0 ± 1.7 vs 6.1 ± 2.0). However, the total score of the Short Form-McGill Pain Questionnaire (SF-MPQ)(15.5 ± 5.2 vs 11.6 ± 5.2; p < 0.01) and both the affective (18.8 ± 6.2 vs 13.4 ± 6.7; p < 0.01) and sensory subscores (34.3 ± 10.7 vs 25.0 ± 9.9; p < 0.01) were significantly higher in patients with NC. The mean pain interference in life activities calculated from the Brief Pain Inventory (BPI) was significantly higher in patients with chronic pain than in patients without it (6.8 ± 1.9 vs 5.9 ± 1.9, p < 0.05). This score was also significantly higher in patients with NC than in those without such a feature (7.2 ± 1.5 vs 6.1 ± 1.9, p < 0.05). Conclusions There exists a specific chronic pain condition associated to CHIKV. Pain with NC seems to be associated with more aggressive clinical picture, more intense impact in QoL and more challenging pharmacological treatment.

Brazilian multicenter study on prevalence of preterm birth and associated factors

Abstract Background The occurrence of preterm birth remains a complex public health condition. It is considered the main cause of neonatal morbidity and mortality, resulting in a high likelihood of sequelae in surviving children. With variable incidence in several countries, it has grown markedly in the last decades. In Brazil, however, there are still difficulties to estimate its real occurrence. Therefore, it is essential to establish the prevalence and causes of this condition in order to propose prevention actions. This study intend to collect information from hospitals nationwide on the prevalence of preterm births, their associated socioeconomic and environmental factors, diagnostic and treatment methods resulting from causes such as spontaneous preterm labor, prelabor rupture of membranes, and therapeutic preterm birth, as well as neonatal results. Methods/Design This proposal is a multicenter cross-sectional study plus a nested case-control study, to be implemented in 27 reference obstetric centers in several regions of Brazil (North: 1; Northeast: 10; Central-west: 1; Southeast: 13; South: 2). For the cross sectional component, the participating centers should perform, during a period of six months, a prospective surveillance of all patients hospitalized to give birth, in order to identify preterm birth cases and their main causes. In the first three months of the study, an analysis of the factors associated with preterm birth will also be carried out, comparing women who have preterm birth with those who deliver at term. For the prevalence study, 37,000 births will be evaluated (at term and preterm), corresponding to approximately half the deliveries of all participating centers in 12 months. For the case-control study component, the estimated sample size is 1,055 women in each group (cases and controls). The total number of preterm births estimated to be followed in both components of the study is around 3,600. Data will be collected through a questionnaire all patients will answer after delivery. The data will then be encoded in an electronic form and sent online by internet to a central database. The data analysis will be carried out by subgroups according to gestational age at preterm birth, its probable causes, therapeutic management, and neonatal outcomes. Then, the respective rates, ratios and relative risks will be estimated for the possible predictors. Discussion These findings will provide information on preterm births in Brazil and their main social and biological risk factors, supporting health policies and the implementation of clinical trials on preterm birth prevention and treatment strategies, a condition with many physical and emotional consequences to children and their families.

Predicting smear negative pulmonary tuberculosis with classification trees and logistic regression: a cross-sectional study

Abstract Background Smear negative pulmonary tuberculosis (SNPT) accounts for 30% of pulmonary tuberculosis cases reported yearly in Brazil. This study aimed to develop a prediction model for SNPT for outpatients in areas with scarce resources. Methods The study enrolled 551 patients with clinical-radiological suspicion of SNPT, in Rio de Janeiro, Brazil. The original data was divided into two equivalent samples for generation and validation of the prediction models. Symptoms, physical signs and chest X-rays were used for constructing logistic regression and classification and regression tree models. From the logistic regression, we generated a clinical and radiological prediction score. The area under the receiver operator characteristic curve, sensitivity, and specificity were used to evaluate the model's performance in both generation and validation samples. Results It was possible to generate predictive models for SNPT with sensitivity ranging from 64% to 71% and specificity ranging from 58% to 76%. Conclusion The results suggest that those models might be useful as screening tools for estimating the risk of SNPT, optimizing the utilization of more expensive tests, and avoiding costs of unnecessary anti-tuberculosis treatment. Those models might be cost-effective tools in a health care network with hierarchical distribution of scarce resources.

The presence of a booster phenomenon among contacts of active pulmonary tuberculosis cases: a retrospective cohort

Abstract Background Assuming a higher risk of latent tuberculosis (TB) infection in the population of Rio de Janeiro, Brazil, in October of 1998 the TB Control Program of Clementino Fraga Filho Hospital (CFFH) routinely started to recommend a two-step tuberculin skin test (TST) in contacts of pulmonary TB cases in order to distinguish a boosting reaction due to a recall of delayed hypersensitivity previously established by infection with Mycobacterium tuberculosis (M.tb) or BCG vaccination from a tuberculin conversion. The aim of this study was to assess the prevalence of boosted tuberculin skin tests among contacts of individuals with active pulmonary tuberculosis (TB). Methods Retrospective cohort of TB contacts ≥ 12 years old who were evaluated between October 1st, 1998 and October 31st 2001. Contacts with an initial TST ≤ 4 mm were considered negative and had a second TST applied after 7–14 days. Boosting reaction was defined as a second TST ≥ 10 mm with an increase in induration ≥ 6 mm related to the first TST. All contacts with either a positive initial or repeat TST had a chest x-ray to rule out active TB disease, and initially positive contacts were offered isoniazid preventive therapy. Contacts that boosted did not receive treatment for latent TB infection and were followed for 24 months to monitor the development of TB. Statistical analysis of dichotomous variables was performed using Chi-square test. Differences were considered significant at a p < 0.05. Results Fifty four percent (572/1060) of contacts had an initial negative TST and 79% of them (455/572) had a second TST. Boosting was identified in 6% (28/455). The mean age of contacts with a boosting reaction was 42.3 ± 21.1 and with no boosting was 28.7 ± 21.7 (p = 0.01). Fifty percent (14/28) of individuals whose test boosted met criteria for TST conversion on the second TST (increase in induration ≥ 10 mm). None of the 28 contacts whose reaction boosted developed TB disease within two years following the TST. Conclusion The low number of contacts with boosting and the difficulty in distinguishing boosting from TST conversion in the second TST suggests that the strategy of two-step TST testing among contacts of active TB cases may not be useful. However, this conclusion must be taken with caution because of the small number of subjects followed.

Acute hepatitis C virus infection assessment among chronic hemodialysis patients in the Southwest Parana State, Brazil

Abstract Background Chronic hemodialysis patients are at higher risk for acquiring hepatitis C virus (HCV). The prevalence varies among different countries and hemodialysis centers. Although guidelines for a comprehensive infection control program exist, the nosocomial transmission still accounts for the new cases of infection. The aim of this study was analyze the follow up of newly acquired acute hepatitis C cases, during the period from January 2002 to May 2005, in the Hemodialysis Center, located in the Southwest region of Parana State, Brazil and to analyze the effectiveness of the measures to restrain the appearance of new cases of acute hepatitis C. Methods Patients were analyzed monthly with anti-HCV tests and ALT measurements. Patients with ALT elevations were monitored for possible acute hepatitis C. Results During this period, 32 new cases were identified with acute hepatitis C virus infection. Blood screening showed variable ALT levels preceding the anti-HCV seroconversion. HCV RNA viremia by PCR analysis was intermittently and even negative in some cases. Ten out of 32 patients received 1 mcg/kg dose of pegylated interferon alfa-2b treatment for 24 weeks. All dialysis personnel were re-trained to strictly follow the regulations and recommendations regarding infection control, proper methods to clean and disinfect equipment were reviewed and HCV-positive patients were isolated. Conclusion Laboratory tests results showed variable ALT preceding anti-HCV seroconversion and intermittent viremia. The applied recommendations contributed importantly to restrain the appearance of new cases of acute hepatitis C in this center and the last case was diagnosed in May 2004.