MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.

Quantized vortices in mixed 3He-4He drops

Using density functional theory, we investigate the structure of mixed 3HeN3-4HeN4 droplets with an embedded impurity (Xe atom or HCN molecule) which pins a quantized vortex line. We find that the dopant+vortex+4HeN4 complex, which in a previous work [F. Dalfovo et al., Phys. Rev. Lett. 85, 1028 (2000)] was found to be energetically stable below a critical size Ncr, is robust against the addition of 3He. While 3He atoms are distributed along the vortex line and on the surface of the 4He drop, the impurity is mostly coated by 4He atoms. Results for N4 = 500 and a number of 3He atoms ranging from 0 to 100 are presented, and the binding energy of the dopant to the vortex line is determined.

Nonrelativistic bound states at finite temperature. II. Muonic hydrogen

We illustrate how to apply modern effective field-theory techniques and dimensional regularization to factorize the various scales, which appear in QED bound states at finite temperature. We focus here on the muonic hydrogen atom. Vacuum polarization effects make the physics of this atom at finite temperature very close to that of heavy quarkonium states. We comment on the implications of our results for these states in the quark gluon plasma. In particular, we estimate the effects of a finite-charm quark mass in the dissociation temperature of bottomonium.

Distorted-wave calculation of cross sections for inner-shell ionization by electron and positron impact

The relativistic distorted-wave Born approximation is used to calculate differential and total cross sections for inner shell ionization of neutral atoms by electron and positron impact. The target atom is described within the independent-electron approximation using the self-consistent Dirac-Fock-Slater potential. The distorting potential for the projectile is also set equal to the Dirac-Fock-Slater potential. For electrons, this guarantees orthogonality of all the orbitals involved and simplifies the calculation of exchange T-matrix elements. The interaction between the projectile and the target electrons is assumed to reduce to the instantaneous Coulomb interaction. The adopted numerical algorithm allows the calculation of differential and total cross sections for projectiles with kinetic energies ranging from the ionization threshold up to about ten times this value. Algorithm accuracy and stability are demonstrated by comparing differential cross sections calculated by our code with the distorting potential set to zero with equivalent results generated by a more robust code that uses the conventional plane-wave Born approximation. Sample calculation results are presented for ionization of K- and L-shells of various elements and compared with the available experimental data.

Phosphorus transformation in poultry litter and litter-treated Oxisol of Brazil assessed by 31P-NMR and wet chemical fractionation

Large quantities of poultry litter are being produced in Brazil, which contain appreciable amounts of phosphorus (P) that could be of environmental concern. To assess the immediate environmental threat, five poultry litters composed of diverse bedding material were incubated for 43 days under greenhouse conditions. The litters consisted of: coffee bean husk (CH); wood chips (WC); rice husk (RH); ground corn cobs (CC) and ground napier grass (NG) (Pennisetum purpureum Schum.), in which the change in forms of soluble P was evaluated using 31P NMR spectroscopy. On average, 80.2 and 19.8 % of the total P in the extract, respectively, accounted for the inorganic and organic forms before incubation and 48 % of the organic P was mineralized to inorganic P in 43 days of incubation. Wide variation in the organic P mineralization rate (from 82 % -WC to 4 % - NG) was observed among litters. Inorganic orthophosphate (99.9 %) and pyrophosphate (0.1 %) were the only inorganic P forms, whereas the organic P forms orthophosphate monoesters (76.3 %) and diester (23.7 %) were detected. Diester P compounds were mineralized almost completely in all litters, except in the CH litter, within the incubation period. Pyrophosphates contributed with less than 0.5% and remained unaltered during the incubation period. Wood-chip litter had a higher organic P (40 %) content and a higher diester: monoester ratio; it was therefore mineralized rapidly, within the first 15 days, achieving steady state by the 29th day. Distinct mineralization patterns were observed in the litter when incubated with a clayey Oxisol. The substantial decrease observed in the organic P fraction (Po) of the litter types followed the order: CH (45 %) > CC (25 %) > RH (13 %) ≈ NG (12 %) > WC (5 %), whereas the Pi fraction increased. Incubation of RH litter in soil slowed down the mineralization of organic P.

