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

Trends in drug consumption and risk of transmission of HIV and hepatitis C virus among injecting drug users in Switzerland, 1993-2006

As a part of the HIV behavioural surveillance system in Switzerland, repeated cross-sectional surveys were conducted in 1993, 1994, 1996, 2000 and 2006 among attenders of all low threshold facilities (LTFs) with needle exchange programmes and/or supervised drug consumption rooms for injection or inhalation in Switzerland. Data were collected in each LTF over five consecutive days, using a questionnaire that was partly completed by an interviewer and partly self administered. The questionnaire was structured around three topics: socio-demographic characteristics, drug consumption, health and risk/preventive behaviour. Analysis was restricted to attenders who had injected drugs during their lifetime (IDUs). Between 1993 and 2006, the median age of IDUs rose by 10 years. IDUs are severely marginalised and their social situation has improved little. The borrowing of used injection equipment (syringe or needle already used by other person) in the last six months decreased (16.5% in 1993, 8.9% in 2006) but stayed stable at around 10% over the past three surveys. Other risk behaviour, such as sharing spoons, cotton or water, was reported more frequently, although also showed a decreasing trend. The reported prevalence of HIV remained fairly stable at around 10% between 1993 and 2006; reported levels of hepatitis C virus (HCV) prevalence were high (56.4% in 2006). In conclusion, the overall decrease in the practice of injection has reduced the potential for transmission of infections. However as HCV prevalence is high this is of particular concern, as the current behaviour of IDUs indicates a potential for further spreading of the infection. Another noteworthy trend is the significant decrease in condom use in the case of paid sex.

Current genetic data do not improve the prediction of type 2 diabetes mellitus: the CoLaus study.

CONTEXT: Several genetic risk scores to identify asymptomatic subjects at high risk of developing type 2 diabetes mellitus (T2DM) have been proposed, but it is unclear whether they add extra information to risk scores based on clinical and biological data. OBJECTIVE: The objective of the study was to assess the extra clinical value of genetic risk scores in predicting the occurrence of T2DM. DESIGN: This was a prospective study, with a mean follow-up time of 5 yr. SETTING AND SUBJECTS: The study included 2824 nondiabetic participants (1548 women, 52 ± 10 yr). MAIN OUTCOME MEASURE: Six genetic risk scores for T2DM were tested. Four were derived from the literature and two were created combining all (n = 24) or shared (n = 9) single-nucleotide polymorphisms of the previous scores. A previously validated clinic + biological risk score for T2DM was used as reference. RESULTS: Two hundred seven participants (7.3%) developed T2DM during follow-up. On bivariate analysis, no differences were found for all but one genetic score between nondiabetic and diabetic participants. After adjusting for the validated clinic + biological risk score, none of the genetic scores improved discrimination, as assessed by changes in the area under the receiver-operating characteristic curve (range -0.4 to -0.1%), sensitivity (-2.9 to -1.0%), specificity (0.0-0.1%), and positive (-6.6 to +0.7%) and negative (-0.2 to 0.0%) predictive values. Similarly, no improvement in T2DM risk prediction was found: net reclassification index ranging from -5.3 to -1.6% and nonsignificant (P ≥ 0.49) integrated discrimination improvement. CONCLUSIONS: In this study, adding genetic information to a previously validated clinic + biological score does not seem to improve the prediction of T2DM.

Data Envelopment Analysis (DEA) : un guide pédagogique à l'intention des décideurs dans le secteur public

Ce guide présente la méthode Data Envelopment Analysis (DEA), une méthode d'évaluation de la performance . Il est destiné aux responsables d'organisations publiques qui ne sont pas familiers avec les notions d'optimisation mathématique, autrement dit de recherche opérationnelle. L'utilisation des mathématiques est par conséquent réduite au minimum. Ce guide est fortement orienté vers la pratique. Il permet aux décideurs de réaliser leurs propres analyses d'efficience et d'interpréter facilement les résultats obtenus. La méthode DEA est un outil d'analyse et d'aide à la décision dans les domaines suivants : - en calculant un score d'efficience, elle indique si une organisation dispose d'une marge d'amélioration ; - en fixant des valeurs-cibles, elle indique de combien les inputs doivent être réduits et les outputs augmentés pour qu'une organisation devienne efficiente ; - en identifiant le type de rendements d'échelle, elle indique si une organisation doit augmenter ou au contraire réduire sa taille pour minimiser son coût moyen de production ; - en identifiant les pairs de référence, elle désigne quelles organisations disposent des best practice à analyser.

Quantitative genetic modeling and inference in the presence of nonignorable missing data.

