TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

TiO2-Photocatalyzed degradation of phenol in saline media in an annular reactor: hydrodynamics, lumped kinetics, intermediates, and acute toxicity

The photocatalytic degradation of phenol in aqueous suspensions of TiO2 under different salt concentrations in an annular reactor has been investigated. In all cases, complete removal of phenol and mineralization degrees above 90% were achieved. The reactor operational parameters were optimized and its hydrodynamics characterized in order to couple mass balance equations with kinetic ones. The photodegradation of the organics followed a Langmuir-Hinshelwood-Hougen-Watson lumped kinetics. From GC/MS analyses, several intermediates formed during oxidation have been identified. The main ones were catechol, hydroquinone, and 3-phenyl-2-propenal, in this order. The formation of negligible concentrations of 4-chlorophenol was observed only in high salinity medium. Acute toxicity was determined by using Artemia sp. as the test organism, which indicated that intermediate products were all less toxic than phenol and a significant abatement of the overall toxicity was accomplished, regardless of the salt concentration.

Observation of the suppression of the flux of cosmic rays above 4x10(19) eV

The energy spectrum of cosmic rays above 2.5 x 10(18) eV, derived from 20 000 events recorded at the Pierre Auger Observatory, is described. The spectral index gamma of the particle flux, J proportional to E(-gamma), at energies between 4 x 10(18) eV and 4 x 10(19) eV is 2.69 +/- 0.02(stat) +/- 0.06(syst), steepening to 4.2 +/- 0.4(stat) +/- 0: 06 (syst) at higher energies. The hypothesis of a single power law is rejected with a significance greater than 6 standard deviations. The data are consistent with the prediction by Greisen and by Zatsepin and Kuz'min.

Upper limit on the diffuse flux of ultrahigh energy tau neutrinos from the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory is sensitive to Earth-skimming tau neutrinos that interact in Earth's crust. Tau leptons from nu(tau) charged-current interactions can emerge and decay in the atmosphere to produce a nearly horizontal shower with a significant electromagnetic component. The data collected between 1 January 2004 and 31 August 2007 are used to place an upper limit on the diffuse flux of nu(tau) at EeV energies. Assuming an E(nu)(-2) differential energy spectrum the limit set at 90% C. L. is E(nu)(2)dN(nu tau)/dE(nu) < 1: 3 x 10(-7) GeV cm(-2) s(-1) sr(-1) in the energy range 2 x 10(17) eV< E(nu) < 2 x 10(19) eV.

Measurement of the depth of maximum of extensive air showers above 10(18) eV

We describe the measurement of the depth of maximum, X(max), of the longitudinal development of air showers induced by cosmic rays. Almost 4000 events above 10(18) eV observed by the fluorescence detector of the Pierre Auger Observatory in coincidence with at least one surface detector station are selected for the analysis. The average shower maximum was found to evolve with energy at a rate of (106 +/- 35-21) g/cm(2)/decade below 10(18.24) +/- (0.05) eV, and d24 +/- 3 g/cm(2)/ecade above this energy. The measured shower-to-shower fluctuations decrease from about 55 to 26 g/cm(2). The interpretation of these results in terms of the cosmic ray mass composition is briefly discussed.

Limit on the diffuse flux of ultrahigh energy tau neutrinos with the surface detector of the Pierre Auger Observatory

Data collected at the Pierre Auger Observatory are used to establish an upper limit on the diffuse flux of tau neutrinos in the cosmic radiation. Earth-skimming nu(tau) may interact in the Earth's crust and produce a tau lepton by means of charged-current interactions. The tau lepton may emerge from the Earth and decay in the atmosphere to produce a nearly horizontal shower with a typical signature, a persistent electromagnetic component even at very large atmospheric depths. The search procedure to select events induced by tau decays against the background of normal showers induced by cosmic rays is described. The method used to compute the exposure for a detector continuously growing with time is detailed. Systematic uncertainties in the exposure from the detector, the analysis, and the involved physics are discussed. No tau neutrino candidates have been found. For neutrinos in the energy range 2x10(17) eV < E(nu)< 2x10(19) eV, assuming a diffuse spectrum of the form E(nu)(-2), data collected between 1 January 2004 and 30 April 2008 yield a 90% confidence-level upper limit of E(nu)(2)dN(nu tau)/dE(nu)< 9x10(-8) GeV cm(-2) s(-1) sr(-1).

