Cluster-based network model for time-course gene expression data

We propose a model-based approach to unify clustering and network modeling using time-course gene expression data. Specifically, our approach uses a mixture model to cluster genes. Genes within the same cluster share a similar expression profile. The network is built over cluster-specific expression profiles using state-space models. We discuss the application of our model to simulated data as well as to time-course gene expression data arising from animal models on prostate cancer progression. The latter application shows that with a combined statistical/bioinformatics analyses, we are able to extract gene-to-gene relationships supported by the literature as well as new plausible relationships.

Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large

Conifers are resistant to attack from a large number of potential herbivores or pathogens. Previous molecular and biochemical characterization of selected conifer defence systems support a model of multigenic, constitutive and induced defences that act on invading insects via physical, chemical, biochemical or ecological (multitrophic) mechanisms. However, the genomic foundation of the complex defence and resistance mechanisms of conifers is largely unknown. As part of a genomics strategy to characterize inducible defences and possible resistance mechanisms of conifers against insect herbivory, we developed a cDNA microarray building upon a new spruce (Picea spp.) expressed sequence tag resource. This first-generation spruce cDNA microarray contains 9720 cDNA elements representing c. 5500 unique genes. We used this array to monitor gene expression in Sitka spruce (Picea sitchensis) bark in response to herbivory by white pine weevils (Pissodes strobi, Curculionidae) or wounding, and in young shoot tips in response to western spruce budworm (Choristoneura occidentalis, Lepidopterae) feeding. Weevils are stem-boring insects that feed on phloem, while budworms are foliage feeding larvae that consume needles and young shoot tips. Both insect species and wounding treatment caused substantial changes of the host plant transcriptome detected in each case by differential gene expression of several thousand array elements at 1 or 2 d after the onset of treatment. Overall, there was considerable overlap among differentially expressed gene sets from these three stress treatments. Functional classification of the induced transcripts revealed genes with roles in general plant defence, octadecanoid and ethylene signalling, transport, secondary metabolism, and transcriptional regulation. Several genes involved in primary metabolic processes such as photosynthesis were down-regulated upon insect feeding or wounding, fitting with the concept of dynamic resource allocation in plant defence. Refined expression analysis using gene-specific primers and real-time PCR for selected transcripts was in agreement with microarray results for most genes tested. This study provides the first large-scale survey of insect-induced defence transcripts in a gymnosperm and provides a platform for functional investigation of plant-insect interactions in spruce. Induction of spruce genes of octadecanoid and ethylene signalling, terpenoid biosynthesis, and phenolic secondary metabolism are discussed in more detail.

Construction, application and validation of selection evaluation model (SEM) for intelligent HVAC control system

Design teams are confronted with the quandary of choosing apposite building control systems to suit the needs of particular intelligent building projects, due to the availability of innumerable ‘intelligent’ building products and a dearth of inclusive evaluation tools. This paper is organised to develop a model for facilitating the selection evaluation for intelligent HVAC control systems for commercial intelligent buildings. To achieve these objectives, systematic research activities have been conducted to first develop, test and refine the general conceptual model using consecutive surveys; then, to convert the developed conceptual framework into a practical model; and, finally, to evaluate the effectiveness of the model by means of expert validation. The results of the surveys are that ‘total energy use’ is perceived as the top selection criterion, followed by the‘system reliability and stability’, ‘operating and maintenance costs’, and ‘control of indoor humidity and temperature’. This research not only presents a systematic and structured approach to evaluate candidate intelligent HVAC control system against the critical selection criteria (CSC), but it also suggests a benchmark for the selection of one control system candidate against another.

Selection-adaption-retention dynamics and variety in organisational routines

The notion of routines as mechanisms for achieving stability and change in organisations is well established in the organisational theory literature (Becker, 2004). However the relationship between the dynamics of selection, adaptation and retention and the increase or decrease in the varieties of routines which are the result of these processes, is not as well established theoretically or empirically. This paper investigates the processes associated with the evolution of an inter-organisational routine over time. The paper contributes to theory by advancing a conceptual clarification between the dynamics of organisational routines which produce variation, and the varieties of routines which are generated as a result of such processes; and an explanation for the relationship between selection, adaptation and retention dynamics and the creation of variety. The research is supported by analysis of empirical data pertaining to the procurement of engineering assets in a large asset intensive organisation.

Inhibitors to optimal project portfolio selection

