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.

Active-control trials with binary data: a comparison of methods for testing superiority or non-inferiority using the odds ratio

This paper considers methods for testing for superiority or non-inferiority in active-control trials with binary data, when the relative treatment effect is expressed as an odds ratio. Three asymptotic tests for the log-odds ratio based on the unconditional binary likelihood are presented, namely the likelihood ratio, Wald and score tests. All three tests can be implemented straightforwardly in standard statistical software packages, as can the corresponding confidence intervals. Simulations indicate that the three alternatives are similar in terms of the Type I error, with values close to the nominal level. However, when the non-inferiority margin becomes large, the score test slightly exceeds the nominal level. In general, the highest power is obtained from the score test, although all three tests are similar and the observed differences in power are not of practical importance. Copyright (C) 2007 John Wiley & Sons, Ltd.

Detecting the node-density artifact in phylogeny reconstruction

The node-density effect is an artifact of phylogeny reconstruction that can cause branch lengths to be underestimated in areas of the tree with fewer taxa. Webster, Payne, and Pagel (2003, Science 301:478) introduced a statistical procedure (the "delta" test) to detect this artifact, and here we report the results of computer simulations that examine the test's performance. In a sample of 50,000 random data sets, we find that the delta test detects the artifact in 94.4% of cases in which it is present. When the artifact is not present (n = 10,000 simulated data sets) the test showed a type I error rate of approximately 1.69%, incorrectly reporting the artifact in 169 data sets. Three measures of tree shape or "balance" failed to predict the size of the node-density effect. This may reflect the relative homogeneity of our randomly generated topologies, but emphasizes that nearly any topology can suffer from the artifact, the effect not being confined only to highly unevenly sampled or otherwise imbalanced trees. The ability to screen phylogenies for the node-density artifact is important for phylogenetic inference and for researchers using phylogenetic trees to infer evolutionary processes, including their use in molecular clock dating. [Delta test; molecular clock; molecular evolution; node-density effect; phylogenetic reconstruction; speciation; simulation.]

Cognitive tests used in chronic adult human randomised controlled trial micronutrient and phytochemical intervention studies

In recent years there has been a rapid growth of interest in exploring the relationship between nutritional therapies and the maintenance of cognitive function in adulthood. Emerging evidence reveals an increasingly complex picture with respect to the benefits of various food constituents on learning, memory and psychomotor function in adults. However, to date, there has been little consensus in human studies on the range of cognitive domains to be tested or the particular tests to be employed. To illustrate the potential difficulties that this poses, we conducted a systematic review of existing human adult randomised controlled trial (RCT) studies that have investigated the effects of 24 d to 36 months of supplementation with flavonoids and micronutrients on cognitive performance. There were thirty-nine studies employing a total of 121 different cognitive tasks that met the criteria for inclusion. Results showed that less than half of these studies reported positive effects of treatment, with some important cognitive domains either under-represented or not explored at all. Although there was some evidence of sensitivity to nutritional supplementation in a number of domains (for example, executive function, spatial working memory), interpretation is currently difficult given the prevailing 'scattergun approach' for selecting cognitive tests. Specifically, the practice means that it is often difficult to distinguish between a boundary condition for a particular nutrient and a lack of task sensitivity. We argue that for significant future progress to be made, researchers need to pay much closer attention to existing human RCT and animal data, as well as to more basic issues surrounding task sensitivity, statistical power and type I error.

Flexible selection of a single treatment incorporating short-term endpoint information in a phase II/III clinical trial

Seamless phase II/III clinical trials in which an experimental treatment is selected at an interim analysis have been the focus of much recent research interest. Many of the methods proposed are based on the group sequential approach. This paper considers designs of this type in which the treatment selection can be based on short-term endpoint information for more patients than have primary endpoint data available. We show that in such a case, the familywise type I error rate may be inflated if previously proposed group sequential methods are used and the treatment selection rule is not specified in advance. A method is proposed to avoid this inflation by considering the treatment selection that maximises the conditional error given the data available at the interim analysis. A simulation study is reported that illustrates the type I error rate inflation and compares the power of the new approach with two other methods: a combination testing approach and a group sequential method that does not use the short-term endpoint data, both of which also strongly control the type I error rate. The new method is also illustrated through application to a study in Alzheimer's disease. © 2015 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.

