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.

Endothelial Nitric Oxide Synthase Haplotypes Associated with Hypertension Do Not Predispose to Cardiac Hypertrophy

Left ventricular hypertrophy (LVH) is a complication that may result from chronic hypertension. While nitric oxide (NO) deficiency has been associated with LVH, inconsistent results have been reported with regards to the association of endothelial NO synthase (eNOS) polymorphisms and LVH in hypertensive patients. This study aims to assess whether eNOS haplotypes are associated with LVH in hypertensive patients. This study included 101 healthy controls and 173 hypertensive patients submitted to echocardiography examination. Genotypes for three eNOS polymorphisms were determined: a single-nucleotide polymorphism in the promoter region (T-786C) and in exon 7 (Glu298Asp), and variable number of tandem repeats in intron 4. We found no significant association between eNOS genotypes and hypertension or with LVH (all p>0.05). However, while we found two eNOS haplotypes associated with variable risk of hypertension (all p<0.05), we found no significant associations between eNOS haplotypes and LVH (all p>0.05), even after adjustment in multiple linear regression analysis. These findings suggest that eNOS haplotypes that have been associated with variable susceptibility to hypertension were not associated with LVH in hypertensive patients. Further studies are necessary to examine whether other genes downstream may interact with eNOS polymorphisms and predispose to LVH in hypertensive patients.

Genetic polymorphisms in TLR4, CR1 and Duffy genes are not associated with malaria resistance in patients from Baixo Amazonas region, Brazil

The main purpose of this research was to analyze the relation of the genetic polymorphisms frequently expressed by antigen-presenting cells, erythrocytes and malaria susceptibility/resistance with the human malaria infection cases. The sample used consisted of 23 Plasmodium vivax ( Pv)- and P. falciparum ( Pf)-infected patients, and 21 healthy individuals as a control group, from the Baixo Amazonas population in Para, Brazil. The Asp299Gly polymorphisms in the Toll-like receptor 4 ( TLR4), and Gly42Asp, Arg89Cys, Ala100Thr, and T-33C in the Duffy gene ( FY) were analyzed by restriction fragment length polymorphism-polymerase chain reaction. The Lys1590Glu and Arg1601Gly polymorphisms in the complement receptor type 1 (CR1) were analyzed by DNA sequencing. According to the results obtained and statistical analysis considering a significance level or alpha = 0.01, we conclude that the low heterozygote frequency (2.27%) for the Asp299Gly mutation, detected in the TLR4 gene, is not related to the Pv and Pf infections in the patients analyzed. Also, the promoter region GATA-1 analysis of the FY gene in the Pv-infected patients showed that the heterozygote frequency for the T-33C mutation (11.36% of the infected patients and 20.45% of the control patients) is not related to infection resistance. Regarding the CR1 gene, the observed heterozygote frequency (9.09%) for the Arg1601Gly mutation in Pf-infected patients when compared to heterozygote frequency in the control group (18.18%) suggests that there is no correlation with infection resistance.

Prognostic value of TP53 Pro47Ser and Arg72Pro single nucleotide polymorphisms and the susceptibility to gliomas in individuals from Southeast Brazil

The TP53 tumor suppressor gene codifies a protein responsible for preventing cells with genetic damage from growing and dividing by blocking cell growth or apoptosis pathways. A common single nucleotide polymorphism (SNP) in TP53 codon 72 (Arg72Pro) induces a 15-fold decrease of apoptosis-inducing ability and has been associated with susceptibility to human cancers. Recently, another TP53 SNP at codon 47 (Pro47Ser) was reported to have a low apoptosis-inducing ability; however, there are no association studies between this SNP and cancer. Aiming to study the role of TP53 Pro47Ser and Arg72Pro on glioma susceptibility and oncologic prognosis of patients, we investigated the genotype distribution of these SNPs in 94 gliomas (81 astrocytomas, 8 ependymomas and 5 oligodendrogliomas) and in 100 healthy subjects by the polymerase chain reaction-restriction fragment length polymorphism approach. Chi-square and Fisher exact test comparisons for genotype distributions and allele frequencies did not reveal any significant difference between patients and control groups. Overall and disease-free survivals were calculated by the Kaplan-Meier method, and the log-rank test was used for comparisons, but no significant statistical difference was observed between the two groups. Our data suggest that TP53 Pro47Ser and Arg72Pro SNPs are not involved either in susceptibility to developing gliomas or in patient survival, at least in the Brazilian population.

Endothelial Nitric Oxide Synthase Polymorphisms and Haplotypes in Amerindians

