Detecting epistasis with restricted response patterns in pairs of biallelic loci.

Well-established examples of genetic epistasis between a pair of loci typically show characteristic patterns of phenotypic distributions in joint genotype tables. However, inferring epistasis given such data is difficult due to the lack of power in commonly used approaches, which decompose the epistatic patterns into main plus interaction effects followed by testing the interaction term. Testing additive-only or all terms may have more power, but they are sensitive to nonepistatic patterns. Alternatively, the epistatic patterns of interest can be enumerated and the best matching one is found by searching through the possibilities. Although this approach requires multiple testing correction over possible patterns, each pattern can be fitted with a regression model with just one degree of freedom and thus the overall power can still be high, if the number of possible patterns is limited. Here we compare the power of the linear decomposition and pattern search methods, by applying them to simulated data generated under several patterns of joint genotype effects with simple biological interpretations. Interaction-only tests are the least powerful; while pattern search approach is the most powerful if the range of possibilities is restricted, but still includes the true pattern.

Regional abundance of on-premise outlets and drinking patterns among Swiss young men: District level analyses and geographic adjustments.

INTRODUCTION AND AIMS: This study investigated the associations of alcohol outlet density with specific alcohol outcomes (consumption and consequences) among young men in Switzerland and assessed the possible geographically related variations. DESIGN AND METHODS: Alcohol consumption and drinking consequences were measured in a 2010-2011 study assessing substance use risk factors (Cohort Study on Substance Use Risk Factors) among 5519 young Swiss men. Outlet density was based on the number of on- and off-premise outlets in the district of residence. Linear regression models were run separately for drinking level, heavy episodic drinking (HED) and drinking consequences. Geographically weighted regression models were estimated when variations were recorded at the district level. RESULTS: No consistent association was found between outlet density and drinking consequences. A positive association between drinking level and HED with on-premise outlet density was found. Geographically weighted regressions were run for drinking level and HED. The predicted values for HED were higher in the southwest part of Switzerland (French-speaking part). DISCUSSION AND CONCLUSIONS: Among Swiss young men, the density of outlets and, in particular, the abundance of bars, clubs and other on-premise outlets was associated with drinking level and HED, even when drinking consequences were not significantly affected. These findings support the idea that outlet density needs to be considered when developing and implementing regional-based prevention initiatives. [Astudillo M, Kuendig H, Centeno-Gil A, Wicki M, Gmel G. Regional abundance of on-premise outlets and drinking patterns among Swiss young men: District level analyses and geographic adjustments. Drug Alcohol Rev 2014;33:526-33].

Training effects with the observer-rated cognitive errors and the coping action patterns scales

The use of observer-rated scales requires that raters be trained until they have become reliable in using the scales. However, few studies properly report how training in using a given rating scale is conducted or indeed how it should be conducted. This study examined progress in interrater reliability over 6 months of training with two observer-rated scales, the Cognitive Errors Rating Scale and the Coping Action Patterns Rating Scale. The evolution of the intraclass correlation coefficients was modeled using hierarchical linear modeling. Results showed an overall training effect as well as effects of the basic training phase and of the rater calibration phase, the latter being smaller than the former. The results are discussed in terms of implications for rater training in psychotherapy research.

Genetic structure of a linear population of Beta vulgaris ssp maritima (sea beet) revealed by isozyme and RFLP analysis

Knowledge of the genetic structure of plant populations is necessary for the understanding of the dynamics of major ecological processes. It also has applications in conservation biology and risk assessment for genetically modified crops. This paper reports the genetic structure of a linear population of sea beet, Beta vulgaris ssp. maritima (the wild relative of sugar beet), on Furzey Island, Poole Harbour. The relative spatial positions of the plants were accurately mapped and the plants were scored for variation at isozyme and RFLP loci. Structure was analysed by repeated subdivision of the population to find the average size of a randomly mating group. Estimates of F-ST between randomly mating units were then made, and gave patterns consistent with the structure of the population being determined largely by founder effects. The implications of these results for the monitoring of transgene spread in wild sea beet populations are discussed.

