” We hypothesized

” We hypothesized ZD1839 purchase that Boston Bowel Preparation Scale (BBPS) scores

could provide a way to standardize the concept of “adequacy”. We performed a retrospective analysis of average-risk screening colonoscopy reports submitted to the Clinical Outcomes Research Initiative (CORI) data repository between 10/2009 and 8/2012. We included only reports documenting a BBPS score and a recommendation for timing of the next colonoscopy and excluded procedures with polyps. We evaluated recommended follow-up intervals stratified by total and segment BBPS scores. We also presented 4 standardized colonoscopy videos with varying degrees of bowel cleanliness to participants of the BBPS Educational Program, a web-based program demonstrating the BBPS, and asked for recommended colonoscopy follow-up intervals. Among 3226 average risk colonoscopies with a BBPS score, 1340 (41.5%) had polyps and 601 (18.6%) lacked follow-up recommendations and were thus excluded. The remaining 1285 procedures, performed by 55 endoscopists, had a median (interquartile range) BBPS score of 8 (7-9). Median recommended follow-up time decreased as BBPS scores decreased, with a sharp drop-off below a BBPS

score of 6 (see Figure). Among reports with total BBPS score of 6 or 7 (n=364), 17 (5%) contained a segment score of 0 or 1 and were associated with shorter median follow-up time compared to reports in which all segment scores were ≥2 (5 vs.10 years, P<.001). Whenever any colonoscopy Clomifene contained a single segment score of 1 (n=55), that segment’s location (right, left, transverse colon) had no impact on recommended follow-up intervals (P=0.955). Video cases were reviewed by 119 endoscopists, including Ruxolitinib cell line 39 CORI users, 51 non-CORI US endoscopists and 29 international endoscopists. Recommended follow-up time decreased as BBPS scores decreased (P<.001; see Table). There was no difference in recommended follow-up time by location

of practice, although more US participants (87%) recommended 10 year follow-up compared to international participants (52%) for Case D (P=.0012). BBPS scores correlate with endoscopist behavior regarding follow-up intervals for colonoscopy. Because BBPS scores have previously been shown to have excellent inter-rater agreement, a total BBPS score ≥ 6 and/or all segment scores ≥ 2 provides a standardized definition of “adequate” when describing bowel cleanliness. Recommended follow-up interval for next colonoscopy for video cases among endoscopists who agreed on the Boston Bowel Preparation Scale score for each case “
“Despite advances in bowel preparation methods, the quality of bowel preparation in patients undergoing colonoscopy remains unsatisfactory. The time point chosen for improvement of education may be important for adequate bowel preparation. To evaluate the effect of telephone re-education on the day before colonoscopy (instead of the day of appointment – regular appointment) on the quality of bowel preparation and colonoscopic findings.

0 and R2013 1, respectively) The POC data product provided by NA

0 and R2013.1, respectively). The POC data product provided by NASA is based on Stramski et al. (2008) algorithm. The full details of the approach used by NASA in standard processing of satellite ocean color data are given at http://oceancolor.gsfc.nasa.gov/.

Spatial resolution of satellite data was about 1.1 km at nadir for the Merged Local Area Coverage (MLAC) SeaWiFS data and 1 km for the Local Area Coverage (LAC) MODIS Aqua data. We also used Global Area Coverage ABT-263 solubility dmso (GAC) SeaWiFS data with effective resolution of about 4.5 km. Satellite POC data have been stored for each pixel containing a coincident in situ data point. Only data pairs with a time difference between in situ measurement and satellite overpass less than 2 h and with a low spatial variability in a 3 × 3 pixel square were used in the analysis. The center pixel in satellite image was the nearest to the in situ measurement. The comparison was carried out if at least 6 of 9 satellite pixels were valid and the average difference between the central pixel and all the other valid pixels was less then 25%. In some cases not one but two overpasses during the same day could have been matched with one in situ measurement. In that case, if both match-ups satisfied the

criteria described above, we have used the one that had the smaller time difference between the satellite and the in situ measurement. These match-up BKM120 clinical trial criteria differ somewhat from those used in Bailey and Werdell (2006). After the compilation of the data using these criteria, the joint satellite and in situ data set included 260 match-ups of POC concentrations. The geographical positions of these data are indicated in Fig. 1. The differences between in situ and satellite-derived POC have been quantified by standard methods (Ostasiewicz et al., 2006): – the absolute

