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.

Steiner trades that give rise to completely decomposable latin interchanges

In this paper we focus on the representation of Steiner trades of volume less than or equal to nine and identify those for which the associated partial latin square can be decomposed into six disjoint latin interchanges.

Pea rms6 mutants exhibit increased basal branching

Our studies on two branching mutants of pea (Pisum sativum L.) have identified a further Ramosus locus, Rms6, with two recessive or partially recessive mutant alleles: rms6-1 (type line S2-271) and rms6-2 (type line K586). Mutants rms6-1 and rms6-2 were derived from dwarf and tall cultivars, Solara and Torsdag, respectively. The rms6 mutants are characterized by increased branching from basal nodes. In contrast, mutants rms1 through rms5 have increased branching from both basal and aerial (upper stem) nodes. Buds at the cotyledonary node of wild-type (WT) plants remain dormant but in rms6 plants these buds were usually released from dormancy. Their growth was either subsequently inhibited, sometimes even prior to emergence above ground, or they grew into secondary stems. The mutant phenotype was strongest for rms6-1 on the dwarf background. Although rms6-2 had a weak single-mutant phenotype, the rms3-1 rms6-2 double mutant showed clear transgression and an additive branching phenotype, with a total lateral length almost 2-fold greater than rms3-1 and nearly 5-fold greater than rms6-2 . Grafting studies between WT and rms6-1 plants demonstrated the primary action of Rms6 may be confined to the shoot. Young WT and rms6-1 shoots had similar auxin levels, and decapitated plants had a similar magnitude of response to applied auxin. Abscisic acid levels were elevated 2-fold at node 2 of young rms6-1 plants. The Rms6 locus mapped to the R to Gp segment of linkage group V (chromosome 3). The rms6 mutants will be useful for basic research and also have possible agronomical value.

Callus induction and plant regeneration from broccoli (Brassica oleracea var. italica) for transformation

We evaluated the efficiency of callus induction and plantlet regeneration from hypocotyl explants of broccoli (Brassica oleracea var. italica). The cultivars were ‘Marathon’, ‘Greenbelt’, and ‘Shogun’. Transformation success was not affected by the presence of tobacco feeder-cell layers on the culture media. The frequency of shoot regeneration was greater from 10-d-old hypocotyls than from 14-d-old hypocotyls. Both ‘Marathon’ and ‘Greenbelt’ had higher potentials for tissue regeneration than did ‘Shogun’. We found that for transformation selection, the optimum concentration was either 50 mg/L kanamycin or 100 mg/L genetkin.

Dry matter production and boron concentrations of vegetative and reproductive tissues of canola and sunflower plants grown in nutrient solution

Canola (Brassica napus L.) and sunflower (Helianthus annuus L.), two important oilseed crops, are sensitive to low boron (B) supply. Symptoms of B deficiency are often more severe during the reproductive stage, but it is not known if this is due to a decreased external B supply with time or an increased sensitivity to low B during this stage. Canola and sunflower were grown for 75 days after transplanting (DAT) in two solution culture experiments using Amberlite (IRA-743) B-specific resin to maintain constant B concentration in solution over the range 0.6 - 53 muM. Initially, the vegetative growth of both crops was good in all treatments. With the onset of the reproductive stage, however, severe B deficiency symptoms developed and growth of canola and sunflower was reduced with less than or equal to 0.9 and less than or equal to 0.7 muM B, respectively. At these concentrations, reproductive parts failed to develop. The critical B concentration (i.e. 90% of maximum shoot dry matter yield) in the youngest opened leaf was 18 mg kg(-1) in canola and 25 mg kg(-1) in sunflower at 75 DAT. The results of this study indicate that the reproductive stage of these two oilseed crops is more sensitive than the vegetative stage to low B supply.

Simulating growth, development, and yield of tillering pearl millet. III. Biomass accumulation and partitioning

