Does Independent Component Analysis Play a Role in Unmixing Hyperspectral Data?

Chapter in Book Proceedings with Peer Review First Iberian Conference, IbPRIA 2003, Puerto de Andratx, Mallorca, Spain, JUne 4-6, 2003. Proceedings

Exploratory analysis of OpenStreetMap for land use classification

In the last years, volunteers have been contributing massively to what we know nowadays as Volunteered Geographic Information. This huge amount of data might be hiding a vast geographical richness and therefore research needs to be conducted to explore their potential and use it in the solution of real world problems. In this study we conduct an exploratory analysis of data from the OpenStreetMap initiative. Using the Corine Land Cover database as reference and continental Portugal as the study area, we establish a possible correspondence between both classification nomenclatures, evaluate the quality of OpenStreetMap polygon features classification against Corine Land Cover classes from level 1 nomenclature, and analyze the spatial distribution of OpenStreetMap classes over continental Portugal. A global classification accuracy around 76% and interesting coverage areas’ values are remarkable and promising results that encourages us for future research on this topic.

CLOCK-controlled polyphonic regulation of circadian rhythms through canonical and noncanonical E-boxes.

In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.

Dynamical modeling and analysis of large cellular regulatory networks.

The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.

Widespread occurrence of non-canonical transcription termination by human RNA polymerase III.

Human RNA polymerase (Pol) III-transcribed genes are thought to share a simple termination signal constituted by four or more consecutive thymidine residues in the coding DNA strand, just downstream of the RNA 3'-end sequence. We found that a large set of human tRNA genes (tDNAs) do not display any T(≥4) stretch within 50 bp of 3'-flanking region. In vitro analysis of tDNAs with a distanced T(≥4) revealed the existence of non-canonical terminators resembling degenerate T(≥5) elements, which ensure significant termination but at the same time allow for the production of Pol III read-through pre-tRNAs with unusually long 3' trailers. A panel of such non-canonical signals was found to direct transcription termination of unusual Pol III-synthesized viral pre-miRNA transcripts in gammaherpesvirus 68-infected cells. Genome-wide location analysis revealed that human Pol III tends to trespass into the 3'-flanking regions of tDNAs, as expected from extensive terminator read-through. The widespread occurrence of partial termination suggests that the Pol III primary transcriptome in mammals is unexpectedly enriched in 3'-trailer sequences with the potential to contribute novel functional ncRNAs.

Detection of matchings in a sequence of underwater images through texture analysis

This paper presents an approach to ameliorate the reliability of the correspondence points relating two consecutive images of a sequence. The images are especially difficult to handle, since they have been acquired by a camera looking at the sea floor while carried by an underwater robot. Underwater images are usually difficult to process due to light absorption, changing image radiance and lack of well-defined features. A new approach based on gray-level region matching and selective texture analysis significantly improves the matching reliability

Analysis of GPCR properties of Frizzled receptors and establishment of the "in vitro" platform for screening of Frizzled agonists/antagonists as potential therapeutic agents

