Mapping and characterization of the Xlsirt P11 RNA cis-acting localization element

In many organisms, polarity of the oocyte is established post-transcriptionally via subcellular RNA localization. Many RNAs are localized during oogenesis in Xenopus laevis, including Xlsirts ( Xenopus laevis short interspersed repeat transcripts) [Kloc, 1993]. Xlsirts constitute a large family defined by highly homologous repeat units 79–81 nucleotides in length. Endogenous Xlsirt RNAs use the METRO (Message Transport Organizer) pathway of localization, where RNAs are transported from the nucleus to the mitochondrial cloud in stage I oocytes. Secondly, RNAs anchor at the vegetal pole in stage II oocytes. Exogenous Xlsirt RNAs can also utilize the Late pathway of localization, which involves localization to the vegetal cortex during stage III of oogenesis and results in RNAs anchored in the cortex of the entire vegetal hemisphere. ^ The Xlsirts localization signal is contained within the repeat region. This study was designed to test the hypothesis that there are cis -acting localization elements in Xlsirts, and that higher order structure plays a role. Results of experiments on Xlsirt P11, a 1700 basepair (bp) family member, led to the conclusion that a 137-bp fragment of the repetitive region is necessary and sufficient for METRO and Late pathway localization. This analysis definitively demonstrates that the Xlsirt localization signal for the METRO and Late pathways reside within the repetitive region and not within the flanking regions. Analysis of Xlsirt linker scanning mutations revealed two METRO-pathway specific subelements, and one Late-pathway specific subelement. Functional, computer, and biochemical evidence relates the higher order structure of this element to its ability to function as a localization element. ^ Xlsirt 137 is 99% identical to the Xlsirt consensus sequence identified in this study, suggesting that it is the localization element for all localized Xlsirt family members. The repeat unit was reframed based on function, rather than arbitrarily based on sequence. This work supports the hypothesis presented in 1981 by George Spohr, who originally isolated the Xlsirts, which stated that the highly conserved repetitive elements must be constrained from variability due to some unknown function of the repeats themselves. These studies shed light on the mechanism of RNA localization, linking structure and function. ^

RNA polymerase II large subunit degradation induced by ecteinascidin 743: Molecular characterization and subsequent rational investigation of antitumor mechanisms in cancers with clinical response

Ecteinascidin 743 (Et-743), which is a novel DNA minor groove alkylator with a unique spectrum of antitumor activity, is currently being evaluated in phase II/III clinical trials. Although the precise molecular mechanisms responsible for the observed antitumor activity are poorly understood, recent data suggests that post-translational modifications of RNA polymerase II Large Subunit (RNAPII LS) may play a central role in the cellular response to this promising anticancer agent. The stalling of an actively transcribing RNAPII LS at Et-743-DNA adducts is the initial cellular signal for transcription-coupled nucleotide excision repair (TC-NER). In this manner, Et-743 poisons TC-NER and produces DNA single strand breaks. Et-743 also inhibits the transcription and RNAPII LS-mediated expression of selected genes. Because the poisoning of TC-NER and transcription inhibition are critical components of the molecular response to Et-743 treatment, we have investigated if changes in RNAPII LS contribute to the disruption of these two cellular pathways. In addition, we have studied changes in RNAPII LS in two tumors for which clinical responses were reported in phase I/II clinical trials: renal cell carcinoma and Ewing's sarcoma. Our results demonstrate that Et-743 induces degradation of the RNAPII LS that is dependent on active transcription, a functional 26S proteasome, and requires functional TC-NER, but not global genome repair. Additionally, we have provided the first experimental data indicating that degradation of RNAPII LS might lead to the inhibition of activated gene transcription. A set of studies performed in isogenic renal carcinoma cells deficient in von Hippel-Lindau protein, which is a ubiquitin-E3-ligase for RNAPII LS, confirmed the central role of RNAPII LS degradation in the sensitivity to Et-743. Finally, we have shown that RNAPII LS is also degraded in Ewing's sarcoma tumors following Et-743 treatment and provide data to suggest that this event plays a role in decreased expression of the Ewing's sarcoma oncoprotein, EWS-Fli1. Altogether, these data implicate degradation of RNAPII LS as a critical event following Et-743 exposure and suggest that the clinical activity observed in renal carcinoma and Ewing's sarcoma may be mediated by disruption of molecular pathways requiring a fully functional RNAPII LS. ^

