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

One person with two DNA profiles: a(nother) case of mosaicism or chimerism.

Nuclear DNA markers, such as short tandem repeats (STR), are widely used for crime investigation and paternity testing. STR were used to determine whether a piece of tissue regurgitated by a dog was part of the penis of a dead, emasculated, man. Unexpectedly, when analyzing the recovered material and a blood sample from the deceased, five out of the 18 loci differed. According to the results, one could have concluded that these samples originated from two different persons. However, taking into account contextual information and data from complementary genetic analyses, the most likely hypothesis was that the deceased was a genetic mosaic or a chimera. Within a forensic genetic context, such genetic peculiarities may prevent associating the perpetrator of an offense with a stain left at a crime scene or lead to false paternity exclusions. Fast recognition of mosaics or chimeras, adapted sampling scheme, as well as careful interpretation of the data should allow avoiding such pitfalls.

Universal features of personality traits from the observer's perspective: Data from 50 cultures

To test hypotheses about the universality of personality traits, college students in 50 cultures identified an adult or college-aged man or woman whom they knew well and rated the 11,985 targets using the 3rd-person version of the Revised NEO Personality Inventory. Factor analyses within cultures showed that the normative American self-report structure was clearly replicated in most cultures and was recognizable in all. Sex differences replicated earlier self-report results, with the most pronounced differences in Western cultures. Cross-sectional age differences for 3 factors followed the pattern identified in self-reports, with moderate rates of change during college age and slower changes after age 40. With a few exceptions, these data support the hypothesis that features of personality traits are common to all human groups.

Plasma and urine profiles of Delta9-tetrahydrocannabinol and its metabolites 11-hydroxy-Delta9-tetrahydrocannabinol and 11-nor-9-carboxy-Delta9-tetrahydrocannabinol after cannabis smoking by male volunteers to estimate recent consumption by athletes.

Since 2004, cannabis has been prohibited by the World Anti-Doping Agency for all sports competitions. In the years since then, about half of all positive doping cases in Switzerland have been related to cannabis consumption. In doping urine analysis, the target analyte is 11-nor-9-carboxy-Delta(9)-tetrahydrocannabinol (THC-COOH), the cutoff being 15 ng/mL. However, the wide urinary detection window of the long-term metabolite of Delta(9)-tetrahydrocannabinol (THC) does not allow a conclusion to be drawn regarding the time of consumption or the impact on the physical performance. The purpose of the present study on light cannabis smokers was to evaluate target analytes with shorter urinary excretion times. Twelve male volunteers smoked a cannabis cigarette standardized to 70 mg THC per cigarette. Plasma and urine were collected up to 8 h and 11 days, respectively. Total THC, 11-hydroxy-Delta(9)-tetrahydrocannabinol (THC-OH), and THC-COOH were determined after hydrolysis followed by solid-phase extraction and gas chromatography/mass spectrometry. The limits of quantitation were 0.1-1.0 ng/mL. Eight puffs delivered a mean THC dose of 45 mg. Plasma levels of total THC, THC-OH, and THC-COOH were measured in the ranges 0.2-59.1, 0.1-3.9, and 0.4-16.4 ng/mL, respectively. Peak concentrations were observed at 5, 5-20, and 20-180 min. Urine levels were measured in the ranges 0.1-1.3, 0.1-14.4, and 0.5-38.2 ng/mL, peaking at 2, 2, and 6-24 h, respectively. The times of the last detectable levels were 2-8, 6-96, and 48-120 h. Besides high to very high THC-COOH levels (245 +/- 1,111 ng/mL), THC (3 +/- 8 ng/mL) and THC-OH (51 +/- 246 ng/mL) were found in 65 and 98% of cannabis-positive athletes' urine samples, respectively. In conclusion, in addition to THC-COOH, the pharmacologically active THC and THC-OH should be used as target analytes for doping urine analysis. In the case of light cannabis use, this may allow the estimation of more recent consumption, probably influencing performance during competitions. However, it is not possible to discriminate the intention of cannabis use, i.e., for recreational or doping purposes. Additionally, pharmacokinetic data of female volunteers are needed to interpret cannabis-positive doping cases of female athletes.

