Biodiversidad de frijol (Phaseolus vulgaris L.) en Honduras. Caracterización agromorfológica, molecular y ecogeográfica

En el presente trabajo se ha llevado a cabo un estudio de la biodiversidad del frijol común (Phaseolus vulgaris L.) en Honduras, que es el segundo de los cultivos de granos básicos en importancia. Dicho estudio se ha realizado mediante una caracterización agromorfológica, molecular y ecogeográfica en una selección de 300 accesiones conservadas en el banco de germoplasma ubicado en la Escuela Agrícola Panamericana (EAP) El Zamorano, y que se colectaron en 13 departamentos del país durante el periodo de 1990 a 1994. Estas accesiones fueron colectadas cuatro años antes del acontecimiento del huracán Mitch, el cual a su paso afectó al 96% del área total cultivable en su momento, lo cual nos hace considerar que la biodiversidad de razas locales (landraces) de frijol común existentes in situ fueron severamente afectadas. Los trabajos dirigidos a analizar la biodiversidad de razas locales de frijol común en Honduras son escasos, y este trabajo se constituye como el primero que incluye una amplia muestra a ser estudiada a través de una caracterización en tres aspectos complementarios (agromorfológico, molecular y ecogeográfico). Se evaluaron 32 caracteres agromorfológicos, 12 cuantitativos y 20 cualitativos, en distintas partes de la planta. Se establecieron las correlaciones entre los caracteres agromorfológicos y se elaboró un dendrograma con los mismos, en el que se formaron ocho grupos, en parte relacionados principalmente con los colores y tamaños de la semilla. Mediante el análisis de componentes principales se estudiaron los caracteres de más peso en cada uno de los tres primeros componentes. Asimismo, se estudiaron las correlaciones entre caracteres, siendo las más altas la longitud y anchura de la hoja, días a madurez y a cosecha y longitud y peso de semilla. Por otra parte, el mapa de diversidad agromorfológica mostró la existencia de tres zonas con mayor diversidad: en el oeste (en los departamentos de Santa Bárbara, Lempira y Copán), en el centro-norte (en los departamentos de Francisco Morazán, Yoro y Atlántida) y en el sur (en el departamento de El Paraíso y al sur de Francisco Morazán). Para la caracterización molecular partimos de 12 marcadores de tipo microsatélite, evaluados en 54 accesiones, que fueron elegidas por constituir grupos que compartían un mismo nombre local. Finalmente, se seleccionaron los cuatro microsatélites (BM53, GATS91, BM211 y PV-AT007) que resultaron ser más polimórficos e informativos para el análisis de las 300 accesiones, con los que se detectaron un total de 119 alelos (21 de ellos únicos o privados de accesión) y 256 patrones alélicos diferentes. Para estudiar la estructura y relaciones genéticas en las 300 accesiones se incluyeron en el análisis tres controles o accesiones de referencia, pertenecientes dos de ellas al acervo genético Andino y una al Mesoamericano. En el dendrograma se obtuvieron 25 grupos de accesiones con idénticas combinaciones de alelos. Al comparar este dendrograma con el de caracteres agromorfológicos se observaron diversos grupos con marcada similitud en ambos. Un total de 118 accesiones resultaron ser homogéneas y homocigóticas, a la vez que representativas del grupo de 300 accesiones, por lo que se analizaron con más detalle. El análisis de la estructura genética definió la formación de dos grupos, supuestamente relacionados con los acervos genéticos Andino (48) y Mesoamericano (61), y un reducido número de accesiones (9) que podrían tener un origen híbrido, debido a la existencia de un cierto grado de introgresión entre ambos acervos. La diferenciación genética entre ambos grupos fue del 13,3%. Asimismo, 66 de los 82 alelos detectados fueron privados de grupo, 30 del supuesto grupo Andino y 36 del Mesoamericano. Con relación al mapa de diversidad molecular, presentó una distribución bastante similar al de la diversidad agromorfológica, detectándose también las zonas de mayor diversidad genética en el oeste (en los departamentos de Lempira y Santa Bárbara), en el centro-norte (en los departamentos de Yoro y Atlántida) y en el sur (en el departamento de El Paraíso y al sur de Francisco Morazán). Para la caracterización ecogeográfica se seleccionaron variables de tipo bioclimático (2), geofísico (2) y edáfico (8), y mediante el método de agrupamiento de partición alrededor de los medoides, la combinación de los grupos con cada uno de los tres tipos de variables definió un total de 32 categorías ecogeográficas en el país, detectándose accesiones en 16 de ellas. La distribución de las accesiones previsiblemente esté relacionada con la existencia de condiciones más favorables al cultivo de frijol. En el mapa de diversidad ecogeográfica, nuevamente, se observaron varias zonas con alta diversidad tanto en el oeste, como en el centro-norte y en el sur del país. Como consecuencia del estudio realizado, se concluyó la existencia de una marcada biodiversidad en el material analizado, desde el punto de vista tanto agromorfológico como molecular. Por lo que resulta de gran importancia plantear la conservación de este patrimonio genético tanto ex situ, en bancos de germoplasma, como on farm, en las propias explotaciones de los agricultores del país, siempre que sea posible. ABSTRACT In the present work we have carried out a study of the biodiversity of the common bean (Phaseolus vulgaris L) in Honduras, which is the second of the basic grain crops in importance. This study was conducted through agro-morphological, molecular and ecogeographical characterization of a selection of 300 accessions conserved in the genebank located in the ‘Escuela Agrícola Panamericana (EAP) El Zamorano’ that were collected in 13 departments of the country during the 1990 to 1994 period. These accessions were collected four years before the occurrence of Mitch hurricane, which affected 96% of the total cultivable area at the time, which makes us to consider that the biodiversity of local landraces of common bean existing in situ were severely affected. The work aimed to analyze the biodiversity of local races of common bean in Honduras are scarce, and this work constitutes the first to include a large sample to be studied through a characterization on three complementary aspects (agromorphological, molecular and ecogeographical). Thirty two agromorphological characters, 12 quantitative and 20 qualitative, in various parts of the plant were evaluated. Correlations between agromorphological characters were established and a dendrogram with them was constructed, in which eight groups were formed, in part mainly related to the colors and sizes of the seeds. By principal component analysis the characters with more weight in each of the first three components were studied. Also, correlations between characters were studied, the highest of them being length and leaf width, days to maturity and harvest, and seed length and weight. Moreover, the map of agromorphological diversity showed the existence of three areas with more diversity: the west (departments of Santa Barbara, Copan and Lempira), the center-north (departments of Francisco Morazán, Yoro and Atlántida) and the south (department of El Paraiso and south of Francisco Morazán). For molecular characterization we started with 12 microsatellite markers, evaluated in 54 accessions, which were chosen because they formed groups that shared the same local name. Finally, four microsatellites (BM53, GATS91, BM211 and PV-AT007) were selected for the analysis of 300 accessions, since they were the most polymorphic and informative. They gave a total of 119 alleles (21 of them unique or private for the accession) and 256 different allelic patterns. To study the structure and genetic relationships in the 300 accessions, three controls or accessions of reference were included in the analysis: two of them belonging to the Andean gene pool and one to the Mesoamerican. In the dendrogram, 25 accession groups with identical allele combinations were obtained. Comparing this dendrogram to the obtained with agromorphological characters, several groups with marked similarity in both were observed. A total of 118 accessions were homozygous and homogeneous, while representing the group of 300 accessions, therefore they were analyzed in more detail. The analysis of the genetic structure defined the formation of two groups, supposedly related to the Andean (48) and the Mesoamerican (61) gene pools, and a small number of accessions (9) which may have a hybrid origin, due to the existence of some degree of introgression between both gene pools. Genetic differentiation between both groups was 13.3%. Also, 66 of the 82 detected alleles were private or unique for the group, 30 of the supposed Andean group and 36 of the Mesoamerican. With relation to the map of molecular diversity, it showed a quite similar distribution to the agromorphological, also detecting the areas of greatest genetic diversity in the west (departments of Lempira and Santa Bárbara), in the center-north (departments Atlántida and Yoro) and in the south (departments of El Paraíso and south of Francisco Morazán). For the ecogeographical characterization, bioclimatic (2), geophysical (2) and edaphic (8) variables were selected, and by the method of clustering partition around the medoids, the combination of the groups to each of the three types of variables defined a total of 32 ecogeographical categories in the country, having accessions in 16 of them. The distribution of accessions is likely related to the existence of more favorable conditions for the cultivation of beans. The map of ecogeographical diversity, again, several areas with high diversity both in the west and in the center-north and in the south of the country were observed. As a result of study, the existence of marked biodiversity in the analyzed material was concluded, both from the agromorphological and from the molecular point of view. Consequently it is very important to propose the conservation of this genetic heritage both ex situ, in genebanks, as on farm, in the holdings of the farmers of the country, whenever possible.

Tratamiento magnético como técnica estimulante de la germinación de semillas de soja

El campo magnético está presente en la naturaleza y afecta al reino vegetal y animal. El bioelectromagnetismo es el estudio de los efectos de los campos electromagnéticos sobre los sistemas biológicos y sus interacciones con los campos magnéticos naturales y artificiales. Estudios previos han determinado que la aplicación de campos magnéticos modifica la velocidad de germinación y el desarrollo de plantas. La soja constituye uno de los cultivos oleaginosos más importantes del mundo: en las últimas décadas su popularidad se ha incrementado por considerar esta leguminosa como una fuente de salud. Este estudio se centra en el efecto que se puede producir en la germinación de la soja (Glycine max L.) al exponerlo a campos magnéticos estacionarios de 125 mT, 250 mT y 300 mT durante diferentes tiempos de exposición (10¿,20¿, 1 hora, 24 horas y de forma permanente). Se observan diferencias muy significativas tanto en la velocidad de germinación como en el porcentaje de semillas germinadas entre las semillas control y las tratadas. Estos resultados ofrecen nuevas posibilidades para una más efectiva producción vegetal.

Estudio de la cadena de valor de la producción hortícola en Loukolela. República del Congo.

El presente proyecto, titulado “Estudio de la cadena de valor de la producción hortícola en Loukolela, República del Congo” fue llevado a cabo durante Febrero de 2014 a Junio de 2014 en la población de Loukolela; una zona muy rural al norte del país. El proyecto se realizó con la colaboración del grupo de cooperación de la ETSI Agrónomos: RISE, la Fundación “Carmen Rodríguez Granda”, la fundación de contraparte en el Congo “Carmen Granda Rodríguez” -que surge de la fundación española con sede en Madrid- y la ayuda de los tutores en España: Isabel de Felipe Boente, y en el Congo: Obeth Harelimana Pascal. Su objetivo principal fue, por un lado, el desarrollo de la comunidad agrícola de Loukolela, utilizando como instrumento el estudio de la cadena de valor de la producción hortícola, localizando los eslabones más débiles y obteniendo información de los distintos mercados de la zona y los distintos agentes de esta. Y por otro, un trabajo sobre el terreno junto a la Cooperativa “Virgen del Carmen” (ACVC) en la que se puso en marcha un huerto piloto cuyo objetivo era diversificar los cultivos presentes en la región con otros que también se dan bien en el terreno; y estudiar sus salidas comerciales y para garantizar la seguridad alimentaria y diversificar la básica dieta de los habitantes de la localidad interesados. Este documento presenta una descripción de las principales actividades que se desarrollaron con el fin de tratar de alcanzar el objetivo planteado y mejorar las condiciones de vida de los habitantes de esta zona.

Integration of physical, chemical and biological tactics against insect pests and virus diseases in horticultural crops

