Machine learning algorithms for spatial data. Case studies: Environmental pollution, natural hazards, renewable resources

This paper presents general problems and approaches for the spatial data analysis using machine learning algorithms. Machine learning is a very powerful approach to adaptive data analysis, modelling and visualisation. The key feature of the machine learning algorithms is that they learn from empirical data and can be used in cases when the modelled environmental phenomena are hidden, nonlinear, noisy and highly variable in space and in time. Most of the machines learning algorithms are universal and adaptive modelling tools developed to solve basic problems of learning from data: classification/pattern recognition, regression/mapping and probability density modelling. In the present report some of the widely used machine learning algorithms, namely artificial neural networks (ANN) of different architectures and Support Vector Machines (SVM), are adapted to the problems of the analysis and modelling of geo-spatial data. Machine learning algorithms have an important advantage over traditional models of spatial statistics when problems are considered in a high dimensional geo-feature spaces, when the dimension of space exceeds 5. Such features are usually generated, for example, from digital elevation models, remote sensing images, etc. An important extension of models concerns considering of real space constrains like geomorphology, networks, and other natural structures. Recent developments in semi-supervised learning can improve modelling of environmental phenomena taking into account on geo-manifolds. An important part of the study deals with the analysis of relevant variables and models' inputs. This problem is approached by using different feature selection/feature extraction nonlinear tools. To demonstrate the application of machine learning algorithms several interesting case studies are considered: digital soil mapping using SVM, automatic mapping of soil and water system pollution using ANN; natural hazards risk analysis (avalanches, landslides), assessments of renewable resources (wind fields) with SVM and ANN models, etc. The dimensionality of spaces considered varies from 2 to more than 30. Figures 1, 2, 3 demonstrate some results of the studies and their outputs. Finally, the results of environmental mapping are discussed and compared with traditional models of geostatistics.

Adaptive strategy for the statistical analysis of connectomes.

We study an adaptive statistical approach to analyze brain networks represented by brain connection matrices of interregional connectivity (connectomes). Our approach is at a middle level between a global analysis and single connections analysis by considering subnetworks of the global brain network. These subnetworks represent either the inter-connectivity between two brain anatomical regions or by the intra-connectivity within the same brain anatomical region. An appropriate summary statistic, that characterizes a meaningful feature of the subnetwork, is evaluated. Based on this summary statistic, a statistical test is performed to derive the corresponding p-value. The reformulation of the problem in this way reduces the number of statistical tests in an orderly fashion based on our understanding of the problem. Considering the global testing problem, the p-values are corrected to control the rate of false discoveries. Finally, the procedure is followed by a local investigation within the significant subnetworks. We contrast this strategy with the one based on the individual measures in terms of power. We show that this strategy has a great potential, in particular in cases where the subnetworks are well defined and the summary statistics are properly chosen. As an application example, we compare structural brain connection matrices of two groups of subjects with a 22q11.2 deletion syndrome, distinguished by their IQ scores.

Automatic front-crawl temporal phase detection using adaptive filtering of inertial signals

This study introduces a novel approach for automatic temporal phase detection and inter-arm coordination estimation in front-crawl swimming using inertial measurement units (IMUs). We examined the validity of our method by comparison against a video-based system. Three waterproofed IMUs (composed of 3D accelerometer, 3D gyroscope) were placed on both forearms and the sacrum of the swimmer. We used two underwater video cameras in side and frontal views as our reference system. Two independent operators performed the video analysis. To test our methodology, seven well-trained swimmers performed three 300 m trials in a 50 m indoor pool. Each trial was in a different coordination mode quantified by the index of coordination. We detected different phases of the arm stroke by employing orientation estimation techniques and a new adaptive change detection algorithm on inertial signals. The difference of 0.2 +/- 3.9% between our estimation and video-based system in assessment of the index of coordination was comparable to experienced operators' difference (1.1 +/- 3.6%). The 95% limits of agreement of the difference between the two systems in estimation of the temporal phases were always less than 7.9% of the cycle duration. The inertial system offers an automatic easy-to-use system with timely feedback for the study of swimming.

Eumelanin-based coloration and fitness parameters in birds: a meta-analysis