Alliance patterns over the course of short-term dynamic psychotherapy: The shape of productive relationships.

The shape of alliance processes over the course of psychotherapy has already been studied in several process-outcome studies on very brief psychotherapy. The present study applies the shape-of-change methodology to short-term dynamic psychotherapies and complements this method with hierarchical linear modeling. A total of 50 psychotherapies of up to 40 sessions were included. Alliance was measured at the end of each session. The results indicate that a linear progression model is most adequate. Three main patterns were found: stable, linear, and quadratic growth. The linear growth pattern, along with the slope parameter, was related to treatment outcome. This study sheds additional light on alliance process research, underscores the importance of linear alliance progression for outcome, and also fosters a better understanding of its limitations.

Simple estimation of incident HIV infection rates in notification cohorts based on window periods of algorithms for evaluation of line-immunoassay result patterns.

BACKGROUND: Tests for recent infections (TRIs) are important for HIV surveillance. We have shown that a patient's antibody pattern in a confirmatory line immunoassay (Inno-Lia) also yields information on time since infection. We have published algorithms which, with a certain sensitivity and specificity, distinguish between incident (< = 12 months) and older infection. In order to use these algorithms like other TRIs, i.e., based on their windows, we now determined their window periods. METHODS: We classified Inno-Lia results of 527 treatment-naïve patients with HIV-1 infection < = 12 months according to incidence by 25 algorithms. The time after which all infections were ruled older, i.e. the algorithm's window, was determined by linear regression of the proportion ruled incident in dependence of time since infection. Window-based incident infection rates (IIR) were determined utilizing the relationship 'Prevalence = Incidence x Duration' in four annual cohorts of HIV-1 notifications. Results were compared to performance-based IIR also derived from Inno-Lia results, but utilizing the relationship 'incident = true incident + false incident' and also to the IIR derived from the BED incidence assay. RESULTS: Window periods varied between 45.8 and 130.1 days and correlated well with the algorithms' diagnostic sensitivity (R(2) = 0.962; P<0.0001). Among the 25 algorithms, the mean window-based IIR among the 748 notifications of 2005/06 was 0.457 compared to 0.453 obtained for performance-based IIR with a model not correcting for selection bias. Evaluation of BED results using a window of 153 days yielded an IIR of 0.669. Window-based IIR and performance-based IIR increased by 22.4% and respectively 30.6% in 2008, while 2009 and 2010 showed a return to baseline for both methods. CONCLUSIONS: IIR estimations by window- and performance-based evaluations of Inno-Lia algorithm results were similar and can be used together to assess IIR changes between annual HIV notification cohorts.

3D seismic facies characterization and geological patterns recognition (Australian North West Shelf)