Biochemical and genetic characterization of the Saccharomyces cerevisiae Hsp110 molecular chaperone Sse1

The Ssel/Hsp110 molecular chaperones are a poorly understood subgroup of the Hsp70 chaperone family. Hsp70 can refold denatured polypeptides via a carboxyl-terminal peptide binding domain (PBD), which is regulated by nucleotide cycling in an amino-terminal ATPase domain. However, unlike Hsp70, both Sse1 and mammalian Hsp110 bind unfolded peptide substrates but cannot refold them. To test the in vivo requirement for interdomain communication, SSE1 alleles carrying amino acid substitutions in the ATPase domain were assayed for their ability to complement sse1Δ phenotypes. Surprisingly, all mutants predicted to abolish ATP hydrolysis complemented the temperature sensitivity of sse1Δ, whereas mutations in predicted ATP binding residues were non-functional. Remarkably, the two domains of Ssel when expressed in trans functionally complement the sse1Δ growth phenotype and interact by coimmunoprecipitation analysis, indicative of a novel type of interdomain communication. ^ Relatively little is known regarding the interactions and cellular functions of Ssel. Through co-immunoprecipitation analysis, we found that Ssel forms heterodimeric complexes with the abundant cytosolic Hsp70s Ssa and Ssb in vivo. Furthermore, these complexes can be efficiently reconstituted in vitro using purified proteins. The ATPase domains of Ssel and the Hsp70s were found to be critical for interaction as inactivating point mutations severely reduced interaction efficiency. Ssel stimulated Ssal ATPase activity synergistically with the co-chaperone Ydj1 via a novel nucleotide exchange activity. Furthermore, FES1, another Ssa nucleotide exchange factor, can functionally substitute for SSE1/2 when overexpressed, suggesting that Hsp70 nucleotide exchange is the fundamental role of the Sse proteins in yeast, and by extension, the Hsp110 homologs in mammals. ^ Cells lacking SSE1 were found to accumulate prepro-α-factor, but not the cotranslationally imported protein Kar2, similar to mutants in the Ssa chaperones. This indicates that the interaction between Ssel and Ssa is functionally significant in vivo. In addition, sse10 cells are compromised for cell wall strength, likely a result of decreased Hsp90 chaperone activity with the cell integrity MAP kinase SIC. Taken together, this work established that the Hsp110 family must be considered an essential component of Hsp70 chaperone biology in the eukaryotic cell.^

Identification and characterization of distinct classes of cAMP-pulse-induced genes in Dictyostelium discoideum development

The unicellular amoeba Dictyostelium discoideum embarks on a developmental program upon starvation. During development, extracellular oscillatory cAMP signaling orchestrates the chemotaxis-mediated aggregation of ∼105 amoebae and is required for optimal induction of so-called pulse-induced genes. This requirement for pulsatile CAMP reflects adaptation of the cAMP-receptor-mediated pathways that regulate these genes. Through examination of a collection of pulse-induced genes, we defined two distinct gene classes based on their induction kinetics and the impact of mutations that impair PKA signaling. The first class (represented by D2 and prtA) is highly dependent on PKA signaling, whereas the second class (represented by carA, gpaB, and acaA) is not. Analysis of expression kinetics revealed that these classes are sequentially expressed with the PKA-independent genes peaking in expression before the PKA-dependent class. Experiments with cycloheximide, an inhibitor of translation, demonstrated that the pulse induction of both classes depends on new protein synthesis early in development. carA and gpaB also exhibit pulse-independent, starvation-induced expression which, unlike their pulse induction, was found to be insensitive to cycloheximide added at the outset of starvation. This result indicates that the mechanism of starvation induction pre-exists in growing cells and is distinct from the pulse induction mechanism for these genes. In order to identify cis-acting elements that are critical for induction of carA, we constructed a GFP reporter controlled by a 914-base-pair portion of its promoter and verified that its expression was PKA-independent, pulse-inducible, and developmentally regulated like the endogenous carA gene. By a combination of truncation, internal deletion, and site-directed mutation, we defined several distinct functional elements within the carA promoter, including a 39-bp region required for pulse induction between base pairs -321 and -282 (relative to the transcription start site), a 131-bp region proximal to the start site that is sufficient for starvation induction, and two separate enhancer domains. Identification of factors that interact with these promoter elements and genetic approaches exploiting the GFP reporter described here should help complete our understanding of the mechanisms regulating these genes, including adaptation mechanisms that likely also govern chemotaxis of Dictyostelium and mammalian cells. ^

