Understanding Subsystems in Biology through Dimensionality Reduction, Graph Partitioning and Analytical Modeling

Biological systems exhibit rich and complex behavior through the orchestrated interplay of a large array of components. It is hypothesized that separable subsystems with some degree of functional autonomy exist; deciphering their independent behavior and functionality would greatly facilitate understanding the system as a whole. Discovering and analyzing such subsystems are hence pivotal problems in the quest to gain a quantitative understanding of complex biological systems. In this work, using approaches from machine learning, physics and graph theory, methods for the identification and analysis of such subsystems were developed. A novel methodology, based on a recent machine learning algorithm known as non-negative matrix factorization (NMF), was developed to discover such subsystems in a set of large-scale gene expression data. This set of subsystems was then used to predict functional relationships between genes, and this approach was shown to score significantly higher than conventional methods when benchmarking them against existing databases. Moreover, a mathematical treatment was developed to treat simple network subsystems based only on their topology (independent of particular parameter values). Application to a problem of experimental interest demonstrated the need for extentions to the conventional model to fully explain the experimental data. Finally, the notion of a subsystem was evaluated from a topological perspective. A number of different protein networks were examined to analyze their topological properties with respect to separability, seeking to find separable subsystems. These networks were shown to exhibit separability in a nonintuitive fashion, while the separable subsystems were of strong biological significance. It was demonstrated that the separability property found was not due to incomplete or biased data, but is likely to reflect biological structure.

Energy partitioning for "fuzzy" atoms

The total energy of molecule in terms of 'fuzzy atoms' presented as sum of one- and two-atomic energy components is described. The divisions of three-dimensional physical space into atomic regions exhibit continuous transition from one to another. The energy components are on chemical energy scale according to proper definitions. The Becke's integration scheme and weight function determines realization of method which permits effective numerical integrations

One- and two-center physical space partitioning of the energy in the density functional theory

A conceptually new approach is introduced for the decomposition of the molecular energy calculated at the density functional theory level of theory into sum of one- and two-atomic energy components, and is realized in the "fuzzy atoms" framework. (Fuzzy atoms mean that the three-dimensional physical space is divided into atomic regions having no sharp boundaries but exhibiting a continuous transition from one to another.) The new scheme uses the new concept of "bond order density" to calculate the diatomic exchange energy components and gives them unexpectedly close to the values calculated by the exact (Hartree-Fock) exchange for the same Kohn-Sham orbitals

Predicting Argentine ant invasion across spatial scales via ecological niche models

La formiga argentina (Linepithema humile) es troba entre les espècies més invasores: originària d'Amèrica del Sud, actualment ha envaït nombroses àrees arreu del món. Aquesta tesi doctoral intenta fer una primera anàlisi integrada i multiescalar de la distribució de la formiga argentina mitjançant l'ús de models de nínxol ecològic. D'acord amb els resultats obtinguts, es preveu que la formiga argentina assoleixi una distribució més àmplia que l'actual. Les prediccions obtingudes a partir dels models concorden amb la distribució actualment coneguda i, a més, indiquen àrees a prop de la costa i dels rius principals com a altament favorables per a l'espècie. Aquests resultats corroboren la idea que la formiga argentina no es troba actualment en equilibri amb el medi. D'altra banda, amb el canvi climàtic, s'espera que la distribució de la formiga argentina s'estengui cap a latituds més elevades en ambdós hemisferis, i sofreixi una retracció en els tròpics a escales globals.

Comparison of subcellular partitioning, distribution, and internal speciation of Cu between Cu-tolerant and naive populations of Dendrodrilus rubidus Savigny

When considering contaminated site ecology and ecological risk assessment a key question is whether organisms that appear unaffected by accumulation of contaminants are tolerant or resistant to those contaminants. A population of Dendrodrilus rubidus Savigny earthworms from the Coniston Copper Mines, an area of former Cu mining, exhibit increased tolerance and accumulation of Cu relative to a nearby non-Cu exposed population. Distribution of total Cu between different body parts (posterior, anterior, body wall) of the two populations was determined after a 14 day exposure to 250 mg Cu kg(-1) in Cu-amended soil. Cu concentrations were greater in Coniston earthworms but relative proportions of Cu in different body parts were the same between populations. Cu speciation was determined using extended X-ray absorption fine structure spectroscopy (EXAFS). Cu was coordinated to 0 atoms in the exposure soil but to S atoms in the earthworms. There was no difference in this speciation between the different earthworm populations. In another experiment earthworms were exposed to a range of Cu concentrations (200-700 mg Cu kg(-1)). Subcellular partitioning of accumulated Cu was determined. Coniston earthworms accumulated more Cu but relative proportions of Cu in the different fractions (cytosol > granular > tissue fragments, cell membranes, and intact cells) were the same between populations. Results suggest that Coniston D. rubidus are able to survive in the Cu-rich Coniston Copper Mines soil through enlargement of the same Cu storage reservoirs that exist in a nearby non-Cu exposed population.

Trace element partitioning and soil particle characterisation around mining and smelting areas at Tharsis, Riotinto and Huelva, SW Spain