Morphogens of the Wnt protein family are the secreted lipoglycoprotein ligands which initiate several pathways heavily involved in the coordination of various developmental stages of organisms in the majority of animal species. Deregulation of these pathways in the adult leads to formation and sustaining of multiple types of cancer. The latter notion is reinforced by the fact that the very discovery of the first Wnt ligand was due to its role as the causative factor of carcinogenic transformation (Nusse and Varmus, 1982). Nowadays our knowledge on Wnt signaling has "moved with the times" and these pathways were identified to be often crucial for tumor formation, its interactions with the microenvironment, and promotion of the metastases (Huang and Du, 2008; Zerlin et al., 2008; Jessen, 2009). Thus the relevance of the pathway as the target for drug development has further increased in the light of modern paradigms of the complex cancer treatments which target also spreading and growth- promoting factors of tumors by specific and highly efficient substances (Pavet et al., 2010). Presently the field of the Wnt-targeting drug research is almost solely dominated by assays based on transcriptional activation induced by the signaling. This approach resulted in development of a number of promising substances (Lee et al., 2011). Despite its effectiveness, the method nevertheless suffers from several drawbacks. Among the major ones is the fact that this approach is prone to identify compounds targeting rather downstream effectors of the pathway, which are indiscriminately used by all the subtypes of the Wnt signaling. Additionally, proteins which are involved in several signaling cascades and not just the Wnt pathway turn out as targets of the new compounds. These issues increase risks of side effects due to off-target interactions and blockade of the pathway in healthy cells. In the present work we put forward a novel biochemical approach for drug development on the Wnt pathway. It targets Frizzleds (Fzs) - a family of 7-transmbembrane proteins which serve as receptors for Wnt ligands. They offer unique properties for the development of highly specific and effective drugs as they control all branches of the Wnt signaling. Recent advances in the understanding of the roles of heterotrimeric G proteins downstream from Fzs (Katanaev et al., 2005; Liu et al., 2005; Jernigan et al., 2010) suggest application of enzymatic properties of these effectors to monitor the receptor-mediated events. We have applied this knowledge in practice and established a specific and efficient method based on utilization of a novel high-throughput format of the GTP-binding assay to follow the activation of Fzs. This type of assay is a robust and well-established technology for the research and screenings on the GPCRs (Harrison and Traynor, 2003). The conventional method of detection involves the radioactively labeled non-hydrolysable GTP analog [35S]GTPyS. Its application in the large-scale screenings is however problematic which promoted development of the novel non-radioactive GTP analog GTP-Eu. The new molecule employs phenomenon of the time-resolved fluorescence to provide sensitivity comparable to the conventional radioactive substance. Initially GTP-Eu was tested only in one of many possible types of GTP-binding assays (Frang et al., 2003). In the present work we expand these limits by demonstrating the general comparability of the novel label with the radioactive method in various types of assays. We provide a biochemical characterization of GTP-Eu interactions with heterotrimeric and small GTPases and a comparative analysis of the behavior of the new label in the assays involving heterotrimeric G protein effectors. These developments in the GTP-binding assay were then applied to monitor G protein activation by the Fz receptors. The data obtained in mammalian cultured cell lines provides for the first time an unambiguous biochemical proof for direct coupling of Fzs with G proteins. The specificity of this interaction has been confirmed by the experiments with the antagonists of Fz and by the pertussis toxin-mediated deactivation. Additionally we have identified the specificity of Wnt3a towards several members of the Fz family and analyzed the properties of human Fz-1 which was found to be the receptor coupled to the Gi/o family of G proteins. Another process playing significant role in the functioning of every GPCR is endocytosis. This phenomenon can also be employed for drug screenings on GPCRs (Bickle, 2010). In the present work we have demonstrated that Drosophila Fz receptors are involved in an unusual for many GPCRs manifestation of the receptor-mediated internalization. Through combination of biochemical approaches and studies on Drosophila as the model organism we have shown that direct interactions of the Fzs and the α-subunit of the heterotrimeric G protein Go with the small GTPase Rab5 regulate internalization of the receptor in early endosomes. We provide data uncovering the decisive role of this self-promoted endocytosis in formation of a proper signaling output in the canonical as well as planar cell polarity (PCP) pathways regulated by Fz. The results of this work thus establish a platform for the high-throughput screening to identify substances active in the cancer-related Wnt pathways. This methodology has been adjusted and applied to provide the important insights in Fz functioning and will be instrumental for further investigations on the Wnt-mediated pathways.

The generation and analysis of combined legless 1 and 2 knock-out mice