Pearl millet landraces from Rajasthan, India, yield significantly less than improved cultivars under optimum growing conditions, but not under stressed conditions. To successfully develop a simulation model for pearl millet, capable of capturing such genotype x environment (G x E) interactions for grain yield, we need to understand the causes of the observed yield interaction. The aim of this paper is to quantify the key parameters that determine the accumulation and partitioning of biomass: the,light extinction coefficient, radiation use efficiency (RUE), pattern of dry matter allocation to the leaf blades, the determination of grain number, and the rate and duration of dry matter accumulation into individual grains. We used data on improved cultivars and landraces, obtained from both published and unpublished sources collected at ICRISAT, Patancheru, India. Where possible, the effects of cultivar and axis (main shoot vs. tillers) on these parameters were analysed, as previous research suggested that G x E interactions for grain yield are associated with differences in tillering habit. Our results indicated there were no cultivar differences in extinction coefficient, RUE, and biomass partitioning before anthesis, and differences between axes in biomass partitioning were negligible. This indicates there was no basis for cultivar differences in the potential grain yield. Landraces, however, produced consistently less grain yield for a given rate of dry matter accumulation at anthesis than did improved cultivars. This was caused by a combination of low grain number and small grain size. The latter was predominantly due to a lower grain growth rate, as genotypic differences in the duration of grain filling were relatively small. Main shoot and tillers also had a similar duration of grain filling. The low grain yield of the landraces was associated with profuse nodal tillering, supporting the hypothesis that grain yield was below the potential yield that could be supported by assimilate availability. We hypothesise this is a survival strategy, which enhances the prospects to escape the effects of stress around anthesis. (C) 2002 E.J. van Oosterom. Published by Elsevier Science B.V. All rights reserved.

Yield response of rice (Oryza sativa L.) genotypes to different types of drought under rainfed lowlands - Part 3. Plant factors contributing to drought resistance

A series of experiments were conducted in drought-prone northeast Thailand to examine the magnitude of yield responses of diverse genotypes to drought stress environments and to identify traits that may confer drought resistance to rainfed lowland rice. One hundred and twenty eight genotypes were grown under non-stress and four different types of drought stress conditions. Under severe drought conditions, the maintenance of PWP of genotypes played a significant role in determining final grain yield. Because of their smaller plant size (lower total dry matter at anthesis) genotypes that extracted less soil water during the early stages of the drought period, tended to maintain higher PWP and had a higher fertile panicle percentage, filled grain percentage and final grain yield than other genotypes. PWP was correlated with delay in flowering (r = -0.387) indicating that the latter could be used as a measure of water potential under stress. Genotypes with well-developed root systems extracted water too rapidly and experienced severe water stress at flowering. RPR which showed smaller coefficient of variation was more useful than root mass density in identifying genotypes with large root system. Under less severe and prolonged drought conditions, genotypes that could achieve higher plant dry matter at anthesis were desirable. They had less delay in flowering, higher grain yield and higher drought response index, indicating the importance of ability to grow during the prolonged stress period. Other shoot characters (osmotic potential, leaf temperature, leaf rolling, leaf death) had little effect on grain yield under different drought conditions. This was associated with a lack of genetic variation and difficulty in estimating trait values precisely. Under mild stress conditions (yield loss less than 50%), there was no significant relationship between the measured drought characters and grain yield. Under these mild drought conditions, yield is determined more by yield potential and phenotype than by drought resistant mechanisms per se. (C) 2002 Elsevier Science B.V. All rights reserved.

Growing Atriplex and Maireana species in saline sodic and waterlogged soils

Utilization of salt affected wasteland by growing forage shrubs has enormous economic and environmental implication for developing countries like Pakistan, where approximately 6.3 million ha of the land is salt affected. Considering the importance of Atriplex and Maireana species, research has been conducted using their different species on the salt affected soils of Faisalabad. Most of Atriplex and Maireana species survived under the environmental conditions of Faisalabad and gave the good yield in the form of forage. Some of these species are woody and can be used for fuel purposes. Sixteen genotypes of Atriplex and Maireana were tested for their tolerance to waterlogging in order to identify halophytic fodder shrubs suitable for growth on secondary salt-affected and waterlogged farmland. The physiological and morphological responses of the species tested were typical of species with a generally poor tolerance to waterlogging. Despite this, some species (e.g., A. amnicola) were surprisingly resistant, surviving up to five months of waterlogging at moderate salinity and high evapotranspirational demand. The most resistant species, A amnicola maintained higher transpiration rates, leaf water potentials and shoot extension rates than most other species during five weeks of waterlogging, and a return to control levels more quickly than other species after plots were drained. Although little morphological adaptation to waterlogged conditions was detected, a shallow and extensive lateral root system and the formation of many short aerenchymatous adventitious roots from procumbent branches appeared to advantage A. amnicola in an environment highly heterogeneous in salinity and low in oxygen concentration. Waterlogging quickly killed shallow fibrous rooted species, although the procumbent branches of some individuals survived as clones if they developed adventitious roots.

Fast neutron mutagenesis of soybean (Glycine soja L.) produces a supernodulating mutant containing a large deletion in linkage group H

