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.

Learning about Bayesian networks for forensic interpretation : An example based on the "the problem of multiple propositions"

Both, Bayesian networks and probabilistic evaluation are gaining more and more widespread use within many professional branches, including forensic science. Notwithstanding, they constitute subtle topics with definitional details that require careful study. While many sophisticated developments of probabilistic approaches to evaluation of forensic findings may readily be found in published literature, there remains a gap with respect to writings that focus on foundational aspects and on how these may be acquired by interested scientists new to these topics. This paper takes this as a starting point to report on the learning about Bayesian networks for likelihood ratio based, probabilistic inference procedures in a class of master students in forensic science. The presentation uses an example that relies on a casework scenario drawn from published literature, involving a questioned signature. A complicating aspect of that case study - proposed to students in a teaching scenario - is due to the need of considering multiple competing propositions, which is an outset that may not readily be approached within a likelihood ratio based framework without drawing attention to some additional technical details. Using generic Bayesian networks fragments from existing literature on the topic, course participants were able to track the probabilistic underpinnings of the proposed scenario correctly both in terms of likelihood ratios and of posterior probabilities. In addition, further study of the example by students allowed them to derive an alternative Bayesian network structure with a computational output that is equivalent to existing probabilistic solutions. This practical experience underlines the potential of Bayesian networks to support and clarify foundational principles of probabilistic procedures for forensic evaluation.

A Robust Multiple Feature Approach To Endpoint Detection In Car Environment Based On Advanced Classifiers

In this paper we propose an endpoint detection system based on the use of several features extracted from each speech frame, followed by a robust classifier (i.e Adaboost and Bagging of decision trees, and a multilayer perceptron) and a finite state automata (FSA). We present results for four different classifiers. The FSA module consisted of a 4-state decision logic that filtered false alarms and false positives. We compare the use of four different classifiers in this task. The look ahead of the method that we propose was of 7 frames, which are the number of frames that maximized the accuracy of the system. The system was tested with real signals recorded inside a car, with signal to noise ratio that ranged from 6 dB to 30dB. Finally we present experimental results demonstrating that the system yields robust endpoint detection.

An allele-specific, stochastic gene expression process controls the expression of multiple Ly49 family genes and generates a diverse, MHC-specific NK cell receptor repertoire.

Mouse NK cells express MHC class I-specific inhibitory Ly49 receptors. Since these receptors display distinct ligand specificities and are clonally distributed, their expression generates a diverse NK cell receptor repertoire specific for MHC class I molecules. We have previously found that the Dd (or Dk)-specific Ly49A receptor is usually expressed from a single allele. However, a small fraction of short-term NK cell clones expressed both Ly49A alleles, suggesting that the two Ly49A alleles are independently and randomly expressed. Here we show that the genes for two additional Ly49 receptors (Ly49C and Ly49G2) are also expressed in a (predominantly) mono-allelic fashion. Since single NK cells can co-express multiple Ly49 receptors, we also investigated whether mono-allelic expression from within the tightly linked Ly49 gene cluster is coordinate or independent. Our clonal analysis suggests that the expression of alleles of distinct Ly49 genes is not coordinate. Thus Ly49 alleles are apparently independently and randomly chosen for stable expression, a process that directly restricts the number of Ly49 receptors expressed per single NK cell. We propose that the Ly49 receptor repertoire specific for MHC class I is generated by an allele-specific, stochastic gene expression process that acts on the entire Ly49 gene cluster.

Familial occurrence of an association of multiple intestinal atresia and choanal atresia: a new syndrome?

We report on two familial cases from a non-consanguineous marriage, presenting multiple intestinal and choanal atresia. Massive hydramnios and dilatation of the bowel were observed at 29 weeks of gestation during routine ultrasound scan of a healthy mother. The fetal karyotype was normal and cystic fibrosis screening was negative. Regular scans were performed throughout the pregnancy. The child was born at 34 weeks gestation. Choanal atresia was diagnosed at birth and abdominal investigations showed multiple atresia interesting both the small bowel and the colon. Further interventions were necessary because of recurrent obstructions. During the following pregnancy, a dilatation of the fetal intestinal tract was detected by ultrasonography at 27 weeks of gestation. Pregnancy was interrupted. Post-mortem examination of the fetus confirmed the stenosis of long segments of the small intestine associated with areas of colonic atresia. In both cases, histology and distribution were consistent with those reported in hereditary multiple intestinal atresia (HMIA). An association between multiple intestinal and choanal atresia has never been reported. We suggest it could correspond to a new autosomal recessive entity for which cytogenetic investigations and high-resolution array CGH revealed no visible anomalies.

Genetic abnormalities and survival in multiple myeloma: the experience of the Intergroupe Francophone du Myélome.

