2

The quality of environmental data analysis and propagation of errors are heavily affected by the representativity of the initial sampling design [CRE 93, DEU 97, KAN 04a, LEN 06, MUL07]. Geostatistical methods such as kriging are related to field samples, whose spatial distribution is crucial for the correct detection of the phenomena. Literature about the design of environmental monitoring networks (MN) is widespread and several interesting books have recently been published [GRU 06, LEN 06, MUL 07] in order to clarify the basic principles of spatial sampling design (monitoring networks optimization) based on Support Vector Machines was proposed. Nonetheless, modelers often receive real data coming from environmental monitoring networks that suffer from problems of non-homogenity (clustering). Clustering can be related to the preferential sampling or to the impossibility of reaching certain regions.

Insights into the biogeographical history of the Lower Guinea Forest Domain: evidence for the role of refugia in the intraspecific differentiation of Aucoumea klaineana.

Determining the biogeographical histories of rainforests is central to our understanding of the present distribution of tropical biodiversity. Ice age fragmentation of central African rainforests strongly influenced species distributions. Elevated areas characterized by higher species richness and endemism have been postulated to be Pleistocene forest refugia. However, it is often difficult to separate the effects of history and of present-day ecological conditions on diversity patterns at the interspecific level. Intraspecific genetic variation could yield new insights into history, because refugia hypotheses predict patterns not expected on the basis of contemporary environmental dynamics. Here, we test geographically explicit hypotheses of vicariance associated with the presence of putative refugia and provide clues about their location. We intensively sampled populations of Aucoumea klaineana, a forest tree sensitive to forest fragmentation, throughout its geographical range. Characterizing variation at 10 nuclear microsatellite loci, we were able to obtain phylogeographic data of unprecedented detail for this region. Using Bayesian clustering approaches, we demonstrated the presence of four differentiated genetic units. Their distribution matched that of forest refugia postulated from patterns of species richness and endemism. Our data also show differences in diversity dynamics at leading and trailing edges of the species' shifting distribution. Our results confirm predictions based on refugia hypotheses and cannot be explained on the basis of present-day ecological conditions.

Methicillin-resistant Staphylococcus aureus: phylogenetic relatedness between European epidemic clones and Swiss sporadic strains.

We have compared the phylogenetic diversity of methicillin-resistant Staphylococcus aureus (MRSA) strains from Switzerland and their phylogenetic relationships with European epidemic clones, using multiprimer random amplification polymorphic DNA (RAPD). Strains included 24 European epidemic clones (59 strains), 66 sporadic strains isolated in Switzerland in 1996-1997, and 15 reference strains of five other Staphylococcus species. Similarity and clustering analysis with the Jaccard's coefficient showed that the maximum genetic distance between MRSA strains was 0.43, whereas the minimum genetic distance between the six Staphylococcus species was 0.97, indicating that the method permits phylogenetic hierarchization. The 24 MRSA clones reported to be epidemic in European countries during the 1990s were distributed into seven different genetic clusters with a maximum distance of 0.29 among them. This clustering pattern was confirmed by the analysis of a subset of MRSA strains by multilocus enzyme electrophoresis at 12 loci. Most of the sporadic Swiss strains were distributed into these seven different genetic clusters, together with the epidemic MRSA clones. This suggests that there is no phylogenetic cluster specific to epidemic clones of MRSA.

Cluster-based active learning for compact image classification

In this paper, we consider active sampling to label pixels grouped with hierarchical clustering. The objective of the method is to match the data relationships discovered by the clustering algorithm with the user's desired class semantics. The first is represented as a complete tree to be pruned and the second is iteratively provided by the user. The active learning algorithm proposed searches the pruning of the tree that best matches the labels of the sampled points. By choosing the part of the tree to sample from according to current pruning's uncertainty, sampling is focused on most uncertain clusters. This way, large clusters for which the class membership is already fixed are no longer queried and sampling is focused on division of clusters showing mixed labels. The model is tested on a VHR image in a multiclass classification setting. The method clearly outperforms random sampling in a transductive setting, but cannot generalize to unseen data, since it aims at optimizing the classification of a given cluster structure.

Aging preferentially affects molecular pathways implicated in development and disease of myelinating glial cells

