The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Unconventional metallic behavior and superconductivity in the K-Mo-O system

Transport properties and magnetization measurements of the K(x)MoO(2-delta) (0 <= x <= 0.25) compound are reported. The compound crystallizes in the oxygen deficient MoO(2) monoclinic structure with potassium atoms occupying interstitial positions. An unconventional metallic behavior with power-law temperature dependence is related to a magnetic ordering. Superconducting transition with small volume fraction is also observed near 7 K for a sample with low potassium composition.

Polyelectrolyte-protein complexation driven by charge regulation

The interplay between the biocolloidal characteristics (especially size and charge), pH, salt concentration and the thermal energy results in a unique collection of mesoscopic forces of importance to the molecular organization and function in biological systems. By means of Monte Carlo simulations and semi-quantitative analysis in terms of perturbation theory, we describe a general electrostatic mechanism that gives attraction at low electrolyte concentrations. This charge regulation mechanism due to titrating amino acid residues is discussed in a purely electrostatic framework. The complexation data reported here for interaction between a polyelectrolyte chain and the proteins albumin, goat and bovine alpha-lactalbumin, beta-lactoglobulin, insulin, k-casein, lysozyme and pectin methylesterase illustrate the importance of the charge regulation mechanism. Special attention is given to pH congruent to pI where ion-dipole and charge regulation interactions could overcome the repulsive ion-ion interaction. By means of protein mutations, we confirm the importance of the charge regulation mechanism, and quantify when the complexation is dominated either by charge regulation or by the ion-dipole term.

Seasonality Role on the Phenolics from Cultivated Baccharis dracunculifolia

Baccharis dracunculifolia is the source of Brazilian green propolis (BGP). Considering the broad spectrum of biological activities attributed to green proplis, B. dracunculifolia has a great potential for the development of new cosmetic and pharmaceutical products. In this work, the cultivation of 10 different populations of native B. dracunculifolia had been undertaken aiming to determine the role of seasonality on its phenolic compounds. For this purpose, fruits of this plant were collected from populations of 10 different regions, and 100 individuals of each population were cultivated in an experimental area of 1800 m(2). With respect to cultivation, the yields of dry plant, essential oil and crude extract were measured monthly resulting in mean values of 399 +/- 80 g, 0.6 +/- 0.1% and 20 +/- 4%, respectively. The HPLC analysis allowed detecting seven phenolic compounds: caffeic acid, ferulic acid, aromadendrin-4'-methyl ether (AME), isosakuranetin, artepillin C, baccharin and 2-dimethyl-6-carboxyethenyl-2H-1-benzopyran acid, which were the major ones throughout the 1-year monthly analysis. Caffeic acid was detected in all cultivated populations with mean of 4.0%. AME displayed the wide variation in relation to other compounds showing means values of 0.65 +/- 0.13% at last quarter. Isosakuranetin and artepillin C showed increasing concentrations with values between 0% and 1.4% and 0% and 1.09%, respectively. The obtained results allow suggesting that the best time for harvesting this plant, in order to obtain good qualitative and quantitative results for these phenolic compounds, is between December and April.

Allosteric Control of Gating Mechanisms Revisited: The Large Conductance Ca(2+)-Activated K(+) Channel

Large-conductance Ca(2+)-activated K(+) channels (BK) play a fundamental role in modulating membrane potential in many cell types. The gating of BK channels and its modulation by Ca(2+) and voltage has been the subject of intensive research over almost three decades, yielding several of the most complicated kinetic mechanisms ever proposed. A large number of open and closed states disposed, respectively, in two planes, named tiers, characterize these mechanisms. Transitions between states in the same plane are cooperative and modulated by Ca(2+). Transitions across planes are highly concerted and voltage-dependent. Here we reexamine the validity of the two-tiered hypothesis by restricting attention to the modulation by Ca(2+). Large single channel data sets at five Ca(2+) concentrations were simultaneously analyzed from a Bayesian perspective by using hidden Markov models and Markov-chain Monte Carlo stochastic integration techniques. Our results support a dramatic reduction in model complexity, favoring a simple mechanism derived from the Monod-Wyman-Changeux allosteric model for homotetramers, able to explain the Ca(2+) modulation of the gating process. This model differs from the standard Monod-Wyman-Changeux scheme in that one distinguishes when two Ca(2+) ions are bound to adjacent or diagonal subunits of the tetramer.

Spin alignment measurements of the K*(0)(892) and phi(1020) vector mesons in heavy ion collisions at root S(NN)=200 GeV