The selection of projects and programs of work is a key function of both public and private sector organisations. Ideally, projects and programs that are selected to be undertaken are consistent with strategic objectives for the organisation; will provide value for money and return on investment; will be adequately resourced and prioritised; will not compete with general operations for resources and not restrict the ability of operations to provide income to the organisation; will match the capacity and capability of the organisation to deliver; and will produce outputs that are willingly accepted by end users and customers. Unfortunately,this is not always the case. Possible inhibitors to optimal project portfolio selection include: processes that are inconsistent with the needs of the organisation; reluctance to use an approach that may not produce predetermined preferences; loss of control and perceived decision making power; reliance on quantitative methods rather than qualitative methods for justification; ineffective project and program sponsorship; unclear project governance, processes and linkage to business strategies; ignorance, taboos and perceived effectiveness; inadequate education and training about the processes and their importance.

Nuclear localization of NPR1 is required for activation of PR gene expression

Systemic acquired resistance (SAR) is a broad-spectrum resistance in plants that involves the upregulation of a battery of pathogenesis-related (PR) genes. NPR1 is a key regulator in the signal transduction pathway that leads to SAR. Mutations in NPR1 result in a failure to induce PR genes in systemic tissues and a heightened susceptibility to pathogen infection, whereas overexpression of the NPR1 protein leads to increased induction of the PR genes and enhanced disease resistance. We analyzed the subcellular localization of NPR1 to gain insight into the mechanism by which this protein regulates SAR. An NPR1–green fluorescent protein fusion protein, which functions the same as the endogenous NPR1 protein, was shown to accumulate in the nucleus in response to activators of SAR. To control the nuclear transport of NPR1, we made a fusion of NPR1 with the glucocorticoid receptor hormone binding domain. Using this steroid-inducible system, we clearly demonstrate that nuclear localization of NPR1 is essential for its activity in inducing PR genes.

Interaction of NPR1 with basic leucine zipper protein transcription factors that bind sequences required for salicylic acid induction of the PR–1 gene

The Arabidopsis thaliana NPR1 has been shown to be a key regulator of gene expression during the onset of a plant disease-resistance response known as systemic acquired resistance. The npr1 mutant plants fail to respond to systemic acquired resistance-inducing signals such as salicylic acid (SA), or express SA-induced pathogenesis-related (PR) genes. Using NPR1 as bait in a yeast two-hybrid screen, we identified a subclass of transcription factors in the basic leucine zipper protein family (AHBP-1b and TGA6) and showed that they interact specifically in yeast and in vitro with NPR1. Point mutations that abolish the NPR1 function in A. thaliana also impair the interactions between NPR1 and the transcription factors in the yeast two-hybrid assay. Furthermore, a gel mobility shift assay showed that the purified transcription factor protein, AHBP-1b, binds specifically to an SA-responsive promoter element of the A. thaliana PR-1 gene. These data suggest that NPR1 may regulate PR-1 gene expression by interacting with a subclass of basic leucine zipper protein transcription factors.

Molecular analysis of the myosin gene family in Arabidopsis thaliana

Myosin is believed to act as the molecular motor for many actin-based motility processes in eukaryotes. It is becoming apparent that a single species may possess multiple myosin isoforms, and at least seven distinct classes of myosin have been identified from studies of animals, fungi, and protozoans. The complexity of the myosin heavy-chain gene family in higher plants was investigated by isolating and characterizing myosin genomic and cDNA clones from Arabidopsis thaliana. Six myosin-like genes were identified from three polymerase chain reaction (PCR) products (PCR1, PCR11, PCR43) and three cDNA clones (ATM2, MYA2, MYA3). Sequence comparisons of the deduced head domains suggest that these myosins are members of two major classes. Analysis of the overall structure of the ATM2 and MYA2 myosins shows that they are similar to the previously-identified ATM1 and MYA1 myosins, respectively. The MYA3 appears to possess a novel tail domain, with five IQ repeats, a six-member imperfect repeat, and a segment of unique sequence. Northern blot analyses indicate that some of the Arabidopsis myosin genes are preferentially expressed in different plant organs. Combined with previous studies, these results show that the Arabidopsis genome contains at least eight myosin-like genes representing two distinct classes.

Genetic resources in wild and cultured stocks of the Asian mud crab, Scylla paramamosain