Interethnic disparities in the distribution of endothelial nitric oxide synthase (eNOS) polymorphisms may affect nitric oxide (NO)-mediated effects of and responses to drugs. While there are differences between black and white subjects there is no information regarding the distribution of eNOS gene alleles and haplotypes in Amerindians. We studied three clinically relevant eNOS polymorphisms (T(-786) C in the promoter, a variable number of tandem repeats in intron 4, and the Glu298Asp in exon 7) and eNOS haplotypes in 170 Amerindians from three tribes of the Brazilian Amazon. The results were compared with previous findings for black and white Brazilians. The Asp298, C(-786), and 4a alleles were much less common in Amerindians (5.0%, 3.2%, and 4.1%, respectively) than in blacks (15.1%, 19.5%, and 32.0%, respectively) or whites (32.8%, 41.9%, and 17.9%, respectively) (p<0.001). The haplotype including the most common alleles for each polymorphism was much more common in Amerindians (89%) than in blacks (45%) or whites (41%). Our findings are consistent with a lower genetic diversity in Amerindians compared with blacks and whites. These striking differences may be of major relevance for case-control association studies focusing on eNOS gene polymorphisms and may explain, at least in part, differences in the responses to cardiovascular drugs.

Endothelial Nitric Oxide Synthase Haplotypes Associated with Aura in Patients with Migraine

There is strong evidence implicating nitric oxide (NO) in the pathophysiology of migraine and aura. Therefore, genetic polymorphisms in the endothelial NO synthase (eNOS) gene have been studied as candidate markers for migraine susceptibility. We compared for the first time the distribution of eNOS haplotypes including the three clinically relevant eNOS polymorphisms (T(-786)C in the promoter, rs2070744; Glu298Asp in exon 7, rs1799983; and a 27 bp variable number of tandem repeats in intron 4) and two additional tagging single-nucleotide polymorphisms (rs3918226 and rs743506) in 178 women with migraine (134 without aura and 44 with aura) and 117 healthy controls (control group). Genotypes were determined by TaqMan allele discrimination assay, real-time polymerase chain reaction, and polymerase chain reaction followed by fragment separation by electrophoresis. The GA (rs743506) genotype was more common in the control group than in women with migraine (odds ratio = 0.47, 95% confidence interval [CI] 0.29-0.78, p<0.01). No significant differences were found in allele distributions for the five eNOS polymorphisms. However, the haplotypes including the variants ""C C a Glu G"" and the variants ""C C b Glu G"" were more common in women with migraine with aura than in women with migraine without aura (odds ratio = 30.71, 95% CI = 1.61-586.4 and odds ratio = 17.26, 95% CI = 1.94-153.4, respectively; both p<0.0015625). These findings suggest that these two eNOS haplotypes affect the susceptibility to the presence of aura in patients with migraine.

Differential expression of E-cadherin gene in human neuroepithelial tumors

Cadherins are cell-to-cell adhesion molecules that play an important role in the establishment of adherent-type junctions by mediating calcium-dependent cellular interactions. The CDH1 gene encodes the transmembrane glycoprotein E-cadherin which is important in maintaining homophilic cell-cell adhesion in epithelial tissues. E-cadherin interacts with catenin proteins to maintain tissue architecture. Structural defects or loss of expression of E-cadherin have been reported as a common feature in several human cancer types. This study aimed to evaluate the expression of E-cadherin and their correlation with clinical features in microdissected brain tumor samples from 81 patients, divided into 62 astrocytic tumors grades I to IV and 19 medulloblastomas, and from 5 white matter non-neoplasic brain tissue samples. E-cadherin (CDH1) gene expression was analyzed by quantitative real-time polymerase chain reaction. Mann-Whitney, Kruskal-Wallis, Kaplan-Meir, and log-rank tests were performed for statistical analyses. We observed a decrease in expression among pathological grades of neuroepithelial tumors. Non-neoplasic brain tissue showed a higher expression level of CDH1 gene than did neuroepithelial tumors. Expression of E-cadherin gene was higher in astrocytic than embryonal tumors (P = 0.0168). Low-grade malignancy astrocytomas (grades I-II) showed higher CDH1 expression than did high-grade malignancy astrocytomas (grades III-IV) and medulloblastomas (P < 0.0001). Non-neoplasic brain tissue showed a higher expression level of CDH1 gene than grade I malignancy astrocytomas, considered as benign tumors (P = 0.0473). These results suggest that a decrease in E-cadherin gene expression level in high-grade neuroepithelial tumors may be a hallmark of malignancy in dedifferentiated tumors and that it may be possibly correlated with their progression and dissemination.

