Biomedical subjects

Grier P Page

Publications and source records attributed to Grier P Page.

At least 19 recordsLinked to original sources

PubMed2025

Aberrant DNA methylation of genes regulating CD4+ T cell HIV-1 reservoir in women with HIV.

BACKGROUND: The HIV-1 reservoir in CD4+ T cells (HRCD4) pose a major challenge to curing HIV, with many of its mechanisms still unclear. HIV-1 DNA integration and immune responses may alter the host's epigenetic landscape, potentially silencing HIV-1 replication. METHODS: This study used bisulphite capture DNA methylation sequencing in CD4+ T cells from the blood of 427 virally suppressed women with HIV to identify differentially methylated sites and regions associated with HRCD4. RESULTS: The average total HRCD4 size was 1409 copies per million cells, with most proviruses defective and only a small proportion intact. The study identified 245 differentially methylated CpG sites and 85 regions linked to HRCD4 size, with 52% of significant sites in intronic regions. Genes associated with HRCD4 were involved in viral replication, HIV-1 latency and cell growth and apoptosis. HRCD4 size was inversely related to DNA methylation of interferon signalling genes and positively associated with methylation at known HIV-1 integration sites. HRCD4-associated genes were enriched on the pathways related to immune defence, transcription repression and host-virus interactions. CONCLUSIONS: These findings suggest that HIV-1 reservoir is linked to aberrant DNA methylation in CD4+ T cells, offering new insights into epigenetic mechanisms of HIV-1 latency and potential molecular targets for eradication strategies. KEY POINTS: Study involved 427 women with HIV. Identified 245 aberrant DNA methylation sites and 85 methylation regions in CD4+ T cells linked to the HIV-1 reservoir. Highlighted genes are involved in viral replication, immune defence, and host genome integration. Findings suggest potential molecular targets for eradication strategies.

Humans↗

PubMed2006

Impact of neurofibromatosis 1 on Quality of Life: a cross-sectional study of 176 American cases.

Neurofibromatosis 1 (NF1), a genetic condition most commonly characterized by the presence of dermal neurofibromas and café au lait macules, has a significant impact on Quality of Life (QoL). There is a wide range of phenotypic variability, so that affected individuals may have either medically devastating or relatively mild manifestations that do not impact their daily lives. In this study, the SF-36 and Skindex questionnaires were used to quantitatively investigate the impact of severity and visibility on QoL in an American population. Participants were recruited primarily through advertisements distributed by The Children's Tumor Foundation (CTF). The majority participated by completing a mailed questionnaire in which they rated themselves using Ablon's visibility and Riccardi's severity scales, and then completed the SF-36 and Skindex questionnaires. Participants with NF1 reported a significant impact in all aspects of skin-disease-specific QoL, but the emotional aspect demonstrated the greatest effect. Participants with more visible signs of NF1 reported significantly greater overall effects on their skin-disease-specific QoL than those whose manifestations were more subtle. All domains of general health QoL were affected as well, especially in participants who reported having severe complications. Interestingly, greater visibility was not found to be associated with a significant decrease in general health QoL. These findings are consistent with those found in a French cohort, and demonstrate the utility of incorporating tools designed for use in both the general population and for patients with skin disease in examining the impact of NF1 on QoL.

Adult↗

PubMed2006

Transcriptional coordination of the metabolic network in Arabidopsis.

