Dalbergia odorifera undergoes massive molecular shifts in response to waterlogging combined with salinity.

Field and greenhouse studies attempting to describe the molecular responses of plant species under waterlogging (WL) combined with salinity (ST) are almost nonexistent. We integrated transcriptional, metabolic, and physiological responses involving several crucial transcripts and common differentially expressed genes and metabolites in fragrant rosewood (Dalbergia odorifera) leaflets to dissect plant-specific molecular responses and patterns under WL combined with ST (SWL). We discovered that the synergistic pattern of the transcriptional response of fragrant rosewood under SWL was exclusively characterized by the number of regulated transcripts. The response patterns under SWL based on transcriptome and metabolome regulation statuses revealed different patterns (additive, dominant, neutral, minor, unilateral, and antagonistic) of transcripts or metabolites that were commonly regulated or expressed uniquely under SWL. Under SWL, the synergistic transcriptional response of several functional gene subsets was positively associated with several metabolomic and physiological responses related to the shutdown of the photosynthetic apparatus and the extensive degradation of starch into saccharides through α-amylase, ß-amylase, and α-glucosidase or plastoglobuli accumulation. The dissimilarity between the regulation status and number of transcripts in plants under combined stresses led to nonsynergistic responses in several physiological and phytohormonal traits. As inferred from the impressive synergistic transcriptional response to morpho-physiological changes, combined stresses exhibited a gradually decreasing effect on the changes observed at the molecular level compared to those in the morphological one. Here, by characterizing the molecular responses and patterns of plant species under SWL, our study considerably improves our understanding of the molecular mechanisms underlying combined stress.

Physio-biochemical and metabolomic responses of the woody plant Dalbergia odorifera to salinity and waterlogging.

BACKGROUND: Trees have developed a broad spectrum of molecular mechanisms to counteract oxidative stress. Secondary metabolites via phenolic compounds emblematized the hidden bridge among plant kingdom, human health, and oxidative stress. Although studies have demonstrated that abiotic stresses can increase the production of medicinal compounds in plants, research comparing the efficiency of these stresses still needs to be explored. Thus, the present research paper provided an exhaustive comparative metabolomic study in Dalbergia odorifera under salinity (ST) and waterlogging (WL). RESULTS: High ST reduced D. odorifera's fresh biomass compared to WL. While WL only slightly affected leaf and vein size, ST had a significant negative impact. ST also caused more significant damage to water status and leaflet anatomy than WL. As a result, WL-treated seedlings exhibited better photosynthesis and an up-regulation of nonenzymatic pathways involved in scavenging reactive oxygen species. The metabolomic and physiological responses of D. odorifera under WL and salinity ST stress revealed an accumulation of secondary metabolites by the less aggressive stress (WL) to counterbalance the oxidative stress. Under WL, more metabolites were more regulated compared to ST. ST significantly altered the metabolite profile in D. odorifera leaflets, indicating its sensitivity to salinity. WL synthesized more metabolites involved in phenylpropanoid, flavone, flavonol, flavonoid, and isoflavonoid pathways than ST. Moreover, the down-regulation of L-phenylalanine correlated with increased p-coumarate, caffeate, and ferulate associated with better cell homeostasis and leaf anatomical indexes under WL. CONCLUSIONS: From a pharmacological and medicinal perspective, WL improved larger phenolics with therapeutic values compared to ST. Therefore, the data showed evidence of the crucial role of medical tree species' adaptability on ROS detoxification under environmental stresses that led to a significant accumulation of secondary metabolites with therapeutic value.

Bayesian Factor-adjusted Sparse Regression.

