Parallelization of whole genome alignment

With the advent of High performance computing, it is now possible to achieve orders of magnitude performance and computation e ciency gains over conventional computer architectures. This thesis explores the potential of using high performance computing to accelerate whole genome alignment. A parallel technique is applied to an algorithm for whole genome alignment, this technique is explained and some experiments were carried out to test it. This technique is based in a fair usage of the available resource to execute genome alignment and how this can be used in HPC clusters. This work is a rst approximation to whole genome alignment and it shows the advantages of parallelism and some of the drawbacks that our technique has. This work describes the resource limitations of current WGA applications when dealing with large quantities of sequences. It proposes a parallel heuristic to distribute the load and to assure that alignment quality is mantained.

Seasonal changes of whole root system conductance by a drought-tolerant grape root system

The role of root systems in drought tolerance is a subject of very limited information compared with above-ground responses. Adjustments to the ability of roots to supply water relative to shoot transpiration demand is proposed as a major means for woody perennial plants to tolerate drought, and is often expressed as changes in the ratios of leaf to root area (AL:AR). Seasonal root proliferation in a directed manner could increase the water supply function of roots independent of total root area (AR) and represents a mechanism whereby water supply to demand could be increased. To address this issue, seasonal root proliferation, stomatal conductance (gs) and whole root system hydraulic conductance (kr) were investigated for a drought-tolerant grape root system (Vitis berlandieri×V. rupestris cv. 1103P) and a non-drought-tolerant root system (Vitis riparia×V. rupestris cv. 101-14Mgt), upon which had been grafted the same drought-sensitive clone of Vitis vinifera cv. Merlot. Leaf water potentials (ψL) for Merlot grafted onto the 1103P root system (–0.91±0.02 MPa) were +0.15 MPa higher than Merlot on 101-14Mgt (–1.06±0.03 MPa) during spring, but dropped by approximately –0.4 MPa from spring to autumn, and were significantly lower by –0.15 MPa (–1.43±0.02 MPa) than for Merlot on 101-14Mgt (at –1.28±0.02 MPa). Surprisingly, gs of Merlot on the drought-tolerant root system (1103P) was less down-regulated and canopies maintained evaporative fluxes ranging from 35–20 mmol vine−1 s−1 during the diurnal peak from spring to autumn, respectively, three times greater than those measured for Merlot on the drought-sensitive rootstock 101-14Mgt. The drought-tolerant root system grew more roots at depth during the warm summer dry period, and the whole root system conductance (kr) increased from 0.004 to 0.009 kg MPa−1 s−1 during that same time period. The changes in kr could not be explained by xylem anatomy or conductivity changes of individual root segments. Thus, the manner in which drought tolerance was conveyed to the drought-sensitive clone appeared to arise from deep root proliferation during the hottest and driest part of the season, rather than through changes in xylem structure, xylem density or stomatal regulation. This information can be useful to growers on a site-specific basis in selecting rootstocks for grape clonal material (scions) grafted to them.

Evaluation of four whole-plant inoculation methods to analyze the pathogenicity of Erwinia amylovora under quarantine conditions

Four methods were tested to assess the fire-blight disease response on grafted pear plants. The leaves of the plants were inoculated with Erwinia amylovora suspensions by pricking with clamps, cutting with scissors, local infiltration, and painting a bacterial suspension onto the leaves with a paintbrush. The effects of the inoculation methods were studied in dose-time-response experiments carried out in climate chambers under quarantine conditions. A modified Gompertz model was used to analyze the disease-time relatiobbnships and provided information on the rate of infection progression (rg) and time delay to the start of symptoms (t0). The disease-pathogen-dose relationships were analyzed according to a hyperbolic saturation model in which the median effective dose (ED50) of the pathogen and maximum disease level (ymax) were determined. Localized infiltration into the leaf mesophile resulted in the early (short t0) but slow (low rg) development of infection whereas in leaves pricked with clamps disease symptoms developed late (long t0) but rapidly (high rg). Paintbrush inoculation of the plants resulted in an incubation period of medium length, a moderate rate of infection progression, and low ymax values. In leaves inoculated with scissors, fire-blight symptoms developed early (short t0) and rapidly (high rg), and with the lowest ED50 and the highest ymax