average Hydroxychloroquine price error (AAE) AAE=1N∑i=1N|Oi−Pi| When comparing the in situ and satellite derived POC concentrations one has to remember that both kinds of POC estimates are subject to errors. In-water POC determinations are subject to several potential sources of errors and there is a continued need for further improvement in the methodology. This issue has been discussed in-depth in Gardner et al. (2003) The causes for the overestimation of POC include potential adsorption of dissolved organic carbon (DOC) onto filters during filtration and contamination of samples during handling. Underestimation of POC can result, for example, from an undersampling of the infrequent large particles, settling of particles below the bottle spigots (Gardner, 1977) or incomplete retention of particles on filters. Therefore the true accuracy of in situ POC determinations remains unspecified. For brevity, in this paper, we refer to in-water POC estimates as ‘measured’ and to the differences between satellite-derived and in-water POC estimates as ‘errors’.

This permits the analysis of more defined antigen specific respon

This permits the analysis of more defined antigen specific responses while reducing the requirement to handle live influenza virus in the laboratory. We have developed a method to potentiate the detection and analysis

of influenza antigen specific T cells utilizing infected selleckchem CKC to present viral peptides in a manner biologically relevant to CD8 T cells. We have demonstrated that our co-culture ELISpot detects greater numbers of antigen specific CD8 T cells than ELISpot with whole virus as an antigen. Our assay can also be adapted to use recombinant viruses to infect CKC, increasing its specificity and reducing the requirement to work with live influenza virus. Our results are the first to demonstrate detection by flow cytometry of influenza-specific IFNγ responses in individual T cells from LPAI infected birds. The ability of our method to detect such large numbers of antigen specific T cells (similar numbers to positive controls with PMA/ionomycin, see example Supplementary Olaparib purchase Fig. 5) likely reflects not only the high promiscuity of the B21 haplotype, but also the fact that our CKC cell line expresses only MHC

class I and presents peptides following a biologically relevant infection process. In ELISpot using whole influenza virus we were able to detect antigen specific responses, although these were much lower (Fig. 1). Although ELISpot has previously been used to measure antiviral responses against other avian viruses, including NDV (Ariaans et al., 2008) and IBV (Ariaans et al., 2009), it has never been employed to analyze avian

responses to influenza. In the present study, Lck following challenge with H7N7 LPAI, the birds became serologically positive and showed specific IFNγ responses, irrespective of whether inactivated or live avian influenza virus was added to endogenous APCs (Fig. 1). Additionally, ELISpot with live virus added to splenocytes from infected birds further reduced SFU counts. It is possible that live virus affects the interactions, and/or the functionality, of cells in vitro (Hinshaw et al., 1994, Oh et al., 2000 and Hao et al., 2008). It was interesting to note that splenocytes from infected birds have greater SFU responses to PMA in our study. PMA does not activate all T cells (Suzawa et al., 1984 and Kim et al., 1986)., It may be that antigen experienced cells (from infected birds) have a lower threshold of activation and are activated more readily by PMA, hence the higher SFU counts in the infected cohort positive control compared with the non-infected. Another possibility is altered lymphocyte subset frequencies in infected birds.

Another strength of the study is that LSI, liver fat Apo A-I and

Another strength of the study is that LSI, liver fat. Apo A-I and R2* increased in parallel showing an internal consistency of the observations. An obvious limitation of the present study is that only female rats were investigated.

As BPA is an estrogenic-acting compound it cannot be taken for granted that different effects would not be seen in males. Unfortunately, we do not have reproducibility data on the methods used in the paper. No detailed histopathological examinations of the livers were performed. The study was performed during 10 weeks of exposure. A longer exposure period might result in effects on the obesity measures used. In the present study we found no evidence that BPA exposure affects fat mass in fructose-fed juvenile Fischer 344 rats. We also suggest that the increase in liver fat infiltration

and apo A-I may result from combination Talazoparib datasheet effects of fructose and BPA exposure, and eventually may lead to more severe metabolic consequences. The present findings would motivate future studies regarding these more long term metabolic consequences. If so, the finding Sirolimus price that fructose fed rats exposed to BPA induced fat infiltration in the liver at dosages close to the current TDI might be of concern given the widespread use of this compound in our environment and since a great proportion of the human population is exposed to both BPA and fructose daily. None declared. We thank Raili Engdahl for excellent Carnitine dehydrogenase technical assistance, Katarina Cvek for expert advice about animal experiments, and Martin Ahlström for assistance with the MR image segmentation and Erik Lampa for statistical support. “
“Carcinogenicity studies have demonstrated that long-term exposure to various respirable micro- and

