Gene-gene interaction and gene-pathway analysis of genome-wide association for schizophrenia

Schizophrenia (SZ) is a complex disorder with high heritability and variable phenotypes that has limited success in finding causal genes associated with the disease development. Pathway-based analysis is an effective approach in investigating the molecular mechanism of susceptible genes associated with complex diseases. The etiology of complex diseases could be a network of genetic factors and within the genes, interaction may occur. In this work we argue that some genes might be of small effect that by itself are neither sufficient nor necessary to cause the disease however, their effect may induce slight changes to the gene expression or affect the protein function, therefore, analyzing the gene-gene interaction mechanism within the disease pathway would play crucial role in dissecting the genetic architecture of complex diseases, making the pathway-based analysis a complementary approach to GWAS technique. ^ In this study, we implemented three novel linkage disequilibrium based statistics, the linear combination, the quadratic, and the decorrelation test statistics, to investigate the interaction between linked and unlinked genes in two independent case-control GWAS datasets for SZ including participants of European (EA) and African (AA) ancestries. The EA population included 1,173 cases and 1,378 controls with 729,454 genotyped SNPs, while the AA population included 219 cases and 288 controls with 845,814 genotyped SNPs. We identified 17,186 interacting gene-sets at significant level in EA dataset, and 12,691 gene-sets in AA dataset using the gene-gene interaction method. We also identified 18,846 genes in EA dataset and 19,431 genes in AA dataset that were in the disease pathways. However, few genes were reported of significant association to SZ. ^ Our research determined the pathways characteristics for schizophrenia through the gene-gene interaction and gene-pathway based approaches. Our findings suggest insightful inferences of our methods in studying the molecular mechanisms of common complex diseases.^

FANCM AND FAAP24 MAINTAIN GENOMIC STABILITY THROUGH COOPERATIVE AND UNIQUE FUNCTIONS

Fanconi anemia (FA) is a rare recessive genetic disease with an array of clinical manifestations including multiple congenital abnormalities, progressive bone marrow failure and profound cancer susceptibility. A hallmark of cells derived from FA patients is hypersensitivity to DNA interstrand crosslinking agents such as mitomycin C (MMC) and cisplatin, suggesting that FA- and FA-associated proteins play important roles in protecting cells from DNA interstrand crosslink (ICL) damage. Two genes involved in the FA pathway, FANCM and FAAP24, are of particular interest because they contain DNA interacting domains. However, there are no definitive patient mutations for these two genes, and the resulting lack of human genetic model system renders their functional studies difficult. In this study, I established isogenic human FANCM- and FAAP24-null mutants through homologous replacement-mediated gene targeting in HCT-116 cells, and systematically investigated the functions of FANCM and FAAP24 inchromosome stability, FA pathway activation, DNA damage checkpoint signaling, and ICL repair. I found that the FANCM-/-/FAAP24-/- double mutant was much more sensitive to DNA crosslinking agents than FANCM-/- and FAAP24-/- single mutants, suggesting that FANCM and FAAP24 possess epistatic as well as unique functions in response to ICL damage. I demonstrated that FANCM and FAAP24 coordinately support the activation of FA pathway by promoting chromatin localization of FA core complex and FANCD2 monoubiqutination. They also cooperatively function to suppress sister chromatid exchange and radial chromosome formation, likely by limiting crossovers in recombination repair. In addition, I defined novel non-overlapping functions of FANCM and FAAP24 in response to ICL damage. FAAP24 plays a major role in activating ICL-induced ATR-dependent checkpoint, which is independent of its interaction with FANCM. On the other hand, FANCM promotes recombination-independent ICL repair independently of FAAP24. Mechanistically, FANCM facilitates recruitment of nucleotide excision repair machinery and lesion bypass factors to ICL damage sites through its translocase activity. Collectively, my studies provide mechanistic insights into how genome integrity is both coordinately and independently protected by FANCM and FAAP24.

INVESTIGATING THE ROLES OF THE p63 ISOFORMS IN THE MICRORNA BIOGENESIS PATHWAY

MicroRNAs play roles in various biological processes like development, tumorigenesis, metastasis and pluripotency. My thesis work has demonstrated roles for p63, a p53 family member, in the upstream regulation of microRNA biogenesis. The p63 gene has a complex gene structure and has multiple isoforms. The TAp63 isoforms contain an acidic transcription activation domain. The ΔNp63 isoforms, lack the TA domain, but have a proline rich region critical for gene transactivation. To understand the functions of these isoforms, the Flores lab generated TAp63 and ΔNp63 conditional knock out mice. Using these mice and tissues and cells from these mice we have found that TAp63 transcriptionally regulates Dicer while ΔNp63 transcriptionally regulates DGCR8. TAp63 -/- mice are highly tumor prone. These mice develop metastatic mammary adenocarcinomas, squamous cell carcinomas, and lung adenocarcinomas to distant sites including the liver, lungs, and brain. I found that TAp63 suppresses metastasis by transcriptionally activating Dicer. TAp63 and Dicer levels were very low or lost in high grade human tumors like mammary adenocarcinomas, squamous cell carcinomas, and lung adenocarcinomas. Expression of Dicer in these tumor cell lines reduced their invasiveness. Using ΔNp63 -/- mice, I found that ΔNp63 transcriptionally activates DGCR8, resulting in a miRNA profile that is critical to reprogram cells to pluripotency. Analysis of epidermal cells derived from ΔNp63 -/- mice revealed that these cells expressed markers of pluripotency, including Sox2, Oct 4 and Nanog; however, genome-wide analysis revealed a novel profile of genes that are common between ΔNp63 -/- epidermal cells and embryonic stem cells. I also found that mouse cells depleted of ΔNp63 form chimeric mice and teratomas in SCID mice, demonstrating that ΔNp63 deficient cells are pluripotent. Further, I found that restoration of DGCR8 in ΔNp63 -/- epidermal cells reduces their pluripotency and induces terminal differentiation. I also demonstrated that iMS (induced multipotent stem) cells could be generated using human keratinocytes by knockdown of ∆Np63 or DGCR8. Taken together, my work has placed p63 and its isoforms at a critical node in controlling miRNA biogenesis.

Signal transduction of {\it HER-2\/}/{\it neu\/} receptor tyrosine kinase