Rational modification of the hierarchy of intermolecular interactions in molecular crystal structures by using tunable halogen bonds

A family of 16 isomolecular salts (3-XpyH)(2)[MX'(4)] (3-XpyH=3-halopyridinium; M=Co, Zn; X=(F), Cl, Br, (I); X'=Cl, Br, I) each containing rigid organic cations and tetrahedral halometallate anions has been prepared and characterized by X-ray single crystal and/or powder diffraction. Their crystal structures reflect the competition and cooperation between non-covalent interactions: N-H center dot center dot center dot X'-M hydrogen bonds, C-X center dot center dot center dot X'-M halogen bonds and pi-pi stacking. The latter are essentially unchanged in strength across the series, but both halogen bonds and hydrogen bonds are modified in strength upon changing the halogens involved. Changing the organic halogen (X) from F to I strengthens the C-X center dot center dot center dot X'-M halogen bonds, whereas an analogous change of the inorganic halogen (X') weakens both halogen bonds and N-H center dot center dot center dot X'-M hydrogen bonds. By so tuning the strength of the putative halogen bonds from repulsive to weak to moderately strong attractive interactions, the hierarchy of the interactions has been modified rationally leading to systematic changes in crystal packing. Three classes of crystal structure are obtained. In type A (C-F center dot center dot center dot X'-M) halogen bonds are absent. The structure is directed by N-H center dot center dot center dot X'-M hydrogen bonds and pi-stacking interactions. In type B structures, involving small organic halogens (X) and large inorganic halogens (X'), long (weak) C-X center dot center dot center dot X'-M interactions are observed with type I halogen-halogen interaction geometries (C-X center dot center dot center dot X' approximate to X center dot center dot center dot X'-M approximate to 155 degrees), but hydrogen bonds still dominate. Thus, minor but quite significant perturbations from the type A structure arise. In type C, involving larger organic halogens (X) and smaller inorganic halogens (X'), stronger halogen bonds are formed with a type II halogen-halogen interaction geometry (C-X center dot center dot center dot X' approximate to 180 degrees; X center dot center dot center dot X'-M approximate to 110 degrees) that is electrostatically attractive. The halogen bonds play a major role alongside hydrogen bonds in directing the type C structures, which as a result are quite different from type A and B.

Rotational spectra induced by vibrations

Mizushima and Venkateswarlu showed in 1953 that certain molecules have the property that excited vibrational states may possess rotational spectra even when the rotational spectrum of the ground vibrational state is forbidden by symmetry. We call such a spectrum a vibrationally induced rotational spectrum, and have made a systematic examination of the point groups which permit such behaviour. We also give formulae for the approximate line frequencies and intensities in these spectra, and discuss some of the problems involved in observing them. The spectra can only arise from degenerate vibrational states, and are of three possible types: i) symmetric top perpendicular spectra, shown by molecules belonging to the point groups Dnh, Dn and Cnh, where n is odd; (ii) symmetric top parallel spectra, shown by molecules belonging to Dnd and S2n, where n is even; and (iii) spherical top spectra, shown by molecules belonging to T or Td. Excited vibrational states of polar molecules of point groups Cnv or Cn, where n is odd, may also possess vibrationally induced perpendicular components of type (i), in addition to their ordinary parallel spectra. In addition to the above limitations on the point groups there are, in general, limitations on the symmetry species of the degenerate vibrational states.

Mutant alleles at the rugosus loci in pea affect seed moisture sorption isotherms and the relations between seed longevity and moisture content