Identification of three distinguishable phenotypes in golden retriever muscular dystrophy

Duchenne muscular dystrophy (DMD) is a human disease characterized by progressive and irreversible skeletal muscle degeneration caused by mutations in genes coding for important muscle proteins. Unfortunately, there is no efficient treatment for this disease; it causes progressive loss of motor and muscular ability until death. The canine model (golden retriever muscular dystrophy) is similar to DMD, showing similar clinical signs. Fifteen dogs were followed from birth and closely observed for clinical signs. Dogs had their disease status confirmed by polymerase chain reaction analysis and genotyping. Clinical observations of musculoskeletal, morphological, gastrointestinal, respiratory, cardiovascular, and renal features allowed us to identify three distinguishable phenotypes in dystrophic dogs: mild (grade I), moderate (grade II) and severe (grade III). These three groups showed no difference in dystrophic alterations of muscle morphology and creatine kinase levels. This information will be useful for therapeutic trials, because DMD also shows significant, inter- and intra-familiar clinical variability. Additionally, being aware of phenotypic differences in this animal model is essential for correct interpretation and understanding of results obtained in pre-clinical trials.

Estimates of genetic trend for carcass traits in a commercial broiler line

Data from the slaughter of 24,001 chickens that were part of a selection program for the production of commercial broilers were used to estimate genetic trend for absolute carcass (CW), breast meat (BRW), and leg (LW) weights, and relative carcass (CY), breast meat (BRY), and leg (LY) weights. The components of (co) variance and breeding values of individuals were obtained by the restricted maximum likelihood method applied to animal models. The relationship matrix was composed of 132,442 birds. The models included as random effects, maternal additive genetic and permanent environmental for CW, BRW, LW, CY, and BRY, and only maternal permanent environmental for LY, besides the direct additive genetic and residual effects, and as fixed effects, hatch week, parents' mating group and sex. The estimates of genetic trend were obtained by average regression of breeding value on generation, and the average genetic trend was estimated by regression coefficients. The genetic trends for CW (+ 6.0336 g/generation), BRW (+ 3.6723 g/generation), LW (+ 1.5846 g/generation), CY (+ 0.1195%/generation), and BRY (+ 0.1388%/generation) were positive, and they were in accordance with the objectives of the selection program for these traits. The genetic trend for LY(-0.0019%/generation) was negative, possibly due to the strong emphasis on selection for BRY and the negative correlations between these two traits.

Newtonian perturbations on models with matter creation

Creation of cold dark matter (CCDM) can macroscopically be described by a negative pressure, and, therefore, the mechanism is capable to accelerate the Universe, without the need of an additional dark energy component. In this framework, we discuss the evolution of perturbations by considering a Neo-Newtonian approach where, unlike in the standard Newtonian cosmology, the fluid pressure is taken into account even in the homogeneous and isotropic background equations (Lima, Zanchin, and Brandenberger, MNRAS 291, L1, 1997). The evolution of the density contrast is calculated in the linear approximation and compared to the one predicted by the Lambda CDM model. The difference between the CCDM and Lambda CDM predictions at the perturbative level is quantified by using three different statistical methods, namely: a simple chi(2)-analysis in the relevant space parameter, a Bayesian statistical inference, and, finally, a Kolmogorov-Smirnov test. We find that under certain circumstances, the CCDM scenario analyzed here predicts an overall dynamics (including Hubble flow and matter fluctuation field) which fully recovers that of the traditional cosmic concordance model. Our basic conclusion is that such a reduction of the dark sector provides a viable alternative description to the accelerating Lambda CDM cosmology.

Confronting dark energy models using galaxy cluster number counts

The mass function of cluster-size halos and their redshift distribution are computed for 12 distinct accelerating cosmological scenarios and confronted to the predictions of the conventional flat Lambda CDM model. The comparison with Lambda CDM is performed by a two-step process. First, we determine the free parameters of all models through a joint analysis involving the latest cosmological data, using supernovae type Ia, the cosmic microwave background shift parameter, and baryon acoustic oscillations. Apart from a braneworld inspired cosmology, it is found that the derived Hubble relation of the remaining models reproduces the Lambda CDM results approximately with the same degree of statistical confidence. Second, in order to attempt to distinguish the different dark energy models from the expectations of Lambda CDM, we analyze the predicted cluster-size halo redshift distribution on the basis of two future cluster surveys: (i) an X-ray survey based on the eROSITA satellite, and (ii) a Sunayev-Zeldovich survey based on the South Pole Telescope. As a result, we find that the predictions of 8 out of 12 dark energy models can be clearly distinguished from the Lambda CDM cosmology, while the predictions of 4 models are statistically equivalent to those of the Lambda CDM model, as far as the expected cluster mass function and redshift distribution are concerned. The present analysis suggests that such a technique appears to be very competitive to independent tests probing the late time evolution of the Universe and the associated dark energy effects.

