Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Levene tests of homogeneity of variance for general block and treatment designs

This article develops a weighted least squares version of Levene's test of homogeneity of variance for a general design, available both for univariate and multivariate situations. When the design is balanced, the univariate and two common multivariate test statistics turn out to be proportional to the corresponding ordinary least squares test statistics obtained from an analysis of variance of the absolute values of the standardized mean-based residuals from the original analysis of the data. The constant of proportionality is simply a design-dependent multiplier (which does not necessarily tend to unity). Explicit results are presented for randomized block and Latin square designs and are illustrated for factorial treatment designs and split-plot experiments. The distribution of the univariate test statistic is close to a standard F-distribution, although it can be slightly underdispersed. For a complex design, the test assesses homogeneity of variance across blocks, treatments, or treatment factors and offers an objective interpretation of residual plot.

Finite, three-dimensional adsorbed phase growth in activated carbon mesopores

The role of PACs (primary adsorption centers) in the mesopore (i.e., transport) region of activated carbons during adsorption of polar species, such as water, is unclear. A classical model of three-dimensional adsorption on finite PACs is presented. The model is a preliminary, theoretical investigation into adsorption on mesopore PACs and is intended to give some insight into the energetic and physical processes at work. Work processes are developed to obtain isotherms and three-dimensional sorbate growth on PACs of varying size and energetic characteristics. The work processes allow two forms of adsorbed phase growth: densification at constant boundary and boundary growth at constant density. Relatively strong sorbate-sorbent interactions and strong surface tension favor adsorbed phase densification over boundary growth. Conversely, relatively weak sorbate-sorbent interactions and weak surface tension favor boundary growth over densification. If sorbate-sorbate interactions are strong compared to sorbate-sorbent interactions, condensation with hysteresis occurs. This can also give rise to delayed boundary growth, where all initial adsorption occurs in the monolayer only. The results indicate that adsorbed phase growth on PACs may be quite complex.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Cadherin-directed actin assembly: E-cadherin physically associates with the Arp2/3 complex to direct actin assembly in nascent adhesive contacts

Cadherin cell adhesion molecules are major determinants of tissue patterning which function in cooperation with the actin cytoskeleton [1-4]. In the context of stable adhesion [1], cadherin/catenin complexes are often envisaged to passively scaffold onto cortical actin filaments. However, cadherins also form dynamic adhesive contacts during wound healing and morphogenesis [2]. Here actin polymerization has been proposed to drive cell surfaces together [5], although F-actin reorganization also occurs as cell contacts mature [6]. The interaction between cadherins and actin is therefore likely to depend on the functional state of adhesion. We sought to analyze the relationship between cadherin homophilic binding and cytoskeletal activity during early cadherin adhesive contacts. Dissecting the specific effect of cadherin ligation alone on actin regulation is difficult in native cell-cell contacts, due to the range of juxtacrine signals that can arise when two cell surfaces adhere [7]. We therefore activated homophilic ligation using a specific functional recombinant protein. We report the first evidence that E-cadherin associates with the Arp2/3 complex actin nucleator and demonstrate that cadherin binding can exert an active, instructive influence on cells to mark sites for actin assembly at the cell surface.

Exploring complex fitness surfaces: Multiple ornamentation and polymorphism in male guppies

The sexual ornamentation used by male guppies to attract females comprises many components, each of which varies considerably among males. Although natural and sexual selection have been shown to contribute to divergence among populations in male sexual ornaments, the role of sexual selection in maintaining polymorphism within populations is less clear. We used both parametric quadratic regression and nonparametric projection pursuit regression techniques to reveal the major axes of non-linear sexual selection on male ornaments. We visualized the fitness surfaces defined by these axes using thin-plate splines to allow a direct comparison of the two methodologies. Identification of the major axes of selection and their visualization was critical in determining the form and strength of nonlinear selection. Both types of analysis revealed fitness surfaces comprising three peaks, suggesting that there is more than one way to make an attractive guppy. Disruptive selection may be an important process underlying the presence of multiple sexual ornaments and may contribute to the maintenance of the high levels of polymorphism in male sexual ornaments found in guppy populations.

A structural basis for the selection of dominant alpha beta T cell receptors in antiviral immunity

We have examined the basis for immunodominant or public TCR usage in an antiviral CTL response. Residues encoded by each of the highly selected genetic elements of an immunodominant clonotype recognizing Epstein-Barr virus were critical to the antigen specificity of the receptor. Upon recognizing antigen the immunodominant TCR undergoes extensive conformational changes in the complementarity determining regions (CDRs), including the disruption of the canonical structures of the germline-encoded CDR1alpha and CDR2alpha loops to produce an enhanced fit with the HLA-peptide complex. TCR ligation induces conformational changes in the TCRalpha constant domain thought to form part of the docking site for CD3epsilon. These findings indicate that TCR immunodominance is associated with structural properties conferring receptor specificity and suggest a novel structural link between TCR ligation and intracellular signaling.

Speciation and phylogeography in Caridina indistincta, a complex of freshwater shrimps from Australian heathland streams

The Caridina indistincta complex is a group of closely related atyid shrimps that inhabit coastal freshwater streams throughout north-eastern Australia. Using mitochondrial DNA sequence data (cytochrome oxidase 1, CO1), we (1) inferred the timing of speciation in the C. indistincta group and (2) examined the intraspecific phylogeographic patterns within the group. Assuming a shrimp-specific rate of CO1 evolution, the level of sequence divergence among species suggests that speciation took place during the Miocene epoch. Within one widespread mainland species, phylogeographic patterns suggest strong geographic 'regionalisation' of mtDNA lineages that are most likely of Pleistocene origin. By contrast, another species comprises two highly divergent mtDNA lineages that occur in sympatry. We suggest that although Pleistocene sea-level regressions appear important in generating population-level phylogeographic patterns, these events were largely unimportant in the formation of species in this group.