Pea (Pisum sativum L.) mutant near-isogenic lines (RRrbrb, rrRbRb, rrrbrb) with lower starch but higher lipid contents, brought about by lesions in the starch biosynthetic pathway, had seed moisture sorption isotherms displaced below that of the wild type (RRRbRb). The negative logarithmic relationship between seed longevity and seed storage moisture content (%, f.wt basis), determined in hermetic storage at 65 degreesC, also differed: longevity in the mutant near-isogenic lines was poorer and less sensitive to moisture content than in the wild type (i.e. C-w was lower). The low-moisture-content limit (m(c)) to this relation also differed, being lower in the mutant near-isogenic lines (5.4-5.9%) than in the wild type (6.1%). In contrast, all four near-isogenic lines showed no difference (P >0.25) in the negative semilogarithmic relationship between equilibrium relative humidity (ERH) and seed longevity. It is concluded that the effect of these alleles at the r and rb loci on seed longevity. was largely indirect; a consequence of their effect on seed composition and hence on moisture sorption isotherms. However, this explanation could not be invoked at moisture contents below mc where differences in longevity remained substantial (RRRbRb double that of rrrbrb). Hence, these mutant alleles affected seed longevity directly at very low moisture contents.

Predation of Bradysia sp (Diptera: Sciaridae), Liriomyza trifolii (Diptera: Agromyzidae) and Bemisia tabaci (Hemiptera: Aleyrodidae) by Coenosia attenuata (Diptera: Muscidae) in greenhouse crops

We studied the predation behaviour of the "hunter fly" (Coenosia attenuata Stein) in the laboratory and greenhouse. In the laboratory, which was conducted at 25 degrees C at 60-80% RH, with a 16L : 8D photoperiod, we examined the functional response of this species to three different pests, namely the sciarid fly (Bradysia sp.), the tobacco whitefly (Bemisia tabaci) and the leaf miner Liriomyza trifolii. In the greenhouse, we studied the population dynamics of the predator and its prey on pepper and water melon crops grown in southern Spain. Adult hunter flies were found to exhibit a type I functional response to adult sciarid flies and whiteflies, but a type II response to adult leaf miners. The type II response was a result of the greater difficulty in capturing and handling leaf miners compared to the other two species. The dynamics of the predator-prey interaction in the greenhouse revealed that the predator specializes mainly on adult sciarids and that the presence of the other prey can be supplemental, but is never essential for survival of the predator; this, however, is crop-dependent. The results oil the dynamics of the predator-prey systems were obtained through a known population dynamics model with modifications.

Concentration of aluminium in breast cyst fluids collected from women affected by gross cystic breast disease

Gross cystic breast disease (GCBD) is the most common benign breast disorder, but the molecular basis of cyst formation remains to be identified. If the use of aluminium-based antiperspirant salts is involved in the etiology of gross breast cyst formation, it might be expected that aluminium would be at elevated levels in human breast cyst fluid (BCF). Aluminium was measured by ICP-MS in 48 samples of BCF, 30 samples of human blood serum and 45 samples of human breast milk at different stages of lactation (colostrum, intermediate, mature). The median level of aluminium in apocrine type I BCF (n:= 27, 150 mu g I-1) was significantly higher than in transuclative type II BCF (n = 21, 32 mu g I-1; P < 0.0001). By comparison, aluminium measurements gave a median concentration of 6 mu g I-1 in human serum and 25 mu g I-1 in human breast milk, with no difference between colostrum, intermediate and mature milk. Levels of aluminium were significantly higher in both types of BCF than in human serum (P < 0.0001). However when compared with human breast milk, aluminium levels were only significantly higher in apocrine type I BCF (P < 0.0001) and not in transudative type II BCF (P = 0.152). It remains to be identified why such high levels of aluminium were found in the apocrine type I BCF and from where the aluminium originated. However, if aluminium-based antiperspirants are found to be the source and to play any causal role in development of breast cysts, then it might become possible to prevent this common breast disorder. Copyright (C) 2008 John Wiley & Sons, Ltd.

NS1 proteins of avian influenza A viruses can act as antagonists of the human alpha/beta interferon response

