Ecological and methodological drivers of non-stationarity in tree growth response to climate

Radial tree growth is sensitive to environmental conditions, making observed growth increments an important indicator of climate change effects on forest growth. However, unprecedented climate variability could lead to non-stationarity, that is, a decoupling of tree growth responses from climate over time, potentially inducing biases in climate reconstructions and forest growth projections. Little is known about whether and to what extent environmental conditions, species, and model type and resolution affect the occurrence and magnitude of non-stationarity. To systematically assess potential drivers of non-stationarity, we compiled tree-ring width chronologies of two conifer species, Picea abies and Pinus sylvestris, distributed across cold, dry, and mixed climates. We analyzed 147 sites across the Europe including the distribution margins of these species as well as moderate sites. We calibrated four numerical models (linear vs. non-linear, daily vs. monthly resolution) to simulate growth chronologies based on temperature and soil moisture data. Climate–growth models were tested in independent verification periods to quantify their non-stationarity, which was assessed based on bootstrapped transfer function stability tests. The degree of non-stationarity varied between species, site climatic conditions, and models. Chronologies of P. sylvestris showed stronger non-stationarity compared with Picea abies stands with a high degree of stationarity. Sites with mixed climatic signals were most affected by non-stationarity compared with sites sampled at cold and dry species distribution margins. Moreover, linear models with daily resolution exhibited greater non-stationarity compared with monthly-resolved non-linear models. We conclude that non-stationarity in climate–growth responses is a multifactorial phenomenon driven by the interaction of site climatic conditions, tree species, and methodological features of the modeling approach. Given the existence of multiple drivers and the frequent occurrence of non-stationarity, we recommend that temporal non-stationarity rather than stationarity should be considered as the baseline model of climate–growth response for temperate forests.     

Large old trees increase growth under shifting climatic constraints: Aligning tree longevity and individual growth dynamics in primary mountain spruce forests

In a world of accelerating changes in environmental conditions driving tree growth, tradeoffs between tree growth rate and longevity could curtail the abundance of large old trees (LOTs), with potentially dire consequences for biodiversity and carbon storage. However, the influence of tree-level tradeoffs on forest structure at landscape scales will also depend on disturbances, which shape tree size and age distribution, and on whether LOTs can benefit from improved growing conditions due to climate warming. We analyzed temporal and spatial variation in radial growth patterns from ~5000 Norway spruce (Picea abies [L.] H. Karst) live and dead trees from the Western Carpathian primary spruce forest stands. We applied mixed-linear modeling to quantify the importance of LOT growth histories and stand dynamics (i.e., competition and disturbance factors) on lifespan. Finally, we assessed regional synchronization in radial growth variability over the 20th century, and modeled the effects of stand dynamics and climate on LOTs recent growth trends. Tree age varied considerably among forest stands, implying an important role of disturbance as an age constraint. Slow juvenile growth and longer period of suppressed growth prolonged tree lifespan, while increasing disturbance severity and shorter time since last disturbance decreased it. The highest age was not achieved only by trees with continuous slow growth, but those with slow juvenile growth followed by subsequent growth releases. Growth trend analysis demonstrated an increase in absolute growth rates in response to climate warming, with late summer temperatures driving the recent growth trend. Contrary to our expectation that LOTs would eventually exhibit declining growth rates, the oldest LOTs (>400 years) continuously increase growth throughout their lives, indicating a high phenotypic plasticity of LOTs for increasing biomass, and a strong carbon sink role of primary spruce forests under rising temperatures, intensifying droughts, and increasing bark beetle outbreaks.     

Syntheses of 1-[2-(Phosphonomethoxy)Alkyl]Thymine Monophosphates and an Evaluation of their Inhibitory Activity Toward Human Thymidine Phosphorylase

A series of new monophosphates of 1-[2-(phosphonomethoxy)alkyl]thymines, such as PMPTp_, 3-MeO-PMPTp, HPMPTp, and FPMPTp, were synthesized and tested for their ability to inhibit human thymidine phosphorylase. Kinetic measurements of enzyme activity were performed using thymidine and inorganic phosphate as the substrates. The data show that some monophosphates provide a considerable increase of the multisubstrate inhibitory effect. The highest inhibitory potency was found with (R)-FPMPTp 4c (K _i ^dT = 4.09 ± 0.47 μM, K _i(P_i) = 2.13 ± 0.29 μM) and (R) 3-MeO-PMPTp 4d (K _i ^dT = 5.78 ± 0.71 μM, K _i(P_i) = 2.71 ± 0.37 μM).     