Patterns of coexpression can reveal networks of functionally related genes and provide deeper understanding of processes requiring multiple gene products. We performed an analysis of coexpression networks for 1,330 genes from the AraCyc database of metabolic pathways in Arabidopsis (Arabidopsis thaliana). We found that genes associated with the same metabolic pathway are, on average, more highly coexpressed than genes from different pathways. Positively coexpressed genes within the same pathway tend to cluster close together in the pathway structure, while negatively correlated genes typically occupy more distant positions. The distribution of coexpression links per gene is highly skewed, with a small but significant number of genes having numerous coexpression partners but most having fewer than 10. Genes with multiple connections (hubs) tend to be single-copy genes, while genes with multiple paralogs are coexpressed with fewer genes, on average, than single-copy genes, suggesting that the network expands through gene duplication, followed by weakening of coexpression links involving duplicate nodes. Using a network-analysis algorithm based on coexpression with multiple pathway members (pathway-level coexpression), we identified and prioritized novel candidate pathway members, regulators, and cross pathway transcriptional control points for over 140 metabolic pathways. To facilitate exploration and analysis of the results, we provide a Web site (http://www.transvar.org/at_coexpress/analysis/web) listing analyzed pathways with links to regression and pathway-level coexpression results. These methods and results will aid in the prioritization of candidates for genetic analysis of metabolism in plants and contribute to the improvement of functional annotation of the Arabidopsis genome.

Arabidopsis↗

PubMed2006

Surfactant dysfunction in SP-A-/- and iNOS-/- mice with mycoplasma infection.

Surfactant dysfunction was studied in C57BL/6 (B6), B6.SP-A(-/-), and B6.iNOS(-/-) mice with pulmonary mycoplasma infection (10(7) colony-forming units). Cell-free bronchoalveolar lavage (BAL) from uninfected B6.SP-A(-/-) versus B6 mice had a reduced content of very large aggregates (VLA) and an increase in intermediate large aggregates (ILA), with no difference in total large aggregates (LA = VLA + ILA). However, LA from uninfected B6.SP-A(-/-) versus B6 mice contained less protein and were more sensitive to inhibition by serum albumin and lysophosphatidylcholine in pulsating bubble studies in vitro. Infection with Mycoplasma pulmonis caused significant lung injury and surfactant abnormalities in B6.SP-A(-/-), B6.iNOS(-/-), and B6 mice at 24, 48, 72 h after infection compared with uninfected mice of the same strain. Analyses of time-pooled data indicated that mycoplasma-infected B6.SP-A(-/-) and B6.iNOS(-/-) mice had significantly lower levels of LA and higher protein/phospholipid ratios in BAL compared with infected B6 mice. Infected B6.iNOS(-/-) versus B6 mice also had increased minimum surface tensions on the pulsating bubble and decreased levels of surfactant protein (SP)-B in BAL. These results indicate that pulmonary mycoplasma infection in vivo causes lung injury and surfactant abnormalities that are dependent in part on iNOS and SP-A. In addition, SP-A deficiency modifies surfactant aggregate content and lowers the inhibition resistance of LA surfactant in vitro compared with congenic normal mice.

Animals↗

PubMed2006

S-adenosylmethionine prevents chronic alcohol-induced mitochondrial dysfunction in the rat liver.

An early event that occurs in response to alcohol consumption is mitochondrial dysfunction, which is evident in changes to the mitochondrial proteome, respiration defects, and mitochondrial DNA (mtDNA) damage. S-adenosylmethionine (SAM) has emerged as a potential therapeutic for treating alcoholic liver disease through mechanisms that appear to involve decreases in oxidative stress and proinflammatory cytokine production as well as the alleviation of steatosis. Because mitochondria are a source of reactive oxygen/nitrogen species and a target for oxidative damage, we tested the hypothesis that SAM treatment during alcohol exposure preserves organelle function. Mitochondria were isolated from livers of rats fed control and ethanol diets with and without SAM for 5 wk. Alcohol feeding caused a significant decrease in state 3 respiration and the respiratory control ratio, whereas SAM administration prevented these alcohol-mediated defects and preserved hepatic SAM levels. SAM treatment prevented alcohol-associated increases in mitochondrial superoxide production, mtDNA damage, and inducible nitric oxide synthase induction, without a significant lessening of steatosis. Accompanying these indexes of oxidant damage, SAM prevented alcohol-mediated losses in cytochrome c oxidase subunits as shown using blue native PAGE proteomics and immunoblot analysis, which resulted in partial preservation of complex IV activity. SAM treatment attenuated the upregulation of the mitochondrial stress chaperone prohibitin. Although SAM supplementation did not alleviate steatosis by itself, SAM prevented several key alcohol-mediated defects to the mitochondria genome and proteome that contribute to the bioenergetic defect in the liver after alcohol consumption. These findings reveal new molecular targets through which SAM may work to alleviate one critical component of alcohol-induced liver injury: mitochondria dysfunction.

Animals↗

PubMed2006

Hedgehog signaling and response to cyclopamine differ in epithelial and stromal cells in benign breast and breast cancer.

The hedgehog pathway regulates epithelial-mesenchymal interactions, differentiation, proliferation and survival during development. Stimulation of hedgehog signaling induces carcinogenesis or promotes cell survival in cancers of multiple organs. Using real-time, quantitative PCR, laser capture microdissection, and immunohistochemistry, distinctive patterns of expression of the hedgehog pathway members patched 1 (PTCH1), smoothened, GLI1, GLI2 and the 3 hedgehog ligands were identified for epithelial cells and stromal fibroblasts in benign breast and breast cancer. Hedgehog ligands were expressed at higher levels in some cancer epithelial cell lines compared to noncancerous epithelial cells. Correspondingly, expression of GLI1, a transcription factor and transcriptional product of hedgehog signaling, was increased 8-fold in cancer epithelial cell lines; however, PTCH1, also a transcriptional target of hedgehog signaling in many cell types, was not increased. GLI1 protein and mRNA, and PTCH1 and sonic hedgehog (SHH) proteins were elevated in 3 of 10 breast cancers; however, PTCH1 transcripts were not consistently increased. Hedgehog-mediated transcription, as indicated by a reporter of GLI-dependent promoter activity and by expression of GLI1 transcripts, was reduced by the hedgehog pathway inhibitor cyclopamine in both MDA-MB-435 cancer epithelial cells and MCF10AT epithelial cells, a cell line derived from benign breast. However, cyclopamine reduced viability of cancer epithelial cell lines, including MDA-MB-435, but did not specifically affect fibroblasts or epithelial cells from benign breast, including MCF10AT. Treatment with sonic hedgehog ligand diminished the cyclopamine-induced reduction in GLI-dependent promoter activity in MCF10AT and MDA-MB-435 and viability of MDA-MB-435. These results demonstrate modulation of GLI-mediated transcription in both cancer and benign-derived epithelial cells by cyclopamine and sonic hedgehog, and further suggest that hedgehog signaling contributes to the survival of only the cancer epithelial cells. Determination as to whether the increase in GLI1 and SHH expression in breast cancer indicates a significant increase in hedgehog signaling will require further evaluation.

Base Sequence↗

PubMed2006

The PowerAtlas: a power and sample size atlas for microarray experimental design and research.

BACKGROUND: Microarrays permit biologists to simultaneously measure the mRNA abundance of thousands of genes. An important issue facing investigators planning microarray experiments is how to estimate the sample size required for good statistical power. What is the projected sample size or number of replicate chips needed to address the multiple hypotheses with acceptable accuracy? Statistical methods exist for calculating power based upon a single hypothesis, using estimates of the variability in data from pilot studies. There is, however, a need for methods to estimate power and/or required sample sizes in situations where multiple hypotheses are being tested, such as in microarray experiments. In addition, investigators frequently do not have pilot data to estimate the sample sizes required for microarray studies. RESULTS: To address this challenge, we have developed a Microrarray PowerAtlas. The atlas enables estimation of statistical power by allowing investigators to appropriately plan studies by building upon previous studies that have similar experimental characteristics. Currently, there are sample sizes and power estimates based on 632 experiments from Gene Expression Omnibus (GEO). The PowerAtlas also permits investigators to upload their own pilot data and derive power and sample size estimates from these data. This resource will be updated regularly with new datasets from GEO and other databases such as The Nottingham Arabidopsis Stock Center (NASC). CONCLUSION: This resource provides a valuable tool for investigators who are planning efficient microarray studies and estimating required sample sizes.

Algorithms↗

PubMed2006

Genome scan for loci predisposing to anxiety disorders using a novel multivariate approach: strong evidence for a chromosome 4 risk locus.

We conducted a 10-centimorgan linkage autosomal genome scan in a set of 19 extended American pedigrees (219 subjects) ascertained through probands with panic disorder. Several anxiety disorders--including social phobia, agoraphobia, and simple phobia--in addition to panic disorder segregate in these families. In previous studies of this sample, linkage analyses were based separately on each of the individual categorical affection diagnoses. Given the substantial comorbidity between anxiety disorders and their probable shared genetic liability, it is clear that this method discards a considerable amount of information. In this article, we propose a new approach that considers panic disorder, simple phobia, social phobia, and agoraphobia as expressions of the same multivariate, putatively genetically influenced trait. We applied the most powerful multipoint Haseman-Elston method, using the grade of membership score generated from a fuzzy clustering of these phenotypes as the dependent variable in Haseman-Elston regression. One region on chromosome 4q31-q34, at marker D4S413 (with multipoint and single-point nominal P values < .00001), showed strong evidence of linkage (genomewide significance at P<.05). The same region is known to be the site of a neuropeptide Y receptor gene, NPY1R (4q31-q32), that was recently connected to anxiolytic-like effects in rats. Several other regions on four chromosomes (4q21.21-22.3, 5q14.2-14.3, 8p23.1, and 14q22.3-23.3) met criteria for suggestive linkage (multipoint nominal P values < .01). Family-by-family analysis did not show any strong evidence of heterogeneity. Our findings support the notion that the major anxiety disorders, including phobias and panic disorder, are complex traits that share at least one susceptibility locus. This method could be applied to other complex traits for which shared genetic-liability factors are thought to be important, such as substance dependencies.

Anxiety↗

PubMed2006

A proposed metric for assessing the measurement quality of individual microarrays.

BACKGROUND: High-density microarray technology is increasingly applied to study gene expression levels on a large scale. Microarray experiments rely on several critical steps that may introduce error and uncertainty in analyses. These steps include mRNA sample extraction, amplification and labeling, hybridization, and scanning. In some cases this may be manifested as systematic spatial variation on the surface of microarray in which expression measurements within an individual array may vary as a function of geographic position on the array surface. RESULTS: We hypothesized that an index of the degree of spatiality of gene expression measurements associated with their physical geographic locations on an array could indicate the summary of the physical reliability of the microarray. We introduced a novel way to formulate this index using a statistical analysis tool. Our approach regressed gene expression intensity measurements on a polynomial response surface of the microarray's Cartesian coordinates. We demonstrated this method using a fixed model and presented results from real and simulated datasets. CONCLUSION: We demonstrated the potential of such a quantitative metric for assessing the reliability of individual arrays. Moreover, we showed that this procedure can be incorporated into laboratory practice as a means to set quality control specifications and as a tool to determine whether an array has sufficient quality to be retained in terms of spatial correlation of gene expression measurements.

Algorithms↗

PubMed2006

Optimal allocation of replicates for measurement evaluation studies.

Optimal experimental design is important for the efficient use of modern high-throughput technologies such as microarrays and proteomics. Multiple factors including the reliability of measurement system, which itself must be estimated from prior experimental work, could influence design decisions. In this study, we describe how the optimal number of replicate measures (technical replicates) for each biological sample (biological replicate) can be determined. Different allocations of biological and technical replicates were evaluated by minimizing the variance of the ratio of technical variance (measurement error) to the total variance (sum of sampling error and measurement error). We demonstrate that if the number of biological replicates and the number of technical replicates per biological sample are variable, while the total number of available measures is fixed, then the optimal allocation of replicates for measurement evaluation experiments requires two technical replicates for each biological replicate. Therefore, it is recommended to use two technical replicates for each biological replicate if the goal is to evaluate the reproducibility of measurements.

Computational Biology↗

PubMed2006

Microarray data analysis: from disarray to consolidation and consensus.

In just a few years, microarrays have gone from obscurity to being almost ubiquitous in biological research. At the same time, the statistical methodology for microarray analysis has progressed from simple visual assessments of results to a weekly deluge of papers that describe purportedly novel algorithms for analysing changes in gene expression. Although the many procedures that are available might be bewildering to biologists who wish to apply them, statistical geneticists are recognizing commonalities among the different methods. Many are special cases of more general models, and points of consensus are emerging about the general approaches that warrant use and elaboration.

Algorithms↗

PubMed2006

Energy restriction lowers the expression of genes linked to inflammation, the cytoskeleton, the extracellular matrix, and angiogenesis in mouse adipose tissue.

Using high-density oligonucleotide microarrays, we examined the actions of energy restriction (ER) on the expression of >11,000 genes in epididymal white adipose tissue (WAT) of 10- to 11-mo-old male C57Bl6 mice. Four groups were studied: controls not subjected to food restriction (CO), food-restricted 18 h before being killed (FR), short-term ER for 23 d (SER), and long-term ER for 9 mo (LER). As we reported previously, compared with CO mice, FR and SER minimally influenced the gene expression profiles; however, 345 transcripts of 6,266 genes determined to be expressed in WAT were significantly altered by LER. We focus here on the 109 (31%) of these genes that were involved in either inflammation (56 genes), cytoskeleton (16 genes), extracellular matrix (23 genes), or angiogenesis (14 genes). Among these 109 genes, 104 transcripts (95%) were down regulated by LER. Western blotting for heat shock protein 47 and osteonectin, and immunohistochemical staining for hypoxia inducible factor (HIF)-1alpha), supported the microarray data that LER down regulated the expressions of these genes. Additionally, a 75% reduction in adipocyte size with LER reflected the change in the expression of genes involved in cell morphology. Our findings provide evidence that LER suppresses the expression of genes encoding inflammatory molecules in WAT while promoting structural remodeling of the cytoskeleton, extracellular matrix, and vasculature. These alterations may play an important role in the protection against WAT-derived inflammation and in lifespan extension by LER.