nanoscale particles (MNP) can induce lung tumors, in particular in the rat model (Saffiotti and Stinson, 1988, Wiessner et al., 1989, Donaldson and Borm, 1998, Muhle et al., 1989, Nikula, 2000 and Roller, 2009). Especially the surface characteristics of poorly soluble particles predominantly determine the carcinogenic potential of MNP (Oberdörster et al., 2005 and Duffin et al., 2007), as they do not act as single molecules, but more likely in a physico-mechanical or physico-chemical way. Different genotoxic modes of action could explain the carcinogenic potential of particles in the lung in non-overload and overload situations. Possible genotoxic mechanisms of MNP in vivo, as summarized earlier by Knaapen et al. (2004), seem to comprise indirect (secondary) mechanisms that are phagocytosis- and/or inflammation-driven, but also directly particle-related (primary) genotoxic modes of action. Release of reactive oxygen (ROS) and nitrogen (RNS) species either by (i) oxidative burst of phagocytes, (ii) disturbance of the respiratory chain, (iii) activation of ROS-/RNS-producing enzyme systems, or (iv) reactive particle surfaces with subsequent oxidative DNA damage is thought to be of principal importance.

All rights reserved http://dx doi org/10 1016/j gde 2012 12 009

All rights reserved. http://dx.doi.org/10.1016/j.gde.2012.12.009 Genomes employ remarkably diverse architectures to store information in DNA sequences and direct all forms of biological function across the tree of life. Information is stored concisely and directly at most bacterial species genomes, where genome evolution favors concise organization and functional specialization. As organisms’ complexity increase, and in particular in multi-cellular eukaryotes, genomes are expanding mildly in terms of new genes, but scale up by two to three orders of

magnitudes in size from millions AP24534 datasheet to billions of bases. Genetic information is then embedded into long and complex DNA sequences in a redundant and indirect fashion. Although the implications of such sparse encoding are widely believed to be profound, it was so far difficult to describe them precisely. Mechanisms that are capable or processing and possibly taking advantage of fragmented and patchy genomic encodings (e.g. RNA splicing) promote the notion that genome sequences are heterogeneous in their information content, ranging from perfectly optimized

elements similar those making up bacterial genomes to ‘junk’-like sequences spanning millions of bases with seemingly no direct function. In contrast, numerous recent studies are utilizing high throughput sequencing to generate rich maps of genomic and epigenomic activity, suggesting that much of the genome Farnesyltransferase is under selection [1 and 2] and involved in gene regulation. Ultimately, understanding Selleckchem Selisistat genome function, and describing how and why metazoan genomes are so large, complex and redundant, must be achieved through physical characterization of genome and chromosome structure. In this short review we survey recent technological

and analytical advances leading to new insight into the structure of complex chromosomes. By mapping chromosomal contacts, we propose, geneticists and epigeneticists are finding vital clues that may lead to integrative, physical and mechanistic models of genome function. Historically, the study of chromosomal architectures relied on structural and biochemical studies of nucleosomes and their modifications at the local level (reviewed in [3]) and on fluorescence-based microcopy (reviewed in [4]) for studying longer range contacts and global chromosomal organization. The development of chromosome conformation capture [5] by Dekker and others and the combination of 3C with powerful genomics approaches [6••, 7••, 8••, 9, 10 and 11] facilitated the quantification of chromatin contacts at unprecedented scale and breadth. 3C is performed through fragmentation of the genome (using, e.g. sequence specific restriction enzymes) followed by re-ligation of DNA fragments that were crosslinked together, owing to physical proximity at the time of nuclei fixation.

All DNA that was sub-diploid in size (sub-G1) was considered to b