Stretching and folding mechanism in foams

We have described the stretching and folding of foams in a vertical Hele-Shaw cell containing air and a surfactant solution, from a sequence of upside-down flips. Besides the firactal dimension of the foam, we have observed the logistic growth for the soap film length. The stretching and folding mechanism is present during the foam formation, and this mechanism is observed even after the foam has reached its respective maximum fractal dimension. Observing the motion of bubbles inside the foam, large bubbles present power spectrum associated with random walk motion in both directions, while the small bubbles are scattered like balls in a Galton board. (C) 2008 Published by Elsevier B.V.

Root Secretion of Defense-related Proteins Is Development-dependent and Correlated with Flowering Time

Proteins found in the root exudates are thought to play a role in the interactions between plants and soil organisms. To gain a better understanding of protein secretion by roots, we conducted a systematic proteomic analysis of the root exudates of Arabidopsis thaliana at different plant developmental stages. In total, we identified 111 proteins secreted by roots, the majority of which were exuded constitutively during all stages of development. However, defense-related proteins such as chitinases, glucanases, myrosinases, and others showed enhanced secretion during flowering. Defense-impaired mutants npr1-1 and NahG showed lower levels of secretion of defense proteins at flowering compared with the wild type. The flowering-defective mutants fca-1, stm-4, and co-1 showed almost undetectable levels of defense proteins in their root exudates at similar time points. In contrast, root secretions of defense-enhanced cpr5-2 mutants showed higher levels of defense proteins. The proteomics data were positively correlated with enzymatic activity assays for defense proteins and with in silico gene expression analysis of genes specifically expressed in roots of Arabidopsis. In conclusion, our results show a clear correlation between defense-related proteins secreted by roots and flowering time.

Consensus Statement: Chromosomal Microarray Is a First-Tier Clinical Diagnostic Test for Individuals with Developmental Disabilities or Congenital Anomalies

Chromosomal microarray (CMA) is increasingly utilized for genetic testing of individuals with unexplained developmental delay/intellectual disability (DD/ID), autism spectrum disorders (ASD), or multiple congenital anomalies (MCA). Performing CMA and G-banded karyotyping on every patient substantially increases the total cost of genetic testing. The International Standard Cytogenomic Array (ISCA) Consortium held two international workshops and conducted a literature review of 33 studies, including 21,698 patients tested by CMA. We provide an evidence-based summary of clinical cytogenetic testing comparing CMA to G-banded karyotyping with respect to technical advantages and limitations, diagnostic yield for various types of chromosomal aberrations, and issues that affect test interpretation. CMA offers a much higher diagnostic yield (15%-20%) for genetic testing of individuals with unexplained DD/ID, ASD, or MCA than a G-banded karyotype (similar to 3%, excluding Down syndrome and other recognizable chromosomal syndromes), primarily because of its higher sensitivity for submicroscopic deletions and duplications. Truly balanced rearrangements and low-level mosaicism are generally not detectable by arrays, but these are relatively infrequent causes of abnormal phenotypes in this population (<1%). Available evidence strongly supports the use of CMA in place of G-banded karyotyping as the first-tier cytogenetic diagnostic test for patients with DD/ID, ASD, or MCA. G-banded karyotype analysis should be reserved for patients with obvious chromosomal syndromes (e.g., Down syndrome), a family history of chromosomal rearrangement, or a history of multiple miscarriages.

Potential sites of triple-helical nucleic acid formation in chromosomes of Rhynchosciara (Diptera: Sciaridae) and Drosophila melanogaster