Overexpression and/or amplification of HER2/neu is frequently detected in many human cancers. Activation of p185 tyrosine kinase can be achieved by point mutation, overexpression, deletion, and heterodimerization with other class I receptors. In this study I investigated the signal transduction pathways mediating the oncogenic signal of the point mutation-activated rat p185. I demonstrated that tyrosine phosphorylation of Shc and formation of Shc/Grb2 complex correlated to the transformation of NIH3T3 cells caused by the point mutation-activated rat HER2/neu. Furthermore, I observed that association with Shc was severely impaired by deletion of most of the major autophosphorylation sites of the point-mutated p185. The truncated p185 product, however, fully retained its ability to transform NIH3T3 cells, induce Shc tyrosine phosphorylation and Shc/Grb2 complex formation. These results suggest that tyrosine phosphorylation of Shc which allows formation of Shc/Grb2 complex may play an important role in cell transformation induced by the point mutation-activated p185, and that stable binding to mutant p185 may not be necessary for Shc to mediate this signaling pathway.^ Recent studies have suggested that formation of the complex containing Sos, Grb2 and Shc is important in coupling receptor tyrosine kinases to the Ras signaling pathway. To clarify the role of this trimer in the oncogenic signaling of the activated p185, I set out to interfere with the protein-protein interactions in Shc/Grb2/Sos complex by introducing Grb2 mutants with deletions in either amino- ($\Delta$N-Grb2) or carboxyl- ($\Delta$C-Grb2) terminal SH3 domains into B104-1-1 cells derived from NIH3T3 cells that express the point mutation-activated HER-2/neu. I found that the transformed phenotypes of the B104-1-1 cells were largely reversed by expression of the $\Delta$N-Grb2. The effect of the $\Delta$C-Grb2 on phenotypic reversion was much weaker. Biochemical analysis showed that the $\Delta$N-Grb2 was able to associate Shc but not the activated p185 nor Sos, while the $\Delta$C-Grb2 bound to Shc, the activated p185, and Sos. The p185-mediated Ras activation was severely inhibited by the $\Delta$N-Grb2 but not the $\Delta$C-Grb2. Taken together, these data demonstrate that interruption of the interaction between Shc and the endogenous Grb2 by the $\Delta$N-Grb2 is able to impair the oncogenic signaling of the mutation-activated p185, indicating that (i) the $\Delta$N-Grb2 functions as a strong dominant-negative mutant, (ii) Shc/Grb2/Sos pathway plays a major role in mediating the oncogenic signal of the mutation-activated p185. Unlike the $\Delta$N-Grb2, the $\Delta$C-Grb2 appears to be a relatively weak dominant-negative mutant, probably due to its ability to largely fulfill the biological functions of the wild-type Grb2. ^

Evaluation of the Hosmer -Lemeshow global goodness -of-fit statistic for mixed variable logistic regression models

The performance of the Hosmer-Lemeshow global goodness-of-fit statistic for logistic regression models was explored in a wide variety of conditions not previously fully investigated. Computer simulations, each consisting of 500 regression models, were run to assess the statistic in 23 different situations. The items which varied among the situations included the number of observations used in each regression, the number of covariates, the degree of dependence among the covariates, the combinations of continuous and discrete variables, and the generation of the values of the dependent variable for model fit or lack of fit.^ The study found that the $\rm\ C$g* statistic was adequate in tests of significance for most situations. However, when testing data which deviate from a logistic model, the statistic has low power to detect such deviation. Although grouping of the estimated probabilities into quantiles from 8 to 30 was studied, the deciles of risk approach was generally sufficient. Subdividing the estimated probabilities into more than 10 quantiles when there are many covariates in the model is not necessary, despite theoretical reasons which suggest otherwise. Because it does not follow a X$\sp2$ distribution, the statistic is not recommended for use in models containing only categorical variables with a limited number of covariate patterns.^ The statistic performed adequately when there were at least 10 observations per quantile. Large numbers of observations per quantile did not lead to incorrect conclusions that the model did not fit the data when it actually did. However, the statistic failed to detect lack of fit when it existed and should be supplemented with further tests for the influence of individual observations. Careful examination of the parameter estimates is also essential since the statistic did not perform as desired when there was moderate to severe collinearity among covariates.^ Two methods studied for handling tied values of the estimated probabilities made only a slight difference in conclusions about model fit. Neither method split observations with identical probabilities into different quantiles. Approaches which create equal size groups by separating ties should be avoided. ^

Molecular mechanisms involved in the regulation of cell polarity and genome stability by the p21 -activated kinase, Shk1, in the fission yeast, Schizosaccharomyces pombe

The essential p21-activated kinase (PAK), Shk1, is a critical component of a Ras/Cdc42/PAK complex required for cell viability, normal cell polarity, proper regulation of cytoskeletal dynamics, and sexual differentiation in the fission yeast, Schizosaccharomyces pombe. While cellular functions of PAKs have been described in eukaryotes from yeasts to mammals, the molecular mechanisms of PAK regulation and function are poorly understood. This study has characterized a novel Shk1 inhibitor, Skb15, and, in addition, identified the cell polarity regulator, Tea1, as a potential biological substrate of Shk1 in S. pombe. Skb15 is a highly conserved WD repeat protein that was discovered from a two-hybrid screen for proteins that interact with the catalytic domain of Shk1. Molecular data indicate that Skb15 negatively regulates Shk1 kinase activity in S. pombe cells. A null mutation in the skb15 gene is lethal and results in deregulation of actin polymerization and localization, microtubule biogenesis, and the cytokinetic machinery, as well as a substantial uncoupling of these processes from the cell cycle. Loss of Skb15 function is suppressed by partial loss of Shk1, demonstrating that negative regulation of Shk1 by Skb15 is required for proper execution of cytoskeletal remodeling and cytokinetic functions. A mouse homolog of Skb15 can substitute for its counterpart in fission yeast, demonstrating that Skb15 protein function has been substantially conserved through evolution. ^ Our laboratory has recently demonstrated that Shk1, in addition to regulating actin cytoskeletal organization, is required for proper regulation of microtubule dynamics in S. pombe cells. The Shk1 protein localizes to interphase and mitotic microtubules, the septum-forming region, and cell ends. This pattern of localization overlaps with that of the cell polarity regulator, Tea1, in S. pombe cells. The tea1 gene was identified by Paul Nurse's laboratory from a screen for genes involved in the control of cell morphogenesis in S. pombe. In contrast to wild type S. pombe cells, which are rod shaped, tea1 null cells are often bent and/or branched in shape. The Tea1 protein localizes to the cell ends, like Shk1, and the growing tips of interphase microtubules. Thus, experiments were performed to investigate whether Tea1 interacts with Shk1. The tea1 null mutation strongly suppresses the loss of function of Skb15, an essential inhibitor of Shk1 function. All defects associated with the skb15 mutation, including defects in F-actin organization, septation, spindle elongation, and chromosome segregation, are suppressed by tea1Δ, suggesting that Tea1 may function in these diverse processes. Consistent with a role for Tea1 in cytokinesis, tea1Δ cells have a modest cell separation defect that is greatly exacerbated by a shk1 mutation and, like Shk1, Tea1 localizes to the septation site. Molecular analyses showed that Tea1 phosphorylation is significantly dependent on Shk1 function in vivo and that bacterially expressed Tea1 protein is directly phosphorylated by recombinant Shk1 kinase in vitro. Taken together, these results identify Tea1 as a potential biological substrate of Shk1 in S. pombe. ^ In summary, this study provides new insights into a conserved regulatory mechanism for PAKs, and also begins to uncover the molecular mechanisms by which the Ras/Cdc42/PAK complex regulates the microtubule and actin cytoskeletons and cell growth polarization in fission yeast. ^