EXECUTIVE SUMMARY This PhD research, funded by the Swiss Sciences Foundation, is principally devoted to enhance the recognition, the visualisation and the characterization of geobodies through innovative 3D seismic approaches. A series of case studies from the Australian North West Shelf ensures the development of reproducible integrated 3D workflows and gives new insight into local and regional stratigraphic as well as structural issues. This project was initiated in year 2000 at the Geology and Palaeontology Institute of the University of Lausanne (Switzerland). Several collaborations ensured the improvement of technical approaches as well as the assessment of geological models. - Investigations into the Timor Sea structural style were carried out at the Tectonics Special Research Centre of the University of Western Australia and in collaboration with Woodside Energy in Perth. - Seismic analysis and attributes classification approach were initiated with Schlumberger Oilfield Australia in Perth; assessments and enhancements of the integrated seismic approaches benefited from collaborations with scientists from Schlumberger Stavanger Research (Norway). Adapting and refining from "linear" exploration techniques, a conceptual "helical" 3D seismic approach has been developed. In order to investigate specific geological issues this approach, integrating seismic attributes and visualisation tools, has been refined and adjusted leading to the development of two specific workflows: - A stratigraphic workflow focused on the recognition of geobodies and the characterization of depositional systems. Additionally, it can support the modelling of the subsidence and incidentally the constraint of the hydrocarbon maturity of a given area. - A structural workflow used to quickly and accurately define major and secondary fault systems. The integration of the 3D structural interpretation results ensures the analysis of the fault networks kinematics which can affect hydrocarbon trapping mechanisms. The application of these integrated workflows brings new insight into two complex settings on the Australian North West Shelf and ensures the definition of astonishing stratigraphic and structural outcomes. The stratigraphic workflow ensures the 3D characterization of the Late Palaeozoic glacial depositional system on the Mermaid Nose (Dampier Subbasin, Northern Carnarvon Basin) that presents similarities with the glacial facies along the Neotethys margin up to Oman (chapter 3.1). A subsidence model reveals the Phanerozoic geodynamic evolution of this area (chapter 3.2) and emphasizes two distinct mode of regional extension for the Palaeozoic (Neotethys opening) and Mesozoic (abyssal plains opening). The structural workflow is used for the definition of the structural evolution of the Laminaria High area (Bonaparte Basin). Following a regional structural characterization of the Timor Sea (chapter 4.1), a thorough analysis of the Mesozoic fault architecture reveals a local rotation of the stress field and the development of reverse structures (flower structures) in extensional setting, that form potential hydrocarbon traps (chapter 4.2). The definition of the complex Neogene structural architecture associated with the fault kinematic analysis and a plate flexure model (chapter 4.3) suggest that the Miocene to Pleistocene reactivation phases recorded at the Laminaria High most probably result from the oblique normal reactivation of the underlying Mesozoic fault planes. This episode is associated with the deformation of the subducting Australian plate. Based on these results three papers were published in international journals and two additional publications will be submitted. Additionally this research led to several communications in international conferences. Although the different workflows presented in this research have been primarily developed and used for the analysis of specific stratigraphic and structural geobodies on the Australian North West Shelf, similar integrated 3D seismic approaches will have applications to hydrocarbon exploration and production phases; for instance increasing the recognition of potential source rocks, secondary migration pathways, additional traps or reservoir breaching mechanisms. The new elements brought by this research further highlight that 3D seismic data contains a tremendous amount of hidden geological information waiting to be revealed and that will undoubtedly bring new insight into depositional systems, structural evolution and geohistory of the areas reputed being explored and constrained and other yet to be constrained. The further development of 3D texture attributes highlighting specific features of the seismic signal, the integration of quantitative analysis for stratigraphic and structural processes, the automation of the interpretation workflow as well as the formal definition of "seismo-morphologic" characteristics of a wide range of geobodies from various environments would represent challenging examples of continuation of this present research. The 21st century will most probably represent a transition period between fossil and other alternative energies. The next generation of seismic interpreters prospecting for hydrocarbon will undoubtedly face new challenges mostly due to the shortage of obvious and easy targets. They will probably have to keep on integrating techniques and geological processes in order to further capitalise the seismic data for new potentials definition. Imagination and creativity will most certainly be among the most important quality required from such geoscientists.

Environment and dispersal paths override life strategies and residence time in determining regional patterns of invasion by alien plants.

