Network-based analysis of genome-wide scans of natural selection

The genomic era brought by recent advances in the next-generation sequencing technology makes the genome-wide scans of natural selection a reality. Currently, almost all the statistical tests and analytical methods for identifying genes under selection was performed on the individual gene basis. Although these methods have the power of identifying gene subject to strong selection, they have limited power in discovering genes targeted by moderate or weak selection forces, which are crucial for understanding the molecular mechanisms of complex phenotypes and diseases. Recent availability and rapid completeness of many gene network and protein-protein interaction databases accompanying the genomic era open the avenues of exploring the possibility of enhancing the power of discovering genes under natural selection. The aim of the thesis is to explore and develop normal mixture model based methods for leveraging gene network information to enhance the power of natural selection target gene discovery. The results show that the developed statistical method, which combines the posterior log odds of the standard normal mixture model and the Guilt-By-Association score of the gene network in a naïve Bayes framework, has the power to discover moderate/weak selection gene which bridges the genes under strong selection and it helps our understanding the biology under complex diseases and related natural selection phenotypes.^

Intrinsic oxidative stress in cancer cells: A biochemical basis for therapeutic selectivity using ROS -generating agents

A major goal of chemotherapy is to selectively kill cancer cells while minimizing toxicity to normal cells. Identifying biological differences between cancer and normal cells is essential in designing new strategies to improve therapeutic selectivity. Superoxide dismutases (SOD) are crucial antioxidant enzymes required for the elimination of superoxide (O2·− ), a free radical produced during normal cellular metabolism. Previous studies in our laboratory demonstrated that 2-methoxyestradiol (2-ME), an estradiol derivative, inhibits the function of SOD and selectively kills human leukemia cells without exhibiting significant cytotoxicity in normal lymphocytes. The present work was initiated to examine the biochemical basis for the selective anticancer activity of 2-ME. Investigations using two-parameter flow cytometric analyses and ROS scavengers established that O2·− is a primary and essential mediator of 2-ME-induced apoptosis in cancer cells. In addition, experiments using SOD overexpression vectors and SOD knockout cells found that SOD is a critical target of 2-ME. Importantly, the administration of 2-ME resulted in the selective accumulation of O 2·− and apoptosis in leukemia and ovarian cancer cells. The preferential activity of 2-ME was found to be due to increased intrinsic oxidative stress in these cancer cells versus their normal counterparts. This intrinsic oxidative stress was associated with the upregulation of the antioxidant enzymes SOD and catalase as a mechanism to cope with the increase in ROS. Furthermore, oxygen consumption experiments revealed that normal lymphocytes decrease their respiration rate in response to 2-ME-induced oxidative stress, while human leukemia cells seem to lack this regulatory mechanism. This leads to an uncontrolled production of O2·−, severe accumulation of ROS, and ultimately ROS-mediated apoptosis in leukemia cells treated with 2-ME. The biochemical differences between cancer and normal cells identified here provide a basis for the development of drug combination strategies using 2-ME with other ROS-generating agents to enhance anticancer activity. The effectiveness of such a combination strategy in killing cancer cells was demonstrated by the use of 2-ME with agents/modalities such as ionizing radiation and doxorubicin. Collectively, the data presented here strongly suggests that 2-ME may have important clinical implications for the selective killing of cancer cells. ^

Data base on heterotrophic dinoflagellates grazing and growth rate as a function of size, prey size and type compiled from the literature

Microzooplankton (the 20 to 200 µm size class of zooplankton) is recognised as an important part of marine pelagic ecosystems. In terms of biomass and abundance heterotrophic dinoflagellates are one of the important groups of organism in microzooplankton. However, their rates - grazing and growth - , feeding behaviour and prey preferences are poorly known and understood. A set of data was assembled in order to derive a better understanding of heterotrophic dinoflagellates rates, in response to parameters such as prey concentration, prey type (size and species), temperature and their own size. With these objectives, literature was searched for laboratory experiments with information on one or more of these parameters effect studied. The criteria for selection and inclusion in the database included: (i) controlled laboratory experiment with a known dinoflagellate feeding on a known prey; (ii) presence of ancillary information about experimental conditions, used organisms - cell volume, cell dimensions, and carbon content. Rates and ancillary information were measured in units that meet the experimenter need, creating a need to harmonize the data units after collection. In addition different units can link to different mechanisms (carbon to nutritive quality of the prey, volume to size limits). As a result, grazing rates are thus available as pg C dinoflagellate-1 h-1, µm3 dinoflagellate-1 h-1 and prey cell dinoflagellate-1 h-1; clearance rate was calculated if not given and growth rate is expressed as the growth rate per day.

Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2013S1, deployed during Antarctic Fast Ice Network 2013 (AFIN 2013)

