An efficient algorithm to perform multiple testing in epistasis screening

Background: Research in epistasis or gene-gene interaction detection for human complex traits has grown over the last few years. It has been marked by promising methodological developments, improved translation efforts of statistical epistasis to biological epistasis and attempts to integrate different omics information sources into the epistasis screening to enhance power. The quest for gene-gene interactions poses severe multiple-testing problems. In this context, the maxT algorithm is one technique to control the false-positive rate. However, the memory needed by this algorithm rises linearly with the amount of hypothesis tests. Gene-gene interaction studies will require a memory proportional to the squared number of SNPs. A genome-wide epistasis search would therefore require terabytes of memory. Hence, cache problems are likely to occur, increasing the computation time. In this work we present a new version of maxT, requiring an amount of memory independent from the number of genetic effects to be investigated. This algorithm was implemented in C++ in our epistasis screening software MBMDR-3.0.3. We evaluate the new implementation in terms of memory efficiency and speed using simulated data. The software is illustrated on real-life data for Crohn’s disease. Results: In the case of a binary (affected/unaffected) trait, the parallel workflow of MBMDR-3.0.3 analyzes all gene-gene interactions with a dataset of 100,000 SNPs typed on 1000 individuals within 4 days and 9 hours, using 999 permutations of the trait to assess statistical significance, on a cluster composed of 10 blades, containing each four Quad-Core AMD Opteron(tm) Processor 2352 2.1 GHz. In the case of a continuous trait, a similar run takes 9 days. Our program found 14 SNP-SNP interactions with a multiple-testing corrected p-value of less than 0.05 on real-life Crohn’s disease (CD) data. Conclusions: Our software is the first implementation of the MB-MDR methodology able to solve large-scale SNP-SNP interactions problems within a few days, without using much memory, while adequately controlling the type I error rates. A new implementation to reach genome-wide epistasis screening is under construction. In the context of Crohn’s disease, MBMDR-3.0.3 could identify epistasis involving regions that are well known in the field and could be explained from a biological point of view. This demonstrates the power of our software to find relevant phenotype-genotype higher-order associations.

Rethinking study and management of agricultural systems for policy design

There is a concern that agriculture will no longer be able to meet, on a global scale, the growing demand for food. Facing such a challenge requires new patterns of thinking in the context of complexity and sustainability sciences. This paper, focused on the social dimension of the study and management of agricultural systems, suggests that rethinking the study of agricultural systems entails analyzing them as complex socio-ecological systems, as well as considering the differing thinking patterns of diverse stakeholders. The intersubjective nature of knowledge, as studied by different philosophical schools, needs to be better integrated into the study and management of agricultural systems than it is done so far, forcing us to accept that there are no simplistic solutions, and to seek a better understanding of the social dimension of agriculture. Different agriculture related problems require different policy and institutional approaches. Finally, the intersubjective nature of knowledge asks for the visualization of different framings and the power relations taking place in the decision-making process. Rethinking management of agricultural systems implies that policy making should be shaped by different principles: learning, flexibility, adaptation, scale-matching, participation, diversity enhancement and precaution hold the promise to significantly improve current standard management procedures.

On the approximation of the subgrid scale for systems of equations

Postprint (published version)

Immunobiology of human papillomavirus infection and vaccination - implications for second generation vaccines.

Prophylactic human papillomavirus (HPV) L1 virus like particle (VLP) vaccines have been shown, in large clinical trials, to be very immunogenic, well-tolerated and highly efficacious against genital disease caused by the vaccine HPV types. However these vaccines, at the present, protect against only two of the 15 oncogenic genital HPV types, they are expensive, delivered by intramuscular injection and require a cold chain. The challenges are to develop cheap, thermo-stable vaccines that can be delivered by non-injectable methods that provide long term (decades) protection at mucosal surfaces to most, if not all, oncogenic HPV types that is as good as the current VLP vaccines. Current approaches include L1 capsomers, L2 protein and peptides, delivery via recombinant L1 bacterial and viral vectors and large-scale VLP production in plants. Rational design and successful development of such vaccines will be based on an understanding of the immune response, and particularly the 'cross talk' between the innate and adaptive responses. This will be central in the development of adjuvants and vaccine formulations that induce the response to provide effective protection.

Psychometric properties of the Marlowe-Crowne Social Desirability Scale in eight African countries and Switzerland