Whole genome analysis of p38 SAPK-mediated gene expression upon stress

Background: Cells have the ability to respond and adapt to environmental changes through activation of stress-activated protein kinases (SAPKs). Although p38 SAPK signalling is known to participate in the regulation of gene expression little is known on the molecular mechanisms used by this SAPK to regulate stress-responsive genes and the overall set of genes regulated by p38 in response to different stimuli.Results: Here, we report a whole genome expression analyses on mouse embryonic fibroblasts (MEFs) treated with three different p38 SAPK activating-stimuli, namely osmostress, the cytokine TNFα and the protein synthesis inhibitor anisomycin. We have found that the activation kinetics of p38α SAPK in response to these insults is different and also leads to a complex gene pattern response specific for a given stress with a restricted set of overlapping genes. In addition, we have analysed the contribution of p38α the major p38 family member present in MEFs, to the overall stress-induced transcriptional response by using both a chemical inhibitor (SB203580) and p38α deficient (p38α-/-) MEFs. We show here that p38 SAPK dependency ranged between 60% and 88% depending on the treatments and that there is a very good overlap between the inhibitor treatment and the ko cells. Furthermore, we have found that the dependency of SAPK varies depending on the time the cells are subjected to osmostress. Conclusions: Our genome-wide transcriptional analyses shows a selective response to specific stimuli and a restricted common response of up to 20% of the stress up-regulated early genes that involves an important set of transcription factors, which might be critical for either cell adaptation or preparation for continuous extra-cellular changes. Interestingly, up to 85% of the up-regulated genes are under the transcriptional control of p38 SAPK. Thus, activation of p38 SAPK is critical to elicit the early gene expression program required for cell adaptation to stress.

Bilinear models for spatio-temporal point distribution analysis: application to extrapolation of left ventricular, biventricular and whole heart cardiac dynamics

In this work we describe the usage of bilinear statistical models as a means of factoring the shape variability into two components attributed to inter-subject variation and to the intrinsic dynamics of the human heart. We show that it is feasible to reconstruct the shape of the heart at discrete points in the cardiac cycle. Provided we are given a small number of shape instances representing the same heart atdifferent points in the same cycle, we can use the bilinearmodel to establish this. Using a temporal and a spatial alignment step in the preprocessing of the shapes, around half of the reconstruction errors were on the order of the axial image resolution of 2 mm, and over 90% was within 3.5 mm. From this, weconclude that the dynamics were indeed separated from theinter-subject variability in our dataset.

Feature Extraction from Whole-Sky Ground-Based Images for Cloud-Type Recognition

Several features that can be extracted from digital images of the sky and that can be useful for cloud-type classification of such images are presented. Some features are statistical measurements of image texture, some are based on the Fourier transform of the image and, finally, others are computed from the image where cloudy pixels are distinguished from clear-sky pixels. The use of the most suitable features in an automatic classification algorithm is also shown and discussed. Both the features and the classifier are developed over images taken by two different camera devices, namely, a total sky imager (TSI) and a whole sky imager (WSC), which are placed in two different areas of the world (Toowoomba, Australia; and Girona, Spain, respectively). The performance of the classifier is assessed by comparing its image classification with an a priori classification carried out by visual inspection of more than 200 images from each camera. The index of agreement is 76% when five different sky conditions are considered: clear, low cumuliform clouds, stratiform clouds (overcast), cirriform clouds, and mottled clouds (altocumulus, cirrocumulus). Discussion on the future directions of this research is also presented, regarding both the use of other features and the use of other classification techniques

Whole-rat protein content estimation: applicability of the Nx6.25 factor

The amino acid composition of the protein from three strains of rat (Wistar, Zucker lean and Zucker obese), subjected to reference and high-fat diets has been used to determine the mean empirical formula, molecular weight and N content of whole-rat protein. The combined whole protein of the rat was uniform for the six experimental groups, containing an estimate of 17.3% N and a mean aminoacyl residue molecular weight of 103.7. This suggests that the appropriate protein factor for the calculation of rat protein from its N content should be 5.77 instead of the classical 6.25. In addition, an estimate of the size of the non-protein N mass in the whole rat gave a figure in the range of 5.5 % of all N. The combination of the two calculations gives a protein factor of 5.5 for the conversion of total N into rat protein.

