BATS: a Bayesian user-friendly software for Analyzing Time Series microarray experiments

BATS is a user-friendly software for Bayesian Analysis of Time Series microarray experiments based on the novel, truly functional and fully Bayesian approach proposed in Angelini et at. (2006). The software is specifically designed for time series data. It allows an user to automatically identify and rank differentially expressed genes and to estimate their expression profiles. BATS successfully manages various technical difficulties which arise in microarray time-course experiments, such as a small number of observations, non-uniform sampling intervals, and presence of missing or multiple data. BATS can carry out analysis with both simulated and real experimental data. It also handles data from different platforms. 1 Availability: BATS is written in Matlab and executable in Windows (Macintosh and Linux version are currently under development). It is freely available upon request from the authors.

A qualitative study of the views of patients with long-term conditions on family doctors in Hong Kong

<b>Background</b> Primary care based management of long-term conditions (LTCs) is high on the international healthcare agenda, including the Asia-Pacific region. Hong Kong has a 'mixed economy' healthcare system with both public and private sectors with a range of types of primary care doctors. Recent Hong Kong Government policy aims to enhance the management of LTCs in primary care possibly based on a 'family doctor' model. Patients' views on this are not well documented and the aim of the present study was to explore the views of patients with LTCs on family doctors in Hong Kong.<p></p> <b>Methods</b> The views of patients (with a variety of LTCs) on family doctors in Hong Kong were explored. Two groups of participants were interviewed; a) those who considered themselves as having a family doctor, b) those who considered themselves as not having a family doctor (either with a regular primary care doctor but not a family doctor or with no regular primary care doctor). In-depth individual semi-structured interviews were carried out with 28 participants (10 with a family doctor, 10 with a regular doctor, and 8 with no regular doctor) and analysed using the constant comparative method.<p></p> <b>Results</b> Participants who did not have a family doctor were familiar with the concept but regarded it as a 'luxury item' for the rich within the private healthcare system. Those with a regular family doctor (all private) regarded having one as important to their and their family's health. Participants in both groups felt that as well as the more usual family medicine specialist or general practitioner, traditional Chinese medicine practitioners also had the potential to be family doctors. However most participants attended the public healthcare system for management of their LTCs whether they had a family doctor or not. Cost, perceived need, quality, trust, and choice were all barriers to the use of family doctors for the management of their LTCs.<p></p> <b>Conclusions</b> Important barriers to the adoption of a 'family doctor' model of management of LTCs exist in Hong Kong. Effective policy implementation seems unlikely unless these complex barriers are addressed

Comparison of coplanar and noncoplanar intensity-modulated radiation therapy and helical tomotherapy for hepatocellular carcinoma

<p>Abstract</p> <p>Background</p> <p>To compare the differences in dose-volume data among coplanar intensity modulated radiotherapy (IMRT), noncoplanar IMRT, and helical tomotherapy (HT) among patients with hepatocellular carcinoma (HCC) and portal vein thrombosis (PVT).</p> <p>Methods</p> <p>Nine patients with unresectable HCC and PVT underwent step and shoot coplanar IMRT with intent to deliver 46 - 54 Gy to the tumor and portal vein. The volume of liver received 30Gy was set to keep less than 30% of whole normal liver (V30 < 30%). The mean dose to at least one side of kidney was kept below 23 Gy, and 50 Gy as for stomach. The maximum dose was kept below 47 Gy for spinal cord. Several parameters including mean hepatic dose, percent volume of normal liver with radiation dose at X Gy (Vx), uniformity index, conformal index, and doses to organs at risk were evaluated from the dose-volume histogram.</p> <p>Results</p> <p>HT provided better uniformity for the planning-target volume dose coverage than both IMRT techniques. The noncoplanar IMRT technique reduces the V10 to normal liver with a statistically significant level as compared to HT. The constraints for the liver in the V30 for coplanar IMRT vs. noncoplanar IMRT vs. HT could be reconsidered as 21% vs. 17% vs. 17%, respectively. When delivering 50 Gy and 60-66 Gy to the tumor bed, the constraints of mean dose to the normal liver could be less than 20 Gy and 25 Gy, respectively.</p> <p>Conclusion</p> <p>Noncoplanar IMRT and HT are potential techniques of radiation therapy for HCC patients with PVT. Constraints for the liver in IMRT and HT could be stricter than for 3DCRT.</p

