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.

Gas permeation in silicon-oxide/polymer (SiOx/PET) barrier films: role of the oxide lattice, nano-defects and macro-defects

We propose a model for permeation in oxide coated gas barrier films. The model accounts for diffusion through the amorphous oxide lattice, nano-defects within the lattice, and macro-defects. The presence of nano-defects indicate the oxide layer is more similar to a nano-porous solid (such as zeolite) than silica glass with respect to permeation properties. This explains why the permeability of oxide coated polymers is much greater, and the activation energy of permeation much lower, than values expected for polymers coated with glass. We have used the model to interpret permeability and activation energies measured for the inert gases (He, Ne and Ar) in evaporated SiOx films of varying thickness (13-70 nm) coated on a polymer substrate. Atomic force and scanning electron microscopy were used to study the structure of the oxide layer. Although no defects could be detected by microscopy, the permeation data indicate that macro-defects (>1 nm), nano-defects (0.3-0.4 nm) and the lattice interstices (<0.3 nm) all contribute to the total permeation. (C) 2002 Elsevier Science B.V. All rights reserved.

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.

Temperature dependence of polaronic transport through single molecules and quantum dots

Motivated by recent experiments on electric transport through single molecules and quantum dots, we investigate a model for transport that allows for significant coupling between the electrons and a boson mode isolated on the molecule or dot. We focus our attention on the temperature-dependent properties of the transport. In the Holstein picture for polaronic transport in molecular crystals the temperature dependence of the conductivity exhibits a crossover from coherent (band) to incoherent (hopping) transport. Here, the temperature dependence of the differential conductance on resonance does not show such a crossover, but is mostly determined by the lifetime of the resonant level on the molecule or dot.

Testing insect movement theory: empirical analysis of pest data routinely collected from agricultural crops

Observations of an insect's movement lead to theory on the insect's flight behaviour and the role of movement in the species' population dynamics. This theory leads to predictions of the way the population changes in time under different conditions. If a hypothesis on movement predicts a specific change in the population, then the hypothesis can be tested against observations of population change. Routine pest monitoring of agricultural crops provides a convenient source of data for studying movement into a region and among fields within a region. Examples of the use of statistical and computational methods for testing hypotheses with such data are presented. The types of questions that can be addressed with these methods and the limitations of pest monitoring data when used for this purpose are discussed. (C) 2002 Elsevier Science B.V. All rights reserved.

Heterogeneous nuclear ribonucleoprotein A3, a novel RNA trafficking response element-binding protein

The cis-acting response element, A2RE, which is sufficient for cytoplasmic mRNA trafficking in oligodendrocytes, binds a small group of rat brain proteins. Predominant among these is heterogeneous nuclear ribonucleoprotein (hnRNP) A2, a trans-acting factor for cytoplasmic trafficking of RNAs bearing A2RE-like sequences. We have now identified the other A2RE-binding proteins as hnRNP A1/A1(B), hnRNP B1, and four isoforms of hnRNP A3. The rat and human hnRNP A3 cDNAs have been sequenced, revealing the existence of alternatively spliced mRNAs. In Western blotting, 38-, 39-, 41 -, and 41.5-kDa components were all recognized by antibodies against a peptide in the glycine-rich region of hnRNP A3, but only the 41- and 41.5-kDa bands bound antibodies to a 15-residue N-terminal peptide encoded by an alternatively spliced part of exon 1. The identities of these four proteins were verified by Edman sequencing and mass spectral analysis of tryptic fragments generated from electrophoretically separated bands. Sequence-specific binding of bacterially expressed hnRNP A3 to A2RE has been demonstrated by biosensor and UV cross-linking electrophoretic mobility shift assays. Mutational analysis and confocal microscopy data support the hypothesis that the hnRNP A3 isoforms have a role in cytoplasmic trafficking of RNA.

Sediment yield and impacts from river catchments to the Great Barrier Reef lagoon

