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.

Mophological analysis of the human vasculature from medical imaging. Application to cerebral and coranary arteries

In this project, we have investigated new ways of modelling and analysis of human vasculature from Medical images. The research was divided in two main areas: cerebral vasculature analysis and coronary arteries modeling. Regarding cerebral vasculature analysis, we have studed cerebral aneurysms, internal carotid and the Circle of Willis (CoW). Aneurysms are abnormal vessel enlargements that can rupture causing important cerebral damages or death. The understanding of this pathology, together with its virtual treatment, and image diagnosis and prognosis, includes identification and detailed measurement of the aneurysms. In this context, we have proposed two automatic aneurysm isolation method, to separate the abnormal part of the vessel from the healthy part, to homogenize and speed-up the processing pipeline usually employed to study this pathology, [Cardenes2011TMI, arrabide2011MedPhys]. The results obtained from both methods have been also compared and validatied in [Cardenes2012MBEC]. A second important task here the analysis of the internal carotid [Bogunovic2011Media] and the automatic labelling of the CoW, Bogunovic2011MICCAI, Bogunovic2012TMI]. The second area of research covers the study of coronary arteries, specially coronary bifurcations because there is where the formation of atherosclerotic plaque is more common, and where the intervention is more challenging. Therefore, we proposed a novel modelling method from Computed Tomography Angiography (CTA) images, combined with Conventional Coronary Angiography (CCA), to obtain realistic vascular models of coronary bifurcations, presented in [Cardenes2011MICCAI], and fully validated including phantom experiments in [Cardene2013MedPhys]. The realistic models obtained from this method are being used to simulate stenting procedures, and to investigate the hemodynamic variables in coronary bifurcations in the works submitted in [Morlachi2012, Chiastra2012]. Additionally, another preliminary work has been done to reconstruct the coronary tree from rotational angiography, and published in [Cardenes2012ISBI].

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})$

The vitamin D receptor gene bAt (CCA) haplotype impairs the response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C patients.

BACKGROUND: Chronic hepatitis C infection is a major cause of end-stage liver disease. Therapy outcome is influenced by 25-OH vitamin D deficiency. To further address this observation, our study investigates the impact of the vitamin D receptor (NR1I1) haplotype and combined effects of plasma vitamin D levels in a well-described cohort of hepatitis C patients. METHODS: A total of 155 chronic hepatitis C patients were recruited from the Swiss Hepatitis C Cohort Study for NR1I1 genotyping and plasma 25-OH vitamin D level measurement. NR1I1 genotype data and combined effects of plasma 25-OH vitamin D level were analysed regarding therapy response (sustained virological response). RESULTS: A strong association was observed between therapy non-response and the NR1I1 CCA (bAt) haplotype consisting of rs1544410 (BsmI) C, rs7975232 (ApaI) C and rs731236 (TaqI) A alleles. Of the HCV patients carrying the CCA haplotype, 50.3% were non-responders (odds ratio [OR] 1.69, 95% CI 1.07, 2.67; P=0.028). A similar association was observed for the combinational CCCCAA genotype (OR 2.94, 95% CI 1.36, 6.37; P=0.007). The combinational CCCCAA genotype was confirmed as an independent risk factor for non-response in multivariate analysis (OR 2.50, 95% CI 1.07, 5.87; P=0.034). Analysing combined effects, a significant impact of low 25-OH vitamin D levels on sustained virological response were only seen in patients with the unfavourable NR1I1 CCA (bAt) haplotype (OR for non-SVR 3.55; 95% CI 1.005, 12.57; P=0.049). CONCLUSIONS: NR1I1 vitamin D receptor polymorphisms influence response to pegylated-interferon/ribavirin-based therapy in chronic hepatitis C and exert an additive genetic predisposition to previously described low 25-OH vitamin D serum levels.

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.

Synthetic calibrators for the analysis of total metanephrines in urine: revisiting the conditions of hydrolysis.