Detection of monoclonal immunoglobulin heavy chain gene rearrangement (FR3) in Thai malignant lymphoma by High Resolution Melting curve analysis

<p>Abstract</p> <p>Malignant lymphoma, especially non-Hodgkin lymphoma, is one of the most common hematologic malignancies in Thailand. The diagnosis of malignant lymphoma is often problematic, especially in early stages of the disease. Detection of antigen receptor gene rearrangement including T cell receptor (TCR) and immunoglobulin heavy chain (IgH) by polymerase chain reaction followed by heteroduplex has currently become standard whereas fluorescent fragment analysis (GeneScan) has been used for confirmation test. In this study, three techniques had been compared: thermocycler polymerase chain reaction (PCR) followed by heteroduplex and polyacrylamide gel electrophoresis, GeneScan analysis, and real time PCR with High Resolution Melting curve analysis (HRM). The comparison was carried out with DNA extracted from paraffin embedded tissues diagnosed as B- cell non-Hodgkin lymphoma. Specific PCR primers sequences for IgH gene variable region 3, including fluorescence labeled IgH primers were used and results were compared with HRM. In conclusion, the detection IgH gene rearrangement by HRM in the LightCycler System showed potential for distinguishing monoclonality from polyclonality in B-cell non-Hodgkin lymphoma.</p> <p>Introduction</p> <p>Malignant lymphoma, especially non-Hodgkin lymphoma, is one of the most common hematologic malignancies in Thailand. The incidence rate as reported by Ministry of Public Health is 3.1 per 100,000 population in female whereas the rate in male is 4.5 per 100,000 population <abbrgrp><abbr bid="B1">1</abbr></abbrgrp>. At Siriraj Hospital, the new cases diagnosed as malignant lymphoma were 214.6 cases/year <abbrgrp><abbr bid="B2">2</abbr></abbrgrp>. The diagnosis of malignant lymphoma is often problematic, especially in early stages of the disease. Therefore, detection of antigen receptor gene rearrangement including T cell receptor (TCR) and immunoglobulin heavy chain (IgH) by polymerase chain reaction (PCR) assay has recently become a standard laboratory test for discrimination of reactive from malignant clonal lymphoproliferation <abbrgrp><abbr bid="B3">3</abbr><abbr bid="B4">4</abbr></abbrgrp>. Analyzing DNA extracted from formalin-fixed, paraffin-embedded tissues by multiplex PCR techniques is more rapid, accurate and highly sensitive. Measuring the size of the amplicon from PCR analysis could be used to diagnose malignant lymphoma with monoclonal pattern showing specific and distinct bands detected on acrylamide gel electrophoresis. However, this technique has some limitations and some patients might require a further confirmation test such as GeneScan or fragment analysis <abbrgrp><abbr bid="B5">5</abbr><abbr bid="B6">6</abbr></abbrgrp>.</p> <p>GeneScan technique or fragment analysis reflects size and peak of DNA by using capillary gel electrophoresis. This technique is highly sensitive and can detect 0.5-1% of clonal lymphoid cells. It measures the amplicons by using various fluorescently labeled primers at forward or reverse sides and a specific size standard. Using a Genetic Analyzer machine and GeneMapper software (Applied Bioscience, USA), the monoclonal pattern revealed one single, sharp and high peak at the specific size corresponding to acrylamide gel pattern, whereas the polyclonal pattern showed multiple and small peak condensed at the same size standard. This technique is the most sensitive and accurate technique; however, it usually requires high technical experience and is also of high cost <abbrgrp><abbr bid="B7">7</abbr></abbrgrp>. Therefore, rapid and more cost effective technique are being sought.</p> <p>LightCycler PCR performs the diagnostic detection of amplicon via melting curve analysis within 2 hours with the use of a specific dye <abbrgrp><abbr bid="B8">8</abbr><abbr bid="B9">9</abbr></abbrgrp>. This dye consists of two types: one known as SYBR-Green I which is non specific and the other named as High Resolution Melting analysis (HRM) which is highly sensitive, more accurate and stable. Several reports demonstrated that this new instrument combined with DNA intercalating dyes can be used to discriminate sequence changes in PCR amplicon without manual handling of PCR product <abbrgrp><abbr bid="B10">10</abbr><abbr bid="B11">11</abbr></abbrgrp>. Therefore, current investigations using melting curve analysis are being developed <abbrgrp><abbr bid="B12">12</abbr><abbr bid="B13">13</abbr></abbrgrp>.</p> <p>In this study, three different techniques were compared to evaluate the suitability of LightCycler PCR with HRM as the clonal diagnostic tool for IgH gene rearrangement in B-cell non-Hogdkin lymphoma, i.e. thermocycler PCR followed by heteroduplex analysis and PAGE, GeneScan analysis and LightCycler PCR with HRM.</p

