Trustworthiness and metrics in visualizing similarity of gene expression

Background

Conventionally, the first step in analyzing the large and high-dimensional data sets measured by microarrays is visual exploration. Dendrograms of hierarchical clustering, self-organizing maps (SOMs), and multidimensional scaling have been used to visualize similarity relationships of data samples. We address two central properties of the methods: (i) Are the visualizations trustworthy, i.e., if two samples are visualized to be similar, are they really similar? (ii) The metric. The measure of similarity determines the result; we propose using a new learning metrics principle to derive a metric from interrelationships among data sets.

Results

The trustworthiness of hierarchical clustering, multidimensional scaling, and the self-organizing map were compared in visualizing similarity relationships among gene expression profiles. The self-organizing map was the best except that hierarchical clustering was the most trustworthy for the most similar profiles. Trustworthiness can be further increased by treating separately those genes for which the visualization is least trustworthy. We then proceed to improve the metric. The distance measure between the expression profiles is adjusted to measure differences relevant to functional classes of the genes. The genes for which the new metric is the most different from the usual correlation metric are listed and visualized with one of the visualization methods, the self-organizing map, computed in the new metric.

Conclusions

The conjecture from the methodological results is that the self-organizing map can be recommended to complement the usual hierarchical clustering for visualizing and exploring gene expression data. Discarding the least trustworthy samples and improving the metric still improves it.

     

Molecular Phylogeny and Systematics of Anoplocephaline Cestodes in Rodents and Lagomorphs

A molecular phylogenetic hypothesis is presented for the anoplocephaline cestodes of placental mammals based on sequence data from the mitochondrial cytochrome c oxidase I (COI) gene, the nuclear-encoded 28S rRNA gene and the internal transcribed spacer region I of rRNA (ITS1). The material consists of 35 species representing nine genera of cestodes, with emphasis on taxa parasitising rodents and lagomorphs in the Holarctic region. The resulting phylogenies show considerable disagreement with earlier systematic and phylogenetic hypotheses derived from morphology. Specifically, the results contradict the view of uterine morphology being the primary determinant of deeper phylogenetic splits within Anoplocephalinae. Also, the role of genital duplication as a means of generic divergence was not found to follow consistently the pattern suggested by earlier hypotheses. Colonisation of novel host lineages has evidently been the predominant mode of diversification in anoplocephaline cestodes of placental mammals; evidence for phyletic co-evolution was obscure. The phylogenies consistently distinguished a large monophyletic group including all species from arvicoline rodents (voles and lemmings), primarily representing the genera Anoplocephaloides Baer, 1923 and Paranoplocephala Lühe, 1910. Phylogenetic relationships within the “arvicoline clade” of cestodes were generally poorly resolved. Consistent support for nodes above and below the unresolved polytomy indicates a rapid radiation involving a nearly simultaneous diversification of many lineages, a scenario also proposed for the arvicoline hosts.     

In vitro experimental studies of sialyl Lewis x and sialyl Lewis a on endothelial and carcinoma cells: crucial glycans on selectin ligands

Extravasation from the blood of malignant tumour cells that form metastasis and leukocytes that go into tissues require contact between selectins and their sialyl Lewis x and sialyl Lewis a (sLex and sLea respectively) decorated ligands. Endothelial cells have been shown to express sLex epitopes in lymph nodes and at sites of inflammation, and this is crucial for the selectin-dependent leukocyte traffic. Besides the ability to synthesize sLex on sialylated N-acetyllactosamine via the action of α(1,3)fucosyltransferase(s), endothelial cells can also degrade sLex to Lewis x through the action of α(2,3)sialidase(s). In addition, several epithelial tumors possess the machinery to synthesize sLex, which facilitates their adhesion to endothelial E- and P-selectin. This revised version was published online in November 2006 with corrections to the Cover Date.     

Characterization of Cyanobacteria by SDS-PAGE of whole-cell proteins and PCR/RFLP of the 16S rRNA gene

Planktonic, filamentous cyanobacterial strains from different genera, both toxic and nontoxic strains, were characterized by SDS-PAGE of whole-cell proteins and PCR/RFLP of the 16S rRNA gene. Total protein pattern analysis revealed the mutual relationships at the genus level. Restriction fragment length polymorphism (RFLP) of the 16S rRNA gene with reference strains proved to be a good method for the cyanobacterial taxonomy. The nonheterocystous strains outgrouped from the nitrogen-fixing ones. With both methods, Aphanizomenon clustered with Anabaena, and Nodularia with Nostoc. In the RFLP study of Anabaena, the neurotoxic strains were identical, but the hepatotoxic ones formed a heterogeneous group. Genetic distances found in the RFLP study were short, confirming that close genotypic relationships underlie considerable diversity among cyanobacterial genera. Received: 16 December 1996 / Accepted: 14 May 1997     

MODE OF SWARMING IN RELATION TO REPRODUCTIVE ISOLATION IN MAYFLIES