We describe for the first time the application of fast neutron mutagenesis to the genetic dissection of root nodulation in legumes. We demonstrate the utility of chromosomal deletion mutations through production of a soybean supernodulation mutant FN37 that lacks the internal autoregulation of nodulation mechanism. After inoculation with microsymbiont Bradyrhizobium japonicum, FN37 forms at least 10 times more nodules than the wild type G. soja parent and has a phenotype identical to that of chemically induced allelic mutants nts382 and nts1007 (NTS-1 locus). Reciprocal grafting of shoots and roots confirmed systemic shoot control of the FN37 nodulation phenotype. RFLP/PCR marker pUTG132a and AFLP marker UQC-IS1 which are tightly linked to NTS-1 allowed the isolation of BAC contigs delineating both ends of the deletion. The genetic/physical distance ratio in the NTS-1 region is 279 kb/cM. The deletion is estimated to be about 460 kb based on the absence of markers and bacterial artificial chromosomes (BAC) ends as well as genetic and physical mapping. Deletion break points were determined physically and placed within flanking BAC contigs.

Identification of research to improve the efficiency of breeding strategies for white clover in Australia - a review

A major challenge faced by today's white clover breeder is how to manage resources within a breeding program. It is essential to utilise these resources with sufficient flexibility to build on past progress from conventional breeding strategies, but also take advantage of emerging opportunities from molecular breeding tools such as molecular markers and transformation. It is timely to review white clover breeding strategies. This background can then be used as a foundation for considering how to continue conventional plant improvement activities and complement them with molecular breeding opportunities. In this review, conventional white clover breeding strategies relevant to the Australian dryland target population environments are considered. Attention is given to: (i) availability of genetic variation, (ii) characterisation of germplasm collections, (iii) quantitative models for estimation of heritability, (iv) the role of multi-environment trials to accommodate genotype-by-environment interactions, (v) interdisciplinary research to understand adaptation to dryland environments, (vi) breeding and selection strategies, and (vii) cultivar structure. Current achievements in biotechnology with specific reference to white clover breeding in Australia are considered, and computer modelling of breeding programs is discussed as a useful integrative tool for the joint evaluation of conventional and molecular breeding strategies and optimisation of resource use in breeding programs. Four areas are identified as future research priorities: (i) capturing the potential genetic diversity among introduced accessions and ecotypes that are adapted to key constraints such as summer moisture stress and the use of molecular markers to assess the genetic diversity, (ii) understanding the underlying physiological/morphological root and shoot mechanisms involved in water use efficiency of white clover, with the objective of identifying appropriate selection criteria, (iii) estimation of quantitative genetic parameters of important morphological/physiological attributes to enable prediction of response to selection in target environments, and (iv) modelling white clover breeding strategies to evaluate the opportunities for integration of molecular breeding strategies with conventional breeding programs.

Modelling High-Dimensional Data by Mixtures of Factor Analyzers

We focus on mixtures of factor analyzers from the perspective of a method for model-based density estimation from high-dimensional data, and hence for the clustering of such data. This approach enables a normal mixture model to be fitted to a sample of n data points of dimension p, where p is large relative to n. The number of free parameters is controlled through the dimension of the latent factor space. By working in this reduced space, it allows a model for each component-covariance matrix with complexity lying between that of the isotropic and full covariance structure models. We shall illustrate the use of mixtures of factor analyzers in a practical example that considers the clustering of cell lines on the basis of gene expressions from microarray experiments. (C) 2002 Elsevier Science B.V. All rights reserved.

Spatial and temporal changes in multiple hormone groups during lateral bud release shortly following apex decapitation of chickpea (Cicer arietinum) seedlings

Although the co-ordination of promotive root-sourced cytokinin (CK) and inhibitory shoot apex-sourced auxin (IAA) is central to all current models on lateral bud dormancy release, control by those hormones alone has appeared inadequate in many studies. Thus it was hypothesized that the IAA : CK model is the central control but that it must be considered within the relevant timeframe leading to lateral bud release and against a backdrop of interactions with other hormone groups. Therefore, IAA and a wide survey of cytokinins (CKs), were examined along with abscisic acid (ABA) and polyamines (PAs) in released buds, tissue surrounding buds and xylem sap at 1 and 4 h after apex removal, when lateral buds of chickpea are known to break dormancy. Three potential lateral bud growth inhibitors, IAA, ABA and cis-zeatin 9-riboside (ZR), declined sharply in the released buds and xylem following decapitation. This is in contrast to potential dormancy breaking CKs like trans-ZR and trans-zeantin 9-riboside 5'phosphate (ZRMP), which represented the strongest correlative changes by increasing 3.5-fold in xylem sap and 22-fold in buds. PAs had not changed significantly in buds or other tissues after 4 h, so they were not directly involved in the breaking of bud dormancy. Results from the xylem and surrounding tissues indicated that bud CK increases resulted from a combination synthesis in the bud and selective loading of CK nucleotides into the xylem from the root.

