Estimating the proportion of true null hypotheses and adaptive false discovery rate control in discrete paradigm

Storey's estimator for the proportion of true null hypotheses, originally proposed under the continuous framework, has been modified in this work under the discrete framework. The modification results in improved estimation of the parameter of interest. The proposed estimator is used to formulate an adaptive version of the Benjamini–Hochberg procedure. Control over the false discovery rate by the proposed adaptive procedure has been proved analytically. The proposed estimate is also used to formulate an adaptive version of the Benjamini–Hochberg–Heyse procedure. Simulation experiments establish the conservative nature of this new adaptive procedure. Substantial amount of gain in power is observed for the new adaptive procedures over the standard procedures. For demonstration of the proposed method, two important real life gene expression data sets, one related to the study of HIV and the other related to methylation study, are used.     

Sample size calculation for multiple testing in microarray data analysis

Microarray technology is rapidly emerging for genome-wide screening of differentially expressed genes between clinical subtypes or different conditions of human diseases. Traditional statistical testing approaches, such as the two-sample t-test or Wilcoxon test, are frequently used for evaluating statistical significance of informative expressions but require adjustment for large-scale multiplicity. Due to its simplicity, Bonferroni adjustment has been widely used to circumvent this problem. It is well known, however, that the standard Bonferroni test is often very conservative. In the present paper, we compare three multiple testing procedures in the microarray context: the original Bonferroni method, a Bonferroni-type improved single-step method and a step-down method. The latter two methods are based on nonparametric resampling, by which the null distribution can be derived with the dependency structure among gene expressions preserved and the family-wise error rate accurately controlled at the desired level. We also present a sample size calculation method for designing microarray studies. Through simulations and data analyses, we find that the proposed methods for testing and sample size calculation are computationally fast and control error and power precisely.     

On sample size and inference for two-stage adaptive designs

Proschan and Hunsberger (1995, Biometrics 51, 1315-1324) proposed a two-stage adaptive design that maintains the Type I error rate. For practical applications, a two-stage adaptive design is also required to achieve a desired statistical power while limiting the maximum overall sample size. In our proposal, a two-stage adaptive design is comprised of a main stage and an extension stage, where the main stage has sufficient power to reject the null under the anticipated effect size and the extension stage allows increasing the sample size in case the true effect size is smaller than anticipated. For statistical inference, methods for obtaining the overall adjusted p-value, point estimate and confidence intervals are developed. An exact two-stage test procedure is also outlined for robust inference.     

The role of conditioning film formation in Pseudomonas aeruginosa PAO1 adhesion to inert surfaces in aquatic environments

Bacterial initial adhesion to inert surfaces in aquatic environments is highly dependent on the surface properties of the substratum, which can be altered significantly by the formation of conditioning films. In this study, the impact of conditioning films formed with extracellular polymeric substances (EPS) on bacterial adhesion was investigated. Adhesion of wild type Pseudomonas aeruginosa PAO1 to slides coated with model EPS components (alginate, humic substances, and bovine serum albumin (BSA)) was examined. Surface roughness of conditioning film coated slides was evaluated by atomic force microscopy (AFM), and its effect on the bacterial initial adhesion was not significant. X-ray photoelectron spectroscopy (XPS) studies were performed to determine the elemental surface compositions of bacterial cells and substrates. Results showed that bacterial adhesion to bare slides and slides coated with alginate and humic substances increased as ionic strength increased. Conversely, BSA coating enhanced bacterial adhesion at low ionic strength but hindered adhesion at higher ionic strength. It was concluded that forces other than hydrophobic and electrostatic interactions were involved in controlling bacterial adhesion to BSA coated surfaces. A steric model for polymer brushes that considers the combined influence of steric effects and DLVO interaction forces was shown to adequately describe the observed bacterial adhesion behaviors.     

Modeling heat stress effects on dairy cattle milk production in a tropical environment using test-day records and random regression models