We describe a novel dissimilarity framework to analyze spatial patterns of species diversity and illustrate it with alien plant invasions in Northern Portugal. We used this framework to test the hypothesis that patterns of alien invasive plant species richness and composition are differently affected by differences in climate, land use and landscape connectivity (i.e. Geographic distance as a proxy and vectorial objects that facilitate dispersal such as roads and rivers) between pairs of localities at the regional scale. We further evaluated possible effects of plant life strategies (Grime's C-S-R) and residence time. Each locality consisted of a 1 km(2) landscape mosaic in which all alien invasive species were recorded by visiting all habitat types. Multi-model inference revealed that dissimilarity in species richness is more influenced by environmental distance (particularly climate), whereas geographic distance (proxies for dispersal limitations) is more important to explain dissimilarity in species composition, with a prevailing role for ecotones and roads. However, only minor differences were found in the responses of the three C-S-R strategies. Some effect of residence time was found, but only for dissimilarity in species richness. Our results also indicated that environmental conditions (e.g. climate conditions) limit the number of alien species invading a given site, but that the presence of dispersal corridors determines the paths of invasion and therefore the pool of species reaching each site. As geographic distances (e.g. ecotones and roads) tend to explain invasion at our regional scale highlights the need to consider the management of alien invasions in the context of integrated landscape planning. Alien species management should include (but not be limited to) the mitigation of dispersal pathways along linear infrastructures. Our results therefore highlight potentially useful applications of the novel multimodel framework to the anticipation and management of plant invasions. (C) 2013 Elsevier GmbH. All rights reserved.

Resting-State Functional Connectivity Emerges from Structurally and Dynamically Shaped Slow Linear Fluctuations.

Brain fluctuations at rest are not random but are structured in spatial patterns of correlated activity across different brain areas. The question of how resting-state functional connectivity (FC) emerges from the brain's anatomical connections has motivated several experimental and computational studies to understand structure-function relationships. However, the mechanistic origin of resting state is obscured by large-scale models' complexity, and a close structure-function relation is still an open problem. Thus, a realistic but simple enough description of relevant brain dynamics is needed. Here, we derived a dynamic mean field model that consistently summarizes the realistic dynamics of a detailed spiking and conductance-based synaptic large-scale network, in which connectivity is constrained by diffusion imaging data from human subjects. The dynamic mean field approximates the ensemble dynamics, whose temporal evolution is dominated by the longest time scale of the system. With this reduction, we demonstrated that FC emerges as structured linear fluctuations around a stable low firing activity state close to destabilization. Moreover, the model can be further and crucially simplified into a set of motion equations for statistical moments, providing a direct analytical link between anatomical structure, neural network dynamics, and FC. Our study suggests that FC arises from noise propagation and dynamical slowing down of fluctuations in an anatomically constrained dynamical system. Altogether, the reduction from spiking models to statistical moments presented here provides a new framework to explicitly understand the building up of FC through neuronal dynamics underpinned by anatomical connections and to drive hypotheses in task-evoked studies and for clinical applications.

New insight in the relationship between regional patterns of knee cartilage thickness, osteoarthritis disease severity, and gait mechanics.

To test if the relationship between knee kinetics during walking and regional patterns of cartilage thickness is influenced by disease severity we tested the following hypotheses in a cross-sectional study of medial compartment osteoarthritis (OA) subjects: (1) the peak knee flexion (KFM) and adduction moments (KAM) during walking are associated with regional cartilage thickness and medial-to-lateral cartilage thickness ratios, and (2) the associations between knee moments and cartilage thickness data are dependent on disease severity. Seventy individuals with medial compartment knee OA were studied. Gait analysis was used to determine the knee moments and cartilage thickness was measured from magnetic resonance imaging. Multiple linear regression analyses tested for associations between cartilage thickness and knee kinetics. Medial cartilage thickness and medial-to-lateral cartilage thickness ratios were lower in subjects with greater KAM for specific regions of the femoral condyle and tibial plateau with no associations for KFM in patients of all disease severities. When separated by severity, the association between KAM and cartilage thickness was found only in patients with more severe OA, and KFM was significantly associated with cartilage thickness only for the less severe OA subjects for specific tibial plateau regions. The results support the idea that the KAM is larger in patients with more severe disease and the KFM has greater influence early in the disease process, which may lessen as pain increases with disease severity. Each component influences different regions of cartilage. Thus the relative contributions of both KAM and KFM should be considered when evaluating gait mechanics and the influence of any intervention for knee OA.

Patterns of oxygen consumption during establishment of cephalocaudal polarity in the early chick embryo.

Oxygen uptake was studied during the establishment of cephalocaudal polarity in the very early chick embryo, i.e., 10 hr before (stage VI) and at laying (stage X). Oxygen fluxes in minute regions of the intact blastoderms were measured in vitro by scanning microspectrophotometry in the presence or absence of glucose. The oxygen consumption of the whole blastoderm remained constant (6 nmol O2 X hr-1) throughout the period studied, although the number of cells increased more than twofold. The regional oxygen fluxes varied from 0.41 to 1.13 nmol O2 X hr-1 X mm-2 at stage VI and from 0.42 to 0.70 nmol O2 X hr-1 X mm-2 at stage X. At stage VI, the oxygen flux in the center of the blastoderm was significantly higher than that in its periphery. This pattern remained evident when the values were corrected for cell number or for cytoplasmic volume. At stage X, there was a tendency for the oxygen fluxes to decrease from the posterior to the anterior regions of the area pellucida. Thus the pattern of oxidative metabolism in the late uterine embryos seems to change from radial to bilateral. This change of symmetry probably reflects the process of formation of the embryonic axis. In addition, the fact that the oxygen uptake was similar in the presence or absence of glucose suggests that early chick embryos metabolize essentially intracellular stores.

Predicting trabecular bone microdamage initiation and accumulation using a non-linear perfect damage model

Studies evaluating the mechanical behavior of the trabecular microstructure play an important role in our understanding of pathologies such as osteoporosis, and in increasing our understanding of bone fracture and bone adaptation. Understanding of such behavior in bone is important for predicting and providing early treatment of fractures. The objective of this study is to present a numerical model for studying the initiation and accumulation of trabecular bone microdamage in both the pre- and post-yield regions. A sub-region of human vertebral trabecular bone was analyzed using a uniformly loaded anatomically accurate microstructural three-dimensional finite element model. The evolution of trabecular bone microdamage was governed using a non-linear, modulus reduction, perfect damage approach derived from a generalized plasticity stress-strain law. The model introduced in this paper establishes a history of microdamage evolution in both the pre- and post-yield regions

Coronary heart disease and cerebrovascular disease mortality in young adults: recent trends in Europe.

Background: Over the last two decades, mortality from coronary heart disease (CHD) and cerebrovascular disease (CVD) declined by about 30% in the European Union (EU). Design: We analyzed trends in CHD (X ICD codes: I20-I25) and CVD (X ICD codes: I60-I69) mortality in young adults (age 35-44 years) in the EU as a whole and in 12 selected European countries, over the period 1980-2007. Methods: Data were derived from the World Health Organization mortality database. With joinpoint regression analysis, we identified significant changes in trends and estimated average annual percent changes (AAPC). Results: CHD mortality rates at ages 35-44 years have decreased in both sexes since the 1980s for most countries, except for Russia (130/100,000 men and 24/100,000 women, in 2005-7). The lowest rates (around 9/100,000 men, 2/100,000 women) were in France, Italy and Sweden. In men, the steepest declines in mortality were in the Czech Republic (AAPC = -6.1%), the Netherlands (-5.2%), Poland (-4.5%), and England and Wales (-4.5%). Patterns were similar in women, though with appreciably lower rates. The AAPC in the EU was -3.3% for men (rate = 16.6/100,000 in 2005-7) and -2.1% for women (rate = 3.5/100,000). For CVD, Russian rates in 2005-7 were 40/100,000 men and 16/100,000 women, 5 to 10-fold higher than in most western European countries. The steepest declines were in the Czech Republic and Italy for men, in Sweden and the Czech Republic for women. The AAPC in the EU was -2.5% in both sexes, with steeper declines after the mid-late 1990s (rates = 6.4/100,000 men and 4.3/100,000 women in 2005-7). Conclusions: CHD and CVD mortality steadily declined in Europe, except in Russia, whose rates were 10 to 15-fold higher than those of France, Italy or Sweden. Hungary and Poland, and also Scotland, where CHD trends were less favourable than in other western European countries, also emerge as priorities for preventive interventions.