Data-driven clustering : new methods and applications

Abstract : This work is concerned with the development and application of novel unsupervised learning methods, having in mind two target applications: the analysis of forensic case data and the classification of remote sensing images. First, a method based on a symbolic optimization of the inter-sample distance measure is proposed to improve the flexibility of spectral clustering algorithms, and applied to the problem of forensic case data. This distance is optimized using a loss function related to the preservation of neighborhood structure between the input space and the space of principal components, and solutions are found using genetic programming. Results are compared to a variety of state-of--the-art clustering algorithms. Subsequently, a new large-scale clustering method based on a joint optimization of feature extraction and classification is proposed and applied to various databases, including two hyperspectral remote sensing images. The algorithm makes uses of a functional model (e.g., a neural network) for clustering which is trained by stochastic gradient descent. Results indicate that such a technique can easily scale to huge databases, can avoid the so-called out-of-sample problem, and can compete with or even outperform existing clustering algorithms on both artificial data and real remote sensing images. This is verified on small databases as well as very large problems. Résumé : Ce travail de recherche porte sur le développement et l'application de méthodes d'apprentissage dites non supervisées. Les applications visées par ces méthodes sont l'analyse de données forensiques et la classification d'images hyperspectrales en télédétection. Dans un premier temps, une méthodologie de classification non supervisée fondée sur l'optimisation symbolique d'une mesure de distance inter-échantillons est proposée. Cette mesure est obtenue en optimisant une fonction de coût reliée à la préservation de la structure de voisinage d'un point entre l'espace des variables initiales et l'espace des composantes principales. Cette méthode est appliquée à l'analyse de données forensiques et comparée à un éventail de méthodes déjà existantes. En second lieu, une méthode fondée sur une optimisation conjointe des tâches de sélection de variables et de classification est implémentée dans un réseau de neurones et appliquée à diverses bases de données, dont deux images hyperspectrales. Le réseau de neurones est entraîné à l'aide d'un algorithme de gradient stochastique, ce qui rend cette technique applicable à des images de très haute résolution. Les résultats de l'application de cette dernière montrent que l'utilisation d'une telle technique permet de classifier de très grandes bases de données sans difficulté et donne des résultats avantageusement comparables aux méthodes existantes.

Phenotypic novelty in experimental hybrids is predicted by the genetic distance between species of cichlid fish.

BACKGROUND: Transgressive segregation describes the occurrence of novel phenotypes in hybrids with extreme trait values not observed in either parental species. A previously experimentally untested prediction is that the amount of transgression increases with the genetic distance between hybridizing species. This follows from QTL studies suggesting that transgression is most commonly due to complementary gene action or epistasis, which become more frequent at larger genetic distances. This is because the number of QTLs fixed for alleles with opposing signs in different species should increase with time since speciation provided that speciation is not driven by disruptive selection. We measured the amount of transgression occurring in hybrids of cichlid fish bred from species pairs with gradually increasing genetic distances and varying phenotypic similarity. Transgression in multi-trait shape phenotypes was quantified using landmark-based geometric morphometric methods. RESULTS: We found that genetic distance explained 52% and 78% of the variation in transgression frequency in F1 and F2 hybrids, respectively. Confirming theoretical predictions, transgression when measured in F2 hybrids, increased linearly with genetic distance between hybridizing species. Phenotypic similarity of species on the other hand was not related to the amount of transgression. CONCLUSION: The commonness and ease with which novel phenotypes are produced in cichlid hybrids between unrelated species has important implications for the interaction of hybridization with adaptation and speciation. Hybridization may generate new genotypes with adaptive potential that did not reside as standing genetic variation in either parental population, potentially enhancing a population's responsiveness to selection. Our results make it conceivable that hybridization contributed to the rapid rates of phenotypic evolution in the large and rapid adaptive radiations of haplochromine cichlids.

Geometric morphometric analyses provide evidence for the adaptive character of the Tanganyikan cichlid fish radiations.