Although melanin is the most common pigment in animal integuments, the adaptive function of variation in melanin-based coloration remains poorly understood. The individual fitness returns associated with melanin pigments can be variable across species as these pigments can have physical and biological protective properties and genes involved in melanogenesis may vary in the intensity of pleiotropic effects. Moreover, dark and pale coloration can also enhance camouflage in alternative habitats and melanin-based coloration can be involved in social interactions. We investigated whether darker or paler individuals achieve a higher fitness in birds, a taxon wherein associations between melanin-based coloration and fitness parameters have been studied in a large number of species. A meta-analysis showed that the degree of melanin-based coloration was not significantly associated with laying date, clutch size, brood size, and survival across 26 species. Similar results were found when restricting the analyses to non-sexually dimorphic birds, colour polymorphic and monomorphic species, in passerines and non-passerines and in species for which inter-individual variation in melanism is due to colour intensity. However, eumelanic coloration was positively associated with clutch and brood size in sexually dimorphic species and those that vary in the size of black patches, respectively. Given that greater extent of melanin-based coloration was positively associated with reproductive parameters and survival in some species but negatively in other species, we conclude that in birds the sign and magnitude of selection exerted on melanin-based coloration is species- or trait-specific.

Control of the adaptive response of the heart to stress via the Notch1 receptor pathway.

In the damaged heart, cardiac adaptation relies primarily on cardiomyocyte hypertrophy. The recent discovery of cardiac stem cells in the postnatal heart, however, suggests that these cells could participate in the response to stress via their capacity to regenerate cardiac tissues. Using models of cardiac hypertrophy and failure, we demonstrate that components of the Notch pathway are up-regulated in the hypertrophic heart. The Notch pathway is an evolutionarily conserved cell-to-cell communication system, which is crucial in many developmental processes. Notch also plays key roles in the regenerative capacity of self-renewing organs. In the heart, Notch1 signaling takes place in cardiomyocytes and in mesenchymal cardiac precursors and is activated secondary to stimulated Jagged1 expression on the surface of cardiomyocytes. Using mice lacking Notch1 expression specifically in the heart, we show that the Notch1 pathway controls pathophysiological cardiac remodeling. In the absence of Notch1, cardiac hypertrophy is exacerbated, fibrosis develops, function is altered, and the mortality rate increases. Therefore, in cardiomyocytes, Notch controls maturation, limits the extent of the hypertrophic response, and may thereby contribute to cell survival. In cardiac precursors, Notch prevents cardiogenic differentiation, favors proliferation, and may facilitate the expansion of a transient amplifying cell compartment.

Split sex ratios in the social Hymenoptera: a meta-analysis

The study of sex allocation in social Hymenoptera (ants, bees, and wasps) provides an excellent opportunity for testing kin-selection theory and studying conflict resolution. A queen-worker conflict over sex allocation is expected because workers are more related to sisters than to brothers, whereas queens are equally related to daughters and sons. If workers fully control sex allocation, split sex ratio theory predicts that colonies with relatively high or low relatedness asymmetry (the relatedness of workers to females divided by the relatedness of workers to males) should specialize in females or males, respectively. We performed a meta-analysis to assess the magnitude of adaptive sex allocation biasing by workers and degree of support for split sex ratio theory in the social Hymenoptera. Overall, variation in relatedness asymmetry (due to mate number or queen replacement) and variation in queen number (which also affects relatedness asymmetry in some conditions) explained 20.9% and 5% of the variance in sex allocation among colonies, respectively. These results show that workers often bias colony sex allocation in their favor as predicted by split sex ratio theory, even if their control is incomplete and a large part of the variation among colonies has other causes. The explanatory power of split sex ratio theory was close to that of local mate competition and local resource competition in the few species of social Hymenoptera where these factors apply. Hence, three of the most successful theories explaining quantitative variation in sex allocation are based on kin selection.

Local-global splitting for spatiotemporal-adaptive multiscale methods

We present a novel spatiotemporal-adaptive Multiscale Finite Volume (MsFV) method, which is based on the natural idea that the global coarse-scale problem has longer characteristic time than the local fine-scale problems. As a consequence, the global problem can be solved with larger time steps than the local problems. In contrast to the pressure-transport splitting usually employed in the standard MsFV approach, we propose to start directly with a local-global splitting that allows to locally retain the original degree of coupling. This is crucial for highly non-linear systems or in the presence of physical instabilities. To obtain an accurate and efficient algorithm, we devise new adaptive criteria for global update that are based on changes of coarse-scale quantities rather than on fine-scale quantities, as it is routinely done before in the adaptive MsFV method. By means of a complexity analysis we show that the adaptive approach gives a noticeable speed-up with respect to the standard MsFV algorithm. In particular, it is efficient in case of large upscaling factors, which is important for multiphysics problems. Based on the observation that local time stepping acts as a smoother, we devise a self-correcting algorithm which incorporates the information from previous times to improve the quality of the multiscale approximation. We present results of multiphase flow simulations both for Darcy-scale and multiphysics (hybrid) problems, in which a local pore-scale description is combined with a global Darcy-like description. The novel spatiotemporal-adaptive multiscale method based on the local-global splitting is not limited to porous media flow problems, but it can be extended to any system described by a set of conservation equations.

