Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

Overfitting Bayesian mixture models with an unknown number of components

This paper proposes solutions to three issues pertaining to the estimation of finite mixture models with an unknown number of components: the non-identifiability induced by overfitting the number of components, the mixing limitations of standard Markov Chain Monte Carlo (MCMC) sampling techniques, and the related label switching problem. An overfitting approach is used to estimate the number of components in a finite mixture model via a Zmix algorithm. Zmix provides a bridge between multidimensional samplers and test based estimation methods, whereby priors are chosen to encourage extra groups to have weights approaching zero. MCMC sampling is made possible by the implementation of prior parallel tempering, an extension of parallel tempering. Zmix can accurately estimate the number of components, posterior parameter estimates and allocation probabilities given a sufficiently large sample size. The results will reflect uncertainty in the final model and will report the range of possible candidate models and their respective estimated probabilities from a single run. Label switching is resolved with a computationally light-weight method, Zswitch, developed for overfitted mixtures by exploiting the intuitiveness of allocation-based relabelling algorithms and the precision of label-invariant loss functions. Four simulation studies are included to illustrate Zmix and Zswitch, as well as three case studies from the literature. All methods are available as part of the R package Zmix, which can currently be applied to univariate Gaussian mixture models.

Thermoanalysis of ammonium perchlorate and sodium nitrate mixture

A binary mixture of ammonium perchlorate-sodium nitrate in molar proportion undergoes partial fusion at 223°C and the transformation of the mixture to sodium perchlorate-ammonium nitrate occurs in the broad endothermic region. The mixture was heated and quenched at various temperatures in a differential thermal analysis assembly. Thermogravimetric analysis, X-ray diffraction, and infrared spectroscopic techniques were used to determine the composition of the quenched sample in order to explain the overall thermal phenomenon. Visual observations of the morphological changes that occur during the course of heating were made using a hot-stage microscope, 30–350°C.

Representations and properties of para-Bose oscillator operators. I. Energy position and momentum eigenstates

Para-Bose commutation relations are related to the SL(2,R) Lie algebra. The irreducible representation [script D]alpha of the para-Bose system is obtained as the direct sum Dbeta[direct-sum]Dbeta+1/2 of the representations of the SL(2,R) Lie algebra. The position and momentum eigenstates are then obtained in this representation [script D]alpha, using the matrix mechanical method. The orthogonality, completeness, and the overlap of these eigenstates are derived. The momentum eigenstates are also derived using the wave mechanical method by specifying the domain of the definition of the momentum operator in addition to giving it a formal differential expression. By a careful consideration in this manner we find that the two apparently different solutions obtained by Ohnuki and Kamefuchi in this context are actually unitarily equivalent. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Strong-coupling expansion for the momentum distribution of the Bose-Hubbard model with benchmarking against exact numerical results

A strong-coupling expansion for the Green's functions, self-energies, and correlation functions of the Bose-Hubbard model is developed. We illustrate the general formalism, which includes all possible (normal-phase) inhomogeneous effects in the formalism, such as disorder or a trap potential, as well as effects of thermal excitations. The expansion is then employed to calculate the momentum distribution of the bosons in the Mott phase for an infinite homogeneous periodic system at zero temperature through third order in the hopping. By using scaling theory for the critical behavior at zero momentum and at the critical value of the hopping for the Mott insulator–to–superfluid transition along with a generalization of the random-phase-approximation-like form for the momentum distribution, we are able to extrapolate the series to infinite order and produce very accurate quantitative results for the momentum distribution in a simple functional form for one, two, and three dimensions. The accuracy is better in higher dimensions and is on the order of a few percent relative error everywhere except close to the critical value of the hopping divided by the on-site repulsion. In addition, we find simple phenomenological expressions for the Mott-phase lobes in two and three dimensions which are much more accurate than the truncated strong-coupling expansions and any other analytic approximation we are aware of. The strong-coupling expansions and scaling-theory results are benchmarked against numerically exact quantum Monte Carlo simulations in two and three dimensions and against density-matrix renormalization-group calculations in one dimension. These analytic expressions will be useful for quick comparison of experimental results to theory and in many cases can bypass the need for expensive numerical simulations.

Modelling predicts positive and negative interactions between three Australian tropical tree species in monoculture and binary mixture