Snow height was measured by the Snow Depth Buoy 2013S1, an autonomous platform, installed close to Neumayer III Base, Antarctic during Antarctic Fast Ice Network 2013 (AFIN 2013). The resulting time series describes the evolution of snow height as a function of place and time between 2013-02-11 and 2013-04-29 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The buoy was installed on the ice shelf. In addition to snow height, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Negative values of snow height occur if surface ablation continues into the sea ice. Thus, these measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). Records without any snow height may still be used for sea ice drift analyses. Note: This data set contains only relative changes in snow height, because no initial readings of absolute snow height are available.

Data base on pelagic ciliates grazing and growth rate as a function of size, prey size and type compiled from the literature

Microzooplankton (the 20 to 200 µm size class of zooplankton) is recognised as an important part of marine pelagic ecosystems. In terms of biomass and abundance pelagic ciliates are one of the important groups of organism in microzooplankton. However, their rates - grazing and growth - , feeding behaviour and prey preferences are poorly known and understood. A set of data was assembled in order to derive a better understanding of pelagic ciliates rates, in response to parameters such as prey concentration, prey type (size and species), temperature and their own size. With these objectives, literature was searched for laboratory experiments with information on one or more of these parameters effect studied. The criteria for selection and inclusion in the database included: (i) controlled laboratory experiment with a known ciliates feeding on a known prey; (ii) presence of ancillary information about experimental conditions, used organisms - cell volume, cell dimensions, and carbon content. Rates and ancillary information were measured in units that meet the experimenter need, creating a need to harmonize the data units after collection. In addition different units can link to different mechanisms (carbon to nutritive quality of the prey, volume to size limits). As a result, grazing rates are thus available as pg C/(ciliate*h), µm**3/(ciliate*h) and prey cell/(ciliate*h); clearance rate was calculated if not given and growth rate is expressed as the growth rate per day.

Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2014S24, deployed during Antarctic Fast Ice Network 2014 (AFIN 2014)

Snow height was measured by the Snow Depth Buoy 2014S24, an autonomous platform, installed close to Neumayer III Base, Antarctic during Antarctic Fast Ice Network 2014 (AFIN 2014). The resulting time series describes the evolution of snow depth as a function of place and time between 2014-03-07 and 2014-05-16 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The buoy was installed on the ice shelf. In addition to snow depth, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Negative values of snow depth occur if surface ablation continues into the sea ice. Thus, these measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). Records without any snow depth may still be used for sea ice drift analyses. Note: This data set contains only relative changes in snow depth, because no initial readings of absolute snow depth are available.

Function of C1A barley peptidases and their inhibitors (cystatins) in barley seed germination : Funcion de las peptidasas C1A de cebada y sus inhibidores (cistatinas) en la germinacion de la semilla