The Validity of Quasi-Steady-State Approximations in Discrete Stochastic Simulations

In biochemical networks, reactions often occur on disparate timescales and can be characterized as either fast or slow. The quasi-steady-state approximation (QSSA) utilizes timescale separation to project models of biochemical networks onto lower-dimensional slow manifolds. As a result, fast elementary reactions are not modeled explicitly, and their effect is captured by nonelementary reaction-rate functions (e.g., Hill functions). The accuracy of the QSSA applied to deterministic systems depends on how well timescales are separated. Recently, it has been proposed to use the nonelementary rate functions obtained via the deterministic QSSA to define propensity functions in stochastic simulations of biochemical networks. In this approach, termed the stochastic QSSA, fast reactions that are part of nonelementary reactions are not simulated, greatly reducing computation time. However, it is unclear when the stochastic QSSA provides an accurate approximation of the original stochastic simulation. We show that, unlike the deterministic QSSA, the validity of the stochastic QSSA does not follow from timescale separation alone, but also depends on the sensitivity of the nonelementary reaction rate functions to changes in the slow species. The stochastic QSSA becomes more accurate when this sensitivity is small. Different types of QSSAs result in nonelementary functions with different sensitivities, and the total QSSA results in less sensitive functions than the standard or the prefactor QSSA. We prove that, as a result, the stochastic QSSA becomes more accurate when nonelementary reaction functions are obtained using the total QSSA. Our work provides an apparently novel condition for the validity of the QSSA in stochastic simulations of biochemical reaction networks with disparate timescales.     

'Genome order index' should not be used for defining compositional constraints in nucleotide sequences - a case study of the Z-curve

Background  

The Z-curve is a three dimensional representation of DNA sequences proposed over a decade ago and has been extensively applied to sequence segmentation, horizontal gene transfer detection, and sequence analysis. Based on the Z-curve, a "genome order index," was proposed, which is defined as S = a ²+ c ²+t ²+g ², where a, c, t, and g are the nucleotide frequencies of A, C, T, and G, respectively. This index was found to be smaller than 1/3 for almost all tested genomes, which was taken as support for the existence of a constraint on genome composition. A geometric explanation for this constraint has been suggested. Each genome was represented by a point P whose distance from the four faces of a regular tetrahedron was given by the frequencies a, c, t, and g. They claimed that an inscribed sphere of radius r = 1/ contains almost all points corresponding to various genomes, implying that S <r ². The distribution of the points P obtained by S was studied using the Z-curve.     

Unidirectional Movement of Cellulose Synthase Complexes in Arabidopsis Seed Coat Epidermal Cells Deposit Cellulose Involved in Mucilage Extrusion,Adherence, and Ray Formation

