MB-MDR: Model-Based Multifactor Dimensionality Reduction for detecting interactions in high-dimensional genomic data

L’anàlisi de l’efecte dels gens i els factors ambientals en el desenvolupament de malalties complexes és un gran repte estadístic i computacional. Entre les diverses metodologies de mineria de dades que s’han proposat per a l’anàlisi d’interaccions una de les més populars és el mètode Multifactor Dimensionality Reduction, MDR, (Ritchie i al. 2001). L’estratègia d’aquest mètode és reduir la dimensió multifactorial a u mitjançant l’agrupació dels diferents genotips en dos grups de risc: alt i baix. Tot i la seva utilitat demostrada, el mètode MDR té alguns inconvenients entre els quals l’agrupació excessiva de genotips pot fer que algunes interaccions importants no siguin detectades i que no permet ajustar per efectes principals ni per variables confusores. En aquest article il•lustrem les limitacions de l’estratègia MDR i d’altres aproximacions no paramètriques i demostrem la conveniència d’utilitzar metodologies parametriques per analitzar interaccions en estudis cas-control on es requereix l’ajust per variables confusores i per efectes principals. Proposem una nova metodologia, una versió paramètrica del mètode MDR, que anomenem Model-Based Multifactor Dimensionality Reduction (MB-MDR). La metodologia proposada té com a objectiu la identificació de genotips específics que estiguin associats a la malaltia i permet ajustar per efectes marginals i variables confusores. La nova metodologia s’il•lustra amb dades de l’Estudi Espanyol de Cancer de Bufeta.

I. Identification and characterisation of epigenetically altered tumor supressor genes by proteomic and genomic approches / II. Studies of genes regualted by MeCP2 union to its promoter regions by proteomic and genmomic approches