Many sparse regression methods are based on the assumption that covariates are weakly correlated, which unfortunately do not hold in many economic and financial datasets. To address this challenge, we model the strongly-correlated covariates by a factor structure: strong correlations among covariates are explained by common factors and the remaining variations are interpreted as idiosyncratic components. We then propose a factor-adjusted sparse regression model with both common factors and idiosyncratic components as decorrelated covariates and develop a semi-Bayesian method. Parameter estimation rate-optimality and model selection consistency are established by non-asymptotic analyses. We show on simulated data that the semi-Bayesian method outperforms its Lasso analogue, manifests insensitivity to the overestimates of the number of common factors, pays a negligible price when covariates are not correlated, scales up well with increasing sample size, dimensionality and sparsity, and converges fast to the equilibrium of the posterior distribution. Numerical results on a real dataset of U.S. bond risk premia and macroeconomic indicators also lend strong supports to the proposed method.

Adaptive Huber Regression on Markov-dependent Data.

High-dimensional linear regression has been intensively studied in the community of statistics in the last two decades. For the convenience of theoretical analyses, classical methods usually assume independent observations and sub-Gaussian-tailed errors. However, neither of them hold in many real high-dimensional time-series data. Recently [Sun, Zhou, Fan, 2019, J. Amer. Stat. Assoc., in press] proposed Adaptive Huber Regression (AHR) to address the issue of heavy-tailed errors. They discover that the robustification parameter of the Huber loss should adapt to the sample size, the dimensionality, and the moments of the heavy-tailed errors. We progress in a vertical direction and justify AHR on dependent observations. Specifically, we consider an important dependence structure - Markov dependence. Our results show that the Markov dependence impacts on the adaption of the robustification parameter and the estimation of regression coefficients in the way that the sample size should be discounted by a factor depending on the spectral gap of the underlying Markov chain.

Simultaneous dimension reduction and adjustment for confounding variation.

Dimension reduction methods are commonly applied to high-throughput biological datasets. However, the results can be hindered by confounding factors, either biological or technical in origin. In this study, we extend principal component analysis (PCA) to propose AC-PCA for simultaneous dimension reduction and adjustment for confounding (AC) variation. We show that AC-PCA can adjust for (i) variations across individual donors present in a human brain exon array dataset and (ii) variations of different species in a model organism ENCODE RNA sequencing dataset. Our approach is able to recover the anatomical structure of neocortical regions and to capture the shared variation among species during embryonic development. For gene selection purposes, we extend AC-PCA with sparsity constraints and propose and implement an efficient algorithm. The methods developed in this paper can also be applied to more general settings. The R package and MATLAB source code are available at https://github.com/linzx06/AC-PCA.

Accommodative responses under different stimulus conditions.

PURPOSE: The method and the sequence of accommodative stimuli are two factors affecting the steady-state accommodative response. This study investigated the possible effects of these two factors on the accommodative stimulus response curves (ASRCs) and compared three indexes to evaluate ASRC in the different stimulus conditions. METHODS: In this study, we investigated the difference in the accommodative stimulus response curves under six stimulus conditions. The ASRCs of 10 emmetropes were measured when the accommodative stimulus (AS) was induced using three methods: negative lenses, positive lenses, and a Badal stimulator; and the stimulus was presented using two sequences: ascending and descending sequences. Three indexes, accommodative error area (AEA), accommodative error index (AEI), and the slope of ASRC, were used for comparison. RESULTS: The accommodative error was larger when accommodative stimulus was presented in the ascending sequence than in the descending sequence (p < 0.05) in the Badal stimulator method whenever AEI, AEA, or the slope of ASRC was used. The ASRCs elicited with positive lenses were more accurate than that elicited with negative lenses or the Badal stimulator when AEA and AEI were used. However, when the slope of the ASRC was used, the result was not significant. CONCLUSIONS: The results suggest that the accommodative response in the steady state may be affected by the method used for inducing AS and the sequence of AS presented. The AEA is suggested to be used in the evaluation of the accuracy of the ASRC in future studies.

Comparison of metagenomic samples using sequence signatures.

