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.

Exact form factors for the Josephson tunneling current and relative particle number fluctuations in a model of two coupled Bose-Einstein condensates

Form factors are derived for a model describing the coherent Josephson tunneling between two coupled Bose-Einstein condensates. This is achieved by studying the exact solution of the model within the framework of the algebraic Bethe ansatz. In this approach the form factors are expressed through determinant representations which are functions of the roots of the Bethe ansatz equations.

Estimating Two-Stage Models for Genetic Influences on Alcohol, Tobacco or Drug Use Initiation and Dependence Vulnerability in Twin and Family Data

Genetic research on risk of alcohol, tobacco or drug dependence must make allowance for the partial overlap of risk-factors for initiation of use, and risk-factors for dependence or other outcomes in users. Except in the extreme cases where genetic and environmental risk-factors for initiation and dependence overlap completely or are uncorrelated, there is no consensus about how best to estimate the magnitude of genetic or environmental correlations between Initiation and Dependence in twin and family data. We explore by computer simulation the biases to estimates of genetic and environmental parameters caused by model misspecification when Initiation can only be defined as a binary variable. For plausible simulated parameter values, the two-stage genetic models that we consider yield estimates of genetic and environmental variances for Dependence that, although biased, are not very discrepant from the true values. However, estimates of genetic (or environmental) correlations between Initiation and Dependence may be seriously biased, and may differ markedly under different two-stage models. Such estimates may have little credibility unless external data favor selection of one particular model. These problems can be avoided if Initiation can be assessed as a multiple-category variable (e.g. never versus early-onset versus later onset user), with at least two categories measurable in users at risk for dependence. Under these conditions, under certain distributional assumptions., recovery of simulated genetic and environmental correlations becomes possible, Illustrative application of the model to Australian twin data on smoking confirmed substantial heritability of smoking persistence (42%) with minimal overlap with genetic influences on initiation.

Genetic variability in cultivated common bean beyond the two major gene pools

It is generally accepted that two major gene pools exist in cultivated common bean (Phaseolus vulgaris L.), a Middle American and an Andean one. Some evidence, based on unique phaseolin morphotypes and AFLP analysis, suggests that at least one more gene pool exists in cultivated common bean. To investigate this hypothesis, 1072 accessions from a common bean core collection from the primary centres of origin, held at CIAT, were investigated. Various agronomic and morphological attributes (14 categorical and 11 quantitative) were measured. Multivariate analyses, consisting of homogeneity analysis and clustering for categorical data, clustering and ordination techniques for quantitative data and nonlinear principal component analysis for mixed data, were undertaken. The results of most analyses supported the existence of the two major gene pools. However, the analysis of categorical data of protein types showed an additional minor gene pool. The minor gene pool is designated North Andean and includes phaseolin types CH, S and T; lectin types 312, Pr, B and K; and mostly A5, A6 and A4 types alpha-amylase inhibitor. Analysis of the combined categorical data of protein types and some plant categorical data also suggested that some other germplasm with C type phaseolin are distinguished from the major gene pools.

A family of textilinin genes, two of which encode proteins with antihaemorrhagic properties

Two peptides, textilinins 1 and 2, isolated from the venom of the Australian common brown snake, Pseudonaja textilis textilis, are effective in preventing blood loss. To further investigate the potential of textilinins as anti-haemorrhagic agents, we cloned cDNAs encoding these proteins. The isolated full-length cDNA (430 bp in size) was shown to code for a 59 amino acid protein, corresponding in size to the native peptide, plus an additional 24 amino acid propeptide. Six such cDNAs were identified, differing in nucleotide sequence in the coding region but with an identical propeptide. All six sequences predicted peptides containing six conserved cysteines common to Kunitz-type serine protease inhibitors. When expressed as glutathione S-transferase (GST) fusion proteins and released by cleavage with thrombin, only those peptides corresponding to textilinin 1 and 2 were active in inhibiting plasmin with K-i values similar to those of their native counterparts and in binding to plasmin less tightly than aprotinin by two orders of magnitude. Similarly, in the mouse tail vein blood loss model only recombinant textilinin 1 and 2 were effective in reducing blood loss. These recombinant textilinins have potential as therapeutic agents for reducing blood loss in humans, obviating the need for reliance on aprotinin, a bovine product with possible risk of transmissible disease, and compromising the fibrinolytic system in a less irreversible manner.

The E(NK) model: extending the NK model to incorporate gene-by-environment interactions and epistasis for diploid genomes