The central and peripheral nervous systems are involved in multiple age-dependent neurological deficits that are often attributed to alterations in function of myelinating glial cells. However, the molecular events that underlie the age-related decline of glial cell function are unknown. We used Schwann cells as a model to study biological processes affected in glial cells by aging. We comprehensively profiled gene expression of the Schwann cellrich mouse sciatic nerve throughout life, from day of birth until senescence (840 days of age). We combined the aging data with the microarray transcriptional data obtained using nerves isolated from Schwann cell-specific neuropathy-inducing mutants MPZCre/+/Lpin1fE2−3/fE2−3 , MPZCre/+/ScapfE1/fE1 and Pmp22-null mice. The majority of age related transcripts were also affected in the analyzed mouse models of neuropathy (54.4%) and in development (59.5%) indicating a high level of overlapping in implicated molecular pathways. We observed that compared to peripheral nerve development, dynamically changing expression profiles in aging have opposite (anticorrelated) orientation while they copy the orientation of transcriptional changes observed in analyzed neuropathy models. Subsequent clustering and biological annotation of dynamically changing transcripts revealed that the processes most significantly deregulated in aging include inflammatory/immune response and lipid biosynthesis/metabolism. Importantly, the changes in these pathways were also observed in myelinated oligodendrocyte-rich optic nerves of aged mice, albeit with lower magnitude. This observation suggests that similar biological processes are affected in aging glial cells in central and peripheral nervous systems, however with different dynamics. Our data, which provide the first comprehensive comparison of molecular changes in glial cells in three distinct biological conditions comprising development, aging and disease, provide not only a new inside into the molecular alterations underlying neural system aging but also identify target pathways for potential therapeutic approaches to prevent or delay complications associated with age-related and inherited forms of neuropathies. *Current address: Department of Physiology, UCSF, San Francisco, CA, USA.

A frequency distribution approach to hotspot identification

We present a new global method for the identification of hotspots in conservation and ecology. The method is based on the identification of spatial structure properties through cumulative relative frequency distributions curves, and is tested with two case studies, the identification of fish density hotspots and terrestrial vertebrate species diversity hotspots. Results from the frequency distribution method are compared with those from standard techniques among local, partially local and global methods. Our approach offers the main advantage to be independent from the selection of any threshold, neighborhood, or other parameter that affect most of the currently available methods for hotspot analysis. The two case studies show how such elements of arbitrariness of the traditional methods influence both size and location of the identified hotspots, and how this new global method can be used for a more objective selection of hotspots.

Weaver syndrome and EZH2 mutations: Clarifying the clinical phenotype.

Weaver syndrome, first described in 1974, is characterized by tall stature, a typical facial appearance, and variable intellectual disability. In 2011, mutations in the histone methyltransferase, EZH2, were shown to cause Weaver syndrome. To date, we have identified 48 individuals with EZH2 mutations. The mutations were primarily missense mutations occurring throughout the gene, with some clustering in the SET domain (12/48). Truncating mutations were uncommon (4/48) and only identified in the final exon, after the SET domain. Through analyses of clinical data and facial photographs of EZH2 mutation-positive individuals, we have shown that the facial features can be subtle and the clinical diagnosis of Weaver syndrome is thus challenging, especially in older individuals. However, tall stature is very common, reported in >90% of affected individuals. Intellectual disability is also common, present in ~80%, but is highly variable and frequently mild. Additional clinical features which may help in stratifying individuals to EZH2 mutation testing include camptodactyly, soft, doughy skin, umbilical hernia, and a low, hoarse cry. Considerable phenotypic overlap between Sotos and Weaver syndromes is also evident. The identification of an EZH2 mutation can therefore provide an objective means of confirming a subtle presentation of Weaver syndrome and/or distinguishing Weaver and Sotos syndromes. As mutation testing becomes increasingly accessible and larger numbers of EZH2 mutation-positive individuals are identified, knowledge of the clinical spectrum and prognostic implications of EZH2 mutations should improve.

The analysis of ecstasy tablets in a forensic drug intelligence perspective

The use by police services and inquiring agencies of forensic data in an intelligence perspective is still fragmentary and to some extent ignored. In order to increase the efficiency of criminal investigation to target illegal drug trafficking organisations and to provide valuable information about their methods, it is necessary to include and interpret objective drug analysis results already during the investigation phase. The value of visual, physical and chemical data of seized ecstasy tablets, as a support for criminal investigation on a strategic and tactical level has been investigated. In a first phase different characteristics of ecstasy tablets have been studied in order to define their relevance, variation, correlation and discriminating power in an intelligence perspective. During 5 years, over 1200 cases of ecstasy seizures (concerning about 150000 seized tablets) coming from different regions of Switzerland (City and Canton of Zurich, Cantons Ticino, Neuchâtel and Geneva) have been systematically recorded. This turned out to be a statistically representative database including large and small cases. During the second phase various comparison and clustering methods have been tested and evaluated, on the type and relevance of tablet characteristics, thus increasing knowledge about synthetic drugs, their manufacturing and trafficking. Finally analytical methodologies have been investigated and formalised, applying traditional intelligence methods. In this context classical tools, which are used in criminal analysis (like the I2 Analyst Notebook, I2 Ibase, ?) have been tested and adapted to address the specific need of forensic drug intelligence. The interpretation of these links provides valuable information about criminal organisations and their trafficking methods. In the final part of this thesis practical examples illustrate the use and value of such information.