In tropical environments, dairy cattle production is constrained by several factors, including climate. The seasonal loss of milk due to heat stress is a recurring challenge for many dairy producers. The objective of this study was to detect heat stress thresholds, milk yield loss and individual animal variations using random regression models for dairy cattle from test-day milk records. Data were obtained from the Kenya Livestock Breeders Organization for the years 2000–2017 and merged with weather data. The weather parameters were grid-interpolated solar and meteorological data obtained from the National Aeronautics and Space Administration/Prediction Of Worldwide Energy Resources (NASA/POWER). After editing, the records comprised 49 993, 45 251 and 36 136 test-day records for first, second, and third lactations, respectively, for the four main dairy breeds: Friesian (68.0%), Ayrshire (21.1%), Jersey (7.6%) and Guernsey (3.3%). Variance components were estimated using Restricted Maximum Likelihood in ASReml software. Random regression models with third-order Legendre polynomials were fitted to the average and individual lactation curves and the reaction norms. An extended factor analytic variance structure for the random cow effects was used to estimate (co)variances between days in milk and thermal load. The daily average temperature (TA) and temperature humidity index (THI) were identified as the most suitable thermal load indicators for assessing milk yield losses. Considering a one day lag, the estimated heat stress thresholds were about 22 °C and 69 index units for TA and THI, respectively. Almost no differences were observed for estimated residual variances between the thermal load indicators, indicating there was no better model fit by TA or THI. The heat stress thresholds and milk loss patterns are important for management of dairy production systems in the tropics with climatic conditions similar to this study. Data recording should be improved as a tool to monitor the expected impacts of climate change and mitigation measures.     

The configuration of the Cu2+ binding region in full-length human prion protein

The cellular prion protein (PrP^C) is a Cu²⁺ binding protein connected to the outer cell membrane. The molecular features of the Cu²⁺ binding sites have been investigated and characterized by spectroscopic experiments on PrP^C-derived peptides and the recombinant human full-length PrP^C (hPrP-[23-231]). The hPrP-[23-231] was loaded with ⁶³Cu under slightly acidic (pH 6.0) or neutral conditions. The PrP^C/Cu²⁺-complexes were investigated by extended X-ray absorption fine structure (EXAFS), electron paramagnetic resonance (EPR), and electron nuclear double resonance (ENDOR). For comparison, peptides from the copper-binding octarepeat domain were investigated in different environments. Molecular mechanics computations were used to select sterically possible peptide/Cu²⁺ structures. The simulated EPR, ENDOR, and EXAFS spectra of these structures were compared with our experimental data. For a stoichiometry of two octarepeats per copper the resulting model has a square planar four nitrogen Cu²⁺ coordination. Two nitrogens belong to imidazole rings of histidine residues. Further ligands are two deprotonated backbone amide nitrogens of the adjacent glycine residues and an axial oxygen of a water molecule. Our complex model differs significantly from those previously obtained for shorter peptides. Sequence context, buffer conditions and stoichiometry of copper show marked influence on the configuration of copper binding to PrP^C. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Indirect effects on fitness between individuals that have never met via an extended phenotype

Interactions between organisms are ubiquitous and have important consequences for phenotypes and fitness. Individuals can even influence those they never meet, if they have extended phenotypes that alter the environments others experience. North American red squirrels (Tamiasciurus hudsonicus) guard food hoards, an extended phenotype that typically outlives the individual and is usually subsequently acquired by non‐relatives. Hoarding by previous owners can, therefore, influence subsequent owners. We found that red squirrels breed earlier and had higher lifetime fitness if the previous hoard owner was a male. This was driven by hoarding behaviour, as males and mid‐aged squirrels had the largest hoards, and these effects persisted across owners, such that if the previous owner was male or died in mid‐age, subsequent occupants had larger hoards. Individuals can, therefore, influence each other's resource‐dependent traits and fitness without ever meeting, such that the past can influence contemporary population dynamics through extended phenotypes.     

On organism: environment buffers and their ecological significance

We consider, from a physical perspective, the case where the interface between an organism and its environment becomes large enough that it acts as a buffer regulating their matter and energy exchanges. We illustrate the physiological and evolutionary role of buffers through the example of lungfish estivation. Then we ponder the relevance of buffers of this kind to the quest for a general definition of concepts like niche construction, the extended phenotype, and related ones, whose meaning is conveyed at present mostly through particular examples. Finally, we comment on the potential significance of buffers to organism—environment codetermination in the sense originally suggested by Lewontin.