A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.

Patterns and levels of endemism in the Australian Wet Tropics rainforest: evidence from flightless insects

This study examines the level and pattern of endemism among 274 flightless rainforest insects found in the Wet Tropics region of Australia. Endemism is measured at two nested scales: (1) those confined to the Wet Tropics, termed 'regional endemics'; and (2) the subset of those species confined to a single subregion of the Wet Tropics, termed 'subregional endemics'. Fifty per cent of the regional endemic flightless insects are also subregional endemics compared with 15% of the known regional endemic vertebrates. The four subregions with the most endemic flightless insect species are the uplands of Mt Finnigan, Carbine, Bellenden-Ker/Bartle Frere and Atherton. Multiple regression suggests that the combination of rainforest area and shape explain the most variance (R-2 = 0.603) in the numbers of species of regional endemic insects. However, subregional endemism is not closely correlated with the size or shape of the subregions in which they occur, or a combination of these factors. Candidate refugial and recolonised subregions are identified, and are consistent with data from palaeoclimatic models and refugia identified using other taxa. We group upland subregions into larger areas of endemism using parsimony analysis of endemism. These groupings are consistent with our understanding of the history of the Wet Tropics rainforests.

A test of the hypothesis of ecological equivalence in an Australian subtropical rain forest

We tested the hypothesis that tree species in a subtropical rain forest in south-east Queensland are ecologically equivalent and therefore have identical environmental requirements for their regeneration. We assessed the evidence that juveniles of species differed in their distributions in treefall gap microsites and along gradients of light availability, soil pH, soil PO4-P availability and soil NO3-N availability. Pairwise comparisons were made on a subset of the common species selected on the basis that they showed a relatively high level of positive association, and would therefore, a priori, be expected to have similar regeneration requirements. Detailed comparisons between the species failed to demonstrate evidence for species differentiation with respect to their tolerance of the disturbance associated with gap microsites or to the gradient of NO3-N availability. However, species differed markedly in their distributions along the soil pH gradient and along the gradients of light availability and soil PO4-P availability. The overall level of ecological differentiation between the species is high: seven out of the 10 possible species pairings showed evidence for ecological differentiation. Such niche differentiation amongst the juveniles of tree species may play an important role in maintaining the species richness of rain-forest communities.

Characterization of the mouse hnRNP A2/B1/B0 gene and identification of processed pseudogenes

The mouse hnRNP A2/B1/B0 gene has been cloned using a PCR-based strategy and sequenced. Analysis of this sequence showed that the gene organization closely follows that of the human orthologue with 12 exons and 11 introns. The hnRNP A2/B1/B0 gene gives rise to four splice variants through alternative splicing of exons 2 and 9. RT-PCR assays indicated that all splice variants were expressed in mouse brain, skin, and stomach tissues of varying ages, although their ratios to one another varied with age and tissue type. We also identified a small subset of all polyadenylated splice variants that included intron 11, which shows 94% sequence identity between human and mouse. Several processed pseudogenes were identified in the mouse genome. A search of the mouse genome databases located five pseudogenes, four of. which are presumed to be non-functional because of the presence of premature stop codons, large deletions or rearrangements within the coding region. The fifth, which possesses putative promoter elements and has a coding sequence identical to that of the hnRNP A2 mRNA, variant, may be functional. (C) 2002 Elsevier Science B.V. All rights reserved.

Alterations in Plasmodium falciparum genotypes during sequential infections suggest the presence of strain specific immunity

Many of the asexual stage Plasmodium falciparum proteins that are the targets of host protective responses are markedly polymorphic. The full repertoire of diversity is not defined for any antigen. Most studies have focused on the genes encoding merozoite surface proteins 1 and 2 (MSP1, MSP2). We explored the extent of diversity of some of the less studied merozoite surface antigens and analyzed the degree of complexity of malaria field isolates by deriving nucleotide sequences of several antigens. We have determined the genotype of apical membrane antigen 1 (AMA1) in a group of 30 field samples, collected over 29 months, from individuals living in an area of intense malaria transmission in Irian Jaya, identifying 14 different alleles. AMA1 genotyping was combined with previously determined MSP2 typing. AMA1 had the greatest power in distinguishing between isolates but methodological problems, especially when mixed infections are present, suggest it is not an ideal typing target. MSP1, MSP3, and glutamate-rich protein genotypes were also determined from a smaller group of samples, and all results were combined to derive an extended antigenic haplotype. Within this subset of 10 patients, nine different genotypes could be discerned; however, five patients were all infected with the same strain. This strain was present in individuals from two separate villages and was still present 12 months later. This strain was predominant at the first time point but had disappeared at the fourth time point. This significant change in malaria genotypes could be due to strain-specific immunity developing in this population.

