Automated Bayesian model development for frequency detection in biological time series

Background

The ability to perform quantitative studies using isotope tracers and metabolic flux analysis (MFA) is critical for detecting pathway bottlenecks and elucidating network regulation in biological systems, especially those that have been engineered to alter their native metabolic capacities. Mathematically, MFA models are traditionally formulated using separate state variables for reaction fluxes and isotopomer abundances. Analysis of isotope labeling experiments using this set of variables results in a non-convex optimization problem that suffers from both implementation complexity and convergence problems.

Results

This article addresses the mathematical and computational formulation of ¹³C MFA models using a new set of variables referred to as fluxomers. These composite variables combine both fluxes and isotopomer abundances, which results in a simply-posed formulation and an improved error model that is insensitive to isotopomer measurement normalization. A powerful fluxomer iterative algorithm (FIA) is developed and applied to solve the MFA optimization problem. For moderate-sized networks, the algorithm is shown to outperform the commonly used 13CFLUX cumomer-based algorithm and the more recently introduced OpenFLUX software that relies upon an elementary metabolite unit (EMU) network decomposition, both in terms of convergence time and output variability.

Conclusions

Substantial improvements in convergence time and statistical quality of results can be achieved by applying fluxomer variables and the FIA algorithm to compute best-fit solutions to MFA models. We expect that the fluxomer formulation will provide a more suitable basis for future algorithms that analyze very large scale networks and design optimal isotope labeling experiments.     

Small tropical forest trees have a greater capacity to adjust carbon metabolism to long-term drought than large canopy trees

The response of small understory trees to long-term drought is vital in determining the future composition, carbon stocks and dynamics of tropical forests. Long-term drought is, however, also likely to expose understory trees to increased light availability driven by drought-induced mortality. Relatively little is known about the potential for understory trees to adjust their physiology to both decreasing water and increasing light availability. We analysed data on maximum photosynthetic capacity (J_max, V_cmax), leaf respiration (R_leaf), leaf mass per area (LMA), leaf thickness and leaf nitrogen and phosphorus concentrations from 66 small trees across 12 common genera at the world's longest running tropical rainfall exclusion experiment and compared responses to those from 61 surviving canopy trees. Small trees increased J_max, V_cmax, R_leaf and LMA (71, 29, 32, 15% respectively) in response to the drought treatment, but leaf thickness and leaf nutrient concentrations did not change. Small trees were significantly more responsive than large canopy trees to the drought treatment, suggesting greater phenotypic plasticity and resilience to prolonged drought, although differences among taxa were observed. Our results highlight that small tropical trees have greater capacity to respond to ecosystem level changes and have the potential to regenerate resilient forests following future droughts.     

The diabetic rat kidney mediates inosituria and selective urinary partitioning of D-chiro-inositol

Diabetic nephropathy is a serious complication of diabetes mellitus with a pressing need for effective metabolic markers to detect renal impairment. Of potential significance are the inositol compounds, myo-inositol (MI), and the less abundant stereoisomer, D-chiro-inositol (DCI), which are excreted at increased levels in the urine in diabetes mellitus, a phenomenon known as inosituria. There is also a selective urinary excretion of DCI compared to MI. As the biological origins of altered inositol metabolism in diabetes mellitus are unknown, the aim of this study was to determine whether the diabetic kidney was directly responsible. Kidneys isolated from four-week streptozotocin-induced diabetic rats were characterized by a 3-fold reduction in glomerular filtration rate (GFR) compared to matched non-diabetic kidneys. When perfused with fixed quantities of MI (50 µM) and DCI (5 µM) under normoglycemic conditions (5 mM glucose), GFR-normalized urinary excretion of MI was increased by 1.7-fold in diabetic vs. non-diabetic kidneys. By comparison, GFR-normalized urinary excretion of DCI was increased by 4-fold. Perfusion conditions replicating hyperglycemia (20 mM glucose) potentiated DCI but not MI urinary excretion in both non-diabetic and diabetic kidneys. Overall, there was a 2.4-fold increase in DCI urinary excretion compared to MI in diabetic kidneys that was independent of glucose ambience. This increased urinary excretion of DCI and MI in diabetic kidneys occurred despite increased renal expression of the inositol transporters, sodium myo-inositol transporter subtype 1 and 2 (SMIT1 and SMIT2). These findings show that the diabetic kidney primarily mediates inosituria and altered urinary partitioning of MI and DCI. Urinary inositol levels might therefore serve as an indicator of impaired renal function in diabetes mellitus with wider implications for monitoring chronic kidney disease.     