Actualmente, la gestión de sistemas de Manejo Integrado de Plagas (MIP) en cultivos hortícolas tiene por objetivo priorizar los métodos de control no químicos en detrimento del consumo de plaguicidas, según recoge la directiva europea 2009/128/CE ‘Uso Sostenible de Plaguicidas’ (OJEC, 2009). El uso de agentes de biocontrol como alternativa a la aplicación de insecticidas es un elemento clave de los sistemas MIP por sus innegables ventajas ambientales que se utiliza ampliamente en nuestro país (Jacas y Urbaneja, 2008). En la región de Almería, donde se concentra el 65% de cultivo en invernadero de nuestro país (47.367 ha), MIP es la principal estrategia en pimiento (MAGRAMA, 2014), y comienza a serlo en otros cultivos como tomate o pepino. El cultivo de pepino, con 8.902 ha (MAGRAMA, 2013), tiene un protocolo semejante al pimiento (Robledo et al., 2009), donde la única especie de pulgón importante es Aphis gossypii Glover. Sin embargo, pese al continuo incremento de la superficie de cultivo agrícola bajo sistemas MIP, los daños originados por virosis siguen siendo notables. Algunos de los insectos presentes en los cultivos de hortícolas son importantes vectores de virus, como los pulgones, las moscas blancas o los trips, cuyo control resulta problemático debido a su elevada capacidad para transmitir virus vegetales incluso a una baja densidad de plaga (Holt et al., 2008; Jacas y Urbaneja, 2008). Las relaciones que se establecen entre los distintos agentes de un ecosistema son complejas y muy específicas. Se ha comprobado que, pese a que los enemigos naturales reducen de manera beneficiosa los niveles de plaga, su incorporación en los sistemas planta-insecto-virus puede desencadenar complicadas interacciones con efectos no deseables (Dicke y van Loon, 2000; Jeger et al., 2011). Así, los agentes de biocontrol también pueden inducir a que los insectos vectores modifiquen su comportamiento como respuesta al ataque y, con ello, el grado de dispersión y los patrones de distribución de las virosis que transmiten (Bailey et al., 1995; Weber et al., 1996; Hodge y Powell, 2008a; Hodge et al., 2011). Además, en ocasiones el control biológico por sí solo no es suficiente para controlar determinadas plagas (Medina et al., 2008). Entre los métodos que se pueden aplicar bajo sistemas MIP están las barreras físicas que limitan la entrada de plagas al interior de los invernaderos o interfieren con su movimiento, como pueden ser las mallas anti-insecto (Álvarez et al., 2014), las mallas fotoselectivas (Raviv y Antignus, 2004; Weintraub y Berlinger, 2004; Díaz y Fereres, 2007) y las mallas impregnadas en insecticida (Licciardi et al., 2008; Martin et al., 2014). Las mallas fotoselectivas reducen o bloquean casi por completo la transmisión de radiación UV, lo que interfiere con la visión de los insectos y dificulta o impide la localización del cultivo y su establecimiento en el mismo (Raviv y Antignus, 2004; Weintraub, 2009). Se ha comprobado cómo su uso puede controlar los pulgones y las virosis en cultivo de lechuga (Díaz et al., 2006; Legarrea et al., 2012a), así como la mosca blanca, los trips y los ácaros, y los virus que estos transmiten en otros cultivos (Costa y Robb, 1999; Antignus et al., 2001; Kumar y Poehling, 2006; Doukas y Payne, 2007a; Legarrea et al., 2010). Sin embargo, no se conoce perfectamente el modo de acción de estas barreras, puesto que existe un efecto directo sobre la plaga y otro indirecto mediado por la planta, cuya fisiología cambia al desarrollarse en ambientes con falta de radiación UV, y que podría afectar al ciclo biológico de los insectos fitófagos (Vänninen et al., 2010; Johansen et al., 2011). Del mismo modo, es necesario estudiar la compatibilidad de esta estrategia con los enemigos naturales de las plagas. Hasta la fecha, los estudios han evidenciado que los agentes de biocontrol pueden realizar su actividad bajo ambientes pobres en radiación UV (Chyzik et al., 2003; Chiel et al., 2006; Doukas y Payne, 2007b; Legarrea et al., 2012c). Otro método basado en barreras físicas son las mallas impregnadas con insecticidas, que se han usado tradicionalmente en la prevención de enfermedades humanas transmitidas por mosquitos (Martin et al., 2006). Su aplicación se ha ensayado en agricultura en ciertos cultivos al aire libre (Martin et al., 2010; Díaz et al., 2004), pero su utilidad en cultivos protegidos para prevenir la entrada de insectos vectores en invernadero todavía no ha sido investigada. Los aditivos se incorporan al tejido durante el proceso de extrusión de la fibra y se liberan lentamente actuando por contacto en el momento en que el insecto aterriza sobre la malla, con lo cual el riesgo medioambiental y para la salud humana es muy limitado. Los plaguicidas que se emplean habitualmente suelen ser piretroides (deltametrina o bifentrín), aunque también se ha ensayado dicofol (Martin et al., 2010) y alfa-cipermetrina (Martin et al., 2014). Un factor que resulta de vital importancia en este tipo de mallas es el tamaño del poro para facilitar una buena ventilación del cultivo, al tiempo que se evita la entrada de insectos de pequeño tamaño como las moscas blancas (Bethke y Paine, 1991; Muñoz et al., 1999). Asimismo, se plantea la necesidad de estudiar la compatibilidad de estas mallas con los enemigos naturales. Es por ello que en esta Tesis Doctoral se plantea la necesidad de evaluar nuevas mallas impregnadas que impidan el paso de insectos de pequeño tamaño al interior de los invernaderos, pero que a su vez mantengan un buen intercambio y circulación de aire a través del poro de la malla. Así, en la presente Tesis Doctoral, se han planteado los siguientes objetivos generales a desarrollar: 1. Estudiar el impacto de la presencia de parasitoides sobre el grado de dispersión y los patrones de distribución de pulgones y las virosis que éstos transmiten. 2. Conocer el efecto directo de ambientes pobres en radiación UV sobre el comportamiento de vuelo de plagas clave de hortícolas y sus enemigos naturales. 3. Evaluar el efecto directo de la radiación UV-A sobre el crecimiento poblacional de pulgones y mosca blanca, y sobre la fisiología de sus plantas hospederas, así como el efecto indirecto de la radiación UV-A en ambas plagas mediado por el crecimiento de dichas planta hospederas. 4. Caracterización de diversas mallas impregnadas en deltametrina y bifentrín con diferentes propiedades y selección de las óptimas para el control de pulgones, mosca blanca y sus virosis asociadas en condiciones de campo. Estudio de su compatibilidad con parasitoides. ABSTRACT Insect vectors of plant viruses are the main agents causing major economic losses in vegetable crops grown under protected environments. This Thesis focuses on the implementation of new alternatives to chemical control of insect vectors under Integrated Pest Management programs. In Spain, biological control is the main pest control strategy used in a large part of greenhouses where horticultural crops are grown. The first study aimed to increase our knowledge on how the presence of natural enemies such as Aphidius colemani Viereck may alter the dispersal of the aphid vector Aphis gossypii Glover (Chapter 4). In addition, it was investigated if the presence of this parasitoid affected the spread of aphid-transmitted viruses Cucumber mosaic virus (CMV, Cucumovirus) and Cucurbit aphid-borne yellows virus (CABYV, Polerovirus) infecting cucumber (Cucumis sativus L). SADIE methodology was used to study the distribution patterns of both the virus and its vector, and their degree of association. Results suggested that parasitoids promoted aphid dispersal in the short term, which enhanced CMV spread, though consequences of parasitism suggested potential benefits for disease control in the long term. Furthermore, A. colemani significantly limited the spread and incidence of the persistent virus CABYV in the long term. The flight activity of pests Myzus persicae (Sulzer), Bemisia tabaci (Gennadius) and Tuta absoluta (Meyrick), and natural enemies A. colemani and Sphaerophoria rueppellii (Weidemann) under UV-deficient environments was studied under field conditions (Chapter 5). One-chamber tunnels were covered with cladding materials with different UV transmittance properties. Inside each tunnel, insects were released from tubes placed in a platform suspended from the ceiling. Specific targets were located at different distances from the platform. The ability of aphids and whiteflies to reach their targets was diminished under UV-absorbing barriers, suggesting a reduction of vector activity under this type of nets. Fewer aphids reached distant traps under UV-absorbing nets, and significantly more aphids could fly to the end of the tunnels covered with non-UV blocking materials. Unlike aphids, differences in B. tabaci captures were mainly found in the closest targets. The oviposition of lepidopteran T. absoluta was also negatively affected by a UV-absorbing cover. The photoselective barriers were compatible with parasitism and oviposition of biocontrol agents. Apart from the direct response of insects to UV radiation, plant-mediated effects influencing insect performance were investigated (Chapter 6). The impact of UV-A radiation on the performance of aphid M. persicae and whitefly B. tabaci, and growth and leaf physiology of host plants pepper and eggplant was studied under glasshouse conditions. Plants were grown inside cages covered by transparent and UV-A-opaque plastic films. Plant growth and insect fitness were monitored. Leaves were harvested for chemical analysis. Pepper plants responded directly to UV-A by producing shorter stems whilst UV-A did not affect the leaf area of either species. UV-A-treated peppers had higher content of secondary metabolites, soluble carbohydrates, free amino acids and proteins. Such changes in tissue chemistry indirectly promoted aphid performance. For eggplants, chlorophyll and carotenoid levels decreased with supplemental UVA but phenolics were not affected. Exposure to supplemental UV-A had a detrimental effect on whitefly development, fecundity and fertility presumably not mediated by plant cues, as compounds implied in pest nutrition were unaltered. Lastly, the efficacy of a wide range of Long Lasting Insecticide Treated Nets (LLITNs) was studied under laboratory and field conditions. This strategy aimed to prevent aphids and whiteflies to enter the greenhouse by determining the optimum mesh size (Chapter 7). This new approach is based on slow release deltamethrin- and bifenthrin-treated nets with large hole sizes that allow improved ventilation of greenhouses. All LLITNs produced high mortality of M. persicae and A. gossypii although their efficacy decreased over time with sun exposure. It was necessary a net with hole size of 0.29 mm2 to exclude B. tabaci under laboratory conditions. The feasibility of two selected nets was studied in the field under a high insect infestation pressure in the presence of CMV- and CABYV-infected cucumber plants. Besides, the compatibility of parasitoid A. colemani with bifenthrin-treated nets was studied in parallel field experiments. Both nets effectively blocked the invasion of aphids and reduced the incidence of both viruses, however they failed to exclude whiteflies. We found that our LLITNs were compatible with parasitoid A. colemani. As shown, the role of natural enemies has to be taken into account regarding the dispersal of insect vectors and subsequent spread of plant viruses. The additional benefits of novel physicochemical barriers, such as photoselective and insecticide-impregnated nets, need to be considered in Integrated Pest Management programs of vegetable crops grown under protected environments.

Cycling Dof Factors: Molecular and functional characterization of Arabidopsis thaliana AtCDF3 and tomato (Solanum lycopersicum L.) SlCDF3 in response to abiotic stress