Acquired genomic aberrations have been shown to significantly impact survival in several hematologic malignancies. We analyzed the prognostic value of the most frequent chromosomal changes in a large series of patients with newly diagnosed symptomatic myeloma prospectively enrolled in homogeneous therapeutic trials. All the 1064 patients enrolled in the IFM99 trials conducted by the Intergroupe Francophone du Myélome benefited from an interphase fluorescence in situ hybridization analysis performed on purified bone marrow plasma cells. They were systematically screened for the following genomic aberrations: del(13), t(11;14), t(4;14), hyperdiploidy, MYC translocations, and del(17p). Chromosomal changes were observed in 90% of the patients. The del(13), t(11;14), t(4;14), hyperdiploidy, MYC translocations, and del(17p) were present in 48%, 21%, 14%, 39%, 13%, and 11% of the patients, respectively. After a median follow-up of 41 months, univariate statistical analyses revealed that del(13), t(4;14), nonhyperdiploidy, and del(17p) negatively impacted both the event-free survival and the overall survival, whereas t(11;14) and MYC translocations did not influence the prognosis. Multivariate analyses on 513 patients annotated for all the parameters showed that only t(4;14) and del(17p) retained prognostic value for both the event-free and overall survivals. When compared with the currently used International Staging System, this prognostic model compares favorably. In myeloma, the genomic aberrations t(4;14) and del(17p), together with beta2-microglobulin level, are important independent predictors of survival. These findings have implications for the design of risk-adapted treatment strategies.

Field Testing of Railroad Flatcar Bridges Volume II: Multiple Spans, TR-498, 2007

Based on the conclusions of IHRB Project TR-444, Demonstration Project Using Railroad Flat Car Bridges for Low Volume Road Bridges, additional research on the use of RRFC bridges was undertaken. This portion of the project investigated the following: (1) Different design and rating procedures; (2) Additional single span configurations plus multiple span configurations; (3) Different mechanisms for connecting adjacent RRFCs and the resulting lateral load distribution factors; (4) Sheet pile abutments; and (5) Behavior RRFCs that had been strengthened so that they could be used on existing abutments. A total of eight RRFC bridges were tested (five single span bridges, two two-span bridges, and one three-span bridge). Based on the results of this study a simplified design and rating procedure has been developed for the economical replacement bridge alternative. In Volume 1, the results from the testing of four single span RRFC bridges are presented, while in Volume 2,this volume, the results from the testing of the strengthened single span bridge plus the three multiple span bridges are presented.

Continuous monitoring and quantification of multiple parameters of daily physical activity in ambulatory Duchenne muscular dystrophy patients.

Multiple motor function and strength assessment tools exist for the evaluation of neuromuscular diseases, but most do not directly assess functional ability in the patients' daily physical activity in their home environment. In this study our aim was to assess: 1) the feasibility and accuracy of physical activity monitoring during two days in a home environment of five DMD patients using a non-commercialized monitor containing a 3D accelerometer and a gyroscope, 2) if a difference in the physical activity parameters could be measured before and one month after starting prednisolone. We reliably quantified the time spend sitting, standing, lying, walking, the number of steps taken, the cadence, the number of walking episodes and their duration as well as how these were distributed over the day. Parameters possibly reflecting endurance, such as the duration of the walking episodes or the succession of two or three walking episodes lasting more than 30 s were the most improved after prednisolone treatment. This degree of detailed determination of physical activity in a home environment has not been previously reported in neuromuscular disorders to our knowledge and some of the reported parameters are potential new outcome measures in clinical trials.

The M-Coffee web server: a meta-method for computing multiple sequence alignments by combining alternative alignment methods.