Characterization of tumor-associated Foxp3+ regulatory T cells and de-novo induction by the tumor microenvironment

Regulatory T cells expressing the fork-head box transcription factor 3 (Foxp3) play a central role in the dominant control of immunological tolerance. Compelling evidence obtained from both animal and clinical studies have now linked the expansion and accumulation of Foxp3+ regulatory T cells associated with tumor lesions to the failure of immune-mediated tumor rejection. However, further progress of the field is hampered by the gap of knowledge regarding their phenotypic, functional, and the developmental origins in which these tumor-associated Foxp3+ regulatory T cells are derived. Here, we have characterized the general properties of tumor-associated Foxp3+ regulatory T cells and addressed the issue of tumor microenvironment mediated de-novo induction by utilizing a well known murine tumor model MCA-205 in combination with our BAC Foxp3-GFP reporter mice and OT-II TCR transgenic mice on the RAG deficient background (RAG OT-II). De-novo induction defines a distinct mechanism of converting non-regulatory precursor cells to Foxp3+ regulatory T cells in the periphery as opposed to the expansion of pre-existing regulatory T cells formed naturally during thymic T cell development. This mechanism is of particularly importance to how tumors induce tumor-antigen-specific suppressor cells to subvert anti-tumor immune responses. Our study has found that tumor-associated Foxp3+ regulatory T cells are highly activated, undergo vigorous proliferation, are more potent by in-vitro suppression assays, and express higher levels of membrane-bound TGF-β1 than non-tumor regulatory T cells. With Foxp3-GFP reporter mice or RAG OT-II TCR transgenic mice, we show that tumor tissue can induce detectable de-novo generation of Foxp3+ regulatory T cells of both polyclonal or antigen specific naïve T cells. This process was not only limited for subcutaneous tumors but for lung tumors as well. Furthermore, this process required the inducing antigen to be co-localized within the tumor tissue. Examination of tumor tissue revealed an abundance of myeloid CD11b+ antigen-presenting cells that were capable of inducing Foxp3+ regulatory T cells. Taken together, these findings elucidate the general attributes and origins of tumor-associated Foxp3+ regulatory T cells in the tumor microenvironment and in their role in the negative regulation of tumor immunity.^

IDENTIFICATION AND CHARACTERIZATION OF DISTINCT POPULATIONS OF CLONOGENIC BONE MARROW STROMAL CELLS CAPABLE OF TRANSFERRING THE HEMATOPOIETIC MICROENVIRONMENT IN VIVO AND SUPPORTING LT-HSCs IN VITRO