Trace elements may present an environmental hazard in the vicinity of mining and smelting activities. However, the factors controlling trace element distribution in soils around ancient and modem mining and smelting areas are not always clear. Tharsis, Riotinto and Huelva are located in the Iberian Pyrite Belt in SW Spain. Tharsis and Riotinto mines have been exploited since 2500 B.C., with intensive smelting taking place. Huelva, established in 1970 and using the Flash Furnace Outokumpu process, is currently one of the largest smelter in the world. Pyrite and chalcopyrite ore have been intensively smelted for Cu. However, unusually for smelters and mines of a similar size, the elevated trace element concentrations in soils were found to be restricted to the immediate vicinity of the mines and smelters, being found up to a maximum of 2 kin from the mines and smelters at Tharsis, Riotinto and Huelva. Trace element partitioning (over 2/3 of trace elements found in the residual immobile fraction of soils at Tharsis) and soil particles examination by SEM-EDX showed that trace elements were not adsorbed onto soil particles, but were included within the matrix of large trace element-rich Fe silicate slag particles (i.e. 1 min circle divide at least 1 wt.% As, Cu and Zn, and 2 wt.% Pb). Slag particle large size (I mm 0) was found to control the geographically restricted trace element distribution in soils at Tharsis, Riotinto and Huelva, since large heavy particles could not have been transported long distances. Distribution and partitioning indicated that impacts to the environment as a result of mining and smelting should remain minimal in the region. (c) 2006 Elsevier B.V. All rights reserved.

A comparative evaluation of functions for partitioning nitrogen and amino acid intake between maintenance and growth in broilers

The results from three types of study with broilers, namely nitrogen (N) balance, bioassays and growth experiments, provided the data used herein. Sets of data on N balance and protein accretion (bioassay studies) were used to assess the ability of the monomolecular equation to describe the relationship between (i) N balance and amino acid (AA) intake and (ii) protein accretion and AA intake. The model estimated the levels of isoleucine, lysine, valine, threonine, methionine, total sulphur AAs and tryptophan resulting in zero balance to be 58, 59, 80, 96, 23, 85 and 32 mg/kg live weight (LW)/day, respectively. These estimates show good agreement with those obtained in previous studies. For the growth experiments, four models, specifically re-parameterized for analysing energy balance data, were evaluated for their ability to determine crude protein (CP) intake at maintenance and efficiency of utilization of CP intake for producing gain. They were: a straight line, two equations representing diminishing returns behaviour (monomolecular and rectangular hyperbola) and one equation describing smooth sigmoidal behaviour with a fixed point of inflexion (Gompertz). The estimates of CP requirement for maintenance and efficiency of utilization of CP intake for producing gain varied from 5.4 to 5.9 g/kg LW/day and 0.60 to 0.76, respectively, depending on the models.

The effects of triazole and strobilurin fungicide programmes on nitrogen uptake, partitioning, remobilization and grain N accumulation in winter wheat cultivars

Field experiments were conducted over 3 years to assess the effect of a triazole fungicide programme, and additions of strobilurin fungicides to it, on nitrogen uptake, accumulation and partitioning in a range of winter wheat cultivars. Commensurate with delayed senescence, fungicide programmes, particularly when including strobilurins, improved grain yield through improvements in both crop biomass and harvest index, although the relationship with green area duration of the flag leaf (GFLAD) depended on year and in some cases, cultivar. In all years fungicide treatments significantly increased the amount of nitrogen in the above-ground biomass, the amount of nitrogen in the grain and the nitrogen harvest index. All these effects could be linearly related to the fungicide effect on GFLAD. These relationships occasionally interacted with cultivar but there was no evidence that fungicide mode of action affected the relationship between GFLAD and yield of nitrogen in the grain. Fungicide treatments significantly reduced the amount of soil mineral N at harvest and when severe disease had been controlled, the net remobilization of N from the vegetation to the grain after anthesis. Fungicide maintained the filling of grain with both dry matter and nitrogen. The proportionate accumulation of nitrogen in the grain was later than that of dry matter and this difference was greater when fungicide had been applied. Effects of fungicide on grain protein concentration and its relationship with GFLAD were inconsistent over year and cultivar. There were several instances where grain protein concentration was unaffected despite large (1(.)5 t/ha) increases in grain yield following fungicide use. Dilution of grain protein concentration following fungicide use, when it did occur, was small compared with what would be predicted by adoption of other yield increasing techniques such as the selection of high yielding cultivars (based on currently available cultivars) or by growing wheat in favourable climates.

Modelling the partitioning of evapotranspiration in a maize-sunflower intercrop

The primary purpose of this study was to model the partitioning of evapotranspiration in a maize-sunflower intercrop at various canopy covers. The Shuttleworth-Wallace (SW) model was extended for intercropping systems to include both crop transpiration and soil evaporation and allowing interaction between the two. To test the accuracy of the extended SW model, two field experiments of maize-sunflower intercrop were conducted in 1998 and 1999. Plant transpiration and soil evaporation were measured using sap flow gauges and lysimeters, respectively. The mean prediction error (simulated minus measured values) for transpiration was zero (which indicated no overall bias in estimation error), and its accuracy was not affected by the plant growth stages, but simulated transpiration during high measured transpiration rates tended to be slightly underestimated. Overall, the predictions for daily soil evaporation were also accurate. Model estimation errors were probably due to the simplified modelling of soil water content, stomatal resistances and soil heat flux as well as due to the uncertainties in characterising the 2 micrometeorological conditions. The SW’s prediction of transpiration was most sensitive to parameters most directly related to the canopy characteristics such as the partitioning of captured solar radiation, canopy resistance, and bulk boundary layer resistance.