El tomate (Solanum lycopersicum L.) es considerado uno de los cultivos hortícolas de mayor importancia económica en el territorio Español. Sin embargo, su producción está seriamente afectada por condiciones ambientales adversas como, salinidad, sequía y temperaturas extremas. Para resolver los problemas que se presentan en condiciones de estrés, se han empleado una serie de técnicas culturales que disminuyen sus efectos negativos, siendo de gran interés el desarrollo de variedades tolerantes. En este sentido la obtención y análisis de plantas transgénicas, ha supuesto un avance tecnológico, que ha facilitado el estudio y la evaluación de genes seleccionados en relación con la tolerancia al estrés. Estudios recientes han mostrado que el uso de genes reguladores como factores de transcripción (FTs) es una gran herramienta para obtener nuevas variedades de tomate con mayor tolerancia a estreses abióticos. Las proteínas DOF (DNA binding with One Finger) son una familia de FTs específica de plantas (Yangisawa, 2002), que están involucrados en procesos fisiológicos exclusivos de plantas como: asimilación del nitrógeno y fijación del carbono fotosintético, germinación de semilla, metabolismo secundario y respuesta al fotoperiodo pero su preciso rol en la tolerancia a estrés abiótico se desconoce en gran parte. El trabajo descrito en esta tesis tiene como objetivo estudiar genes reguladores tipo DOF para incrementar la tolerancia a estrés abiotico tanto en especies modelo como en tomate. En el primer capítulo de esta tesis se muestra la caracterización funcional del gen CDF3 de Arabidopsis, así como su papel en la respuesta a estrés abiótico y otros procesos del desarrollo. La expresión del gen AtCDF3 es altamente inducido por sequía, temperaturas extremas, salinidad y tratamientos con ácido abscísico (ABA). La línea de inserción T-DNA cdf3-1 es más sensible al estrés por sequía y bajas temperaturas, mientras que líneas transgénicas de Arabidopsis 35S::AtCDF3 aumentan la tolerancia al estrés por sequía, osmótico y bajas temperaturas en comparación con plantas wild-type (WT). Además, estas plantas presentan un incremento en la tasa fotosintética y apertura estomática. El gen AtCDF3 se localiza en el núcleo y que muestran una unión específica al ADN con diferente afinidad a secuencias diana y presentan diversas capacidades de activación transcripcional en ensayos de protoplastos de Arabidopsis. El dominio C-terminal de AtCDF3 es esencial para esta localización y su capacidad activación, la delección de este dominio reduce la tolerancia a sequía en plantas transgénicas 35S::AtCDF3. Análisis por microarray revelan que el AtCDF3 regula un set de genes involucrados en el metabolismo del carbono y nitrógeno. Nuestros resultados demuestran que el gen AtCDF3 juega un doble papel en la regulación de la respuesta a estrés por sequía y bajas temperaturas y en el control del tiempo de floración. En el segundo capítulo de este trabajo se lleva a cabo la identificación de 34 genes Dof en tomate que se pueden clasificar en base a homología de secuencia en cuatro grupos A-D, similares a los descritos en Arabidopsis. Dentro del grupo D se han identificado cinco genes DOF que presentan características similares a los Cycling Dof Factors (CDFs) de Arabidopsis. Estos genes son considerados ortólogos de Arabidopsis CDF1-5, y han sido nombrados como Solanum lycopersicum CDFs o SlCDFs. Los SlCDF1-5 son proteínas nucleares que muestran una unión específica al ADN con diferente afinidad a secuencias diana y presentan diversas capacidades de activación transcripcional in vivo. Análisis de expresión de los genes SlCDF1-5 muestran diferentes patrones de expresión durante el día y son inducidos de forma diferente en respuesta a estrés osmótico, salino, y de altas y bajas temperaturas. Plantas de Arabidopsis que sobre-expresan SlCDF1 y SlCDF3 muestran un incremento de la tolerancia a la sequía y salinidad. Además, de la expresión de varios genes de respuesta estrés como AtCOR15, AtRD29A y AtERD10, son expresados de forma diferente en estas líneas. La sobre-expresión de SlCDF3 en Arabidopsis promueve un retardo en el tiempo de floración a través de la modulación de la expresión de genes que controlan la floración como CONSTANS (CO) y FLOWERING LOCUS T (FT). En general, nuestros datos demuestran que los SlCDFs están asociados a funciones aun no descritas, relacionadas con la tolerancia a estrés abiótico y el control del tiempo de floración a través de la regulación de genes específicos y a un aumento de metabolitos particulares. ABSTRACT Tomato (Solanum lycopersicum L.) is one of the horticultural crops of major economic importance in the Spanish territory. However, its production is being affected by adverse environmental conditions such as salinity, drought and extreme temperatures. To resolve the problems triggered by stress conditions, a number of agricultural techniques that reduce the negative effects of stress are being frequently applied. However, the development of stress tolerant varieties is of a great interest. In this direction, the technological progress in obtaining and analysis of transgenic plants facilitated the study and evaluation of selected genes in relation to stress tolerance. Recent studies have shown that a use of regulatory genes such as transcription factors (TFs) is a great tool to obtain new tomato varieties with greater tolerance to abiotic stresses. The DOF (DNA binding with One Finger) proteins form a family of plant-specific TFs (Yangisawa, 2002) that are involved in the regulation of particular plant processes such as nitrogen assimilation, photosynthetic carbon fixation, seed germination, secondary metabolism and flowering time bur their precise roles in abiotic stress tolerance are largely unknown. The work described in this thesis aims at the study of the DOF type regulatory genes to increase tolerance to abiotic stress in both model species and the tomato. In the first chapter of this thesis, we present molecular characterization of the Arabidopsis CDF3 gene as well as its role in the response to abiotic stress and in other developmental processes. AtCDF3 is highly induced by drought, extreme temperatures, salt and abscisic acid (ABA) treatments. The cdf3-1 T-DNA insertion mutant was more sensitive to drought and low temperature stresses, whereas the AtCDF3 overexpression enhanced the tolerance of transgenic plants to drought, cold and osmotic stress comparing to the wild-type (WT) plants. In addition, these plants exhibit increased photosynthesis rates and stomatal aperture. AtCDF3 is localized in the nuclear region, displays specific binding to the canonical DNA target sequences and has a transcriptional activation activity in Arabidopsis protoplast assays. In addition, the C-terminal domain of AtCDF3 is essential for its localization and activation capabilities and the deletion of this domain significantly reduces the tolerance to drought in transgenic 35S::AtCDF3 overexpressing plants. Microarray analysis revealed that AtCDF3 regulated a set of genes involved in nitrogen and carbon metabolism. Our results demonstrate that AtCDF3 plays dual roles in regulating plant responses to drought and low temperature stress and in control of flowering time in vegetative tissues. In the second chapter this work, we carried out to identification of 34 tomato DOF genes that were classified by sequence similarity into four groups A-D, similar to the situation in Arabidopsis. In the D group we have identified five DOF genes that show similar characteristics to the Cycling Dof Factors (CDFs) of Arabidopsis. These genes were considered orthologous to the Arabidopsis CDF1 - 5 and were named Solanum lycopersicum CDFs or SlCDFs. SlCDF1-5 are nuclear proteins that display specific binding to canonical DNA target sequences and have transcriptional activation capacities in vivo. Expression analysis of SlCDF1-5 genes showed distinct diurnal expression patterns and were differentially induced in response to osmotic, salt and low and high temperature stresses. Arabidopsis plants overexpressing SlCDF1 and SlCDF3 showed increased drought and salt tolerance. In addition, various stress-responsive genes, such as AtCOR15, AtRD29A and AtERD10, were expressed differently in these lines. The overexpression of SlCDF3 in Arabidopsis also results in the late flowering phenotype through the modulation of the expression of flowering control genes such CONSTANS (CO) and FLOWERING LOCUS T (FT). Overall, our data connet SlCDFs to undescribed functions related to abiotic stress tolerance and flowering time through the regulation of specific target genes and an increase in particular metabolites.

Self-organizing cultured neural networks: Image analysis techniques for longitudinal tracking and modeling of the underlying network structure

Esta tesis estudia la evolución estructural de conjuntos de neuronas como la capacidad de auto-organización desde conjuntos de neuronas separadas hasta que forman una red (clusterizada) compleja. Esta tesis contribuye con el diseño e implementación de un algoritmo no supervisado de segmentación basado en grafos con un coste computacional muy bajo. Este algoritmo proporciona de forma automática la estructura completa de la red a partir de imágenes de cultivos neuronales tomadas con microscopios de fase con una resolución muy alta. La estructura de la red es representada mediante un objeto matemático (matriz) cuyos nodos representan a las neuronas o grupos de neuronas y los enlaces son las conexiones reconstruidas entre ellos. Este algoritmo extrae también otras medidas morfológicas importantes que caracterizan a las neuronas y a las neuritas. A diferencia de otros algoritmos hasta el momento, que necesitan de fluorescencia y técnicas inmunocitoquímicas, el algoritmo propuesto permite el estudio longitudinal de forma no invasiva posibilitando el estudio durante la formación de un cultivo. Además, esta tesis, estudia de forma sistemática un grupo de variables topológicas que garantizan la posibilidad de cuantificar e investigar la progresión de las características principales durante el proceso de auto-organización del cultivo. Nuestros resultados muestran la existencia de un estado concreto correspondiente a redes con configuracin small-world y la emergencia de propiedades a micro- y meso-escala de la estructura de la red. Finalmente, identificamos los procesos físicos principales que guían las transformaciones morfológicas de los cultivos y proponemos un modelo de crecimiento de red que reproduce el comportamiento cuantitativamente de las observaciones experimentales. ABSTRACT The thesis analyzes the morphological evolution of assemblies of living neurons, as they self-organize from collections of separated cells into elaborated, clustered, networks. In particular, it contributes with the design and implementation of a graph-based unsupervised segmentation algorithm, having an associated very low computational cost. The processing automatically retrieves the whole network structure from large scale phase-contrast images taken at high resolution throughout the entire life of a cultured neuronal network. The network structure is represented by a mathematical object (a matrix) in which nodes are identified neurons or neurons clusters, and links are the reconstructed connections between them. The algorithm is also able to extract any other relevant morphological information characterizing neurons and neurites. More importantly, and at variance with other segmentation methods that require fluorescence imaging from immunocyto- chemistry techniques, our measures are non invasive and entitle us to carry out a fully longitudinal analysis during the maturation of a single culture. In turn, a systematic statistical analysis of a group of topological observables grants us the possibility of quantifying and tracking the progression of the main networks characteristics during the self-organization process of the culture. Our results point to the existence of a particular state corresponding to a small-world network configuration, in which several relevant graphs micro- and meso-scale properties emerge. Finally, we identify the main physical processes taking place during the cultures morphological transformations, and embed them into a simplified growth model that quantitatively reproduces the overall set of experimental observations.

Self-organizing cultured neural networks: Image analysis techniques for longitudinal tracking and modeling of the underlying network structure (versión actualizada)