The mud crab (Scylla spp.) aquaculture industry has expanded rapidly in recent years in many countries in the Indo - West Pacific (IWP) region as an alternative to marine shrimp culture because of significant disease outbreaks and associated failures of many shrimp culture industries in the region. Currently, practices used to produce and manage breeding crabs in hatcheries may compromise levels of genetic diversity, ultimately compromising growth rates, disease resistance and stock productivity. Therefore, to avoid “genetic pollution” and its harmful effects and to promote further development of mud crab aquaculture and fisheries in a sustainable way, a greater understanding of the genetic attributes of wild and cultured mud crab stocks is required. Application of these results can provide benefits for managing wild and cultured Asian mud crab populations for multiple purposes including for commercial production, recreation and conservation and to increase profitability and sustainability of newly emerging crab culture industries. Phylogeographic patterns and the genetic structure of Asian mud crab populations across the IWP were assessed to determine if they were concordant with those of other widespread taxa possessing pelagic larvae of relatively long duration. A 597 bp fragment of the mitochondrial DNA COI gene was amplified and screened for variation in a total of 297 individuals of S. paramamosain from six sampling sites across the species’ natural geographical distribution in the IWP and 36 unique haplotypes were identified. Haplotype diversities per site ranged from 0.516 to 0.879. Nucleotide diversity estimates among haplotypes were 0.11% – 0.48%. Maximum divergence observed among S. paramamosain samples was 1.533% and samples formed essentially a single monophyletic group as no obvious clades were related to geographical location of sites. A weak positive relationship was observed however, between genetic distance and geographical distance among sites. Microsatellite markers were then used to assess contemporary gene flow and population structure in Asian mud crab populations sampled across their natural distribution in the IWP. Eight microsatellite loci were screened in sampled S. paramamosain populations and all showed high allelic diversity at all loci in sampled populations. In total, 344 individuals were analysed, and 304 microsatellite alleles were found across the 8 loci. The mean number of alleles per locus at each site ranged from 20.75 to 28.25. Mean allelic richness per site varied from 17.2 to 18.9. All sites showed high levels of heterozygosity as average expected heterozygosities for all loci ranged from 0.917 – 0.953 while mean observed heterozygosity ranged from 0.916 – 0.959. Allele diversities were similar at all sites and across all loci. The results did not show any evidence for major differences in allele frequencies among sites and patterns of allele frequencies were very similar in all populations across all loci. Estimates of population differentiation (FST) were relatively low and most probably largely reflect intra – individual variation for very highly variable loci. Results from nDNA analysis showed evidence for only very limited population genetic structure among sampled S. paramamosain, and a positive and significant association for genetic and geographical distance among sample sites. Microsatellite markers were then employed to determine if adequate levels of genetic diversity has been captured in crab hatcheries for the breeding cycle. The results showed that all microsatellite loci were polymorphic in hatchery samples. Culture populations were in general, highly genetically depauperate, compared with comparable wild populations, with only 3 to 8 alleles recorded for the same loci set per population. In contrast, very high numbers of alleles per locus were found in reference wild S. paramamosain populations, which ranged from 18 to 46 alleles per locus per population. In general, this translates into a 3 to 10 fold decline in mean allelic richness per locus in all culture stocks compared with wild reference counterparts. Furthermore, most loci in all cultured S. paramamosain samples showed departures from HWE equilibrium. Allele frequencies were very different in culture samples from that present in comparable wild reference samples and this in particular, was reflected in a large decline in allele diversity per locus. The pattern observed was best explained by significant impacts of breeding practices employed in hatcheries rather than natural differentiation among wild populations used as the source of brood stock. Recognition of current problems and management strategies for the species both for the medium and long-term development of the new culture industry are discussed. The priority research to be undertaken over the medium term for S. paramamosain should be to close the life cycle fully to allow individuals to be bred on demand and their offspring equalised to control broodstock reproductive contributions. Establishing a broodstock register and pedigree mating system will be required before any selection program is implemented. This will ensure that sufficient genetic variation will be available to allow genetic gains to be sustainably achieved in a future stock improvement program. A fundamental starting point to improve hatchery practices will be to encourage farmers and hatchery managers to spawn more females in their hatcheries as it will increase background genetic diversity in culture stocks. Combining crablet cohorts from multiple hatcheries into a single cohort for supply to farmers or rotation of breeding females regularly in hatcheries will help to address immediate genetic diversity problems in culture stocks. Application of these results can provide benefits for managing wild and cultured Asian mud crab populations more efficiently. Over the long-term, application of data on genetic diversity in wild and cultured stocks of Asian mud crab will contribute to development of sustainable and productive culture industries in Vietnam and other countries in the IWP and can contribute towards conservation of wild genetic resources.