BACKGROUND: The quantification of total (free+sulfated) metanephrines in urine is recommended to diagnose pheochromocytoma. Urinary metanephrines include metanephrine itself, normetanephrine and methoxytyramine, mainly in the form of sulfate conjugates (60-80%). Their determination requires the hydrolysis of the sulfate ester moiety to allow electrochemical oxidation of the phenolic group. Commercially available urine calibrators and controls contain essentially free, unhydrolysable metanephrines which are not representative of native urines. The lack of appropriate calibrators may lead to uncertainty regarding the completion of the hydrolysis of sulfated metanephrines, resulting in incorrect quantification. METHODS: We used chemically synthesized sulfated metanephrines to establish whether the procedure most frequently recommended for commercial kits (pH 1.0 for 30 min over a boiling water bath) ensures their complete hydrolysis. RESULTS: We found that sulfated metanephrines differ in their optimum pH to obtain complete hydrolysis. Highest yields and minimal variance were established for incubation at pH 0.7-0.9 during 20 min. CONCLUSION: Urinary pH should be carefully controlled to ensure an efficient and reproducible hydrolysis of sulfated metanephrines. Synthetic sulfated metanephrines represent the optimal material for calibrators and proficiency testing to improve inter-laboratory accuracy.

Comparative analysis of the transcriptome across distant species.

The transcriptome is the readout of the genome. Identifying common features in it across distant species can reveal fundamental principles. To this end, the ENCODE and modENCODE consortia have generated large amounts of matched RNA-sequencing data for human, worm and fly. Uniform processing and comprehensive annotation of these data allow comparison across metazoan phyla, extending beyond earlier within-phylum transcriptome comparisons and revealing ancient, conserved features. Specifically, we discover co-expression modules shared across animals, many of which are enriched in developmental genes. Moreover, we use expression patterns to align the stages in worm and fly development and find a novel pairing between worm embryo and fly pupae, in addition to the embryo-to-embryo and larvae-to-larvae pairings. Furthermore, we find that the extent of non-canonical, non-coding transcription is similar in each organism, per base pair. Finally, we find in all three organisms that the gene-expression levels, both coding and non-coding, can be quantitatively predicted from chromatin features at the promoter using a 'universal model' based on a single set of organism-independent parameters.

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background: Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results: Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions: Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Multi-scale effects of nestling diet on breeding performance in a terrestrial top predator inferred from stable isotope analysis

Inter-individual diet variation within populations is likely to have important ecological and evolutionary implications. The diet-fitness relationships at the individual level and the emerging population processes are, however, poorly understood for most avian predators inhabiting complex terrestrial ecosystems. In this study, we use an isotopic approach to assess the trophic ecology of nestlings in a long-lived raptor, the Bonelli"s eagle Aquila fasciata, and investigate whether nestling dietary breath and main prey consumption can affect the species" reproductive performance at two spatial scales: territories within populations and populations over a large geographic area. At the territory level, those breeding pairs whose nestlings consumed similar diets to the overall population (i.e. moderate consumption of preferred prey, but complemented by alternative prey categories) or those disproportionally consuming preferred prey were more likely to fledge two chicks. An increase in the diet diversity, however, related negatively with productivity. The age and replacements of breeding pair members had also an influence on productivity, with more fledglings associated to adult pairs with few replacements, as expected in long-lived species. At the population level, mean productivity was higher in those population-years with lower dietary breadth and higher diet similarity among territories, which was related to an overall higher consumption of preferred prey. Thus, we revealed a correspondence in diet-fitness relationships at two spatial scales: territories and populations. We suggest that stable isotope analyses may be a powerful tool to monitor the diet of terrestrial avian predators on large spatio-temporal scales, which could serve to detect potential changes in the availability of those prey on which predators depend for breeding. We encourage ecologists and evolutionary and conservation biologists concerned with the multi-scale fitness consequences of inter-individual variation in resource use to employ similar stable isotope-based approaches, which can be successfully applied to complex ecosystems such as the Mediterranean.