Bone marrow (BM) stromal cells are ascribed two key functions, 1) stem cells for non-hematopoietic tissues (MSC) and 2) as components of the hematopoietic stem cell niche. Current approaches studying the stromal cell system in the mouse are complicated by the low yield of clonogenic progenitors (CFU-F). Given the perivascular location of MSC in BM, we developed an alternative methodology to isolate MSC from mBM. An intact ‘plug’ of bone marrow is expelled from bones and enzymatically disaggregated to yield a single cell suspension. The recovery of CFU-F (1917.95+199) reproducibly exceeds that obtained using the standard BM flushing technique (14.32+1.9) by at least 2 orders of magnitude (P<0.001; N = 8) with an accompanying 196-fold enrichment of CFU-F frequency. Purified BM stromal and vascular endothelial cell populations are readily obtained by FACS. A detailed immunophenotypic analysis of lineage depleted BM identified PDGFRαβPOS stromal cell subpopulations distinguished by their expression of CD105. Both subpopulations retained their original phenotype of CD105 expression in culture and demonstrate MSC properties of multi-lineage differentiation and the ability to transfer the hematopoietic microenvironment in vivo. To determine the capacity of either subpopulation to support long-term multi-lineage reconstituting HSCs, we fractionated BM stromal cells into either the LinNEGPDGFRαβPOSCD105POS and LINNEGPDGFRαβPOSCD105LOW/- populations and tested their capacity to support LT-HSC by co-culturing each population with either 1 or 10 HSCs for 10 days. Following the 10 day co-culture period, both populations supported transplantable HSCs from 10 cells/well co-cultures demonstrating high levels of donor repopulation with an average of 65+23.6% chimerism from CD105POS co-cultures and 49.3+19.5% chimerism from the CD105NEG co-cultures. However, we observed a significant difference when mice were transplanted with the progeny of a single co-cultured HSC. In these experiments, CD105POS co-cultures (100%) demonstrated long-term multi- lineage reconstitution, while only 4 of 8 mice (50%) from CD105NEG -single HSC co-cultures demonstrated long-term reconstitution, suggesting a more limited expansion of functional stem cells. Taken together, these results demonstrate that the PDGFRαβCD105POS stromal cell subpopulation is distinguished by a unique capacity to support the expansion of long-term reconstituting HSCs in vitro.

BIOCHEMICAL CHARACTERIZATION OF BINDING PARTNERS OF TWO HSP70 CO-CHAPERONES IN SACCHAROMYCES CEREVISIAE

Cells are exposed to a variety of environmental and physiological changes including temperature, pH and nutrient availability. These changes cause stress to cells, which results in protein misfolding and altered cellular protein homeostasis. How proteins fold into their three-dimensional functional structure is a fundamental biological process with important relevance to human health. Misfolded and aggregated proteins are linked to multiple neurodegenerative diseases, cardiovascular disease and cystic fibrosis. To combat proteotoxic stress, cells deploy an array of molecular chaperones that assist in the repair or removal of misfolded proteins. Hsp70, an evolutionarily conserved molecular chaperone, promotes protein folding and helps maintain them in a functional state. Requisite co-chaperones, including nucleotide exchange factors (NEFs) strictly regulate and serve to recruit Hsp70 to distinct cellular processes or locations. In yeast and human cells, three structurally non-related cytosolic NEFs are present: Sse1 (Hsp110), Fes1 (HspBP1) and Snl1 (Bag-1). Snl1 is unique among the cytosolic NEFs as it is localized at the ER membrane with its Hsp70 binding (BAG) domain exposed to the cytosol. I discovered that Snl1 distinctly interacts with assembled ribosomes and several lines of evidence indicate that this interaction is both independent of and concurrent with binding to Hsp70 and is not dependent on membrane localization. The ribosome-binding site is identified as a short lysine-rich motif within the amino terminus of the Snl1 BAG domain distinct from the Hsp70 interaction region. In addition, I demonstrate ribosome association with the Snl1 homolog in the pathogenic fungus, Candida albicans and localize this putative NEF to a perinuclear/ER membrane, suggesting functional conservation in fungal BAG domain-containing proteins. As a first step in determining specific domain architecture in fungal BAG proteins, I present the preliminary steps of protein purification and analysis of the minimal Hsp70 binding region in in both S.cerevisiae and C. albicans Snl1. Contrary to previous in vitro evidence which showed the Fes1 NEF to interact with both cytosolic Hsp70s, Ssa and Ssb, Fes1 is shown to interact specifically with Ssa when expressed under normal cellular conditions in S. cerevisiae. This is the first reported evidence of Hsp70 binding selectivity for a cytosolic NEF, and suggests a possible mechanism to achieve specificity in Hsp70-dependent functions. Taken together, the work presented in this dissertation highlights the striking divergence among Hsp70 co-chaperones in selecting binding partners, which may correlate with their specific roles in the cell.