Cadmium levels in the lung, liver, kidney cortex, and urine samples from Australians without occupational exposure to metals

The authors undertook this study to assess levels of cadmium exposure in the general population. Samples of lung, liver, and kidney were obtained from 61 cadavers (43 males, 18 females; 2-89 yr of age, mean age = 38.5 yr) who died from accidental causes and who were subject to postmortem examinations at the John Tonge Centre for Forensic Sciences, Queensland Health Scientific Services, Brisbane, Australia, in 1997 and 1998. Samples of bladder urine were also obtained from 22 cadavers. Tissue and urine samples were analyzed for cadmium, zinc, and copper with inductively coupled plasm (ICP) mass spectrometry. The overall mean values for cadmium in the lung, liver, and kidney cortex samples were 0.13, 0.95, and 15.45 mug/gm wet tissue weight. The average renal cadmium level in subjects with high lung-cadmium levels (n = 13) was 6 mug/gm wet tissue weight higher than that of similarly aged subjects who had medium lung-cadmium levels (n = 30). In females, the average level of cadmium in the liver was 74% greater than in males, and the average liver cadmium in females with high lung-cadmium levels was 100% higher than in males in the same age range who had the same high lung-cadmium levels. Renal cadmium accumulation tended to be greater in females than in males who were in the same age range and who had similar lung-cadmium levels, a result that suggested that there was a higher absorption rate of cadmium in females. The mean value for a urinary cadmium excretion of 2.30 mug/gm creatinine was found in a subset of samples that had a mean age of 39 yr and a renal cortex cadmium concentration of 18.6 mug/gm wet tissue weight. Urinary cadmium excretion rates were correlated more strongly with lung and kidney cadmium content than with age or liver cadmium levels. The results suggest that urinary cadmium excretion may be increased in smokers and could provide some estimate of body cadmium burdens in future Australian epidemiological studies.

Expression of specific glycoconjugates in both primary and secondary olfactory pathways in BALB/C mice

Binding of cell surface carbohydrates to their receptors specifically promotes axon growth and synaptogenesis in select regions of the developing nervous system. In some cases these interactions depend upon cell-cell adhesion mediated by the same glycoconjugates present on the surface of apposing cells or their processes. We have previously shown that the plant lectin Dolichos biflorus agglutinin (DBA) binds to: a subpopulation of mouse primary olfactory neurons whose axons selectively fasciculate prior to terminating in the olfactory bulb. In the present study, we investigated whether these glycoconjugates were also expressed by postsynaptic olfactory neurons specifically within the olfactory pathway. We show here for the first time that DBA ligands were expressed both by a subset of primary olfactory neurons as well as by the postsynaptic mitral/tufted cells in BALB/C mice. These glycoconjugates were first detected on mitral/tufted cell axons during the early postnatal period, at a time when there is considerable synaptogenesis and synaptic remodelling in the primary olfactory cortex. This is one of the few examples of the selective expression of molecules in contiguous axon tracts in the mammalian nervous system. These results suggest that glycoconjugates recognized by DBA may have a specific role in the formation and maintenance of neural connections within a select functional pathway in the brain. J. Comp. Neurol. 443:213-225, 2002. (C) 2002 Wiley-Liss, Inc.

Expression and role of roundabout-1 in embryonic Xenopus forebrain

The receptor Roundabout-1 (Robo1) and its ligand Slit are known to influence axon guidance and central nervous system (CNS) patterning in both vertebrate and nonvertebrate systems. Although Robo-Slit interactions mediate axon guidance in the Drosophila CNS, their role in establishing the early axon scaffold in the embryonic vertebrate brain remains unclear. We report here the identification and expression of a Xenopus Robo1 orthologue that is highly homologous to mammalian Robo1. By using overexpression studies and immunohistochemical and in situ hybridization techniques, we have investigated the role of Robo1 in the development of a subset of neurons and axon tracts in the Xenopus forebrain. Robo1 is expressed in forebrain nuclei and in neuroepithelial cells underlying the main axon tracts. Misexpression of Robo1 led to aberrant development of axon tracts as well as the ectopic differentiation of forebrain neurons. These results implicate Robo1 in both neuronal differentiation and axon guidance in embryonic vertebrate forebrain. (C) 2002 Wiley-Liss, Inc.