CELLULOSE SYNTHASE5 (CESA5) synthesizes cellulose necessary for seed mucilage adherence to seed coat epidermal cells of Arabidopsis (Arabidopsis thaliana). The involvement of additional CESA proteins in this process and details concerning the manner in which cellulose is deposited in the mucilage pocket are unknown. Here, we show that both CESA3 and CESA10 are highly expressed in this cell type at the time of mucilage synthesis and localize to the plasma membrane adjacent to the mucilage pocket. The isoxaben resistant1-1 and isoxaben resistant1-2 mutants affecting CESA3 show defects consistent with altered mucilage cellulose biosynthesis. CESA3 can interact with CESA5 in vitro, and green fluorescent protein-tagged CESA5, CESA3, and CESA10 proteins move in a linear, unidirectional fashion around the cytoplasmic column of the cell, parallel with the surface of the seed, in a pattern similar to that of cortical microtubules. Consistent with this movement, cytological evidence suggests that the mucilage is coiled around the columella and unwinds during mucilage extrusion to form a linear ray. Mutations in CESA5 and CESA3 affect the speed of mucilage extrusion and mucilage adherence. These findings imply that cellulose fibrils are synthesized in an ordered helical array around the columella, providing a distinct structure to the mucilage that is important for both mucilage extrusion and adherence.The epidermal cells of Arabidopsis (Arabidopsis thaliana) seed coats produce two distinct secondary cell walls: pectin-rich mucilage and cellulose-rich columellae (Western et al., 2000). When seeds are hydrated, mucilage expands rapidly, rupturing the outer tangential cell wall and forming a mucilage capsule that surrounds the seed. Seed coat mucilage is composed primarily of rhamnogalacturonan I (RG I) and also contains homogalacturonan (HG), hemicelluloses (such as xylans and glucomannans), and cellulose (for review, see Haughn and Western, 2012). Extruded mucilage consists of an outer, nonadherent fraction and an inner, adherent fraction (Western et al., 2000, 2001; Macquet et al., 2007a). The adherent and nonadherent mucilage layers differ in the amount of methylesterified HG (Rautengarten et al., 2008; Saez-Aguayo et al., 2013; Voiniciuc et al., 2013), galactans (Dean et al., 2007; Macquet et al., 2007b), arabinans (Arsovski et al., 2009), mannans (Yu et al., 2014), and cellulose (Harpaz-Saad et al., 2011; Mendu et al., 2011; Sullivan et al., 2011), all of which influence the physical properties of the layers.Adherent mucilage has a distinct structure, which can be examined using cell wall dyes and antibodies. When treated with cellulose-specific dyes, densely stained rays extend from the top of each columella to the outer edge of the adherent layer, many cell lengths above the seed surface (Mendu et al., 2011; Sullivan et al., 2011). Cytological evidence indicates that cellulose, pectins, and mannans are components of the ray (Haughn and Western, 2012; Griffiths et al., 2014; North et al., 2014; Yu et al., 2014), although the exact manner in which they are assembled is unknown.Cellulose is abundant in mucilage rays and mediates adherence. Loss-of-function mutations in CELLULOSE SYNTHASE5 (CESA5) result in reduced cellulose levels and increased detachment of mucilage from the seed (Harpaz-Saad et al., 2011; Mendu et al., 2011; Sullivan et al., 2011; Griffiths et al., 2014). How a reduction in cellulose results in a loss of adherence is still unknown, but it likely involves interaction with other mucilage components such as pectin and arabinogalactan proteins (Griffiths et al., 2014). Since cesa5 mutants still have some cellulose in the rays of the adherent mucilage halo (Mendu et al., 2011; Sullivan et al., 2011), additional cellulose synthases must be involved in mucilage cellulose biosynthesis.The Arabidopsis genome encodes 10 different CESAs (Delmer, 1999; Richmond and Somerville, 2000). Multiple lines of evidence suggest that three different CESAs are required to form one active cellulose synthase complex (CSC; for review, see Somerville, 2006). CSCs are membrane-bound protein complexes that synthesize cellulose microfibrils in the apoplast (for review, see Somerville, 2006; Endler and Persson, 2011; Lei et al., 2012). CESA1, CESA3, and CESA6 are considered the core components of the primary wall CSC (Desprez et al., 2007; Persson et al., 2007). CESA2, CESA5, and CESA9 are partially redundant to CESA6 in primary wall biosynthesis, and genetic evidence suggests that each of these CESA polypeptides can form a functional CSC with CESA3 and CESA1 (Desprez et al., 2007; Persson et al., 2007). CESA10 is expressed in young plants, stems, floral tissue, and the base of rosette leaves (Beeckman et al., 2002; Doblin et al., 2002), but its function in cellulose biosynthesis is unclear. Other cesa mutant lines have been examined for altered mucilage phenotypes (cesa1, radially swollen1 [Burn et al., 2002; Sullivan et al., 2011], cesa2, cesa6, and cesa9 [Mendu et al., 2011]; CESA3, je5 [Sullivan et al., 2011] and cesa10-1 [Sullivan et al., 2011]); to date, only CESA5 has been shown to be required for cellulose biosynthesis during mucilage deposition.Two mutant alleles of CESA3, isoxaben resistant1-1 (ixr1-1) and ixr1-2, were isolated in a screen for resistance to the herbicide isoxaben (Scheible et al., 2001). Isoxaben inhibits the incorporation of Glc into the emerging cellulose polymer and is considered a potent and specific inhibitor of cellulose biosynthesis (Heim et al., 1990). Homozygous ixr1-1 and ixr1-2 lines show increased resistance to the herbicide, and the mutations causing this resistance were mapped to the genomic locus of CESA3 (Heim et al., 1990; Scheible et al., 2001). The ixr1-1 and ixr1-2 mutations cause amino acid substitutions near the C terminus of the CESA3 protein. ixr1-1 causes a Gly-to-Asn substitution (G998A) located in a transmembrane domain, while ixr1-2 contains a Thr-to-Ile substitution (T942I) in an apoplastic region of the protein between two transmembrane domains (Scheible et al., 2001). Recently, the ixr1-2 allele was shown to affect the velocity of CSCs in the plasma membrane, which consequently modifies cellulose crystallinity in the cell wall (Harris et al., 2012). It is not exactly clear how the ixr1-1 mutation affects cellulose biosynthesis. The effects of either of these mutations on seed coat mucilage have not been investigated.Since mucilage is composed primarily of pectins with smaller amounts of cellulose, seed coat epidermal cells represent an excellent system to study cellulose biosynthesis and interactions between cellulose and other wall components in muro. In this study, we investigated how cellulose is synthesized and deposited in seed coat epidermal cells. We show that at least three different CESA proteins are highly expressed in the seed coat epidermis during mucilage biosynthesis. These CESAs are oriented and move in a linear fashion around the cytoplasmic column of each cell in an identical pattern to cortical microtubules. In addition, we provide evidence that the adherent mucilage has a helical structure that expands and unwinds during extrusion to form the mucilage ray. We propose that during seed coat epidermal cell development, the biosynthesis of cellulose predetermines the structure of rays in the adherent mucilage layer.     