Many viruses, including human influenza A virus, have developed strategies for counteracting the host type I interferon (IFN) response. We have explored whether avian influenza viruses were less capable of combating the type I IFN response in mammalian cells, as this might be a determinant of host range restriction. A panel of avian influenza viruses isolated between 1927 and 1997 was assembled. The selected viruses showed variation in their ability to activate the expression of a reporter gene under the control of the IFN-beta promoter and in the levels of IFN induced in mammalian cells. Surprisingly, the avian NS1 proteins expressed alone or in the genetic background of a human influenza virus controlled IFN-beta induction in a manner similar to the NS1 protein of human strains. There was no direct correlation between the IFN-beta induction and replication of avian influenza viruses in human A549 cells. Nevertheless, human cells deficient in the type I IFN system showed enhanced replication of the avian viruses studied, implying that the human type I IFN response limits avian influenza viruses and can contribute to host range restriction.

Bone morphogenetic proteins (BMP)-4,-6, and-7 potently suppress basal and luteinizing hormone-induced androgen production by bovine theca interna cells in primary culture: Could ovarian hyperandrogenic dysfunction be caused by a defect in thecal BMP signaling?

We reported recently that bovine theca interna cells in primary culture express several type-I and type-II receptors for bone morphogenetic proteins (BMPs). The same cells express at least two potential ligands for these receptors (BMP-4 and - 7), whereas bovine granulosa cells and oocytes express BMP-6. Therefore, BMPs of intrafollicular origin may exert autocrine/paracrine actions to modulate theca cell function. Here we report that BMP-4, - 6, and - 7 potently suppress both basal ( P < 0.0001; respective IC50 values, 0.78, 0.30, and 1.50 ng/ml) and LH-induced ( P < 0.0001; respective IC50 values, 5.00, 0.55, and 4.55 ng/ml) androgen production by bovine theca cells while having only a moderate effect on progesterone production and cell number. Semiquantitative RT-PCR showed that all three BMPs markedly reduced steady-state levels of mRNA for P450c17. Levels of mRNA encoding steroidogenic acute regulatory protein, P450scc, and 3 beta-hydroxysteroid dehydrogenase were also reduced but to a much lesser extent. Immunocytochemistry confirmed a marked reduction in cellular content of P450c17 protein after BMP treatment ( P < 0.001). Exposure to BMPs led to cellular accumulation of phosphorylated Smad1, but not Smad2, confirming that the receptors signal via a Smad1 pathway. The specificity of the BMP response was further explored by coincubating cells with BMPs and several potential BMP antagonists, chordin, gremlin, and follistatin. Gremlin and chordin were found to be effective antagonists of BMP-4 and - 7, respectively, and the observation that both antagonists enhanced ( P < 0.01) androgen production in the absence of exogenous BMP suggests an autocrine/paracrine role for theca-derived BMP- 4 and - 7 in modulating androgen production. Collectively, these data indicate that an intrafollicular BMP signaling pathway contributes to the negative regulation of thecal androgen production and that ovarian hyperandrogenic dysfunction could be a result of a defective autoregulatory pathway involving thecal BMP signaling.

Combined affinity labelling and mass spectrometry analysis of differential cell surface protein expression in normal and prostate cancer cells

Differences in the expression of cell surface proteins between a normal prostate epithelial (1542-NP2TX) and a prostate cancer cell line (1542-CP3TX) derived from the same patient were investigated. A combination of affinity chromatographic purification of biotin-tagged surface proteins with mass spectrometry analysis identified 26 integral membrane proteins and 14 peripheral surface proteins. The findings confirm earlier reports of altered expression in prostate cancer for several cell surface proteins, including ALCAM/CD166, the Ephrin type A receptor, EGFR and the prostaglandin F2 receptor regulatory protein. In addition, several novel findings of differential expression were made, including the voltage-dependent anion selective channel proteins Porin 1 and 2, ecto-5'-nucleotidase (CD73) and Scavenger receptor B1. Cell surface protein expression changed both qualitatively and quantitatively when the cells were grown in the presence of either or both interferon INFalpha and INFgamma. Costimulation with type I and II interferons had additive or synergistic effects on the membrane density of several, mainly peripherally attached surface proteins. Concerted upregulation of surface exposed antigens may be of benefit in immuno-adjuvant-based treatment of interferon-responsive prostate cancer. In conclusion, this study demonstrates that differences in the expression of membrane proteins between normal and prostate cancer cells are reproducibly detectable following vectorial labelling with biotin, and that detailed analysis of extracellular-induced surface changes can be achieved by combining surface-specific labelling with high-resolution two-dimensional gel electrophoresis and mass spectrometry.