BACKGROUND: Sequence signatures, as defined by the frequencies of k-tuples (or k-mers, k-grams), have been used extensively to compare genomic sequences of individual organisms, to identify cis-regulatory modules, and to study the evolution of regulatory sequences. Recently many next-generation sequencing (NGS) read data sets of metagenomic samples from a variety of different environments have been generated. The assembly of these reads can be difficult and analysis methods based on mapping reads to genes or pathways are also restricted by the availability and completeness of existing databases. Sequence-signature-based methods, however, do not need the complete genomes or existing databases and thus, can potentially be very useful for the comparison of metagenomic samples using NGS read data. Still, the applications of sequence signature methods for the comparison of metagenomic samples have not been well studied. RESULTS: We studied several dissimilarity measures, including d2, d2* and d2S recently developed from our group, a measure (hereinafter noted as Hao) used in CVTree developed from Hao's group (Qi et al., 2004), measures based on relative di-, tri-, and tetra-nucleotide frequencies as in Willner et al. (2009), as well as standard lp measures between the frequency vectors, for the comparison of metagenomic samples using sequence signatures. We compared their performance using a series of extensive simulations and three real next-generation sequencing (NGS) metagenomic datasets: 39 fecal samples from 33 mammalian host species, 56 marine samples across the world, and 13 fecal samples from human individuals. Results showed that the dissimilarity measure d2S can achieve superior performance when comparing metagenomic samples by clustering them into different groups as well as recovering environmental gradients affecting microbial samples. New insights into the environmental factors affecting microbial compositions in metagenomic samples are obtained through the analyses. Our results show that sequence signatures of the mammalian gut are closely associated with diet and gut physiology of the mammals, and that sequence signatures of marine communities are closely related to location and temperature. CONCLUSIONS: Sequence signatures can successfully reveal major group and gradient relationships among metagenomic samples from NGS reads without alignment to reference databases. The d2S dissimilarity measure is a good choice in all application scenarios. The optimal choice of tuple size depends on sequencing depth, but it is quite robust within a range of choices for moderate sequencing depths.

Effects of palm oil blended with oxidized fish oil on growth performances, hematology, and several immune parameters in juvenile Japanese sea bass, Lateolabrax japonicas.

A 60-day feeding trial was conducted to determine the effects of palm oil blended with oxidized and non-oxidized fish oil on growth performances, hematology, and non-specific immune response in juvenile Japanese sea bass, Lateolabrax japonicas. Japanese sea bass (1.73 ± 0.01 g) were fed seven experimental diets containing 100 g/kg of dietary lipid in forms of palm oil (10P), fish oil (10F), fish oil blended with palm oil at different ratios, 6:4 (6F4P) and 4:6 (4F6P), oxidized fish oil (10OF), and oxidized fish oil blended with palm oil at different ratios, 6:4 (6OF4P) and 4:6 (4OF6P). After the feeding trial, the following results were illustrated. No significant effects were observed in survival, feed conversion ratio, condition factor, and hematocrit after feeding with experimental diets for 60 days. The relatively higher specific growth rate and hematology were observed in 6F4P. Furthermore, both palm oil and oxidized fish oil acted as a negatively on serum lysozyme activity (P < 0.05). This study suggested that a ration of 6F4P is recommended as an innocuous ratio for Japanese sea bass. Furthermore, according to the present investigation, palm oil seems to have the ability to improve the protein efficiency when added to oxidized fish diets as well as a positive trend to the growth performance (P > 0.05).

Terahertz Rectangular Waveguides by UV-LIGA with Megasonic Agitation.

This paper researches the fabrication of a WR2.8 terahertz rectangular waveguide operating at the frequency ranging from 260 GHz to 400 GHz via UV-LIGA technology (UV-lithography, electroplating, and molding). In the process, megasonic agitation is applied to improve the mechanical properties and internal surface roughness of the WR2.8 rectangular waveguide. The effects of process parameters on the properties of structures are discussed, and optimized parameters are obtained to achieve accurate geometry dimensions. In addition, the highly crosslinked SU-8 is reliably removed from structures without damage through a synthesis method. The accuracy of the height and width of the WR2.8 rectangular waveguide is 5 µm and 2 µm, respectively, and the measured internal surface roughness is 79.6 nm. Moreover, experimental measurements and numerical simulations of the waveguide are conducted, and the difference between the two is also highlighted.