Heterosis in papaya: inter and intragroup analysis

Papaya (Carica papaya L.) is a typical crop of tropical areas, and Brazil is one of the leading world producers. In recent decades, papaya culture has expanded to different regions of the country, but the number of cultivars available is still limited. In the present study, a complete diallel cross was carried out using eight accessions of papaya from the UENF/Caliman germplasm bank. Four genotypes belong to the Formosa heterotic group and four, to the Solo group. This study aimed to evaluate the occurrence and viability of exploring heterosis in heterotic intragroup hybrids. Fifty-six hybrid progenies were generated and evaluated. Among the Formosa intragroup hybrids, two hybrid combinations (MR x J4 and MR x SK) showed heterosis for all traits, as well as good average total fruit production. Among the Solo intragroup hybrids, three hybrid combinations (WM x GG, WM x SS and WM x SM) stand out for fruit production and high content of soluble solids. In Formosa x Solo hybrids, all hybrid combinations with the parent JS (JS x WM, JS x GG, JS x SS and JS x SM) showed high fruit quality and good average for fruit production. The heterotic profile of the hybrids tested allowed the identification of promising hybrids within Formosa and Solo heterotic groups. The analysis of the canonical variables also allowed the visualization of distinct groups of hybrids, depending on the provenance of the parents.

Virulence and rep-PCR analysis of Pyricularia grisea isolates from two Brazilian upland rice cultivars

The phenotypic and genetic diversity of 77 isolates of Pyricularia grisea collected from two upland rice cultivars, Maravilha and Primavera, was studied. Isolates exhibiting compatible reaction to cv.Primavera were incompatible to cv.Maravilha and vice versa, with the exception of six isolates that were compatible to both cultivars. The virulence of isolates from cv. Maravilha on 32 test genotypes of rice was significantly higher (t = 9.09, p < 0.0001) than the isolates from cv.Primavera. A phenogram constructed from virulence data showed two main groups, one constituted mainly of isolates from cv.Primavera (97.6%) and the other of isolates from cv.Maravilha (91.17%). Rep-PCR analysis of isolates using two primers designed from sequences of Pot2 showed that isolates could be clustered broadly into two groups. The average similarity within a cluster of isolates from cv.Primavera was significantly greater than the average similarity among the isolates of cv.Maravilha (t = 5.37, p < 0.0001). There was close correspondence between clusters based on PCR and virulence data (r = 0.48, p < 0.011). The results showed that isolates of P. grisea were cultivar specific and had low phenotypic and genetic diversity.

COMPARISON BETWEEN WOOD DRYING DEFECT SCORES: SPECIMEN TESTING X ANALYSIS OF KILN-DRIED BOARDS

It is important to develop drying technologies for Eucalyptus grandis lumber, which is one of the most planted species of this genus in Brazil and plays an important role as raw material for the wood industry. The general aim of this work was to assess the conventional kiln drying of juvenile wood of three clones of Eucalyptus grandis. The specific aims were to compare the behavior between: i) drying defects indicated by tests with wood specimens and conventional kiln-dried boards; and ii) physical properties and the drying quality. Five 11-year-old trees of each clone were felled, and only flatsawn boards of the first log were used. Basic density and total shrinkage were determined, and the drying test with wood specimens at 100 °C was carried out. Kiln drying of boards was performed, and initial and final moisture content, moisture gradient in thickness, drying stresses and drying defects were assessed. The defect scoring method was used to verify the behavior between the defects detected by specimen testing and the defects detected in kiln-dried boards. As main results, the drying schedule was too severe for the wood, resulting in a high level of boards with defects. The behavior between the defects in the drying test with specimens and the defects of kiln-dried boards was different, there was no correspondence, according to the defect scoring method.