Important prognostic factors for the long-term survival of lung cancer subjects in Taiwan

<p>Abstract</p> <p>Background</p> <p>This study used a large-scale cancer database in determination of prognostic factors for the survival of lung cancer subjects in Taiwan.</p> <p>Methods</p> <p>Total of 24,910 subjects diagnosed with lung cancer was analysed. Survival estimates by Kaplan-Meier methods. Cox proportional-hazards model estimated the death risk (hazard ratio (HR)) for various prognostic factors.</p> <p>Results</p> <p>The prognostic indicators associated with a higher risk of lung cancer deaths are male gender (males versus females; HR = 1.07, 95% confidence intervals (CI): 1.03–1.11), males diagnosed in later periods (shown in 1991–1994 versus 1987–1990; HR = 1.13), older age at diagnosis, large cell carcinoma (LCC)/small cell carcinoma (SCC), and supportive care therapy over chemotherapy. The overall 5-year survival rate for lung cancer death was significantly poorer for males (21.3%) than females (23.6%). Subjects with squamous cell carcinoma (SQCC) and treatment by surgical resection alone had better prognosis. We find surgical resections to markedly increase 5-year survival rate from LCC, decreased risk of death from LCC, and no improved survival from SCC.</p> <p>Conclusion</p> <p>Gender and clinical characteristics (i.e. diagnostic period, diagnostic age, histological type and treatment modality) play important roles in determining lung cancer survival.</p

The efficacy of four-slice helical CT in evaluating pancreatic trauma: a single institution experience

<p>Abstract</p> <p>Study objective</p> <p>To assess the efficacy of computed tomography (CT) in evaluating patients with pancreatic trauma.</p> <p>Methods</p> <p>We undertook a retrospective review of all blunt trauma patients admitted to the Chi-Mei Medical Center from January 2004 to June 2006. Every patients underwent abdominal CT scan in emergency department and the CT scans were obtained with a four-slice helical CT. Diagnosis of a pancreatic injury in these patients was by surgical observation or by CT findings. Radiographic pancreatic injuries were classified as deep or superficial lesions. Deep lesions were defined as the hematomas or lacerations >50% thickness of the pancreas. Superficial lesions were described as the hematomas or lacerations <50% thickness of the pancreas; pancreatic edema; and focal fluid accumulation around the pancreas</p> <p>Results</p> <p>Nineteen patients with pancreatic trauma, fourteen males and five females, average age 40.6 ± 21.4 years, were included. Most patients (73.7%) with pancreatic trauma had associated organ injuries. CT was performed in all patients and laparotomy in 14 patients. CT was 78.9% sensitive in detecting pancreatic trauma. All deep pancreatic lesions revealed on CT required surgical treatment, and complication was discovered in two patients undergoing delayed surgery. Superficial lesions were managed conservatively.</p> <p>Conclusion</p> <p>Four-slice helical CT can detect most pancreatic trauma and provide practical therapeutic guidance. Delayed operation might result in complications and is associated with prolonged hospital stays.</p