Statistical discrimination of steroid profiles in doping control with support vector machines.

Due to their performance enhancing properties, use of anabolic steroids (e.g. testosterone, nandrolone, etc.) is banned in elite sports. Therefore, doping control laboratories accredited by the World Anti-Doping Agency (WADA) screen among others for these prohibited substances in urine. It is particularly challenging to detect misuse with naturally occurring anabolic steroids such as testosterone (T), which is a popular ergogenic agent in sports and society. To screen for misuse with these compounds, drug testing laboratories monitor the urinary concentrations of endogenous steroid metabolites and their ratios, which constitute the steroid profile and compare them with reference ranges to detect unnaturally high values. However, the interpretation of the steroid profile is difficult due to large inter-individual variances, various confounding factors and different endogenous steroids marketed that influence the steroid profile in various ways. A support vector machine (SVM) algorithm was developed to statistically evaluate urinary steroid profiles composed of an extended range of steroid profile metabolites. This model makes the interpretation of the analytical data in the quest for deviating steroid profiles feasible and shows its versatility towards different kinds of misused endogenous steroids. The SVM model outperforms the current biomarkers with respect to detection sensitivity and accuracy, particularly when it is coupled to individual data as stored in the Athlete Biological Passport.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Profiles of drug users in Switzerland and effects of early-onset intensive use of alcohol, tobacco and cannabis on other illicit drug use.

QUESTIONS UNDER STUDY / PRINCIPLES: The main aim of this study was to investigate profiles of drug users, with a particular focus on illicit drugs other than cannabis, and to explore the effect of early-onset intensive use (drunkenness, daily smoking, high on cannabis) on profiles of illicit drug use. METHODS: Baseline data from a representative sample of 5,831 young Swiss men in the ongoing Cohort Study on Substance Use Risk Factors were used. Substance use (alcohol, tobacco, cannabis and 15 types of other illicit drug) and age of onset of intensive use were assessed. The Item Response Theory (IRT) and prevalence rates at different ages of onset were used to reveal different profiles of illicit drug use. RESULTS: In addition to cannabis, there were two profiles of other illicit drug use: (a) "softer" drug users (uppers, hallucinogens and inhaled drugs), among which ecstasy had the highest discriminatory potential (IRT slope = 4.68, standard error (SE) = 0.48; p <0.001); and (b) "harder" drug users (heroin, ketamine, gamma-hydroxybutyrate/gamma-hydroxylactone, research chemicals, crystal meth and spice), among which ketamine had the highest discriminatory potential (slope = 4.05; SE = 0.63; p <0.001). Onset of intensive use at the age of 12 years or younger also discriminated between these two profiles. CONCLUSION: Both the IRT model and the effect of onset of intensive use enabled two groups of illicit drugs to be identified. In particular, very early onset (at 12 years or younger) intensive use of any substance was a marker for later use of the second group of drugs.

Inferring transport characteristics in a fractured rock aquifer by combining single-hole ground-penetrating radar reflection monitoring and tracer test data

Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.

Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures.

INTRODUCTION: Breast cancer subtyping and prognosis have been studied extensively by gene expression profiling, resulting in disparate signatures with little overlap in their constituent genes. Although a previous study demonstrated a prognostic concordance among gene expression signatures, it was limited to only one dataset and did not fully elucidate how the different genes were related to one another nor did it examine the contribution of well-known biological processes of breast cancer tumorigenesis to their prognostic performance. METHOD: To address the above issues and to further validate these initial findings, we performed the largest meta-analysis of publicly available breast cancer gene expression and clinical data, which are comprised of 2,833 breast tumors. Gene coexpression modules of three key biological processes in breast cancer (namely, proliferation, estrogen receptor [ER], and HER2 signaling) were used to dissect the role of constituent genes of nine prognostic signatures. RESULTS: Using a meta-analytical approach, we consolidated the signatures associated with ER signaling, ERBB2 amplification, and proliferation. Previously published expression-based nomenclature of breast cancer 'intrinsic' subtypes can be mapped to the three modules, namely, the ER-/HER2- (basal-like), the HER2+ (HER2-like), and the low- and high-proliferation ER+/HER2- subtypes (luminal A and B). We showed that all nine prognostic signatures exhibited a similar prognostic performance in the entire dataset. Their prognostic abilities are due mostly to the detection of proliferation activity. Although ER- status (basal-like) and ERBB2+ expression status correspond to bad outcome, they seem to act through elevated expression of proliferation genes and thus contain only indirect information about prognosis. Clinical variables measuring the extent of tumor progression, such as tumor size and nodal status, still add independent prognostic information to proliferation genes. CONCLUSION: This meta-analysis unifies various results of previous gene expression studies in breast cancer. It reveals connections between traditional prognostic factors, expression-based subtyping, and prognostic signatures, highlighting the important role of proliferation in breast cancer prognosis.