We present the first spin alignment measurements for the K*(0)(892) and phi(1020) vector mesons produced at midrapidity with transverse momenta up to 5 GeV/c at root s(NN) = 200 GeV at RHIC. The diagonal spin-density matrix elements with respect to the reaction plane in Au+Au collisions are rho(00) = 0.32 +/- 0.04 (stat) +/- 0.09 (syst) for the K*(0) (0.8 < p(T) < 5.0 GeV/c) and rho(00) = 0.34 +/- 0.02 (stat) +/- 0.03 (syst) for the phi (0.4 < p(T) < 5.0 GeV/c) and are constant with transverse momentum and collision centrality. The data are consistent with the unpolarized expectation of 1/3 and thus no evidence is found for the transfer of the orbital angular momentum of the colliding system to the vector-meson spins. Spin alignments for K(*0) and phi in Au+Au collisions were also measured with respect to the particle's production plane. The phi result, rho(00) = 0.41 +/- 0.02 (stat) +/- 0.04 (syst), is consistent with that in p+p collisions, rho(00) = 0.39 +/- 0.03 (stat) +/- 0.06 (syst), also measured in this work. The measurements thus constrain the possible size of polarization phenomena in the production dynamics of vector mesons.

Changes in temporomandibular joint disc position and form following Herbst and fixed orthodontic treatment

Objective: To determine the changes in the position and form of the temporomandibular joint articular disc in adolescents with Class II division 1 malocclusion and mandibular retrognathism treated with the Herbst appliance (phase I) and fixed orthodontic appliance (phase II). Materials and Methods: Thirty-two consecutive adolescents went through phase I of treatment and 23 completed phase II. The temporomandibular joints were evaluated qualitatively by means of magnetic resonance images at the beginning of treatment (T1), during phase I (T2), at the end of phase I (T3), and at the end of phase II (T4). Results: Significant changes in disc position were not observed with the mouth closed between T1 X T3 (P = .317), T3 X T4 (P = .287), or T1 X T4 (P = .261). At T2, on average, the disc was positioned regressively. With the mouth open, no difference was observed between T1 X T3 (P = .223) or T1 X T4 (P = .082). We did observe a significant difference between T3 X T4 (P < .05). Significant changes in the disc form were found with the mouth closed between T1 X T2 (P < .001) and T2 X T3 (P < .001). Conclusions: At the end of the two-phase treatment, in general terms, the position and form of the initial articular discs were maintained; however, in some temporomandibular joints some seemingly adverse effects were observed at T4. (Angle Orthod. 2010;80:843-852.)

Alpha element abundances and gradients in the Milky Way bulge from FLAMES-GIRAFFE spectra of 650 K giants

Aims. We present the analysis of the [alpha/Fe] abundance ratios for a large number of stars at several locations in the Milky Way bulge with the aim of constraining its formation scenario. Methods. We obtained FLAMES-GIRAFFE spectra (R = 22 500) at the ESO Very Large Telescope for 650 bulge red giant branch (RGB) stars and performed spectral synthesis to measure Mg, Ca, Ti, and Si abundances. This sample is composed of 474 giant stars observed in 3 fields along the minor axis of the Galactic bulge and at latitudes b = -4 degrees, b = -6 degrees, b = -12 degrees. Another 176 stars belong to a field containing the globular cluster NGC 6553, located at b = -3 degrees and 5 degrees away from the other three fields along the major axis. Stellar parameters and metallicities for these stars were presented in Zoccali et al. (2008, A&A, 486, 177). We have also re-derived stellar parameters and abundances for the sample of thick and thin disk red giants analyzed in Alves-Brito et al. (2010, A&A, 513, A35). Therefore using a homogeneous abundance database for the bulge, thick and thin disk, we have performed a differential analysis minimizing systematic errors, to compare the formation scenarios of these Galactic components. Results. Our results confirm, with large number statistics, the chemical similarity between the Galactic bulge and thick disk, which are both enhanced in alpha elements when compared to the thin disk. In the same context, we analyze [alpha/Fe] vs. [Fe/H] trends across different bulge regions. The most metal rich stars, showing low [alpha/Fe] ratios at b = -4 degrees disappear at higher Galactic latitudes in agreement with the observed metallicity gradient in the bulge. Metal-poor stars ([Fe/H] < -0.2) show a remarkable homogeneity at different bulge locations. Conclusions. We have obtained further constrains for the formation scenario of the Galactic bulge. A metal-poor component chemically indistinguishable from the thick disk hints for a fast and early formation for both the bulge and the thick disk. Such a component shows no variation, neither in abundances nor kinematics, among different bulge regions. A metal-rich component showing low [alpha/Fe] similar to those of the thin disk disappears at larger latitudes. This allows us to trace a component formed through fast early mergers (classical bulge) and a disk/bar component formed on a more extended timescale.

A FIRST CONSTRAINT ON THE THICK DISK SCALE LENGTH: DIFFERENTIAL RADIAL ABUNDANCES IN K GIANTS AT GALACTOCENTRIC RADII 4, 8, AND 12 kpc