Land use intensification is estimated to result in an overall increase in sediment delivery to the Great Barrier Reef lagoon by a factor of approximately four. Modelling suggests that, following land use intensification, croplands cause the greatest increase of sediment yield and sediment concentration, whereas erosion of grazing land is the main contemporary source of sediments, primarily owing to the large spatial extent of this land use. The spatial pattern of sediment yield to the coast after land use intensification is strongly correlated with the pattern under natural conditions, although the greatest increase is estimated to have occurred in the wet-dry catchments. Sediment transport and resuspension processes have led to the development of a strongly sediment-partitioned shelf, with modern mud-rich sediments almost exclusively restricted to the inner and inner-middle shelf, northward-facing embayments and in the lee of headlands. Elevated sediment concentrations increase the potential transport rates of nutrients and other pollutants. Whether increased sediment supply to the coastal zone has impacted on reefs remains a point of contention. More sediment load data need to be collected and analysed in order to make detailed estimates of catchment yields and establish the possible sediment impact on the Great Barrier Reef.

Adipose tissue as an endocrine organ

Adipose tissue is a highly active endocrine organ secreting a range of soluble products with both local and distant actions. These hormones have important roles in metabolism, reproduction, cardiovascular function and immunity. It is now evident that adipose endocrine function directly influences other organ systems, including the brain, liver and skeletal muscle. The endocrine function of adipose tissue is significantly regulated by nutritional status, and both are inextricably linked to the energy storage role of adipose tissue. This chapter highlights the endocrinology of adipose tissue by concentrating on functional aspects of the secreted products. The data of particular relevance to humans are highlighted, and areas in need of future research are suggested.

Shear stress and sediment transport calculations for swash zone modelling

It is shown that the observed difference in sediment transporting efficiency by the swash uprush, compared with the downrush, could be mainly due to greater bed shear stress for a given velocity in the more abruptly accelerated uprush. The bed shear stress generated by an arbitrary free stream velocity time series is modelled in terms of usual wave boundary layer models plus a phase lead (phi(tau) of the bed shear stress compared with the free stream velocity at the peak frequency. With this approach, the total transport amounts in uprush and downrush can be modelled satisfactorily with the same sediment transport formula, without the need for different uprush and downrush coefficients. While the adaptation of sediment transport formulae from steady flow can thus lead to the right total amounts of sediment moved by this method, the timing of the instantaneous sediment transport rates are probably not accurately modelled due to the highly unsteady nature of the swash and the presence of pre-suspended sediment in the uprush. Nevertheless, the proposed method is a useful intermediate step before we have a complete understanding of sediment transport under very rapid accelerations and of the relative contribution of pre-suspended sediment to the onshore sediment transport in swash zones. (C) 2002 Published by Elsevier Science B.V.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.

Vertical fluxes of sediment in oscillatory sheet flow

Time series of vertical sediment fluxes are derived from concentration time series in sheet flow under waves. While the concentrations C(z,t) vary very little with time for \z\ < 10d(50), the measured vertical sediment fluxes Q(zs)(z,t) vary strongly with time in this vertical band and their time variation follows, to some extent, the variation of the grain roughness Shields parameter 02,5(t). Thus, sediment distribution models based on the pickup function boundary condition are in some qualitative agreement with the measurements. However, the pickup function models are only able to model the upward bursts of sediment during the accelerating phases of the flow. They are, so far, unable to model the following strong downward sediment fluxes, which are observed during the periods of flow deceleration. Classical pickup functions, which essentially depend on the Shields parameter, are also incapable of modelling the secondary entrainment fluxes, which sometimes occur at free stream velocity reversal. The measured vertical fluxes indicate that the effective sediment settling velocity in the high [(0.3 < C(z,t) < 0.4] concentration area is typically only a few percent of the clear water settling velocity, while the measurements of Richardson and Jeronimo [Chem. Eng. Sci. 34 (1979) 1419], from a different physical setting, lead to estimates of the order 20%. The data does not support gradient diffusion as a model for sediment entrainment from the bed. That is, detailed modelling of the observed near-bed fluxes would require diffusivities that go negative during periods of flow deceleration. An observed general trend for concentration variability to increase with elevation close to the bed is also irreconcilable with diffusion models driven by a bottom boundary condition. (C) 2002 Published by Elsevier Science B.V.