Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

Development and mechanical testing of a short intramedullary nail for fixation of femoral rotational osteotomy in cerebral palsy patients

Background: Rotational osteotomy is frequently indicated to correct excessive femoral anteversion in cerebral palsy patients. Angled blade plate is the standard fixation device used when performed in the proximal femur, but extensile exposure is required for plate accommodation. The authors developed a short locked intramedullary nail to be applied percutaneously in the fixation of femoral rotational osteotomies in children with cerebral palsy and evaluated its mechanical properties. Methods: The study was divided into three stages. In the first part, a prototype was designed and made based on radiographic measurements of the femoral medullary canal of ten-year-old patients. In the second, synthetic femoral models based on rapid-prototyping of 3D reconstructed images of patients with cerebral palsy were obtained and were employed to adjust the nail prototype to the morphological changes observed in this disease. In the third, rotational osteotomies were simulated using synthetic femoral models stabilized by the nail and by the AO-ASIF fixed-angle blade plate. Mechanical testing was done comparing both devices in bending-compression and torsion. Results: The authors observed proper adaptation of the nail to normal and morphologically altered femoral models, and during the simulated osteotomies. Stiffness in bending-compression was significantly higher in the group fixed by the plate (388.97 +/- 57.25 N/mm) than in that fixed by the nail (268.26 +/- 38.51 N/mm) as torsional relative stiffness was significantly higher in the group fixed by the plate (1.07 +/- 0.36 Nm/degrees) than by the nail (0.35 +/- 0.13 Nm/degrees). Conclusions: Although the device presented adequate design and dimension to fit into the pediatric femur, mechanical tests indicated that the nail was less stable than the blade plate in bending-compression and torsion. This may be a beneficial property, and it can be attributed to the more flexible fixation found in intramedullary devices.

Effects of Er:YAG Laser Irradiation on the Microtensile Bond Strength to Bleached Enamel

Objective: The objective of this study was to evaluate the influence of different Er:YAG laser (lambda = 2.94 mu m) energy parameters on the microtensile bond strength (mu TBS) and superficial morphology of bovine enamel bleached with 16% carbamide peroxide. Background: Laser irradiation could improve adhesion to bleached enamel surfaces. Methods: Sixty bovine enamel blocks (7x3x3 mm(3)) were randomly assigned to six groups according to enamel preparation procedures (n = 10): G1-bleaching and Er:YAG laser irradiation with 25.52 J/cm(2) (laser A, LA); G2-bleaching and Er:YAG laser irradiation with 4.42J/cm(2) (laser B, LB); G3-bleaching; G4-Er:YAG laser irradiation with 25.52 J/cm(2); G5-Er:YAG laser irradiation with 4.42J/cm(2); G6-control, no treatment. G1 to G3 were bleached for 6 h during 21 days. Afterwards, enamel surfaces in all groups were slightly abraded with 600-grit SiC papers and G1, G2, G4 and G5 were irradiated according to each protocol. Enamel blocks were then restored with an etch-and-rinse adhesive system and a 4-mm thick composite buildup was made in two increments (n = 9). After 24 h, restored blocks were serially sectioned with a cross-section area of similar to 1 mm(2) at the bonded interface and tested in tension in a universal testing machine (1 mm/min). Failure mode was determined at a magnification of x100 using a stereomicroscope. One treated block of each group was selected for scanning electron microscopy (SEM) analysis. mu TBS data were analyzed by two-way ANOVA and no statistical differences were observed among groups. Results: Mean bond strengths (SD) in MPa were: G1-30.4(6.2); G2-27.9(8.5); G3-32.3(3.9); G4-23.7(5.8); G5-29.3(6.0); G6-29.1(6.1). A large number of adhesive failures was recorded for bleached and irradiated enamel surfaces. Conclusions: Bleached enamel surfaces mu TBS values were not significantly different from those of unbleached enamel. Even though Er:YAG laser irradiation with both parameters had no influence on mu TBS for bleached and unbleached enamel, SEM analysis revealed that Er:YAG laser irradiation with 25.52J/cm(2) should not be recommended, as enamel ablation was observed, whereas irradiation with 4.42J/cm(2) did not promote any remarkable changes on enamel surface.

Increased Occurrence of Dental Anomalies Associated with Second-Premolar Agenesis

Objective: To evaluate the prevalence of dental anomalies in patients with agenesis of second premolars and compare the findings with the prevalence of these anomalies in the general population. Materials and Methods: A Brazilian sample of 203 patients aged 8 to 22 years was selected. All patients presented agenesis of at least one second premolar. Panoramic and periapical radiographs and dental casts were used to analyze the presence of other associated dental anomalies, including agenesis of other permanent teeth, ectopia of unerupted permanent teeth, infraocclusion of deciduous molars, microdontia of maxillary lateral incisors, and supernumerary teeth. The occurrence of these anomalies was compared with occurrence data previously reported for the general population. Statistical testing was performed using the chi-square test (P < .05) and the odds ratio. Results: The sample with agenesis of at least one second premolar presented a significantly increased prevalence rate of permanent tooth agenesis (21%), excluding third molars. Among the sample segment aged 14 years or greater (N = 77), occurrence of third-molar agenesis (48%) exceeded twice its normal frequency. Significant increases in occurrence of microdontia of maxillary lateral incisors (20.6%), infraocclusion of deciduous molars (24.6%), and distoangulation of mandibular second premolars (7.8%) were observed. Palatally displaced canine anomaly was also significantly elevated (8.1%). Conclusion: The results provide evidence that agenesis of other permanent teeth, microdontia, deciduous molar infraocclusion, and certain dental ectopias are the products of the same genetic mechanisms that cause second-premolar agenesis. (Angle Orthod. 2009;79:436-441.)