Aging preferentially affects molecular pathways implicated in development and disease of myelinating glial cells

The central and peripheral nervous systems are involved in multiple agedependent neurological deficits that are often attributed to alterations in function of myelinating glial cells. However, the molecular events that underlie the age-related decline of glial cell function are unknown. We used Schwann cells as a model to study biological processes affected in glial cells by aging. We comprehensively profiled gene expression of the Schwann cell-rich mouse sciatic nerve throughout life, from day of birth until senescence (840 days of age). We combined the aging data with the microarray transcriptional data obtained using nerves isolated from Schwann cell-specific neuropathy-inducing mutants MPZCre/þ/Lpin1fE2-3/fE2-3, MPZCre/þ/ScapfE1/fE1 and Pmp22-null mice. A majority of age related transcripts were also affected in the analyzed mouse models of neuropathy (54.4%) and in development (59.5%) indicating a high level of overlapping in implicated molecular pathways. We observed that compared to peripheral nerve development, dynamically changing expression profiles in aging have opposite (anticorrelated) orientation while they copy the orientation of transcriptional changes observed in analyzed neuropathy models. Subsequent clustering and biological annotation of dynamically changing transcripts revealed that the processes most significantly deregulated in aging include inflammatory/ immune response and lipid biosynthesis/metabolism. Importantly, the changes in these pathways were also observed in myelinated oligodendrocyte- rich optic nerves of aged mice, albeit with lower magnitude. This observation suggests that similar biological processes are affected in aging glial cells in central and peripheral nervous systems, however with different dynamics. Our data, which provide the first comprehensive comparison of molecular changes in glial cells in three distinct biological conditions comprising development, aging and disease, provide not only a new inside into the molecular alterations underlying neural system aging but also identify target pathways for potential therapeutical approaches to prevent or delay complications associated with age-related and inherited forms of neuropathies.

Detection of clusters using space-time scan statistics

This paper aims at detecting spatio-temporal clustering in fire sequences using space?time scan statistics, a powerful statistical framework for the analysis of point processes. The methodology is applied to active fire detection in the state of Florida (US) identified by MODIS (Moderate Resolution Imaging Spectroradiometer) during the period 2003?06. Results of the present study show that statistically significant clusters can be detected and localized in specific areas and periods of the year. Three out of the five most likely clusters detected for the entire frame period are localized in the north of the state, and they cover forest areas; the other two clusters cover a large zone in the south, corresponding to agricultural land and the prairies in the Everglades. In order to analyze if the wildfires recur each year during the same period, the analyses have been performed separately for the 4 years: it emerges that clusters of forest fires are more frequent in hot seasons (spring and summer), while in the southern areas, they are widely present during the whole year. The recognition of overdensities of events and the ability to locate them in space and in time can help in supporting fire management and focussing on prevention measures.

Identifying the index lesion with template prostate mapping biopsies.

PURPOSE: The natural history of prostate cancer might be driven by the index lesion. We determined the percent of men in whom the index lesion could be defined using transperineal template prostate mapping biopsies. MATERIALS AND METHODS: Included in study were consecutive men undergoing transperineal template prostate mapping biopsies with biopsies grouped into 20 zones. Men with clinically significant disease in only 1 prostate area were considered to have an identifiable index lesion. We evaluated the impact of using 2 definitions of clinically significant disease (Gleason grade pattern 4 and/or lesion volume 0.5 cc or greater) and 2 clustering rules (stringent and tolerant) to define the index lesion. RESULTS: Included in study were 391 men with a median age of 62 years (IQR 58-67) and a median prostate specific antigen of 6.9 ng/ml (IQR 4.8-10.0). Of the men 269 (69%) were previously diagnosed with prostate cancer. By deploying a median of 1.2 cores per ml (IQR 0.9-1.7) cancer was diagnosed in 82.9% of the men (324 of 391) with a median of 6 positive cores (IQR 2-9), a median maximum cancer core length of 5 mm (IQR 3-8) and a total cancer core length per zone of 7 mm (IQR 3-13). Insignificant disease was found in 26.3% to 42.9% of cases. When a stringent spatial relationship was used to define individual lesions, 44.4% to 54.6% of patients had 1 index lesion and 12.7% to 19.1% had more than 1 area with clinically significant disease. These proportions changed to 46.6% to 59.2% and 10.5% to 14.5%, respectively, when less stringent spatial clustering was applied. CONCLUSIONS: Transperineal template prostate mapping biopsies enable the index lesion to be localized in most men with clinically significant disease. This information may be important to select appropriate candidates for targeted therapy and to plan a tailored treatment strategy in men undergoing radical therapy.

In this issue.