Colony-specific investigations reveal highly variable responses among individual corals to ocean acidification and warming

As anthropogenic climate change is an ongoing concern, scientific investigations on its impacts on coral reefs are increasing. Although impacts of combined ocean acidification (OA) and temperature stress (T) on reef-building scleractinian corals have been studied at the genus, species and population levels, there are little data available on how individual corals respond to combined OA and anomalous temperatures. In this study, we exposed individual colonies of Acropora digitifera, Montipora digitata and Porites cylindrica to four pCO2-temperature treatments including 400 µatm-28 °C, 400 µatm-31 °C, 1000 µatm-28 °C and 1000 µatm-31 °C for 26 days. Physiological parameters including calcification, protein content, maximum photosynthetic efficiency, Symbiodinium density, and chlorophyll content along with Symbiodinium type of each colony were examined. Along with intercolonial responses, responses of individual colonies versus pooled data to the treatments were investigated. The main results were: 1) responses to either OA or T or their combination were different between individual colonies when considering physiological functions; 2) tolerance to either OA or T was not synonymous with tolerance to the other parameter; 3) tolerance to both OA and T did not necessarily lead to tolerance of OA and T combined (OAT) at the same time; 4) OAT had negative, positive or no impacts on physiological functions of coral colonies; and 5) pooled data were not representative of responses of all individual colonies. Indeed, the pooled data obscured actual responses of individual colonies or presented a response that was not observed in any individual. From the results of this study we recommend improving experimental designs of studies investigating physiological responses of corals to climate change by complementing them with colony-specific examinations.

(Table A1) Species list of the high Antarctic Weddell Sea food web

Human-induced habitat destruction, overexploitation, introduction of alien species and climate change are causing species to go extinct at unprecedented rates, from local to global scales. There are growing concerns that these kinds of disturbances alter important functions of ecosystems. Our current understanding is that key parameters of a community (e.g. its functional diversity, species composition, and presence/absence of vulnerable species) reflect an ecological network's ability to resist or rebound from change in response to pressures and disturbances, such as species loss. If the food web structure is relatively simple, we can analyse the roles of different species interactions in determining how environmental impacts translate into species loss. However, when ecosystems harbour species-rich communities, as is the case in most natural systems, then the complex network of ecological interactions makes it a far more challenging task to perceive how species' functional roles influence the consequences of species loss. One approach to deal with such complexity is to focus on the functional traits of species in order to identify their respective roles: for instance, large species seem to be more susceptible to extinction than smaller species. Here, we introduce and analyse the marine food web from the high Antarctic Weddell Sea Shelf to illustrate the role of species traits in relation to network robustness of this complex food web. Our approach was threefold: firstly, we applied a new classification system to all species, grouping them by traits other than body size; secondly, we tested the relationship between body size and food web parameters within and across these groups and finally, we calculated food web robustness. We addressed questions regarding (i) patterns of species functional/trophic roles, (ii) relationships between species functional roles and body size and (iii) the role of species body size in terms of network robustness. Our results show that when analyzing relationships between trophic structure, body size and network structure, the diversity of predatory species types needs to be considered in future studies.

Stable isotope record of planktonic foraminifera from the western Iberian margin during marine isotope stage 2

The surface water hydrography along the western Iberian margin, as part of the North Atlantic's eastern boundary upwelling system, consists of a complex, seasonally variable system of equatorward and poleward surface and subsurface currents and seasonal upwelling. Not much information exists to ascertain if the modern current and productivity patterns subsisted under glacial climate conditions, such as during marine isotope stage (MIS) 2, and how North Atlantic meltwater events, especially Heinrich events, affected them. To help answer these questions we are combining stable isotope records of surface to subsurface dwelling planktonic foraminifer species with sea surface temperature and export productivity data for four cores distributed along the western and southwestern Iberian margin (MD95-2040, MD95-2041, MD99-2336, and MD99-2339). The records reveals that with the exception of the Heinrich events and Greenland Stadial (GS) 4 hydrographic conditions along the western Iberian margin were not much different from the present. During the Last Glacial Maximum (LGM), subtropical surface and subsurface waters penetrated poleward to at least 40.6°N (site MD95-2040). Export productivity was, in general, high on the western margin during the LGM and low in the central Gulf of Cadiz, in agreement with the modern situation. During the Heinrich events and GS 4, on the other hand, productivity was high in the Gulf of Cadiz and suppressed in the upwelling regions along the western margin where a strong halocline inhibited upwelling. Heinrich event 1 had the strongest impact on the hydrography and productivity off Iberia and was the only period when subarctic surface waters were recorded in the central Gulf of Cadiz. South of Lisbon (39°N), the impact of the other Heinrich events was diminished, and not all of them led to a significant cooling in the surface waters. Thus, climatic impacts of Heinrich events highly varied with latitude and the prevailing hydrographic conditions in this region.

Implicaciones de la accesibilidad configuracional en la movilidad peatonal. El caso de Madrid = The implications of configurational accessibility in pedestrian mobility. The case of Madrid