The purpose of this study was to assess the cross-cultural validity of the Marlowe-Crowne Social Desirability scale short form C, in a large sample of French-speaking participants from eight African countries and Switzerland. Exploratory and confirmatory analyses suggested retaining a two-factor structure. Item bias detection according to country was conducted for all 13 items and effect was calculated with R2. For the two-factor solution, 9 items were associated with a negligible effect size, 3 items with a moderate one, and 1 item with a large one. A series of analyses of covariance considering the acquiescence variable as a covariate showed that the acquiescence tendency does not contribute to the bias at item level. This research indicates that the psychometric properties of this instrument do not reach a scalar equivalence but that a culturally reliable measurement of social desirability could be developed.

Scale Effects in Hydraulic Model Tests of Rock Protected Structures, HR-119, 1970

A laboratory investigation was undertaken to determine the limiting model Reynolds number above which the scour behavior of rock protected structures can be reproduced in hydraulic models scaled according to the Froude criterion. A submerged jet was passed over an initially full scour pocket containing uniform glass spheres and the rate of scour was measured as a function of time. The dimensions of the scour pocket and jet and the particle diameters were varied as needed to maintain strict geometric similarity. For each of two different Froude numbers the Reynolds number was varied over a wide range. The normalized scour rate was found to be practically independent of the Reynolds number, R, (based on the jet velocity and particle diameter) at values of R above about 2.5 x 10^3, and to decrease with Rat smaller values. A grid placed in the jet was found to have a very strong effect on the scour rate. In an attempt to explain the effect of R on the scour behavior, turbulent pressure and velocity fluctuations were measured in air flows and water flows, respectively, over rigid scour pockets having the same geometry as those formed in the scour experiments. The normalized spectra of the fluctuations were found to be nearly independent of R, but the flow pattern was found to be very sensitive to the inlet condition, the jet deflecting upward or downward in a not wholly explainable manner. This indicates that scour behavior can be modeled only if the approach flow is also accurately modeled.

A Ten Year Performance Summary of Fibrous PC Concrete Overlay Research in Greene County, Iowa, HR-165, 1984

The Greene County, Iowa, overlay project, completed in October, 1973, was evaluated in October, 1978, after five years of service and most recently in October, 1983, after ten years of service. The 33 fibrous concrete sections, four CRCP sections, two mesh reinforced and two plain concrete sections with doweled reinforcement were rated relative to each other on a scale of 0 to 100. The rating was conducted by original members of the Project Planning Committee, Iowa DOT, Iowa County, Federal Highway Administration, University of Illinois and industry representatives. In all, there were 23 and 24 representatives who rated the project in 1978 and 1983 respectively. The 23 or 24 values were then averaged to provide a final rating number for each section or variable. All experimental overlay sections had performed quite well in the period from five through 10 years, experiencing only limited additional deterioration. Based upon this relatively good performance through 10 years, the sections will be maintained for further research with another evaluation at 15 years. The 4" thick nonfibrous mesh reinforced continuous reinforced concrete pavement overlay sections provided the best performance in this research project. Another nonfibrous 5" thick bar reinforced overlay section performed almost as well. The best performance of a fibrous reinforced concrete section was obtained with 160 pounds of fiber per cubic yard.

The link between land-use management and fluvial flood risk : a chaotic conception?

There is much policy interest in the possible linkages that might exist between land use and downstream fluvial flood risk. On the one hand, this position is sustained by observations from plot- and field-scale studies that suggest land management does affect runoff. On the other, upscaling these effects to show that land-management activities impact upon flood risk at larger catchment scales has proved to be elusive. This review considers the reasons for why this upscaling is problematic. We argue that, rather than it reflecting methodological challenges associated with the difficulties of modelling hydrological processes over very large areas and during extreme runoff events, it reflects the fact that any linkage between land management and flood risk cannot be generalized and taken out of its specific spatial (catchment) and temporal (flood event) context. We use Sayer's (1992) notion of a `chaotic conception' to describe the belief that there is a simple and general association between land management and downstream flood risk rather than the impacts of land management being spatially and temporally contingent in relation to the particular geographical location, time period and scale being considered. Our argument has important practical consequences because it implies that land-management activities to reduce downstream flood risk will be different to traditional flood-reduction interventions such as levees. The purpose of demonstration projects then needs careful consideration such that conclusions made for one project are not transferred uncritically to other scales of analysis or geographical locations.

The birth of the Rheic Ocean - Early Palaeozoic subsidence patterns and subsequent tectonic plate scenarios