Granulin epithelin precursor: a bone morphogenic protein 2-inducible growth factor that activates Erk1/2 signaling and JunB transcription factor in chondrogenesis.

Granulin epithelin precursor (GEP) has been implicated in development, tissue regeneration, tumorigenesis, and inflammation. Herein we report that GEP stimulates chondrocyte differentiation from mesenchymal stem cells in vitro and endochondral ossification ex vivo, and GEP-knockdown mice display skeleton defects. Similar to bone morphogenic protein (BMP) 2, application of the recombinant GEP accelerates rabbit cartilage repair in vivo. GEP is a key downstream molecule of BMP2, and it is required for BMP2-mediated chondrocyte differentiation. We also show that GEP activates chondrocyte differentiation through Erk1/2 signaling and that JunB transcription factor is one of key downstream molecules of GEP in chondrocyte differentiation. Collectively, these findings reveal a novel critical role of GEP growth factor in chondrocyte differentiation and the molecular events both in vivo and in vitro.

The effects of spherical aberration on static accommodative responses in emmetropes and myopes.

PURPOSE: To compare the effects of spherical aberration (SA) on accommodative responses between emmetropes (EMs), stable myopes (SMs), and progressing myopes (PMs). METHODS: The accommodative responses (AR) induced by negative lenses from 0 to -5 D in 1 D steps were measured with a Complete Ophthalmic Analysis System (COAS) for 10 emmetropes (EMs), 13 stable myopes (SMs), and 9 progressing myopes (PMs). Three types of AR were calculated from wavefront aberration data based on ocular foci within a paraxial area, natural pupil area, and a 3.5 mm-diameter annular pupil area. The accommodative stimulus response curves (ASRCs) were compared between the three groups and between the three AR conditions. RESULTS: The difference in the slope of ASRC was significant between the three AR conditions (p < 0.0001) but it was not significant between the three refractive groups. In each group, the slope of ASRC for the accommodative response measured within the paraxial area of the pupil (paraxial AR) was significantly steeper than that measured for the natural pupil (total AR) and also for that estimated for the 3.5 mm-diameter annular pupil (p < 0.01). The difference in the y-intercept of ASRC between the three AR conditions was also significant (p = 0.02) but not between the three refractive groups. The y-intercept of ASRC for paraxial AR was more positive than that for total AR in SMs (p < 0.01) and in PMs (p < 0.05); the y-intercept of ASRC for the 3.5 mm peripheral AR was more positive than that for total AR in SMs (p < 0.05). CONCLUSIONS: The spherical aberration of the eye has a significant effect on the slope of ASRC and is less in PMs than in EMs and SMs.

Hoeffding's inequality for general Markov chains with its applications to statistical learning.

This paper establishes Hoeffding's lemma and inequality for bounded functions of general-state-space and not necessarily reversible Markov chains. The sharpness of these results is characterized by the optimality of the ratio between variance proxies in the Markov-dependent and independent settings. The boundedness of functions is shown necessary for such results to hold in general. To showcase the usefulness of the new results, we apply them for non-asymptotic analyses of MCMC estimation, respondent-driven sampling and high-dimensional covariance matrix estimation on time series data with a Markovian nature. In addition to statistical problems, we also apply them to study the time-discounted rewards in econometric models and the multi-armed bandit problem with Markovian rewards arising from the field of machine learning.

Hard thresholding regression.

In this paper, we propose the hard thresholding regression (HTR) for estimating high-dimensional sparse linear regression models. HTR uses a two-stage convex algorithm to approximate the ℓ 0-penalized regression: The first stage calculates a coarse initial estimator, and the second stage identifies the oracle estimator by borrowing information from the first one. Theoretically, the HTR estimator achieves the strong oracle property over a wide range of regularization parameters. Numerical examples and a real data example lend further support to our proposed methodology.