The Role of Imported Cases and Favorable Meteorological Conditions in the Onset of Dengue Epidemics

Dengue/dengue hemorrhagic fever is the world's most widely spread mosquito-borne arboviral disease and threatens more than two-thirds of the world's population. Cases are mainly distributed in tropical and subtropical areas in accordance with vector habitats for Aedes aegypti and Ae. albopictus. However, the role of imported cases and favorable meteorological conditions has not yet been quantitatively assessed. This study verified the correlation between the occurrence of indigenous dengue and imported cases in the context of weather variables (temperature, rainfall, relative humidity, etc.) for different time lags in southern Taiwan. Our findings imply that imported cases have a role in igniting indigenous outbreaks, in non-endemics areas, when favorable weather conditions are present. This relationship becomes insignificant in the late phase of local dengue epidemics. Therefore, early detection and case management of imported cases through timely surveillance and rapid laboratory-diagnosis may avert large scale epidemics of dengue/dengue hemorrhagic fever. An early-warning surveillance system integrating meteorological data will be an invaluable tool for successful prevention and control of dengue, particularly in non-endemic countries

Identifying Modules of Coexpressed Transcript Units and Their Organization of Saccharopolyspora erythraea from Time Series Gene Expression Profiles

BACKGROUND: The Saccharopolyspora erythraea genome sequence was released in 2007. In order to look at the gene regulations at whole transcriptome level, an expression microarray was specifically designed on the S. erythraea strain NRRL 2338 genome sequence. Based on these data, we set out to investigate the potential transcriptional regulatory networks and their organization. METHODOLOGY/PRINCIPAL FINDINGS: In view of the hierarchical structure of bacterial transcriptional regulation, we constructed a hierarchical coexpression network at whole transcriptome level. A total of 27 modules were identified from 1255 differentially expressed transcript units (TUs) across time course, which were further classified in to four groups. Functional enrichment analysis indicated the biological significance of our hierarchical network. It was indicated that primary metabolism is activated in the first rapid growth phase (phase A), and secondary metabolism is induced when the growth is slowed down (phase B). Among the 27 modules, two are highly correlated to erythromycin production. One contains all genes in the erythromycin-biosynthetic (ery) gene cluster and the other seems to be associated with erythromycin production by sharing common intermediate metabolites. Non-concomitant correlation between production and expression regulation was observed. Especially, by calculating the partial correlation coefficients and building the network based on Gaussian graphical model, intrinsic associations between modules were found, and the association between those two erythromycin production-correlated modules was included as expected. CONCLUSIONS: This work created a hierarchical model clustering transcriptome data into coordinated modules, and modules into groups across the time course, giving insight into the concerted transcriptional regulations especially the regulation corresponding to erythromycin production of S. erythraea. This strategy may be extendable to studies on other prokaryotic microorganisms

A simple approach to ranking differentially expressed gene expression time courses through Gaussian process regression.

BACKGROUND: The analysis of gene expression from time series underpins many biological studies. Two basic forms of analysis recur for data of this type: removing inactive (quiet) genes from the study and determining which genes are differentially expressed. Often these analysis stages are applied disregarding the fact that the data is drawn from a time series. In this paper we propose a simple model for accounting for the underlying temporal nature of the data based on a Gaussian process. RESULTS: We review Gaussian process (GP) regression for estimating the continuous trajectories underlying in gene expression time-series. We present a simple approach which can be used to filter quiet genes, or for the case of time series in the form of expression ratios, quantify differential expression. We assess via ROC curves the rankings produced by our regression framework and compare them to a recently proposed hierarchical Bayesian model for the analysis of gene expression time-series (BATS). We compare on both simulated and experimental data showing that the proposed approach considerably outperforms the current state of the art. CONCLUSIONS: Gaussian processes offer an attractive trade-off between efficiency and usability for the analysis of microarray time series. The Gaussian process framework offers a natural way of handling biological replicates and missing values and provides confidence intervals along the estimated curves of gene expression. Therefore, we believe Gaussian processes should be a standard tool in the analysis of gene expression time series