The cichlids of East Africa are renowned as one of the most spectacular examples of adaptive radiation. They provide a unique opportunity to investigate the relationships between ecology, morphological diversity, and phylogeny in producing such remarkable diversity. Nevertheless, the parameters of the adaptive radiations of these fish have not been satisfactorily quantified yet. Lake Tanganyika possesses all of the major lineages of East African cichlid fish, so by using geometric morphometrics and comparative analyses of ecology and morphology, in an explicitly phylogenetic context, we quantify the role of ecology in driving adaptive speciation. We used geometric morphometric methods to describe the body shape of over 1000 specimens of East African cichlid fish, with a focus on the Lake Tanganyika species assemblage, which is composed of more than 200 endemic species. The main differences in shape concern the length of the whole body and the relative sizes of the head and caudal peduncle. We investigated the influence of phylogeny on similarity of shape using both distance-based and variance partitioning methods, finding that phylogenetic inertia exerts little influence on overall body shape. Therefore, we quantified the relative effect of major ecological traits on shape using phylogenetic generalized least squares and disparity analyses. These analyses conclude that body shape is most strongly predicted by feeding preferences (i.e., trophic niches) and the water depths at which species occur. Furthermore, the morphological disparity within tribes indicates that even though the morphological diversification associated with explosive speciation has happened in only a few tribes of the Tanganyikan assemblage, the potential to evolve diverse morphologies exists in all tribes. Quantitative data support the existence of extensive parallelism in several independent adaptive radiations in Lake Tanganyika. Notably, Tanganyikan mouthbrooders belonging to the C-lineage and the substrate spawning Lamprologini have evolved a multitude of different shapes from elongated and Lamprologus-like hypothetical ancestors. Together, these data demonstrate strong support for the adaptive character of East African cichlid radiations.

Programming of marginal zone B-cell fate by basic Kruppel-like factor (BKLF/KLF3).

Splenic marginal zone (MZ) B cells are a lineage distinct from follicular and peritoneal B1 B cells. They are located next to the marginal sinus where blood is released. Here they pick up antigens and shuttle the load onto follicular dendritic cells inside the follicle. On activation, MZ B cells rapidly differentiate into plasmablasts secreting antibodies, thereby mediating humoral immune responses against blood-borne type 2 T-independent antigens. As Krüppel-like factors are implicated in cell differentiation/function in various tissues, we studied the function of basic Krüppel-like factor (BKLF/KLF3) in B cells. Whereas B-cell development in the bone marrow of KLF3-transgenic mice was unaffected, MZ B-cell numbers in spleen were increased considerably. As revealed in chimeric mice, this occurred cell autonomously, increasing both MZ and peritoneal B1 B-cell subsets. Comparing KLF3-transgenic and nontransgenic follicular B cells by RNA-microarray revealed that KLF3 regulates a subset of genes that was similarly up-regulated/down-regulated on normal MZ B-cell differentiation. Indeed, KLF3 expression overcame the lack of MZ B cells caused by different genetic alterations, such as CD19-deficiency or blockade of B-cell activating factor-receptor signaling, indicating that KLF3 may complement alternative nuclear factor-κB signaling. Thus, KLF3 is a driving force toward MZ B-cell maturation.

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae).

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.

Programming Hippocampal Neural Stem/Progenitor Cells into Oligodendrocytes Enhances Remyelination in the Adult Brain after Injury.

Demyelinating diseases are characterized by a loss of oligodendrocytes leading to axonal degeneration and impaired brain function. Current strategies used for the treatment of demyelinating disease such as multiple sclerosis largely rely on modulation of the immune system. Only limited treatment options are available for treating the later stages of the disease, and these treatments require regenerative therapies to ameliorate the consequences of oligodendrocyte loss and axonal impairment. Directed differentiation of adult hippocampal neural stem/progenitor cells (NSPCs) into oligodendrocytes may represent an endogenous source of glial cells for cell-replacement strategies aiming to treat demyelinating disease. Here, we show that Ascl1-mediated conversion of hippocampal NSPCs into mature oligodendrocytes enhances remyelination in a diphtheria-toxin (DT)-inducible, genetic model for demyelination. These findings highlight the potential of targeting hippocampal NSPCs for the treatment of demyelinated lesions in the adult brain.

Fetal programming and epigenetic mechanisms in arterial hypertension.

