Positional bias of general and tissue-specific regulatory motifs in mouse gene promoters

Background: The arrangement of regulatory motifs in gene promoters, or promoterarchitecture, is the result of mutation and selection processes that have operated over manymillions of years. In mammals, tissue-specific transcriptional regulation is related to the presence ofspecific protein-interacting DNA motifs in gene promoters. However, little is known about therelative location and spacing of these motifs. To fill this gap, we have performed a systematic searchfor motifs that show significant bias at specific promoter locations in a large collection ofhousekeeping and tissue-specific genes.Results: We observe that promoters driving housekeeping gene expression are enriched inparticular motifs with strong positional bias, such as YY1, which are of little relevance in promotersdriving tissue-specific expression. We also identify a large number of motifs that show positionalbias in genes expressed in a highly tissue-specific manner. They include well-known tissue-specificmotifs, such as HNF1 and HNF4 motifs in liver, kidney and small intestine, or RFX motifs in testis,as well as many potentially novel regulatory motifs. Based on this analysis, we provide predictionsfor 559 tissue-specific motifs in mouse gene promoters.Conclusion: The study shows that motif positional bias is an important feature of mammalianproximal promoters and that it affects both general and tissue-specific motifs. Motif positionalconstraints define very distinct promoter architectures depending on breadth of expression andtype of tissue.

BadiRate: Estimating Family Turnover Rates by Likelihood-Based Methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Interrogating eleven fast-evolving genes for signatures of recent positive selection in worldwide human populations

Different signatures of natural selection persist over varying time scales in our genome, revealing possible episodes of adaptative evolution during human history. Here, we identify genes showing signatures of ancestral positive selection in the human lineage and investigate whether some of those genes have been evolving adaptatively in extant human populations. Specifically, we compared more than 11,000 human genes with their orthologs inchimpanzee, mouse, rat and dog and applied a branch-site likelihood method to test for positive selection on the human lineage. Among the significant cases, a robust set of 11 genes were then further explored for signatures of recent positive selection using SNP data. We genotyped 223 SNPs in 39 worldwide populations from the HGDP Diversity panel and supplemented this information with available genotypes for up to 4,814 SNPs distributed along 2 Mb centered on each gene. After exploring the allele frequency spectrum, population differentiation and the maintainance of long unbroken haplotypes, we found signals of recent adaptative phenomena in only one of the 11 candidate gene regions. However, the signal ofrecent selection in this region may come from a different, neighbouring gene (CD5) ratherthan from the candidate gene itself (VPS37C). For this set of positively-selected genes in thehuman lineage, we find no indication that these genes maintained their rapid evolutionarypace among human populations. Based on these data, it therefore appears that adaptation forhuman-specific and for population-specific traits may have involved different genes.

BadiRate: estimating family turnover rates by likelihood-based methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Antipsychotic polypharmacy in a regional health service: a population-based study

