Sequentially testing for a gene-drug interaction in a genomewide analysis

Assaying a large number of genetic markers from patients in clinical trials is now possible in order to tailor drugs with respect to efficacy. The statistical methodology for analysing such massive data sets is challenging. The most popular type of statistical analysis is to use a univariate test for each genetic marker, once all the data from a clinical study have been collected. This paper presents a sequential method for conducting an omnibus test for detecting gene-drug interactions across the genome, thus allowing informed decisions at the earliest opportunity and overcoming the multiple testing problems from conducting many univariate tests. We first propose an omnibus test for a fixed sample size. This test is based on combining F-statistics that test for an interaction between treatment and the individual single nucleotide polymorphism (SNP). As SNPs tend to be correlated, we use permutations to calculate a global p-value. We extend our omnibus test to the sequential case. In order to control the type I error rate, we propose a sequential method that uses permutations to obtain the stopping boundaries. The results of a simulation study show that the sequential permutation method is more powerful than alternative sequential methods that control the type I error rate, such as the inverse-normal method. The proposed method is flexible as we do not need to assume a mode of inheritance and can also adjust for confounding factors. An application to real clinical data illustrates that the method is computationally feasible for a large number of SNPs. Copyright (c) 2007 John Wiley & Sons, Ltd.

Amplified fragment length polymorphism genetic fingerprinting challenges the taxonomic status of the near-endemic species Chara curta Nolte ex Kutz. (Characeae)

Amplified fragment length polymorphism (AFLP) genetic fingerprinting of 14 accessions of Chara curta and Chara aspera Willd., sampled across a range of habitats and morphologies in Britain, suggests that these taxa are part of the variation within a single species complex. Two primer combinations generating 397 fragments (97% of which were polymorphic), analysed by Jaccard's similarity coefficient and principal co-ordinate analysis, did not recover groups which reflect the current taxonomy. By contrast with the genetic study, a Gower general similarity coefficient and principal co-ordinate analysis of 52 morphological characters recovered the currently recognized species groups. A Mantel test showed no significant correlation between the genetic data and the morphological data, supporting the hypothesis that phenotypic variability in Chara L. is either to some extent environmentally induced or represents developmental stages. Implications for the conservation status of C. curta in Britain are discussed. (c) 2007 The Linnean Society of London, Botanical Journal of the Linnean Society, 2007, 155, 467-476.

The expected performance of single nucleotide polymorphism loci in paternity testing

We discuss the utility of single nucleotide polymorphism loci for full trio and mother-unavailable paternity testing cases, in the presence of population substructure and relatedness of putative and actual fathers. We focus primarily on the expected number of loci required to gain specified probabilities of mismatches, and report the expected proportion of paternity indices greater than three threshold values for these loci. (c) 2004 Elsevier Ireland Ltd. All rights reserved.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.

Single-trial event-related potential analysis for brain-computer interfaces

Abstract. Different types of mental activity are utilised as an input in Brain-Computer Interface (BCI) systems. One such activity type is based on Event-Related Potentials (ERPs). The characteristics of ERPs are not visible in single-trials, thus averaging over a number of trials is necessary before the signals become usable. An improvement in ERP-based BCI operation and system usability could be obtained if the use of single-trial ERP data was possible. The method of Independent Component Analysis (ICA) can be utilised to separate single-trial recordings of ERP data into components that correspond to ERP characteristics, background electroencephalogram (EEG) activity and other components with non- cerebral origin. Choice of specific components and their use to reconstruct “denoised” single-trial data could improve the signal quality, thus allowing the successful use of single-trial data without the need for averaging. This paper assesses single-trial ERP signals reconstructed using a selection of estimated components from the application of ICA on the raw ERP data. Signal improvement is measured using Contrast-To-Noise measures. It was found that such analysis improves the signal quality in all single-trials.

Algorithmic statistical analysis of electrophysiological data for the investigation of structure-activity relationship in single neurons

We are developing computational tools supporting the detailed analysis of the dependence of neural electrophysiological response on dendritic morphology. We approach this problem by combining simulations of faithful models of neurons (experimental real life morphological data with known models of channel kinetics) with algorithmic extraction of morphological and physiological parameters and statistical analysis. In this paper, we present the novel method for an automatic recognition of spike trains in voltage traces, which eliminates the need for human intervention. This enables classification of waveforms with consistent criteria across all the analyzed traces and so it amounts to reduction of the noise in the data. This method allows for an automatic extraction of relevant physiological parameters necessary for further statistical analysis. In order to illustrate the usefulness of this procedure to analyze voltage traces, we characterized the influence of the somatic current injection level on several electrophysiological parameters in a set of modeled neurons. This application suggests that such an algorithmic processing of physiological data extracts parameters in a suitable form for further investigation of structure-activity relationship in single neurons.