Mayfly males swarm, that is they fly in a fixed pattern by a specific object, the swarm marker. Females orientate to the same markers. Leptophlebia marginata mayflies were observed to orientate to two kinds of objects in a single locality in central Finland: to trees and to horizontal pale objects on the ground; when dispersed or moved to the other type of marker, they returned to their former orientation. Tree swarming is by far the most common mode of swarming, but some horizontally orientating populations were found. Sympatric populations are genetically and morphologically distinct, whereas other populations appear to have some gene flow between the swarming types. The tree-swarming mode appears to be primitive and the horizontal mode derived; wind rather than predation is the factor favoring swarming close to the ground. Swarming constitutes an effective mechanism of premating isolation in mayflies.     

Interfertility and genetic variability among European and North American isolates of the basidiomycete fungus Chondrostereum purpureum

The conspecificity of Finnish and western Canadian isolates of the decay fungus Chondrostereum purpureum was investigated by several approaches, including the assessment of genetic variability, mating and progeny analysis, and the analysis of selected phenotypic traits. Eight second-generation single spore strains per fungal isolate pairing were investigated with specific genetic markers developed for both Finnish and Canadian parental isolates. Tests of linkage disequilibrium were used to analyze whether these markers assorted independently among single spore strains. This procedure was similarly applied to the third-generation spore progeny. Finally, global non-metric multidimensional scaling was used to analyze independent random amplified microsatellite marker data to assess the genetic variability of the parental Finnish and Canadian isolates, and their second- and third-generation progeny. Our results revealed that the parental isolates from Finland and western Canada were genetically divergent, but no interfertility barriers were identified between these geographically distant fungi. Furthermore, parental genetic markers used in mating studies demonstrated that second- and third-generation spore progenies underwent normal meiosis and genetic recombination without linkage disequilibrium. Based on this work, the studied C. purpureum isolates from Finland and Canada can be considered as belonging to a single biological species, although genetic and limited phenotypic differentiation was observed.     

Microsatellite-based species identification method for Drosophila virilis group species

Species of the D. virilis group are widely used in evolutionary research, but the individuals of different species are difficult to distinguish from each other morphologically. We constructed a fast and easy microsatellite-based identification method for the species of the group occurring sympatrically in northern Europe. The neighbor joining tree based on 14 microsatellite loci also gave a good resolution of the species divergence pattern in the whole group.     

Genetic architecture of resistance in Daphnia hosts against two species of host-specific parasites

Understanding the genetic architecture of host resistance is key for understanding the evolution of host–parasite interactions. Evolutionary models often assume simple genetics based on few loci and strong epistasis. It is unknown, however, whether these assumptions apply to natural populations. Using a quantitative trait loci (QTL) approach, we explore the genetic architecture of resistance in the crustacean Daphnia magna to two of its natural parasites: the horizontally transmitted bacterium Pasteuria ramosa and the horizontally and vertically transmitted microsporidium Hamiltosporidium tvaerminnensis. These two systems have become models for studies on the evolution of host–parasite interactions. In the QTL panel used here, Daphnia''s resistance to P. ramosa is controlled by a single major QTL (which explains 50% of the observed variation). Resistance to H. tvaerminnensis horizontal infections shows a signature of a quantitative trait based in multiple loci with weak epistatic interactions (together explaining 38% variation). Resistance to H. tvaerminnensis vertical infections, however, shows only one QTL (explaining 13.5% variance) that colocalizes with one of the QTLs for horizontal infections. QTLs for resistance to Pasteuria and Hamiltosporidium do not colocalize. We conclude that the genetics of resistance in D. magna are drastically different for these two parasites. Furthermore, we infer that based on these and earlier results, the mechanisms of coevolution differ strongly for the two host–parasite systems. Only the Pasteuria–Daphnia system is expected to follow the negative frequency-dependent selection (Red Queen) model. How coevolution works in the Hamiltosporidium–Daphnia system remains unclear.     

The elemental composition, the microfibril angle distribution and the shape of the cell cross-section in Norway spruce xylem

Relationships between the elemental composition, the microfibril angle (MFA) distribution and the average shape of the cell cross-section of irrigated-fertilised and untreated Norway spruce (Picea abies [L.] Karst.) earlywood were studied. Sample material was obtained from Flakaliden, Sweden. The elemental composition was studied by determining the relative mass fractions of the elements P, S, Cl, K, Ca and Mn by X-ray fluorescence and by determining the mass absorption coefficients for X-rays. X-ray diffraction was used to determine the MFA distribution and the average shape of the cell cross-section. The latter was also determined by light microscopy. In transition from juvenile wood to mature wood, a decrease of the mode of the MFA distribution from 13°–24° to 3°–6° was connected to a change in the shape of the cell cross-section from circular to rectangular. The irrigation-fertilisation treatment caused no change in the MFA distribution or in the shape of the cell cross-section, whereas the mass absorption coefficient was higher and the density was smaller in irrigated-fertilised wood. Larger proportion of the elements S, Cl and K, but smaller proportion of the element Mn, were observed due to the treatment. The results indicate that the shape of the cell cross-section or the MFA distribution are not directly linked to the growth rate of tracheids or to the nutrient-element content in the xylem and only show notable changes as a function of the cambial age.