Esta tesis estudia la evolución estructural de conjuntos de neuronas como la capacidad de auto-organización desde conjuntos de neuronas separadas hasta que forman una red (clusterizada) compleja. Esta tesis contribuye con el diseño e implementación de un algoritmo no supervisado de segmentación basado en grafos con un coste computacional muy bajo. Este algoritmo proporciona de forma automática la estructura completa de la red a partir de imágenes de cultivos neuronales tomadas con microscopios de fase con una resolución muy alta. La estructura de la red es representada mediante un objeto matemático (matriz) cuyos nodos representan a las neuronas o grupos de neuronas y los enlaces son las conexiones reconstruidas entre ellos. Este algoritmo extrae también otras medidas morfológicas importantes que caracterizan a las neuronas y a las neuritas. A diferencia de otros algoritmos hasta el momento, que necesitan de fluorescencia y técnicas inmunocitoquímicas, el algoritmo propuesto permite el estudio longitudinal de forma no invasiva posibilitando el estudio durante la formación de un cultivo. Además, esta tesis, estudia de forma sistemática un grupo de variables topológicas que garantizan la posibilidad de cuantificar e investigar la progresión de las características principales durante el proceso de auto-organización del cultivo. Nuestros resultados muestran la existencia de un estado concreto correspondiente a redes con configuracin small-world y la emergencia de propiedades a micro- y meso-escala de la estructura de la red. Finalmente, identificamos los procesos físicos principales que guían las transformaciones morfológicas de los cultivos y proponemos un modelo de crecimiento de red que reproduce el comportamiento cuantitativamente de las observaciones experimentales. ABSTRACT The thesis analyzes the morphological evolution of assemblies of living neurons, as they self-organize from collections of separated cells into elaborated, clustered, networks. In particular, it contributes with the design and implementation of a graph-based unsupervised segmentation algorithm, having an associated very low computational cost. The processing automatically retrieves the whole network structure from large scale phase-contrast images taken at high resolution throughout the entire life of a cultured neuronal network. The network structure is represented by a mathematical object (a matrix) in which nodes are identified neurons or neurons clusters, and links are the reconstructed connections between them. The algorithm is also able to extract any other relevant morphological information characterizing neurons and neurites. More importantly, and at variance with other segmentation methods that require fluorescence imaging from immunocyto- chemistry techniques, our measures are non invasive and entitle us to carry out a fully longitudinal analysis during the maturation of a single culture. In turn, a systematic statistical analysis of a group of topological observables grants us the possibility of quantifying and tracking the progression of the main networks characteristics during the self-organization process of the culture. Our results point to the existence of a particular state corresponding to a small-world network configuration, in which several relevant graphs micro- and meso-scale properties emerge. Finally, we identify the main physical processes taking place during the cultures morphological transformations, and embed them into a simplified growth model that quantitatively reproduces the overall set of experimental observations.

Evaluación del uso y de productividad del agua de la Comunidad de Regantes “Río Adaja”

El agua de riego en España se ha reducido del 80 % al 70% tras la rehabilitación de los sistemas tradicionales de riego y el incremento de riegos a presión. La política española ha favorecido la creación de nuevos regadíos con fines sociales, para asentar a la población rural en zonas con disponibilidad de recursos hídricos. Este contexto es aplicable a la Comunidad de Regantes “Rio Adaja” (CCRR), que comenzó a funcionar en 2010 por lo que se la ha elegido para evaluar el uso y productividad del agua y manejo del riego en CCRR modernizadas de la cuenca del Duero. El estudio del manejo del riego se realizó con evaluaciones de campo, el primer año de funcionamiento, en una muestra de sistemas de riego (pivotes centrales, ramales de avance frontal, cobertura total) representativa de los sistemas predominantes en la CCRR. Además, se analizó la carta de riego propuesta por el fabricante de los pivotes centrales, considerando una distribución de caudal continua a lo largo del ramal, y se propuso una nueva carta con emisores de riego que mejoraban la uniformidad de aplicación del agua. El uso del agua en la CCRR se evaluó considerando tanto los indicadores de eficiencia del riego: suministro relativo de riego (anual relative irrigation supply, ARIS), suministro relativo del agua (anual relativewater supply, ARWS), suministro relativo de precipitación (rainfall relative supply, RRS) como los de productividad: productividad del agua (water productivity, WP) productividad del agua de riego (irrigation water productivity, IWP) y productividad de la evapotranspiración (evapotranspiration water productivity, ETWP). Primero, se determinaron: las necesidades hídricas de los cultivos para mantener un contenido de humedad óptimo en su zona radical, el coeficiente dual del cultivo, el agua disponible total (ADP) y agua fácilmente aprovechable (AFA). Después, se estimaron las necesidades hídricas de los cultivos considerando tres años tipo: húmedo, normal y seco correspondientes a la probabilidad de disponibilidad de la precipitación del 20, 50 y 80%, respectivamente. Así mismo, se realizó una encuesta a los regantes de la CCRR para conocer la dosis de riego y rendimiento anual de los cultivos principales durante sus tres años de funcionamiento: 2010-2011, 2011-2012 y 2012-2013.Finalmente, se simuló el efecto del riego y su manejo en la producción de los cultivos y en la productividad del agua. Además, el modelo de simulación AQUACROP (Geerts et al., 2010) se ha utilizado para estudiar la mejora del uso del agua de los cultivos de la CCRR. Dado que este modelo requiere de calibración específica para cada cultivo y cada zona y dado que, de todos los cultivos de la CCRR, sólo el girasol cumplía el requisito, este cultivo fue elegido para estudiar si la estrategia de riego deficitario mejoraría el uso del agua. Los resultados obtenidos indican que el 90% de los sistemas de riego evaluados distribuye el agua con una uniformidad adecuada (CUC≥75%). La simulación de la distribución del agua con las cartas de riego propuestas por el fabricante en pivotes centrales resultó en coeficientes CUC< 75% y sus valores mejoraban al eliminar el aspersor distal. La uniformidad del riego mejoraría si se trabajase con la carta de riego propuesta y que se compone por emisores de riego seleccionados en este estudio. En la mayoría de los cultivos, se aplicó riegos deficitarios (ARIS < 1 en los dos primeros años de funcionamiento de la CCRR y riegos excedentarios (ARIS > 1) el tercer año siendo significativas las diferencias observadas. El manejo del riego fue bueno (0,9 ≤ ARWS ≤1,2) en la mayoría de los cultivos. Así mismo, los indicadores de productividad del agua (WP e IWP (€.m-3)) varió entre cultivos y años estudiados y, destacan los valores observados en: cebolla, patata, zanahoria y cebada. En general, la productividad del agua en los riegos deficitarios fue mayor observándose además, que los índices de productividad mayores correspondieron al año con precipitación mayor aunque, las diferencias entre sus valores medios no fueron significativas en las tres campañas de riego estudiadas. Los resultados apuntan a que la metodología del balance hídrico y las herramientas presentadas en este trabajo (uniformidad de distribución de agua, indicadores de eficiencia del uso de agua y de su productividad) son adecuadas para estudiar el manejo del agua en CCRR. En concreto, la uniformidad en la aplicación del agua de la CCRR mejoraría seleccionando emisores de riego que proporcionen una mayor uniformidad de distribución del agua, lo que conllevaría a cambiar el diámetro de la boquilla de los emisores y/o eliminar el aspersor distal. Así mismo, puede ser de interés adoptar estrategias de riego deficitario para incrementar la productividad en el uso del agua, y las rentas de los regantes, para lo cual se propone utilizar un patrón de cultivos de referencia. Finalmente, el riego deficitario puede ser una estrategia para mejorar la eficiencia y productividad en el uso del agua de la CCRR siempre que lleve asociado un manejo del riego adecuado que resulta, relativamente, más fácil cuando se dispone de sistemas de riego con una uniformidad de aplicación alta. Sin embargo su aplicación no sería aconsejable en los cultivos de remolacha azucarera, regado con sistemas de riego con un coeficiente de uniformidad de Christiansen CUC < 75%, y maíz, regado con sistemas de riego con un coeficiente de uniformidad de Christiansen CUC < 65%. ABSTRACT The irrigation scheme modernization and the increase of sprinkler irrigation area have reduced the irrigation water use from 80 to 70%. The national irrigation policy favored the creation of new irrigation schemes with the purpose to settle the rural population in areas with availability in water resources. Within this context, the irrigation district “Río Adaja” (CCRR) started in 2010 so, it has been chosen as a case study to evaluate the water use and the irrigation management in a modernized CCRR. Several field evaluations were carried out during the first operation year, in a sample of irrigation systems (center pivot, moving lateral and solid set) selected among all the systems in the CCRR. Likewise, the manufacturer irrigation chart for the center pivot systems has been considered and the pressure and discharge distribution along the pivot have been estimated, assuming a continuous flow along the pipe. Then; the sprinkler nozzles were selected order to increase the uniformity on water application. The water use in the CCRR has been assessed by considering the water use efficiency indicators: annual relative irrigation supply (ARIS), annual relative water supply (ARWS), relative rainfall supply (RRS) and also the productivity indicators: water productivity (WP), irrigation water productivity (IWP) and evapotranspiration water productivity (ETWP). On the one hand, it has been determined the crop water requirement (to maintain the optimal soil water content in the rooting zone), the dual crop coefficient, the total available water and the readily available water. The crop water requirement was estimated by considering the typical wet, normal and dry years which correspond to the probability of effective precipitation exceedance of 20, 50 and 80%, respectively. On the other hand, the irrigation depth and crop yield by irrigation campaign have been considered for the main crops in the area. This information was obtained from a farmer’s survey in 2010-2011, 2011-2012 and 2012-2013. For sunflower, the irrigation effect and its management on the crop yield and water productivity have been simulated. Also a deficit irrigation strategy, which improves the water resources, has been determined by means of AQUACROP (FAO). The results showed that 90% of the evaluated irrigation systems have adequate irrigation water application uniformity (CUC ≥ 75%). The CUC values in center pivots, which were calculated using the manufacturer irrigation chart, are below < 75% . However, these values would increase with the change of emitter nozzle to the proposed nozzles selection. The results on water use showed a deficit irrigation management (ARIS < 1), in most of crops during the first two operation years, and an excess in irrigation for the third year (ARIS > 1) although non-significant difference was observed. In most cases, the management of irrigation is adequate (0,9≤ ARWS≤ 1,2) although there are differences among crops. Likewise, the productivity indicators (WP and IWP (€.m-3)) varied among crops and with irrigation events. The highest values corresponded to onion, potato, carrot and barley. The values for deficit irrigation were the highest and the productivity indicators increased the year with the highest effective precipitation. Nevertheless, the differences between the average values of these indicators by irrigation campaign were non-significant. This study highlights that the soil water balance methodology and other tools used in the methodology are adequate to study the use and productivity of water in the irrigation district. In fact, the water use in this CCRR can be improved if the irrigation systems were designed with higher water distribution uniformity what would require the change of sprinkler nozzles and/or eliminate the end gun. Likewise, it is advisable to set up deficit irrigation strategies to increase the water productivity taking into account certain limits on water application uniformities. In this respect, a reference cropping pattern has been proposed and the limits for water uniformity have been calculated for several crops.

Compromisos en la adaptación a distintos huéspedes en virus de plantas

Los patógenos han desarrollado estrategias para sobrevivir en su entorno, infectar a sus huéspedes, multiplicarse dentro de estos y posteriormente transmitirse a otros huéspedes. Todos estos componentes hacen parte de la eficacia biológica de los patógenos, y les permiten ser los causantes de enfermedades infecciosas tanto en hombres y animales, como en plantas. El proceso de infección produce efectos negativos en la eficacia biológica del huésped y la gravedad de los efectos, dependerá de la virulencia del patógeno. Por su parte, el huésped ha desarrollado mecanismos de respuesta en contra del patógeno, tales como la resistencia, por la que reduce la multiplicación del patógeno, o la tolerancia, por la que disminuye el efecto negativo de la infección. Estas respuestas del huésped a la infección producen efectos negativos en la eficacia biológica del patógeno, actuando como una presión selectiva sobre su población. Si la presión selectiva sobre el patógeno varía según el huésped, se predice que un mismo patógeno no podrá aumentar su eficacia biológica en distintos huéspedes y estará más adaptado a un huésped y menos a otro, disminuyendo su gama de huéspedes. Esto supone que la adaptación de un patógeno a distintos huéspedes estará a menudo dificultada por compromisos (trade-off) en diferentes componentes de la eficacia biológica del patógeno. Hasta el momento, la evidencia de compromisos de la adaptación del patógeno a distintos huéspedes no es muy abundante, en lo que se respecta a los virus de plantas. En las últimas décadas, se ha descrito un aumento en la incidencia de virus nuevos o previamente descritos que producen enfermedades infecciosas con mayor gravedad y/o diferente patogenicidad, como la infección de huéspedes previamente resistentes. Esto se conoce como la emergencia de enfermedades infecciosas y está causada por patógenos emergentes, que proceden de un huésped reservorio donde se encuentran adaptados. Los huéspedes que actúan como reservorios pueden ser plantas silvestres, que a menudo presentan pocos síntomas o muy leves a pesar de estar infectados con diferentes virus, y asimismo se encuentran en ecosistemas con ninguna o poca intervención humana. El estudio de los factores ecológicos y biológicos que actúan en el proceso de la emergencia de enfermedades infecciosas, ayudará a entender sus causas para crear estrategias de prevención y control. Los virus son los principales patógenos causales de la emergencia de enfermedades infecciosas en humanos, animales y plantas y un buen modelo para entender los procesos de la emergencia. Asimismo, las plantas a diferencia de los animales, son huéspedes fáciles de manipular y los virus que las afectan, más seguros para el trabajo en laboratorio que los virus de humanos y animales, otros modelos también usados en la investigación. Por lo tanto, la interacción virus – planta es un buen modelo experimental para el estudio de la emergencia de enfermedades infecciosas. El estudio de la emergencia de virus en plantas tiene también un interés particular, debido a que los virus pueden ocasionar pérdidas económicas en los cultivos agrícolas y poner en riesgo la durabilidad de la resistencia de plantas mejoradas, lo que supone un riesgo en la seguridad alimentaria con impactos importantes en la sociedad, comparables con las enfermedades infecciosas de humanos y animales domésticos. Para que un virus se convierta en un patógeno emergente debe primero saltar desde su huésped reservorio a un nuevo huésped, segundo adaptarse al nuevo huésped hasta que la infección dentro de la población de éste se vuelva independiente del reservorio y finalmente debe cambiar su epidemiología. En este estudio, se escogió la emergencia del virus del mosaico del pepino dulce (PepMV) en el tomate, como modelo experimental para estudiar la emergencia de un virus en una nueva especie de huésped, así como las infecciones de distintos genotipos del virus del moteado atenuado del pimiento (PMMoV) en pimiento, para estudiar la emergencia de un virus que aumenta su patogenicidad en un huésped previamente resistente. El estudio de ambos patosistemas nos permitió ampliar el conocimiento sobre los factores ecológicos y evolutivos en las dos primeras fases de la emergencia de enfermedades virales en plantas. El PepMV es un patógeno emergente en cultivos de tomate (Solanum lycopersicum) a nivel mundial, que se describió primero en 1980 infectando pepino dulce (Solanum muricatum L.) en Perú, y casi una década después causando una epidemia en cultivos de tomate en Holanda. La introducción a Europa posiblemente fue a través de semillas infectadas de tomate procedentes de Perú, y desde entonces se han descrito nuevos aislados que se agrupan en cuatro cepas (EU, LP, CH2, US1) que infectan a tomate. Sin embargo, el proceso de su emergencia desde pepino dulce hasta tomate es un interrogante de gran interés, porque es uno de los virus emergentes más recientes y de gran importancia económica. Para la emergencia de PepMV en tomate, se recolectaron muestras de tomate silvestre procedentes del sur de Perú, se analizó la presencia y diversidad de aislados de PepMV y se caracterizaron tanto biológicamente (gama de huéspedes), como genéticamente (secuencias genomicas). Se han descrito en diferentes regiones del mundo aislados de PMMoV que han adquirido la capacidad de infectar variedades previamente resistentes de pimiento (Capsicum spp), es decir, un típico caso de emergencia de virus que implica la ampliación de su gama de huéspedes y un aumento de patogenicidad. Esto tiene gran interés, ya que compromete el uso de variedades resistentes obtenidas por mejora genética, que es la forma de control de virus más eficaz que existe. Para estudiar la emergencia de genotipos altamente patogénicos de PMMoV, se analizaron clones biológicos de PMMoV procedentes de aislados de campo cuya patogenicidad era conocida (P1,2) y por mutagénesis se les aumentó la patogenicidad (P1,2,3 y P1,2,3,4), introduciendo las mutaciones descritas como responsables de estos fenotipos. Se analizó si el aumento de la patogenicidad conlleva un compromiso en la eficacia biológica de los genotipos de PMMoV. Para ello se evaluaron diferentes componentes de la eficacia biológica del virus en diferentes huéspedes con distintos alelos de resistencia. Los resultados de esta tesis demuestran: i). El potencial de las plantas silvestres como reservorios de virus emergentes, en este caso tomates silvestres del sur de Perú, así como la existencia en estas plantas de aislados de PepMV de una nueva cepa no descrita que llamamos PES. ii) El aumento de la gama de huéspedes no es una condición estricta para la emergencia de los virus de plantas. iii) La adaptación es el mecanismo más probable en la emergencia de PepMV en tomate cultivado. iv) El aumento de la patogenicidad tiene un efecto pleiotrópico en distintos componentes de la eficacia biológica, así mismo el signo y magnitud de este efecto dependerá del genotipo del virus, del huésped y de la interacción de estos factores. ABSTRACT host Pathogens have evolved strategies to survive in their environment, infecting their hosts, multiplying inside them and being transmitted to other hosts. All of these components form part of the pathogen fitness, and allow them to be the cause of infectious diseases in humans, animals, and plants. The infection process produces negative effects on the host fitness and the effects severity will depend on the pathogen virulence. On the other hand, hosts have developed response mechanisms against pathogens such as resistance, which reduces the growth of pathogens, or tolerance, which decreases the negative effects of infection. T he se responses of s to infection cause negative effects on the pathogen fitness, acting as a selective pressure on its population. If the selective pressures on pathogens va ry according to the host s , probably one pathogen cannot increase its fitness in different hosts and will be more adapted to one host and less to another, decreasing its host range. This means that the adaptation of one pathogen to different hosts , will be often limited by different trade - off components of biological effectiveness of pathogen. Nowadays , trade - off evidence of pathogen adaptation to different hosts is not extensive, in relation with plant viruses. In last decades, an increase in the incidence of new or previously detected viruses has been described, causing infectious diseases with increased severity and/or different pathogenicity, such as the hosts infection previously resistants. This is known as the emergence of infectious diseases and is caused by emerging pathogens that come from a reservoir host where they are adapted. The hosts which act as reservoirs can be wild plants, that often have few symptoms or very mild , despite of being infected with different viruses, and being found in ecosystems with little or any human intervention. The study of ecological and biological factors , acting in the process of the infectious diseases emergence will help to understand its causes to create strategies for its prevention and control. Viruses are the main causative pathogens of the infectious diseases emergence in humans, animals and plants, and a good model to understand the emergency processes. Likewise, plants in contrast to animals are easy host to handle and viruses that affect them, safer for laboratory work than viruses of humans and animals, another models used in research. Therefore, the interaction plant-virus is a good experimental model for the study of the infectious diseases emergence. The study of virus emergence in plants also has a particular interest, because the viruses can cause economic losses in agricultural crops and threaten the resistance durability of improved plants, it suppose a risk for food security with significant impacts on society, comparable with infectious diseases of humans and domestic animals. To become an emerging pathogen, a virus must jump first from its reservoir host to a new host, then adapt to a new host until the infection within the population becomes independent from the reservoir, and finally must change its epidemiology. In this study, the emergence of pepino mosaic virus (PepMV) in tomato, was selected as experimental model to study the emergence of a virus in a new host specie, as well as the infections of different genotypes of pepper mild mottle virus (PMMoV) in pepper, to study the emergence of a virus that increases its pathogenicity in a previously resistant host. The study of both Pathosystems increased our knowledge about the ecological and evolutionary factors in the two first phases of the emergence of viral diseases in plants. The PepMV is an emerging pathogen in tomato (Solanum lycopersicum L.) in the world, which was first described in 1980 by infecting pepino (Solanum muricatum L.) in Peru, and almost after a decade caused an epidemic in tomato crops in Netherlands. The introduction to Europe was possibly through infected tomato seeds from Peru, and from then have been described new isolates that are grouped in four strains (EU, LP, CH2, US1) that infect tomato. However, the process of its emergence from pepino up tomato is a very interesting question, because it is one of the newest emerging viruses and economically important. For the PepMV emergence in tomato, wild tomato samples from southern Peru were collected, and the presence and diversity of PepMV isolates were analyzed and characterized at biological (host range) and genetics (genomic sequences) levels. Isolates from PMMoV have been described in different world regions which have acquired the ability to infect pepper varieties that were previously resistants (Capsicum spp), it means, a typical case of virus emergence which involves the host range extension and an increased pathogenicity. This is of great interest due to involve the use of resistant varieties obtained by breeding, which is the most effective way to control virus. To study the emergence of highly pathogenic genotypes of PMMoV, biological clones from field isolates whose pathogenicity was known were analyzed (P1,2) and by mutagenesis we increased its pathogenicity (P1,2,3 and P1,2, 3,4), introducing the mutations described as responsible for these phenotypes. We analyzed whether the increased pathogenicity involves a trade-off in fitness of PMMoV genotypes. For this aim, different components of virus fitness in different hosts with several resistance alleles were evaluated. The results of this thesis show: i). The potential of wild plants as reservoirs of emerging viruses, in this case wild tomatoes in southern Peru, and the existence in these plants of PepMV isolates of a new undescribed strain that we call PES. ii) The host range expansion is not a strict condition for the plant virus emergence. iii) The adaptation is the most likely mechanism in the PepMV emergence in cultivated tomato. iv) The increased pathogenicity has a pleiotropic effect on several fitness components, besides the sign and magnitude of this effect depends on the virus genotype, the host and the interaction of both.