Inversion of multiple intersecting high-resolution crosshole GPR profiles for hydrological characterization at the Boise Hydrogeophysical Research Site

The integration of geophysical data into the subsurface characterization problem has been shown in many cases to significantly improve hydrological knowledge by providing information at spatial scales and locations that is unattainable using conventional hydrological measurement techniques. In particular, crosshole ground-penetrating radar (GPR) tomography has shown much promise in hydrology because of its ability to provide highly detailed images of subsurface radar wave velocity, which is strongly linked to soil water content. Here, we develop and demonstrate a procedure for inverting together multiple crosshole GPR data sets in order to characterize the spatial distribution of radar wave velocity below the water table at the Boise Hydrogeophysical Research Site (BHRS) near Boise, Idaho, USA. Specifically, we jointly invert 31 intersecting crosshole GPR profiles to obtain a highly resolved and consistent radar velocity model along the various profile directions. The model is found to be strongly correlated with complementary neutron porosity-log data and is further corroborated by larger-scale structural information at the BHRS. This work is an important prerequisite to using crosshole GPR data together with existing hydrological measurements for improved groundwater flow and contaminant transport modeling.

Beta-catenin Status in P?aediatric Medulloblastomas: correlation of immunohistochemical expression with mutational status, genetic profiles, and clinical characteristics.

Medulloblastoma is the most frequent malignant paediatric brain tumour. The activation of the Wnt/beta-catenin pathway occurs in 10-15% of medulloblastomas and has been recently described as a marker for favourable patient outcome. We report a series of 72 paediatric medulloblastomas evaluated for beta-catenin protein expression, CTNNB1 mutations, and comparative genomic hybridization. Gene expression profiles were also available in a subset of 40 cases. Immunostaining of beta-catenin showed extensive nuclear staining (>50% of the tumour cells) in six cases and focal nuclear staining (<10% of cells) in three cases. The other cases either exhibited a signal strictly limited to the cytoplasm (58 cases) or were negative (five cases). CTNNB1 mutations were detected in all beta-catenin extensively nucleopositive cases. The expression profiles of these cases documented strong activation of the Wnt/beta-catenin pathway. Remarkably, five out of these six tumours showed a complete loss of chromosome 6. In contrast, cases with focal nuclear beta-catenin staining, as well as tumours with negative or cytoplasmic staining, never demonstrated CTNNB1 mutation, Wnt/beta-catenin pathway activation or chromosome 6 loss. Patients with extensive nuclear staining were significantly older at diagnosis and were in continuous complete remission after a mean follow-up of 75.7 months (range 27.5-121.2 months) from diagnosis. All three patients with focal nuclear staining of beta-catenin died within 36 months from diagnosis. Altogether, these data confirm and extend previous observations that CTNNB1-mutated tumours represent a distinct molecular subgroup of medulloblastomas with favourable outcome, indicating that therapy de-escalation should be considered. International consensus on the definition criteria of this distinct medulloblastoma subgroup should be achieved.

Genome-wide RNA polymerase II profiles and RNA accumulation reveal kinetics of transcription and associated epigenetic changes during diurnal cycles.