PURPOSE OF REVIEW: To provide an overview of available evidence of the potential role of epigenetics in the pathogenesis of hypertension and vascular dysfunction. RECENT FINDINGS: Arterial hypertension is a highly heritable condition. Surprisingly, however, genetic variants only explain a tiny fraction of the phenotypic variation and the term 'missing heritability' has been coined to describe this phenomenon. Recent evidence suggests that phenotypic alteration that is unrelated to changes in DNA sequence (thereby escaping detection by classic genetic methodology) offers a potential explanation. Here, we present some basic information on epigenetics and review recent work consistent with the hypothesis of epigenetically induced arterial hypertension. SUMMARY: New technologies that enable the rigorous assessment of epigenetic changes and their phenotypic consequences may provide the basis for explaining the missing heritability of arterial hypertension and offer new possibilities for treatment and/or prevention.

The role of geography and ecology in shaping repeated patterns of morphological and genetic differentiation between European minnows (Phoxinus phoxinus) from the Pyrenees and the Alps

Neutral and selective processes c an drive repeated patterns of evolu tion in dif ferent groups of populationsexp eriencing similar ecol ogica l gradients. In this paper, we used a combinat ion of nucl ear and mitochondrialDNA markers, as well as geometric morphometrics, to investigate repeated patterns of morphological andgenetic divergence of E uropean minnows in two mountain ranges : the Pyrenees and the Al ps. Europeanminnows (Phoxinus phoxinus) are cyprinid fish i nha bitin g most freshwater bodies in Europe, including those indifferent mountain r anges that could act as major geographical barriers to gene flow. We explored patterns ofP. phoxinus phenotypic and genetic di versi fication along a gradi ent of alti tude common to the two mountainranges, and tested for isolation by distance (IBD), isolation by environment (IBE) and isolation by adaptation(IBA). The results indicated that populations from the Pyr enees a nd the Alps bel ong to two well differentiated,reciprocally monophyletic mt DNA lineages. Substantial genetic differentiation due to geographical isolationwithin and between populations from the Pyrenees and the Alps was also found using rapidly evolving AFLPsmarkers (isolation by distance or IBD), as well as morphological differences between mountain ranges. Als o,morphology varied strong ly with elevation and so did genetic differentiation to a lower extent. Despitemoderate evidence for IBE and IBA, and therefore of repeated evolution, substantial population heterogeneitywas found at the genetic level, suggesting that selection and population specific genetic drift act in concert toaffect genetic divergence.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Impact of genetic SLC28 transporter and ITPA variants on ribavirin 21 serum level, hemoglobin drop and therapeutic response in patients with HCV infection

Background: T reatment o f chronic hepatitis C i s evolving, a nd direct acting antivirals ( DAAs) are now a dded to p egylated interferon-α ( Peg- INF-α) and ribavirin (RBV) for the treatment o f hepatitis C v irus ( HCV) genotype 1 infection. DAAs c ause d ifferent side effects and can even worsen RBV induced hemolytic anemia. T herefore, identifying host genetic d eterminants of R BV bioavailability and therapeutic e fficacy will remain crucial for individualized treatment. Recent d ata showed associations between R BV induced h emolytic anemia and genetic polymorphisms o f concentrative nucleoside transporters s uch as C NT3 (SLC28A3) and i nosine t riphosphatase (ITPA). T o analyze t he association of genetic variants of SLC28 transporters and ITPA with RBV induced hemolytic anemia and treatment o utcome. Methods: I n our study, 173 patients f rom t he S wiss Hepatitis C C ohort Study and 2 2 patients from Swiss Association for the Study of the Liver study 24 (61% HCV g enotype 1, 3 9% genotypes 2 o r 3) were analyzed for SLC28A2 single nucleotide p olymorphism (SNP) rs11854484, SLC28A3 rs56350726 and SLC28A3 rs10868138 as well as ITPA SNPs rs1127354 and rs7270101. RBV serum levels during treatment were measured in 49 patients. Results: SLC28A2 r s11854484 genotype TT was associated with significantly higher dosage- and body weight-adjusted RBV levels as compared to genotypes TC and CC (p=0.04 and p=0.02 at weeks 4 and 8, respectively). ITPA SNPs rs1127354 and rs7270101 were associated with h emolytic a nemia both in genotype as w ell as i n allelic a nalyses. SLC28A3 rs56350726 genotype TT (vs. AT/AA, RR=2.1; 95% CI 1.1-4.1) as well as the T allele (vs. A; RR=1.8, 95% CI 1.1-3.2) were associated with increased SVR rates. The combined analysis of overall ITPA activity and SLC28 v ariants together revealed n o significant a dditive effects on either treatment-related anemia or SVR. Conclusions: T he newly identified association between RBV serum levels a nd SLC28A2 rs11854484 genotype as well as the replicated association of ITPA and SLC28A3 g enetic p olymorphisms w ith RBV induced hemolytic anemia and treatment r esponse underpin the need for further studies on host genetic d eterminants of R BV bioavailability and therapeutic e fficacy f or individualized treatment of chronic hepatitis C.