Transport and metabolism of xylem cytokinins during lateral bud release in decapitated chickpea (Cicer arietinum) seedlings

Although cytokinins (CKs) are widely thought to have a role in promoting shoot branching, there is little data supporting a causative or even a correlative relationship between endogenous CKs and timing of bud outgrowth. We previously showed that lateral bud CK content increased rapidly following shoot decapitation. However, it is not known whether roots are the source of this CK. Here, we have used shoot decapitation to instantaneously induce lateral bud release in chickpea seedlings. This treatment rapidly alters rate and direction of solvent and solute (including CK) trafficking, which may be a passive signalling mechanism central to initiation of lateral bud release. To evaluate changes in xylem transport, intact and decapitated plants were infiltrated with [H-3]zeatin riboside ([H-3]ZR), a water-soluble blue dye or [H-3]H2O by injection into the hypocotyl. All three tracers were recovered in virtually all parts of the shoot within I h of injection. In intact plants, solute accumulation in the lateral bud at node 1 was significantly less than in the adjacent stipule and nodal tissue. In decapitated plants, accumulation of [H-3]ZR and of blue dye in the same bud position was increased 3- to 10-fold relative to intact plants, whereas content of [H-3]H2O was greatly reduced indicating an increased solvent throughput. The stipule and cut stem, predicted to have high evapotranspiration rates, also showed increased solute content accompanied by enhanced depletion of [H-3]H2O. To assess whether metabolism modifies quantities of active CK reaching the buds, we followed the metabolic fate of [H-3]ZR injected at physiological concentrations. Within 1 h, 80-95% of [H-3]ZR was converted to other active CKs (mainly zeatin riboside-5'phosphate (ZRMP) and zeatin (Z)), other significant, but unconfirmed metabolites some of which may be active (O-acetylZR, O-acetylZRMP and a compound correlated with sites of high CK-concentrations) and inactive catabolites (adenosine, adenine, 5'AMP and water). Despite rapid metabolic degradation, the total active label, which was indicative of CK concentration in buds, increased rapidly following decapitation. It can be inferred that xylem sap CKs represent one source of active CKs appearing in lateral buds after shoot decapitation.

Do increases in xylem sap pH and/or ABA concentration mediate stomatal closure following nitrate deprivation?

Stomatal conductance (g(s)) of pepper (Capsicum annuum L.) plants decreased during the second photoperiod (day 2) after withholding nitrate (N). Stomatal closure of N-deprived plants was not associated with a decreased shoot water potential (Psi(shoot)); conversely Psi(shoot) was lower in N-supplied plants. N deprivation transiently (days 2 and 3) alkalized (0.2-0.3 pH units) xylem sap exuded from de-topped root systems under root pressure, and xylem sap expressed from excised shoots by pressurization. The ABA concentration of expressed sap increased 3-4-fold when measured on days 2 and 4. On day 2, leaves detached from N-deprived and N-supplied plants showed decreased transpiration rates when fed an alkaline (pH 7) artificial xylem (AX) solution, independent of the ABA concentration (10-100 nM) supplied. Thus changes in xylem sap composition following N deprivation can potentially close stomata. However, the lower transpiration rate of detached N-deprived leaves relative to N-supplied leaves shows that factors residing within N-deprived leaves also mediate stomatal closure, and that these factors assume greater importance as the duration of N deprivation increases.

Leaf area development of ABA-deficient and wild-type peas at two levels of nitrogen supply

The ABA-deficient wilty pea (Pisum sativum L.) and its wild-type (WT) were grown at two levels of nitrogen supply (0.5 and 5.0 mM) for 5-6 weeks from sowing, to determine whether leaf ABA status altered the leaf growth response to N deprivation. Plants were grown at high relative humidity to prevent wilting of the wilty peas. Irrespective of N supply, expanding wilty leaflets had ca 50% less ABA than WT leaflets but similar ethylene evolution rates. Fully expanded wilty leaflets had lower relative water contents (RWC) and were 10-60% smaller in area (according to the node of measurement) than WT leaflets. However, there were no genotypic differences in plant relative leaf expansion rate (RLER). Growth of both genotypes at 0.5 mM N increased the RWC of fully expanded leaflets, but did not alter ethylene evolution or ABA concentration of expanding leaflets. Plants grown at 0.5 mM N showed a 20-30% reduction in RLER, which was similar in magnitude in both wilty and WT peas. Thus, leaf ABA status did not alter the leaf growth response to N deprivation.