Analysis of disputed single-parent/child and sibling relationships using 16 STR loci

This study describes the validation of short tandem repeat (STR) systems for the resolution of cases of disputed parentage where only a single parent is available for testing or where the claimed relationship of both parents is in doubt and also cases where sibship must be tested. Three separate multiplex systems the Second Generation Multiplex, Powerplex 1.2 and FFFL have been employed, giving a total of 16 STR loci. Both empirical and theoretical approaches to the validation have been adopted. Appropriate equations have been derived to calculate likelihood ratios for different relationships, incorporating a correction for subpopulation effects. An F(ST) point estimate of 1% has been applied throughout. Empirically, 101 cases of alleged father, alleged mother and child where analysed using six SLP systems and also using the three multiplex STR systems. Of the 202 relationships tested, 197 were independently resolved by both systems, providing either clear evidence of non-parentage or strong support for the relationship.

NMR studies of the conformational effect of single and double 3'-S-phosphorothiolate substitutions within deoxythymidine trinucleotides

NMR spectroscopy has been used to investigate the conformational effects of single and two consecutive 3′-S-phosphorothiolate modifications within a deoxythymidine trinucleotide. The presence of a single 3′-phosphorothioate modification shifts the conformation of the sugar ring it is attached to, from a mainly south to north pucker; this effect is also transmitted to the 3′-neighbour deoxyribose. This transmission is thought to be caused by favourable stacking of the heterocyclic bases. Similar observations have been made previously by this group. When two adjacent modifications are present, the conformations of the attached deoxyribose rings are again shifted almost completely to the north, however, there is no transmission to the 3′ deoxyribose ring. Base proton chemical shift analysis and molecular modelling have been used to aid elucidation of the origin of this feature. The observation for the dimodified sequence is consistent with our previously reported results for a related system in which spaced modifications are more thermodynamically stable than consecutive ones.

Single-cell proteomic analysis of glucosinolate-rich S-cells in Arabidopsis thaliana

Single-cell analysis is essential for understanding the processes of cell differentiation and metabolic specialisation in rare cell types. The amount of single proteins in single cells can be as low as one copy per cell and is for most proteins in the attomole range or below; usually considered as insufficient for proteomic analysis. The development of modern mass spectrometers possessing increased sensitivity and mass accuracy in combination with nano-LC-MS/MS now enables the analysis of single-cell contents. In Arabidopsis thaliana, we have successfully identified nine unique proteins in a single-cell sample and 56 proteins from a pool of 15 single-cell samples from glucosinolate-rich S-cells by nanoLC-MS/MS proteomic analysis, thus establishing the proof-of-concept for true single-cell proteomic analysis. Dehydrin (ERD14_ARATH), two myrosinases (BGL37_ARATH and BGL38_ARATH), annexin (ANXD1_ARATH), vegetative storage proteins (VSP1_ARATH and VSP2_ARATH) and four proteins belonging to the S-adenosyl-l-methionine cycle (METE_ARATH, SAHH1_ARATH, METK4_ARATH and METK1/3_ARATH) with associated adenosine kinase (ADK1_ARATH), were amongst the proteins identified in these single-S-cell samples. Comparison of the functional groups of proteins identified in S-cells with epidermal/cortical cells and whole tissue provided a unique insight into the metabolism of S-cells. We conclude that S-cells are metabolically active and contain the machinery for de novo biosynthesis of methionine, a precursor for the most abundant glucosinolate glucoraphanine in these cells. Moreover, since abundant TGG2 and TGG1 peptides were consistently found in single-S-cell samples, previously shown to have high amounts of glucosinolates, we suggest that both myrosinases and glucosinolates can be localised in the same cells, but in separate subcellular compartments. The complex membrane structure of S-cells was reflected by the presence of a number of proteins involved in membrane maintenance and cellular organisation.

A detailed single chain model for molten poly(tetrafluoroethylene) from a novel structure refinement technique on neutron scattering data

We present a new methodology that couples neutron diffraction experiments over a wide Q range with single chain modelling in order to explore, in a quantitative manner, the intrachain organization of non-crystalline polymers. The technique is based on the assignment of parameters describing the chemical, geometric and conformational characteristics of the polymeric chain, and on the variation of these parameters to minimize the difference between the predicted and experimental diffraction patterns. The method is successfully applied to the study of molten poly(tetrafluoroethylene) at two different temperatures, and provides unambiguous information on the configuration of the chain and its degree of flexibility. From analysis of the experimental data a model is derived with CC and CF bond lengths of 1.58 and 1.36 Å, respectively, a backbone valence angle of 110° and a torsional angle distribution which is characterized by four isometric states, namely a split trans state at ± 18°, giving rise to a helical chain conformation, and two gauche states at ± 112°. The probability of trans conformers is 0.86 at T = 350°C, which decreases slightly to 0.84 at T = 400°C. Correspondingly, the chain segments are characterized by long all-trans sequences with random changes in sign, rather anisotropic in nature, which give rise to a rather stiff chain. We compare the results of this quantitative analysis of the experimental scattering data with the theoretical predictions of both force fields and molecular orbital conformation energy calculations.