Antibodies to specific nucleic acid conformations are amongst the methods that have allowed the study of non-canonical (Watson-Crick) DNA structures in higher organisms. In this work, the structural limitations for the immunological detection of DNA.RNA hybrid duplexes were examined using specific RNA homopolymers as probes for homopolymer polydeoxyadenylic acid (poly(dA)).polydeoxythymidylic acid (poly(dT))-rich regions of Rhynchosciara americana (Diptera: Sciaridae) chromosomes. Anti-DNA.RNA duplexes did not react with the complex formed between chromosomal poly(dA) and exogenous polyuridylic acid (poly(rU)). Additionally, poly(rU) prevented the detection of polyadenylic acid.poly(dT) hybrid duplexes preformed in situ. These results raised the possibility that three-stranded structures rather than duplexes were formed in chromosomal sites. To test this hypothesis, the specificity of antibodies to triple-helical nucleic acids was reassessed employing distinct nucleic acid configurations. These antibodies were raised to the poly(dA).poly(rU).poly(rU) complex and have been used here for the first time in immunocytochemistry. Anti-triplex antibodies recognised the complex poly(dA).poly(rU).poly(rU) assembled with poly(rU) in poly(dA).poly(dT)-rich homopolymer regions of R. americana chromosomes. The antibodies could not detect short triplex stretches, suggesting the existence of constraints for triple-helix detection, probably related to triplex tract length. In addition, anti-poly(dA).poly(rU).poly(rU) antibodies reacted with the pericentric heterochromatin of RNase-treated polytene chromosomes of R. americana and Drosophila melanogaster. In apparent agreement with data obtained in cell types from other organisms, the results of this work suggest that significant triple-helix DNA extensions can be formed in pericentric regions of these species.

Upper limit on the cosmic-ray photon flux above 10(19) eV using the surface detector of the Pierre Auger Observatory

A method is developed to search for air showers initiated by photons using data recorded by the surface detector of the Auger Observatory. The approach is based on observables sensitive to the longitudinal shower development, the signal risetime and the curvature of the shower front. Applying this method to the data, tipper limits on the flux of photons of 3.8 x 10(-3), 2.5 x 10(-3), and 2.2 x 10(-3) km(-2) sr(-1) yr(-1) above 10(19) eV, 2 x 10(19) eV, and 4 x 10(19) eV are derived, with corresponding limits on the fraction of photons being 2.0%, 5.1%, and 31% (all limits at 95% c.l.). These photon limits disfavor certain exotic models of sources of cosmic rays. The results also show that the approach adopted by the Auger Observatory to calibrate the shower energy is not strongly biased by a contamination from photons. (C) 2008 Elsevier B.V. All rights reserved.

Correlation of the highest-energy cosmic rays with the positions of nearby active galactic nuclei

Data collected by the Pierre Auger Observatory provide evidence for anisotropy in the arrival directions of the cosmic rays with the highest-energies, which are correlated with the positions of relatively nearby active galactic nuclei (AGN) [Pierre Auger Collaboration, Science 318 (2007) 938]. The correlation has maximum significance for cosmic rays with energy greater than similar to 6 x 10(19) eV and AGN at a distance less than similar to 75 Mpc. We have confirmed the anisotropy at a confidence level of more than 99% through a test with parameters specified a priori, using an independent data set. The observed correlation is compatible with the hypothesis that cosmic rays with the highest-energies originate from extra-galactic sources close enough so that their flux is not significantly attenuated by interaction with the cosmic background radiation (the Greisen-Zatsepin-Kuz`min effect). The angular scale of the correlation observed is a few degrees, which suggests a predominantly light composition unless the magnetic fields are very weak outside the thin disk of our galaxy. Our present data do not identify AGN as the sources of cosmic rays unambiguously, and other candidate sources which are distributed as nearby AGN are not ruled out. We discuss the prospect of unequivocal identification of individual sources of the highest-energy cosmic rays within a few years of continued operation of the Pierre Auger Observatory. (C) 2008 Elsevier B.V. All rights reserved.

Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs ""radio-hybrid"" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluorescence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request. (C) 2011 Elsevier B.V. All rights reserved.

Trigger and aperture of the surface detector array of the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above 3 x 10(18) eV, for all zenith angles between 0 degrees and 60 degrees, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance. (C) 2009 Elsevier B.V. All rights reserved.

Atmospheric effects on extensive air showers observed with the surface detector of the Pierre Auger observatory

Atmospheric parameters, Such as pressure (P), temperature (T) and density (rho proportional to P/T), affect the development of extensive air showers initiated by energetic cosmic rays. We have Studied the impact of atmospheric variations on extensive air showers by means of the surface detector of the Pierre Auger Observatory. The rate of events shows a similar to 10% seasonal modulation and similar to 2% diurnal one. We find that the observed behaviour is explained by a model including the effects associated with the variations of P and rho. The former affects the longitudinal development of air showers while the latter influences the Moliere radius and hence the lateral distribution of the shower particles. The model is validated with full simulations of extensive air showers using atmospheric profiles measured at the site of the Pierre Auger Observatory. (C) 2009 Elsevier B.V. All rights reserved.