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.

DNA vaccine coding for the full-length infectious Kunjin virus RNA protects mice against the New York strain of West Nile virus

A plasmid DNA directing transcription of the infectious full-length RNA genome of Kunjin (KUN) virus in vivo from a mammalian expression promoter was used to vaccinate mice intramuscularly. The KUN viral cDNA encoded in the plasmid contained the mutation in the NS1 protein (Pro-250 to Leu) previously shown to attenuate KUN virus in weanling mice. KUN virus was isolated from the blood of immunized mice 3-4 days after DNA inoculation, demonstrating that infectious RNA was being transcribed in vivo; however, no symptoms of virus-induced disease were observed. By 19 days postimmunization, neutralizing antibody was detected in the serum of immunized animals. On challenge with lethal doses of the virulent New York strain of West Nile (WN) or wild-type KUN virus intracerebrally or intraperitoneally, mice immunized with as little as 0.1-1 mug of KUN plasmid DNA were solidly protected against disease. This finding correlated with neutralization data in vitro showing that serum from KUN DNA-immunized mice neutralized KUN and WN,viruses with similar efficiencies. The results demonstrate that delivery of an attenuated but replicating KUN virus via a plasmid DNA vector may provide an effective vaccination strategy against virulent strains of WN virus.

Zygosity Diagnosis in the Absence of Genotypic Data: An Approach Using Latent Class Analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Systematic characterization of mutations in yeast acetohydroxyacid synthase: Interpretation of herbicide-resistance data

Acetohydroxyacid synthase (AHAS, EC 4.1.3.18) catalyses the first step in branched-chain amino acid biosynthesis and is the target for sulfonylurea and imidazolinone herbicides, which act as potent and specific inhibitors. Mutants of the enzyme have been identified that are resistant to particular herbicides. However, the selectivity of these mutants towards various sulfonylureas and imidazolinones has not been determined systematically. Now that the structure of the yeast enzyme is known, both in the absence and presence of a bound herbicide, a detailed understanding of the molecular interactions between the enzyme and its inhibitors becomes possible. Here we construct 10 active mutants of yeast AHAS, purify the enzymes and determine their sensitivity to six sulfonylureas and three imidazolinones. An additional three active mutants were constructed with a view to increasing imidazolinone sensitivity. These three variants were purified and tested for their sensitivity to the imidazolinones only. Substantial differences are observed in the sensitivity of the 13 mutants to the various inhibitors and these differences are interpreted in terms of the structure of the herbicide-binding site on the enzyme.

The characteristics of six species of living hollow-bearing trees and their importance for arboreal marsupials in the dry sclerophyll forests of southeast Queensland, Australia

Six species of trees located in the dry sclerophyll forests of southeast Queensland were studied to ascertain which was most suitable to be retained as hollow-bearing trees for nesting and denning by arboreal marsupials. Generally for all tree species, the number of entrances to hollows was positively correlated with the diameter at breast height (DBH) and the growth stage, and entrance diameters also increased in trees with a larger DBH. However, there were differences between the species; Corymbia citriodora had few hollows until the individuals were very large while Eucalyptus crebra had low numbers of hollows throughout its entire size range. It was concluded that a mixture of tree species provided a range of hollow sizes and positions that would be suitable for nesting and denning by arboreal marsupials in those forests. There were large differences between tree species in the relationship between tree size and estimated age. Five of the tree species took between 186 and 230 years to begin to produce hollows while E. crebra took up to 324 years. This suggests that tree species other than E. crebra may be the most preferred for retention in areas where hollow-bearing tree densities are lower than the prescribed level. Other data also suggests there are likely to be enough trees in larger size classes that would begin to form hollows within the next 50 years to compensate for an expected loss of hollow-bearing stags during that same period. In terms of forest operation, the retention of six hollow-bearing trees/ha would represent an estimated loss of 7.3-15% wood production. (C) 2003 Elsevier Science B.V. All rights reserved.

Obesity is associated with worse peritoneal dialysis outcomes in the Australia and New Zealand patient populations

Although obesity is associated with increased risks of morbidity and death in the general population, a number of studies of patients undergoing hemodialysis have demonstrated that increasing body mass index (BMI) is correlated with decreased mortality risk. Whether this association holds true among patients treated with peritoneal dialysis (PD) has been less well studied. The aim of this investigation was to examine the association between BMI and outcomes among new PD patients in a large cohort, with long-term follow-up monitoring. Using data from the Australia and New Zealand Dialysis and Transplant Registry, an analysis of all new adult patients (n = 9679) who underwent an episode of PD treatment in Australia or New Zealand between April 1, 1991, and March 31, 2002, was performed. Patients were classified as obese (BMI of greater than or equal to30 kg/m(2)), overweight (BMI of 25.0 to 29.9 kg/m(2)), normal weight (BMI of 20 to 24.9 kg/m(2)), or underweight (BMI of

Potential for early involvement of CYP isoforms in aspects of human cadmium toxicity

This paper investigates the possible link between non-workplace cadmium (Cd) exposure, cytochrome P450 expression and hypertension. We present results of our investigation into the relationships between liver and kidney Cd burdens and the abundance of the CYP isoform 4A11. Our data show associations between non-workplace Cd exposure and changes in the abundance of hepatic and renal cortical CYP4A11. In liver the levels of immunochemically detectable CYP4A11 were positively correlated with tissue Cd content while in contrast CYP4A11 abundance was inversely correlated with kidney Cd burden. These differences are most likely related to the different Cd burden of the tissues. These observations suggest the potential for involvement of Cd as a mediator of CYP4A11 expression in kidney cortex and indicate that elevations in kidney Cd content may be involved in hypertension via alteration of the expression of this particular isoform. Potential mechanisms by which Cd may alter CYP4A11 expression are discussed briefly. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Solubility of selected esters in supercritical carbon dioxide

The solubility of ethyl propionate, ethyl butyrate, and ethyl isovalerate in supercritical carbon dioxide was measured at temperature ranging from 308.15 to 333.15 K and pressure ranging from 85 to 195 bar. At the same temperature, the solubility of these compounds increases with pressure. The crossover pressure region was also observed in this study. The experimental data were correlated by the semi-empirical Chrastil equation and Peng-Robinson equation of state (EOS) using several mixing rules. The Peng-Robinson EOS gives better solubility prediction than the empirical Chrastil equation. (C) 2002 Elsevier Science B.V. All rights reserved.