A modular approach for integrative analysis of large-scale gene-expression and drug-response data

High-throughput technologies are now used to generate more than one type of data from the same biological samples. To properly integrate such data, we propose using co-modules, which describe coherent patterns across paired data sets, and conceive several modular methods for their identification. We first test these methods using in silico data, demonstrating that the integrative scheme of our Ping-Pong Algorithm uncovers drug-gene associations more accurately when considering noisy or complex data. Second, we provide an extensive comparative study using the gene-expression and drug-response data from the NCI-60 cell lines. Using information from the DrugBank and the Connectivity Map databases we show that the Ping-Pong Algorithm predicts drug-gene associations significantly better than other methods. Co-modules provide insights into possible mechanisms of action for a wide range of drugs and suggest new targets for therapy

Meta-analysis of genome-wide association data identifies a risk locus for major mood disorders on 3p21.1.

The major mood disorders, which include bipolar disorder and major depressive disorder (MDD), are considered heritable traits, although previous genetic association studies have had limited success in robustly identifying risk loci. We performed a meta-analysis of five case-control cohorts for major mood disorder, including over 13,600 individuals genotyped on high-density SNP arrays. We identified SNPs at 3p21.1 associated with major mood disorders (rs2251219, P = 3.63 x 10(-8); odds ratio = 0.87; 95% confidence interval, 0.83-0.92), with supportive evidence for association observed in two out of three independent replication cohorts. These results provide an example of a shared genetic susceptibility locus for bipolar disorder and MDD.

Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data.

Despite the central role of quantitative PCR (qPCR) in the quantification of mRNA transcripts, most analyses of qPCR data are still delegated to the software that comes with the qPCR apparatus. This is especially true for the handling of the fluorescence baseline. This article shows that baseline estimation errors are directly reflected in the observed PCR efficiency values and are thus propagated exponentially in the estimated starting concentrations as well as 'fold-difference' results. Because of the unknown origin and kinetics of the baseline fluorescence, the fluorescence values monitored in the initial cycles of the PCR reaction cannot be used to estimate a useful baseline value. An algorithm that estimates the baseline by reconstructing the log-linear phase downward from the early plateau phase of the PCR reaction was developed and shown to lead to very reproducible PCR efficiency values. PCR efficiency values were determined per sample by fitting a regression line to a subset of data points in the log-linear phase. The variability, as well as the bias, in qPCR results was significantly reduced when the mean of these PCR efficiencies per amplicon was used in the calculation of an estimate of the starting concentration per sample.

Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Alcohol-related and hepatocellular cancer deaths by country of birth in England and Wales: analysis of mortality and census data.

BACKGROUND: The incidence of and mortality from alcohol-related conditions, liver disease and hepatocellular cancer (HCC) are increasing in the UK. We compared mortality rates by country of birth to explore potential inequalities and inform clinical and preventive care. DESIGN: Analysis of mortality for people aged 20 years and over using the 2001 Census data and death data from 1999 and 2001-2003. SETTING: England and Wales. MAIN OUTCOME MEASURES: Standardized mortality ratios (SMRs) for alcohol-related deaths and HCC. RESULTS: Mortality from alcohol-related deaths (23 502 deaths) was particularly high for people born in Ireland (SMR for men [M]: 236, 95% confidence interval [CI]: 219-254; SMR for women [F]: 212, 95% CI: 191-235) and Scotland (SMR-M: 187, CI: 173-213; SMR-F 182, CI: 163-205) and men born in India (SMR-M: 161, CI: 144-181). Low alcohol-related mortality was found in women born in other countries and men born in Bangladesh, Middle East, West Africa, Pakistan, China and Hong Kong, and the West Indies. Similar mortality patterns were observed by country of birth for alcoholic liver disease and other liver diseases. Mortality from HCC (8266 deaths) was particularly high for people born in Bangladesh (SMR-M: 523, CI: 380-701; SMR-F: 319, CI: 146-605), China and Hong Kong (SMR-M: 492, CI: 168-667; SMR-F: 323, CI: 184-524), West Africa (SMR-M: 440, CI, 308-609; SMR-F: 319, CI: 165-557) and Pakistan (SMR-M: 216, CI: 113-287; SMR-F: 215, CI: 133-319). CONCLUSIONS: These findings show persistent differences in mortality by country of birth for both alcohol-related and HCC deaths and have important clinical and public health implications. New policy, research and practical action are required to address these differences.This resource was contributed by The National Documentation Centre on Drug Use.