Adaptive filtering methods for identifying cross-frequency couplings in human EEG.

Oscillations have been increasingly recognized as a core property of neural responses that contribute to spontaneous, induced, and evoked activities within and between individual neurons and neural ensembles. They are considered as a prominent mechanism for information processing within and communication between brain areas. More recently, it has been proposed that interactions between periodic components at different frequencies, known as cross-frequency couplings, may support the integration of neuronal oscillations at different temporal and spatial scales. The present study details methods based on an adaptive frequency tracking approach that improve the quantification and statistical analysis of oscillatory components and cross-frequency couplings. This approach allows for time-varying instantaneous frequency, which is particularly important when measuring phase interactions between components. We compared this adaptive approach to traditional band-pass filters in their measurement of phase-amplitude and phase-phase cross-frequency couplings. Evaluations were performed with synthetic signals and EEG data recorded from healthy humans performing an illusory contour discrimination task. First, the synthetic signals in conjunction with Monte Carlo simulations highlighted two desirable features of the proposed algorithm vs. classical filter-bank approaches: resilience to broad-band noise and oscillatory interference. Second, the analyses with real EEG signals revealed statistically more robust effects (i.e. improved sensitivity) when using an adaptive frequency tracking framework, particularly when identifying phase-amplitude couplings. This was further confirmed after generating surrogate signals from the real EEG data. Adaptive frequency tracking appears to improve the measurements of cross-frequency couplings through precise extraction of neuronal oscillations.

Detection of ingested cocaine-filled packets: comparison of filtered back projection CT with adaptive statistical iterative reconstructed images

Purpose: To evaluate the diagnostic value and image quality of CT with filtered back projection (FBP) compared with adaptive statistical iterative reconstructed images (ASIR) in body stuffers with ingested cocaine-filled packets.Methods and Materials: Twenty-nine body stuffers (mean age 31.9 years, 3 women) suspected for ingestion of cocaine-filled packets underwent routine-dose 64-row multidetector CT with FBP (120kV, pitch 1.375, 100-300 mA and automatic tube current modulation (auto mA), rotation time 0.7sec, collimation 2.5mm), secondarily reconstructed with 30 % and 60 % ASIR. In 13 (44.83%) out of the body stuffers cocaine-filled packets were detected, confirmed by exact analysis of the faecal content including verification of the number (range 1-25). Three radiologists independently and blindly evaluated anonymous CT examinations (29 FBP-CT and 68 ASIR-CT) for the presence and number of cocaine-filled packets indicating observers' confidence, and graded them for diagnostic quality, image noise, and sharpness. Sensitivity, specificity, area under the receiver operating curve (ROC) Az and interobserver agreement between the 3 radiologists for FBP-CT and ASIR-CT were calculated.Results: The increase of the percentage of ASIR significantly diminished the objective image noise (p<0.001). Overall sensitivity and specificity for the detection of the cocaine-filled packets were 87.72% and 76.15%, respectively. The difference of ROC area Az between the different reconstruction techniques was significant (p= 0.0101), that is 0.938 for FBP-CT, 0.916 for 30 % ASIR-CT, and 0.894 for 60 % ASIR-CT.Conclusion: Despite the evident image noise reduction obtained by ASIR, the diagnostic value for detecting cocaine-filled packets decreases, depending on the applied ASIR percentage.

Monitoring of innate and adaptive immune responses in the context of human immunotherapy