Propuesta de implementación del cultivo del manzano (Malus domestica) en la Vall d'Alinyà : análisis de viabilidad ambiental-económica y certificación de créditos voluntarios de carbono

El presente proyecto tiene por objetivo realizar una propuesta de cultivo de manzano (Malus domestica) en la Vall d’Alinyà para el desarrollo de un producto agroalimentario que contribuya a la revitalización de esta área rural y a la mitigación del calentamiento global. Para ello, se pretende proporcionar una metodología para calcular, con mayor precisión, el comportamiento de los cultivos como sumidero de carbono a partir del balance neto de gases con efecto invernadero (GEI) del sistema. Con tal de obtener las emisiones de GEI generadas por los cultivos se desarrolla un análisis del ciclo de vida (ACV), y para obtener la fijación de los cultivos, se realiza un análisis de su capacidad de absorción de dióxido de carbono en la biomasa de la plantación. Además, mediante el uso de sistemas de información geográfica (SIG), se determinaron los terrenos potenciales para el cultivo del manzano en la zona de estudio. Posteriormente, se verifica la certificación de créditos en el mercado de carbono voluntario y se analiza la viabilidad económica del proyecto, obteniendo así un producto (la manzana) con valor ambiental añadido. Como conclusión se obtuvo una viabilidad positiva de la verificación de créditos, puesto que el balance neto de carbono fue positivo, absorbiéndose 234,54 t CO2 en todo el territorio potencialmente cultivable (16,92 ha) y durante 15 años de actividad agrícola. Al mismo tiempo, la propuesta resultó viable económicamente, generándose unos beneficios ligados a la venta de producción frutícola y de los créditos de carbono de 79.484 € durante los 15 años de actividad productiva.

Assessing farmers' vulnerability to climate change : a case study in Karnataka, India

In the context of observed climate change impacts and their effect on agriculture and crop production, this study intends to assess the vulnerability of rural livelihoods through a study case in Karnataka, India. The social approach of climate change vulnerability in this study case includes defining and exploring factors that determine farmers’ vulnerability in four villages. Key informant interviews, farmer workshops and structured household interviews were used for data collection. To analyse the data, we adapted and applied three vulnerability indices: Livelihood Vulnerability Index (LVI), LVI-IPCC and the Livelihood Effect Index (LEI), and used descriptive statistical methods. The data was analysed at two scales: whole sample-level and household level. The results from applying the indices for the whole-sample level show that this community's vulnerability to climate change is moderate, whereas the household-level results show that most of the households' vulnerability is high-very high, while 15 key drivers of vulnerability were identified. Results and limitations of the study are discussed under the rural livelihoods framework, in which the indices are based, allowing a better understanding of the social behaviouraltrends, as well as an holistic and integrated view of the climate change, agriculture, and livelihoods processes shaping vulnerability. We conclude that these indices, although a straightforward method to assess vulnerability, have limitations that could account for inaccuracies and inability to be standardised for benchmarking, therefore we stress the need for further research.

The genome of the recently domesticated crop plant sugar beet (Beta vulgaris)

Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world's annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714-758 megabases and shares an ancient genome triplication with other eudicot plants. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant's potential in energy biotechnology.

A germination study of herbicide-resistant field poppies in Spain