Adipose Tissue↗

PubMed2006

Familial and genetic susceptibility to major neonatal morbidities in preterm twins.

BACKGROUND: Intraventricular hemorrhage, necrotizing enterocolitis, and bronchopulmonary dysplasia remain significant causes of morbidity and mortality in preterm newborns. OBJECTIVES: Our goal was to assess the familial and genetic susceptibility to intraventricular hemorrhage, necrotizing enterocolitis, and bronchopulmonary dysplasia. METHODS: Mixed-effects logistic-regression and latent variable probit model analysis were used to assess the contribution of several covariates in a multicenter retrospective study of 450 twin pairs born at < or =32 weeks of gestation. To determine the genetic contribution, concordance rates in a subset of 252 monozygotic and dizygotic twin pairs were compared. RESULTS: The study population had a mean gestational age of 29 weeks and birth weight of 1286 g. After controlling for effects of covariates, the twin data showed that 41.3%, 51.9%, and 65.2%, respectively, of the variances in liability for intraventricular hemorrhage, necrotizing enterocolitis, and bronchopulmonary dysplasia could be accounted for by genetic and shared environmental factors. Among the 63 monozygotic twin pairs, the observed concordance for bronchopulmonary dysplasia was significantly higher than the expected concordance; 12 of 18 monozygotic twin pairs with > or =1 affected member had both members affected versus 3.69 expected. After controlling for covariates, genetic factors accounted for 53% of the variance in liability for bronchopulmonary dysplasia. CONCLUSIONS: Twin analyses show that intraventricular hemorrhage, necrotizing enterocolitis, and bronchopulmonary dysplasia are familial in origin. These data demonstrate, for the first time, the significant genetic susceptibility for bronchopulmonary dysplasia in preterm infants.