Risk prediction of prevalent diabetes in a Swiss population using a weighted genetic score--the CoLaus Study.

AIMS/HYPOTHESIS: Several susceptibility genes for type 2 diabetes have been discovered recently. Individually, these genes increase the disease risk only minimally. The goals of the present study were to determine, at the population level, the risk of diabetes in individuals who carry risk alleles within several susceptibility genes for the disease and the added value of this genetic information over the clinical predictors. METHODS: We constructed an additive genetic score using the most replicated single-nucleotide polymorphisms (SNPs) within 15 type 2 diabetes-susceptibility genes, weighting each SNP with its reported effect. We tested this score in the extensively phenotyped population-based cross-sectional CoLaus Study in Lausanne, Switzerland (n = 5,360), involving 356 diabetic individuals. RESULTS: The clinical predictors of prevalent diabetes were age, BMI, family history of diabetes, WHR, and triacylglycerol/HDL-cholesterol ratio. After adjustment for these variables, the risk of diabetes was 2.7 (95% CI 1.8-4.0, p = 0.000006) for individuals with a genetic score within the top quintile, compared with the bottom quintile. Adding the genetic score to the clinical covariates improved the area under the receiver operating characteristic curve slightly (from 0.86 to 0.87), yet significantly (p = 0.002). BMI was similar in these two extreme quintiles. CONCLUSIONS/INTERPRETATION: In this population, a simple weighted 15 SNP-based genetic score provides additional information over clinical predictors of prevalent diabetes. At this stage, however, the clinical benefit of this genetic information is limited.

Genetic variation in IL28B is associated with chronic hepatitis C and treatment failure: a genome-wide association study.

BACKGROUND & AIMS: Hepatitis C virus (HCV) induces chronic infection in 50% to 80% of infected persons; approximately 50% of these do not respond to therapy. We performed a genome-wide association study to screen for host genetic determinants of HCV persistence and response to therapy. METHODS: The analysis included 1362 individuals: 1015 with chronic hepatitis C and 347 who spontaneously cleared the virus (448 were coinfected with human immunodeficiency virus [HIV]). Responses to pegylated interferon alfa and ribavirin were assessed in 465 individuals. Associations between more than 500,000 single nucleotide polymorphisms (SNPs) and outcomes were assessed by multivariate logistic regression. RESULTS: Chronic hepatitis C was associated with SNPs in the IL28B locus, which encodes the antiviral cytokine interferon lambda. The rs8099917 minor allele was associated with progression to chronic HCV infection (odds ratio [OR], 2.31; 95% confidence interval [CI], 1.74-3.06; P = 6.07 x 10(-9)). The association was observed in HCV mono-infected (OR, 2.49; 95% CI, 1.64-3.79; P = 1.96 x 10(-5)) and HCV/HIV coinfected individuals (OR, 2.16; 95% CI, 1.47-3.18; P = 8.24 x 10(-5)). rs8099917 was also associated with failure to respond to therapy (OR, 5.19; 95% CI, 2.90-9.30; P = 3.11 x 10(-8)), with the strongest effects in patients with HCV genotype 1 or 4. This risk allele was identified in 24% of individuals with spontaneous HCV clearance, 32% of chronically infected patients who responded to therapy, and 58% who did not respond (P = 3.2 x 10(-10)). Resequencing of IL28B identified distinct haplotypes that were associated with the clinical phenotype. CONCLUSIONS: The association of the IL28B locus with natural and treatment-associated control of HCV indicates the importance of innate immunity and interferon lambda in the pathogenesis of HCV infection.

A simple genetic basis for complex social behaviour mediates widespread gene expression differences.