[Significance of Bone Marrow Microvessel Density and Vascular- Related Factors in Multiple Myeloma].

OBJECTIVE: To investigate the clinical significance of bone marrow microvessel density(MVD) and angiogenesis related factors in multipic myeloma(MM). METHODS: Twenty cases of MM and 20 cases of simple fracture were selected and enrolled in MM group and control group respectively. The clinical data and results of laboratorial tests were collected; the bone marrow MVD of patients was detected by using the modified plastic-embedded pathologic sections of bone marrow tissue and histochemistry staining, the expression levels of amgiogenesis-related factors including VEGF, TNF-α, HGF, TGF-α, TGF-ß1, bFGF, Ang-Ⅰ, Ang-Ⅱ in bone marrow supernatant were detected by ELISA; the mRNA expression levels of above-mentioned cytokines in bone marrow mononuclear cells were detected by real time-PCR; the pearson correlation analysis was used to analyze the correlation of MVD with VEGF, HGF and bFGF levels. RESULTS: The MVD in MM group was significantly higher than that in control group (P＜0.001); the mRNA expression of VEGF, TGF-α, TGF-ß1 and HGF in bone marrow mononuclear cells of MM group was higher than that of control group(P＜0.001); the levels of VEGF, HGF, bFGF and THF-α in bone marrow supernatant of MM group were higher than those in control group(P＜0.05), moreover, the MVD positively correlated with levels of VEGF, HGF and bFGF in bone marrow(r=0.488, 0.472 and 0.457). CONCLUSION: The MVD and levels vessel-related factors in bone marrow supernatant of MM patients increase, among which the levels of VEGF and HGF in bone marrow supernatant are consistant with those mRNA expression level in bone marrow mononuclear cells, moreover, the MVD possitively cerrelates with levels of VEGF, HGF and bFGF in bone marrow supernatant, suggesting that the changes of bone marrow microenvironment vassel-related factors play an important role in angiogenesis and pathogenesis of multiple myeloma.

[Effect of HIF-1α on Angiogenesis-Related Factors in K562 Cells].

OBJECTIVE: To explore the mechanisms of angiogenesis in chronic myeloid leukemia (CML) through detecting the levels of angiogenesis-related factors secreted from K562 cells after overexpression and interference of HIF-1α gene in K562 cells. METHODS: The K562 cells were transfected by lentiviruses carried and interfered HIF-1α gene, then the transtected K562 cells with carried and interfered with HIF-1α gene were enrolled in overexpression and interference groups respectively, at the same time the K562 cells transfected by the empty virus were enrolled in control group. The cells were harvested after culture for 72 hours under normoxid condition. The transfection efficient in 3 groups was detected by fluorescence microscopy; the mRNA expression of HIF-1α gene and angiogenesis-related factors was detected by RT-PCR; the concentration of angiogenesis-related factors in the caltured supernatant was detected by ELISA. RESULTS: The optimal MOI of K562 cells transfected with lentivirus was 10 and the transfection efficiency was about 50%. The positive rate of transfection after screening by puromycin was more than 90%. The mRNA expression of ANG-I, ANG-II, TGF-α and VEGF in the interference group was lower than that in the over-expression group, and the TGF-ß1 mRNA expression in the interference group was higher than in the over-expression group. The mRNA expression of ANG-I and VEGF in the interference group was lower than that in the control group. TGF-αdid not could be detected, and the culture supernatant concentration of ANG-I and TNF-α in the interference group was lower than in the over-expression group, while the VEGF concentration in the interference group was higher than that in the over-expression group. All of the above-mentioned differences were statistically significant (P＜0.05). CONCLUSION: The positive K562 cells transfected with leutivirus have been harvested by screening with puromycin. The HIF-1α mRNA positively regulates the mRNA expression of ANG-1, ANG-2, TGF-α, VEGF in K562 cells, promotes the antocrine ability of ANG-1 and TNF-α, moreover not stimulates the autocrine of TGF-α, the up-regulation of HIF-1α expression can inhibit the expression TGF-ß1 in K562 cells and the autocrine of TGF-ß1.