Plant proteolysis is a metabolic process where specific enzymes called peptidases degrade proteins. In plants, this complex process involves broad metabolic networks and different sub-cellular compartments. Several types of peptidases take part in the proteolytic process, mainly cysteine-, serine-, aspartyl- and metallo- peptidases. Among the cysteine-peptidases, the papain-like or C1A peptidases (family C1, clan CA) are extensively present in land plants and are classified into catepsins L-, B-, H- and Flike. The catalytic mechanism of these C1A peptidases is highly conserved and involves the three amino acids Cys, His and Asn in the catalytic triad, and a Gln residue which seems essential for maintaining an active enzyme conformation. These proteins are synthesized as inactive precursors, which comprise an N-terminal signal peptide, a propeptide, and the mature protein. In barley, we have identified 33 cysteine-peptidases from the papain-like family, classifying them into 8 different groups. Five of them corresponded to cathepsins L-like (5 subgroups), 1 cathepsin B-like group, 1 cathepsin F-like group and 1 cathepsin H-like group. Besides, C1A peptidases are the specific targets of the plant proteinaceous inhibitors known as phytocystatins (PhyCys). The cystatin inhibitory mechanism is produced by a tight and reversible interaction with their target enzymes. In barley, the cystatin gene family is comprised by 13 members. In this work we have tried to elucidate the role of the C1A cysteine-peptidases and their specific inhibitors (cystatins) in the germination process of the barley grain. Therefore, we selected a representative member of each group/subgroup of C1A peptidases (1 cathepsin B-like, 1 cathepsin F-like, 1 cathepsin H-like and 5 cathepsins L-like). The molecular characterization of the cysteine-peptidases was done and the peptidase-inhibitor interaction was analyzed in vitro and in vivo. A study in the structural basis for specificity of pro-peptide/enzyme interaction in barley C1A cysteine-peptidases has been also carried out by inhibitory assays and the modeling of the three-dimensional structures. The barley grain maturation produces the accumulation of storage proteins (prolamins) in the endosperm which are mobilized during germination to supply the required nutrients until the photosynthesis is fully established. In this work, we have demonstrated the participation of the cysteine-peptidases and their inhibitors in the degradation of the different storage protein fractions (hordeins, albumins and globulins) present in the barley grain. Besides, transgenic barley plants overexpressing or silencing cysteine-peptidases or cystatins were obtained by Agrobacterium-mediated transformation of barley immature embryos to analyze their physiological function in vivo. Preliminary assays were carried out with the T1 grains of several transgenic lines. Comparing the knock-out and the overexpressing lines with the WT, alterations in the germination process were detected and were correlated with their grain hordein content. These data will be validated with the homozygous grains that are being produced through the double haploid technique by microspore culture. Resumen La proteólisis es un proceso metabólico por el cual se lleva a cabo la degradación de las proteínas de un organismo a través de enzimas específicas llamadas proteasas. En plantas, este complejo proceso comprende un entramado de rutas metabólicas que implican, además, diferentes compartimentos subcelulares. En la proteólisis participan numerosas proteasas, principalmente cisteín-, serín-, aspartil-, y metalo-proteasas. Dentro de las cisteín-proteasas, las proteasas tipo papaína o C1A (familia C1, clan CA) están extensamente representadas en plantas terrestres, y se clasifican en catepsinas tipo L, B, H y F. El mecanismo catalítico de estas proteasas está altamente conservado y la triada catalítica formada por los aminoácidos Cys, His y Asn, y a un aminoácido Gln, que parece esencial para el mantenimiento de la conformación activa de la proteína. Las proteasas C1A se sintetizan como precursores inactivos y comprenden un péptido señal en el extremo N-terminal, un pro-péptido y la proteína madura. En cebada hemos identificado 33 cisteín-proteasas de tipo papaína y las hemos clasificado filogenéticamente en 8 grupos diferentes. Cinco de ellos pertenecen a las catepsinas tipo L (5 subgrupos), un grupo a las catepsinas tipo-B, otro a las catepsinas tipo-F y un último a las catepsinas tipo-H. Las proteasas C1A son además las dianas específicas de los inhibidores protéicos de plantas denominados fitocistatinas. El mecanismo de inhibición de las cistatinas está basado en una fuerte interacción reversible. En cebada, se conoce la familia génica completa de las cistatinas, que está formada por 13 miembros. En el presente trabajo se ha investigado el papel de las cisteín-proteasas de cebada y sus inhibidores específicos en el proceso de la germinación de la semilla. Para ello, se seleccionó una proteasa representante de cada grupo/subgrupo (1 catepsina tipo- B, 1 tipo-F, 1 tipo-H, y 5 tipo-L, una por cada subgrupo). Se ha llevado a cabo su caracterización molecular y se ha analizado la interacción enzima-inhibidor tanto in vivo como in vitro. También se han realizado estudios sobre las bases estructurales que demuestran la especificidad en la interacción enzima/propéptido en las proteasas C1A de cebada, mediante ensayos de inhibición y la predicción de modelos estructurales de la interacción. Finalmente, y dado que durante la maduración de la semilla se almacenan proteínas de reserva (prolaminas) en el endospermo que son movilizadas durante la germinación para suministrar los nutrientes necesarios hasta que la nueva planta pueda realizar la fotosíntesis, en este trabajo se ha demostrado la participación de las cisteínproteasas y sus inhibidores en la degradación de las diferentes tipos de proteínas de reserva (hordeinas, albúmins y globulinas) presentes en el grano de cebada. Además, se han obtenido plantas transgénicas de cebada que sobre-expresan o silencian cistatinas y cisteín-proteasas con el fin de analizar la función fisiológica in vivo. Se han realizado análisis preliminares en las semillas T1 de varias líneas tránsgenicas de cebada y al comparar las líneas knock-out y las líneas de sobre-expresión con las silvestres, se han detectado alteraciones en la germinación que están además correlacionadas con el contenido de hordeinas de las semillas. Estos datos serán validados en las semillas homocigotas que se están generando mediante la técnica de dobles haploides a partir del cultivo de microesporas.