Interactions of cell-autonomous circadian oscillators with diurnal cycles govern the temporal compartmentalization of cell physiology in mammals. To understand the transcriptional and epigenetic basis of diurnal rhythms in mouse liver genome-wide, we generated temporal DNA occupancy profiles by RNA polymerase II (Pol II) as well as profiles of the histone modifications H3K4me3 and H3K36me3. We used these data to quantify the relationships of phases and amplitudes between different marks. We found that rhythmic Pol II recruitment at promoters rather than rhythmic transition from paused to productive elongation underlies diurnal gene transcription, a conclusion further supported by modeling. Moreover, Pol II occupancy preceded mRNA accumulation by 3 hours, consistent with mRNA half-lives. Both methylation marks showed that the epigenetic landscape is highly dynamic and globally remodeled during the 24-hour cycle. While promoters of transcribed genes had tri-methylated H3K4 even at their trough activity times, tri-methylation levels reached their peak, on average, 1 hour after Pol II. Meanwhile, rhythms in tri-methylation of H3K36 lagged transcription by 3 hours. Finally, modeling profiles of Pol II occupancy and mRNA accumulation identified three classes of genes: one showing rhythmicity both in transcriptional and mRNA accumulation, a second class with rhythmic transcription but flat mRNA levels, and a third with constant transcription but rhythmic mRNAs. The latter class emphasizes widespread temporally gated posttranscriptional regulation in the mouse liver.

Methadone maintenance treatment in the Swiss Canton of Vaud: demographic and clinical data on 1,782 ambulatory patients.

Using data from the Public Health Service, we studied the demographic and clinical characteristics of 1,782 patients enrolled in methadone maintenance treatment (MMT) during 2001 in the Swiss Canton of Vaud, comparing our findings with the results of a previous study from 1976 to 1986. In 2001, most patients (76.9%) were treated in general practice. Mortality is low in this MMT population (1%/year). While patient age and sex profiles were similar to those found in the earlier study, we did observe a substantial increase in the number of patients and the number of practitioners treating MMT patients, probably reflecting the low-threshold governmental policies and the creation of specialized centers. In conclusion, easier access to MMT enhances the number of patients, but new concerns about the quality of management emerge: benzodiazepine as a concomitant prescription; low rates of screening for hepatitis B, C and HIV, and social and psychiatric preoccupations.

Disease profiles of detainees in the Canton of Vaud in Switzerland: gender and age differences in substance abuse, mental health and chronic health conditions.

BACKGROUND: Literature on the disease profile of prisoners that differentiates by age and gender remains sparse. This study aimed to describe the health of correctional inmates in terms of substance abuse problems and mental and somatic health conditions, and compare it by gender and age. METHODS: This study examined cross-sectional data from the Canton of Vaud in Switzerland on the health conditions of detainees who were in prison on January 1, 2011 or entered prison in 2011. Health conditions validated by physician examination were reported using the International Classification of Diseases (ICD) version 10. The analyses were descriptive by groups of prisoners: the entire sample (All), Men, Older adults and Women. RESULTS: A total of 1,664 individuals were included in the analysis. Men comprised 91.5 % of the sample and had a mean age of 33 years. The other 8.5 % were women and had an average age of 39. Older adults (i.e., age 50 and older) represented 7 % of the total sample. Overall, 80 % of inmates were non-Swiss citizens, but the proportion of Swiss prisoners was higher among the older adults (51 %) and women (29 %). Overall, 41 % of inmates self-reported substance abuse problems. Of those, 27 % were being treated by psychiatrists for behavioral disorders related to substance abuse. Chronic infectious diseases were found in 9 % of the prison population. In addition, 27 % of detainees suffered from serious mental health conditions. Gender and age had an influence on the disease profile of this sample: compared to the entire prison population, the older inmates were less likely to misuse illegal drugs and to suffer from communicable infections but exhibited more problems with alcohol and a higher burden of chronic health conditions. Female prisoners were more disposed to mental health problems (including drug abuse) and infectious diseases. In terms of chronic diseases, women suffered from the same conditions as men, but the diseases were more prevalent in women. CONCLUSION: It is important to understand the different disease profiles of prisoners by gender and age, as it helps identify the needs of different groups and tailor age-and gender-specific interventions.