Optimal dioptric value of near addition lenses intended to slow myopic progression.

PURPOSE: The purpose of this study was to determine the optimal power value of near addition lenses, which would create the least error in accommodative and vergence responses. METHODS: We evaluated accommodative response, phoria, and fixation disparity when the subject viewed through various addition lenses at three working distances for 30 young adults (11 emmetropic, 17 myopic, and 2 hyperopic). Accommodative response was determined with a Canon R-1 infrared optometer under binocular viewing conditions, phoria was determined by the alternating cover test with prism neutralization, and fixation disparity was measured with a Sheedy disparometer. RESULTS: We found that the optimal powers of near addition lenses for the young adult subjects associated with zero retinal defocus were +0.92 D, +1.04 D, and +1.28 D at three viewing distances, 50 cm, 40 cm, and 30 cm, respectively. The optimal powers associated with -3 prism diopters (Delta) near phoria were +0.58 D, +0.35 D, and +0.20 D at the three distances, 50 cm, 40 cm, and 30 cm, respectively. In addition, we found high correlations between the initial accommodative error and the optimal power of the near addition lenses and between the initial near phoria and the optimal power of the near addition lenses. CONCLUSIONS: The results suggest that when the effects of near addition lenses on the accommodative and vergence systems are both considered, the optimal dioptric power of the near addition lens is in a range between +0.20 D and +1.28 D for the three viewing distances. Using progressive lenses to delay the progression of myopia may have promising results if each subject's prescription is customized based on establishing a balance between the accommodative and vergence systems. Formulas derived from this study provide a basis for such considerations.

Changes in accommodative and vergence responses when viewing through near addition lenses.

It has been proposed that near addition lenses reduce the lag of accommodation and may slow myopia progression. In this study, we investigated the differences in accommodative response and near phoria in subjects with and without near addition lenses. Fourteen subjects (7 emmetropic and 7 myopic) participated in the study. Monocular and binocular accommodative responses to a target at 40 cm were measured with and without +2.00 diopter (D) lenses using a Canon R-1 optometer (Canon Europe N.V., Amsterdam, The Netherlands). Near dissociated phoria was measured using 3 testing methods: Maddox rod, cover test, and Von Graefe technique. The differences in accommodative response and near phoria between the 2 viewing conditions with and without the near addition lens were significant (P < 0.0001). No significant differences were revealed in accommodative response and near phoria between refractive error groups. The average accommodative responses of all subjects were 2.03 +/- 0.06 (SE) D (without +2.00 D lens) and 0.60 +/- 0.07 D (with +2.00 D lens) under monocular viewing conditions and 2.11 +/- 0.06 D (without +2.00 D lenses) and 0.77 +/- 0.07 D (with +2.00 D lenses) under binocular viewing conditions. The average near phorias of all subjects were -2.08 +/- 0.69 prism diopters (PD, without +2.00 D lenses) and -7.90 +/- 0.68 PD (with +2.00 D lenses). We discuss the effect of near addition lenses on the defocus of the retinal image, accommodative response, and near phoria. In addition, we propose that the observed difference between binocular and monocular accommodation, when viewing through near addition lenses, is caused by an increased vergence accommodation.

The adaptive effect of narrowing the interocular separation on the AC/A ratio.