A new issue, once again a bouquet of attractive papers. First of all the paper by Droit-Dupré et al. (10.1007/s00428-015-1724-9). The group studied colonic adenocarcinomas, not otherwise specified, by immunohistochemistry for the expression of markers of intestinal epithelial cell differentiation. Hierarchical clustering analysis identified a major cluster of two thirds of the case series, expressing cytokeratin 20, CDX2 and MUC2 and invariably mismatch repair competent, which they called crypt-like. In stage III colon cancer, the crypt-like cluster had a better prognosis. The paper is a relatively simple example of what is happening in cancer classification beyond morphology: multiparameter differentiation and (epi)genomic markers defining new subtypes of cancer with potential clinical significance in clinical decision making.

Extra Forces induced by wide-pulse, high-frequency electrical stimulation: Occurrence, magnitude, variability and underlying mechanisms.

OBJECTIVE: In contrast to conventional (CONV) neuromuscular electrical stimulation (NMES), the use of "wide-pulse, high-frequencies" (WPHF) can generate higher forces than expected by the direct activation of motor axons alone. We aimed at investigating the occurrence, magnitude, variability and underlying neuromuscular mechanisms of these "Extra Forces" (EF). METHODS: Electrically-evoked isometric plantar flexion force was recorded in 42 healthy subjects. Additionally, twitch potentiation, H-reflex and M-wave responses were assessed in 13 participants. CONV (25Hz, 0.05ms) and WPHF (100Hz, 1ms) NMES consisted of five stimulation trains (20s on-90s off). RESULTS: K-means clustering analysis disclosed a responder rate of almost 60%. Within this group of responders, force significantly increased from 4% to 16% of the maximal voluntary contraction force and H-reflexes were depressed after WPHF NMES. In contrast, non-responders showed neither EF nor H-reflex depression. Twitch potentiation and resting EMG data were similar between groups. Interestingly, a large inter- and intrasubject variability of EF was observed. CONCLUSION: The responder percentage was overestimated in previous studies. SIGNIFICANCE: This study proposes a novel methodological framework for unraveling the neurophysiological mechanisms involved in EF and provides further evidence for a central contribution to EF in responders.

Analysis of the genetic basis of periodic fever with aphthous stomatitis, pharyngitis, and cervical adenitis (PFAPA) syndrome.

PFAPA syndrome is the most common autoinflammatory syndrome in children from Western countries. In spite of its strong familial clustering, its genetic basis and inheritance pattern are still unknown. We performed a comprehensive genetic study on 68 individuals from 14 families. Linkage analysis suggested a susceptibility locus on chromosome 8, but direct molecular sequencing did not support this initial statistical finding. Exome sequencing revealed the absence of any gene that was mutated in all patients. Exhaustive screening of genes involved in other autoinflammatory syndromes or encoding components of the human inflammasome showed no DNA variants that could be linked to PFAPA molecular pathology. Among these, the previously-reported missense mutation V198M in the NLRP3 gene was clearly shown not to co-segregate with PFAPA. Our results on this relatively large cohort indicate that PFAPA syndrome is unlikely to be a monogenic condition. Moreover, none of the several genes known to be involved in inflammation or in autoinflammatory disorders seem to be relevant, alone, to its etiology, suggesting that PFAPA results from oligogenic or complex inheritance of variants in multiple disease genes and/or non-genetic factors.

Mutations in LONP1, a mitochondrial matrix protease, cause CODAS syndrome.

Cerebral, ocular, dental, auricular, skeletal anomalies (CODAS) syndrome (MIM 600373) was first described and named by Shehib et al, in 1991 in a single patient. The anomalies referred to in the acronym are as follows: cerebral-developmental delay, ocular-cataracts, dental-aberrant cusp morphology and delayed eruption, auricular-malformations of the external ear, and skeletal-spondyloepiphyseal dysplasia. This distinctive constellation of anatomical findings should allow easy recognition but despite this only four apparently sporadic patients have been reported in the last 20 years indicating that the full phenotype is indeed very rare with perhaps milder or a typical presentations that are allelic but without sufficient phenotypic resemblance to permit clinical diagnosis. We performed exome sequencing in three patients (an isolated case and a brother and sister sib pair) with classical features of CODAS. Sanger sequencing was used to confirm results as well as for mutation discovery in a further four unrelated patients ascertained via their skeletal features. Compound heterozygous or homozygous mutations in LONP1 were found in all (8 separate mutations; 6 missense, 1 nonsense, 1 small in-frame deletion) thus establishing the genetic basis of CODAS and the pattern of inheritance (autosomal recessive). LONP1 encodes an enzyme of bacterial ancestry that participates in protein turnover within the mitochondrial matrix. The mutations cluster at the ATP-binding and proteolytic domains of the enzyme. Biallelic inheritance and clustering of mutations confirm dysfunction of LONP1 activity as the molecular basis of CODAS but the pathogenesis remains to be explored.