DNA methylation has an important impact on normal cell physiology, thus any defects in this mechanism may be related to the development of various diseases In this project we are interested in identifying epigeneticaliy modified genes, in general controlled by processes related to the DNA methylation, by means of a new strategy combining protomic and genomic analyses. First, the two Dimensional-Difference Gel Electrophoresis (2-DIGE) protein analyses of extracts obtained from HCT-116 wt and double knockout for DNMT1 and DNMT3b (DKO) cells revealed 34 proteins overexpressed in the condition of DNMTs depletion. From five genes with higher transcript lavels in DKO cells, comparing with HCT-116 wt. oniy AKR1B1, UCHLl and VIM are melhylated in HCT-116. As expected. the DNA methvlation 1s lost in DKO cells. The rneth,vl ation of VIM and UCHLl promoters in some cancer samples has already been repaired, thus further studies has been focused on AKRlBI. AKR1B1 expression due lo DNA methyiaton of promoter region seems to occur specilfically in the colon cancer cell Iines. which was confirmed in the DNA rnethylation status and expression analyses. performed on 32 different cancer cell lines (including colon, breast, lymphoma, leukemia, neuroblastoma, glioma and lung cancer cell Iines) as well as normal colon and normal lymphocytes samples. AKRIBI expression after treatments with DNA demethvlating agent (AZA) was rescued in 5 coloncancer cell lines (including genetic regulation of the candidate gene. The methylation status of the rest of the genes identified in proteomic analysis was checked by methylation specific PCR (MSP) experiment and all appeared to be unmethylated. The similar research has been done also bv means of Mecp2-null mouse model For 14 selected candidate genes the analyses of expression leveis, methylation Status and MeCP2 interaction with promoters are currently being performed.

Behavioural Flexibility: evolutionary past and its role in a changing world

A través de la historia de la vida, gran parte de los organismos han desarrollado estrategias para responder a un mundo en constante cambio. Hoy en día, las actividades humanas producen cambios ambientales a una velocidad sin precedentes, lo cual se traduce en grandes desafíos para la persistencia de biodiversidad. Esta investigación evalúa las respuesta de los animales a los cambios ambientales enfocándose en la flexibilidad del comportamiento como estrategia adaptativa. En una primera aproximación a una escala evolutiva, se otorgan evidencias del vínculo hasta ahora tenue entre la cognición e historias de vida, entregando un claro apoyo a la relación entre longevidad, vida reproductiva y el tamaño del cerebro en mamíferos. La longevidad es el centro de muchas hipótesis respecto a las ventajas de desarrollar un cerebro grande, como por ejemplo en la hipótesis del buffer cognitivo y las respuestas flexibles frente a nuevos ambientes. En un segundo nivel, se abordan factores extrínsecos e intrínsecos que podrían explicar las diferencias individuales en innovación, un componente clave en la flexibilidad del comportamiento. Por medio de una aproximación experimental, se evalúan potenciales escenarios que podrían conducir a consistentes diferencias individuales en uno de los principales factores subyacentes a la innovación (i.e. la motivación), y el potencial control endocrino sobre estos escenarios. Posteriormente, con el objetivo de evaluar la respuesta de los animales frente a los cambios ambientales actuales, se explora la respuesta de los animales frente a una de las actividades humanas mas disruptivas sobre los ecosistemas, la urbanización. Por medio de un analisis filogenetico comparativo a nivel global en aves se abordan los mecanismos implicados en la perdida de biodiversidad observada en ambientes urbanos. Los resultados entregan evidencias sobre la importancia de procesos de dispersión local junto con el papel clave de los rasgos de historia de vida, pero en un sentido diferente al clasicamente pensado. Finalmente por medio de una revisión bibliográfica se entregan evidencias teóricas y empíricas que respaldan el rol clave de la flexibilidad del comportamiento en confrontar los desafíos de una vida urbana. La integración de estos resultados muestra cómo el pasado evolutivo contribuye a hacer frente a los retos ambientales actuales, y pone de relieve posibles consecuencias ante un planeta más cambiante que nunca.

Behavioural Flexibility: evolutionary past and its role in a changing world

A través de la historia de la vida, gran parte de los organismos han desarrollado estrategias para responder a un mundo en constante cambio. Hoy en día, las actividades humanas producen cambios ambientales a una velocidad sin precedentes, lo cual se traduce en grandes desafíos para la persistencia de biodiversidad. Esta investigación evalúa las respuesta de los animales a los cambios ambientales enfocándose en la flexibilidad del comportamiento como estrategia adaptativa. En una primera aproximación a una escala evolutiva, se otorgan evidencias del vínculo hasta ahora tenue entre la cognición e historias de vida, entregando un claro apoyo a la relación entre longevidad, vida reproductiva y el tamaño del cerebro en mamíferos. La longevidad es el centro de muchas hipótesis respecto a las ventajas de desarrollar un cerebro grande, como por ejemplo en la hipótesis del buffer cognitivo y las respuestas flexibles frente a nuevos ambientes. En un segundo nivel, se abordan factores extrínsecos e intrínsecos que podrían explicar las diferencias individuales en innovación, un componente clave en la flexibilidad del comportamiento. Por medio de una aproximación experimental, se evalúan potenciales escenarios que podrían conducir a consistentes diferencias individuales en uno de los principales factores subyacentes a la innovación (i.e. la motivación), y el potencial control endocrino sobre estos escenarios. Posteriormente, con el objetivo de evaluar la respuesta de los animales frente a los cambios ambientales actuales, se explora la respuesta de los animales frente a una de las actividades humanas mas disruptivas sobre los ecosistemas, la urbanización. Por medio de un analisis filogenetico comparativo a nivel global en aves se abordan los mecanismos implicados en la perdida de biodiversidad observada en ambientes urbanos. Los resultados entregan evidencias sobre la importancia de procesos de dispersión local junto con el papel clave de los rasgos de historia de vida, pero en un sentido diferente al clasicamente pensado. Finalmente por medio de una revisión bibliográfica se entregan evidencias teóricas y empíricas que respaldan el rol clave de la flexibilidad del comportamiento en confrontar los desafíos de una vida urbana. La integración de estos resultados muestra cómo el pasado evolutivo contribuye a hacer frente a los retos ambientales actuales, y pone de relieve posibles consecuencias ante un planeta más cambiante que nunca.

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Assembling genes from predicted exons in linear time with dynamic programming

In a number of programs for gene structure prediction in higher eukaryotic genomic sequences, exon prediction is decoupled from gene assembly: a large pool of candidate exons is predicted and scored from features located in the query DNA sequence, and candidate genes are assembled from such a pool as sequences of nonoverlapping frame-compatible exons. Genes are scored as a function of the scores of the assembled exons, and the highest scoring candidate gene is assumed to be the most likely gene encoded by the query DNA sequence. Considering additive gene scoring functions, currently available algorithms to determine such a highest scoring candidate gene run in time proportional to the square of the number of predicted exons. Here, we present an algorithm whose running time grows only linearly with the size of the set of predicted exons. Polynomial algorithms rely on the fact that, while scanning the set of predicted exons, the highest scoring gene ending in a given exon can be obtained by appending the exon to the highest scoring among the highest scoring genes ending at each compatible preceding exon. The algorithm here relies on the simple fact that such highest scoring gene can be stored and updated. This requires scanning the set of predicted exons simultaneously by increasing acceptor and donor position. On the other hand, the algorithm described here does not assume an underlying gene structure model. Indeed, the definition of valid gene structures is externally defined in the so-called Gene Model. The Gene Model specifies simply which gene features are allowed immediately upstream which other gene features in valid gene structures. This allows for great flexibility in formulating the gene identification problem. In particular it allows for multiple-gene two-strand predictions and for considering gene features other than coding exons (such as promoter elements) in valid gene structures.

Identifying protein-coding genes in genomic sequences

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

A genomic analysis identifies a novel component in the genetic structure of sub-Saharan African populations

Studies of large sets of SNP data have proven to be a powerful tool in the analysis of the genetic structure of human populations. In this work, we analyze genotyping data for 2,841 SNPs in 12 Sub-Saharan African populations, including a previously unsampled region of south-eastern Africa (Mozambique). We show that robust results in a world-wide perspective can be obtained when analyzing only 1,000 SNPs. Our main results both confirm the results of previous studies, and show new and interesting features in Sub-Saharan African genetic complexity. There is a strong differentiation of Nilo-Saharans, much beyond what would be expected by geography. Hunter-gatherer populations (Khoisan and Pygmies) show a clear distinctiveness with very intrinsic Pygmy (and not only Khoisan) genetic features. Populations of the West Africa present an unexpected similarity among them, possibly the result of a population expansion. Finally, we find a strong differentiation of the south-eastern Bantu population from Mozambique, which suggests an assimilation of a pre-Bantu substrate by Bantu speakers in the region.

Design and evaluation of a panel os single-nucleotide polymorphisms in microRNA genomic regions for association studies in human disease

MicroRNAs (miRNA) are recognized posttranscriptional gene repressors involved in the control of almost every biological process. Allelic variants in these regions may be an important source of phenotypic diversity and contribute to disease susceptibility. We analyzed the genomic organization of 325 human miRNAs (release 7.1, miRBase) to construct a panel of 768 single-nucleotide polymorphisms (SNPs) covering approximately 1 Mb of genomic DNA, including 131 isolated miRNAs (40%) and 194 miRNAs arranged in 48 miRNA clusters, as well as their 5-kb flanking regions. Of these miRNAs, 37% were inside known protein-coding genes, which were significantly associated with biological functions regarding neurological, psychological or nutritional disorders. SNP coverage analysis revealed a lower SNP density in miRNAs compared with the average of the genome, with only 24 SNPs located in the 325 miRNAs studied. Further genotyping of 340 unrelated Spanish individuals showed that more than half of the SNPs in miRNAs were either rare or monomorphic, in agreement with the reported selective constraint on human miRNAs. A comparison of the minor allele frequencies between Spanish and HapMap population samples confirmed the applicability of this SNP panel to the study of complex disorders among the Spanish population, and revealed two miRNA regions, hsa-mir-26a-2 in the CTDSP2 gene and hsa-mir-128-1 in the R3HDM1 gene, showing geographical allelic frequency variation among the four HapMap populations, probably because of differences in natural selection. The designed miRNA SNP panel could help to identify still hidden links between miRNAs and human disease.

Genetic and genomic analysis modeling of germline c-MYC overexpression and cancer susceptibility

Background: Germline genetic variation is associated with the differential expression of many human genes. The phenotypic effects of this type of variation may be important when considering susceptibility to common genetic diseases. Three regions at 8q24 have recently been identified to independently confer risk of prostate cancer. Variation at 8q24 has also recently been associated with risk of breast and colorectal cancer. However, none of the risk variants map at or relatively close to known genes, with c-MYC mapping a few hundred kilobases distally. Results: This study identifies cis-regulators of germline c-MYC expression in immortalized lymphocytes of HapMap individuals. Quantitative analysis of c-MYC expression in normal prostate tissues suggests an association between overexpression and variants in Region 1 of prostate cancer risk. Somatic c-MYC overexpression correlates with prostate cancer progression and more aggressive tumor forms, which was also a pathological variable associated with Region 1. Expression profiling analysis and modeling of transcriptional regulatory networks predicts a functional association between MYC and the prostate tumor suppressor KLF6. Analysis of MYC/Myc-driven cell transformation and tumorigenesis substantiates a model in which MYC overexpression promotes transformation by down-regulating KLF6. In this model, a feedback loop through E-cadherin down-regulation causes further transactivation of c-MYC.Conclusion: This study proposes that variation at putative 8q24 cis-regulator(s) of transcription can significantly alter germline c-MYC expression levels and, thus, contribute to prostate cancer susceptibility by down-regulating the prostate tumor suppressor KLF6 gene.

Drivers and sources of supply flexibility: An exploratory study

Purpose - There has been much research on manufacturing flexibility, but supply chain flexibility is still an under-investigated area. This paper focuses on supply flexibility, the aspects of flexibility related to the upstream supply chain. Our purpose is to investigate why and how firms increase supply flexibility.Methodology/Approach An exploratory multiple case study was conducted. We analyzed seven Spanish manufacturers from different sectors (automotive, apparel, electronics and electrical equipment).Findings - The results show that there are some major reasons why firms need supply flexibility (manufacturing schedule fluctuations, JIT purchasing, manufacturing slack capacity, low level of parts commonality, demand volatility, demand seasonality and forecast accuracy), and that companies increase this type of flexibility by implementing two main strategies: to increase suppliers responsiveness capability and flexible sourcing . The results also suggest that the supply flexibility strategy selected depends on two factors: the supplier searching and switching costs and the type of uncertainty (mix, volume or delivery).Research limitations - This paper has some limitations common to all case studies, such as the subjectivity of the analysis, and the questionable generalizability of results (since the sample of firms is not statistically significant).Implications - Our study contributes to the existing literature by empirically investigating which are the main reasons for companies needing to increase supply flexibility, how they increase this flexibility, and suggesting some factors that could influence the selection of a particular supply flexibility strategy.

Flexibility vs. rigidity: Does Spain have the worst of both worlds?

In this paper we study the structure of labor market flows in Spain and compare them with France and the US. We characterize a number of empirical regularities and stylized facts. One striking result is that the job finding rate is slightly higher than in France, while the jon loss rate is much higher, putting Spain half-way between France and the US. This suggests that while Spain has borne the full cost of its labor market reforms in terms of job precarity, the benefits in terms of job creation have been quite modest. We hypothesize that this has been due to the reform s credibility being imperfect, which leads toexpectation of reversal.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

"Next-generation" Sequencing (NGS): The new genomic revolution

In the past 5 years "Next-generation" Sequencing (NGS) technologies have transformed genomics by delivering fast, inexpensive and accurate genomeinformation changing the way we think about scientific approaches in basic,applied and clinical research. The inexpensive production of large volumes ofsequence data is the main advantage over the automated Sanger method,making this new technology useful for many applications. In this chapter, a brieftechnical review of NGS technologies is given, along with the keys to NGSsuccess and a broad range of applications for NGS technologies.

Understanding the gains from wage flexibility: The exchange rate connection

We study the gains from increased wage flexibility and their dependence on exchange rate policy, using a small open economy model with staggered price andwage setting. Two results stand out: (i) the impact of wage adjustments on employment is smaller the more the central bank seeks to stabilize the exchange rate,and (ii) an increase in wage flexibility often reduces welfare, and more likely so ineconomies under an exchange rate peg or an exchange rate-focused monetary policy.Our findings call into question the common view that wage flexibility is particularlydesirable in a currency union.