Spatial distribution of schistosomiasis foci on Itamaracá Island, Pernambuco, Brazil

Acute cases of schistosomiasis have been found on the coastal area of Pernambuco, Brazil, due to environmental disturbances and disorderly occupation of the urban areas. This study identifies and spatially marks the main foci of the snail host species, Biomphalaria glabrata on Itamaracá Island. The chaotic occupation of the beach resorts has favoured the emergence of transmission foci, thus exposing residents and tourists to the risk of infection. A database covering five years of epidemiological investigation on snails infected by Schistosoma mansoni in the island was produced with information from the geographic positioning of the foci, number of snails collected, number of snails tested positive, and their infection rate. The spatial position of the foci were recorded through the Global Positioning System (GPS), and the geographical coordinates were imported by AutoCad. The software packages ArcView and Spring were used for data processing and spatial analysis. AutoCad 2000 was used to plot the pairs of coordinates obtained from GPS. Between 1998 and 2002 5009 snails, of which 12.2% were positive for S. mansoni, were collected in Forte Beach. A total of 27 foci and areas of environmental risk were identified and spatially analyzed allowing the identification of the areas exposed to varying degrees of risk.

Diabetes in the United Kingdom: Analysis of QRESEARCH Data

Report describing 10 year trends in incidence, prevalence and mortality of diabetes and obesity by age, sex, deprivation and ethnicity from this Primary Care database.

Capillaria spp. eggs in Patagonian archaeological sites: statistical analysis of morphometric data