Stichorchis subtriquetrus in European beaver from Croatia: first report

Since 1996, European beavers (Castor fiber) have been reintroduced from Germany to Croatia. However, little is known about the health status of the established population. Necropsy of a 2-year-old female European beaver, which died from a car accident, revealed 28 adult trematodes in the preserved fragments of the colon and in the peritoneal cavity. All of them were identified as Stichorchis subtriquetrus. The typical location of these parasites is the large intestine, and the finding of several specimens in the peritoneal cavity can be explained by the severe trauma. We assume that the trematode was introduced to Croatia along with the beavers. Furthermore, it is very probable that the life cycle can be completed as the intermediate hosts are part of the local fauna. This is the first record of S. subtriquetrus in Croatia.     

The unsaturated acyclic nucleoside analogues bearing a sterically constrained (Z)-4′-benzamido-2′-butenyl moiety: Synthesis,X-ray crystal structure study,cytostatic and antiviral activity evaluations

A series of the novel acyclic unsaturated pyrimidine (1–12) and adenine (13) nucleoside analogues bearing conformationally restricted (Z)-2′-butenyl moiety were synthesized and evaluated for their antiviral and cytostatic activity potency against malignant tumor cell lines and normal human fibroblast (WI38). The N-1 and/or N-3 acyclic side chain substitution in pyrimidine ring in N-3 substituted 5-trifluoromethyluracil derivative (11), N-1, N-3 disubstituted 5-fluorouracil derivative (12) and adenine derivative (13) was deduced from their ¹H and ¹³C NMR spectra and confirmed by single crystal X-ray structure analysis. The X-ray crystal structure analysis 11–13 revealed also supramolecular self-assemblies, in which infinite chains or dimers built two- and three-dimensional networks. The results of the in vitro cytostatic activity evaluations of 1–13 indicate that the majority of the compounds tested exhibited a non-specific and moderate antiproliferative effect at the highest concentration (100 μM). Of all evaluated compounds on the cell lines tested only the N-1 4″-fluoro-substituted-benzamide uracil derivative (7) showed rather marked and selective inhibitory activity against the growth of MCF-7 cells at a concentration of 2.7 μM and no cytotoxic effect on normal fibroblasts WI38. This compound can be therefore considered as a potential antitumor lead compound for further synthetic structure modification.     

Biophysical methods for biochip analysis. Use of wide field digital fluorescent microscopy

This paper discusses an issue on the development of biophysical methods for biochip analysis. A scheme and construction of a biochip analyzer based on wide-field digital fluorescence microscopy are described. The analyzer is designed to register images of biological microchips labeled with fluorescent dyes. The device developed is useful for high-sensitive throughput recording analyses by biochips after interaction of immobilized probes with fluorescently labeled sample molecules as well as it provides the higher rate of the analysis compared to laser scanning devices. With this analyzer a scope where biological microchips can be applied becomes wider, the development of new protocols of the analyses is possible and standard analyses run faster with the use of biochips, the expenses for the analysis performance can be reduced.