Temporal and spatial variability in nearshore bacterioplankton communities of Lake Michigan

The spatial and temporal variability of bacterial communities were determined for the nearshore waters of Lake Michigan, an oligotrophic freshwater inland sea. A freshwater estuary and nearshore sites were compared six times during 2006 using denaturing gradient gel electrophoresis (DGGE). Bacterial composition clustered by individual site and date rather than by depth. Seven 16S rRNA gene clone libraries were constructed, yielding 2717 bacterial sequences. Spatial variability was detected among the DGGE banding patterns and supported by clone library composition. The clone libraries from deep waters and the estuary environment revealed highest overall bacterial diversity. Betaproteobacteria sequence types were the most dominant taxa, comprising 40.2–67.7% of the clone libraries. BAL 47 was the most abundant freshwater cluster of Betaproteobacteria , indicating widespread distribution of this cluster in the nearshore waters of Lake Michigan. Incertae sedis 5 and Oxalobacteraceae sequence types were prevalent in each clone library, displaying more diversity than previously described in other freshwater environments. Among the Oxalobacteraceae sequences, a globally distributed freshwater cluster was determined. The nearshore waters of Lake Michigan are a dynamic environment that experience forces similar to the coastal ocean environment and share common bacterial diversity with other freshwater habitats.     

Hexavalent chromium reduction in <Emphasis Type="Italic">Desulfovibrio vulgaris</Emphasis> Hildenborough causes transitory inhibition of sulfate reduction and cell growth

Desulfovibrio vulgaris Hildenborough is a well-studied sulfate reducer that can reduce heavy metals and radionuclides [e.g., Cr(VI) and U(VI)]. Cultures grown in a defined medium had a lag period of approximately 30 h when exposed to 0.05 mM Cr(VI). Substrate analyses revealed that although Cr(VI) was reduced within the first 5 h, growth was not observed for an additional 20 h. The growth lag could be explained by a decline in cell viability; however, during this time small amounts of lactate were still utilized without sulfate reduction or acetate formation. Approximately 40 h after Cr exposure (0.05 mM), sulfate reduction occurred concurrently with the accumulation of acetate. Similar amounts of hydrogen were produced by Cr-exposed cells compared to control cells, and lactate was not converted to glycogen during non-growth conditions. D. vulgaris cells treated with a reducing agent and then exposed to Cr(VI) still experienced a growth lag, but the addition of ascorbate at the time of Cr(VI) addition prevented the lag period. In addition, cells grown on pyruvate displayed more tolerance to Cr(VI) compared to lactate-grown cells. These results indicated that D. vulgaris utilized lactate during Cr(VI) exposure without the reduction of sulfate or production of acetate, and that ascorbate and pyruvate could protect D. vulgaris cells from Cr(VI)/Cr(III) toxicity. J.D. Wall and M.W. Fields are both affiliated to the Virtual Institute of Microbial Stress and Survival (). M.E. Clark and S.B. Thieman contributed equally to this work.     

The NFP locus of Medicago truncatula controls an early step of Nod factor signal transduction upstream of a rapid calcium flux and root hair deformation