All DNA that was sub-diploid in size (sub-G1) was considered to be caused by internucleosomal DNA fragmentation. Table 1 indicates the cell cycle distribution obtained. After a 12-h incubation, the ATZD treated with AC-4, AC-7 and AC-10 (2.5 μg/ml) caused a small PLX4032 supplier increase in the number of cells in the G2/M phase compared with the negative control (15.7%, p < 0.05). For the ATZD-treated cells, the percentage of cells in the G2/M phase were 19.7%, 19.2% and 19.9%, for AC-4, AC-7 and AC-10, respectively. After a 24-h incubation, the cells in the G0/G1 and S phases remained mostly

unchanged; however, there were fewer cells in the G2/M phase. Additionally, all ATZD caused significant internucleosomal DNA fragmentation at all of the concentrations tested (p < 0.05), which implies that ATZD preferentially caused cells from the G2/M phase to transition into sub-G1. Cells treated with m-AMSA served as the positive control, and had an increased number of cells in the G2/M interval and a

significant amount of internucleosomal DNA fragmentation. After 12- and 24-h incubations, the effects of ATZD were evaluated based on cell morphology using hematoxylin–eosin and acridine orange/ethidium bromide staining. The integrity of the cell membrane and selleck the mitochondrial membrane potential were also determined by flow cytometry. Additionally, after a 24-h incubation, phosphatidylserine externalisation and caspase 3/7 activation were measured by flow cytometry. After a 12-h incubation, HCT-8 cells either treated or untreated with ATZD, were tested at all concentrations and presented

slight morphological changes (data not shown). On the other hand, after a 24-h incubation, morphological examination of HCT-8 cells showed severe drug-mediated changes. The hematoxylin–eosin stained HCT-8 cells treated with ATZD presented a morphology consistent with apoptosis, including a reduction in cell volume, chromatin condensation and nuclei fragmentation (Fig. 4). The acridine orange/ethidium bromide stained and treated cells also displayed a morphology consistent with apoptosis, in a time- and concentration-dependent manner (p < 0.05, Fig. 5). m-AMSA, served as the positive control, which also induced morphological changes consistent with apoptosis. The integrity of the cell Resveratrol membrane is a parameter of cell viability that differs between apoptotic and necrotic cells. After 12- or 24-h of exposure, ATZD induced a slight disruption in the plasmatic membrane, which was only observed at the higher concentrations tested (Figs. 6A, B). As cited above, the internucleosomal DNA fragmentation was markedly increased in ATZD-treated cells (p < 0.05, Table 1). Both of these modifications are characteristics of apoptotic cells. In addition, ATZD induced mitochondrial depolarisation in a time- and concentration-dependent manner (p < 0.05, Figs. 6C, D).

One of the known literature formulas for estimating Chl

One of the known literature formulas for estimating Chl AZD0530 ic50 a obtained for the Baltic Sea environment is the one given by Siegel et al. (1994). It uses the green-to-red reflectance

ratio (but at wavelengths slightly shifted compared to the wavelengths already analysed in this work) and takes the following form: Chl a = 31.05(Rrs (510)/Rrs(670))− 2.115. If we used the modelled reflectance spectra obtained in this work, the equivalent formula would take the form Chl a = 32.3(Rrs(510)/Rrs(670))− 1.24 (n = 82; r2 = 0.7; X = 1.54). As can be seen in Figure 10a, these two last formulas would agree only in the ranges of the relatively low values of the Rrs(510)/Rrs(670) ratio (which corresponds to Chl a concentrations

of the order of 10 mg m− 3 and higher). For high values of that green-to-red reflectance ratio, the latter formula would predict Chl a values several times higher than the one given by Siegel et al. (1994). www.selleckchem.com/products/PLX-4032.html The other formula known from the literature is the one from the paper by Darecki et al. (2005). It uses the green-to-orange ratio of Rrs(550)/Rrs(590) and after simple transformation takes the form Chl a = 5.47 (Rrs(550)/Rrs(590))− 4.681. Based on the modelling results obtained in the present work, the equivalent formula using the same reflectance ratio would be Chl a = 30 (Rrs(550)/Rrs(590))− 3.33 (n = 82; r2 = 0.76;

Y-27632 solubility dmso X = 1.48). Figure 10b shows that these last two formulas would exhibit distinct differences. Both formulas are relatively steep functions of the green-to-orange reflectance ratio but for the same values of this, the predicted ranges of Chl a would differ by about one order of magnitude. However, in view of the results of the latter comparison, it has to be emphasised that the 590 nm reflectance band taken for that additional test lies relatively far from the modelling input data on the light absorption coefficient an(λ) (we recall that the nearest an input data bands were at 555 and 650 nm). As a consequence, the modelled values of Rrs at 590 nm band should be treated with a relatively low level of confidence. Nevertheless, the last two additional quantitative comparisons of the relationships between Chl a and different colour ratios should warn the potential user that all the results of the simplified modelling performed here, and in effect, all the semi-empirical (reflectance-based) formulas presented in this work, should be treated as qualitative rather than quantitative. Finally, let us comment on the comparison of all the statistical parameters obtained here for different variants of both empirical (see Table 1 and Table 2) and semi-empirical formulas (see Table 3 and Table 4).