Angiotensinogen genotype, plasma protein and mRNA concentration in isolated systolic hypertension

Isolated systolic hypertension (ISH) occurs predominantly in the elderly, with a considerable morbidity and mortality. Its etiology is unknown but is likely to involve a significant genetic component. The aim of this study was to examine the angiotensinogen gene in ISH. The M235T and G(- 6)A polymorphisms were genotyped by polymerase chain reaction (PCR) in 86 ISH patients and 120 normotensive controls. Plasma angiotensinogen concentration was determined in 198 subjects by an indirect radioimmunoassay technique. Angiotensinogen mRNA concentration was determined by quantitative competitive reverse transcription (RT)-PCR in subcutaneous adipose tissue from a subset of these patients (n = 8) and controls (n = 6). Both the M235T (p = 0.0015) and G(- 6)A (p = 0.029) polymorphisms were associated with ISH. Plasma angiotensinogen concentration was higher in patients than controls (p < 0.0001), but was not associated with genotype. Angiotensinogen mRNA concentration in adipose tissue from ISH subjects was significantly lower than in adipose tissue from normotensive subjects (p = 0.033). The association of angiotensinogen gene variants with ISH and the elevation of plasma angiotensinogen concentration in these patients suggests a role of the angiotensinogen gene in this form of hypertension. Angiotensinogen gene expression may be altered in ISH, but this requires further examination.

T cell regulation of the immune response to infection in periodontal diseases

Although T cells have been implicated in the pathogenesis and are considered to be central both in progression and control of the chronic inflammatory periodontal diseases, the precise contribution of T cells to the regulation of tissue destruction has not been fully elucidated. Current dogma suggests that immunity to infection is controlled by distinct T helper 1 (Th1) and T helper 2 (Th2) subsets of T cells classified on the basis of their cytokine profile. Further, a subset of T cells with immunosuppressive function and cytokine profile distinct from Th1 or Th2 has been described and designated as regulatory T cells. Although these regulatory T cells have been considered to maintain self-tolerance resulting in the suppression of auto-immune responses, recent data suggest that these cells may also play a role in preventing infection-induced immunopathology. In this review, the role of functional and regulatory T cells in chronic inflammatory periodontal diseases will be summarized. This should not only provide an insight into the relationship between the immune response to periodontopathic bacteria and disease but should also highlight areas of development for potentially new therapeutic modalities.

Implicit vector equilibrium problems with applications

Let X and Y be Hausdorff topological vector spaces, K a nonempty, closed, and convex subset of X, C: K--> 2(Y) a point-to-set mapping such that for any x is an element of K, C(x) is a pointed, closed, and convex cone in Y and int C(x) not equal 0. Given a mapping g : K --> K and a vector valued bifunction f : K x K - Y, we consider the implicit vector equilibrium problem (IVEP) of finding x* is an element of K such that f (g(x*), y) is not an element of - int C(x) for all y is an element of K. This problem generalizes the (scalar) implicit equilibrium problem and implicit variational inequality problem. We propose the dual of the implicit vector equilibrium problem (DIVEP) and establish the equivalence between (IVEP) and (DIVEP) under certain assumptions. Also, we give characterizations of the set of solutions for (IVP) in case of nonmonotonicity, weak C-pseudomonotonicity, C-pseudomonotonicity, and strict C-pseudomonotonicity, respectively. Under these assumptions, we conclude that the sets of solutions are nonempty, closed, and convex. Finally, we give some applications of (IVEP) to vector variational inequality problems and vector optimization problems. (C) 2003 Elsevier Science Ltd. All rights reserved.

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.

Localization of a novel melanoma susceptibility locus to 1p22