Structure function and Multifractal spectrum applied to Digital Elevation Model

A Digital Elevation Model (DEM) provides the information basis used for many geographic applications such as topographic and geomorphologic studies, landscape through GIS (Geographic Information Systems) among others. The DEM capacity to represent Earth?s surface depends on the surface roughness and the resolution used. Each DEM pixel depends on the scale used characterized by two variables: resolution and extension of the area studied. DEMs can vary in resolution and accuracy by the production method, although there are statistical characteristics that keep constant or very similar in a wide range of scales. Based on this property, several techniques have been applied to characterize DEM through multiscale analysis directly related to fractal geometry: multifractal spectrum and the structure function. The comparison of the results by both methods is discussed. The study area is represented by a 1024 x 1024 data matrix obtained from a DEM with a resolution of 10 x 10 m each point, which correspond with a region known as ?Monte de El Pardo? a property of Spanish National Heritage (Patrimonio Nacional Español) of 15820 Ha located to a short distance from the center of Madrid. Manzanares River goes through this area from North to South. In the southern area a reservoir is found with a capacity of 43 hm3, with an altitude of 603.3 m till 632 m when it is at the highest capacity. In the middle of the reservoir the minimum altitude of this area is achieved.

Unveiling Protein Functions through the Dynamics of the Interaction Network

Protein interaction networks have become a tool to study biological processes, either for predicting molecular functions or for designing proper new drugs to regulate the main biological interactions. Furthermore, such networks are known to be organized in sub-networks of proteins contributing to the same cellular function. However, the protein function prediction is not accurate and each protein has traditionally been assigned to only one function by the network formalism. By considering the network of the physical interactions between proteins of the yeast together with a manual and single functional classification scheme, we introduce a method able to reveal important information on protein function, at both micro- and macro-scale. In particular, the inspection of the properties of oscillatory dynamics on top of the protein interaction network leads to the identification of misclassification problems in protein function assignments, as well as to unveil correct identification of protein functions. We also demonstrate that our approach can give a network representation of the meta-organization of biological processes by unraveling the interactions between different functional classes

Topology of a percolating soil pore network.

A connectivity function defined by the 3D-Euler number, is a topological indicator and can be related to hydraulic properties (Vogel and Roth, 2001). This study aims to develop connectivity Euler indexes as indicators of the ability of soils for fluid percolation. The starting point was a 3D grey image acquired by X-ray computed tomography of a soil at bulk density of 1.2 mg cm-3. This image was used in the simulation of 40000 particles following a directed random walk algorithms with 7 binarization thresholds. These data consisted of 7 files containing the simulated end points of the 40000 random walks, obtained in Ruiz-Ramos et al. (2010). MATLAB software was used for computing the frequency matrix of the number of particles arriving at every end point of the random walks and their 3D representation.

The Phase Space as a New Representation of the Dynamical Behaviour of Temperature and Enthalpy in a Reefer monitored with a Multidistributed Sensors Network

The study of temperature gradients in cold stores and containers is a critical issue in the food industry for the quality assurance of products during transport, as well as forminimizing losses. The objective of this work is to develop a new methodology of data analysis based on phase space graphs of temperature and enthalpy, collected by means of multidistributed, low cost and autonomous wireless sensors and loggers. A transoceanic refrigerated transport of lemons in a reefer container ship from Montevideo (Uruguay) to Cartagena (Spain) was monitored with a network of 39 semi-passive TurboTag RFID loggers and 13 i-button loggers. Transport included intermodal transit from transoceanic to short shipping vessels and a truck trip. Data analysis is carried out using qualitative phase diagrams computed on the basis of Takens?Ruelle reconstruction of attractors. Fruit stress is quantified in terms of the phase diagram area which characterizes the cyclic behaviour of temperature. Areas within the enthalpy phase diagram computed for the short sea shipping transport were 5 times higher than those computed for the long sea shipping, with coefficients of variation above 100% for both periods. This new methodology for data analysis highlights the significant heterogeneity of thermohygrometric conditions at different locations in the container.

Protein isolate and inulin: chemical, rheological and structural basis

