Stomatal frequency adjustment of four conifer species to historical changes in atmospheric CO2

The species-specific inverse relation between atmospheric CO(2) concentration and stomatal frequency for many woody angiosperm species is being used increasingly with fossil leaves to reconstruct past atmospheric CO(2) levels. To extend our limited knowledge of the responsiveness of conifer needles to CO(2) fluctuations, the stomatal frequency response of four native North American conifer species (Tsuga heterophylla, Picea glauca, Picea mariana, and Larix laricina) to a range of historical CO(2) mixing ratios (290 to 370 ppmV) was analyzed. Because of the specific mode of leaf development and the subsequent stomatal patterning in conifer needles, the stomatal index of these species was not affected by CO(2). In contrast, a new measure of stomatal frequency, based on the number of stomata per millimeter of needle length, decreased significantly with increasing CO(2). For Tsuga heterophylla, the stomatal frequency response to CO(2) changes in the last century is validated through assessment of the influence of other biological and environmental variables. Because of their sensitive response to CO(2), combined with a high preservation capacity, fossil needles of Tsuga heterophylla, Picea glauca, P. mariana, and Larix laricina have great potential for detecting and quantifying past atmospheric CO(2) fluctuations.     

Interpreting scratch assays using pair density dynamics and approximate Bayesian computation

Quantifying the impact of biochemical compounds on collective cell spreading is an essential element of drug design, with various applications including developing treatments for chronic wounds and cancer. Scratch assays are a technically simple and inexpensive method used to study collective cell spreading; however, most previous interpretations of scratch assays are qualitative and do not provide estimates of the cell diffusivity, D, or the cell proliferation rate, λ. Estimating D and λ is important for investigating the efficacy of a potential treatment and provides insight into the mechanism through which the potential treatment acts. While a few methods for estimating D and λ have been proposed, these previous methods lead to point estimates of D and λ, and provide no insight into the uncertainty in these estimates. Here, we compare various types of information that can be extracted from images of a scratch assay, and quantify D and λ using discrete computational simulations and approximate Bayesian computation. We show that it is possible to robustly recover estimates of D and λ from synthetic data, as well as a new set of experimental data. For the first time, our approach also provides a method to estimate the uncertainty in our estimates of D and λ. We anticipate that our approach can be generalized to deal with more realistic experimental scenarios in which we are interested in estimating D and λ, as well as additional relevant parameters such as the strength of cell-to-cell adhesion or the strength of cell-to-substrate adhesion.     

Statistical analysis of an RNA titration series evaluates microarray precision and sensitivity on a whole-array basis

Background

It is one of the ultimate goals for modern biological research to fully elucidate the intricate interplays and the regulations of the molecular determinants that propel and characterize the progression of versatile life phenomena, to name a few, cell cycling, developmental biology, aging, and the progressive and recurrent pathogenesis of complex diseases. The vast amount of large-scale and genome-wide time-resolved data is becoming increasing available, which provides the golden opportunity to unravel the challenging reverse-engineering problem of time-delayed gene regulatory networks.

Results

In particular, this methodological paper aims to reconstruct regulatory networks from temporal gene expression data by using delayed correlations between genes, i.e., pairwise overlaps of expression levels shifted in time relative each other. We have thus developed a novel model-free computational toolbox termed TdGRN (Time-delayed Gene Regulatory Network) to address the underlying regulations of genes that can span any unit(s) of time intervals. This bioinformatics toolbox has provided a unified approach to uncovering time trends of gene regulations through decision analysis of the newly designed time-delayed gene expression matrix. We have applied the proposed method to yeast cell cycling and human HeLa cell cycling and have discovered most of the underlying time-delayed regulations that are supported by multiple lines of experimental evidence and that are remarkably consistent with the current knowledge on phase characteristics for the cell cyclings.

Conclusion

We established a usable and powerful model-free approach to dissecting high-order dynamic trends of gene-gene interactions. We have carefully validated the proposed algorithm by applying it to two publicly available cell cycling datasets. In addition to uncovering the time trends of gene regulations for cell cycling, this unified approach can also be used to study the complex gene regulations related to the development, aging and progressive pathogenesis of a complex disease where potential dependences between different experiment units might occurs.     

Corticosteroids prevent myofibroblast accumulation and airway remodeling in mice

At present there are conflicting results from studies investigating the role of corticosteroids in inhibiting airway remodeling in asthma. We have used a mouse model to determine whether administration of corticosteroids prevents the development of allergen-induced structural features of airway remodeling. Mice treated with corticosteroids were subjected to repetitive ovalbumin (OVA) challenge for 3 mo, at which time levels of peribronchial fibrosis and the thickness of the peribronchial smooth muscle layer were assessed by immunohistology, levels of transforming growth factor (TGF)-beta1 by ELISA, and the number of alpha-smooth muscle actin+/Col-1+ peribronchial myofibroblasts by immunohistochemistry. Corticosteroids significantly reduced allergen-induced increases in peribronchial collagen deposition and levels of total lung collagen but did not reduce allergen-induced increases in the thickness of the peribronchial smooth muscle layer. Levels of lung TGF-beta1 were significantly reduced in mice treated with systemic corticosteroids, and this was associated with a significant decrease in the number of peribronchial inflammatory cells that expressed TGF-beta1, including eosinophils and mononuclear cells. Corticosteroids also significantly reduced the number of peribronchial myofibroblasts. Overall, these studies demonstrate that administration of corticosteroids significantly reduces levels of allergen-induced peribronchial fibrosis. The reduction in peribronchial fibrosis mediated by corticosteroids is likely to be due to several mechanisms including inhibition of expression of TGF-beta1, a reduction in the number of peribronchial inflammatory cells expressing TGF-beta1 (eosinophils, macrophages), as well as by corticosteroids reducing the accumulation of peribronchial myofibroblasts that contribute to collagen expression.     