Natural selection is typically exerted at some specific life stages. If natural selection takes place before a trait can be measured, using conventional models can cause wrong inference about population parameters. When the missing data process relates to the trait of interest, a valid inference requires explicit modeling of the missing process. We propose a joint modeling approach, a shared parameter model, to account for nonrandom missing data. It consists of an animal model for the phenotypic data and a logistic model for the missing process, linked by the additive genetic effects. A Bayesian approach is taken and inference is made using integrated nested Laplace approximations. From a simulation study we find that wrongly assuming that missing data are missing at random can result in severely biased estimates of additive genetic variance. Using real data from a wild population of Swiss barn owls Tyto alba, our model indicates that the missing individuals would display large black spots; and we conclude that genes affecting this trait are already under selection before it is expressed. Our model is a tool to correctly estimate the magnitude of both natural selection and additive genetic variance.

HIV behavioural surveillance in Switzerland: trends in drug consumption, injection and syringe distribution, 1993-2011

Background: As part of the second generation surveillance system for HIV/Aids in Switzerland, repeated cross-sectional surveys were conducted in 1993, 1994, 1996, 2000, 2006 and 2011 among attenders of all low threshold facilities (LTFs) with needle exchange programmes and/or supervised drug consumption rooms for injection or inhalation. The number of syringes distributed to the injectors has also been measured annually since 2000. Distribution in other settings, such as pharmacies, is also monitored nationally. Methods: Periodic surveys of LTFs have been conducted using an interviewer/self-administered questionnaire structured along four themes: socio-demographic characteristics, drug consumption, risk/preventive behaviour and health. Analysis is restricted to attenders who had injected drugs during their lifetime (IDU´s). Pearson's chi-square test and trend analysis were conducted on annual aggregated data. Trend significance was assessed using Stata's non parametric test nptrend. Results: Median age of IDU´s increased from 26 years in 1993 to 40 in 2011; most are men (78%). Total yearly number of syringes distributed by LTFs has decreased by 44% in 10 years. Use of cocaine has increased (Table 1). Injection, regular use of heroin and borrowing of syringes/needles have decreased, while sharing of other material remains stable. There are fewer new injectors; more IDU´s report substitution treatment. Most attenders had ever been tested for HIV (90% in 1993, 94% in 2011). Reported prevalence of HIV remained stable around 10%; that of HCV decreased from 62% in 2000 to 42% in 2011. Conclusions: Overall, findings indicate a decrease in injection as a means of drug consumption in that population. This interpretation is supported by data from other sources, such as a national decrease in distribution from other delivery points. Switzerland's behavioural surveillance system is sustainable and allows the HIV epidemic to be monitored among this hard-to-reach population, providing information for planning and evaluation.

Regional-scale integration of multi-scale hydrological and geophysical data using a two-step Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale for the purpose of improving predictions of groundwater flow and solute transport. However, extending corresponding approaches to the regional scale still represents one of the major challenges in the domain of hydrogeophysics. To address this problem, we have developed a regional-scale data integration methodology based on a two-step Bayesian sequential simulation approach. Our objective is to generate high-resolution stochastic realizations of the regional-scale hydraulic conductivity field in the common case where there exist spatially exhaustive but poorly resolved measurements of a related geophysical parameter, as well as highly resolved but spatially sparse collocated measurements of this geophysical parameter and the hydraulic conductivity. To integrate this multi-scale, multi-parameter database, we first link the low- and high-resolution geophysical data via a stochastic downscaling procedure. This is followed by relating the downscaled geophysical data to the high-resolution hydraulic conductivity distribution. After outlining the general methodology of the approach, we demonstrate its application to a realistic synthetic example where we consider as data high-resolution measurements of the hydraulic and electrical conductivities at a small number of borehole locations, as well as spatially exhaustive, low-resolution estimates of the electrical conductivity obtained from surface-based electrical resistivity tomography. The different stochastic realizations of the hydraulic conductivity field obtained using our procedure are validated by comparing their solute transport behaviour with that of the underlying ?true? hydraulic conductivity field. We find that, even in the presence of strong subsurface heterogeneity, our proposed procedure allows for the generation of faithful representations of the regional-scale hydraulic conductivity structure and reliable predictions of solute transport over long, regional-scale distances.

Clustering and hot spot detection in socio-economic spatio-temporal data

Distribution of socio-economic features in urban space is an important source of information for land and transportation planning. The metropolization phenomenon has changed the distribution of types of professions in space and has given birth to different spatial patterns that the urban planner must know in order to plan a sustainable city. Such distributions can be discovered by statistical and learning algorithms through different methods. In this paper, an unsupervised classification method and a cluster detection method are discussed and applied to analyze the socio-economic structure of Switzerland. The unsupervised classification method, based on Ward's classification and self-organized maps, is used to classify the municipalities of the country and allows to reduce a highly-dimensional input information to interpret the socio-economic landscape. The cluster detection method, the spatial scan statistics, is used in a more specific manner in order to detect hot spots of certain types of service activities. The method is applied to the distribution services in the agglomeration of Lausanne. Results show the emergence of new centralities and can be analyzed in both transportation and social terms.