Biofortificación de cereales con micronutrientes

La malnutrición en humanos por falta de micronutrientes afecta a más de 3.000 millones de personas en todo el mundo. Las deficiencias más frecuentes son de zinc (Zn), hierro (Fe), vitamina A y yodo (I) y son particularmente evidentes en las mujeres y los niños de países en desarrollo. Las carencias de Zn y Fe afectan a alrededor de un tercio de la población del mundo y se están extendiendo en los países desarrollados, especialmente en Europa. La biofortificación agronómica, incrementando la concentración de micronutriente Zn en la planta, es un medio para evitar la deficiencia de Zn en animales y humanos y consiste en aumentar el contenido de micronutrientes en las plantas mediante la fertilización de los cultivos.

Contribución al estudio y clasificación de las funciones de implicación borrosas

La tesis doctoral "Contribución al estudio y clasificación de las funciones de implicación borrosas" constituye a la vez una revisión y un conjunto de aportaciones a la modelización de enunciados condicionales, o enunciados del tipo "Si P, entonces Q". Dentro del marco de la lógica borrosa, tradicionalmente, se considera que la mayor parte de las funciones de implicación, que modelizan los enunciados condicionales, deberían generalizar la implicación material booleana. En esta memoria se apoya el argumento de que la implicación material booleana no es siempre el modelo más adecuado para la modelización de los enunciados condicionales, por lo que se definen y estudian clases o grupos de funciones de implicación que se adecúen a las necesidades de cada caso de aplicación. Así pues, tras un capítulo introductorio, en el capítulo 2 se plantean clases de funciones de implicación que sirvan de apoyo en la definición de aplicaciones de carácter borroso más allá de las funciones de implicación borrosas derivadas de la implicación material booleana. Así se llega a la conclusión de que son necesarios, por lo menos, cuatro clases de funciones de implicación de las que se estudian sus propiedades. En el tercer capítulo se aborda el estudio de la verificación tanto de la regla del Modus Ponens como de la regla del Modus Tollens por parte de las funciones de implicación pertenecientes a cada una de las clases planteadas. Por último, en el capítulo 4 se presenta un estudio de los operadores citados como implicaciones borrosas en la literatura bajo el enfoque presentado en esta memoria, ampliándose estudios realizados por otros autores. ---ABSTRACT--- The Ph.D. Thesis "Contribution al estudio y clasificacion de las funciones de implication borrosas" is a revision and a set of contributions to the conditional statements or statements of the type "If P then Q" modelization. It has been usually considered in Fuzzy Logic, that implication functions that model conditional statements should generalize the boolean material implication. The present memory supports the argument that boolean material implication is not always the most appropiated model for conditional statements, so implication function groups or clusters are denned in order to fit in each application. Therefore, after an introductory chapter, in chapter 2 groups or clusters of implication functions are given in order to bear out fuzzy applications definitions further on from implication functions coming from boolean material implication. Then, it has been reached the conclusion that, at least, four types of implication functions are needed so their properties are studied.