- Resumen La hipótesis que anima esta tesis doctoral es que algunas de las características del entorno urbano, en particular las que describen la accesibilidad de su red de espacio público, podrían estar relacionadas con la proporción de viajes a pie o reparto modal, que tiene cada zona o barrio de Madrid. Uno de los puntos de partida de dicha hipótesis que el entorno urbano tiene una mayor influencia sobre los viaje a pie que en sobre otros modos de transporte, por ejemplo que en los viajes de bicicleta o en transporte público; y es que parece razonable suponer que estos últimos van a estar más condicionadas por ejemplo por la disponibilidad de vías ciclistas, en el primer caso, o por la existencia de un servicio fiable y de calidad, en el segundo. Otra de las motivaciones del trabajo es que la investigación en este campo de la accesibilidad del espacio público, en concreto la denominada “Space Syntax”, ha probado en repetidas ocasiones la influencia de la red de espacio público en cómo se distribuye la intensidad del tráfico peatonal por la trama urbana, pero no se han encontrado referencias de la influencia de dicho elemento sobre el reparto modal. De acuerdo con la hipótesis y con otros trabajos anteriores se propone una metodología basada en el análisis empírico y cuantitativo. Su objetivo es comprobar si la red de espacio público, independientemente de otras variables como los usos del suelo, incluso de las variables de ajenas entorno no construido, como las socioeconómicas, está o no relacionada estadísticamente con la proporción de peatones viajes en las zonas urbanas. Las técnicas estadísticas se utilizan para comprobar sistemáticamente la asociación de las variables del entorno urbano, denominadas variables independientes, con el porcentaje de viajes a pie, la variable dependiente. En términos generales, la metodología es similar a la usada en otros trabajos en este campo como los de CERVERÓ y KOCKLEMAN (1997), CERVERÓ y DUNCAN (2003), o para los que se utilizan principalmente en la revisión general de TRB (2005) o, más recientemente, en ZEGRAS (2006) o CHATMAN (2009). Otras opciones metodológicas, como los métodos de preferencias declaradas (ver LOUVIERE, HENSHER y SWAIT, 2000) o el análisis basado en agentes (PENN & TURNER, 2004) fueron descartados, debido a una serie de razones, demasiado extensas para ser descritas aquí. El caso de estudio utilizado es la zona metropolitana de Madrid, abarcándola hasta la M-50, es decir en su mayor parte, con un tamaño aproximado de 31x34 Km y una población de 4.132.820 habitantes (aproximadamente el 80% de la población de la región). Las principales fuentes de datos son la Encuesta Domiciliaria de Movilidad de 2004 (EDM04), del Consorcio Regional de Transportes de Madrid que es la última disponible (muestra: > 35.000 familias,> 95.000 personas), y un modelo espacial del área metropolitana, integrando el modelo para calcular los índices de Space Syntax y un Sistema de Información Geográfica (SIG). La unidad de análisis, en este caso las unidades espaciales, son las zonas de transporte (con una población media de 7.063 personas) y los barrios (con una población media de 26.466 personas). Las variables del entorno urbano son claramente el centro del estudio. Un total de 20 índices (de 21) se seleccionan de entre los más relevantes encontrados en la revisión de la producción científica en este campo siendo que, al mismo tiempo, fueran accesibles. Nueve de ellos se utilizan para describir las características de los usos del suelo, mientras que otros once se usan para describir la red de espacios públicos. Estos últimos incluyen las variables de accesibilidad configuracional, que son, como se desprende de su título, el centro del estudio propuesto. La accesibilidad configuracional es un tipo especial de accesibilidad que se basa en la configuración de la trama urbana, según esta fue definida por HILLIER (1996), el autor de referencia dentro de esta línea de investigación de Space Syntax. Además se incluyen otras variables de la red de espacio público más habituales en los estudios de movilidad, y que aquí se denominan características geométricas de los elementos de la red, tales como su longitud, tipo de intersección, conectividad, etc. Por último se incluye además una variable socioeconómica, es decir ajena al entorno urbano, para evaluar la influencia de los factores externos, pues son varios los que pueden tener un impacto en la decisión de caminar (edad, género, nivel de estudios, ingresos, tasa de motorización, etc.). La asociación entre las variables se han establecido usando análisis de correlación (bivariante) y modelos de análisis multivariante. Las primeras se calculan entre por pares entre cada una de las 21 variables independientes y la dependiente, el porcentaje de viajes a pie. En cuanto a los segundos, se han realizado tres tipos de estudios: modelo multivariante general lineal, modelo multivariante general curvilíneo y análisis discriminante. Todos ellos son capaces de generar modelos de asociación entre diversas variables, pudiéndose de esta manera evaluar con bastante precisión en qué medida cada modelo reproduce el comportamiento de la variable dependiente, y además, el peso o influencia de cada variable en el modelo respecto a las otras. Los resultados fundamentales del estudio se expresan en dos modelos finales alternativos, que demuestran tener una significativa asociación con el porcentaje de viajes a pie (R2 = 0,6789, p <0,0001), al explicar las dos terceras partes de su variabilidad. En ellos, y en general en todo el estudio realizado, se da una influencia constante de tres índices en particular, que quedan como los principales. Dos de ellos, de acuerdo con muchos de los estudios previos, corresponden a la densidad y la mezcla de usos del suelo. Pero lo más novedoso de los resultados obtenidos es que el tercero es una medida de la accesibilidad de la red de espacio público, algo de lo que no había referencias hasta ahora. Pero, ¿cuál es la definición precisa y el peso relativo de cada uno en el modelo, es decir, en la variable independiente? El de mayor peso en la mayor parte de los análisis realizados es el índice de densidad total (n º residentes + n º puestos de trabajo + n º alumnos / Ha). Es decir, una densidad no sólo de población, sino que incluye algunas de las actividades más importantes que pueden darse una zona para generar movilidad a pie. El segundo que mayor peso adquiere, llegando a ser el primero en alguno de los análisis estadísticos efecturados, es el índice de accesibuilidad configuracional denominado integración de radio 5. Se trata de una medida de la accesibilidad de la zona, de su centralidad, a la escala de, más un menor, un distrito o comarca. En cuanto al tercero, obtiene una importancia bastante menor que los anteriores, y es que representa la mezcla de usos. En concreto es una medida del equilibrio entre los comercios especializados de venta al por menor y el número de residentes (n º de tiendas especializadas en alimentación, bebidas y tabaco / n º de habitantes). Por lo tanto, estos resultados confirman buena parte de los de estudios anteriores, especialmente los relativas a los usos del suelo, pero al mismo tiempo, apuntan a que la red de espacio público podría tener una influir mayor de la comprobada hasta ahora en la proporción de peatones sobre el resto de modos de transportes. Las razones de por qué esto puede ser así, se discuten ampliamente en las conclusiones. Finalmente se puede precisar que dicha conclusión principal se refiere a viajes de una sola etapa (no multimodales) que se dan en los barrios y zonas del área metropolitana de Madrid. Por supuesto, esta conclusión tiene en la actualidad, una validez limitada, ya que es el resultado de un solo caso — Abstract The research hypothesis for this Ph.D. Thesis is that some characteristics of the built environment, particularly those describing the accessibility of the public space network, could be associated with the proportion of pedestrians in all trips (modal split), found in the different parts of a city. The underlying idea is that walking trips are more sensitive to built environment than those by other transport modes, such as for example those by bicycle or by public transport, which could be more conditioned by, e.g. infrastructure availability or service frequency and quality. On the other hand, it has to be noted that the previously research on this field, in particular within Space Syntax’s where this study can be referred, have tested similar hypothesis using pedestrian volumes as the dependent variable, but never against modal split. According to such hypothesis, research methodology is based primarily on empirical quantitative analysis, and it is meant to be able to assess whether public space network, no matter other built environment and non-built environment variables, could have a relationship with the proportion of pedestrian trips in urban areas. Statistical techniques are used to check the association of independent variables with the percentage of walking in all trips, the dependent one. Broadly speaking this methodology is similar to that of previous studies in the field such as CERVERO&KOCKLEMAN (1997), CERVERO & DUNCAN (2003), or to those used mainly in the general review of T.R.B. (2005) or, more recently in ZEGRAS (2006) or CHATMAN (2009). Other methodological options such as stated choice methods (see LOUVIERE, HENSHER & SWAIT, 2000) or agent based analysis (PENN & TURNER, 2004), were discarded, due to a number of reasons, too long to be described here. The case study is not the entire Madrid’s metropolitan area, but almost (4.132.820 inhabitants, about 80% of region´s population). Main data sources are the Regional Mobility Home Based Survey 2004 (EDM04), which is the last available (sample: >35.000 families, > 95.000 individuals), and a spatial model of the metropolitan area, developed using Space Syntax and G.I.S. techniques. The analysis unit, in this case spatial units, are both transport zones (mean population = 7.063) and neighborhoods (mean population = 26.466). The variables of the built environment are clearly the core of the study. A total of 20 (out of 21) are selected from among those found in the literature while, at the same time, being accessible. Nine out of them are used to describe land use characteristics while another eleven describe the network of public spaces. Latter ones include configurational accessibility or Space Syntax variables. This is a particular sort of accessibility related with the concept of configuration, by HILLIER (1996), one of the main authors of Space Syntax, But it also include more customary variables used in mobility research to describe the urban design or spatial structure (here public space network), which here are called geometric characteristics of the such as its length, type of intersection, conectivity, density, etc. Finally a single socioeconomic variable was included in order to assess the influence non built environment factors that also may have an impact on walking (age, income, motorization rate, etc.). The association among variables is worked out using bi-variate correlation analysis and multivariate-analysis. Correlations are calculated among the 21 independent variables and the dependent one, the percentage of walking trips. Then, three types of multi-variate studies are run: general linear, curvilinear and discriminant multi-variate analysis. The latter are fully capable of generating complex association models among several variables, assessing quite precisely to what extent each model reproduces the behavior of the dependent variable, and also the weight or influence of each variable in the model. This study’s results show a consistent influence of three particular indexes in the two final alternative models of the multi-variate study (best, R2=0,6789, p<0,0000). Not surprisingly, two of them correspond to density and mix of land uses. But perhaps more interesting is that the third one is a measure of the accessibility of the public space network, a variable less important in the literature up to now. Additional precisions about them and their relative weight could also be of some interest. The density index is not only about population but includes most important activities in an area (nº residents + nº jobs+ nº students/Ha). The configurational index (radius 5 integration) is a measure of the accessibility of the area, i.e. centrality, at the scale of, more a less, a district. Regarding the mix of land uses index, this one is a measure of the balance between retail, in fact local basic retail, and the number of residents (nº of convenience shops / nº of residents). Referring to their weights, configurational index (radius 5 integration) gets the higher standardized coefficient of the final equation. However, in the final equations, there are a higher number of indexes coming from the density or land use mix categories than from public space network enter. Therefore, these findings seem to support part of the field’s knowledge, especially those concerning land uses, but at the same time they seem to bring in the idea that the configuration of the urban grid could have an influence in the proportion of walkers (as a part of total trips on any transport mode) that do single journey trips in the neighborhoods of Madrid, Spain. Of course this conclusion has, at present, a limited validity since it’s the result of a single case. The reasons of why this can be so, are discussed in the last part of the thesis.