Based on high-resolution spectra obtained with the MIKE spectrograph on the Magellan telescopes, we present detailed elemental abundances for 20 red giant stars in the outer Galactic disk, located at Galactocentric distances between 9 and 13 kpc. The outer disk sample is complemented with samples of red giants from the inner Galactic disk and the solar neighborhood, analyzed using identical methods. For Galactocentric distances beyond 10 kpc, we only find chemical patterns associated with the local thin disk, even for stars far above the Galactic plane. Our results show that the relative densities of the thick and thin disks are dramatically different from the solar neighborhood, and we therefore suggest that the radial scale length of the thick disk is much shorter than that of the thin disk. We make a first estimate of the thick disk scale length of L(thick) = 2.0 kpc, assuming L(thin) = 3.8 kpc for the thin disk. We suggest that radial migration may explain the lack of radial age, metallicity, and abundance gradients in the thick disk, possibly also explaining the link between the thick disk and the metal-poor bulge.

A binary engine fuelling HD 87643's complex circumstellar environment Determined using AMBER/VLTI imaging

Context. The star HD 87643, exhibiting the ""B[e] phenomenon"", has one of the most extreme infrared excesses for this object class. It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty material was tested through a combination of multi-wavelength high spatial resolution observations. Methods. We observed HD 87643 with high spatial resolution techniques, using the near-IR AMBER/VLTI interferometer with baselines ranging from 60 m to 130 m and the mid-IR MIDI/VLTI interferometer with baselines ranging from 25 m to 65 m. These observations are complemented by NACO/VLT adaptive-optics-corrected images in the K and L-bands, and ESO-2.2m optical Wide-Field Imager large-scale images in the B, V and R-bands. Results. We report the direct detection of a companion to HD 87643 by means of image synthesis using the AMBER/VLTI instrument. The presence of the companion is confirmed by the MIDI and NACO data, although with a lower confidence. The companion is separated by similar to 34 mas with a roughly north-south orientation. The period must be large (several tens of years) and hence the orbital parameters are not determined yet. Binarity with high eccentricity might be the key to interpreting the extreme characteristics of this system, namely a dusty circumstellar envelope around the primary, a compact dust nebulosity around the binary system and a complex extended nebula suggesting past violent ejections.

The B0.5 IVe CoRoT target HD 49330 II. Spectroscopic ground-based observations

Context. We present spectroscopic ground-based observations of the early Be star HD 49330 obtained simultaneously with the CoRoT-LRA1 run just before the burst observed in the CoRoT data. Aims. Ground-based spectroscopic observations of the early Be star HD 49330 obtained during the precursor phase and just before the start of an outburst allow us to disantangle stellar and circumstellar contributions and identify modes of stellar pulsations in this rapidly rotating star. Methods. Time series analysis (TSA) is performed on photospheric line profiles of He I and Si III by means of the least squares method. Results. We find two main frequencies f1 = 11.86 c d(-1) and f2 = 16.89 c d(-1) which can be associated with high order p-mode pulsations. We also detect a frequency f3 = 1.51 c d(-1) which can be associated with a low order g-mode. Moreover we show that the stellar line profile variability changed over the spectroscopic run. These results are in agreement with the results of the CoRoT data analysis, as shown in Huat et al. (2009). Conclusions. Our study of mid-and short-term spectroscopic variability allows the identification of p-and g-modes in HD 49330. It also allows us to display changes in the line profile variability before the start of an outburst. This brings new constraints for the seimic modelling of this star.

Expression of an activating mutation in the gene encoding the K(ATP) channel subunit Kir6.2 in mouse pancreatic beta cells recapitulates neonatal diabetes

Neonatal diabetes is a rare monogenic form of diabetes that usually presents within the first six months of life. It is commonly caused by gain-of-function mutations in the genes encoding the Kir6.2 and SUR1 subunits of the plasmalemmal ATP-sensitive K(+) (K(ATP)) channel. To better understand this disease, we generated a mouse expressing a Kir6.2 mutation (V59M) that causes neonatal diabetes in humans and we used Cre-lox technology to express the mutation specifically in pancreatic beta cells. These beta-V59M mice developed severe diabetes soon after birth, and by 5 weeks of age, blood glucose levels were markedly increased and insulin was undetectable. Islets isolated from beta-V59M mice secreted substantially less insulin and showed a smaller increase in intracellular calcium in response to glucose. This was due to a reduced sensitivity of K(ATP) channels in pancreatic beta cells to inhibition by ATP or glucose. In contrast, the sulfonylurea tolbutamide, a specific blocker of K(ATP) channels, closed K(ATP) channels, elevated intracellular calcium levels, and stimulated insulin release in beta-V59M beta cells, indicating that events downstream of K(ATP) channel closure remained intact. Expression of the V59M Kir6.2 mutation in pancreatic beta cells alone is thus sufficient to recapitulate the neonatal diabetes observed in humans. beta-V59M islets also displayed a reduced percentage of beta cells, abnormal morphology, lower insulin content, and decreased expression of Kir6.2, SUR1, and insulin mRNA. All these changes are expected to contribute to the diabetes of beta-V59M mice. Their cause requires further investigation.