Background To analyse the extent and profile of outpatient regular dispensation of antipsychotics, both in combination and monotherapy, in the Barcelona Health Region (Spain), focusing on the use of clozapine and long-acting injections (LAI). Methods Antipsychotic drugs dispensed for people older than 18 and processed by the Catalan Health Service during 2007 were retrospectively reviewed. First and second generation antipsychotic drugs (FGA and SGA) from the Anatomical Therapeutic Chemical classification (ATC) code N05A (except lithium) were included. A patient selection algorithm was designed to identify prescriptions regularly dispensed. Variables included were age, gender, antipsychotic type, route of administration and number of packages dispensed. Results A total of 117,811 patients were given any antipsychotic, of whom 71,004 regularly received such drugs. Among the latter, 9,855 (13.9%) corresponded to an antipsychotic combination, 47,386 (66.7%) to monotherapy and 13,763 (19.4%) to unspecified combinations. Of the patients given antipsychotics in association, 58% were men. Olanzapine (37.1%) and oral risperidone (36.4%) were the most common dispensations. Analysis of the patients dispensed two antipsychotics (57.8%) revealed 198 different combinations, the most frequent being the association of FGA and SGA (62.0%). Clozapine was dispensed to 2.3% of patients. Of those who were receiving antipsychotics in combination, 6.6% were given clozapine, being clozapine plus amisulpride the most frequent association (22.8%). A total of 3.800 patients (5.4%) were given LAI antipsychotics, and 2.662 of these (70.1%) were in combination. Risperidone was the most widely used LAI. Conclusions The scant evidence available regarding the efficacy of combining different antipsychotics contrasts with the high number and variety of combinations prescribed to outpatients, as well as with the limited use of clozapine. Background To analyse the extent and profile of outpatient regular dispensation of antipsychotics, both in combination and monotherapy, in the Barcelona Health Region (Spain), focusing on the use of clozapine and long-acting injections (LAI). Methods Antipsychotic drugs dispensed for people older than 18 and processed by the Catalan Health Service during 2007 were retrospectively reviewed. First and second generation antipsychotic drugs (FGA and SGA) from the Anatomical Therapeutic Chemical classification (ATC) code N05A (except lithium) were included. A patient selection algorithm was designed to identify prescriptions regularly dispensed. Variables included were age, gender, antipsychotic type, route of administration and number of packages dispensed. Results A total of 117,811 patients were given any antipsychotic, of whom 71,004 regularly received such drugs. Among the latter, 9,855 (13.9%) corresponded to an antipsychotic combination, 47,386 (66.7%) to monotherapy and 13,763 (19.4%) to unspecified combinations. Of the patients given antipsychotics in association, 58% were men. Olanzapine (37.1%) and oral risperidone (36.4%) were the most common dispensations. Analysis of the patients dispensed two antipsychotics (57.8%) revealed 198 different combinations, the most frequent being the association of FGA and SGA (62.0%). Clozapine was dispensed to 2.3% of patients. Of those who were receiving antipsychotics in combination, 6.6% were given clozapine, being clozapine plus amisulpride the most frequent association (22.8%). A total of 3.800 patients (5.4%) were given LAI antipsychotics, and 2.662 of these (70.1%) were in combination. Risperidone was the most widely used LAI. Conclusions The scant evidence available regarding the efficacy of combining different antipsychotics contrasts with the high number and variety of combinations prescribed to outpatients, as well as with the limited use of clozapine.

Effect of growth selection on morphology in Pirenaica cattle

Twelve variables were identified to define morphology in 109 calves sired by eight Pirenaica bulls widely used in artificial insemination. The effect of selection for weight at 210 days of age (W210) was detected from the regression coefficient between the 12 variables and the selection index for W210 used in the breeding plan. Unbiased estimates of the genetic correlated responses on the morphology were obtained without previously estimating the genetic correlation among traits. In Pirenaica cattle, selection for W210 increased slaughter live weight, chest depth and corporal length. The expected changes on height (withers, back and rump), rump width and thoracic perimeter were lower but positive and statistically significant.

Review. Characterization and selection of hexaploid wheats containing resistance to Heterodera avenae or Mayetiola destructor introgressed from Aegilops

Cereal cyst nematode (CCN, Heterodera avenae) and Hessian fly (HF, Mayetiola destructor) are two major pests affecting wheat crops worldwide including important cereal areas of Spain. Aegilops ventricosa and Ae. triuncialis were used as donors in a strategy to introduce resistance genes (RG) for these pests in hexaploid wheat (Triticum aestivum L.). Two 42 chromosomes introgression lines have been derived from Ae. ventricosa: H-93-8 and H-93-33 carrying genes Cre2 and H27 conferring resistance to CCN and HF, respectively. Line TR-3531 with 42 chromosomes has been derived from Ae. triuncialis and carries RGs conferring resistance for CCN (Cre7) and for HF (H30). Alien material has been incorporated in lines H-93 by chromosomal substitution and recombination, while in line TR-3531 homoeologous recombination affecting small DNA fragments has played a major role. It has been demonstrated that Cre2, Cre7, H27 and H30 are major single dominant genes and not allelic of other previously described RGs. Biochemical and molecular-biology studies of the defense mechanism triggered by Cre2 and Cre7 have revealed specific induction of peroxidase and other antioxidant enzymes. In parallel to these basic studies advanced lines carrying resistance genes for CNN and/or HF have been developed. Selection was done using molecular markers for eventually «pyramiding» resistance genes. Several isozyme and RAPD markers have been described and, currently, new markers based on transposable elements and NBS-LRR sequences are being developed. At present, two advanced lines have already been included at the Spanish Catalogue of Commercial Plant Varieties.