Bronchopulmonary Dysplasia↗

PubMed2005

COGA phenotypes and linkages on chromosome 2.

An initial linkage analysis of the alcoholism phenotype as defined by DSM-III-R criteria and alcoholism defined by DSM-IV criteria showed many, sometimes striking, inconsistencies. These inconsistencies are greatly reduced by making the definition of alcoholism more specific. We defined new phenotypes combining the alcoholism definitions and the latent variables, defining an individual as affected if that individual is alcoholic under one of the definitions (either DSM-III-R or DSM-IV), and indicated having a symptom defined by one of the latent variables. This was done for each of the two alcoholism definitions and five latent variables, selected from a canonical discriminant analyses indicating they formed significant groupings using the electrophysiological variables. We found that linkage analyses utilizing these latent variables were much more robust and consistent than the linkage results based on DSM-III-R or DSM-IV criteria for definition of alcoholism. We also performed linkage analyses on two first principal components derived phenotypes, one derived from the electrophysiological variables, and the other derived from the latent variables. A region on chromosome 2 at 250 cM was found to be linked to both of these derived phenotypes. Further examination of the SNPs in this region identified several haplotypes strongly associated with these derived phenotypes.

Alcoholism↗

PubMed2005

DCDC2 is associated with reading disability and modulates neuronal development in the brain.