Soy protein isolate is typical vegetable protein with health-enhancing activities. Inulin, a prebiotic no digestible carbohydrate, has functional properties. A mashed potato serving of 200 g with added soy protein isolate and inulin concentrations of 15?60 g kg provides from 3 to 12 g of soy protein isolate and/or inulin, respectively. Currently, no information is available about the possible texture-modifying effect of this non-ionizable polar carbohydrate in different soy-based food systems. In this study, the effect of the addition of soy protein isolate and inulin blends at different soy protein isolate: inulin ratios on the degree of inulin polymerization and the rheological and structural properties of fresh mashed and frozen/thawed mashed potatoes were evaluated. The inulin chemical structure remained intact throughout the various treatments, and soy protein isolate did not affect inulin composition being a protein compatible with this fructan. Small-strain rheology showed that both ingredients behaved like soft fillers. In the frozen/thawed mashed potatoes samples,0 addition of 30 : 30 and 15 : 60 blend ratios significantly increased elasticity (G value) compared with 0 : 0 control, consequently reducing the freeze/thaw stability conferred by the cryoprotectants. Inulin crystallites caused a significant strengthening effect on soy protein isolate gel. Micrographs revealed that soy protein isolate supports the inulin structure by building up a second fine-stranded network. Thereby, possibility of using soy protein isolate and inulin in combination with mashed potatoes to provide a highly nutritious and healthy product is promising.

Algoritmos para mineria de datos con redes de neuronas

Este trabajo propone una serie de algoritmos con el objetivo de extraer información de conjuntos de datos con redes de neuronas. Se estudian dichos algoritmos con redes de neuronas Enhenced Neural Networks (ENN), debido a que esta arquitectura tiene algunas ventajas cuando se aproximan funciones mediante redes neuronales. En la red ENN los pesos de la matriz principal varián con cada patrón, por lo que se comete un error menor en la aproximación. Las redes de neuronas ENN reúnen la información en los pesos de su red auxiliar, se propone un método para obtener información de la red a través de dichos pesos en formas de reglas y asignando un factor de certeza de dichas reglas. La red ENN obtiene un error cuadrático medio menor que el error teórico de una aproximación matemática por ejemplo mediante polinomios de Taylor. Se muestra como una red ENN, entrenada a partir un conjunto de patrones obtenido de una función de variables reales, sus pesos asociados tienen unas relaciones similares a las que se veri_can con las variables independientes con dicha función de variables reales. Las redes de neuronas ENN aproximan polinomios, se extrae conocimiento de un conjunto de datos de forma similar a la regresión estadística, resolviendo de forma más adecuada el problema de multicolionalidad en caso de existir. Las relaciones a partir de los pesos asociados de la matriz de la red auxiliar se obtienen similares a los coeficientes de una regresión para el mismo conjunto numérico. Una red ENN entrenada a partir de un conjunto de datos de una función boolena extrae el conocimiento a partir de los pesos asociados, y la influencia de las variables de la regla lógica de la función booleana, queda reejada en esos pesos asociados a la red auxiliar de la red ENN. Se plantea una red de base radial (RBF) para la clasificación y predicción en problemas forestales y agrícolas, obteniendo mejores resultados que con el modelo de regresión y otros métodos. Los resultados con una red RBF mejoran al método de regresión si existe colinealidad entre los datos que se dispone y no son muy numerosos. También se detecta que variables tienen más importancia en virtud de la variable pronóstico. Obteniendo el error cuadrático medio con redes RBF menor que con otros métodos, en particular que con el modelo de regresión. Abstract A series of algorithms is proposed in this study aiming at the goal of producing information about data groups with a neural network. These algorithms are studied with Enheced Neural Networks (ENN), owing to the fact that this structure shows sever advantages when the functions are approximated by neural networks. Main matrix weights in th ENN vary on each pattern; so, a smaller error is produced when approximating. The neural network ENN joins the weight information contained in their auxiliary network. Thus, a method to obtain information on the network through those weights is proposed by means of rules adding a certainty factor. The net ENN obtains a mean squared error smaller than the theorical one emerging from a mathematical aproximation such as, for example, by means of Taylor's polynomials. This study also shows how in a neural network ENN trained from a set of patterns obtained through a function of real variables, its associated weights have relationships similar to those ones tested by means of the independent variables connected with such functions of real variables. The neural network ENN approximates polynomials through it information about a set of data may be obtained in a similar way than through statistical regression, solving in this way possible problems of multicollinearity in a more suitable way. Relationships emerging from the associated weights in the auxiliary network matrix obtained are similar to the coeficients corresponding to a regression for the same numerical set. A net ENN trained from a boolean function data set obtains its information from its associated weights. The inuence of the variables of the boolean function logical rule are reected on those weights associated to the net auxiliar of the ENN. A radial basis neural networks (RBF) for the classification and prediction of forest and agricultural problems is proposed. This scheme obtains better results than the ones obtained by means of regression and other methods. The outputs with a net RBF better the regression method if the collineality with the available data and their amount is not very large. Detection of which variables are more important basing on the forecast variable can also be achieved, obtaining a mean squared error smaller that the ones obtained through other methods, in special the one produced by the regression pattern.