New plate-tectonic reconstructions of the Gondwana margin suggest that the location of Gondwana-derived terranes should not only be guided by the models, but should also consider the possible detrital input from some Asian blocks (Hunia), supposed to have been located along the Cambrian Gondwana margin, and accreted in the Silurian to the North-Chinese block. Consequently, the Gondwana margin has to be subdivided into a more western domain, where the future Avalonian blocks will be separated from Gondwana by the opening Rheic Ocean, whereas in its eastern continuation, hosting the future basement areas of Central Europe, different periods of crustal extension should be distinguished. Instead of applying a rather cylindrical model, it is supposed that crustal extension follows a much more complex pattern, where local back-arcs or intra-continental rifts are involved. Guided by the age data of magmatic rocks and the pattern of subsidence curves, the following extensional events can be distinguished: During the early to middle Cambrian, a back-arc setting guided the evolution at the Gondwana margin. Contemporaneous intra-continental rift basins developed at other places related to a general post-PanAfrican extensional phase affecting Africa Upper Cambrian formation of oceanic crust is manifested in the Chamrousse area, and may have lateral cryptic relics preserved in other places. This is regarded as the oceanisation of some marginal basins in a context of back-arc rifting. These basins were closed in a mid-Ordovician tectonic phase, related to the subduction of buoyant material (mid-ocean ridge?) Since the Early Ordovician, a new phase of extension is observed, accompanied by a large-scale volcanic activity, erosion of the rift shoulders generated detritus (Armorican Quartzite) and the rift basins collected detrital zircons from a wide hinterland. This phase heralded the opening of Palaeotethys, but it failed due to the Silurian collision (Eo-Variscan phase) of an intra-oceanic arc with the Gondwana margin. During this time period, at the eastern wing of the Gondwana margin begins the drift of the future Hunia microcontinents, through the opening of an eastern prolongation of the already existing Rheic Ocean. The passive margin of the remaining Gondwana was composed of the Galatian superterranes, constituents of the future Variscan basement areas. Remaining under the influence of crustal extension, they will start their drift to Laurussia since the earliest Devonian during the opening of the Palaeotethys Ocean. (C) 2008 Elsevier B.V. All rights reserved.

A non-linear Bayesian approach for upscaling local-scale hydraulic conductivity measurements based on local- and regional-scale geophysical data

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed an upscaling procedure based on a Bayesian sequential simulation approach. This method is then applied to the stochastic integration of low-resolution, regional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this upscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Climate change and plant distribution: local models predict high-elevation persistence