The M-Coffee server is a web server that makes it possible to compute multiple sequence alignments (MSAs) by running several MSA methods and combining their output into one single model. This allows the user to simultaneously run all his methods of choice without having to arbitrarily choose one of them. The MSA is delivered along with a local estimation of its consistency with the individual MSAs it was derived from. The computation of the consensus multiple alignment is carried out using a special mode of the T-Coffee package [Notredame, Higgins and Heringa (T-Coffee: a novel method for fast and accurate multiple sequence alignment. J. Mol. Biol. 2000; 302: 205-217); Wallace, O'Sullivan, Higgins and Notredame (M-Coffee: combining multiple sequence alignment methods with T-Coffee. Nucleic Acids Res. 2006; 34: 1692-1699)] Given a set of sequences (DNA or proteins) in FASTA format, M-Coffee delivers a multiple alignment in the most common formats. M-Coffee is a freeware open source package distributed under a GPL license and it is available either as a standalone package or as a web service from www.tcoffee.org.

Reduced left hemisphere dominance for multiple versus single anagram solutions and its modulation by individuals' schizotypy

Problem solving (including insight, divergent thinking) seems to rely on the right hemisphere (RH). These functions are difficult to assess behaviorally. We propose anagram resolution as a suitable paradigm. University students (n=32) performed three tachistoscopic lateralized visual half-field experiments (stimulus presentation 150ms). In Experiment 1, participants recalled four-letter strings. Subsequently, participants provided solutions for four-letter anagrams (one solution in Experiment 2; two solutions in Experiment 3). Additionally, participants completed a schizotypy questionnaire (O-LIFE). Results showed a right visual field advantage in Experiment 1 and 2, but no visual field advantage in Experiment 3. In Experiment 1, increasing positive schizotypy associated with a RH performance shift. Problem solving seems to require increasingly the RH when facing several rather than one solution. This result supports previous studies on the RH's role in remote associative, metaphor and discourse processing. The more complex language requirements, the less personality traits seem to matter.

Optogenetic Interrogation of Dopaminergic Modulation of the Multiple Phases of Reward-Seeking Behavior.

Phasic activation of dopaminergic neurons is associated with reward-predicting cues and supports learning during behavioral adaptation. While noncontingent activation of dopaminergic neurons in the ventral tegmental are (VTA) is sufficient for passive behavioral conditioning, it remains unknown whether the phasic dopaminergic signal is truly reinforcing. In this study, we first targeted the expression of channelrhodopsin-2 to dopaminergic neurons of the VTA and optimized optogenetically evoked dopamine transients. Second, we showed that phasic activation of dopaminergic neurons in freely moving mice causally enhances positive reinforcing actions in a food-seeking operant task. Interestingly, such effect was not found in the absence of food reward. We further found that phasic activation of dopaminergic neurons is sufficient to reactivate previously extinguished food-seeking behavior in the absence of external cues. This was also confirmed using a single-session reversal paradigm. Collectively, these data suggest that activation of dopaminergic neurons facilitates the development of positive reinforcement during reward-seeking and behavioral flexibility.

Targeting vascular NADPH oxidase 1 blocks tumor angiogenesis through a PPARα mediated mechanism.

Reactive oxygen species, ROS, are regulators of endothelial cell migration, proliferation and survival, events critically involved in angiogenesis. Different isoforms of ROS-generating NOX enzymes are expressed in the vasculature and provide distinct signaling cues through differential localization and activation. We show that mice deficient in NOX1, but not NOX2 or NOX4, have impaired angiogenesis. NOX1 expression and activity is increased in primary mouse and human endothelial cells upon angiogenic stimulation. NOX1 silencing decreases endothelial cell migration and tube-like structure formation, through the inhibition of PPARα, a regulator of NF-κB. Administration of a novel NOX-specific inhibitor reduced angiogenesis and tumor growth in vivo in a PPARα dependent manner. In conclusion, vascular NOX1 is a critical mediator of angiogenesis and an attractive target for anti-angiogenic therapies.

Sustained NKG2D engagement induces cross-tolerance of multiple distinct NK cell activation pathways

NKG2D is a multisubunit activation receptor that allows natural killer (NK) cells to detect and eliminate stressed, infected, and transformed host cells. However, the chronic exposure of NK cells to cell-bound NKG2D ligands has been shown to impair NKG2D function both in vitro and in vivo. Here we have tested whether continuous NKG2D engagement selectively impacted NKG2D function or whether heterologous NK cell activation pathways were also affected. We found that sustained NKG2D engagement induced cross-tolerization of several unrelated NK cell activation receptors. We show that receptors that activate NK cells via the DAP12/KARAP and DAP10 signaling adaptors, such as murine NKG2D and Ly49D, cross-tolerize preferentially NK cell activation pathways that function independent of DAP10/12, such as antibody-dependent cell-mediated cytotoxicity and missing-self recognition. Conversely, DAP10/12-independent pathways are unable to cross-tolerize unrelated NK cell activation receptors such as NKG2D or Ly49D. These data define a class of NK cell activation receptors that can tolerize mature NK cells. The reversible suppression of the NK cells' cytolytic function probably reduces the NK cells' efficacy to control endogenous and exogenous stress yet may be needed to limit tissue damage

Multiple processes regulate long-term population dynamics of sea urchins on Mediterranean rocky reefs

We annually monitored the abundance and size structure of herbivorous sea urchin populations (Paracentrotus lividus and Arbacia lixula) inside and outside a marine reserve in the Northwestern Mediterranean on two distinct habitats (boulders and vertical walls) over a period of 20 years, with the aim of analyzing changes at different temporal scales in relation to biotic and abiotic drivers. P. lividus exhibited significant variability in density over time on boulder bottoms but not on vertical walls, and temporal trends were not significantly different between the protection levels. Differences in densities were caused primarily by variance in recruitment, which was less pronounced inside the MPA and was correlated with adult density, indicating density-dependent recruitment under high predation pressure, as well as some positive feedback mechanisms that may facilitate higher urchin abundances despite higher predator abundance. Populations within the reserve were less variable in abundance and did not exhibit the hyper-abundances observed outside the reserve, suggesting that predation effects maybe more subtle than simply lowering the numbers of urchins in reserves. A. lixula densities were an order of magnitude lower than P. lividus densities and varied within sites and over time on boulder bottoms but did not differ between protection levels. In December 2008, an exceptionally violent storm reduced sea urchin densities drastically (by 50% to 80%) on boulder substrates, resulting in the lowest values observed over the entire study period, which remained at that level for at least two years (up to the present). Our results also showed great variability in the biological and physical processes acting at different temporal scales. This study highlights the need for appropriate temporal scales for studies to fully understand ecosystem functioning, the concepts of which are fundamental to successful conservation and management.