Cryptographic properties of Boolean functions defining elementary cellular automata

In this work, the algebraic properties of the local transition functions of elementary cellular automata (ECA) were analysed. Specifically, a classification of such cellular automata was done according to their algebraic degree, the balancedness, the resiliency, nonlinearity, the propagation criterion and the existence of non-zero linear structures. It is shown that there is not any ECA satisfying all properties at the same time.

Analysis of service-oriented computing systems

La computación basada en servicios (Service-Oriented Computing, SOC) se estableció como un paradigma ampliamente aceptado para el desarollo de sistemas de software flexibles, distribuidos y adaptables, donde las composiciones de los servicios realizan las tareas más complejas o de nivel más alto, frecuentemente tareas inter-organizativas usando los servicios atómicos u otras composiciones de servicios. En tales sistemas, las propriedades de la calidad de servicio (Quality of Service, QoS), como la rapídez de procesamiento, coste, disponibilidad o seguridad, son críticas para la usabilidad de los servicios o sus composiciones en cualquier aplicación concreta. El análisis de estas propriedades se puede realizarse de una forma más precisa y rica en información si se utilizan las técnicas de análisis de programas, como el análisis de complejidad o de compartición de datos, que son capables de analizar simultáneamente tanto las estructuras de control como las de datos, dependencias y operaciones en una composición. El análisis de coste computacional para la composicion de servicios puede ayudar a una monitorización predictiva así como a una adaptación proactiva a través de una inferencia automática de coste computacional, usando los limites altos y bajos como funciones del valor o del tamaño de los mensajes de entrada. Tales funciones de coste se pueden usar para adaptación en la forma de selección de los candidatos entre los servicios que minimizan el coste total de la composición, basado en los datos reales que se pasan al servicio. Las funciones de coste también pueden ser combinadas con los parámetros extraídos empíricamente desde la infraestructura, para producir las funciones de los límites de QoS sobre los datos de entrada, cuales se pueden usar para previsar, en el momento de invocación, las violaciones de los compromisos al nivel de servicios (Service Level Agreements, SLA) potenciales or inminentes. En las composiciones críticas, una previsión continua de QoS bastante eficaz y precisa se puede basar en el modelado con restricciones de QoS desde la estructura de la composition, datos empiricos en tiempo de ejecución y (cuando estén disponibles) los resultados del análisis de complejidad. Este enfoque se puede aplicar a las orquestaciones de servicios con un control centralizado del flujo, así como a las coreografías con participantes multiples, siguiendo unas interacciones complejas que modifican su estado. El análisis del compartición de datos puede servir de apoyo para acciones de adaptación, como la paralelización, fragmentación y selección de los componentes, las cuales son basadas en dependencias funcionales y en el contenido de información en los mensajes, datos internos y las actividades de la composición, cuando se usan construcciones de control complejas, como bucles, bifurcaciones y flujos anidados. Tanto las dependencias funcionales como el contenido de información (descrito a través de algunos atributos definidos por el usuario) se pueden expresar usando una representación basada en la lógica de primer orden (claúsulas de Horn), y los resultados del análisis se pueden interpretar como modelos conceptuales basados en retículos. ABSTRACT Service-Oriented Computing (SOC) is a widely accepted paradigm for development of flexible, distributed and adaptable software systems, in which service compositions perform more complex, higher-level, often cross-organizational tasks using atomic services or other service compositions. In such systems, Quality of Service (QoS) properties, such as the performance, cost, availability or security, are critical for the usability of services and their compositions in concrete applications. Analysis of these properties can become more precise and richer in information, if it employs program analysis techniques, such as the complexity and sharing analyses, which are able to simultaneously take into account both the control and the data structures, dependencies, and operations in a composition. Computation cost analysis for service composition can support predictive monitoring and proactive adaptation by automatically inferring computation cost using the upper and lower bound functions of value or size of input messages. These cost functions can be used for adaptation by selecting service candidates that minimize total cost of the composition, based on the actual data that is passed to them. The cost functions can also be combined with the empirically collected infrastructural parameters to produce QoS bounds functions of input data that can be used to predict potential or imminent Service Level Agreement (SLA) violations at the moment of invocation. In mission-critical applications, an effective and accurate continuous QoS prediction, based on continuations, can be achieved by constraint modeling of composition QoS based on its structure, known data at runtime, and (when available) the results of complexity analysis. This approach can be applied to service orchestrations with centralized flow control, and choreographies with multiple participants with complex stateful interactions. Sharing analysis can support adaptation actions, such as parallelization, fragmentation, and component selection, which are based on functional dependencies and information content of the composition messages, internal data, and activities, in presence of complex control constructs, such as loops, branches, and sub-workflows. Both the functional dependencies and the information content (described using user-defined attributes) can be expressed using a first-order logic (Horn clause) representation, and the analysis results can be interpreted as a lattice-based conceptual models.