Cloning and characterization of the full length cDNA and promoter of mouse tissue transglutaminase

Transglutaminases are a family of enzymes that catalyze the covalent cross-linking of proteins through the formation of $\varepsilon$-($\gamma$-glutaminyl)-lysyl isopeptide bonds. Tissue transglutaminase (Tgase) is an intracellular enzyme which is expressed in terminally differentiated and senescent cells and also in cells undergoing apoptotic cell death. To characterize this enzyme and examine its relationship with other members of the transglutaminase family, cDNAs, the first two exons of the gene and 2 kb of the 5$\sp\prime$ flanking region, including the promoter, were isolated. The full length Tgase transcript consists of 66 bp of 5$\sp\prime$-UTR (untranslated) sequence, an open reading frame which encodes 686 amino acids and 1400 bp of 3$\sp\prime$-UTR sequence. Alignment of the deduced Tgase protein sequence with that of other transglutaminases revealed regions of strong homology, particularly in the active site region.^ The Tgase cDNA was used to isolate and characterize a genomic clone encompassing the 5$\sp\prime$ end of the mouse Tgase gene. The transcription start site was defined using genomic and cDNA clones coupled with S1 protection analysis and anchored PCR. This clone includes 2.3 kb upstream of the transcription start site and two exons that contain the first 256 nucleotides of the mouse Tgase cDNA sequence. The exon intron boundaries have been mapped and compared with the exon intron boundaries of three members of the transglutaminase family: human factor XIIIa, the human keratinocyte transglutaminase and human erythrocyte band 4.1. Tissue Tgase exon II is similar to comparable exons of these genes. However, exon I bears no resemblance with any of the other transglutaminase amino terminus exons.^ Previous work in our laboratory has shown that the transcription of the Tgase gene is directly controlled by retinoic acid and retinoic acid receptors. To identify the region of the Tgase gene responsible for regulating its expression, fragments of the Tgase promoter and 5$\sp\prime$-flanking region were cloned into the chloramphenicol actetyl transferase (CAT) reporter constructs. Transient transfection experiments with these constructs demonstrated that the upstream region of Tgase is a functional promoter which contains a retinoid response element within a 1573 nucleotide region spanning nucleotides $-$252 to $-$1825. ^

Mechanism of mRNA stability control mediated by AU -rich element: Characterization of cis-elements and trans-factors

Regulation of cytoplasmic deadenylation, the first step in mRNA turnover, has direct impact on the fate of gene expression. AU-rich elements (AREs) found in the 3′ untranslated regions of many labile mRNAs are the most common RNA-destabilizing elements known in mammalian cells. Based on their sequence features and functional properties, AREs can be divided into three classes. Class I or class III ARE directs synchronous deadenylation, whereas class II ARE directs asynchronous deadenylation with the formation of poly(A)-intermediates. Through systematic mutagenesis study, we found that a cluster of five or six copies of AUUUA motifs forming various degrees of reiteration is the key feature dictating the choice between asynchronous versus synchronous deadenylation. A 20–30 nt AU-rich sequence immediately 5 ′ to this cluster of AUUUA motifs can greatly enhance its destabilizing ability and is an integral part of the AREs. These two features are the defining characteristics of class II AREs. ^ To better understand the decay mechanism of AREs, current methods have several limitations. Taking the advantage of tetracycline-regulated promoter, we developed a new transcriptional pulse strategy, Tet-system. By controlling the time and the amount of Tet addition, a pulse of RNA could be generated. Using this new system, we showed that AREs function in both growth- and density-arrested cells. The new strategy offers for the first time an opportunity to investigate control of mRNA deadenylation and decay kinetics in mammalian cells that exhibit physiologically relevant conditions. ^ As a member of heterogeneous nuclear RNA-binding protein, hnRNP D 0/AUF1 displays specific affinities for ARE sequences in vitro . But its in vivo function in ARE-mediated mRNA decay is unclear. AUF1/hnRNP D0 is composed of at least four isoforms derived by alternative RNA splicing. Each isoform exhibits different affinity for ARE sequence in vitro. Here, we examined in vivo effect of AUF1s/hnRNP D0s on degradation of ARE-containing mRNA. Our results showed that all four isoforms exhibit various RNA stabilizing effects in NIH3T3 cells, which are positively correlated with their binding affinities for ARE sequences. Further experiments indicated that AUF1/hnRNP D0 has a general role in modulating the stability of cytoplasmic mRNAs in mammalian cells. ^