The constancy of global regulation across a species: the concentrations of ppGpp and RpoS are strain-specific in Escherichia coli

Background: Sigma factors and the alarmone ppGpp control the allocation of RNA polymerase to promoters under stressful conditions. Both ppGpp and the sigma factor sigma(S) (RpoS) are potentially subject to variability across the species Escherichia coli. To find out the extent of strain variation we measured the level of RpoS and ppGpp using 31 E. coli strains from the ECOR collection and one reference K-12 strain. Results: Nine ECORs had highly deleterious mutations in rpoS, 12 had RpoS protein up to 7-fold above that of the reference strain MG1655 and the remainder had comparable or lower levels. Strain variation was also evident in ppGpp accumulation under carbon starvation and spoT mutations were present in several low-ppGpp strains. Three relationships between RpoS and ppGpp levels were found: isolates with zero RpoS but various ppGpp levels, strains where RpoS levels were proportional to ppGpp and a third unexpected class in which RpoS was present but not proportional to ppGpp concentration. High-RpoS and high-ppGpp strains accumulated rpoS mutations under nutrient limitation, providing a source of polymorphisms. Conclusions: The ppGpp and sigma(S) variance means that the expression of genes involved in translation, stress and other traits affected by ppGpp and/or RpoS are likely to be strain-specific and suggest that influential components of regulatory networks are frequently reset by microevolution. Different strains of E. coli have different relationships between ppGpp and RpoS levels and only some exhibit a proportionality between increasing ppGpp and RpoS levels as demonstrated for E. coli K-12.

Chain motifs: The tails and handles of complex networks

A great part of the interest in complex networks has been motivated by the presence of structured, frequently nonuniform, connectivity. Because diverse connectivity patterns tend to result in distinct network dynamics, and also because they provide the means to identify and classify several types of complex network, it becomes important to obtain meaningful measurements of the local network topology. In addition to traditional features such as the node degree, clustering coefficient, and shortest path, motifs have been introduced in the literature in order to provide complementary descriptions of the network connectivity. The current work proposes a different type of motif, namely, chains of nodes, that is, sequences of connected nodes with degree 2. These chains have been subdivided into cords, tails, rings, and handles, depending on the type of their extremities (e.g., open or connected). A theoretical analysis of the density of such motifs in random and scale-free networks is described, and an algorithm for identifying these motifs in general networks is presented. The potential of considering chains for network characterization has been illustrated with respect to five categories of real-world networks including 16 cases. Several interesting findings were obtained, including the fact that several chains were observed in real-world networks, especially the world wide web, books, and the power grid. The possibility of chains resulting from incompletely sampled networks is also investigated.

Detailed study of the K(L)->pi(0)pi(0)pi(0) Dalitz plot

Using a sample of 68.3x10(6) K(L)->pi(0)pi(0)pi(0) decays collected in 1996-1999 by the KTeV (E832) experiment at Fermilab, we present a detailed study of the K(L)->pi(0)pi(0)pi(0) Dalitz plot density. We report the first observation of interference from K(L)->pi(+)pi(-)pi(0) decays in which pi(+)pi(-) rescatters to pi(0)pi(0) in a final-state interaction. This rescattering effect is described by the Cabibbo-Isidori model, and it depends on the difference in pion scattering lengths between the isospin I=0 and I=2 states, a(0)-a(2). Using the Cabibbo-Isidori model, and fixing (a(0)-a(2))m(pi)(+)=0.268 +/- 0.017 as measured by the CERN-NA48 collaboration, we present the first measurement of the K(L)->pi(0)pi(0)pi(0) quadratic slope parameter that accounts for the rescattering effect: h(000)=(+0.59 +/- 0.20(stat)+/- 0.48(syst)+/- 1.06(ext))x10(-3), where the uncertainties are from data statistics, KTeV systematic errors, and external systematic errors. Fitting for both h(000) and a(0)-a(2), we find h(000)=(-2.09 +/- 0.62(stat)+/- 0.72(syst)+/- 0.28(ext))x10(-3), and m(pi)(+)(a(0)-a(2))=0.215 +/- 0.014(stat)+/- 0.025(syst)+/- 0.006(ext); our value for a(0)-a(2) is consistent with that from NA48.