Correlation of the EPR properties of perchlorotriphenylmethyl radical and their efficiency as DNP polarizers

Water soluble perchlorinated trityl (PTM) radicals were found to be effective 95 GHz DNP (dynamic nuclear polarization) polarizers in ex situ (dissolution) 13C DNP (Gabellieri et al., Angew Chem., Int. Ed. 2010, 49, 3360). The degree of the nuclear polarization obtained was reported to be dependent on the position of the chlorine substituents on the trityl skeleton. In addition, on the basis of the DNP frequency sweeps it was suggested that the 13C NMR signal enhancement is mediated by the Cl nuclei. To understand the DNP mechanism of the PTM radicals we have explored the 95 GHz EPR characteristics of these radicals that are relevant to their performance as DNP polarizers. The EPR spectra of the radicals revealed axially symmetric g-tensors. A comparison of the spectra with the 13C DNP frequency sweeps showed that although the solid effect mechanism is operational the DNP frequency sweeps reveal some extra width suggesting that contributions from EPR forbidden transitions involving 35,37Cl nuclear flips are likely. This was substantiated experimentally by ELDOR (electron-electron double resonance) detected NMR measurements, which map the EPR forbidden transitions, and ELDOR experiments that follow the depolarization of the electron spin upon irradiation of the forbidden EPR transitions. DFT (density functional theory) calculations helped to assign the observed transitions and provided the relevant spin Hamiltonian parameters. These results show that the 35,37Cl hyperfine and nuclear quadrupolar interactions cause a considerable nuclear state mixing at 95 GHz thus facilitating the polarization of the Cl nuclei upon microwave irradiation. Overlap of Cl nuclear frequencies and the 13C Larmor frequency further facilitates the polarization of the 13C nuclei by spin diffusion. Calculation of the 13C DNP frequency sweep based on the Cl nuclear polarization showed that it does lead to an increase in the width of the spectra, improving the agreement with the experimental sweeps, thus supporting the existence of a new heteronuclear assisted DNP mechanism.

Effects of P300-based BCI use on reported presence in a virtual environment

Brain-computer interfaces (BCIs) are becoming more and more popular as an input device for virtual worlds and computer games. Depending on their function, a major drawback is the mental workload associated with their use and there is significant effort and training required to effectively control them. In this paper, we present two studies assessing how mental workload of a P300-based BCI affects participants" reported sense of presence in a virtual environment (VE). In the first study, we employ a BCI exploiting the P300 event-related potential (ERP) that allows control of over 200 items in a virtual apartment. In the second study, the BCI is replaced by a gaze-based selection method coupled with wand navigation. In both studies, overall performance is measured and individual presence scores are assessed by means of a short questionnaire. The results suggest that there is no immediate benefit for visualizing events in the VE triggered by the BCI and that no learning about the layout of the virtual space takes place. In order to alleviate this, we propose that future P300-based BCIs in VR are set up so as require users to make some inference about the virtual space so that they become aware of it,which is likely to lead to higher reported presence.

Disentangling the formation of contrasting tree line physiognomies combining model selection and Bayesian parameterization for simulation models.

Alpine tree-line ecotones are characterized by marked changes at small spatial scales that may result in a variety of physiognomies. A set of alternative individual-based models was tested with data from four contrasting Pinus uncinata ecotones in the central Spanish Pyrenees to reveal the minimal subset of processes required for tree-line formation. A Bayesian approach combined with Markov chain Monte Carlo methods was employed to obtain the posterior distribution of model parameters, allowing the use of model selection procedures. The main features of real tree lines emerged only in models considering nonlinear responses in individual rates of growth or mortality with respect to the altitudinal gradient. Variation in tree-line physiognomy reflected mainly changes in the relative importance of these nonlinear responses, while other processes, such as dispersal limitation and facilitation, played a secondary role. Different nonlinear responses also determined the presence or absence of krummholz, in agreement with recent findings highlighting a different response of diffuse and abrupt or krummholz tree lines to climate change. The method presented here can be widely applied in individual-based simulation models and will turn model selection and evaluation in this type of models into a more transparent, effective, and efficient exercise.