The B0.5IVe CoRoT target HD 49330 I. Photometric analysis from CoRoT data

Context. Be stars undergo outbursts producing a circumstellar disk from the ejected material. The beating of non-radial pulsations has been put forward as a possible mechanism of ejection. Aims. We analyze the pulsational behavior of the early B0.5IVe star HD 49330 observed during the first CoRoT long run towards the Galactical anticenter (LRA1). This Be star is located close to the lower edge of the beta Cephei instability strip in the HR diagram and showed a 0.03 mag outburst during the CoRoT observations. It is thus an ideal case for testing the aforementioned hypothesis. Methods. We analyze the CoRoT light curve of HD 49330 using Fourier methods and non-linear least square fitting. Results. In this star, we find pulsation modes typical of beta Cep stars (p modes) and SPB stars (g modes) with amplitude variations along the run directly correlated with the outburst. These results provide new clues about the origin of the Be phenomenon as well as strong constraints on the seismic modelling of Be stars.

Cyclic variability of the circumstellar disk of the Be star zeta Tauri II. Testing the 2D global disk oscillation model

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the HI emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star zeta Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for zeta Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the H alpha and Br gamma lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

Testing Evolutionary and Dispersion Scenarios for the Settlement of the New World

Background: Discussion surrounding the settlement of the New World has recently gained momentum with advances in molecular biology, archaeology and bioanthropology. Recent evidence from these diverse fields is found to support different colonization scenarios. The currently available genetic evidence suggests a ""single migration'' model, in which both early and later Native American groups derive from one expansion event into the continent. In contrast, the pronounced anatomical differences between early and late Native American populations have led others to propose more complex scenarios, involving separate colonization events of the New World and a distinct origin for these groups. Methodology/Principal Findings: Using large samples of Early American crania, we: 1) calculated the rate of morphological differentiation between Early and Late American samples under three different time divergence assumptions, and compared our findings to the predicted morphological differentiation under neutral conditions in each case; and 2) further tested three dispersal scenarios for the colonization of the New World by comparing the morphological distances among early and late Amerindians, East Asians, Australo-Melanesians and early modern humans from Asia to geographical distances associated with each dispersion model. Results indicate that the assumption of a last shared common ancestor outside the continent better explains the observed morphological differences between early and late American groups. This result is corroborated by our finding that a model comprising two Asian waves of migration coming through Bering into the Americas fits the cranial anatomical evidence best, especially when the effects of diversifying selection to climate are taken into account. Conclusions: We conclude that the morphological diversity documented through time in the New World is best accounted for by a model postulating two waves of human expansion into the continent originating in East Asia and entering through Beringia.

Beyond the Fragmentation Threshold Hypothesis: Regime Shifts in Biodiversity Across Fragmented Landscapes

Ecological systems are vulnerable to irreversible change when key system properties are pushed over thresholds, resulting in the loss of resilience and the precipitation of a regime shift. Perhaps the most important of such properties in human-modified landscapes is the total amount of remnant native vegetation. In a seminal study Andren proposed the existence of a fragmentation threshold in the total amount of remnant vegetation, below which landscape-scale connectivity is eroded and local species richness and abundance become dependent on patch size. Despite the fact that species patch-area effects have been a mainstay of conservation science there has yet to be a robust empirical evaluation of this hypothesis. Here we present and test a new conceptual model describing the mechanisms and consequences of biodiversity change in fragmented landscapes, identifying the fragmentation threshold as a first step in a positive feedback mechanism that has the capacity to impair ecological resilience, and drive a regime shift in biodiversity. The model considers that local extinction risk is defined by patch size, and immigration rates by landscape vegetation cover, and that the recovery from local species losses depends upon the landscape species pool. Using a unique dataset on the distribution of non-volant small mammals across replicate landscapes in the Atlantic forest of Brazil, we found strong evidence for our model predictions - that patch-area effects are evident only at intermediate levels of total forest cover, where landscape diversity is still high and opportunities for enhancing biodiversity through local management are greatest. Furthermore, high levels of forest loss can push native biota through an extinction filter, and result in the abrupt, landscape-wide loss of forest-specialist taxa, ecological resilience and management effectiveness. The proposed model links hitherto distinct theoretical approaches within a single framework, providing a powerful tool for analysing the potential effectiveness of management interventions.

Quantum statistical correlations in thermal field theories: Boundary effective theory