The Lake Tekapo experiment (LTEX): An investigation of atmospheric boundary layer processes in complex terrain

A research program on atmospheric boundary layer processes and local wind regimes in complex terrain was conducted in the vicinity of Lake Tekapo in the southern Alps of New Zealand, during two 1-month field campaigns in 1997 and 1999. The effects of the interaction of thermal and dynamic forcing were of specific interest, with a particular focus on the interaction of thermal forcing of differing scales. The rationale and objectives of the field and modeling program are described, along with the methodology used to achieve them. Specific research aims include improved knowledge of the role of surface forcing associated with varying energy balances across heterogeneous terrain, thermal influences on boundary layer and local wind development, and dynamic influences of the terrain through channeling effects. Data were collected using a network of surface meteorological and energy balance stations, radiosonde and pilot balloon soundings, tethered balloon and kite-based systems, sodar, and an instrumented light aircraft. These data are being used to investigate the energetics of surface heat fluxes, the effects of localized heating/cooling and advective processes on atmospheric boundary layer development, and dynamic channeling. A complementary program of numerical modeling includes application of the Regional Atmospheric Modeling System (RAMS) to case studies characterizing typical boundary layer structures and airflow patterns observed around Lake Tekapo. Some initial results derived from the special observation periods are used to illustrate progress made to date. In spite of the difficulties involved in obtaining good data and undertaking modeling experiments in such complex terrain, initial results show that surface thermal heterogeneity has a significant influence on local atmospheric structure and wind fields in the vicinity of the lake. This influence occurs particularly in the morning. However, dynamic channeling effects and the larger-scale thermal effect of the mountain region frequently override these more local features later in the day.

Genomovar diversity amongst Burkholderia cepacia complex isolates from an Australian adult cystic fibrosis unit

In this study, a combination of recA-based PCR assays and 16S rDNA restriction fragment length polymorphism (RFLP) analysis was used to determine the genomovar diversity of clinical Burkholderia cepacia complex isolates. Twenty-eight isolates were prospectively collected from patients attending a large Australian adult cystic fibrosis (CF) unit, 22 isolates were referred from other Australian CF units and a further eight isolates originated from patients without CF. The 28 prospectively collected isolates were distributed amongst the following genomovars: Burkholderia cepacia genomovar I (28.6%), Burkholderia multivorans (21.4%), Burkholderia cepacia genomovar III (39.3%), Burkholderia vietnamiensis (3.6%) and Burkholderia ambifaria (7.1%). The results of this study highlight the usefulness of 16S rDNA RFLP typing for the identification of other Burkholderia spp. and non-fermenting gram-negative bacteria.

Phylogeny of the Sporobolus indicus complex, based on internal transcribed spacer (ITS) sequences

The entire internal transcribed spacer ( ITS) region, including the 5.8S subunit of the nuclear ribosomal DNA ( rDNA), was sequenced by direct double-stranded sequencing of polymerase chain reaction (PCR) amplified fragments. The study included 40 Sporobolus ( Family Poaceae, subfamily Chloridoideae) seed collections from 14 putative species ( all 11 species from the S. indicus complex and three Australian native species). These sequences, along with those from two out-group species [ Pennisetum alopecuroides ( L.) Spreng. and Heteropogon contortus ( L.) P. Beauv. ex Roemer & Schultes, Poaceae, subfamily Panicoideae], were analysed by the parsimony method (PAUP; version 4.0b4a) to infer phylogenetic relationships among these species. The length of the ITS1, 5.8S subunit and ITS2 region were 222, 164 and 218 base pairs ( bp), respectively, in all species of the S. indicus complex, except for the ITS2 region of S. diandrus P. Beauv. individuals, which was 217 bp long. Of the 624 characters included in the analysis, 245 ( 39.3%) of the 330 variable sites contained potential phylogenetic information. Differences in sequences among the members of the S. pyramidalis P. Beauv., S. natalensis (Steud.) Dur & Schinz and S. jacquemontii Kunth. collections were 0%, while differences ranged from 0 to 2% between these and other species of the complex. Similarly, differences in sequences among collections of S. laxus B. K. Simon, S. sessilis B. K. Simon, S. elongatus R. Br. and S. creber De Nardi were 0%, compared with differences of 1-2% between these four species and the rest of the complex. When comparing S. fertilis ( Steud.) Clayton and S. africanus (Poir.) Robyns & Tourney, differences between collections ranged from 0 to 1%. Parsimony analysis grouped all 11 species of the S. indicus complex together, indicating a monophyletic origin. For the entire data set, pair-wise distances among members of the S. indicus complex varied from 0.00 to 1.58%, compared with a range of 20.08-21.44% among species in the complex and the Australian native species studied. A strict consensus phylogenetic tree separated 11 species of the S. indicus complex into five major clades. The phylogeny, based on ITS sequences, was found to be congruent with an earlier study on the taxonomic relationship of the weedy Sporobolus grasses revealed from random amplified polymorphic DNA ( RAPD). However, this cladistic analysis of the complex was not in agreement with that created on past morphological analyses and therefore gives a new insight into the phylogeny of the S. indicus complex.