The haploid NK model developed by Kauffman can be extended to diploid genomes and to incorporate gene-by-environment interaction effects in combination with epistasis. To provide the flexibility to include a wide range of forms of gene-by-environment interactions, a target population of environment types (TPE) is defined. The TPE consists of a set of E different environment types, each with their own frequency of occurrence. Each environment type conditions a different NK gene network structure or series of gene effects for a given network structure, providing the framework for defining gene-by-environment interactions. Thus, different NK models can be partially or completely nested within the E environment types of a TPE, giving rise to the E(NK) model for a biological system. With this model it is possible to examine how populations of genotypes evolve in context with properties of the environment that influence the contributions of genes to the fitness values of genotypes. We are using the E(NK) model to investigate how both epistasis and gene-by-environment interactions influence the genetic improvement of quantitative traits by plant breeding strategies applied to agricultural systems. © 2002 Wiley Periodicals, Inc.

Genes induced by growth hormone in a model of adipogenic differentiation

A substantial number of GH regulated genes have been reported in mature hepatocytes. but genes involved in GH-initiated cell differentiation have not yet been identified. Here we have studied a, ell-characterised model of GH-dependent differentiation, adipogenesis of 3T3-F442A preadipocytes, to identify genes rapidly induced by GH. Using the suppression subtractive hybridisation technique, we have identified eight genes induced within 60 min of GH treatment, and verified these by northern analysis. Six were identifiable as Stat 2. Stat 3, thrombospondin-1. oncostatin M receptor beta chain. a DEAD box RNA helicase. and muscleblind. a developmental transcription factor. Bioinformatic approaches assigned one of the two remaining unknown genes as a novel 436 residue serine,threonine kinase. As each of the identified genes hake important developmental roles. they may be important in initiating GH-induced adipogenesis. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Structural characterization of respiratory syncytial virus fusion inhibitor escape mutants: homology model of the F protein and a syncytium formation assay

Respiratory syncytial virus (RSV) is a ubiquitous human pathogen and the leading cause of lower respiratory tract infections in infants. Infection of cells and subsequent formation of syncytia occur through membrane fusion mediated by the RSV fusion protein (RSV-F). A novel in vitro assay of recombinant RSV-F function has been devised and used to characterize a number of escape mutants for three known inhibitors of RSV-F that have been isolated. Homology modeling of the RSV-F structure has been carried out on the basis of a chimera derived from the crystal structures of the RSV-F core and a fragment from the orthologous fusion protein from Newcastle disease virus (NDV). The structure correlates well with the appearance of RSV-F in electron micrographs, and the residues identified as contributing to specific binding sites for several monoclonal antibodies are arranged in appropriate solvent-accessible clusters. The positions of the characterized resistance mutants in the model structure identify two promising regions for the design of fusion inhibitors. (C) 2003 Elsevier Science (USA). All rights reserved.

The application of wave induced forces to a two-dimensional finite element long wave hydrodynamic model

This paper describes the modification of a two-dimensional finite element long wave hydrodynamic model in order to predict the net current and water levels attributable to the influences of waves. Tests examine the effects of the application of wave induced forces, including comparisons to a physical experiment. An example of a real river system is presented with comparisons to measured data, which demonstrate the importance of simulating the combined effects of tides and waves upon hydrodynamic behavior. (C) 2002 Elsevier Science Ltd. All rights reserved.

A two-Higgs doublet model with remarkable CP properties

We analyze generalized CP symmetries of two-Higgs doublet models, extending them from the scalar to the fermion sector of the theory. We show that, other than the usual CP transformation, there is only one of those symmetries which does not imply massless charged fermions. That single model which accommodates a fermionic mass spectrum compatible with experimental data possesses a remarkable feature. Through a soft breaking of the symmetry it displays a new type of spontaneous CP violation, which does not occur in the scalar sector responsible for the symmetry breaking mechanism but, rather, in the fermion sector.

Basis invariant conditions for supersymmetry in the two-Higgs-doublet model

The minimal supersymmetric standard model involves a rather restrictive Higgs potential with two Higgs fields. Recently, the full set of classes of symmetries allowed in the most general two-Higgs-doublet model was identified; these classes do not include the supersymmetric limit as a particular class. Thus, a physically meaningful definition of the supersymmetric limit must involve the interaction of the Higgs sector with other sectors of the theory. Here we show how one can construct basis invariant probes of supersymmetry involving both the Higgs sector and the gaugino-Higgsino-Higgs interactions.

Two-Higgs-doublet model in light of the standard model H -> tau(+)tau(-) search at the LHC