Over the past 20 years, the incidence of cutaneous malignant melanoma (CMM) has increased dramatically worldwide. A positive family history of the disease is among the most established risk factors for CMM; it is estimated that 10% of CMM cases result from an inherited predisposition. Although mutations in two genes, CDKN2A and CDK4, have been shown to confer an increased risk of CMM, they account for only 20%-25% of families with multiple cases of CMM. Therefore, to localize additional loci involved in melanoma susceptibility, we have performed a genomewide scan for linkage in 49 Australian pedigrees containing at least three CMM cases, in which CDKN2A and CDK4 involvement has been excluded. The highest two-point parametric LOD score (1.82; recombination fraction [theta] 0.2) was obtained at D1S2726, which maps to the short arm of chromosome 1 (1p22). A parametric LOD score of 4.65 (theta = 0) and a nonparametric LOD score of 4.19 were found at D1S2779 in nine families selected for early age at onset. Additional typing yielded seven adjacent markers with LOD scores 13 in this subset, with the highest parametric LOD score, 4.95 (theta = 0) ( nonparametric LOD score 5.37), at D1S2776. Analysis of 33 additional multiplex families with CMM from several continents provided further evidence for linkage to the 1p22 region, again strongest in families with the earliest mean age at diagnosis. A nonparametric ordered sequential analysis was used, based on the average age at diagnosis in each family. The highest LOD score, 6.43, was obtained at D1S2779 and occurred when the 15 families with the earliest ages at onset were included. These data provide significant evidence of a novel susceptibility gene for CMM located within chromosome band 1p22.

A molecular mechanism for mRNA trafficking in neuronal dendrites

Specific neuronal mRNAs are localized in dendrites, often concentrated in dendritic spines and spine synapses, where they are translated. The molecular mechanism of localization is mostly unknown. Here we have explored the roles of A2 response element (A2RE), a cis-acting signal for oligodendrocyte RNA trafficking, and its cognate trans-acting factor, heterogeneous nuclear ribonucleoprotein ( hnRNP) A2, in neurons. Fluorescently labeled chimeric RNAs containing A2RE were microinjected into hippocampal neurons, and RNA transport followed using confocal laser scanning microscopy. These RNA molecules, but not RNA lacking the A2RE sequence, were transported in granules to the distal neurites. hnRNP A2 protein was implicated as the cognate trans-acting factor: it was colocalized with RNA in cytoplasmic granules, and RNA trafficking in neurites was compromised by A2RE mutations that abrogate hnRNP A2 binding. Coinjection of antibodies to hnRNP A2 halved the number of trafficking cells, and treatment of neurons with antisense oligonucleotides also disrupted A2RE - RNA transport. Colchicine inhibited trafficking, whereas cells treated with cytochalasin were unaffected, implicating involvement of microtubules rather than microfilaments. A2RE-like sequences are found in a subset of dendritically localized mRNAs, which, together with these results, suggests that a molecular mechanism based on this cis-acting sequence may contribute to dendritic RNA localization.

The t-SNARE syntaxin 4 is regulated during macrophage activation to function in membrane traffic and cytokine secretion

Activation of macrophages with lipopolysaccharide (LPS) induces the rapid synthesis and secretion of proinflammatory cytokines, such as tumor necrosis factor (TNFalpha), for priming the immune response [1, 2]. TNFalpha plays a key role in inflammatory disease [3]; yet, little is known of the intracellular trafficking events leading to its secretion. In order to identify molecules involved in this secretory pathway, we asked whether any of the known trafficking proteins are regulated by LPS. We found that the levels of SNARE proteins were rapidly and significantly up- or downregulated during macrophage activation. A subset of t-SNAREs (Syntaxin 4/SNAP23/Munc18c) known to control regulated exocytosis in other cell types [4, 5] was substantially increased by LPS in a temporal pattern coinciding with peak TNFalpha secretion. Syntaxin 4 formed a complex with Munc18c at the cell surface of macrophages. Functional studies involving the introduction of Syntaxin 4 cDNA or peptides into macrophages implicate this t-SNARE in a rate-limiting step of TNFalpha secretion and in membrane ruffling during macrophage activation. We conclude that in macrophages, SNAREs are regulated in order to accommodate the rapid onset of cytokine secretion and for membrane traffic associated with the phenotypic changes of immune activation. This represents a novel regulatory role for SNAREs in regulated secretion and in macrophage-mediated host defense.