Novel O-palmitolylated beta-E1 subunit of pyruvate dehydrogenase is phosphorylated during ischemia/reperfusion injury

Background  

During and following myocardial ischemia, glucose oxidation rates are low and fatty acids dominate as a source of oxidative metabolism. This metabolic phenotype is associated with contractile dysfunction during reperfusion. To determine the mechanism of this reliance on fatty acid oxidation as a source of ATP generation, a functional proteomics approach was utilized.     

Differences in the response sensitivity of stomatal index to atmospheric CO2 among four genera of Cupressaceae conifers

Background and Aims

The inverse relationship between stomatal density (SD: number of stomata per mm² leaf area) and atmospheric concentration of CO₂ ([CO₂]) permits the use of plants as proxies of palaeo-atmospheric CO₂. Many stomatal reconstructions of palaeo-[CO₂] are based upon multiple fossil species. However, it is unclear how plants respond to [CO₂] across genus, family or ecotype in terms of SD or stomatal index (SI: ratio of stomata to epidermal cells). This study analysed the stomatal numbers of conifers from the ancient family Cupressaceae, in order to examine the nature of the SI–[CO₂] relationship, and potential implications for stomatal reconstructions of palaeo-[CO₂].

Methods

Stomatal frequency measurements were taken from historical herbarium specimens of Athrotaxis cupressoides, Tetraclinis articulata and four Callitris species, and live A. cupressoides grown under CO₂-enrichment (370, 470, 570 and 670 p.p.m. CO₂).

Key Results

T. articulata, C. columnaris and C. rhomboidea displayed significant reductions in SI with rising [CO₂]; by contrast, A. cupressoides, C. preissii and C. oblonga show no response in SI. However, A. cupressoides does reduce SI to increases in [CO₂] above current ambient (approx. 380 p.p.m. CO₂). This dataset suggests that a shared consistent SI–[CO₂] relationship is not apparent across the genus Callitris.

Conclusions

The present findings suggest that it is not possible to generalize how conifer species respond to fluctuations in [CO₂] based upon taxonomic relatedness or habitat. This apparent lack of a consistent response, in conjunction with high variability in SI, indicates that reconstructions of absolute palaeo-[CO₂] based at the genus level, or upon multiple species for discrete intervals of time are not as reliable as those based on a single or multiple temporally overlapping species.     

A mathematical model of tumour-induced capillary growth

The corneal limbal vessels of an animal host respond to the presence of a source of Tumour Angiogenesis Factor (TAF) implanted in the cornea by the formation of new capillaries which grow towards the source. This neovasculature can be easily seen and studied and this paper describes a mathematical model of some of the important features of the growth. The model includes the diffusion of TAF, the formation of sprouts from pre-existing vessels and models the movement of these sprouts to form new capillaries as a chemotactic response to the presence of TAF. Numerical results are produced for various values of the parameters which characterize the model and it is suggested that the model might form the framework for further theoretical work on related phenomena such as wound healing or to develop strategies for the investigation of anti-angiogenesis.     

Computational Approaches to Solving Equations Arising from Wound Healing

In the wound healing process, the cell movement associated with chemotaxis generally outweighs the movement associated with random motion, leading to advection-dominated mathematical models of wound healing. The equations in these models must be solved with care, but often inappropriate approaches are adopted. Two one-dimensional test problems arising from advection-dominated models of wound healing are solved using four algorithms—MATLAB’s inbuilt routine pdepe.m, the Numerical Algorithms Group routine d03pcf.f, and two finite volume methods. The first finite volume method is based on a first-order upwinding treatment of chemotaxis terms and the second on a flux limiting approach. The first test problem admits an analytic solution which can be used to validate the numerical results by analyzing two measures of the error for each method: the average absolute difference and a mass balance error. These criteria as well as the visual comparison between the numerical methods and the exact solution lead us to conclude that flux limiting is the best approach to solving advection-dominated wound healing problems numerically in one dimension. The second test problem is a coupled nonlinear three species model of wound healing angiogenesis. Measurement of the mass balance error for this test problem further confirms our hypothesis that flux limiting is the most appropriate method for solving advection-dominated governing equations in wound healing models. We also consider two two-dimensional test problems arising from wound healing, one that admits an analytic solution and a more complicated problem of blood vessels growth into a devascularized wound bed. The results from the two-dimensional test problems also demonstrate that the flux limiting treatment of advective terms is ideal for an advection-dominated problem.     

Statistical analysis of an RNA titration series evaluates microarray precision and sensitivity on a whole-array basis

Background  

Concerns are often raised about the accuracy of microarray technologies and the degree of cross-platform agreement, but there are yet no methods which can unambiguously evaluate precision and sensitivity for these technologies on a whole-array basis.