Optimisation of algorithms to compute information theoretic indexes

Analysis of big amount of data is a field with many years of research. It is centred in getting significant values, to make it easier to understand and interpret data. Being the analysis of interdependence between time series an important field of research, mainly as a result of advances in the characterization of dynamical systems from the signals they produce. In the medicine sphere, it is easy to find many researches that try to understand the brain behaviour, its operation mode and its internal connections. The human brain comprises approximately 1011 neurons, each of which makes about 103 synaptic connections. This huge number of connections between individual processing elements provides the fundamental substrate for neuronal ensembles to become transiently synchronized or functionally connected. A similar complex network configuration and dynamics can also be found at the macroscopic scales of systems neuroscience and brain imaging. The emergence of dynamically coupled cell assemblies represents the neurophysiological substrate for cognitive function such as perception, learning, thinking. Understanding the complex network organization of the brain on the basis of neuroimaging data represents one of the most impervious challenges for systems neuroscience. Brain connectivity is an elusive concept that refers to diferent interrelated aspects of brain organization: structural, functional connectivity (FC) and efective connectivity (EC). Structural connectivity refers to a network of physical connections linking sets of neurons, it is the anatomical structur of brain networks. However, FC refers to the statistical dependence between the signals stemming from two distinct units within a nervous system, while EC refers to the causal interactions between them. This research opens the door to try to resolve diseases related with the brain, like Parkinson’s disease, senile dementia, mild cognitive impairment, etc. One of the most important project associated with Alzheimer’s research and other diseases are enclosed in the European project called Blue Brain. The center for Biomedical Technology (CTB) of Universidad Politecnica de Madrid (UPM) forms part of the project. The CTB researches have developed a magnetoencephalography (MEG) data processing tool that allow to visualise and analyse data in an intuitive way. This tool receives the name of HERMES, and it is presented in this document. Analysis of big amount of data is a field with many years of research. It is centred in getting significant values, to make it easier to understand and interpret data. Being the analysis of interdependence between time series an important field of research, mainly as a result of advances in the characterization of dynamical systems from the signals they produce. In the medicine sphere, it is easy to find many researches that try to understand the brain behaviour, its operation mode and its internal connections. The human brain comprises approximately 1011 neurons, each of which makes about 103 synaptic connections. This huge number of connections between individual processing elements provides the fundamental substrate for neuronal ensembles to become transiently synchronized or functionally connected. A similar complex network configuration and dynamics can also be found at the macroscopic scales of systems neuroscience and brain imaging. The emergence of dynamically coupled cell assemblies represents the neurophysiological substrate for cognitive function such as perception, learning, thinking. Understanding the complex network organization of the brain on the basis of neuroimaging data represents one of the most impervious challenges for systems neuroscience. Brain connectivity is an elusive concept that refers to diferent interrelated aspects of brain organization: structural, functional connectivity (FC) and efective connectivity (EC). Structural connectivity refers to a network of physical connections linking sets of neurons, it is the anatomical structur of brain networks. However, FC refers to the statistical dependence between the signals stemming from two distinct units within a nervous system, while EC refers to the causal interactions between them. This research opens the door to try to resolve diseases related with the brain, like Parkinson’s disease, senile dementia, mild cognitive impairment, etc. One of the most important project associated with Alzheimer’s research and other diseases are enclosed in the European project called Blue Brain. The center for Biomedical Technology (CTB) of Universidad Politecnica de Madrid (UPM) forms part of the project. The CTB researches have developed a magnetoencephalography (MEG) data processing tool that allow to visualise and analyse data in an intuitive way. This tool receives the name of HERMES, and it is presented in this document.