Uso de microalgas inmovilizadas en biofilm para tratamiento de aguas residuales

Los sistemas desarrollados para el cultivo de las microalgas pueden clasificarse en dos grandes grupos, en función de la forma de desarrollo de las poblaciones algales: cultivos en suspensión y cultivos inmovilizados. A su vez, ambos sistemas de cultivo pueden ser abiertos o cerrados, según que las microalgas y el medio de cultivo estén en contacto directo con el aire atmosférico o no. Los cultivos de microalgas inmovilizadas en una superficie que actúa de soporte, forman un 'biofilm' contínuo sobre dicha superficie formando un ecosistema específico de agregación denominado 'biofilm', normalmente constituido por varias especies, incluyendo bacterias, aunque puede haber alguna predominante. Las principales ventajas que tienen los fotobior reactores de biofilm sobre los de células suspendidas son: a) Facilidad para la cosecha, b) Mayor concentración por unidad de volumen de medio, c) reducción o ausencia de células en el efluente, d) reducción en la necesidad de energía y e) mayor eficiencia en el empleo de la radiación recibida por unidad de superficie de suelo. El uso de fotobiorreactores abiertos de microalgas inmovilizadas es relativamente reciente, dirigiéndose los desarrollos en dos direcciones: una que considera el soporte del biofilm fijo, siendo el medio acuoso el que se mueve sobre él y otra que considera el soporte del biofilm móvil en un medio acuático estático o de lento movimiento. Al primero de los tipos pertenece el fotobiorreactor laminar desarrollado por el Grupo de Agroenergética de la UPM (GA-UPM) (patente ES 2 347 515 B2 ), especialmente concebido para la captación del CO2 procedente de emisiones y para la producción de biomasa algal. Otro desarrollo del GA-UPM perteneciente al segundo tipo ha sido un fotobiorreactor de biodisco de eje flotante, que elimina los inconvenientes de los contactores de biodisco tradicionales, ya que no necesita cojinetes de apoyo ni eje rígido y al estar dotado de un propulsor hidrobárico consume muy poca energía, sin limitación técnica en cuanto a la longitud total del eje. (Modelo de Utilidad ES 1 099 707 U). Este fotobiorreactor está especialmente diseñado para la eliminación de nutrientes (N y P) de los efluentes del secundario de las EDAR y también puede utilizarse como un contactor tradicional de biodiscos para eliminación de la materia orgánica de las aguas residuales con indudables ventajas de reducción de costes de instalación y funcionamiento.

Development of somatic embryogenesis for cloning and conservation of mature holm oak trees (Quercus ilex L.)