Field poppy, Papaver rhoeas L., is a very common weed in winter cereals in North-Eastern Spain. Its control is becoming difficult due to expanding herbicide resistance. To control field poppies there are alternative strategies such as non-chemical control that take into account the weed emergence period. However, there is a lack of knowledge of P. rhoeas emergence patterns in semi-arid conditions. Thus, here we conducted pot experiments on the emergence of P. rhoeas. We aimed to describe the emergence period and to quantify the emergence of a susceptible and of a herbicide-resistant P. rhoeas population at two locations in Catalonia, Spain, from 1998 to 2001 and until 2004 at one of them. Therefore, pots containing seeds of both populations were established at the two locations and emergence was recorded monthly. We studied the origin of the population, the sowing location, the effect of cultivation and the sowing year. First, we found that the main emergence peaks in our experiments occurred in autumn, accounting for between 65.7 and 98.5% of the annual emergence from October to December, and only little emergence was recorded in spring. This emergence pattern is different from those found in the literature corresponding to Northern European countries, where in some cases main flushes occur only in autumn, in spring and winter or only in spring. The emergence was mainly affected by cultivation, but the effect of light stimulus was observed several months later. As a consequence, cultivation should be done in early autumn, promoting emergence during the whole autumn and winter so that emerged seedlings can be controlled before sowing a spring crop. Second, most experiments showed that the emergence was significantly higher in the first autumn than in the following seasons, e.g. 4.1% emergence in the first year and only 2.1, 2.3, 0.5 and 0.6% new emergence at one of the locations for the second, third, fourth and fifth years. Thus, after having a severe P. rhoeas infestation causing a big seed rain, emergence should be stimulated by autumn cultivation in the following season and seedlings controlled by trying to deplete the soil seed bank as much as possible. Despite the fact that emergence will be staggered throughout several years and that there was a significant relationship between rainfall and emergence, so that dry years will cause a smaller emergence rate of the weed, these findings define a cultural management strategy to reduce P. rhoeas infestations and to contribute to integrated weed management strategies combining it with other tools.

Effects of crop and weed densities on the interactions between barley and Lolium rigidum in several Mediterranean locations

The effects of both barley and Lolium rigidum densities on weed growth and spike production and on crop yield were examined in five field experiments carried out in the Mediterranean drylands of Spain and Western Australia. The aim was to check the consistency of the competitiveness of the crop in different environmental and management conditions. L. rigidum reduced barley yields in most of the experiments (between 0 and 85%), the number of ears per m2 being the most affected. It was found that increasing the barley seeding rate did not reduce the crop losses but did limit weed biomass (between 5 and 61%) and spike production (between 24 and 85%). The variability observed in crop yield losses between sites and seasons was related to rainfall at the beginning of the season. The most sensitive component of yield to weed competition was the number of ears per plant.

A priori parameterisation of the CERES soil-crop models and tests against several European data sets

Mechanistic soil-crop models have become indispensable tools to investigate the effect of management practices on the productivity or environmental impacts of arable crops. Ideally these models may claim to be universally applicable because they simulate the major processes governing the fate of inputs such as fertiliser nitrogen or pesticides. However, because they deal with complex systems and uncertain phenomena, site-specific calibration is usually a prerequisite to ensure their predictions are realistic. This statement implies that some experimental knowledge on the system to be simulated should be available prior to any modelling attempt, and raises a tremendous limitation to practical applications of models. Because the demand for more general simulation results is high, modellers have nevertheless taken the bold step of extrapolating a model tested within a limited sample of real conditions to a much larger domain. While methodological questions are often disregarded in this extrapolation process, they are specifically addressed in this paper, and in particular the issue of models a priori parameterisation. We thus implemented and tested a standard procedure to parameterize the soil components of a modified version of the CERES models. The procedure converts routinely-available soil properties into functional characteristics by means of pedo-transfer functions. The resulting predictions of soil water and nitrogen dynamics, as well as crop biomass, nitrogen content and leaf area index were compared to observations from trials conducted in five locations across Europe (southern Italy, northern Spain, northern France and northern Germany). In three cases, the model’s performance was judged acceptable when compared to experimental errors on the measurements, based on a test of the model’s root mean squared error (RMSE). Significant deviations between observations and model outputs were however noted in all sites, and could be ascribed to various model routines. In decreasing importance, these were: water balance, the turnover of soil organic matter, and crop N uptake. A better match to field observations could therefore be achieved by visually adjusting related parameters, such as field-capacity water content or the size of soil microbial biomass. As a result, model predictions fell within the measurement errors in all sites for most variables, and the model’s RMSE was within the range of published values for similar tests. We conclude that the proposed a priori method yields acceptable simulations with only a 50% probability, a figure which may be greatly increased through a posteriori calibration. Modellers should thus exercise caution when extrapolating their models to a large sample of pedo-climatic conditions for which they have only limited information.