We show that the one-loop effective action at finite temperature for a scalar field with quartic interaction has the same renormalized expression as at zero temperature if written in terms of a certain classical field phi(c), and if we trade free propagators at zero temperature for their finite-temperature counterparts. The result follows if we write the partition function as an integral over field eigenstates (boundary fields) of the density matrix element in the functional Schrodinger field representation, and perform a semiclassical expansion in two steps: first, we integrate around the saddle point for fixed boundary fields, which is the classical field phi(c), a functional of the boundary fields; then, we perform a saddle-point integration over the boundary fields, whose correlations characterize the thermal properties of the system. This procedure provides a dimensionally reduced effective theory for the thermal system. We calculate the two-point correlation as an example.

Testing Lorentz Invariance and CPT Conservation with NuMI Neutrinos in the MINOS Near Detector

A search for a sidereal modulation in the MINOS near detector neutrino data was performed. If present, this signature could be a consequence of Lorentz and CPT violation as predicted by the effective field theory called the standard-model extension. No evidence for a sidereal signal in the data set was found, implying that there is no significant change in neutrino propagation that depends on the direction of the neutrino beam in a sun-centered inertial frame. Upper limits on the magnitudes of the Lorentz and CPT violating terms in the standard-model extension lie between 10(-4) and 10(-2) of the maximum expected, assuming a suppression of these signatures by a factor of 10(-17).

Phase transitions and spatially ordered counterion association in ionic-lipid membranes: A statistical model

We propose a statistical model to account for the gel-fluid anomalous phase transitions in charged bilayer- or lamellae-forming ionic lipids. The model Hamiltonian comprises effective attractive interactions to describe neutral-lipid membranes as well as the effect of electrostatic repulsions of the discrete ionic charges on the lipid headgroups. The latter can be counterion dissociated (charged) or counterion associated (neutral), while the lipid acyl chains may be in gel (low-temperature or high-lateral-pressure) or fluid (high-temperature or low-lateral-pressure) states. The system is modeled as a lattice gas with two distinct particle types-each one associated, respectively, with the polar-headgroup and the acyl-chain states-which can be mapped onto an Ashkin-Teller model with the inclusion of cubic terms. The model displays a rich thermodynamic behavior in terms of the chemical potential of counterions (related to added salt concentration) and lateral pressure. In particular, we show the existence of semidissociated thermodynamic phases related to the onset of charge order in the system. This type of order stems from spatially ordered counterion association to the lipid headgroups, in which charged and neutral lipids alternate in a checkerboard-like order. Within the mean-field approximation, we predict that the acyl-chain order-disorder transition is discontinuous, with the first-order line ending at a critical point, as in the neutral case. Moreover, the charge order gives rise to continuous transitions, with the associated second-order lines joining the aforementioned first-order line at critical end points. We explore the thermodynamic behavior of some physical quantities, like the specific heat at constant lateral pressure and the degree of ionization, associated with the fraction of charged lipid headgroups.

Statistical models of mixtures with a biaxial nematic phase

We consider a simple Maier-Saupe statistical model with the inclusion of disorder degrees of freedom to mimic the phase diagram of a mixture of rodlike and disklike molecules. A quenched distribution of shapes leads to a phase diagram with two uniaxial and a biaxial nematic structure. A thermalized distribution, however, which is more adequate to liquid mixtures, precludes the stability of this biaxial phase. We then use a two-temperature formalism, and assume a separation of relaxation times, to show that a partial degree of annealing is already sufficient to stabilize a biaxial nematic structure.

Testing for large extra dimensions with neutrino oscillations

We consider a model where sterile neutrinos can propagate in a large compactified extra dimension giving rise to Kaluza-Klein (KK) modes and the standard model left-handed neutrinos are confined to a 4-dimensional spacetime brane. The KK modes mix with the standard neutrinos modifying their oscillation pattern. We examine former and current experiments such as CHOOZ, KamLAND, and MINOS to estimate the impact of the possible presence of such KK modes on the determination of the neutrino oscillation parameters and simultaneously obtain limits on the size of the largest extra dimension. We found that the presence of the KK modes does not essentially improve the quality of the fit compared to the case of the standard oscillation. By combining the results from CHOOZ, KamLAND, and MINOS, in the limit of a vanishing lightest neutrino mass, we obtain the stronger bound on the size of the extra dimension as similar to 1.0(0.6) mu m at 99% C.L. for normal (inverted) mass hierarchy. If the lightest neutrino mass turns out to be larger, 0.2 eV, for example, we obtain the bound similar to 0.1 mu m. We also discuss the expected sensitivities on the size of the extra dimension for future experiments such as Double CHOOZ, T2K, and NO nu A.

Testing the interaction between dark energy and dark matter via the latest observations

Cosmological analyses based on currently available observations are unable to rule out a sizeable coupling between dark energy and dark matter. However, the signature of the coupling is not easy to grasp, since the coupling is degenerate with other cosmological parameters, such as the dark energy equation of state and the dark matter abundance. We discuss possible ways to break such degeneracy. Based on the perturbation formalism, we carry out the global fitting by using the latest observational data and get a tight constraint on the interaction between dark sectors. We find that the appropriate interaction can alleviate the coincidence problem.

Systematic investigation of a family of gradient-dependent functionals for solids

Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. On a test set of 60 solids, we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.