Discriminant analysis was used to identify eggs of Capillaria spp. at specific level found in organic remains from an archaeological site in Patagonia, Argentina, dated of 6,540 ± 110 years before present. In order to distinguish eggshell morphology 149 eggs were measured and grouped into four arbitrary subsets. The analysis used on egg width and length discriminated them into different morphotypes (Wilks' lambda = 0.381, p < 0.05). The correlation analysis suggests that width was the most important variable to discriminate among the Capillaria spp. egg morphotypes (Pearson coefficient = 0.950, p < 0.05). The study of eggshell patterns, the relative frequency in the sample, and the morphometric data allowed us to correlate the four morphotypes with Capillaria species.

Thermomechanical analysis of strain localization in a ductile detachment zone

At the latitude of the Thor-Odin dome (British Columbia) the Columbia River Detachment defines the eastern margin of the Shuswap metamorphic core complex and localizes in a 1 km thick muscovite-bearing quartzite mylonite. We present a combined Ar-40/Ar-39, (micro) structural, and oxygen isotope study of the deformation history in the detachment and evaluate the spatial and temporal relationships between microstructure formation and localization of strain. High-precision Ar-40/Ar-39 geochronology from different levels in the mylonite delineates a pattern of increasingly younger (49.0 to 47.9 Ma) deformation ages in deeper levels of the mylonitic footwall. The correlation of Ar-40/Ar-39 ages with decreasing deformation temperatures (similar to 550 degrees-400 degrees C) in the top 200 m of the mylonite indicates that deformation migrated downward from the contact with the hanging wall. Strain localization was diachronous in progressively deeper levels of the footwall and was likely controlled by fluid-assisted strain hardening due to advective heat removal and contemporaneous reaction weakening due to dissolution-reprecipitation of white mica. The observed constant high-stress microstructures across the entire detachment indicate that flow stress was buffered by the interplay of strain rate and temperature, where high strain rates at elevated temperature produced the same microstructure as lower strain rates under decreasing temperature conditions. The combined data suggest that the complex interplay among temporally nonuniform rates of footwall exhumation, heat advection, and embrittlement by meteoric fluids strongly determines the thermomechanical behavior of extensional detachments.

Cluster analysis of behavioural and event-related potentials during a contingent negative variation paradigm in remitting-relapsing and benign forms of multiple sclerosis

Background: Event-related potentials (ERPs) may be used as a highly sensitive way of detecting subtle degrees of cognitive dysfunction. On the other hand, impairment of cognitive skills is increasingly recognised as a hallmark of patients suffering from multiple sclerosis (MS). We sought to determine the psychophysiological pattern of information processing among MS patients with the relapsing-remitting form of the disease and low physical disability considered as two subtypes: 'typical relapsing-remitting' (RRMS) and 'benign MS' (BMS). Furthermore, we subjected our data to a cluster analysis to determine whether MS patients and healthy controls could be differentiated in terms of their psychophysiological profile.Methods: We investigated MS patients with RRMS and BMS subtypes using event-related potentials (ERPs) acquired in the context of a Posner visual-spatial cueing paradigm. Specifically, our study aimed to assess ERP brain activity in response preparation (contingent negative variation -CNV) and stimuli processing in MS patients. Latency and amplitude of different ERP components (P1, eN1, N1, P2, N2, P3 and late negativity -LN) as well as behavioural responses (reaction time -RT; correct responses -CRs; and number of errors) were analyzed and then subjected to cluster analysis. Results: Both MS groups showed delayed behavioural responses and enhanced latency for long-latency ERP components (P2, N2, P3) as well as relatively preserved ERP amplitude, but BMS patients obtained more important performance deficits (lower CRs and higher RTs) and abnormalities related to the latency (N1, P3) and amplitude of ERPs (eCNV, eN1, LN). However, RRMS patients also demonstrated abnormally high amplitudes related to the preparation performance period of CNV (cCNV) and post-processing phase (LN). Cluster analyses revealed that RRMS patients appear to make up a relatively homogeneous group with moderate deficits mainly related to ERP latencies, whereas BMS patients appear to make up a rather more heterogeneous group with more severe information processing and attentional deficits. Conclusions: Our findings are suggestive of a slowing of information processing for MS patients that may be a consequence of demyelination and axonal degeneration, which also seems to occur in MS patients that show little or no progression in the physical severity of the disease over time.

Compositional data analysis of archaeological glass: problems and possible solutions

Presentation in CODAWORK'03, session 4: Applications to archeometry

Compositional ideas in the bayesian analysis of categorical data with application to dose finding clinical trials

Compositional random vectors are fundamental tools in the Bayesian analysis of categorical data.Many of the issues that are discussed with reference to the statistical analysis of compositionaldata have a natural counterpart in the construction of a Bayesian statistical model for categoricaldata.This note builds on the idea of cross-fertilization of the two areas recommended by Aitchison (1986)in his seminal book on compositional data. Particular emphasis is put on the problem of whatparameterization to use

Clinical data and molecular analysis of Mycobacterium tuberculosis isolates from drug-resistant tuberculosis patients in Goiás, Brazil

Drug resistance is one of the major concerns regarding tuberculosis (TB) infection worldwide because it hampers control of the disease. Understanding the underlying mechanisms responsible for drug resistance development is of the highest importance. To investigate clinical data from drug-resistant TB patients at the Tropical Diseases Hospital, Goiás (GO), Brazil and to evaluate the molecular basis of rifampin (R) and isoniazid (H) resistance in Mycobacterium tuberculosis. Drug susceptibility testing was performed on 124 isolates from 100 patients and 24 isolates displayed resistance to R and/or H. Molecular analysis of drug resistance was performed by partial sequencing of the rpoB and katGgenes and analysis of the inhA promoter region. Similarity analysis of isolates was performed by 15 loci mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing. The molecular basis of drug resistance among the 24 isolates from 16 patients was confirmed in 18 isolates. Different susceptibility profiles among the isolates from the same individual were observed in five patients; using MIRU-VNTR, we have shown that those isolates were not genetically identical, with differences in one to three loci within the 15 analysed loci. Drug-resistant TB in GO is caused by M. tuberculosis strains with mutations in previously described sites of known genes and some patients harbour a mixed phenotype infection as a consequence of a single infective event; however, further and broader investigations are needed to support our findings.