Analysis of 16S rDNA sequences from pathogenic Leptospira serovars and use of single nucleotide polymorphisms for rapid speciation by D-HPLC

Leptospira have a worldwide distribution and include important zoonotic pathogens yet diagnosis and differentiation still tend to rely on traditional bacteriological and serological approaches. In this study a 1.3 kb fragment of the rrs gene (16S rDNA) was sequenced from a panel of 22 control strains, representing serovars within the pathogenic species Leptospira interrogans, Leptospira borgpetersenii, and Leptospira kirschneri, to identify single nucleotide polymorphisms (SNPs). These were identified in the 5' variable region of the 16S sequence and a 181 bp PCR fragment encompassing this region was used for speciation by Denaturing High Performance Liquid Chromatography (D-HPLC). This method was applied to eleven additional species, representing pathogenic, non-pathogenic and intermediate species and was demonstrated to rapidly differentiate all but 2 of the non-pathogenic Leptospira species. The method was applied successfully to infected tissues from field samples proving its value for diagnosing leptospiral infections found in animals in the UK. Crown Copyright (C) 2010 Published by Elsevier Ltd. All rights reserved.

The Glu298Asp single nucleotide polymorphism in the endothelial nitric oxide synthase gene differentially effects the vascular response to acute consumption of fruit and vegetable puree-based drinks

Scope Diets low in fruits and vegetables (FV) are responsible for 2.7 million deaths from cardiovascular diseases (CVD) and certain cancers annually. Many FV and their juices contain flavonoids, some of which increase endothelial nitric oxide synthase (eNOS) activity. A single nucleotide polymorphism in the eNOS gene, where thymine (T) replaces guanine (G) at position 894 predicting substitution of glutamate for aspartate at codon 298 (Glu298Asp), has been associated with increased CVD risk due to effects on nitric oxide synthesis and subsequently vascular reactivity. Individuals can be homozygous for guanine (GG), thymine (TT) or heterozygous (GT). Methods and results We investigated the effects of acute ingestion of a FV-puree-based-drink (FVPD) on vasodilation and antioxidant status in subjects retrospectively genotyped for this polymorphism. Healthy volunteers (n = 24; 11 GG, 11 GT, 2 TT) aged 30–70 were recruited to a randomized, controlled, crossover, acute study. We showed that acute consumption of 400 mL FVPD differentially affected individuals depending on their genotype. There was a significant genotype interaction for endothelium-dependent vasodilation measured by laser Doppler imaging with iontophoresis (P < 0.05) and ex vivo low-density lipoproteins (LDL) oxidation (P = 0.002). GG subjects had increased endothelium-dependent vasodilation 180 min (P = 0.028) and reduced ex vivo LDL oxidation (P = 0.013) after 60 min after FVPD compared with control, no differences were observed in GT subjects. Conclusion eNOS Glu298Asp genotype differentially affects vasodilation and ex vivo LDL oxidation after consumption of FV in the form of a puree-based drink.

Restriction fragment length polymorphism of polymerase chain reaction products applied to the differentiation of poultry campylobacters for epidemiological investigations

A technique for subtyping Camplobacter jejuni isolates has been developed by using the restriction fragment length polymorphism (Rnp) of polymerase chain reaction (PCR) products of the fluA and flaB genes. The technique was validated by using strains representing 28 serotypes of C jejuni and it may also be applied to C coli. From these strains 12 distinct RFLP profiles were observed but there was no direct relationship between the RFLP profile and the serotype. One hundred and thirty-five campylobacter isolates from 15 geographically distinct broiler flocks were investigated. All the isolates could be subtyped by using the RFLP method. Isolates from most of the flocks had a single RFLP profile despite data indicating that several serotypes were involved. Although it is possible that further restriction analysis may have demonstrated profile variations in these strains, it is more likely that antigenic variation can occur within genotypically related campylobacters. As a result, serotyping may give conflicting information for veterinary epidemiological purposes. This RFLP typing scheme appears to provide a suitable tool for the investigation of the sources and routes of transmission of campylobacters in chickens.

Analysis and modeling of ants’ behavior from single to multi-body

In order to shed light on the collective behavior of social insects, we analyzed the behavior of ants from single to multi-body. In an experimental set-up, ants are placed in hemisphere without a nest and food. Trajectory of ants is recorded. From this bottom-up approach, we found that collective behavior of ants as follows: 1. Activity of single ant increases and decreases periodically. 2. Spontaneous meeting process is observed between two ants and meeting spot of two ants is localized in hemisphere. 3. Result on division of labor is obtained between two ants.