Isolation and characterization of a novel human Mix -like homeobox gene MIXL

Molecular events involved in specification of early hematopoietic system are not well known. In Xenopus, a paired-box homeodomain family (Mix.1–4) has been implicated in this process. Although Mix-like homeobox genes have been isolated from zebrafish (bon), chicken (CMIX) and mice (MmI/MIXL1), isolation of a human Mix-like gene has remained elusive. ^ We have recently isolated and characterized a novel human Mix-like homeobox gene with a predicted open reading frame of 232 amino acids designated the Mix.1 homeobox (Xenopus laevis)-like gene (MIXL). The overall identity of this novel protein to CMIX and MmI/MIXL1 is 41% and 69%, respectively. However, the identity in the homeodomain is 66% to that of Xenopus Mix.1, 79% to that of CMIX, and 94% to that of MmI/MIXL1. In normal hematopoiesis, MIXL expression appears to be restricted immature B and T lymphoid cells. Several acute leukemic cell lines of B, T and myeloid lineages express MIXL suggesting a survival/block in differentiation advantage. Furthermore, Xenopus animal cap assay revealed that MIXL could induce expression of the α-globin gene, suggesting a functional conservation of the homeodomain. ^ Biochemical analysis revealed that MIXL proteins are phosphorylated at multiple sites. Immunoprecipitation and immunoblotting confirmed that MIXL is tyrosine phosphorylated. Mutational analysis determined that Tyr20 appears to be the site for phosphorylation. However, deletion analysis preliminarily showed that the proline-rich domain appears not to be necessary for tyrosine phosphorylation. The novel finding will help us make a deeper understanding of the regulation on homeodomain proteins by rarely reported tyrosine phosphorylation. ^ Taken together, isolation of the MIXL gene is the first step toward understanding novel regulatory circuits in early hematopoietic differentiation and malignant transformation. ^

Plant Functional Type (PFT) map over Siberia derived from GlobCover 2005 dataset, link to dataset in NetCDF format

High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics.

HERMES: Towards an Integrated Toolbox to Characterize Functional and Effective Brain Connectivity