Summary The Wnt signaling pathway plays an important role during development and also for maintaining tissue homeostasis due to its function in proliferation, differentiation and cell fate decisions. Wnt ligands bind to Frizzled receptors and activate a signaling cascade that results in the stabilization of β-Catenin, a key component of the pathway. β-Catenin translocates to the nucleus, where, together with a transcription factor of the Tcf/Lef family, it activates the expression of target genes. Legless and Pygopus are two recently discovered essential components of the Wnt pathway in Drosophila, which may mediate the nuclear import and retention of beta-Catenin and/or contribute directly to the activation of Wnt target genes. To address the function of Legless in the mouse, we have generated compound constitutive and conditional knockout alleles of the two homologues legless 'I (bc1-9) and 2. We have induced the deletion of legless in self-renewing tissues such as the gastrointestinal tract, the mammary gland and the skin during adulthood and constitutively in the embryo. The present thesis focused on the consequences of the inactivation of legless in epithelial homeostasis as well as in a regeneration model and its comparison to pygopus. Deletion of neither legless nor pygopus in the adult small intestine resulted in any apparent anomaly, contrasting expectations from the phenotype caused by over-expression of Dickkopf, a Wnt inhibitor (Pinto et al., 2003). These observations indicate that canonical Wnt signaling might not be indispensable for normal gastrointestinal epithelium homeostasis, or that, in this context, Legless and Pygopus are not essential components of the Wnt pathway. However, the regeneration of the colonic epithelium after DSS induced damage was markedly impaired in legless, but not in pygopus deficient mice. Thus, unlike in Drosophila, deletion of mammalian legless and pygopus resulted in different phenotypes, suggesting that Legless might interact with as yet unidentified partners in addition to Pygopus. Resumé La voie de signalisation Wnt joue un rôle important au cours du développement ainsi que pour le maintien de l' homéostase tissulaire due à sa fonction durant la prolifération, la différentiation et les décisions sur l'avenir des cellules. Les ligands de Wnt se lient aux récepteurs Frizzled et activent une cascade de signalisation résultant en la stabilisation de β-Catenin, un composant central de cette voie. β-Catenin est transloquée dans le noyau ou, avec l'aide des facteurs de transcription de la famille Tcf/lef, elle active la transcription des gènes cibles. Legless et Pygopus sont deux composants récemment découverts et essentiels de la voie de signalisation Wnt chez la Drosophile qui pourraient être des médiateurs de l'import et de la rétention nucléaire de bêta-catenin et/ou contribuer directement a l'activation des gènes cibles. Afin de comprendre la fonction de Legless chez la souris, nous avons généré simultanément les allèles « knock-out » constitutifs et conditionnels des deux homologues legless 1 (bc1-9) et 2. Nous avons induit la délétion de legless dans des tissus capables de s'auto renouveler comme le tract gastro-intestinal, la glande mammaire et la peau chez l'adulte et nous avons supprimé constitutivement legless chez l'embryon. La présente thèse est concentrée sur les conséquences de l'inactivation de legless au cours de l' homéostase épithéliale ainsi que dans un modèle de régénération et sur sa comparaison avec pygopus. Ni la délétion de legless ni celle de pygopus dans l'intestin adulte n'ont résulté en quelque anomalie, contrastant nos attentes provenant des phénotypes causes par la surexpression de Dickkpof, un inhibiteur de Wnt (Pinto et al., 2003). Ces observations indiquent que la voie de signalisation Wnt/β-Catenin pourrait ne pas être indispensable à l' homéostase normale du tract gastro-intestinal, ou que, dans ce contexte, Legless et Pygopus ne sont pas des composants essentiels de la vole Wnt. Cependant, la régénération de l'épithélium du colon après induction de son endommagement au DSS fut dramatiquement diminuée chez legless mais pas chez les souris mutantes pour pygopus. Ainsi, a la différence de chez la Drosophile, la délétion de legless et pygopus chez les mammifères a résulté en des phénotypes différents, suggérant que Legless pourrait interagir avec d'autres partenaires, encore non identifies, que Pygopus.

Application of discriminant function analysis and change-point problem in dating volcanic ashes

The application of Discriminant function analysis (DFA) is not a new idea in the studyof tephrochrology. In this paper, DFA is applied to compositional datasets of twodifferent types of tephras from Mountain Ruapehu in New Zealand and MountainRainier in USA. The canonical variables from the analysis are further investigated witha statistical methodology of change-point problems in order to gain a betterunderstanding of the change in compositional pattern over time. Finally, a special caseof segmented regression has been proposed to model both the time of change and thechange in pattern. This model can be used to estimate the age for the unknown tephrasusing Bayesian statistical calibration

Multisensor change detection with nonlinear canonical correlation