Chemical evolution models for spiral disks: the Milky Way, M31, and M33

Context. The distribution of chemical abundances and their variation with time are important tools for understanding the chemical evolution of galaxies. In particular, the study of chemical evolution models can improve our understanding of the basic assumptions made when modelling our Galaxy and other spirals. Aims. We test a standard chemical evolution model for spiral disks in the Local Universe and study the influence of a threshold gas density and different efficiencies in the star formation rate (SFR) law on radial gradients of abundance, gas, and SFR. The model is then applied to specific galaxies. Methods. We adopt a one-infall chemical evolution model where the Galactic disk forms inside-out by means of infall of gas, and we test different thresholds and efficiencies in the SFR. The model is scaled to the disk properties of three Local Group galaxies (the Milky Way, M31 and M33) by varying its dependence on the star formation efficiency and the timescale for the infall of gas onto the disk. Results. Using this simple model, we are able to reproduce most of the observed constraints available in the literature for the studied galaxies. The radial oxygen abundance gradients and their time evolution are studied in detail. The present day abundance gradients are more sensitive to the threshold than to other parameters, while their temporal evolutions are more dependent on the chosen SFR efficiency. A variable efficiency along the galaxy radius can reproduce the present day gas distribution in the disk of spirals with prominent arms. The steepness in the distribution of stellar surface density differs from massive to lower mass disks, owing to the different star formation histories. Conclusions. The most massive disks seem to have evolved faster (i.e., with more efficient star formation) than the less massive ones, thus suggesting a downsizing in star formation for spirals. The threshold and the efficiency of star formation play a very important role in the chemical evolution of spiral disks. For instance, an efficiency varying with radius can be used to regulate the star formation. The oxygen abundance gradient can steepen or flatten in time depending on the choice of this parameter.

Automated supervised classification of variable stars in the CoRoT programme Method and application to the first four exoplanet fields

Aims. In this work, we describe the pipeline for the fast supervised classification of light curves observed by the CoRoT exoplanet CCDs. We present the classification results obtained for the first four measured fields, which represent a one-year in-orbit operation. Methods. The basis of the adopted supervised classification methodology has been described in detail in a previous paper, as is its application to the OGLE database. Here, we present the modifications of the algorithms and of the training set to optimize the performance when applied to the CoRoT data. Results. Classification results are presented for the observed fields IRa01, SRc01, LRc01, and LRa01 of the CoRoT mission. Statistics on the number of variables and the number of objects per class are given and typical light curves of high-probability candidates are shown. We also report on new stellar variability types discovered in the CoRoT data. The full classification results are publicly available.

Spatial patterns in the brood combs of Nannotrigona testaceicornis (Hymenoptera: Meliponinae): male clusters

Genetic models of sex and caste determination in eusocial stingless bees suggest specific patterns of male, worker and gyne cell distribution in the brood comb. Conflict between queen and laying workers over male parentage and center-periphery gradients of conditions, such as food and temperature, could also contribute to non-random spatial configuration. We converted the positions of the hexagonal cells in a brood comb to Cartesian coordinates, labeled by sex or caste of the individuals inside. To detect and locate clustered patterns, the mapped brood combs were evaluated by indexes of dispersion (MMC, mean distance of cells of a given category from their centroid) and eccentricity (DMB, distance between this centroid and the overall brood comb centroid) that we developed. After randomizing the labels and recalculating the indexes, we calculated probabilities that the original values had been generated by chance. We created sets of binary brood combs in which males were aggregated, regularly or randomly distributed among females. These stylized maps were used to describe the power of MMC and DMB, and they were applied to evaluate the male distribution in the sampled Nannotrigona testaceicornis brood combs. MMC was very sensitive to slight deviations from a perfectly rounded clump; DMB detected any asymmetry in the location of these compact to fuzzy clusters. Six of the 82 brood combs of N. testaceicornis that we analyzed had more than nine males, distributed according to variations in spatial patterns, as indicated by the two indexes.

Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence

Background: Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei. Results: The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infraordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes. Conclusions: Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.