La encina (Quercus ilex L.) es una de las especies forestales mediterráneas más importantes. Constituye gran parte del estrato arbóreo de dehesas o montados, produce bellota como alimento del ganado y establece simbiosis con hongos micorrizógenos de gran valor económico. La encina está considerada como una especie recalcitrante en términos de conservación de semillas y capacidad morfogénica, lo que dificulta los programas de conservación de recursos genéticos y la mejora de la especie. La propagación vegetativa es una potente herramienta de los programas de mejora, por lo que es preciso desarrollar protocolos de regeneración somática en encina. La embriogénesis somática está considerada como la modalidad más adecuada de regeneración basada en técnicas de cultivo de tejidos vegetales utilizada en biotecnología forestal. Este trabajo se centra en el estudio de determinados aspectos de la embriogénesis somática para la regeneración clonal de encinas adultas. La memoria de esta tesis se ha dividido en capítulos que se corresponden con diferentes aspectos del sistema embriogénico. La embriogénesis somática se indujo en tegumentos maternos de óvulos en desarrollo procedentes de bellotas inmaduras de encinas adultas. A pesar de las bajas frecuencias de inducción, las líneas embriogénicas generadas se amplificaron mediante embriogénesis secundaria observándose cierta pérdida de la capacidad de diferenciación con el tiempo. Tanto el genotipo como la formulación del medio de cultivo influyeron en la respuesta embriogénica, concluyendo que la formulación de macronutrientes de Schenk y Hildebrant del medio sin reguladores de crecimiento fue la combinación más efectiva en la inducción. Los resultados sugirieron la existencia de una ventana en el desarrollo del óvulo más sensible a la inducción. El genotipo in[luyó en la capacidad proliferativa de los cultivos y en la conversión de los embriones somáticos, que se incrementó suplementando el medio con ácido indol-3-butírico y 6-benciladenina. El cultivo en medio líquido de líneas embriogénicas en condiciones de inmersión transitoria incrementó el crecimiento, dependiendo del genotipo, con respecto al cultivo en medio semisólido. Sin embargo, no mejoró la capacidad de diferenciar embriones cotiledonares aislados. Se estableció un protocolo de inicio y mantenimiento de cultivos en suspensión para varias líneas embriogénicas mediante inoculación en alta densidad de agregados embrionarios procedentes del medio semisólido. Para evitar la pérdida de vigor y la capacidad morfogénica debida al cultivo prolongado se desarrolló un protocolo de crioconservación de líneas embriogénicas mediante vitrificación. Al determinar la influencia de los agentes crioprotectores antes y después de su inmersión en nitrógeno líquido se concluyó que las respuestas de capacidad de crecimiento y de diferenciación del material embriogénico son independientes, además de estar bajo influencia del genotipo y el tipo de material crioconservado. La combinación de sacarosa y PVS2 previa a la inmersión en nitrógeno líquido proporcionó la mayor tasa de recuperación. Cuando las líneas fueron crioconservadas 30 días la capacidad de diferenciación se perdió en todas ellas. El análisis de SSR detectó variación somaclonal en el material crioconservado a corto plazo. SSR y RAPD mostraron importantes diferencias genéticas entre los árboles donantes y el material embriogénico que dependieron del genotipo. El grado de detección dependió del marcador empleado. Ambos marcadores revelaron baja inestabilidad intraclonal. Los RAPD revelaron variación genética intra-individuo en las encinas donantes. Se discuten la variación genética pre-existente en encina, su aparición durante las primeras fases de la inducción de embriogénesis, y la presencia de tejidos provenientes de la fertilización en el explanto materno. Esto hace preciso definir la identidad genética del material donante y acometer ensayos de detección precoz de variación somaclonal. ABSTRACT Holm oak (Quercus ilex L.) is one of the most important Mediterranean forest species. It conforms the tree layer of dehesas or montados, it produces acorns to feed the livestock and it establishes symbiosis with profitable mycorrhizal fungi. Holm oak is considered as recalcitrant species in terms of seed conservation and morphogenic capacities, which complicates the development of genetic conservation and improvement programs. Vegetative propagation is one of the mightiest tools for breeding programs therefore; developing protocols for clonal regeneration of holm oak is essential. Somatic embryogenesis is considered the best tissue culture-based way of plant regeneration in forest biotechnology. The present study is focused on the study of certain aspects of somatic embryogenesis for clonal regeneration of mature holm oak. This thesis manuscript is divided into several chapters that match with different aspects of the embryogenic system. Somatic embryogenesis induction was achieved on maternal teguments of developing ovules from immature acorns of adult holm oak trees. Despite the low induction frequencies, the generated embryogenic lines were amplified by secondary embryogenesis. A decline in the differentiation capacity over time was also observed. It was concluded that both genotype and culture media formulation influenced the embryogenic response, being the Schenk and Hildebrandt´s macronutrients formulation from culture medium and the lack of plant growth regulators the most effective combination for the induction of the embryogenic response. It has been suggested the existence of a developmental window in which ovules are prone to induction. Genotype influenced the proliferation capacity and the plant conversion of somatic embryos, which was also favoured by the presence of indol-3-butyric acid and 6-bencyladenine. The use of temporary immersion systems as proliferation in liquid culture of the embryogenic lines increased the growth depending on genotype, when compared to semisolid cultures. However, it did not improve the differentiation of single cotyledonary embryos. A protocol for the initiation and maintenance of embryogenic suspension cultures was established for several embryogenic lines with highly dense inoculi of embryogenic clusters from proliferating semisolid cultures. In order to avoid the loss of vigour and morphogenic ability of embryogenic lines due to prolonged cultures, a cryopreservation protocol for embryogenic lines of holm oak has been developed. During the determination of the influence of cryoprotective agents on the growth and differentiation capacities before and after liquid nitrogen immersion, it was concluded that both responses were independent from each other and also under the influence of genotype and the type of cryopreserved material. The combination of sucrose and PVS2 prior liquid nitrogen immersion provided higher recovery rates. When the same embryogenic lines were cryopreserved for 30 days, none was able to differentiate. The SSRs analysis of the short-term cryopreserved material detected somaclonal variation. Both SSR and RAPD markers showed high sensitivity to detect genetic differences between the donor trees and the generated embryogenic material. Nevertheless, the degree of instability detection depended on the marker. The SSR analysis indicated a relationship between genotype, the studied loci and the located polymorphisms. Also, both markers revealed low intraclonal genetic variation. The RAPD detected genetic variation within the donor trees. The presence of pre-existent genetic variation within mature trees, in addition to its occurrence during the early stages of the embryogenic induction, and the presence of tissues of fertilisation origin within the maternal explants are all discussed. Nonetheless, the determination of the genetic identity of donor material is required, in addition to early detection methods of somaclonal variation.

Development of somatic embryogenesis for cloning and conservation of mature holm oak trees (Quercus ilex L.)

La encina (Quercus ilex L.) es una de las especies forestales mediterráneas más importantes. Constituye gran parte del estrato arbóreo de dehesas o montados, produce bellota como alimento del ganado y establece simbiosis con hongos micorrizógenos de gran valor económico. La encina está considerada como una especie recalcitrante en términos de conservación de semillas y capacidad morfogénica, lo que dificulta los programas de conservación de recursos genéticos y la mejora de la especie. La propagación vegetativa es una potente herramienta de los programas de mejora, por lo que es preciso desarrollar protocolos de regeneración somática en encina. La embriogénesis somática está considerada como la modalidad más adecuada de regeneración basada en técnicas de cultivo de tejidos vegetales utilizada en biotecnología forestal. Este trabajo se centra en el estudio de determinados aspectos de la embriogénesis somática para la regeneración clonal de encinas adultas. La memoria de esta tesis se ha dividido en capítulos que se corresponden con diferentes aspectos del sistema embriogénico. La embriogénesis somática se indujo en tegumentos maternos de óvulos en desarrollo procedentes de bellotas inmaduras de encinas adultas. A pesar de las bajas frecuencias de inducción, las líneas embriogénicas generadas se amplificaron mediante embriogénesis secundaria observándose cierta pérdida de la capacidad de diferenciación con el tiempo. Tanto el genotipo como la formulación del medio de cultivo influyeron en la respuesta embriogénica, concluyendo que la formulación de macronutrientes de Schenk y Hildebrant del medio sin reguladores de crecimiento fue la combinación más efectiva en la inducción. Los resultados sugirieron la existencia de una ventana en el desarrollo del óvulo más sensible a la inducción. El genotipo in[luyó en la capacidad proliferativa de los cultivos y en la conversión de los embriones somáticos, que se incrementó suplementando el medio con ácido indol-3-butírico y 6-benciladenina. El cultivo en medio líquido de líneas embriogénicas en condiciones de inmersión transitoria incrementó el crecimiento, dependiendo del genotipo, con respecto al cultivo en medio semisólido. Sin embargo, no mejoró la capacidad de diferenciar embriones cotiledonares aislados. Se estableció un protocolo de inicio y mantenimiento de cultivos en suspensión para varias líneas embriogénicas mediante inoculación en alta densidad de agregados embrionarios procedentes del medio semisólido. Para evitar la pérdida de vigor y la capacidad morfogénica debida al cultivo prolongado se desarrolló un protocolo de crioconservación de líneas embriogénicas mediante vitrificación. Al determinar la influencia de los agentes crioprotectores antes y después de su inmersión en nitrógeno líquido se concluyó que las respuestas de capacidad de crecimiento y de diferenciación del material embriogénico son independientes, además de estar bajo influencia del genotipo y el tipo de material crioconservado. La combinación de sacarosa y PVS2 previa a la inmersión en nitrógeno líquido proporcionó la mayor tasa de recuperación. Cuando las líneas fueron crioconservadas 30 días la capacidad de diferenciación se perdió en todas ellas. El análisis de SSR detectó variación somaclonal en el material crioconservado a corto plazo. SSR y RAPD mostraron importantes diferencias genéticas entre los árboles donantes y el material embriogénico que dependieron del genotipo. El grado de detección dependió del marcador empleado. Ambos marcadores revelaron baja inestabilidad intraclonal. Los RAPD revelaron variación genética intra-individuo en las encinas donantes. Se discuten la variación genética pre-existente en encina, su aparición durante las primeras fases de la inducción de embriogénesis, y la presencia de tejidos provenientes de la fertilización en el explanto materno. Esto hace preciso definir la identidad genética del material donante y acometer ensayos de detección precoz de variación somaclonal. ABSTRACT Holm oak (Quercus ilex L.) is one of the most important Mediterranean forest species. It conforms the tree layer of dehesas or montados, it produces acorns to feed the livestock and it establishes symbiosis with profitable mycorrhizal fungi. Holm oak is considered as recalcitrant species in terms of seed conservation and morphogenic capacities, which complicates the development of genetic conservation and improvement programs. Vegetative propagation is one of the mightiest tools for breeding programs therefore; developing protocols for clonal regeneration of holm oak is essential. Somatic embryogenesis is considered the best tissue culture-based way of plant regeneration in forest biotechnology. The present study is focused on the study of certain aspects of somatic embryogenesis for clonal regeneration of mature holm oak. This thesis manuscript is divided into several chapters that match with different aspects of the embryogenic system. Somatic embryogenesis induction was achieved on maternal teguments of developing ovules from immature acorns of adult holm oak trees. Despite the low induction frequencies, the generated embryogenic lines were amplified by secondary embryogenesis. A decline in the differentiation capacity over time was also observed. It was concluded that both genotype and culture media formulation influenced the embryogenic response, being the Schenk and Hildebrandt´s macronutrients formulation from culture medium and the lack of plant growth regulators the most effective combination for the induction of the embryogenic response. It has been suggested the existence of a developmental window in which ovules are prone to induction. Genotype influenced the proliferation capacity and the plant conversion of somatic embryos, which was also favoured by the presence of indol-3-butyric acid and 6-bencyladenine. The use of temporary immersion systems as proliferation in liquid culture of the embryogenic lines increased the growth depending on genotype, when compared to semisolid cultures. However, it did not improve the differentiation of single cotyledonary embryos. A protocol for the initiation and maintenance of embryogenic suspension cultures was established for several embryogenic lines with highly dense inoculi of embryogenic clusters from proliferating semisolid cultures. In order to avoid the loss of vigour and morphogenic ability of embryogenic lines due to prolonged cultures, a cryopreservation protocol for embryogenic lines of holm oak has been developed. During the determination of the influence of cryoprotective agents on the growth and differentiation capacities before and after liquid nitrogen immersion, it was concluded that both responses were independent from each other and also under the influence of genotype and the type of cryopreserved material. The combination of sucrose and PVS2 prior liquid nitrogen immersion provided higher recovery rates. When the same embryogenic lines were cryopreserved for 30 days, none was able to differentiate. The SSRs analysis of the short-term cryopreserved material detected somaclonal variation. Both SSR and RAPD markers showed high sensitivity to detect genetic differences between the donor trees and the generated embryogenic material. Nevertheless, the degree of instability detection depended on the marker. The SSR analysis indicated a relationship between genotype, the studied loci and the located polymorphisms. Also, both markers revealed low intraclonal genetic variation. The RAPD detected genetic variation within the donor trees. The presence of pre-existent genetic variation within mature trees, in addition to its occurrence during the early stages of the embryogenic induction, and the presence of tissues of fertilisation origin within the maternal explants are all discussed. Nonetheless, the determination of the genetic identity of donor material is required, in addition to early detection methods of somaclonal variation.

Climate variability effects on field crops in Iberian Peninsula. Predictability

La presente Tesis constituye un avance en el conocimiento de los efectos de la variabilidad climática en los cultivos en la Península Ibérica (PI). Es bien conocido que la temperatura del océano, particularmente de la región tropical, es una de las variables más convenientes para ser utilizado como predictor climático. Los océanos son considerados como la principal fuente de almacenamiento de calor del planeta debido a la alta capacidad calorífica del agua. Cuando se libera esta energía, altera los regímenes globales de circulación atmosférica por mecanismos de teleconexión. Estos cambios en la circulación general de la atmósfera afectan a la temperatura, precipitación, humedad, viento, etc., a escala regional, los cuales afectan al crecimiento, desarrollo y rendimiento de los cultivos. Para el caso de Europa, esto implica que la variabilidad atmosférica en una región específica se asocia con la variabilidad de otras regiones adyacentes y/o remotas, como consecuencia Europa está siendo afectada por los patrones de circulaciones globales, que a su vez, se ven afectados por patrones oceánicos. El objetivo general de esta tesis es analizar la variabilidad del rendimiento de los cultivos y su relación con la variabilidad climática y teleconexiones, así como evaluar su predictibilidad. Además, esta Tesis tiene como objetivo establecer una metodología para estudiar la predictibilidad de las anomalías del rendimiento de los cultivos. El análisis se centra en trigo y maíz como referencia para otros cultivos de la PI, cultivos de invierno en secano y cultivos de verano en regadío respectivamente. Experimentos de simulación de cultivos utilizando una metodología en cadena de modelos (clima + cultivos) son diseñados para evaluar los impactos de los patrones de variabilidad climática en el rendimiento y su predictibilidad. La presente Tesis se estructura en dos partes: La primera se centra en el análisis de la variabilidad del clima y la segunda es una aplicación de predicción cuantitativa de cosechas. La primera parte está dividida en 3 capítulos y la segundo en un capitulo cubriendo los objetivos específicos del presente trabajo de investigación. Parte I. Análisis de variabilidad climática El primer capítulo muestra un análisis de la variabilidad del rendimiento potencial en una localidad como indicador bioclimático de las teleconexiones de El Niño con Europa, mostrando su importancia en la mejora de predictibilidad tanto en clima como en agricultura. Además, se presenta la metodología elegida para relacionar el rendimiento con las variables atmosféricas y oceánicas. El rendimiento de los cultivos es parcialmente determinado por la variabilidad climática atmosférica, que a su vez depende de los cambios en la temperatura de la superficie del mar (TSM). El Niño es el principal modo de variabilidad interanual de la TSM, y sus efectos se extienden en todo el mundo. Sin embargo, la predictibilidad de estos impactos es controversial, especialmente aquellos asociados con la variabilidad climática Europea, que se ha encontrado que es no estacionaria y no lineal. Este estudio mostró cómo el rendimiento potencial de los cultivos obtenidos a partir de datos de reanálisis y modelos de cultivos sirve como un índice alternativo y más eficaz de las teleconexiones de El Niño, ya que integra las no linealidades entre las variables climáticas en una única serie temporal. Las relaciones entre El Niño y las anomalías de rendimiento de los cultivos son más significativas que las contribuciones individuales de cada una de las variables atmosféricas utilizadas como entrada en el modelo de cultivo. Además, la no estacionariedad entre El Niño y la variabilidad climática europea se detectan con mayor claridad cuando se analiza la variabilidad de los rendimiento de los cultivos. La comprensión de esta relación permite una cierta predictibilidad hasta un año antes de la cosecha del cultivo. Esta predictibilidad no es constante, sino que depende tanto la modulación de la alta y baja frecuencia. En el segundo capítulo se identifica los patrones oceánicos y atmosféricos de variabilidad climática que afectan a los cultivos de verano en la PI. Además, se presentan hipótesis acerca del mecanismo eco-fisiológico a través del cual el cultivo responde. Este estudio se centra en el análisis de la variabilidad del rendimiento de maíz en la PI para todo el siglo veinte, usando un modelo de cultivo calibrado en 5 localidades españolas y datos climáticos de reanálisis para obtener series temporales largas de rendimiento potencial. Este estudio evalúa el uso de datos de reanálisis para obtener series de rendimiento de cultivos que dependen solo del clima, y utilizar estos rendimientos para analizar la influencia de los patrones oceánicos y atmosféricos. Los resultados muestran una gran fiabilidad de los datos de reanálisis. La distribución espacial asociada a la primera componente principal de la variabilidad del rendimiento muestra un comportamiento similar en todos los lugares estudiados de la PI. Se observa una alta correlación lineal entre el índice de El Niño y el rendimiento, pero no es estacionaria en el tiempo. Sin embargo, la relación entre la temperatura del aire y el rendimiento se mantiene constante a lo largo del tiempo, siendo los meses de mayor influencia durante el período de llenado del grano. En cuanto a los patrones atmosféricos, el patrón Escandinavia presentó una influencia significativa en el rendimiento en PI. En el tercer capítulo se identifica los patrones oceánicos y atmosféricos de variabilidad climática que afectan a los cultivos de invierno en la PI. Además, se presentan hipótesis acerca del mecanismo eco-fisiológico a través del cual el cultivo responde. Este estudio se centra en el análisis de la variabilidad del rendimiento de trigo en secano del Noreste (NE) de la PI. La variabilidad climática es el principal motor de los cambios en el crecimiento, desarrollo y rendimiento de los cultivos, especialmente en los sistemas de producción en secano. En la PI, los rendimientos de trigo son fuertemente dependientes de la cantidad de precipitación estacional y la distribución temporal de las mismas durante el periodo de crecimiento del cultivo. La principal fuente de variabilidad interanual de la precipitación en la PI es la Oscilación del Atlántico Norte (NAO), que se ha relacionado, en parte, con los cambios en la temperatura de la superficie del mar en el Pacífico Tropical (El Niño) y el Atlántico Tropical (TNA). La existencia de cierta predictibilidad nos ha animado a analizar la posible predicción de los rendimientos de trigo en la PI utilizando anomalías de TSM como predictor. Para ello, se ha utilizado un modelo de cultivo (calibrado en dos localidades del NE de la PI) y datos climáticos de reanálisis para obtener series temporales largas de rendimiento de trigo alcanzable y relacionar su variabilidad con anomalías de la TSM. Los resultados muestran que El Niño y la TNA influyen en el desarrollo y rendimiento del trigo en el NE de la PI, y estos impactos depende del estado concurrente de la NAO. Aunque la relación cultivo-TSM no es igual durante todo el periodo analizado, se puede explicar por un mecanismo eco-fisiológico estacionario. Durante la segunda mitad del siglo veinte, el calentamiento (enfriamiento) en la superficie del Atlántico tropical se asocia a una fase negativa (positiva) de la NAO, que ejerce una influencia positiva (negativa) en la temperatura mínima y precipitación durante el invierno y, por lo tanto, aumenta (disminuye) el rendimiento de trigo en la PI. En relación con El Niño, la correlación más alta se observó en el período 1981 -2001. En estas décadas, los altos (bajos) rendimientos se asocian con una transición El Niño - La Niña (La Niña - El Niño) o con eventos de El Niño (La Niña) que están finalizando. Para estos eventos, el patrón atmosférica asociada se asemeja a la NAO, que también influye directamente en la temperatura máxima y precipitación experimentadas por el cultivo durante la floración y llenado de grano. Los co- efectos de los dos patrones de teleconexión oceánicos ayudan a aumentar (disminuir) la precipitación y a disminuir (aumentar) la temperatura máxima en PI, por lo tanto el rendimiento de trigo aumenta (disminuye). Parte II. Predicción de cultivos. En el último capítulo se analiza los beneficios potenciales del uso de predicciones climáticas estacionales (por ejemplo de precipitación) en las predicciones de rendimientos de trigo y maíz, y explora métodos para aplicar dichos pronósticos climáticos en modelos de cultivo. Las predicciones climáticas estacionales tienen un gran potencial en las predicciones de cultivos, contribuyendo de esta manera a una mayor eficiencia de la gestión agrícola, seguridad alimentaria y de subsistencia. Los pronósticos climáticos se expresan en diferentes formas, sin embargo todos ellos son probabilísticos. Para ello, se evalúan y aplican dos métodos para desagregar las predicciones climáticas estacionales en datos diarios: 1) un generador climático estocástico condicionado (predictWTD) y 2) un simple re-muestreador basado en las probabilidades del pronóstico (FResampler1). Los dos métodos se evaluaron en un caso de estudio en el que se analizaron los impactos de tres escenarios de predicciones de precipitación estacional (predicción seco, medio y lluvioso) en el rendimiento de trigo en secano, sobre las necesidades de riego y rendimiento de maíz en la PI. Además, se estimó el margen bruto y los riesgos de la producción asociada con las predicciones de precipitación estacional extremas (seca y lluviosa). Los métodos predWTD y FResampler1 usados para desagregar los pronósticos de precipitación estacional en datos diarios, que serán usados como inputs en los modelos de cultivos, proporcionan una predicción comparable. Por lo tanto, ambos métodos parecen opciones factibles/viables para la vinculación de los pronósticos estacionales con modelos de simulación de cultivos para establecer predicciones de rendimiento o las necesidades de riego en el caso de maíz. El análisis del impacto en el margen bruto de los precios del grano de los dos cultivos (trigo y maíz) y el coste de riego (maíz) sugieren que la combinación de los precios de mercado previstos y la predicción climática estacional pueden ser una buena herramienta en la toma de decisiones de los agricultores, especialmente en predicciones secas y/o localidades con baja precipitación anual. Estos métodos permiten cuantificar los beneficios y riesgos de los agricultores ante una predicción climática estacional en la PI. Por lo tanto, seríamos capaces de establecer sistemas de alerta temprana y diseñar estrategias de adaptación del manejo del cultivo para aprovechar las condiciones favorables o reducir los efectos de condiciones adversas. La utilidad potencial de esta Tesis es la aplicación de las relaciones encontradas para predicción de cosechas de la próxima campaña agrícola. Una correcta predicción de los rendimientos podría ayudar a los agricultores a planear con antelación sus prácticas agronómicas y todos los demás aspectos relacionados con el manejo de los cultivos. Esta metodología se puede utilizar también para la predicción de las tendencias futuras de la variabilidad del rendimiento en la PI. Tanto los sectores públicos (mejora de la planificación agrícola) como privados (agricultores, compañías de seguros agrarios) pueden beneficiarse de esta mejora en la predicción de cosechas. ABSTRACT The present thesis constitutes a step forward in advancing of knowledge of the effects of climate variability on crops in the Iberian Peninsula (IP). It is well known that ocean temperature, particularly the tropical ocean, is one of the most convenient variables to be used as climate predictor. Oceans are considered as the principal heat storage of the planet due to the high heat capacity of water. When this energy is released, it alters the global atmospheric circulation regimes by teleconnection1 mechanisms. These changes in the general circulation of the atmosphere affect the regional temperature, precipitation, moisture, wind, etc., and those influence crop growth, development and yield. For the case of Europe, this implies that the atmospheric variability in a specific region is associated with the variability of others adjacent and/or remote regions as a consequence of Europe being affected by global circulations patterns which, in turn, are affected by oceanic patterns. The general objective of this Thesis is to analyze the variability of crop yields at climate time scales and its relation to the climate variability and teleconnections, as well as to evaluate their predictability. Moreover, this Thesis aims to establish a methodology to study the predictability of crop yield anomalies. The analysis focuses on wheat and maize as a reference crops for other field crops in the IP, for winter rainfed crops and summer irrigated crops respectively. Crop simulation experiments using a model chain methodology (climate + crop) are designed to evaluate the impacts of climate variability patterns on yield and its predictability. The present Thesis is structured in two parts. The first part is focused on the climate variability analyses, and the second part is an application of the quantitative crop forecasting for years that fulfill specific conditions identified in the first part. This Thesis is divided into 4 chapters, covering the specific objectives of the present research work. Part I. Climate variability analyses The first chapter shows an analysis of potential yield variability in one location, as a bioclimatic indicator of the El Niño teleconnections with Europe, putting forward its importance for improving predictability in both climate and agriculture. It also presents the chosen methodology to relate yield with atmospheric and oceanic variables. Crop yield is partially determined by atmospheric climate variability, which in turn depends on changes in the sea surface temperature (SST). El Niño is the leading mode of SST interannual variability, and its impacts extend worldwide. Nevertheless, the predictability of these impacts is controversial, especially those associated with European climate variability, which have been found to be non-stationary and non-linear. The study showed how potential2 crop yield obtained from reanalysis data and crop models serves as an alternative and more effective index of El Niño teleconnections because it integrates the nonlinearities between the climate variables in a unique time series. The relationships between El Niño and crop yield anomalies are more significant than the individual contributions of each of the atmospheric variables used as input in the crop model. Additionally, the non-stationarities between El Niño and European climate variability are more clearly detected when analyzing crop-yield variability. The understanding of this relationship allows for some predictability up to one year before the crop is harvested. This predictability is not constant, but depends on both high and low frequency modulation. The second chapter identifies the oceanic and atmospheric patterns of climate variability affecting summer cropping systems in the IP. Moreover, hypotheses about the eco-physiological mechanism behind crop response are presented. It is focused on an analysis of maize yield variability in IP for the whole twenty century, using a calibrated crop model at five contrasting Spanish locations and reanalyses climate datasets to obtain long time series of potential yield. The study tests the use of reanalysis data for obtaining only climate dependent time series of simulated crop yield for the whole region, and to use these yield to analyze the influences of oceanic and atmospheric patterns. The results show a good reliability of reanalysis data. The spatial distribution of the leading principal component of yield variability shows a similar behaviour over all the studied locations in the IP. The strong linear correlation between El Niño index and yield is remarkable, being this relation non-stationary on time, although the air temperature-yield relationship remains on time, being the highest influences during grain filling period. Regarding atmospheric patterns, the summer Scandinavian pattern has significant influence on yield in IP. The third chapter identifies the oceanic and atmospheric patterns of climate variability affecting winter cropping systems in the IP. Also, hypotheses about the eco-physiological mechanism behind crop response are presented. It is focused on an analysis of rainfed wheat yield variability in IP. Climate variability is the main driver of changes in crop growth, development and yield, especially for rainfed production systems. In IP, wheat yields are strongly dependent on seasonal rainfall amount and temporal distribution of rainfall during the growing season. The major source of precipitation interannual variability in IP is the North Atlantic Oscillation (NAO) which has been related in part with changes in the Tropical Pacific (El Niño) and Atlantic (TNA) sea surface temperature (SST). The existence of some predictability has encouraged us to analyze the possible predictability of the wheat yield in the IP using SSTs anomalies as predictor. For this purpose, a crop model with a site specific calibration for the Northeast of IP and reanalysis climate datasets have been used to obtain long time series of attainable wheat yield and relate their variability with SST anomalies. The results show that El Niño and TNA influence rainfed wheat development and yield in IP and these impacts depend on the concurrent state of the NAO. Although crop-SST relationships do not equally hold on during the whole analyzed period, they can be explained by an understood and stationary ecophysiological mechanism. During the second half of the twenty century, the positive (negative) TNA index is associated to a negative (positive) phase of NAO, which exerts a positive (negative) influence on minimum temperatures (Tmin) and precipitation (Prec) during winter and, thus, yield increases (decreases) in IP. In relation to El Niño, the highest correlation takes place in the period 1981-2001. For these decades, high (low) yields are associated with an El Niño to La Niña (La Niña to El Niño) transitions or to El Niño events finishing. For these events, the regional associated atmospheric pattern resembles the NAO, which also influences directly on the maximum temperatures (Tmax) and precipitation experienced by the crop during flowering and grain filling. The co-effects of the two teleconnection patterns help to increase (decrease) the rainfall and decrease (increase) Tmax in IP, thus on increase (decrease) wheat yield. Part II. Crop forecasting The last chapter analyses the potential benefits for wheat and maize yields prediction from using seasonal climate forecasts (precipitation), and explores methods to apply such a climate forecast to crop models. Seasonal climate prediction has significant potential to contribute to the efficiency of agricultural management, and to food and livelihood security. Climate forecasts come in different forms, but probabilistic. For this purpose, two methods were evaluated and applied for disaggregating seasonal climate forecast into daily weather realizations: 1) a conditioned stochastic weather generator (predictWTD) and 2) a simple forecast probability resampler (FResampler1). The two methods were evaluated in a case study where the impacts of three scenarios of seasonal rainfall forecasts on rainfed wheat yield, on irrigation requirements and yields of maize in IP were analyzed. In addition, we estimated the economic margins and production risks associated with extreme scenarios of seasonal rainfall forecasts (dry and wet). The predWTD and FResampler1 methods used for disaggregating seasonal rainfall forecast into daily data needed by the crop simulation models provided comparable predictability. Therefore both methods seem feasible options for linking seasonal forecasts with crop simulation models for establishing yield forecasts or irrigation water requirements. The analysis of the impact on gross margin of grain prices for both crops and maize irrigation costs suggests the combination of market prices expected and the seasonal climate forecast can be a good tool in farmer’s decision-making, especially on dry forecast and/or in locations with low annual precipitation. These methodologies would allow quantifying the benefits and risks of a seasonal weather forecast to farmers in IP. Therefore, we would be able to establish early warning systems and to design crop management adaptation strategies that take advantage of favorable conditions or reduce the effect of adverse conditions. The potential usefulness of this Thesis is to apply the relationships found to crop forecasting on the next cropping season, suggesting opportunity time windows for the prediction. The methodology can be used as well for the prediction of future trends of IP yield variability. Both public (improvement of agricultural planning) and private (decision support to farmers, insurance companies) sectors may benefit from such an improvement of crop forecasting.