DYX2 on 6p22 is the most replicated reading disability (RD) locus. By saturating a previously identified peak of association with single nucleotide polymorphism markers, we identified a large polymorphic deletion that encodes tandem repeats of putative brain-related transcription factor binding sites in intron 2 of DCDC2. Alleles of this compound repeat are in significant disequilibrium with multiple reading traits. RT-PCR data show that DCDC2 localizes to the regions of the brain where fluent reading occurs, and RNA interference studies show that down-regulation alters neuronal migration. The statistical and functional studies are complementary and are consistent with the latest clinical imaging data for RD. Thus, we propose that DCDC2 is a candidate gene for RD.

Adult↗

PubMed2005

TDT-association analysis of EKN1 and dyslexia in a Colorado twin cohort.

A candidate gene, EKN1, was recently described in a cohort from Finland for the dyslexia locus on chromosome 15q, DYX1. This report described a (2;15) (q11;21) translocation disrupting EKN1 that cosegregated with dyslexia in a two-generation family. It also characterized a sequence polymorphism in the 5' untranslated region and a missense mutation that showed significant association in 109 dyslexics compared to 195 controls (p=0.002 and p=0.006, respectively). To confirm these results we interrogated the same polymorphisms in a cohort of 150 nuclear families with dyslexia ascertained through the Colorado Learning Disabilities Research Center. Using QTDT analysis with nine individual quantitative tasks and two composite measures of reading performance, we could not replicate the reported association. We conclude that the polymorphisms identified in the Finland sample are unlikely to be functional DNA changes contributing to dyslexia, and that if variation in EKN1 is causal such changes are more likely to be in regulatory regions that were not sequenced in this study. Alternatively, the published findings of association with markers in EKN1 may reflect linkage disequilibrium with variation in another gene(s) in the region.

