CDS: a fold-change based statistical test for concomitant identification of distinctness and similarity in gene expression analysis

The problem of identifying differential activity such as in gene expression is a major defeat in biostatistics and bioinformatics.Equally important,however much less frequently studied,is the question of similar activity from one biological condition to another.The foldchange,or ratio,is usually considered a relevant criterion for stating difference and similarity between measurements.Importantly,no statistical method for concomitant evaluation of similarity and distinctness currently exists for biological applications.Modem microarray,digital PCR(dPCR),and Next-Generation Sequencing(NGS) technologies frequently provide a means of coefficient of variation estimation for individual measurements.Using fold-change,and by making the assumption that measurements are normally distributed with known variances,we designed a novel statistical test that allows us to detect concomitantly,thus using the same formalism,differentially and similarly expressed genes(http://cds.ihes.fr).Given two sets of gene measurements in different biological conditions,the probabilities of making type I and type II errors in stating that a gene is differentially or similarly expressed from one condition to the other can be calculated.Furthermore,a confidence interval for the fold-change can be delineated.Finally,we demonstrate that the assumption of normality can be relaxed to consider arbitrary distributions numerically.The Concomitant evaluation of Distinctness and Similarity(CDS) statistical test correctly estimates similarities and differences between measurements of gene expression.The implementation,being time and memory efficient,allows the use of the CDS test in high-throughput data analysis such as microarray,dPCR,and NGS experiments.Importantly,the CDS test can be applied to the comparison of single measurements(N:1) provided the variance(or coefficient of variation) of the signals is known,making CDS a valuable tool also in biomedical analysis where typically a single measurement per subject is available.     

Evaluating alternative stem cell hypotheses for adult corneal epithelial maintenance

In this review we evaluate evidence for three different hypotheses that explain how the corneal epithelium is maintained. The limbal epithelial stem cell(LESC)hypothesis is most widely accepted. This proposes that stem cells in the basal layer of the limbal epithelium, at the periphery of the cornea, maintain themselves and also produce transient(or transit) amplifying cells(TACs). TACs then move centripetally to the centre of the cornea in the basal layer of the corneal epithelium and also replenish cells in the overlying suprabasal layers. The LESCs maintain the corneal epithelium during normal homeostasis and become more active to repair significant wounds. Second, the corneal epithelial stem cell(CESC) hypothesis postulates that, during normal homeostasis, stem cells distributed throughout the basal corneal epithelium, maintain the tissue. According to this hypothesis, LESCs are present in the limbus but are only active during wound healing. We also consider a third possibility, that the corneal epithelium is maintained during normal homeostasis by proliferation of basal corneal epithelial cells without any input from stem cells. After reviewing the published evidence, we conclude that the LESC and CESC hypotheses are consistent with more of the evidence than the third hypothesis, so we do not consider this further. The LESC and CESC hypotheses each have difficulty accounting for one main type of evidence so we evaluate the two key lines of evidence that discriminate between them. Finally, we discuss how lineage-tracing experiments have begun to resolve the debate in favour of the LESC hypothesis. Nevertheless, it also seems likely that some basal corneal epithelial cells can act as long-term progenitors if limbal stem cell function is compromised. Thus, this aspect of the CESC hypothesis may have a lasting impact on our understanding of corneal epithelial maintenance, even if it is eventually shown that stem cells are restricted to the limbus as proposed by the LESC hypothesis.     

Repetitive Sequences in Plant Nuclear DNA:Types,Distribution, Evolution and Function

Repetitive DNA sequences are a major component of eukaryotic genomes and may account for up to 90% of the genome size. They can be divided into minisatellite, microsatellite and satellite sequences. Satellite DNA sequences are considered to be a fast-evolving component of eukaryotic genomes, comprising tandemly-arrayed, highly-repetitive and highly-conserved monomer sequences. The monomer unit of satellite DNA is 150–400 base pairs(bp) in length.Repetitive sequences may be species- or genus-specific, and may be centromeric or subtelomeric in nature. They exhibit cohesive and concerted evolution caused by molecular drive, leading to high sequence homogeneity. Repetitive sequences accumulate variations in sequence and copy number during evolution, hence they are important tools for taxonomic and phylogenetic studies, and are known as ‘‘tuning knobs' ' in the evolution. Therefore, knowledge of repetitive sequences assists our understanding of the organization, evolution and behavior of eukaryotic genomes. Repetitive sequences have cytoplasmic, cellular and developmental effects and play a role in chromosomal recombination. In the post-genomics era, with the introduction of next-generation sequencing technology, it is possible to evaluate complex genomes for analyzing repetitive sequences and deciphering the yet unknown functional potential of repetitive sequences.     

What is the sound of fear？ Behavioral responses of white-crowned sparrows Zonotrichia leucophrys to synthesized nonlinear acoustic phenomena

Fear and anxiety may be adaptive responses to life-threatening situations, and animals may communicate fear to others vocally. A fundamental understanding of fear inducing sounds is important for both wildlife conservation and management because it helps us understand how to design repellents and also how （and why） animals may be negatively impacted by anthropogenic sounds. Nonlinear phenomena--sounds produced by the desynchronization of vibrations in a sound production system-are commonly found in stress-induced animal vocalizations, such as in alarm calls, mobbing calls, and fear screams. There are several functional hypotheses for these nonlinear phenomena. One specific hypothesis is the unpredictability hypothesis, which suggests that because nonlinear phenomena are more variable and somewhat unpredictable, animals are less likely to habituate to them. Animals should, therefore, have a prolonged response to sounds with nonlinear phenomena than sounds without them. Most of the studies involving nonlinear phenomena have used mammalian subjects and conspecific stimuli. Our study fo- cused on white-crowned sparrows （Zonotrichia leucophrys ssp. oriantha） and used synthesized acoustic stimuli to investigate behavioral responses to stimuli with and without nonlinear phenomena. We predicted that birds would be less relaxed after hearing a stimulus with a nonlinear component. We calculated the difference from baseline of proportion of time spent in relaxed behaviors and performed pair-wise comparisons between a pure tone control stimulus and each of three experimental stimuli, including a frequency jump up, a frequency jump down, and white noise. These comparisons showed that in the 30q50 s after the playback experiment, birds were significantly less relaxed after hearing noise or an abrupt frequency jump down an octave but not an abrupt frequency jump up an octave or a pure tone. Nonlinear phenomena, therefore, may be generally arousing to animals and may explain why these acoustic properties are commonly fou     

Ecological Niche Divergence between Trapelus ruderatus(Olivier,1807) and T.persicus(Blanford,1881)(Sauria:Agamidae) in the Middle East

Modeling the potential distribution areas for a given species is important in understanding the relationship between the actual distribution and the most suitable habitat for a species. In this study, we obtained all available records of Trapelus ruderatus and Trapelus persicus from museums, literature and fieldwork and used them with environmental layers in the Maximum Entropy algorithm to predict highly suitable habitat areas. The distribution model of T. ruderatus and T. persicus showed excellent performance for both models(T.ruderatus AUC =0.964 ± 0.001 and T.persicus AUC = 0.996 ± 0.003), and predicted suitable regions in Iran, Turkey, Iraq and Syria. Niche overlap was measured between the two groups by ENMtools and 13% overlapped. We used a niche identity test to determine differences between the niches of the two species. Finally, by comparing our null hypothesis to the true niche overlap of the two species, we were able to reject our null hypothesis of no difference between the niches. Due to the sympatric distribution pattern of these species, we do not need a background test for niche divergence.     

Electronic Nose with an Air Sensor Matrix for Detecting Beef Freshness

China is one of the largest meat producing countries in the wodd. With the growing concern for food safety more attention has been paid to meat quality. The application of conventional test methods for meat quality is limited by many factors, and subjectiveness, such as longer time to prepare samples and to test. A sensor matrix was constructed with several separate air sensors, and tests were conducted to detect the freshness of the beef. The results show that the air sensors TGS2610, TGS2600, TGS2611, TGS2620 and TGS2602 made by Tianjin Figaro Electronic Co, Ltd could be used to determine the degree of freshness but TGS2442 is not suitable. This study provides a foundation for designing and making an economical and practical detector for beef freshness.     

Construction of brewing-wine Aspergillus oryzaepyrG- mutant bypyrG gene deletion and its application in homology transformation

pyrG- host cells are indispensable for pyrG- based transformation system. Isolations ofpyrG- host cells by random mutations are limited by time-consuming, unclear genetic background and potential interferences of homogenous recombination. The purpose of this study was to construct brewing-wine Aspergillus oryzae pyrG- mutant by sitedirected mutation of pyrG gene deletion which would be used as a host for further transformation, pMD-pyrGAB, a vector carrying pyrG deletion cassette, was used to con- struct pyrG- mutant of A. oryzae. Three stable pyrG deletion mutants of A. oryzae were isolated by resistant to 5-fluoroorotic acid and confirmed by polymerase chain reaction analysis, indicating thatpyrG was completely excised. The ApyrG mutants were applied as pyrG- host cells to disrupt xdh gene encoding xylitol dehydrogenase, which involves in xylitol production ofA. Oryzae. The xdh disrup- tion mutants were efficiently constructed by transforming a pMD-pyrC-xdh disruption plasmid carrying pyrG, and the produced xylitol concentration of the Axdh mutant was three times as much as that of the ApyrG recipient. Site-directed pyrG gene deletion is thus an effective way for the isolation of pyrG- host cells, and the established host-vector system could be applied in further functional genomics analysis and molecular breeding ofA. Oryzae.     

Applying Methods from Differential Geometry to Devise Stable and Persistent Air Layers Attached to Objects Immersed in Water

We describe a few mathematical tools which allow to investigate whether air-water interfaces exist(under prescribedconditions)and are mechanically stable and temporally persistent.In terms of physics,air-water interfaces are governed by theYoung-Laplace equation.Mathematically they are surfaces of constant mean curvature which represent solutions of a nonlinearelliptic partial differential equation.Although explicit solutions of this equation can be obtained only in very special cases,it is-under moderately special circumstances-possible to establish the existence of a solution without actually solving thedifferential equation.We also derive criteria for mechanical stability and temporal persistence of an air layer.Furthermorewe calculate the lifetime of a non-persistent air layer.Finally,we apply these tools to two examples which exhibit the symmetriesof 2D lattices.These examples can be viewed as abstractions of the biological model represented by the aquatic fernSalvinia.     

Quantitative genetic analysis of chlorophyll a fluorescence parameters in maize in the field environments

Chlorophyll fluorescence transient from initial to maximum fluorescence("P" step) throughout two intermediate steps("J" and "I")(JIP‐test) is considered a reliable early quantitative indicator of stress in plants. The JIP‐test is particularly useful for crop plants when applied in variable field environments. The aim of the present study was to conduct a quantitative trait loci(QTL) analysis for nine JIP‐test parameters in maize during flowering in four field environments differing in weather conditions. QTL analysis and identification of putative candidate genes might help to explain the genetic relationship between photosynthesis and different field scenarios in maize plants. The JIP‐test parameters were analyzed in the intermated B73 Mo17(IBM) maize population of 205 recombinant inbred lines. A set of 2,178 molecular markers across the whole maize genome was used for QTL analysis revealing 10 significant QTLs for seven JIP‐test parameters, of which five were co‐localized when combinedover the four environments indicating polygenic inheritance and pleiotropy. Our results demonstrate that QTL analysis of chlorophyll fluorescence parameters was capable of detecting one pleiotropic locus on chromosome 7, coinciding with the gene gst23 that may be associated with efficient photosynthesis under different field scenarios.     

Bioconvection in Casson nanofluid flow with Gyrotactic microorganisms and variable surface heat flux

This paper presents a two-dimensional unsteady laminar boundary layer mixed convection flow heat and mass transfer along a vertical plate filled with Casson nanofluid located in a porous quiescent medium that contains both nanoparticles and gyrotactic microorganisms. This permeable vertical plate is assumed to be moving in the same direction as the free stream velocity. The flow is subject to a variable heat flux, a zero nanoparticle flux and a constant density of motile microorganisms on the surface. The free stream velocity is time-dependent resulting in a non-similar solution. The transport equations are solved using the bivariate spectral quasilinearization method. A grid independence test for the validity of the result is given. The significance of the inclusion of motile microorganisms to heat transfer processes is discussed. We show, inter alia, that introducing motile microorganisms into the flow reduces the skin friction coefficient and that the random motion of the nanoparticles improves the rate of transfer of the motile microorganisms.     

A note on the Wilcoxon-Mann-Whitney test for 2 X kappa ordered tables

Biological and medical investigations often use ordered categorical data. When two groups are to be compared and the data for the groups fall in three or more ordered categories, the Wilcoxon-Mann-Whitney (WMW) test uses information in the ordering to give a test that is usually powerful against shift alternatives. However, such applications of WMW often involve distributions for which extensive ties play an important role. Newly available computer programs for performing exact tests give deeper insights into the characteristics of the exact WMW distributions and the suitability of normal approximations. We offer practical advice, based on experience with published biomedical data sets and on numerical studies of hypothetical ordered tables, for the use of WMW and its normal approximations.     

Power and sample size estimation for the Wilcoxon rank sum test with application to comparisons of C statistics from alternative prediction models

Summary .  The Wilcoxon Mann-Whitney (WMW) U test is commonly used in nonparametric two-group comparisons when the normality of the underlying distribution is questionable. There has been some previous work on estimating power based on this procedure ( Lehmann, 1998 , Nonparametrics ). In this article, we present an approach for estimating type II error, which is applicable to any continuous distribution, and also extend the approach to handle grouped continuous data allowing for ties. We apply these results to obtaining standard errors of the area under the receiver operating characteristic curve (AUROC) for risk-prediction rules under H ₁ and for comparing AUROC between competing risk prediction rules applied to the same data set. These results are based on SAS -callable functions to evaluate the bivariate normal integral and are thus easily implemented with standard software.     

Statistical power when testing for genetic differentiation

A variety of statistical procedures are commonly employed when testing for genetic differentiation. In a typical situation two or more samples of individuals have been genotyped at several gene loci by molecular or biochemical means, and in a first step a statistical test for allele frequency homogeneity is performed at each locus separately, using, e.g. the contingency chi-square test, Fisher's exact test, or some modification thereof. In a second step the results from the separate tests are combined for evaluation of the joint null hypothesis that there is no allele frequency difference at any locus, corresponding to the important case where the samples would be regarded as drawn from the same statistical and, hence, biological population. Presently, there are two conceptually different strategies in use for testing the joint null hypothesis of no difference at any locus. One approach is based on the summation of chi-square statistics over loci. Another method is employed by investigators applying the Bonferroni technique (adjusting the P-value required for rejection to account for the elevated alpha errors when performing multiple tests simultaneously) to test if the heterogeneity observed at any particular locus can be regarded significant when considered separately. Under this approach the joint null hypothesis is rejected if one or more of the component single locus tests is considered significant under the Bonferroni criterion. We used computer simulations to evaluate the statistical power and realized alpha errors of these strategies when evaluating the joint hypothesis after scoring multiple loci. We find that the 'extended' Bonferroni approach generally is associated with low statistical power and should not be applied in the current setting. Further, and contrary to what might be expected, we find that 'exact' tests typically behave poorly when combined in existing procedures for joint hypothesis testing. Thus, while exact tests are generally to be preferred over approximate ones when testing each particular locus, approximate tests such as the traditional chi-square seem preferable when addressing the joint hypothesis.     

An Exact Paired Rank Test

An exact rank test for two dependent samples based on overall mid‐ranks is discussed which can be applied to metric as well as to ordinal data. The exact conditional distribution of the test statistic given the observed vector of rank differences is determined. A recursion formula is given as well as a fast shift algorithm in SAS/IML code. Moreover, it is demonstrated that the paired rank test can be more powerful than other tests for paired samples by means of a simulation study. Finally, the test is applied to a psychiatric trial with longitudinal ordinal data.     

Signed-Rank Tests in the Presence of Ties (with extended tables)

We begin with a review of the areas of application of the signed-rank tests (SRTs) and we conclude that the results are exact only if no ties of non-null differences exist. In order to apply the SRTs according to WILCOXON and according to PRATT also in the presence of ties, by assigning midranks, we derive their null distributions. As special cases the null distributions for the problem without ties are obtained. In order to save the practising statistician the time-consuming calculations of the distribution functions, we compute tables of critical values (for reasons of volume they will be published as part of the reprints only). For N₀ = 0 (1) 5 null differences and M = = 1(1) 10 non-null differences the critical values of all distributions with all possible tie vectors are calculated. Instructions are provided and an example serves to illustrate the use of the table. The extension of the tables are obtained by means of counting formulas given in the text. Approximations are provided in order to make the application of tests possible for larger samples as well. It is shown that the approximation of the null distribution in the presence of ties by the null distributions under the assumption of no ties in some cases overstates and sometimes understates the exact rejection probability. For N₀ = 0 (1) 10 and M = 1 (1) 10 all distributions with all possible tie vectors for the SRTs with WILCOXON and PRATT ranking are examined with respect to the lattice type of the test statistic. The result is given in table 6. It is evident that the portion of PRATT -distributions with lattice character decreases as the number of null differences increases. Continuity corrections are obtained for the asymptotic normal distribution which take into account the lattice character of the distribution of the test statistic.     

The Exact Distribution of the Anderson-Kannemann Statistic

The Anderson-Kannemann test is a rank test for treatment effects in a randomized block design with K treatments and N blocks. In this paper, an algorithm for computing the exact distribution of the Anderson-Kannemann test statistic under the null hypothesis is deduced. Then, the exact distribution is compared with the asymptotic χ²-distribution, and it is shown that the exact distribution is approximated fairly well by the asymptotic distribution. Tables of the exact distribution are given for K = 3, N = 3(1)15; K = 4, N = 3(1)11; K = 5, N = 3(1)7; and K = 6, N = 3(1)5.     

A Comparison of the Three Conditional Exact Tests in Two‐way Contingency Tables Using the Unconditional Exact Power

The conditional exact tests of homogeneity of two binomial proportions are often used in small samples, because the exact tests guarantee to keep the size under the nominal level. The Fisher's exact test, the exact chi‐squared test and the exact likelihood ratio test are popular and can be implemented in software StatXact. In this paper we investigate which test is the best in small samples in terms of the unconditional exact power. In equal sample cases it is proved that the three tests produce the same unconditional exact power. A symmetry of the unconditional exact power is also found. In unequal sample cases the unconditional exact powers of the three tests are computed and compared. In most cases the Fisher's exact test turns out to be best, but we characterize some cases in which the exact likelihood ratio test has the highest unconditional exact power. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Efficient exact p-value computation for small sample, sparse, and surprising categorical data.

A major obstacle in applying various hypothesis testing procedures to datasets in bioinformatics is the computation of ensuing p-values. In this paper, we define a generic branch-and-bound approach to efficient exact p-value computation and enumerate the required conditions for successful application. Explicit procedures are developed for the entire Cressie-Read family of statistics, which includes the widely used Pearson and likelihood ratio statistics in a one-way frequency table goodness-of-fit test. This new formulation constitutes a first practical exact improvement over the exhaustive enumeration performed by existing statistical software. The general techniques we develop to exploit the convexity of many statistics are also shown to carry over to contingency table tests, suggesting that they are readily extendible to other tests and test statistics of interest. Our empirical results demonstrate a speed-up of orders of magnitude over the exhaustive computation, significantly extending the practical range for performing exact tests. We also show that the relative speed-up gain increases as the null hypothesis becomes sparser, that computation precision increases with increase in speed-up, and that computation time is very moderately affected by the magnitude of the computed p-value. These qualities make our algorithm especially appealing in the regimes of small samples, sparse null distributions, and rare events, compared to the alternative asymptotic approximations and Monte Carlo samplers. We discuss several established bioinformatics applications, where small sample size, small expected counts in one or more categories (sparseness), and very small p-values do occur. Our computational framework could be applied in these, and similar cases, to improve performance.     

Unconditional Exact Tests for Equivalence or Noninferiority for Paired Binary Endpoints

Problems of establishing equivalence or noninferiority between two medical diagnostic procedures involve comparisons of the response rates between correlated proportions. When the sample size is small, the asymptotic tests may not be reliable. This article proposes an unconditional exact test procedure to assess equivalence or noninferiority. Two statistics, a sample-based test statistic and a restricted maximum likelihood estimation (RMLE)-based test statistic, to define the rejection region of the exact test are considered. We show the p-value of the proposed unconditional exact tests can be attained at the boundary point of the null hypothesis. Assessment of equivalence is often based on a comparison of the confidence limits with the equivalence limits. We also derive the unconditional exact confidence intervals on the difference of the two proportion means for the two test statistics. A typical data set of comparing two diagnostic procedures is analyzed using the proposed unconditional exact and asymptotic methods. The p-value from the unconditional exact tests is generally larger than the p-value from the asymptotic tests. In other words, an exact confidence interval is generally wider than the confidence interval obtained from an asymptotic test.     

A cautionary note on exact unconditional inference for a difference between two independent binomial proportions

Fisher's exact test for comparing response proportions in a randomized experiment can be overly conservative when the group sizes are small or when the response proportions are close to zero or one. This is primarily because the null distribution of the test statistic becomes too discrete, a partial consequence of the inference being conditional on the total number of responders. Accordingly, exact unconditional procedures have gained in popularity, on the premise that power will increase because the null distribution of the test statistic will presumably be less discrete. However, we caution researchers that a poor choice of test statistic for exact unconditional inference can actually result in a substantially less powerful analysis than Fisher's conditional test. To illustrate, we study a real example and provide exact test size and power results for several competing tests, for both balanced and unbalanced designs. Our results reveal that Fisher's test generally outperforms exact unconditional tests based on using as the test statistic either the observed difference in proportions, or the observed difference divided by its estimated standard error under the alternative hypothesis, the latter for unbalanced designs only. On the other hand, the exact unconditional test based on the observed difference divided by its estimated standard error under the null hypothesis (score statistic) outperforms Fisher's test, and is recommended. Boschloo's test, in which the p-value from Fisher's test is used as the test statistic in an exact unconditional test, is uniformly more powerful than Fisher's test, and is also recommended.