The analysis of multi-modal and multi-sensor images is nowadays of paramount importance for Earth Observation (EO) applications. There exist a variety of methods that aim at fusing the different sources of information to obtain a compact representation of such datasets. However, for change detection existing methods are often unable to deal with heterogeneous image sources and very few consider possible nonlinearities in the data. Additionally, the availability of labeled information is very limited in change detection applications. For these reasons, we present the use of a semi-supervised kernel-based feature extraction technique. It incorporates a manifold regularization accounting for the geometric distribution and jointly addressing the small sample problem. An exhaustive example using Landsat 5 data illustrates the potential of the method for multi-sensor change detection.

TCRep 3D: an automated in silico approach to study the structural properties of TCR repertoires.

TCRep 3D is an automated systematic approach for TCR-peptide-MHC class I structure prediction, based on homology and ab initio modeling. It has been considerably generalized from former studies to be applicable to large repertoires of TCR. First, the location of the complementary determining regions of the target sequences are automatically identified by a sequence alignment strategy against a database of TCR Vα and Vβ chains. A structure-based alignment ensures automated identification of CDR3 loops. The CDR are then modeled in the environment of the complex, in an ab initio approach based on a simulated annealing protocol. During this step, dihedral restraints are applied to drive the CDR1 and CDR2 loops towards their canonical conformations, described by Al-Lazikani et. al. We developed a new automated algorithm that determines additional restraints to iteratively converge towards TCR conformations making frequent hydrogen bonds with the pMHC. We demonstrated that our approach outperforms popular scoring methods (Anolea, Dope and Modeller) in predicting relevant CDR conformations. Finally, this modeling approach has been successfully applied to experimentally determined sequences of TCR that recognize the NY-ESO-1 cancer testis antigen. This analysis revealed a mechanism of selection of TCR through the presence of a single conserved amino acid in all CDR3β sequences. The important structural modifications predicted in silico and the associated dramatic loss of experimental binding affinity upon mutation of this amino acid show the good correspondence between the predicted structures and their biological activities. To our knowledge, this is the first systematic approach that was developed for large TCR repertoire structural modeling.

Some spectral properties of the canonical solution operator to $\bar\partial$ on weighted Fock spaces

We characterize the Schatten class membership of the canonical solution operator to $\overline{\partial}$ acting on $L^2(e^{-2\phi})$, where $\phi$ is a subharmonic function with $\Delta\phi$ a doubling measure. The obtained characterization is in terms of $\Delta\phi$. As part of our approach, we study Hankel operators with anti-analytic symbols acting on the corresponding Fock space of entire functions in $L^2(e^{-2\phi})$

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of p H and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups.

A Canonical Genetic Algorithm for Blind Inversion of Linear Channels

It is well known the relationship between source separation and blind deconvolution: If a filtered version of an unknown i.i.d. signal is observed, temporal independence between samples can be used to retrieve the original signal, in the same manner as spatial independence is used for source separation. In this paper we propose the use of a Genetic Algorithm (GA) to blindly invert linear channels. The use of GA is justified in the case of small number of samples, where other gradient-like methods fails because of poor estimation of statistics.

Genetic diversity between improved banana diploids using canonical variables and the Ward-MLM method

The objective of this work was to estimate the genetic diversity of improved banana diploids using data from quantitative analysis and from simple sequence repeats (SSR) marker, simultaneously. The experiment was carried out with 33 diploids, in an augmented block design with 30 regular treatments and three common ones. Eighteen agronomic characteristics and 20 SSR primers were used. The agronomic characteristics and the SSR were analyzed simultaneously by the Ward-MLM, cluster, and IML procedures. The Ward clustering method considered the combined matrix obtained by the Gower algorithm. The Ward-MLM procedure identified three ideal groups (G1, G2, and G3) based on pseudo-F and pseudo-t² statistics. The dendrogram showed relative similarity between the G1 genotypes, justified by genealogy. In G2, 'Calcutta 4' appears in 62% of the genealogies. Similar behavior was observed in G3, in which the 028003-01 diploid is the male parent of the 086079-10 and 042079-06 genotypes. The method with canonical variables had greater discriminatory power than Ward-MLM. Although reduced, the genetic variability available is sufficient to be used in the development of new hybrids.