Mountain ecosystems will likely be affected by global warming during the 21st century, with substantial biodiversity loss predicted by species distribution models (SDMs). Depending on the geographic extent, elevation range and spatial resolution of data used in making these models, different rates of habitat loss have been predicted, with associated risk of species extinction. Few coordinated across-scale comparisons have been made using data of different resolution and geographic extent. Here, we assess whether climate-change induced habitat losses predicted at the European scale (10x10' grid cells) are also predicted from local scale data and modeling (25x25m grid cells) in two regions of the Swiss Alps. We show that local-scale models predict persistence of suitable habitats in up to 100% of species that were predicted by a European-scale model to lose all their suitable habitats in the area. Proportion of habitat loss depends on climate change scenario and study area. We find good agreement between the mismatch in predictions between scales and the fine-grain elevation range within 10x10' cells. The greatest prediction discrepancy for alpine species occurs in the area with the largest nival zone. Our results suggest elevation range as the main driver for the observed prediction discrepancies. Local scale projections may better reflect the possibility for species to track their climatic requirement toward higher elevations.

Mesurer la stigmatisation perçue chez les personnes souffrant de troubles psychiques : traduction française, validation et adaptation de la Stigma Scale [French translation, validation and adaptation of the Stigma Scale].

AIM: People suffering from mental illness are exposed to stigma. However, only few tools are available to assess stigmatization as perceived from the patient's perspective. The aim of this study is to adapt and validate a French version of the Stigma Scale (King et al., 2007 [8]). This self-report questionnaire has a three-factor structure: discrimination, disclosure and positive aspects of mental illness. Discrimination subscale refers to perceived negative reactions of others. Disclosure subscale refers mainly to managing disclosure to avoid discrimination and finally positive aspects subscale taps into how patients are becoming more accepting, more understanding toward their illness. METHOD: In the first step, internal consistency, convergent validity and test-retest reliability of the French adaptation of the 28-item scale were assessed in a sample of 183 patients. Results of confirmatory factor analyses (CFA) did not confirm the hypothesized structure. In the light of the failed attempts to validate the original version, an alternative 9-item short-form version of the Stigma Scale, maintaining the integrity of the original model, was developed based on results of exploratory factor analyses in the first sample and cross-validated in a new sample of 234 patients. RESULTS: Results of CFA did not confirm that the data fitted well to the three-factor model of the 28-item Stigma Scale (χ(2)/df=2.02, GFI=0.77, AGFI=0.73, RMSEA=0.07, CFI=0.77 and NNFI=0.75). Cronbach's α was excellent for discrimination (0.84) and disclosure (0.83) subscales but poor for potential positive aspects (0.46). External validity was satisfactory. Overall Stigma Scale total score was negatively correlated with the score on Rosenberg's Self-Esteem Scale (r=-0.49), and each subscale was significantly correlated with a visual analogue scale that referred to the specific aspect of stigma (0.43≤|r|≤0.60). Intraclass correlation coefficients between 0.68 and 0.89 indicated good test-retest reliability. The results of the CFA demonstrated that the items chosen for the short version of the Stigma Scale have the expected fit properties (χ(2)/df=1.02, GFI=0.98, AGFI=0.98, RMSEA=0.01, CFI=1.0 and NNFI=1.0). Considering the small number (three) of items in each subscale of the short version of the Stigma Scale, α coefficients for discrimination (0.57), disclosure (0.80) and potential positive aspects subscales (0.62) are considered as good. CONCLUSION: Our results suggest that the 9-item French short version of the Stigma Scale is a useful, reliable and valid self-report questionnaire to assess perceived stigmatization in people suffering from mental illness. The time of completion is really short and questions are well understood and accepted by the patients.

Introduction of snow and geomorphic disturbance variables into predictive models of alpine plant distribution in the Western Swiss Alps

Indirect topographic variables have been used successfully as surrogates for disturbance processes in plant species distribution models (SDM) in mountain environments. However, no SDM studies have directly tested the performance of disturbance variables. In this study, we developed two disturbance variables: a geomorphic index (GEO) and an index of snow redistribution by wind (SNOW). These were developed in order to assess how they improved both the fit and predictive power of presenceabsence SDM based on commonly used topoclimatic (TC) variables for 91 plants in the Western Swiss Alps. The individual contribution of the disturbance variables was compared to TC variables. Maps of models were prepared to spatially test the effect of disturbance variables. On average, disturbance variables significantly improved the fit but not the predictive power of the TC models and their individual contribution was weak (5.6% for GEO and 3.3% for SNOW). However their maximum individual contribution was important (24.7% and 20.7%). Finally, maps including disturbance variables (i) were significantly divergent from TC models in terms of predicted suitable surfaces and connectivity between potential habitats, and (ii) were interpreted as more ecologically relevant. Disturbance variables did not improve the transferability of models at the local scale in a complex mountain system, and the performance and contribution of these variables were highly species-specific. However, improved spatial projections and change in connectivity are important issues when preparing projections under climate change because the future range size of the species will determine the sensitivity to changing conditions.

Longitudinal Joint Forming in PCC Pavements, MLR-00-05, Construction Report, 2003

In conventional construction practices, a longitudinal joint is sawed in a PCC (Portland Cement Concrete) pavement to control concrete shrinkage cracking between two lanes of traffic. Sawing a joint in hardened concrete is an expensive and time consuming operation. The longitudinal joint is not a working joint (in comparison to a transverse joint) as it is typically tied with a tie bar at 30 inch spacing. The open joint reservoir, left by the saw blade, typically is filled or sealed with a durable crack sealant to keep incompressibles and water from getting into the joint reservoir. An experimental joint forming knife has been developed. It is installed under the paving machine to form the longitudinal joint in the wet concrete as a part of the paving process. Through this research method, forming a very narrow longitudinal joint during the paving process, two conventional paving operations can be eliminated. Joint forming eliminates the need of the joint sawing operation in the hard concrete, and as the joint that is formed does not leave a wide-open reservoir, but only a hairline crack, it does not need the joint filling or sealing operation. Therefore, the two conventional longitudinal joint sawing and sealing operations are both being eliminated by this innovation. A laboratory scale prototype joint forming knife was built and tested, initially forming joints in small concrete beams. The results were positive so the method was proposed for field testing. Initial field tests were done in the construction season of 2001, limited to one paving contractor. A number of modifications were made to the knife throughout the field tests. About 3000 feet of longitudinal joint was formed in 2001. Additional testing was done in the 2002 construction season, working with the same contractor. About 150,000 feet of longitudinal joint was formed in 2002. Evaluations of the formed joints were done to determine longitudinal joint hairline crack development rate and appearance. Additional tests will be done in the next construction season to improve or perfect the longitudinal joint forming technique.