In addition to the presentation of IOP-based relationships for th

In addition to the presentation of IOP-based relationships for the two satellite light wavelengths of 443 and 555 nm, the statistical analyses are supplemented with examples of analogous relationships but determined at the optimal bands chosen from among those original light wavelengths for which the HydroScat-4 and AC-9 instruments performed in situ measurements. To derive statistical formulas for biogeochemical properties of suspended matter as functions of remote-sensing reflectance values, the available dataset has to be extended

with the aid of radiative transfer modelling. It has been common practice in much optical modelling work that the average values of the constituent-specific optical coefficient multiplied by the assumed selleck products concentrations of these constituents learn more give the modelled absolute values of these optical coefficients, which are then used as further inputs for radiative transfer modelling. But because the very large variability of constituent-specific optical coefficients of suspended matter in the southern Baltic Sea were documented in an earlier work by S. B. Woźniak et al. (2011), it was decided not

to use averaged values as the modelling input. Instead, a different approach to the problem is taken: in each separate modelling case the real, measured optical coefficients (i.e. the values of the coefficients an(λ), cn(λ) and bbp(λ)) are used as modelling input and the corresponding and actually measured values of biogeochemical properties are also used in the subsequent statistical analyses. From the available empirical material a subset of 83 cases was selected (see the stations denoted by grey dots in Figure 2), which

consists of only those cases for which all the biogeochemical properties of the relevant particulate matter (i.e. concentrations of SPM, POM, Cisplatin cost POC and Chl a) and all the seawater IOPs (i.e. values of an(λ), cn(λ) and bb(λ)) required for further modelling were measured at the same time. For this particular data subset, the hypothetical spectra of the remote-sensing reflectance Rrs [sr− 1] were then determined on the basis of radiative transfer numerical simulations. The Hydrolight-Ecolight 5.0 (Sequoia Scientific, Inc.) model was applied with a set of simplifying assumptions. The modelled hypothetical water bodies were chosen to be infinitely deep, and all the IOPs of the modelled waters were chosen to be constant with depth. This assumption is obviously a significant simplification, but it most likely represents quite well a common situation in the Baltic Sea, where the relatively shallow subsurface layer of water penetrated by sunlight is mixed as a result of wave action and turbulence caused by surface wind stress. Another simplification was the assumption that no inelastic scattering (no Raman scattering, or chlorophyll or CDOM fluorescence) and no internal sources (no bioluminescence) were taken into account.

Using the simulation parameters in Table 1, the linear stability

Using the simulation parameters in Table 1, the linear stability analysis was insensitive to setting νvνv to this smaller value, selleckchem so for the purpose of this modeling exercise the smaller viscosity/diffusivity sufficed. One consequence of varying N2N2 and M2M2 is that the dynamics may become sensitive to whether the hydrostatic approximation is employed. Because the balanced Richardson number can be tuned by adjusting

the values of M2,N2M2,N2, and f  , the individual parameters for each set are chosen to fix the hydrostatic parameter ( Marshall et al., 1997) equation(25) η=γ2Ri,where γ=h/Lγ=h/L is the aspect ratio of the motion. For η≪1η≪1 it is appropriate to use the hydrostatic approximation to the vertical momentum equation. The parameter γγ is estimated according to the initial M2M2 and N2N2 from the simulations. Because the unstable modes lie in an arc symmetric about the isopycnal, the mean aspect ratio of the motions can be taken as γ=M2/N2γ=M2/N2, and simple algebra gives equation(26) η=f2N21Ri2.The parameter choices in Table 1 are chosen so that η=0.1η=0.1 for the “hydrostatic” parameters and η=10η=10 for the “nonhydrostatic” parameters. Note that in both cases, the fully nonhydrostatic equations are solved. To check whether the results are sensitive to whether a model is run in hydrostatic mode, a parallel this website set of the η=0.1η=0.1 simulations was