The purpose of this study was to determine whether the response AC/A ratio could be altered when the subject's interpupillary distance (IPD) was optically halved. We measured the changes in the AC/A ratio for 10 subjects after using the optical device for 30 min. Accommodative response was measured using a Canon R-1 optometer, and vergence response was measured with an ASL 210 Eye Movement Monitor. The average AC/A ratios were 1.20+/-0.35 (SD) (MA/D) and 0.84+/-0.39 (MA/D) before and after wearing the device, respectively. The decrease in AC/A ratio was statistically significant (p=0.01). This was mainly caused by a reduction in the slope of the accommodative vergence. The results of this study suggest that the AC/A ratio can be decreased if an IPD-narrowing device is used. A possible application of this mechanism in the study of myopia is discussed.

Myopic progression and dark focus variation in optometric students during the first academic year.

BACKGROUND: The aim of this research was to investigate the change in refractive error (RE) of optometric students during their first academic year and whether these changes relate to changes in their dark focus (DF). METHOD: The RE and DF of 64 students were measured objectively every three months during the first academic year, a total of four times, using a Canon R-1 infrared optometer. Thirty-five of the 64 students had an additional RE and DF measurement three weeks immediately after their Summer vacation. Students completed a survey regarding the near work demands they experienced during the Winter break and the teaching semesters. RESULTS: Over nine months, the average RE of the students changed significantly from -2.22 +/- 1.93 (SD) D to -2.50 +/- 2.05 D (p = 0.0002). The rate of myopic progression averaged -0.37 dioptres per year. Inclusion of measurements taken on 35 students immediately after the Summer vacation showed that their change in RE during the Summer vacation was not significant (p = 0.79). For these subjects, the DF measured immediately after the vacation was significantly lower than the DF measured before the vacation (p = 0.007). The reduction in the DF after the vacations corresponded to a period of relative myopic stability in these subjects. CONCLUSIONS: The results of this study suggest that optometric students performing extensive near work are at risk of developing myopia. The variation of their DF values indicates the changing demand for near work during different periods of the year. After Winter and Summer vacations, the DF was lower and the myopic progression was suspended. These findings further support the notion that myopic progression is related to high near work demands and suggest that this progression can be slowed by a period of reduced near work, for example, vacation periods.

[Modification of titanium surface with calcium and phosphorus ions using micro-arc oxidation and its effect on osteoblast attachment].

OBJECTIVE: To study the method for modifying titanium surface with calcium and phosphorus ions using micro-arc oxidation technique and observe osteoblast attachment to the modified surface. METHODS: TA(2) titanium discs were treated with micro-arc oxidation in electrolyte solution containing Ca(2+) and PO(4)(3-). The influence of Ca(2+) and PO(4)(3-) concentrations in the solution and the electrical parameters of the micro-arc oxidization on the content of calcium and phosphorus ions incorporated into the surface of titanium was investigated using energy-dispersive X-ray (EDX) analysis. Scanning electron microscopy was employed for morphological observation of the ceramic coating on the metal surface. The binding strength of ceramic coating with titanium was tested by shear bonding experiment. MC-3T3-E1 osteoblast-like cells were then cultured on the treated surface of titanium discs to investigate the influence of the ceramic coating on osteoblast attachment. RESULTS: Micro-arc oxidation treatment produced a layer of porous TiO(2) coating on the surface of titanium discs, with the average pore size of 2 to 10 mum. EDX analysis revealed that the ceramic coating contained Ca and P elements, whose content had close correlation with Ca(2+) and PO(4)(3-) concentrations in the electrolyte solution and voltage, duty cycle and frequency of micro-arc oxidation. The average bonding strength of the ceramic coating with titanium was 22+/-3 MPa, and TiO(2) coating promoted attachment and spread of osteoblast-like cells on the metal surface as demonstrated by cell culture. CONCLUSIONS: Porous TiO(2) coating can be constructed on the surface of titanium using micro-arc oxidation, and Ca(2+) and PO(4)(3-) incorporated into the coating can improve the biocompatibility of titanium. The content of Ca(2+) and PO(4)(3-) in the coating can be modulated by adjusting the concentrations of Ca(2+) and PO(4)(3-) in the electrolyte solution and the electrical parameters of the micro-arc oxidation.