Isotope determination of sulfur by mass spectrometry in soil samples

Sulphur plays an essential role in plants and is one of the main nutrients in several metabolic processes. It has four stable isotopes (32S, 33S, 34S, and 36S) with a natural abundance of 95.00, 0.76, 4.22, and 0.014 in atom %, respectively. A method for isotopic determination of S by isotope-ratio mass spectrometry (IRMS) in soil samples is proposed. The procedure involves the oxidation of organic S to sulphate (S-SO4(2-)), which was determined by dry combustion with alkaline oxidizing agents. The total S-SO4(2-) concentration was determined by turbidimetry and the results showed that the conversion process was adequate. To produce gaseous SO2 gas, BaSO4 was thermally decomposed in a vacuum system at 900 ºC in the presence of NaPO3. The isotope determination of S (atom % 34S atoms) was carried out by isotope ratio mass spectrometry (IRMS). In this work, the labeled material (K2(34)SO4) was used to validate the method of isotopic determination of S; the results were precise and accurate, showing the viability of the proposed method.

Microscopic calculations of the excitation spectrum of one 3He impurity in liquid 4He

We calculate the chemical potential ¿0 and the effective mass m*/m3 of one 3He impurity in liquid 4He. First a variational wave function including two- and three-particle dynamical correlations is adopted. Triplet correlations bring the computed values of ¿0 very close to the experimental results. The variational estimate of m*/m3 includes also backflow correlations between the 3He atom and the particles in the medium. Different approximations for the three-particle distribution function give almost the same values for m*/m3. The variational approach underestimates m*/m3 by ~10% at all of the considered densities. Correlated-basis perturbation theory is then used to improve the wave function to include backflow around the particles of the medium. The perturbative series built up with one-phonon states only is summed up to infinite order and gives results very close to the variational ones. All the perturbative diagrams with two independent phonons have then been summed to compute m*/m3. Their contribution depends to some extent on the form used for the three-particle distribution function. When the scaling approximation is adopted, a reasonable agreement with the experimental results is achieved.

Charged-particle interaction with liquids: Ripplon excitations

We calculate the ripplon field contribution to the self-energy of an electron exterior to a liquid for planar and spherical geometries. We compare the full dielectric calculation of the electron-liquid interaction with the simpler alternative method consisting of integrating the electron-atom static-induced-dipolar potential through the whole liquid volume. We obtain good agreement between both methods for a nonpolar liquid such as 4He but differences up to 40% for a polar liquid such as water. We study the conditions under which the ripplon contribution to the self-energy is a perturbation. For an electron moving parallel to a planar liquid surface, we calculate the ripplon contribution to its stopping power. For this dynamical case, we conclude that the alternative method is a good approximation even for polar liquids.

Calorimetry of hydrogen desorption from a-Si nanoparticles

The process of hydrogen desorption from amorphous silicon (a-Si) nanoparticles grown by plasma-enhanced chemical vapor deposition (PECVD) has been analyzed by differential scanning calorimetry (DSC), mass spectrometry, and infrared spectroscopy, with the aim of quantifying the energy exchanged. Two exothermic peaks centered at 330 and 410 C have been detected with energies per H atom of about 50 meV. This value has been compared with the results of theoretical calculations and is found to agree with the dissociation energy of Si-H groups of about 3.25 eV per H atom, provided that the formation energy per dangling bond in a-Si is about 1.15 eV. It is shown that this result is valid for a-Si:H films, too.

Evolución: de la especie humana al cyborg