Base Sequence↗

PubMed2005

Sources of variation in Affymetrix microarray experiments.

BACKGROUND: A typical microarray experiment has many sources of variation which can be attributed to biological and technical causes. Identifying sources of variation and assessing their magnitude, among other factors, are important for optimal experimental design. The objectives of this study were: (1) to estimate relative magnitudes of different sources of variation and (2) to evaluate agreement between biological and technical replicates. RESULTS: We performed a microarray experiment using a total of 24 Affymetrix GeneChip arrays. The study included 4th mammary gland samples from eight 21-day-old Sprague Dawley CD female rats exposed to genistein (soy isoflavone). RNA samples from each rat were split to assess variation arising at labeling and hybridization steps. A general linear model was used to estimate variance components. Pearson correlations were computed to evaluate agreement between technical and biological replicates. CONCLUSION: The greatest source of variation was biological variation, followed by residual error, and finally variation due to labeling when *.cel files were processed with dChip and RMA image processing algorithms. When MAS 5.0 or GCRMA-EB were used, the greatest source of variation was residual error, followed by biology and labeling. Correlations between technical replicates were consistently higher than between biological replicates.

Animals↗

PubMed2005

Reproducible clusters from microarray research: whither?

MOTIVATION: In cluster analysis, the validity of specific solutions, algorithms, and procedures present significant challenges because there is no null hypothesis to test and no 'right answer'. It has been noted that a replicable classification is not necessarily a useful one, but a useful one that characterizes some aspect of the population must be replicable. By replicable we mean reproducible across multiple samplings from the same population. Methodologists have suggested that the validity of clustering methods should be based on classifications that yield reproducible findings beyond chance levels. We used this approach to determine the performance of commonly used clustering algorithms and the degree of replicability achieved using several microarray datasets. METHODS: We considered four commonly used iterative partitioning algorithms (Self Organizing Maps (SOM), K-means, Clutsering LARge Applications (CLARA), and Fuzzy C-means) and evaluated their performances on 37 microarray datasets, with sample sizes ranging from 12 to 172. We assessed reproducibility of the clustering algorithm by measuring the strength of relationship between clustering outputs of subsamples of 37 datasets. Cluster stability was quantified using Cramer's v2 from a kXk table. Cramer's v2 is equivalent to the squared canonical correlation coefficient between two sets of nominal variables. Potential scores range from 0 to 1, with 1 denoting perfect reproducibility. RESULTS: All four clustering routines show increased stability with larger sample sizes. K-means and SOM showed a gradual increase in stability with increasing sample size. CLARA and Fuzzy C-means, however, yielded low stability scores until sample sizes approached 30 and then gradually increased thereafter. Average stability never exceeded 0.55 for the four clustering routines, even at a sample size of 50. These findings suggest several plausible scenarios: (1) microarray datasets lack natural clustering structure thereby producing low stability scores on all four methods; (2) the algorithms studied do not produce reliable results and/or (3) sample sizes typically used in microarray research may be too small to support derivation of reliable clustering results. Further research should be directed towards evaluating stability performances of more clustering algorithms on more datasets specially having larger sample sizes with larger numbers of clusters considered.

Cluster Analysis↗