The UniProtKB/Swiss-Prot Tox-Prot program: A central hub of integrated venom protein data.

Animal toxins are of interest to a wide range of scientists, due to their numerous applications in pharmacology, neurology, hematology, medicine, and drug research. This, and to a lesser extent the development of new performing tools in transcriptomics and proteomics, has led to an increase in toxin discovery. In this context, providing publicly available data on animal toxins has become essential. The UniProtKB/Swiss-Prot Tox-Prot program (http://www.uniprot.org/program/Toxins) plays a crucial role by providing such an access to venom protein sequences and functions from all venomous species. This program has up to now curated more than 5000 venom proteins to the high-quality standards of UniProtKB/Swiss-Prot (release 2012_02). Proteins targeted by these toxins are also available in the knowledgebase. This paper describes in details the type of information provided by UniProtKB/Swiss-Prot for toxins, as well as the structured format of the knowledgebase.

Interpretation of array comparative genome hybridization data: a major challenge.

The advent and application of high-resolution array-based comparative genome hybridization (array CGH) has led to the detection of large numbers of copy number variants (CNVs) in patients with developmental delay and/or multiple congenital anomalies as well as in healthy individuals. The notion that CNVs are also abundantly present in the normal population challenges the interpretation of the clinical significance of detected CNVs in patients. In this review we will illustrate a general clinical workflow based on our own experience that can be used in routine diagnostics for the interpretation of CNVs.

The effect of patient, provider and financing regulations on the intensity of ambulatory physical therapy episodes : a multilevel analysis based on routinely available data

Background: Many studies have found considerable variations in the resource intensity of physical therapy episodes. Although they have identified several patient-and provider-related factors, few studies have examined their relative explanatory power. We sought to quantify the contribution of patients and providers to these differences and examine how effective Swiss regulations are (nine-session ceiling per prescription and bonus for first treatments). Methods: Our sample consisted of 87,866 first physical therapy episodes performed by 3,365 physiotherapists based on referrals by 6,131 physicians. We modeled the number of visits per episode using a multilevel log linear regression with crossed random effects for physiotherapists and physicians and with fixed effects for cantons. The three-level explanatory variables were patient, physiotherapist and physician characteristics. Results: The median number of sessions was nine (interquartile range 6-13). Physical therapy use increased with age, women, higher health care costs, lower deductibles, surgery and specific conditions. Use rose with the share of nine-session episodes among physiotherapists or physicians, but fell with the share of new treatments. Geographical area had no influence. Most of the variance was explained at the patient level, but the available factors explained only 4% thereof. Physiotherapists and physicians explained only 6% and 5% respectively of the variance, although the available factors explained most of this variance. Regulations were the most powerful factors. Conclusion: Against the backdrop of abundant physical therapy supply, Swiss financial regulations did not restrict utilization. Given that patient-related factors explained most of the variance, this group should be subject to closer scrutiny. Moreover, further research is needed on the determinants of patient demand.

Evolution of life history traits in a natural population of the shrew Crocidura russula

Abstract Life history traits encompass all the decisions concerning fitness an individual is faced with during his life. The study of these traits is crucial to understand the factors shaping the biology of living organisms. Up until now, most of the information on the evolution of life history traits comes from laboratory studies. While these studies are interesting to test the effect of specific parameters, their conclusions are difficult to extrapolate to natural populations. Investigating the evolution of life history traits in natural populations is of great interest. This may be tricky because it requires information on reproduction, survival and morphology of individuals. Mark-recapture methods allow most of this information to be obtained. However, when direct observations of a species are not possible due to its ecology, indirect methods must be used to infer lifetime reproductive success. In this case, molecular markers are particularly helpful in assessing the genetic relationships between individuals and allow the construction of a pedigree. This thesis focuses on a natural population of a small insectivorous mammal, the greater white-toothed shrew, Crocidura russula. Because of its hidden lifestyle, the two complementary techniques mentioned above were combined to gather information on this population. The data were used to explore diverse aspects of evolutionary biology. We demonstrated that the high genetic variance displayed by the species was not maintained by its mating system because this shrew was less monogamous than previously thought. The large genetic diversity was most likely promoted by gene flow from the neighborhood. Dispersal was thus a central topic in this thesis. We showed that dispersal was not driven by inbreeding avoidance. In addition, we did not find any inbreeding depression in the population. Dispersal was promoted by a high number of vacant territories in the population for both sexes, meaning that territory acquisition played an important role in driving dispersal. Moreover, dispersal propensity was shown to have a genetic basis and, once achieved, to have no effect on individual fitness. Body mass was found to be a life history trait strongly influenced by sexual and viability selection in both sexes. Larger individuals had higher access to reproduction through territory acquisition and defense than lighter ones. By contrast, intermediate size individuals were favored by viability selection presumably because of ecological constraints and metabolic costs. Finally, we demonstrated that the majority of the life history traits in our shrew population has the potential to evolve because they maintained substantial amounts of additive genetic variance. Nonetheless, life history traits had no significant heritability due to their high level of nonadditive or environmental variance. Résumé Les traits d'histoire de vie comprennent toutes les décisions auxquelles un individu est confronté au cours de sa vie et qui concernent sa valeur adaptative. L'étude de ces traits est cruciale pour comprendre les facteurs qui façonnent la biologie des êtres vivants. Jusqu'à ce jour, la majorité des informations sur l'évolution des traits d'histoire de vie provient d'études réalisées en laboratoire. Alors que ces études sont intéressantes pour tester l'effet de paramètres spécifiques, leurs conclusions sont difficilement extrapolables aux populations naturelles. Il est particulièrement intéressant d'étudier l'évolution des traits d'histoire de vie dans des populations naturelles. Toutefois, ces études peuvent se révéler difficiles parce qu'elles requièrent des informations sur la reproduction, la survie et la morphologie des individus. Des méthodes de marquage-recapture permettent d'obtenir ces informations. Cependant, lorsque l'écologie de l'espèce rend les obervations directes impossibles, des méthodes indirectes doivent être utilisées pour obtenir le succès reproducteur des individus. Dans ce cas, les marqueurs moléculaires sont particulièrement utiles pour évaluer les relations génétiques entre individus et permettre la construction d'un pedigree. Cette thèse porte sur une population naturelle d'un petit mammifère insectivore, la musaraigne musette, Crocidura russula. Parce que cette espèce présente un mode de vie souterrain, les deux techniques complémentaires mentionnées ci-dessus ont été combinées pour acquérir les informations nécessaires. Les données ont été utilisées pour explorer divers aspects de biologie evolutive. Nous avons montré que la grande quantité de variance génétique trouvée chez cette espèce n'est pas maintenue par son système d'appariement. Celle-ci s'est en effet avérée être moins monogame que ce qui était admis jusqu'ici. Sa grande diversité génétique est plutôt entretenue par le flux de gènes provenant du voisinage. La dispersion a donc été un sujet phare dans cette thèse. Nous avons montré qu'elle n'est pas provoquée par un évitement de la consanguinité et nous n'avons pas trouvé de dépression de consanguité dans notre population. L'acquisition d'un territoire joue par contre un rôle important dans la dispersion. En outre, la dispersion possède une base génétique chez cette espèce. De plus, une fois qu'ils ont dispersé, les individus n'ont pas une valeur adaptative differente d'individus philopatriques. Le poids s'est avéré être un trait d'histoire de vie fortement influencé par la sélection sexuelle et de viabilité chez les deux sexes. Les gros individus ont accès à la reproduction parce qu'ils acquièrent et défendent un territoire plus facilement que les plus légers. Au contraire, les individus de taille intermédiaire sont favorisés par la sélection de viabilité, certainement à cause de contraintes écologiques et de coûts métaboliques. Finalement, nous avons montré que la majorité des traits d'histoire de vie dans notre population a le potentiel d'évoluer parce qu'elle maintient des quantités considérables de variance génétique additive. Néanmoins, l'héritabilité de ces traits d'histoire de vie n'est pas significative à cause de la grande quantité de variance non-additive ou environmentale associée à ces traits.

Technical note: Implementation of a geodatabase of published and unpublished data on the catastrophic Vaiont landslide

On 9 October 1963 a catastrophic landslide suddenly occurred on the southern slope of the Vaiont dam reservoir. A mass of approximately 270 million m3 collapsed into the reservoir generating a wave that overtopped the dam and hit the town of Longarone and other villages nearby. Several investigations and interpretations of the slope collapse have been carried out during the last 45 years, however, a comprehensive explanation of both the triggering and the dynamics of the phenomenon has yet to be provided. In order to re-evaluate the currently existing information on the slide, an electronic bibliographic database and an ESRI-geodatabase have been developed. The chronology of the collected documentation showed that most of the studies for re-evaluating the failure mechanisms were conducted in the last decade, as a consequence of knowledge, methods and techniques recently acquired. The current contents of the geodatabase will improve definition of the structural setting that influenced the slide and led to the the propagation of the displaced rock mass. The objectives, structure and contents of the e-bibliography and Geodatabase are indicated, together with a brief description on the possible use of the alphanumeric and spatial contents of the databases.