Bone morphogenetic protein (BMP) ligands and receptors in bovine ovarian follicle cells: actions of BMP-4,-6 and-7 on granulosa cells and differential modulation of Smad-1 phosphorylation by follistatin

Given the paucity of information on the potential roles of bone morphogenetic proteins (BMPs) in the ruminant ovary we conducted immunolocalization and functional studies on cells isolated from bovine antral follicles. Immunocytochemistry revealed expression of BMP-4 and -7 in isolated theca cells whereas granulosa cells and oocytes selectively expressed RMP-6. All three cell types expressed a range of BMP-responsive type-I (BMPRIB, ActRI) and type-II (BMPRII, ActRII, ActRIIB) receptors supporting autocrine/paracrine roles within the follicle. This was reinforced by functional experiments on granulosa cells which showed that BMP-4, -6 and -7 promoted cellular accumulation of phosphorylated Smad-1 but not Smad-2 and enhanced 'basal' and IGF-stimulated secretion of oestradiol (E2), inhibin-A, activin-A and follistatin (FS). Concomitantly, each BMP suppressed 'basal' and IGF-stimulated progesterone secretion, consistent with an action to prevent or delay atresia and/or luteinization. BMPs also increased viable cell number under 'basal' (BMP-4 and -7) and IGF-stimulated (BMP-4, -6 and -7) conditions. Since FS, a product of bovine granulosa cells, has been shown to bind several BMPs, we used the Biacore technique to compare its binding affinities for activin-A (prototype FS ligand) and BMP-4, -6 and -7. Compared with activin-A (K-d 0.28 +/- 0.02 nM; 100%), the relative affinities of FS for BMP-4, -6 and -7 were 10, 5 and 1% respectively. Moreover, studies on granulosa cells showed that preincubation of ligand with excess FS abolished activin-A-induced phosphorylation of Smad-2 and BMP-4-induced phosphorylation of Smad-1. However, FS only partially reversed BMP-6-induced Smad-1 phosphorylation and had no inhibitory effect on BMP-7-induced Smad-1 phosphorylation. These findings support functional roles for BMP-4, -6 and -7 as paracrine/autocrine modulators of granulosa cell steroidogenesis, peptide secretion and proliferation in bovine antral follicles. The finding that FS can differentially modulate BMP-induced receptor activation and that this correlates with the relative binding affinity of FS for each BMP type implicates FS as a potential modulator of BMP action in the ovary.

Mutant alleles at the rugosus loci in pea affect seed moisture sorption isotherms and the relations between seed longevity and moisture content

Pea (Pisum sativum L.) mutant near-isogenic lines (RRrbrb, rrRbRb, rrrbrb) with lower starch but higher lipid contents, brought about by lesions in the starch biosynthetic pathway, had seed moisture sorption isotherms displaced below that of the wild type (RRRbRb). The negative logarithmic relationship between seed longevity and seed storage moisture content (%, f.wt basis), determined in hermetic storage at 65 °C, also differed: longevity in the mutant near-isogenic lines was poorer and less sensitive to moisture content than in the wild type (i.e. CW was lower). The low-moisture-content limit (mc) to this relation also differed, being lower in the mutant near-isogenic lines (5.4–5.9%) than in the wild type (6.1%). In contrast, all four near-isogenic lines showed no difference (P >0.25) in the negative semi-logarithmic relationship between equilibrium relative humidity (ERH) and seed longevity. It is concluded that the effect of these alleles at the r and rb loci on seed longevity was largely indirect; a consequence of their effect on seed composition and hence on moisture sorption isotherms. However, this explanation could not be invoked at moisture contents below mc where differences in longevity remained substantial (RRRbRb double that of rrrbrb). Hence, these mutant alleles affected seed longevity directly at very low moisture contents.