Functional and effective connectivity in MEG. Application to the study of epilepsy

Nuestro cerebro contiene cerca de 1014 sinapsis neuronales. Esta enorme cantidad de conexiones proporciona un entorno ideal donde distintos grupos de neuronas se sincronizan transitoriamente para provocar la aparición de funciones cognitivas, como la percepción, el aprendizaje o el pensamiento. Comprender la organización de esta compleja red cerebral en base a datos neurofisiológicos, representa uno de los desafíos más importantes y emocionantes en el campo de la neurociencia. Se han propuesto recientemente varias medidas para evaluar cómo se comunican las diferentes partes del cerebro a diversas escalas (células individuales, columnas corticales, o áreas cerebrales). Podemos clasificarlos, según su simetría, en dos grupos: por una parte, la medidas simétricas, como la correlación, la coherencia o la sincronización de fase, que evalúan la conectividad funcional (FC); mientras que las medidas asimétricas, como la causalidad de Granger o transferencia de entropía, son capaces de detectar la dirección de la interacción, lo que denominamos conectividad efectiva (EC). En la neurociencia moderna ha aumentado el interés por el estudio de las redes funcionales cerebrales, en gran medida debido a la aparición de estos nuevos algoritmos que permiten analizar la interdependencia entre señales temporales, además de la emergente teoría de redes complejas y la introducción de técnicas novedosas, como la magnetoencefalografía (MEG), para registrar datos neurofisiológicos con gran resolución. Sin embargo, nos hallamos ante un campo novedoso que presenta aun varias cuestiones metodológicas sin resolver, algunas de las cuales trataran de abordarse en esta tesis. En primer lugar, el creciente número de aproximaciones para determinar la existencia de FC/EC entre dos o más señales temporales, junto con la complejidad matemática de las herramientas de análisis, hacen deseable organizarlas todas en un paquete software intuitivo y fácil de usar. Aquí presento HERMES (http://hermes.ctb.upm.es), una toolbox en MatlabR, diseñada precisamente con este fin. Creo que esta herramienta será de gran ayuda para todos aquellos investigadores que trabajen en el campo emergente del análisis de conectividad cerebral y supondrá un gran valor para la comunidad científica. La segunda cuestión practica que se aborda es el estudio de la sensibilidad a las fuentes cerebrales profundas a través de dos tipos de sensores MEG: gradiómetros planares y magnetómetros, esta aproximación además se combina con un enfoque metodológico, utilizando dos índices de sincronización de fase: phase locking value (PLV) y phase lag index (PLI), este ultimo menos sensible a efecto la conducción volumen. Por lo tanto, se compara su comportamiento al estudiar las redes cerebrales, obteniendo que magnetómetros y PLV presentan, respectivamente, redes más densamente conectadas que gradiómetros planares y PLI, por los valores artificiales que crea el problema de la conducción de volumen. Sin embargo, cuando se trata de caracterizar redes epilépticas, el PLV ofrece mejores resultados, debido a la gran dispersión de las redes obtenidas con PLI. El análisis de redes complejas ha proporcionado nuevos conceptos que mejoran caracterización de la interacción de sistemas dinámicos. Se considera que una red está compuesta por nodos, que simbolizan sistemas, cuyas interacciones se representan por enlaces, y su comportamiento y topología puede caracterizarse por un elevado número de medidas. Existe evidencia teórica y empírica de que muchas de ellas están fuertemente correlacionadas entre sí. Por lo tanto, se ha conseguido seleccionar un pequeño grupo que caracteriza eficazmente estas redes, y condensa la información redundante. Para el análisis de redes funcionales, la selección de un umbral adecuado para decidir si un determinado valor de conectividad de la matriz de FC es significativo y debe ser incluido para un análisis posterior, se convierte en un paso crucial. En esta tesis, se han obtenido resultados más precisos al utilizar un test de subrogadas, basado en los datos, para evaluar individualmente cada uno de los enlaces, que al establecer a priori un umbral fijo para la densidad de conexiones. Finalmente, todas estas cuestiones se han aplicado al estudio de la epilepsia, caso práctico en el que se analizan las redes funcionales MEG, en estado de reposo, de dos grupos de pacientes epilépticos (generalizada idiopática y focal frontal) en comparación con sujetos control sanos. La epilepsia es uno de los trastornos neurológicos más comunes, con más de 55 millones de afectados en el mundo. Esta enfermedad se caracteriza por la predisposición a generar ataques epilépticos de actividad neuronal anormal y excesiva o bien síncrona, y por tanto, es el escenario perfecto para este tipo de análisis al tiempo que presenta un gran interés tanto desde el punto de vista clínico como de investigación. Los resultados manifiestan alteraciones especificas en la conectividad y un cambio en la topología de las redes en cerebros epilépticos, desplazando la importancia del ‘foco’ a la ‘red’, enfoque que va adquiriendo relevancia en las investigaciones recientes sobre epilepsia. ABSTRACT There are about 1014 neuronal synapses in the human brain. This huge number of connections provides the substrate for neuronal ensembles to become transiently synchronized, producing the emergence of cognitive functions such as perception, learning or thinking. Understanding the complex brain network organization on the basis of neuroimaging data represents one of the most important and exciting challenges for systems neuroscience. Several measures have been recently proposed to evaluate at various scales (single cells, cortical columns, or brain areas) how the different parts of the brain communicate. We can classify them, according to their symmetry, into two groups: symmetric measures, such as correlation, coherence or phase synchronization indexes, evaluate functional connectivity (FC); and on the other hand, the asymmetric ones, such as Granger causality or transfer entropy, are able to detect effective connectivity (EC) revealing the direction of the interaction. In modern neurosciences, the interest in functional brain networks has increased strongly with the onset of new algorithms to study interdependence between time series, the advent of modern complex network theory and the introduction of powerful techniques to record neurophysiological data, such as magnetoencephalography (MEG). However, when analyzing neurophysiological data with this approach several questions arise. In this thesis, I intend to tackle some of the practical open problems in the field. First of all, the increase in the number of time series analysis algorithms to study brain FC/EC, along with their mathematical complexity, creates the necessity of arranging them into a single, unified toolbox that allow neuroscientists, neurophysiologists and researchers from related fields to easily access and make use of them. I developed such a toolbox for this aim, it is named HERMES (http://hermes.ctb.upm.es), and encompasses several of the most common indexes for the assessment of FC and EC running for MatlabR environment. I believe that this toolbox will be very helpful to all the researchers working in the emerging field of brain connectivity analysis and will entail a great value for the scientific community. The second important practical issue tackled in this thesis is the evaluation of the sensitivity to deep brain sources of two different MEG sensors: planar gradiometers and magnetometers, in combination with the related methodological approach, using two phase synchronization indexes: phase locking value (PLV) y phase lag index (PLI), the latter one being less sensitive to volume conduction effect. Thus, I compared their performance when studying brain networks, obtaining that magnetometer sensors and PLV presented higher artificial values as compared with planar gradiometers and PLI respectively. However, when it came to characterize epileptic networks it was the PLV which gives better results, as PLI FC networks where very sparse. Complex network analysis has provided new concepts which improved characterization of interacting dynamical systems. With this background, networks could be considered composed of nodes, symbolizing systems, whose interactions with each other are represented by edges. A growing number of network measures is been applied in network analysis. However, there is theoretical and empirical evidence that many of these indexes are strongly correlated with each other. Therefore, in this thesis I reduced them to a small set, which could more efficiently characterize networks. Within this framework, selecting an appropriate threshold to decide whether a certain connectivity value of the FC matrix is significant and should be included in the network analysis becomes a crucial step, in this thesis, I used the surrogate data tests to make an individual data-driven evaluation of each of the edges significance and confirmed more accurate results than when just setting to a fixed value the density of connections. All these methodologies were applied to the study of epilepsy, analysing resting state MEG functional networks, in two groups of epileptic patients (generalized and focal epilepsy) that were compared to matching control subjects. Epilepsy is one of the most common neurological disorders, with more than 55 million people affected worldwide, characterized by its predisposition to generate epileptic seizures of abnormal excessive or synchronous neuronal activity, and thus, this scenario and analysis, present a great interest from both the clinical and the research perspective. Results revealed specific disruptions in connectivity and network topology and evidenced that networks’ topology is changed in epileptic brains, supporting the shift from ‘focus’ to ‘networks’ which is gaining importance in modern epilepsy research.

From plant surface to plant metabolism: the uncertain fate of foliar-applied nutrients

The application of agrochemical sprays to the aerial parts of crop plants is an important agricultural practice world-wide. While variable effectiveness is often seen in response to foliar treatments, there is abundant evidence showing the beneficial effect of foliar fertilizers in terms of improving the metabolism, quality, and yields of crops. This mini-review is focused on the major bottlenecks associated with the uptake and translocation of foliar-applied nutrient solutions. A better understanding of the complex scenario surrounding the ultimate delivery of foliar-applied nutrients to sink cells and organs is essential for improving the effectiveness and performance of foliar fertilizers.

Real-time determination of airspace wind speed field based on FDR = Determinación en tiempo real del campo de velocidad del viento basado en registradores de datos de vuelo

El estudio del comportamiento de la atmósfera ha resultado de especial importancia tanto en el programa SESAR como en NextGen, en los que la gestión actual del tránsito aéreo (ATM) está experimentando una profunda transformación hacia nuevos paradigmas tanto en Europa como en los EE.UU., respectivamente, para el guiado y seguimiento de las aeronaves en la realización de rutas más eficientes y con mayor precisión. La incertidumbre es una característica fundamental de los fenómenos meteorológicos que se transfiere a la separación de las aeronaves, las trayectorias de vuelo libres de conflictos y a la planificación de vuelos. En este sentido, el viento es un factor clave en cuanto a la predicción de la futura posición de la aeronave, por lo que tener un conocimiento más profundo y preciso de campo de viento reducirá las incertidumbres del ATC. El objetivo de esta tesis es el desarrollo de una nueva técnica operativa y útil destinada a proporcionar de forma adecuada y directa el campo de viento atmosférico en tiempo real, basada en datos de a bordo de la aeronave, con el fin de mejorar la predicción de las trayectorias de las aeronaves. Para lograr este objetivo se ha realizado el siguiente trabajo. Se han descrito y analizado los diferentes sistemas de la aeronave que proporcionan las variables necesarias para obtener la velocidad del viento, así como de las capacidades que permiten la presentación de esta información para sus aplicaciones en la gestión del tráfico aéreo. Se ha explorado el uso de aeronaves como los sensores de viento en un área terminal para la estimación del viento en tiempo real con el fin de mejorar la predicción de las trayectorias de aeronaves. Se han desarrollado métodos computacionalmente eficientes para estimar las componentes horizontales de la velocidad del viento a partir de las velocidades de las aeronaves (VGS, VCAS/VTAS), la presión y datos de temperatura. Estos datos de viento se han utilizado para estimar el campo de viento en tiempo real utilizando un sistema de procesamiento de datos a través de un método de mínima varianza. Por último, se ha evaluado la exactitud de este procedimiento para que esta información sea útil para el control del tráfico aéreo. La información inicial proviene de una muestra de datos de Registradores de Datos de Vuelo (FDR) de aviones que aterrizaron en el aeropuerto Madrid-Barajas. Se dispuso de datos de ciertas aeronaves durante un periodo de más de tres meses que se emplearon para calcular el vector viento en cada punto del espacio aéreo. Se utilizó un modelo matemático basado en diferentes métodos de interpolación para obtener los vectores de viento en áreas sin datos disponibles. Se han utilizado tres escenarios concretos para validar dos métodos de interpolación: uno de dos dimensiones que trabaja con ambas componentes horizontales de forma independiente, y otro basado en el uso de una variable compleja que relaciona ambas componentes. Esos métodos se han probado en diferentes escenarios con resultados dispares. Esta metodología se ha aplicado en un prototipo de herramienta en MATLAB © para analizar automáticamente los datos de FDR y determinar el campo vectorial del viento que encuentra la aeronave al volar en el espacio aéreo en estudio. Finalmente se han obtenido las condiciones requeridas y la precisión de los resultados para este modelo. El método desarrollado podría utilizar los datos de los aviones comerciales como inputs utilizando los datos actualmente disponibles y la capacidad computacional, para proporcionárselos a los sistemas ATM donde se podría ejecutar el método propuesto. Estas velocidades del viento calculadas, o bien la velocidad respecto al suelo y la velocidad verdadera, se podrían difundir, por ejemplo, a través del sistema de direccionamiento e informe para comunicaciones de aeronaves (ACARS), mensajes de ADS-B o Modo S. Esta nueva fuente ayudaría a actualizar la información del viento suministrada en los productos aeronáuticos meteorológicos (PAM), informes meteorológicos de aeródromos (AIRMET), e información meteorológica significativa (SIGMET). ABSTRACT The study of the atmosphere behaviour is been of particular importance both in SESAR and NextGen programs, where the current air traffic management (ATM) system is undergoing a profound transformation to the new paradigms both in Europe and the USA, respectively, to guide and track aircraft more precisely on more efficient routes. Uncertainty is a fundamental characteristic of weather phenomena which is transferred to separation assurance, flight path de-confliction and flight planning applications. In this respect, the wind is a key factor regarding the prediction of the future position of the aircraft, so that having a deeper and accurate knowledge of wind field will reduce ATC uncertainties. The purpose of this thesis is to develop a new and operationally useful technique intended to provide adequate and direct real-time atmospheric winds fields based on on-board aircraft data, in order to improve aircraft trajectory prediction. In order to achieve this objective the following work has been accomplished. The different sources in the aircraft systems that provide the variables needed to derivate the wind velocity have been described and analysed, as well as the capabilities which allow presenting this information for air traffic management applications. The use of aircraft as wind sensors in a terminal area for real-time wind estimation in order to improve aircraft trajectory prediction has been explored. Computationally efficient methods have been developed to estimate horizontal wind components from aircraft velocities (VGS, VCAS/VTAS), pressure, and temperature data. These wind data were utilized to estimate a real-time wind field using a data processing approach through a minimum variance method. Finally, the accuracy of this procedure has been evaluated for this information to be useful to air traffic control. The initial information comes from a Flight Data Recorder (FDR) sample of aircraft landing in Madrid-Barajas Airport. Data available for more than three months were exploited in order to derive the wind vector field in each point of the airspace. Mathematical model based on different interpolation methods were used in order to obtain wind vectors in void areas. Three particular scenarios were employed to test two interpolation methods: a two-dimensional one that works with both horizontal components in an independent way, and also a complex variable formulation that links both components. Those methods were tested using various scenarios with dissimilar results. This methodology has been implemented in a prototype tool in MATLAB © in order to automatically analyse FDR and determine the wind vector field that aircraft encounter when flying in the studied airspace. Required conditions and accuracy of the results were derived for this model. The method developed could be fed by commercial aircraft utilizing their currently available data sources and computational capabilities, and providing them to ATM system where the proposed method could be run. Computed wind velocities, or ground and true airspeeds, would then be broadcasted, for example, via the Aircraft Communication Addressing and Reporting System (ACARS), ADS-B out messages, or Mode S. This new source would help updating the wind information furnished in meteorological aeronautical products (PAM), meteorological aerodrome reports (AIRMET), and significant meteorological information (SIGMET).