Summary1 SummaryCancer patients have a better clinical outcome when their tumours display marked infiltration by memory Τ cells. Moreover, the overrepresentation of Th1 gene signatures in primary tumours correlates with favourable prognosis. Thus, vaccination to induce Τ cells capable of infiltrating and eradicating the tumour seems a promising strategy for the treatment of cancer. Here, I monitored CD4 Τ cell responses in melanoma patients vaccinated with the long synthetic peptides Melan- A16-35(A27L) and NY-ESO-179.108. Most of the patients developed strong and diverse peptide antigen specific CD4 Τ cell responses. Analysis of the fine specificity of CD4 Τ cell antigen recognition led to the identification of two new epitopes. The peptide Melan-A16_35(A27L) was delivered by virus-like particles (VLPs) derived from bacteriophage Οβ, which themselves displayed strong immunogenicity. I show evidence for induction of Οβ- and Melan-A specific CD4 Τ cell responses that developed a Th1 functional profile after repeated vaccination cycles. They also specifically released the chemokines CCL-3 and CCL-4, which play important roles in attracting CD8 Τ cells to the APC surface for priming and formation of Τ cell memory. We further found induction of robust humoral IgG responses upon VLP vaccination, and the lgG1-lgG4 isotype composition depended on the adjuvant used. Since heavy chain class switching largely dépends on the presence of CD4 Τ cell help, this result suggests that the adjuvant can influence the differentiation of elicited CD4 Τ cells, thereby contributing to the quality and function of both Β cells and CD8 Τ cells. The nature of the inflammatory processes in the tumour microenvironment can modulate CD8 Τ cell function. A collaboration was established for the investigation regulation of inflammasome activation in human primary monocytes. We identified IL- 4 and TGF-β as strong inhibitors of IL-1 β secretion, Indicating some level of regulation from effector Th2 and Treg responses. We further found a potent inhibition of inflammasome activation by type I interferon, and demonstrated in vivo inhibition of IL-1 β responses in monocytes from active multiple sclerosis patients under IFN-β therapy. This finding further offers a possible explanation for its success, which mechanism of action is still largely unclear. Interestingly, type I interferon is also being used as adjuvant treatment for tumour free metastatic cutaneous melanoma patients. While its clinical benefit has remained controversial, recent data suggest that the subset of patients with ulcerated primary melanoma lesions can benefit from this therapy. Future investigations will shed light on the implication of the inflammasome in this context, and may offer new strategies for improved adjuvant treatments of melanoma.2 RésuméLes patients atteints de cancer ont une meilleure chance de survie si leurs tumeurs s'avèrent être largement infiltrées par des cellules Τ mémoires. De plus, la surreprésentation d'une signature génique Th1 est en corrélation avec un pronostic favorable. Ainsi, la vaccination visant à induire des cellules Τ capables d'infiltrer et de détruire la tumeur parait être une stratégie prometteuse pour le traitement du cancer. Dans ce travail, j'ai procédé au monitoring de la réponse des cellules Τ CD4 dans des patients atteints de mélanome vaccinés avec les longs peptides synthétiques Melan-A16_35(A27L) et NY-ESO-179_108. Ces peptides représentent des antigènes tumoraux reconnus par des lymphocytes T. La majorité des patients a développé une réponse forte et diversifiée des cellules Τ CD4 spécifiques contre les peptides. L'analyse de la spécificité fine de la reconnaissance antigénique des cellules Τ CD4 nous a conduits à l'identification de deux nouveaux épitopes. Le peptide Melan-Aie. 35(A27L) a été délivré par des particules de type viral (VLPs) dérivés de bactériophages Qβ, qui ont eux-mêmes démontré une forte immunogénicité. Mon travail montre les preuves d'une induction de réponses spécifiques des cellules Τ CD4 contre les Qβ et Melan-A développant un profil fonctionnel Th1 après plusieurs cycles de vaccination. Elles secrètent aussi spécifiquement les chimiokines CCL-3 et CCL-4, qui jouent un rôle important dans l'attraction des cellules Τ CD8 à la surface des cellules présentatrices d'antigènes et contribuent ainsi à induire et former la mémoire cellulaire Τ CD8. Nous avons également remarqué une induction de fortes réponses humorales IgG après vaccination avec les VLPs, et que la composition des isotypes lgG1-lgG4 dépendait de l'adjuvant utilisé. Etant donné qu'une commutation de classe de la chaîne lourde dépend largement ùie l'aide des cellules Τ CD4, ce résultat suggère que l'adjuvant puisse influencer la différeritiation de cellules Τ CD4 en différent types, contribuant ainsi à la qualité et à la fonction des cellules Β et des cellules Τ CD8.La nature des processus d'inflammation dans le microenvironnement tumoral peut moduler la fonction des cellules Τ CD8. Une collaboration a été établie pour investiguer la régulation de l'activation de l'inflammasome dans des monocytes primaires humains. Nous avons identifié l'IL-4 et le TGF-β comme étant de puissants inhibiteurs de la sécrétion de IL-Ιβ, indiquant une certaine régulation de la réponse inflammatoire induite par les cellules Th2 et Τ régulatrices. Nous avons également trouvé une forte inhibition de l'activation de l'inflammasome par l'interféron type I, et nous avons démontré une inhibition in vivo de la réponse IL-1 β dans des monocytes de patients atteints d'une sclérose en plaque active sous traitement IFN-β. Ce résultat nous offre une possible explication du succès de cette thérapie, dont le mécanisme reste à ce jour encore largement obscur. Il est intéressant de noter que l'interféron de type I est également utilisé pour le traitement de patients atteints de mélanome cutané métastasique sans tumeurs. Bien que le bénéfice clinique de ce traitement reste controversé, des études récentes montrent qu'une partie des patients atteints de mélanome primaire ulcéré peut tirer bénéfice de cette thérapie. De futures investigations pourront mieux nous renseigner sur l'implication de l'inflammasome dans ce contexte et offrir de nouvelles stratégies pour améliorer les traitements adjuvants du mélanome.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Adaptive sequence evolution in a color gene involved in the formation of the characteristic egg-dummies of male haplochromine cichlid fishes.