EduBuilding: material selection from life cycle costing sensitivity

Life Cycle Cost Analysis provides a form of synopsis of the initial and consequential costs of building related decisions. These cost figures may be implemented to justify higher investments, for example, in the quality or flexibility of building solutions through a long term cost reduction. The emerging discipline of asset mnagement is a promising approach to this problem, because it can do things that techniques such as balanced scorecards and total quantity cannot. Decisions must be made about operating and maintaining infrastructure assets. An injudicious sensitivity of life cycle costing is that the longer something lasts, the less it costs over time. A life cycle cost analysis will be used as an economic evaluation tool and collaborate with various numbers of analyses. LCCA quantifies incurring costs commonly overlooked (by property and asset managers and designs) as replacement and maintenance costs. The purpose of this research is to examine the Life Cycle Cost Analysis on building floor materials. By implementing the life cycle cost analysis, the true cost of each material will be computed projecting 60 years as the building service life and 5.4% as the inflation rate percentage to classify and appreciate the different among the materials. The analysis results showed the high impact in selecting the floor materials according to the potential of service life cycle cost next.

Improving the performance of model-order selection criteria by partial-model selection search

The traditional searching method for model-order selection in linear regression is a nested full-parameters-set searching procedure over the desired orders, which we call full-model order selection. On the other hand, a method for model-selection searches for the best sub-model within each order. In this paper, we propose using the model-selection searching method for model-order selection, which we call partial-model order selection. We show by simulations that the proposed searching method gives better accuracies than the traditional one, especially for low signal-to-noise ratios over a wide range of model-order selection criteria (both information theoretic based and bootstrap-based). Also, we show that for some models the performance of the bootstrap-based criterion improves significantly by using the proposed partial-model selection searching method. Index Terms— Model order estimation, model selection, information theoretic criteria, bootstrap 1. INTRODUCTION Several model-order selection criteria can be applied to find the optimal order. Some of the more commonly used information theoretic-based procedures include Akaike’s information criterion (AIC) [1], corrected Akaike (AICc) [2], minimum description length (MDL) [3], normalized maximum likelihood (NML) [4], Hannan-Quinn criterion (HQC) [5], conditional model-order estimation (CME) [6], and the efficient detection criterion (EDC) [7]. From a practical point of view, it is difficult to decide which model order selection criterion to use. Many of them perform reasonably well when the signal-to-noise ratio (SNR) is high. The discrepancies in their performance, however, become more evident when the SNR is low. In those situations, the performance of the given technique is not only determined by the model structure (say a polynomial trend versus a Fourier series) but, more importantly, by the relative values of the parameters within the model. This makes the comparison between the model-order selection algorithms difficult as within the same model with a given order one could find an example for which one of the methods performs favourably well or fails [6, 8]. Our aim is to improve the performance of the model order selection criteria in cases where the SNR is low by considering a model-selection searching procedure that takes into account not only the full-model order search but also a partial model order search within the given model order. Understandably, the improvement in the performance of the model order estimation is at the expense of additional computational complexity.

Genetic Diversity of Rice Tungro Spherical Virus in Tungro-Endemic Provinces of the Philippines and Indonesia

The two adjacent genes of coat protein 1 and 2 of rice tungro spherical virus (RTSV) were amplified from total RNA extracts of serologically indistinguishable field isolates from the Philippines and Indonesia, using reverse transcriptase polymerase chain reaction (RT-PCR). Digestion with HindIII and BstYI restriction endonucleases differentiated the amplified DNA products into eight distinct coat protein genotypes. These genotypes were then used as indicators of virus diversity in the field. Inter- and intra-site diversities were determined over three cropping seasons. At each of the sites surveyed, one or two main genotypes prevailed together with other related minor or mixed genotypes that did not replace the main genotype over the sampling time. The cluster of genotypes found at the Philippines sites was significantly different from the one at the Indonesia sites, suggesting geographic isolation for virus populations. Phylogenetic studies based on the nucleotide sequences of 38 selected isolates confirm the spatial distribution of RTSV virus populations but show that gene flow may occur between populations. Under the present conditions, rice varieties do not seem to exert selective pressure on the virus populations. Based on the selective constraints in the coat protein amino acid sequences and the virus genetic composition per site, a negative selection model followed by random-sampling events due to vector transmissions is proposed to explain the inter-site diversity observed