run using the MITgcm (Marshall et al., 1997) in hydrostatic mode and with identical initial conditions. The hydrostatic MITgcm gave nearly identical results (not shown) as long as the grid spacing ΔxΔx was less than half the wavelength of the most unstable mode; when ΔxΔx was set above this threshold the MITgcm was prone to numerical instability which eventually led to the simulation crashing. This numerical instability influenced 4��8C the choice to use the nonhydrostatic solver for these simulations over the MITgcm. Nonetheless, previous work by Mahadevan (2006) suggests that the average vertical fluxes at the length scales in these simulations should be similar regardless of whether the model is run hydrostatically or nonhydrostatically, so it is likely that the results from

the nonhydrostatic solver are robust for the η=0.1η=0.1 simulations at all resolutions. The simulation parameters in Table 1 were chosen specifically to demonstrate cases of grid-arrested restratification (Sets A and C) and completed restratification (B and D) by varying νhνh. The amount of restratification that takes place is not uniquely dependent on the parameter choices in each set; all of the parameters can be varied in relation to one another to change the anticipated final value of Ri  . Fig. 4 shows the growth rate plots for each parameter set. In each case the horizontal viscosity damps the highest wavenumber modes, so that increasing the resolution beyond a certain point does not permit extra modes to become resolved or further restratification to occur.

In contrast, G8, G6, and G9 were among the genotypes with the low

In contrast, G8, G6, and G9 were among the genotypes with the lowest stability and with higher (G8) and lower (G6 and G9) mean yield performances than the overall mean. The yield, stability and yield–stability ranks for 20 tested genotypes in 24 environments based on each of the statistical methods mentioned above are given in Table 3. Comparison of the statistical methods PI3K targets based on the yield ranks showed that the methods generally gave similar results in the ranking of genotypes. For example, the five top-ranked

genotypes based on AMMI were G4, followed by G10, G19, G1, and G17; based on the GGE biplot were G4 followed by G10, G1, G19, and G17; based on JRA were G8, G4, G1 = G12, and G10; and based on the YSi statistic were G4, G10, G19, G1, and G17. With respect to stability ranks, genotypes G2, G15, G12, G11, and G17 were found to be stable based http://www.selleckchem.com/products/r428.html on AMMI distance, whereas the five top-ranked genotypes based on the GGE biplot were G18 = G12 = G2, G14, and G3, showing that AMMI and the GGE biplot gave similar results in identifying two of the five top-ranking genotypes as stable. According to JRA the most desirable genotypes based on stability ranks were G2, G17, G10,

G16, and G3, and based on the YSi statistic the most stable genotypes were G2, G17, G3, G16, and G18. Similar stability ranks were assigned by the JRA method and YSi statistic, as they identified four of the five top-ranking genotypes as stable. For yield–stability, the AMMI analysis identified G10 followed by G17, G3, G15, and G12 as the top-ranking high-yielding and stable genotypes; whereas G18 followed by G17 = G12 and G4 = G10 were characterized

by the GGE biplot as high-yielding and stable. According to JRA, the top-ranking high-yielding and stable genotypes were G10, followed by G4, G12, G17, and G3, and based on the almost YSi statistic the highest-ranking genotypes were G4 = G10, G17, G19, and G1 = G18. All four methods identified G10 and G17 as among the five top-ranking high-yielding and stable genotypes. Significant rank correlations were found between the statistical methods in the ranking of genotypes for yield, stability and yield–stability (Table 4). With respect to yield, the statistical methods were significantly correlated (P < 0.01) in the ranking of genotypes. The correlations varied from 0.72 (JRA–YSi; P < 0.01) to 0.99 (GGE–AMMI; P < 0.01) indicating that AMMI and the GGE biplot agreed most closely in ranking genotypes for yield. The statistical methods were positively correlated in identifying stable genotypes. The Spearman’s rank correlations for stability indices ranged from 0.53 (GGE–YSi; P < 0.05) to 0.97 (JRA–YSi). The AMMI distance (AMMID) was highly correlated with the stability indices in JRA (r = 0.83; P < 0.01) and YSi (r = 0.86; P < 0.01).