BACKGROUND: The exceptionally diverse species flocks of cichlid fishes in East Africa are prime examples of parallel adaptive radiations. About 80% of East Africa's more than 1 800 endemic cichlid species, and all species of the flocks of Lakes Victoria and Malawi, belong to a particularly rapidly evolving lineage, the haplochromines. One characteristic feature of the haplochromines is their possession of egg-dummies on the males' anal fins. These egg-spots mimic real eggs and play an important role in the mating system of these maternal mouthbrooding fish. RESULTS: Here, we show that the egg-spots of haplochromines are made up of yellow pigment cells, xanthophores, and that a gene coding for a type III receptor tyrosine kinase, colony-stimulating factor 1 receptor a (csf1ra), is expressed in egg-spot tissue. Molecular evolutionary analyses reveal that the extracellular ligand-binding and receptor-interacting domain of csf1ra underwent adaptive sequence evolution in the ancestral lineage of the haplochromines, coinciding with the emergence of egg-dummies. We also find that csf1ra is expressed in the egg-dummies of a distantly related cichlid species, the ectodine cichlid Ophthalmotilapia ventralis, in which markings with similar functions evolved on the pelvic fin in convergence to those of the haplochromines. CONCLUSION: We conclude that modifications of existing signal transduction mechanisms might have evolved in the haplochromine lineage in association with the origination of anal fin egg-dummies. That positive selection has acted during the evolution of a color gene that seems to be involved in the morphogenesis of a sexually selected trait, the egg-dummies, highlights the importance of further investigations of the comparative genomic basis of the phenotypic diversification of cichlid fishes.

Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation.

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.

The adaptive dynamics of niche constructing traits in spatially subdivided populations: evolving posthumous extended phenotypes.

Niche construction, by which organisms modify the environment in which they live, has been proposed to affect the evolution of many phenotypic traits. But what about the evolution of a niche constructing trait itself, whose expression changes the pattern of natural selection to which the trait is exposed in subsequent generations? This article provides an inclusive fitness analysis of selection on niche constructing phenotypes, which can affect their environment from local to global scales in arbitrarily spatially subdivided populations. The model shows that phenotypic effects of genes extending far beyond the life span of the actor can be affected by natural selection, provided they modify the fitness of those individuals living in the future that are likely to have inherited the niche construction lineage of the actor. Present benefits of behaviors are thus traded off against future indirect costs. The future costs will generally result from a complicated interplay of phenotypic effects, population demography and environmental dynamics. To illustrate these points, I derive the adaptive dynamics of a trait involved in the consumption of an abiotic resource, where resource abundance in future generations feeds back to the evolutionary dynamics of the trait.

A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion

An epidemic model is formulated by a reactionâeuro"diffusion system where the spatial pattern formation is driven by cross-diffusion. The reaction terms describe the local dynamics of susceptible and infected species, whereas the diffusion terms account for the spatial distribution dynamics. For both self-diffusion and cross-diffusion, nonlinear constitutive assumptions are suggested. To simulate the pattern formation two finite volume formulations are proposed, which employ a conservative and a non-conservative discretization, respectively. An efficient simulation is obtained by a fully adaptive multiresolution strategy. Numerical examples illustrate the impact of the cross-diffusion on the pattern formation.