We study the implications of the searches based on H -> tau(+)tau-by the ATLAS and CMS collaborations on the parameter space of the two-Higgs-doublet model (2HDM). In the 2HDM, the scalars can decay into a tau pair with a branching ratio larger than the SM one, leading to constraints on the 2HDM parameter space. We show that in model II, values of tan beta > 1.8 are definitively excluded if the pseudoscalar is in the mass range 110 GeV < m(A) < 145 GeV. We have also discussed the implications for the 2HDM of the recent dimuon search by the ATLAS collaboration for a CP-odd scalar in the mass range 4-12 GeV.

Identificação de compostos que modulam a patogénese da doença de Machado-Joseph em C.elegans: autofagia como alvo terapêutico: identification of small compounds that modulate pathogenesis of Machado-Joseph disease in a C.elegans model: targeting autophagy

A doença de Machado-Joseph (DMJ) ou ataxia espinocerebelosa do tipo 3 (SCA3), conhecida por ser a mais comum das ataxias hereditárias dominantes em todo o mundo, é uma doença neurodegenerativa autossómica dominante que leva a uma grande incapacidade motora, embora sem alterar o intelecto, culminando com a morte do doente. Atualmente não existe nenhum tratamento eficaz para esta doença. A DMJ é resultado de uma alteração genética causada pela expansão de uma sequência poliglutamínica (poliQ), na região C-terminal do gene que codifica a proteína ataxina-3 (ATXN3). Os mecanismos celulares das doenças de poliglutaminas que provocam toxicidade, bem como a função da ATXN3, não são ainda totalmente conhecidos. Neste trabalho, usamos, pela sua simplicidade e potencial genético, um pequeno animal invertebrado, o nemátode C. elegans, com o objetivo de identificar fármacos eficazes para o combate contra a patogénese da DMJ, analisando simultaneamente o seu efeito na agregação da ATXN3 mutante nas células neuronais in vivo e o seu impacto no comportamento motor dos animais. Este pequeno invertebrado proporciona grandes vantagens no estudo dos efeitos tóxicos de proteínas poliQ nos neurónios, uma vez que a transparência das suas 959 células (das quais 302 são neurónios) facilita a deteção de proteínas fluorescentes in vivo. Para além disso, esta espécie tem um ciclo de vida curto, é económica e de fácil manutenção. Neste trabalho testámos no nosso modelo transgénico da DMJ com 130Qs em C.elegans dois compostos potencialmente moduladores da agregação da ATXN3 mutante e da resultante disfunção neurológica, atuando pela via da autofagia. De modo a validar a possível importância terapêutica da ativação da autofagia os compostos candidatos escolhidos foram o Litío e o análogo da Rapamicina CCI-779, testados independentemente e em combinação. A neuroproteção conferida pelo Litío e pelo CCI-779 independentemente sugere que o uso destes fármacos possa ser considerado uma boa estratégia como terapia para a DMJ, a testar em organismos evolutivamente mais próximos do humano. A manipulação da autofagia, segundo vários autores, parece ser benéfica e pode ser a chave para o desenvolvimento de novos tratamentos para várias doenças relacionadas com a agregação proteica e o envelhecimento.

Abelian symmetries in the two-Higgs-doublet model with fermions

We classify all possible implementations of an Abelian symmetry in the two-Higgs-doublet model with fermions. We identify those symmetries which are consistent with nonvanishing quark masses and a Cabibbo-Kobayashi-Maskawa quark-mixing matrix (CKM), which is not block-diagonal. Our analysis takes us from a plethora of possibilities down to 246 relevant cases, requiring only 34 distinct matrix forms. We show that applying Z(n) with n >= 4 to the scalar sector leads to a continuous U(1) symmetry in the whole Lagrangian. Finally, we address the possibilities of spontaneous CP violation and of natural suppression of the flavor-changing neutral currents. We explain why our work is relevant even for non-Abelian symmetries.

Geometric picture of generalized-CP and Higgs-family transformations in the two-Higgs-doublet model

In the two-Higgs-doublet model (THDM), generalized-CP transformations (phi(i) -> X-ij phi(*)(j) where X is unitary) and unitary Higgs-family transformations (phi(i) -> U-ij phi(j)) have recently been examined in a series of papers. In terms of gauge-invariant bilinear functions of the Higgs fields phi(i), the Higgs-family transformations and the generalized-CP transformations possess a simple geometric description. Namely, these transformations correspond in the space of scalar-field bilinears to proper and improper rotations, respectively. In this formalism, recent results relating generalized CP transformations with Higgs-family transformations have a clear geometric interpretation. We will review what is known regarding THDM symmetries, as well as derive new results concerning those symmetries, namely how they can be interpreted geometrically as applications of several CP transformations.