Computer modelling promises to be an important tool for analysing and predicting interactions between trees within mixed species forest plantations. This study explored the use of an individual-based mechanistic model as a predictive tool for designing mixed species plantations of Australian tropical trees. The `spatially explicit individually based-forest simulator' (SeXI-FS) modelling system was used to describe the spatial interaction of individual tree crowns within a binary mixed-species experiment. The three-dimensional model was developed and verified with field data from three forest tree species grown in tropical Australia. The model predicted the interactions within monocultures and binary mixtures of Flindersia brayleyana, Eucalyptus pellita and Elaeocarpus grandis, accounting for an average of 42% of the growth variation exhibited by species in different treatments. The model requires only structural dimensions and shade tolerance as species parameters. By modelling interactions in existing tree mixtures, the model predicted both increases and reductions in the growth of mixtures (up to +/-50% of stem volume at 7 years) compared to monocultures. This modelling approach may be useful for designing mixed tree plantations.

NOx Reduction and Desorption Studies Under Electric Discharge Plasma Using a Simulated Gas Mixture: A Case Study on the Effect of Corona Electrodes

In this study, reduction and desorption of oxides of nitrogen (NOx) were conducted using an electrical discharge plasma technique. The study was carried out using a simulated gas mixture to explore the possibility of re-generation of used adsorbents by a nonthermal plasma desorption technique. Three different types of corona electrodes, namely, pipe, helical wire, and straight wire, were used for analyzing their effectiveness in NOx reduction/desorption. The pipe-type corona electrode exhibited a nitric oxide (NO) conversion of 50%, which is 1.5 times that of the straight-wire-type electrode at an energy density of 175J/L. The helical-wire-type corona electrode exhibited a NOx desorption efficiency almost 4 times that of the pipe-type electrode,indicating the possibility that corona-generated species play a crucial role in desorption.

Comparative productivity of irrigated short-term ryegrass (Lolium multiflorum) pasture receiving nitrogen, grown alone or in a mixture with white (Trifolium repens) and Persian (T. resupinatum) clovers

Dairy farms in subtropical Australia use irrigated, annually sown short-term ryegrass (Lolium multiflorum) or mixtures of short-term ryegrass and white (Trifolium repens) and Persian (shaftal) (T. resupinatum) clover during the winter-spring period in all-year-round milk production systems. A series of small plot cutting experiments was conducted in 3 dairying regions (tropical upland, north Queensland, and subtropical southeast Queensland and northern New South Wales) to determine the most effective rate and frequency of application of nitrogen (N) fertiliser. The experiments were not grazed, nor was harvested material returned to the plots, after sampling. Rates up to 100 kg N/ha.month (as urea or calcium ammonium nitrate) and up to 200 kg N/ha every 2 months (as urea) were applied to pure stands of ryegrass in 1991. In 1993 and 1994, urea, at rates up to 150 kg N/ha.month and to 200 kg N/ha every 2 months, was applied to pure stands of ryegrass; urea, at rates up to 50 kg N/ha.month, was also applied to ryegrass-clover mixtures. The results indicate that applications of 50-85 kg N/ha.month can be recommended for short-term ryegrass pastures throughout the subtropics and tropical uplands of eastern Australia, irrespective of soil type. At this rate, dry matter yields will reach about 90% of their potential, forage nitrogen concentration will be increased, there is minimal risk to stock from nitrate poisoning and there will be no substantial increase in soil N. The rate of N for ryegrass-clover pastures is slightly higher than for pure ryegrass but, at these rates, the clover component will be suppressed. However, increased ryegrass yields and higher forage nitrogen concentrations will compensate for the reduced clover component. At application rates up to 100 kg N/ha.month, build-up of NO3--N and NH4+-N in soil was generally restricted to the surface layers (0-20 cm) of the soil, but there was a substantial increase throughout the soil profile at 150 kg N/ha.month. The build-up of NO3--N and NH4+-N was greater and was found at lower rates on the lighter soil compared with heavy clays. Generally, most of the soil N was in the NO3--N form and most was in the top 20 cm.

The performance of irrigated mixtures of tall fescue, ryegrass and white clover in subtropical Australia. 1. The effects of sowing mixture combinations, nitrogen and oversowing on establishment, productivity, botanical composition and persistence.

In the subtropics of Australia, the ryegrass component of irrigated perennial ryegrass (Lolium perenne) - white clover (Trifolium repens) pastures declines by approximately 40% in the summer following establishment, being replaced by summer-active C4 grasses. Tall fescue (Festuca arundinacea) is more persistent than perennial ryegrass and might resist this invasion, although tall fescue does not compete vigorously as a seedling. This series of experiments investigated the influence of ryegrass and tall fescue genotype, sowing time and sowing mixture as a means of improving tall fescue establishment and the productivity and persistence of tall fescue, ryegrass and white clover-based mixtures in a subtropical environment. Tall fescue frequency at the end of the establishment year decreased as the number of companion species sown in the mixture increased. Neither sowing mixture combinations nor sowing rates influenced overall pasture yield (of around 14 t/ha) in the establishment year but had a significant effect on botanical composition and component yields. Perennial ryegrass was less competitive than short-rotation ryegrass, increasing first-year yields of tall fescue by 40% in one experiment and by 10% in another but total yield was unaffected. The higher establishment-year yield (3.5 t/ha) allowed Dovey tall fescue to compete more successfully with the remaining pasture components than Vulcan (1.4 t/ha). Sowing 2 ryegrass cultivars in the mixture reduced tall fescue yields by 30% compared with a single ryegrass (1.6 t/ha), although tall fescue alone achieved higher yields (7.1 t/ha). Component sowing rate had little influence on composition or yield. Oversowing the ryegrass component into a 6-week-old sward of tall fescue and white clover improved tall fescue, white clover and overall yields in the establishment year by 83, 17 and 11%, respectively, but reduced ryegrass yields by 40%. The inclusion of red (T. pratense) and Persian (T. resupinatum) clovers and chicory (Cichorium intybus) increased first-year yields by 25% but suppressed perennial grass and clover components. Yields were generally maintained at around 12 t/ha/yr in the second and third years, with tall fescue becoming dominant in all 3 experiments. The lower tall fescue seeding rate used in the first experiment resulted in tall fescue dominance in the second year following establishment, whereas in Experiments 2 and 3 dominance occurred by the end of the first year. Invasion by the C4 grasses was relatively minor (<10%) even in the third year. As ryegrass plants died, tall fescue and, to a lesser extent, white clover increased as a proportion of the total sward. Treatment effects continued into the second, but rarely the third, year and mostly affected the yield of one of the components rather than total cumulative yield. Once tall fescue became dominant, it was difficult to re-introduce other pasture components, even following removal of foliage and moderate renovation. Severe renovation (reducing the tall fescue population by at least 30%) seems a possible option for redressing this situation.

Thermoanalysis of ammonium perchlorate and sodium nitrate mixture

A binary mixture of ammonium perchlorate-sodium nitrate in molar proportion undergoes partial fusion at 223°C and the transformation of the mixture to sodium perchlorate-ammonium nitrate occurs in the broad endothermic region. The mixture was heated and quenched at various temperatures in a differential thermal analysis assembly. Thermogravimetric analysis, X-ray diffraction, and infrared spectroscopic techniques were used to determine the composition of the quenched sample in order to explain the overall thermal phenomenon. Visual observations of the morphological changes that occur during the course of heating were made using a hot-stage microscope, 30–350°C.

Insulating phase of a disordered system: Fermi glass versus electron glass

The authors examine the critical divergence of the low-frequency conductivity of the noninteracting Fermi glass and interacting electron glass models of the insulating phase of a disordered system as the metallic phase is approached. Results for the two are found to be rather different, which can be tested experimentally. In particular, for the electron glass, there exists a nonvanishing contribution to the dielectric constants from the low-frequency (hopping) conductivity even at low temperatures, which scales with the high-frequency (optical) contribution, and diverges with the same exponent at the insulator-metal transition.

Representation and properties of para-Bose oscillator operators. II. Coherent states and the minimum uncertainty states

The energy, position, and momentum eigenstates of a para-Bose oscillator system were considered in paper I. Here we consider the Bargmann or the analytic function description of the para-Bose system. This brings in, in a natural way, the coherent states ||z;alpha> defined as the eigenstates of the annihilation operator ?. The transformation functions relating this description to the energy, position, and momentum eigenstates are explicitly obtained. Possible resolution of the identity operator using coherent states is examined. A particular resolution contains two integrals, one containing the diagonal basis ||z;alpha><−z;alpha||. We briefly consider the normal and antinormal ordering of the operators and their diagonal and discrete diagonal coherent state approximations. The problem of constructing states with a minimum value of the product of the position and momentum uncertainties and the possible alpha dependence of this minimum value is considered. Journal of Mathematical Physics is copyrighted by The American Institute of Physics.

Analysis-By-Synthesis Based Switched Transform Domain Split Vq Using Gaussian Mixture Model

Using analysis-by-synthesis (AbS) approach, we develop a soft decision based switched vector quantization (VQ) method for high quality and low complexity coding of wideband speech line spectral frequency (LSF) parameters. For each switching region, a low complexity transform domain split VQ (TrSVQ) is designed. The overall rate-distortion (R/D) performance optimality of new switched quantizer is addressed in the Gaussian mixture model (GMM) based parametric framework. In the AbS approach, the reduction of quantization complexity is achieved through the use of nearest neighbor (NN) TrSVQs and splitting the transform domain vector into higher number of subvectors. Compared to the current LSF quantization methods, the new method is shown to provide competitive or better trade-off between R/D performance and complexity.

Mixture Effects of Benzene, Toluene, Ethylbenzene and Xylenes (BTEX) to Lung Carcinoma Cells via a Hanging Drop Air Exposure System

A recently developed hanging drop air exposure system for toxicity studies of volatile chemicals was applied to evaluate the cell viability of lung carcinoma A549 cells after 1 h and 24 h of exposure to benzene, toluene, ethylbenzene and xylenes (BTEX) as individual compounds and mixtures of 4 or 6 components. The cellular chemical concentrations causing 50% reduction of cell viability (EC50) were calculated use a mass balance model and came to 17, 12, 11, 9, 4 and 4 mmol/kg cell dry weight for benzene, toluene, ethylbenzene, m-xylene, o-xylene and p-xylene respectively after 1 h of exposure. The EC50 decreased by a factor of four after 24 h of exposure. All mixture effects were best described by the mixture toxicity model of concentration addition, which is valid for chemicals with the same mode of action. Good agreement with the model predictions were found for benzene, toluene, ethylbenzene and m-xylene at four different representative fixed concentration ratios after 1 h of exposure but lower agreement to mixture prediction was obtained after 24 h of exposure. A recreated car exhaust mixture, which involved the contribution of the more toxic p-xylene and o-xylene, yielded an acceptable but lower quality prediction as well.