The analysis of the interdependence between time series has become an important field of research in the last years, mainly as a result of advances in the characterization of dynamical systems from the signals they produce, the introduction of concepts such as generalized and phase synchronization and the application of information theory to time series analysis. In neurophysiology, different analytical tools stemming from these concepts have added to the ‘traditional’ set of linear methods, which includes the cross-correlation and the coherency function in the time and frequency domain, respectively, or more elaborated tools such as Granger Causality. This increase in the number of approaches to tackle the existence of functional (FC) or effective connectivity (EC) between two (or among many) neural networks, along with the mathematical complexity of the corresponding time series analysis tools, makes it desirable to arrange them into a unified-easy-to-use software package. The goal is to allow neuroscientists, neurophysiologists and researchers from related fields to easily access and make use of these analysis methods from a single integrated toolbox. Here we present HERMES (http://hermes.ctb.upm.es), a toolbox for the Matlab® environment (The Mathworks, Inc), which is designed to study functional and effective brain connectivity from neurophysiological data such as multivariate EEG and/or MEG records. It includes also visualization tools and statistical methods to address the problem of multiple comparisons. We believe that this toolbox will be very helpful to all the researchers working in the emerging field of brain connectivity analysis.

Measurement and characterization of interoperability between devices with SpaceWire communications

This Master Final Project is intended to show the process developed to the functional and electrical characterization between different devices that use the SpaceWire space communications standard integrated into an evaluation board designed for this purpose. In order to carry out this characterization, firstly, a study to understand the SpaceWire standard is done. After that, another study for the understanding of the demonstration board with its different interfaces and IPs of SpW is done. According to this, it is expected to find out how the SpW devices are structured, especially at FPGA level, and how is the communication between them. Based on the knowledge obtained about SpaceWire and the SpW devices integrated into the evaluation board, the set of measurements and the strategy to validate electrical interoperability between the different devices are defined, as well as to perform functional checks required to ensure its proper understanding. Furthermore, it will let check whether the standard is met and search the limit of operation within a communication system representative of existing equipment in a satellite. Once finished the test plan and implemented on the representative hardware, the board will be considered characterized at SpW level and a report with the conclusions reached about the operation of the SpW interfaces in the board and constraints found will be done. RESUMEN. El presente Trabajo Fin de Máster pretende mostrar el proceso realizado para la caracterización eléctrica y funcional entre distintos dispositivos que utilizan el estándar de comunicaciones espaciales SpaceWire integrados en una tarjeta de evaluación diseñada para tal efecto. Para poder llevar a cabo dicha caracterización, en primer lugar, se realiza un estudio para el conocimiento del estándar SpaceWire. A continuación, se lleva a cabo otro estudio para el conocimiento de la tarjeta de demostración en la que se encuentran los distintos interfaces e IPs de SpW. Con esto último, se pretende conocer como están estructurados los dispositivos SpW, sobre todo a nivel de FPGA, y como se realiza la comunicación entre ellos. En base a los conocimientos adquiridos acerca de SpaceWire y los dispositivos SpW de la tarjeta de evaluación, se definen el conjunto de medidas y la estrategia a seguir para validar eléctricamente la interoperabilidad entre los distintos dispositivos, así como para realizar las comprobaciones funcionales necesarias para asegurar su correcto entendimiento. Además, con ello se podrá comprobar si se cumple el estándar y se podrá también buscar el límite de operación dentro de un sistema de comunicaciones representativo de los equipos existentes en un satélite. Realizado el plan de pruebas y aplicado sobre el hardware representativo se podrá dar por caracterizada la tarjeta a nivel SpW y realizar un informe con las conclusiones alcanzadas acerca del funcionamiento de los interfaces SpW de la tarjeta y las limitaciones encontradas.

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

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

HERMES: towards an integrated toolbox to characterize functional and effective brain connectivity

The analysis of the interdependence between time series has become an important field of research in the last years, mainly as a result of advances in the characterization of dynamical systems from the signals they produce, the introduction of concepts such as generalized and phase synchronization and the application of information theory to time series analysis. In neurophysiology, different analytical tools stemming from these concepts have added to the ?traditional? set of linear methods, which includes the cross-correlation and the coherency function in the time and frequency domain, respectively, or more elaborated tools such as Granger Causality. This increase in the number of approaches to tackle the existence of functional (FC) or effective connectivity (EC) between two (or among many) neural networks, along with the mathematical complexity of the corresponding time series analysis tools, makes it desirable to arrange them into a unified, easy-to-use software package. The goal is to allow neuroscientists, neurophysiologists and researchers from related fields to easily access and make use of these analysis methods from a single integrated toolbox. Here we present HERMES (http://hermes.ctb.upm.es), a toolbox for the Matlab® environment (The Mathworks, Inc), which is designed to study functional and effective brain connectivity from neurophysiological data such as multivariate EEG and/or MEG records. It includes also visualization tools and statistical methods to address the problem of multiple comparisons. We believe that this toolbox will be very helpful to all the researchers working in the emerging field of brain connectivity analysis.