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.

Surface-functionalized ultrasmall superparamagnetic nanoparticles as magnetic delivery vectors for camptothecin.

Drug-nanoparticle conjugates: The anticancer drug camptothecin (CPT) was covalently linked at the surface of ultrasmall superparamagnetic iron oxide nanoparticles (USPIOs) via a linker, allowing drug release by cellular esterases. Nanoparticles were hierarchically built to achieve magnetically-enhanced drug delivery to human cancer cells and antiproliferative activity.The linking of therapeutic drugs to ultrasmall superparamagnetic iron oxide nanoparticles (USPIOs) allowing intracellular release of the active drug via cell-specific mechanisms would achieve tumor-selective magnetically-enhanced drug delivery. To validate this concept, we covalently attached the anticancer drug camptothecin (CPT) to biocompatible USPIOs (iron oxide core, 9-10 nm; hydrodynamic diameter, 52 nm) coated with polyvinylalcohol/polyvinylamine (PVA/aminoPVA). A bifunctional, end-differentiated dicarboxylic acid linker allowed the attachment of CPT to the aminoPVA as a biologically labile ester substrate for cellular esterases at one end, and as an amide at the other end. These CPT-USPIO conjugates exhibited antiproliferative activity in vitro against human melanoma cells. The intracellular localization of CPT-USPIOs was confirmed by transmission electron microscopy (iron oxide core), suggesting localization in lipid vesicles, and by fluorescence microscopy (CPT). An external static magnetic field applied during exposure increased melanoma cell uptake of the CPT-USPIOs.

Prediction of Childhood Asthma Using Conditional Probability and Discrete Event Simulation

Abstract: Asthma prevalence in children and adolescents in Spain is 10-17%. It is the most common chronic illness during childhood. Prevalence has been increasing over the last 40 years and there is considerable evidence that, among other factors, continued exposure to cigarette smoke results in asthma in children. No statistical or simulation model exist to forecast the evolution of childhood asthma in Europe. Such a model needs to incorporate the main risk factors that can be managed by medical authorities, such as tobacco (OR = 1.44), to establish how they affect the present generation of children. A simulation model using conditional probability and discrete event simulation for childhood asthma was developed and validated by simulating realistic scenario. The parameters used for the model (input data) were those found in the bibliography, especially those related to the incidence of smoking in Spain. We also used data from a panel of experts from the Hospital del Mar (Barcelona) related to actual evolution and asthma phenotypes. The results obtained from the simulation established a threshold of a 15-20% smoking population for a reduction in the prevalence of asthma. This is still far from the current level in Spain, where 24% of people smoke. We conclude that more effort must be made to combat smoking and other childhood asthma risk factors, in order to significantly reduce the number of cases. Once completed, this simulation methodology can realistically be used to forecast the evolution of childhood asthma as a function of variation in different risk factors.

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.

Graphical probabilistic analysis of the combination of items of evidence

Unlike the evaluation of single items of scientific evidence, the formal study and analysis of the jointevaluation of several distinct items of forensic evidence has to date received some punctual, ratherthan systematic, attention. Questions about the (i) relationships among a set of (usually unobservable)propositions and a set of (observable) items of scientific evidence, (ii) the joint probative valueof a collection of distinct items of evidence as well as (iii) the contribution of each individual itemwithin a given group of pieces of evidence still represent fundamental areas of research. To somedegree, this is remarkable since both, forensic science theory and practice, yet many daily inferencetasks, require the consideration of multiple items if not masses of evidence. A recurrent and particularcomplication that arises in such settings is that the application of probability theory, i.e. the referencemethod for reasoning under uncertainty, becomes increasingly demanding. The present paper takesthis as a starting point and discusses graphical probability models, i.e. Bayesian networks, as frameworkwithin which the joint evaluation of scientific evidence can be approached in some viable way.Based on a review of existing main contributions in this area, the article here aims at presentinginstances of real case studies from the author's institution in order to point out the usefulness andcapacities of Bayesian networks for the probabilistic assessment of the probative value of multipleand interrelated items of evidence. A main emphasis is placed on underlying general patterns of inference,their representation as well as their graphical probabilistic analysis. Attention is also drawnto inferential interactions, such as redundancy, synergy and directional change. These distinguish thejoint evaluation of evidence from assessments of isolated items of evidence. Together, these topicspresent aspects of interest to both, domain experts and recipients of expert information, because theyhave bearing on how multiple items of evidence are meaningfully and appropriately set into context.

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.