Resumen: Gracias a la evolución biológica y cultural el ser humano ha sido capaz de adaptarse al ambiente. Una opción para mejorar las prestaciones humanas podría ser la obtención de un cyborg o ser humano biónico en el cual la parte inorgánica fuese un chip computacional conectado directamente al cerebro. Esta y otras posibilidades se debaten en el presente artículo. Palabras clave: Evolución biológica, selección natural, evolución cultural, cyborg, biónico, inteligencia.

dKDM5/LID regulates H3K4me3 dynamics at the transcription-start site (TSS) of actively transcribed developmental genes

H3K4me3 is a histone modification that accumulates at the transcription-start site (TSS) of active genes and is known to be important for transcription activation. The way in which H3K4me3 is regulated at TSS and the actual molecular basis of its contribution to transcription remain largely unanswered. To address these questions, we have analyzed the contribution of dKDM5/LID, the main H3K4me3 demethylase in Drosophila, to the regulation of the pattern of H3K4me3. ChIP-seq results show that, at developmental genes, dKDM5/LID localizes at TSS and regulates H3K4me3. dKDM5/LID target genes are highly transcribed and enriched in active RNApol II and H3K36me3, suggesting a positive contribution to transcription. Expression-profiling show that, though weakly, dKDM5/LID target genes are significantly downregulated upon dKDM5/LID depletion. Furthermore, dKDM5/LID depletion results in decreased RNApol II occupancy, particularly by the promoter-proximal Pol lloser5 form. Our results also show that ASH2, an evolutionarily conserved factor that locates at TSS and is required for H3K4me3, binds and positively regulates dKDM5/LID target genes. However, dKDM5/LID and ASH2 do not bind simultaneously and recognize different chromatin states, enriched in H3K4me3 and not, respectively. These results indicate that, at developmental genes, dKDM5/LID and ASH2 coordinately regulate H3K4me3 at TSS and that this dynamic regulation contributes to transcription.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

Rapid fluidic exchange microsystem for recording of fast ion channel kinetics in Xenopus oocytes.

We present a new lab-on-a-chip system for electrophysiological measurements on Xenopus oocytes. Xenopus oocytes are widely used host cells in the field of pharmacological studies and drug development. We developed a novel non-invasive technique using immobilized non-devitellinized cells that replaces the traditional "two-electrode voltage-clamp" (TEVC) method. In particular, rapid fluidic exchange was implemented on-chip to allow recording of fast kinetic events of exogenous ion channels expressed in the cell membrane. Reducing fluidic exchange times of extracellular reagent solutions is a great challenge with these large millimetre-sized cells. Fluidic switching is obtained by shifting the laminar flow interface in a perfusion channel under the cell by means of integrated poly-dimethylsiloxane (PDMS) microvalves. Reagent solution exchange times down to 20 ms have been achieved. An on-chip purging system allows to perform complex pharmacological protocols, making the system suitable for screening of ion channel ligand libraries. The performance of the integrated rapid fluidic exchange system was demonstrated by investigating the self-inhibition of human epithelial sodium channels (ENaC). Our results show that the response time of this ion channel to a specific reactant is about an order of magnitude faster than could be estimated with the traditional TEVC technique.

X-ray photoelectron spectroscopy of oxygen adsorbates on Al(111): Theory experiment

We present the result of polar angle resolved x¿ray photoemission spectroscopy on Al(111)/O and cluster calculations of the O(1s) binding energy (BE) for various model situations. In the experimental data two O(1s) peaks are observed, separated by 1.3 eV. The angular behavior (depth¿resolution) could indicate that the lower BE peak is associated with an O atom under the surface, and the higher BE peak with an O atom above the surface. Equally, it could indicate oxygen islands on the surface where the perimeter atoms have a higher O(1s) BE than the interior atoms. The cluster calculations show that the former interpretation cannot be correct, since an O ads below the surface has a higher calculated O(1s) BE than one above. Cluster calculations simulating oxygen islands are, however, consistent with the experimental data.