Establishment of the Rhizobium-legume symbiosis depends on a molecular dialogue, in which rhizobial nodulation (Nod) factors act as symbiotic signals, playing a key role in the control of specificity of infection and nodule formation. Using nodulation-defective (Nod-) mutants of Medicago truncatula to study the mechanisms controlling Nod factor perception and signalling, we have previously identified five genes that control components of a Nod factor-activated signal transduction pathway. Characterisation of a new M. truncatula Nod- mutant led to the identification of the Nod Factor Perception (NFP) locus. The nfp mutant has a novel phenotype among Nod- mutants of M. truncatula, as it does not respond to Nod factors by any of the responses tested. The nfp mutant thus shows no rapid calcium flux, the earliest detectable Nod factor response of wild-type plants, and no root hair deformation. The nfp mutant is also deficient in Nod factor-induced calcium spiking and early nodulin gene expression. While certain genes controlling Nod factor signal transduction also control the establishment of an arbuscular mycorrhizal symbiosis, the nfp mutant shows a wild-type mycorrhizal phenotype. These data indicate that the NFP locus controls an early step of Nod factor signal transduction, upstream of previously identified genes and specific to nodulation.     

Regulation of Inflammation by IL-17A and IL-17F Modulates Non-Alcoholic Fatty Liver Disease Pathogenesis

Non-alcoholic fatty liver disease (NAFLD) has become the most common chronic liver disease worldwide. While it is well-accepted that inflammation is central to NAFLD pathogenesis, the immune pathway(s) orchestrating disease progression are poorly defined. Notably, IL-17RA signaling, via IL-17A, plays an important role in obesity-driven NAFLD pathogenesis. However, the role of the IL-17F, another IL-17RA ligand, in NAFLD pathogenesis has not been examined. Further, the cell types expressing IL-17RA and producing IL-17RA ligands in the pathogenesis of NAFLD have not been defined. Here, IL-17RA^-/-, IL-17A^-/-, IL-17F^-/- and wild-type (WT) mice were fed either standard chow diet or methionine and choline deficient diet (MCDD)—a diet known to induce steatosis and hepatic inflammation through beta-oxidation dysfunction—and hepatic inflammation and NAFLD progression were subsequently quantified. MCDD feeding augmented hepatic IL-17RA expression and significantly increased hepatic infiltration of macrophages and IL-17A and IL-17F producing CD4⁺ and CD8⁺ T cells in WT mice. In contrast, IL-17RA^-/-, IL-17A^-/-, and IL-17F^-/- mice, despite increased steatosis, exhibited significant protection from hepatocellular damage compared to WT controls. Protection from hepatocellular damage correlated with decreased levels of hepatic T-cell and macrophage infiltration and decreased expression of inflammatory mediators associated with NAFLD. In sum, our results indicate that the IL-17 axis also plays a role in a MCDD-induced model of NAFLD pathogenesis. Further, we show for the first time that IL-17F, and not only IL-17A, plays an important role in NAFLD driven inflammation.     

Unexpected oocyte growth after follicular antrum formation in four marsupial species.

During examination of maturing preovulatory marsupial oocytes we noted that oocyte diameters were invariably about 50% greater than the figures reported in earlier histological studies. As all previous investigations were limited to small follicles (at most 25% the size of the ovulating follicle), the present study was initiated to examine oocyte growth during the whole period of follicular development. Oocyte and follicle diameters were measured for three Australian (Trichosurus vulpecula, Macropus eugenii and Bettongia penicillata--fresh nonfixed material) and one American marsupial species (Monodelphis domestica--histological sections) in which multiple follicle development had been induced by exogenous gonadotrophin treatment. In all species oocytes were obtained from follicles ranging from pre-antral to immediately pre-ovulatory (maximum follicle sizes obtained were: T. vulpecula, 4.5 mm; M. eugenii, 4.3 mm; B. penicillata, 2.5 mm; M. domestica, 0.7 mm). In two of the species (T. vulpecula and B. penicillata) ovulated oocytes were also examined. In T. vulpecula and M. eugenii oocytes were found to achieve much greater diameters than previously reported from histological studies of small follicles (< 0.8 mm) and similar patterns of growth were found in the other two species. In the four species oocytes reached diameters about two to three times that found for eutherian mammals. It was concluded that the marsupial oocyte continued to grow after formation of the follicular antrum and that, although the rate of oocyte growth slowed in larger follicles, it continued into the period immediately before ovulation. In B. penicillata the largest oocytes were obtained after ovulation.(ABSTRACT TRUNCATED AT 250 WORDS)