A remarkable social polymorphism is controlled by a single Mendelian factor in the fire ant Solenopsis invicta. A genomic element marked by the gene Gp-9 determines whether workers tolerate one or many fertile queens in their colony. Gp-9 was recently shown to be part of a supergene with two nonrecombining variants, SB and Sb. SB/SB and SB/Sb queens differ in how they initiate new colonies, and in many physiological traits, for example odour and maturation rate. To understand how a single genetic element can affect all these traits, we used a microarray to compare gene expression patterns between SB/SB and SB/Sb queens of three different age classes: 1-day-old unmated queens, 11-day-old unmated queens and mated, fully reproductive queens collected from mature field colonies. The number of genes that were differentially expressed between SB/SB and SB/Sb queens of the same age class was smallest in 1-day-old queens, maximal in 11-day-old queens and intermediate in reproductive queens. Gene ontology analysis showed that SB/SB queens upregulate reproductive genes faster than SB/Sb queens. For all age classes, genes inside the supergene were overrepresented among the differentially expressed genes. Consistent with the hypothesized greater number of transposons in the Sb supergene, 13 transposon genes were upregulated in SB/Sb queens. Viral genes were also upregulated in SB/Sb mature queens, consistent with the known greater parasite load in colonies headed by SB/Sb queens compared with colonies headed by SB/SB queens. Eighteen differentially expressed genes between reproductive queens were involved in chemical signalling. Our results suggest that many genes in the supergene are involved in regulating social organization and queen phenotypes in fire ants.

Evolutionary switch and genetic convergence on rbcL following the evolution of C4 photosynthesis.

Rubisco is responsible for the fixation of CO2 into organic compounds through photosynthesis and thus has a great agronomic importance. It is well established that this enzyme suffers from a slow catalysis, and its low specificity results into photorespiration, which is considered as an energy waste for the plant. However, natural variations exist, and some Rubisco lineages, such as in C4 plants, exhibit higher catalytic efficiencies coupled to lower specificities. These C4 kinetics could have evolved as an adaptation to the higher CO2 concentration present in C4 photosynthetic cells. In this study, using phylogenetic analyses on a large data set of C3 and C4 monocots, we showed that the rbcL gene, which encodes the large subunit of Rubisco, evolved under positive selection in independent C4 lineages. This confirms that selective pressures on Rubisco have been switched in C4 plants by the high CO2 environment prevailing in their photosynthetic cells. Eight rbcL codons evolving under positive selection in C4 clades were involved in parallel changes among the 23 independent monocot C4 lineages included in this study. These amino acids are potentially responsible for the C4 kinetics, and their identification opens new roads for human-directed Rubisco engineering. The introgression of C4-like high-efficiency Rubisco would strongly enhance C3 crop yields in the future CO2-enriched atmosphere.

Host factors and genetic susceptibility to infections due to intracellular bacteria and fastidious organisms.

While genetic polymorphisms play a paramount role in tuberculosis (TB), less is known about their contribution to the severity of diseases caused by other intracellular bacteria and fastidious microorganisms. We searched electronic databases for observational studies reporting on host factors and genetic predisposition to infections caused by intracellular fastidious bacteria published up to 30 May 2014. The contribution of genetic polymorphisms was documented for TB. This includes genetic defects in the mononuclear phagocyte/T helper cell type 1 (Th1) pathway contributing to disseminated TB disease in children and genome-wide linkage analysis (GWAS) in reactivated pulmonary TB in adults. Similarly, experimental studies supported the role of host genetic factors in the clinical presentation of illnesses resulting from other fastidious intracellular bacteria. These include IL-6 -174G/C or low mannose-binding (MBL) polymorphisms, which are incriminated in chronic pulmonary conditions triggered by C. pneumoniae, type 2-like cytokine secretion polymorphisms, which are correlated with various clinical patterns of M. pneumoniae infections, and genetic variation in the NOD2 gene, which is an indicator of tubal pathology resulting from Chamydia trachomatis infections. Monocyte/macrophage migration and T lymphocyte recruitment defects are corroborated to ineffective granuloma formation observed among patients with chronic Q fever. Similar genetic polymorphisms have also been suggested for infections caused by T. whipplei although not confirmed yet. In conclusion, this review supports the paramount role of genetic factors in clinical presentations and severity of infections caused by intracellular fastidious bacteria. Genetic predisposition should be further explored through such as exome sequencing.