Spontaneous Speech and Emotional Response modeling based on One-class classifier oriented to Alzheimer Disease diagnosis

The purpose of our project is to contribute to earlier diagnosis of AD and better estimates of its severity by using automatic analysis performed through new biomarkers extracted from non-invasive intelligent methods. The methods selected in this case are speech biomarkers oriented to Sponta-neous Speech and Emotional Response Analysis. Thus the main goal of the present work is feature search in Spontaneous Speech oriented to pre-clinical evaluation for the definition of test for AD diagnosis by One-class classifier. One-class classifi-cation problem differs from multi-class classifier in one essen-tial aspect. In one-class classification it is assumed that only information of one of the classes, the target class, is available. In this work we explore the problem of imbalanced datasets that is particularly crucial in applications where the goal is to maximize recognition of the minority class as in medical diag-nosis. The use of information about outlier and Fractal Dimen-sion features improves the system performance.

Regression-based techniques for statistical decision making in single-case designs

The present study evaluates the performance of four methods for estimating regression coefficients used to make statistical decisions regarding intervention effectiveness in single-case designs. Ordinary least squares estimation is compared to two correction techniques dealing with general trend and one eliminating autocorrelation whenever it is present. Type I error rates and statistical power are studied for experimental conditions defined by the presence or absence of treatment effect (change in level or in slope), general trend, and serial dependence. The results show that empirical Type I error rates do not approximate the nominal ones in presence of autocorrelation or general trend when ordinary and generalized least squares are applied. The techniques controlling trend show lower false alarm rates, but prove to be insufficiently sensitive to existing treatment effects. Consequently, the use of the statistical significance of the regression coefficients for detecting treatment effects is not recommended for short data series.

Phase mask selection in wavefront coding systems: an adaptive approach

A method for optimizing the strength of a parametric phase mask for a wavefront coding imaging system is presented. The method is based on an optimization process that minimizes a proposed merit function. The goal is to achieve modulation transfer function invariance while quantitatively maintaining nal image delity. A parametric lter that copes with the noise present in the captured images is used to obtain the nal images, and this lter is optimized. The whole process results in optimum phase mask strength and optimal parameters for the restoration lter. The results for a particular optical system are presented and tested experimentally in the labo- ratory. The experimental results show good agreement with the simulations, indicating that the procedure is useful.

Approach to the Lower Pliocene marine-continental correlation from southern Spain. The micrommamal site of Alhaurín el Grande-1 (Málaga Basin, Betic Cordillera, Spain)

A new micromammal site at Alhaurín el Grande (Málaga, southern Spain) located above early Pliocene marine deposits allows an approach to the marine-continental correlation for this age. The early Pliocene marine filling throughout the Málaga Basin is developed in three transgressive-regressive sequences (Pl-1, Pl-2, and Pl-3 units) bounded by discontinuities. At the top of the intermediate sequence Pl-2, peaty sediments have yielded fossils of Rodentia, Lagomorpha, Insectivora, and Crocodylia. The presence of Cricetus barrieri Mein & Michaux, 1970 in combination with murids, both of primitive morphology, such as Apodemus gudrunae Van de Weerd, 1976, and more advanced forms (i.e. Occitanomys brailloni Michaux, 1969 and Stephanomys donnezani cordii Ruiz Bustos, 1986), points to an early Ruscinian age (MN 14 biozone). Based on the planktonic foraminifers, the biostratigraphic data indicate that marine sediments just below the micromammal beds belong to the MPl-2 biozone of the early Zanclean. Available paleomagnetic data from the marine sediments show that the micromammal bed must be located between the normal geomagnetic subchron C3n3n (4.89-4.80 Ma) and the subchron C3n2n (4.63-4.49 Ma), limiting the age of this site to the late part of the early Zanclean.