[The value of DCE-MRI in predicting IDH gene mutation of high-grade gliomas]

Objective: To investigate the value of quantitative and semiquantitative parameters of DCE-MRI in predicting IDH gene mutation of high-grade gliomas before the operation. Methods: Twenty-six individuals with surgically and pathologically proved WHO Ⅲ-Ⅳ gliomas collected from April 2016 to June 2019 in First People's Hospital of Changzhou, were divided into two groups, IDH mutation group (7 cases, 27-67 years, 3 males and 4 females,) and IDH gene wild group (19 cases, 42-75 years, 12 males and 7 females) according to the results of molecular pathology. All individuals underwent conventional plain (T(1)WI, T(2)WI), enhanced MR scanning (T(1)WI) and dynamic contrast enhancement (DCE). Four quantitative parameters:volume transfer constant (K(trans)), ratio constant of back flux (Kep), extravascular extracellular space fractional volume (Ve), and blood plasma fractional volume (Vp), and four semiquantitative parameters: time to peak (TTP), maximum concentration (MAX Conc), initial area under the gadolinium concentration-time curve (IAUC) and maximum slope of decrease (MAX Slope) were measured. The independent samples t test (normal distribution and homogeneity of variance) or Mann-Whitney rank sum test (abnormal distribution or heterogeneity of variance) were used to compare the differences of quantitative and semiquantitative parameters between IDH gene mutation group and IDH gene wild type group. Receiver operating characteristic (ROC) curve was used to evaluate the efficiency of quantitative and semiquantitative parameters in predicting IDH gene mutation of high-grade gliomas. Results: The value of K(trans),TTP in IDH mutated group were 0.096 (0.080,0.135)/min and (3.95±0.34) s, respectively. The value of K(trans), TTP in IDH wild type group were 0.168 (0.132, 0.337)/min and (2.58±1.15) s, respectively. The value of K(trans) in IDH mutated group was significantly less than the value of K(trans) in IDH gene wild type group (Z value was -2.168, P value was 0.030). The value of K(trans) in IDH mutated group was significantly greater than the value of K(trans) in IDH gene wild type (Z value was -2.630, P value was 0.007). The area under the ROC curve (AUC) of K(trans) and TTP in predicting IDH gene mutation of high-grade gliomas was 0.782 and 0.842, respectively. The specificity of K(trans) was higher (73.7%), The sensitivity of TTP was the higher (100.0%). Combined K(trans)and TTP were the best for predicting IDH gene mutation of high-grade gliomas, AUC was 0.865. Conclusion: Quantitative and semiquantitative parameters of DCE-MRI can help to predict IDH gene mutation of high-grade gliomas before the operation.

Targeted genetic screening in mice through haploid embryonic stem cells identifies critical genes in bone development

Mutagenic screening is powerful for identifying key genes involved in developmental processes. However, such screens are successful only in lower organisms. Here, we develop a targeted genetic screening approach in mice through combining androgenetic haploid embryonic stem cells (AG-haESCs) and clustered regularly interspaced palindromic repeats/CRISPR-associated protein 9 (CRISPR-Cas9) technology. We produced a mutant semi-cloned (SC) mice pool by oocyte injection of AG-haESCs carrying constitutively expressed Cas9 and an single guide RNA (sgRNA) library targeting 72 preselected genes in one step and screened for bone-development-related genes through skeletal analysis at birth. This yielded 4 genes: Zic1 and Clec11a, which are required for bone development, and Rln1 and Irx5, which had not been previously considered. Whereas Rln1-/- mice exhibited small skeletal size only at birth, Irx5-/- mice showed skeletal abnormalities both in postnatal and adult phases due to decreased bone mass and increased bone marrow adipogenesis. Mechanistically, iroquois homeobox 5 (IRX5) promotes osteoblastogenesis and inhibits adipogenesis by suppressing peroxisome proliferator activated receptor γ (PPARγ) activation. Thus, AG-haESC-mediated functional mutagenic screening opens new avenues for genetic interrogation of developmental processes in mice.

[Feasibility study of blood oxygen level-dependent magnetic resonance imaging in evaluating the response of metastatic lymph nodes of rabbit VX2 tumor to radiotherapy]

Objective: To investigate the feasibility of blood oxygen level-dependent magnetic resonance imaging (MRI) in evaluating the response of metastatic lymph nodes of rabbit VX2 tumor to radiotherapy. Methods: Twenty-eight healthy New Zealand white rabbits which were provided by the Laboratory Animal Center of Soochow University, male or female, 2 to 3 months, weighing 2 to 3 kg, were used to establish the animal model of VX2 tumor popliteal fossa metastatic lymph node, and then were divided into either the radiotherapy group (n=16) or the control group (n=12). The radiotherapy group received a 20 Gy radiotherapy per rabbit, the control group received sham radiotherapy. All rabbits underwent MRI scan on four time points, including before (0 day), 3rd, 7th and 14th days after radiotherapy. The two parameters of size and R(2*) value (s(-1)) of lymph node were measured. At each time point,two rabbits in each group were sacrificed randomly to resect lymph nodes for pathological examination, and two parameters of microvessel density (MVD, strip/HP) and apoptosis index (AI, %) were analyzed. The parameters among the four time points in each group or between the two groups were compared. The correlation of lymph node size and R(2*) value with MVD or AI was analyzed, respectively. Results: A significant size difference was neither between the two groups or among the each time points in each group (P>0.05). The R(2*) of lymph node in the radiotherapy group was (29.6±1.7),(36.8±2.6),(44.8±5.8) and (57.7±6.2) s(-1) at the time points of 0, 3, 7 and 14 days, respectively, showing a gradual increase trend; MVD was (52.3±2.5),(41.0±3.6),(34.0±3.6) and (22.7±2.5) strip/HP respectively, showing a decreasing trend; AI was 12.8%±0.5%,14.9%±0.6%,20.6%±0.5% and 27.5%±0.7% respectively, showing a gradual increase trend (all P<0.05). In the control group, both R(2*) value and AI among the four time points did not change statistically (all P>0.05), but MVD showed a gradual increase trend,(50.0±3.0),(53.0±1.7),(60.3±2.5) and (70.0±2.0) strip/HP, respectively, P<0.05. There were significant differences in R(2*) and MVD at 3, 7 and 14 days, in AI at 7 and 14 days between the two groups (all P<0.05). There was a linear correlation of R(2*) value, but not of size, with MVD and AI (r=-0.87 and 0.94, respectively). Conclusion: Blood oxygen level-dependent MRI can indirectly reflect the hypoxic status of metastatic lymph nodes after radiotherapy, and has potential value in evaluating the response of metastatic lymph nodes to radiotherapy.

A Surrogate Modelling Approach Based on Nonlinear Dimension Reduction for Uncertainty Quantification in Groundwater Flow Models

In this paper, we develop a surrogate modelling approach for capturing the output field (e.g. the pressure head) from groundwater flow models involving a stochastic input field (e.g. the hydraulic conductivity). We use a Karhunen–Loève expansion for a log-normally distributed input field and apply manifold learning (local tangent space alignment) to perform Gaussian process Bayesian inference using Hamiltonian Monte Carlo in an abstract feature space, yielding outputs for arbitrary unseen inputs. We also develop a framework for forward uncertainty quantification in such problems, including analytical approximations of the mean of the marginalized distribution (with respect to the inputs). To sample from the distribution, we present Monte Carlo approach. Two examples are presented to demonstrate the accuracy of our approach: a Darcy flow model with contaminant transport in 2-d and a Richards equation model in 3-d.

Lkb1 deletion in periosteal mesenchymal progenitors induces osteogenic tumors through mTORC1 activation

Bone osteogenic sarcoma has a poor prognosis as the exact cell of origin and the signaling pathways underling tumor formation remain undefined. Here, we report an osteogenic tumor mouse model based on the conditional knockout of liver kinase b1 (Lkb1; also known as Stk11) in Cathepsin K (Ctsk)-Cre expressing cells. Lineage tracing studies demonstrated that Ctsk-Cre could label a population of periosteal cells. The cells functioned as mesenchymal progenitors with regard to markers and functional properties. LKB1 deficiency increased proliferation and osteoblast differentiation of Ctsk+ periosteal cells, while downregulation of mTORC1 activity, using Raptor genetic mouse model or mTORC1 inhibitor treatment, ameliorated tumor progression of Ctsk-Cre Lkb1fllfl mice. Xenograft mouse models, using human osteosarcoma cell lines, also demonstrated that LKB1 deficiency promoted tumor formation, while mTOR inhibition suppressed xenograft tumor growth. In summary, we identified periosteum-derived Ctsk-Cre expressing cells as a cell of origin for osteogenic tumor and suggested the LKB1-mTORC1 pathway as a promising target for treatment of osteogenic tumor.

Abstract P6-02-14: Not presented

This abstract was not presented at the conference.

Citation Format: Xing W, Li Q, Sun G, Cao R, Chen B, Jiang C, Ma L, Wang K. Not presented [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P6-02-14.

Abstract P2-01-25: Not presented

This abstract was not presented at the conference.

Citation Format: Xing W, Li Q, Cao R, Sun G, Jiang C. Not presented [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr P2-01-25.

Phenotype Extraction Based on Word Embedding to Sentence Embedding Cascaded Approach

As a significant determinant in the development of named entity recognition, phenotypic descriptions are normally presented differently in biomedical literature with the use of complicated semantics. In this paper, a novel approach has been proposed to identify plant phenotypes by adopting word embedding to sentence embedding cascaded approach. We make use of a word embedding method to find high-frequency phenotypes with original sentences used as input in a sentence embedding method. In doing so, a variety of complicated phenotypic expressions can be recognized accurately. Besides, the state-of-the-art word representation models have been compared and among them, skip-gram with negative sampling was selected with the best performance. To evaluate the performance of our approach, we applied it to the dataset composed of 56 748 PubMed abstracts of model organism Arabidopsis thaliana. The experiment results showed that our approach yielded the best performance, as it achieved a 2.588-fold increase in terms of the number of new phenotypic descriptions when compared to the original phenotype ontology.

A gene-phenotype relationship extraction pipeline from the biomedical literature using a representation learning approach

Motivation

The fundamental challenge of modern genetic analysis is to establish gene-phenotype correlations that are often found in the large-scale publications. Because lexical features of gene are relatively regular in text, the main challenge of these relation extraction is phenotype recognition. Due to phenotypic descriptions are often study- or author-specific, few lexicon can be used to effectively identify the entire phenotypic expressions in text, especially for plants.

Results

We have proposed a pipeline for extracting phenotype, gene and their relations from biomedical literature. Combined with abbreviation revision and sentence template extraction, we improved the unsupervised word-embedding-to-sentence-embedding cascaded approach as representation learning to recognize the various broad phenotypic information in literature. In addition, the dictionary- and rule-based method was applied for gene recognition. Finally, we integrated one of famous information extraction system OLLIE to identify gene-phenotype relations. To demonstrate the applicability of the pipeline, we established two types of comparison experiment using model organism Arabidopsis thaliana. In the comparison of state-of-the-art baselines, our approach obtained the best performance (F1-Measure of 66.83%). We also applied the pipeline to 481 full-articles from TAIR gene-phenotype manual relationship dataset to prove the validity. The results showed that our proposed pipeline can cover 70.94% of the original dataset and add 373 new relations to expand it.

Availability and implementation

The source code is available at http://www.wutbiolab.cn: 82/Gene-Phenotype-Relation-Extraction-Pipeline.zip.

Supplementary information

Supplementary data are available at Bioinformatics online.

Periodontal infection with Porphyromonas gingivalis induces preterm birth and lower birth weight in rats

Preterm birth (PTB), accompanied by low birth weight (LBW) or not, is a syndrome with tremendous risk factors and long-term health consequences for children. In recent decades, overwhelming studies have shown that periodontitis contributes to prematurity and LBW. This study was conducted to determine the link between maternal periodontitis and the pathogenesis of PTB and/or LBW through a rat infection model induced by Porphyromonas gingivalis, an important periodontopathic bacterium. The murine model was established by surgically ligating the left mandibular first molars and inoculating with P. gingivalis, and then all female rats initiated mating 6 weeks post infection. The gestational day and birth weight were recorded, and blood, amniotic fluid, and placental specimens were collected. Rats with a PTB and LBW newborns were observed in the P. gingivalis-infected group. Additionally, P. gingivalis infection significantly increased the maternal serum levels of interferon-γ and interleukin-1β, whereas no significant difference in the cytokine response was observed in the amniotic fluid. Moreover, with the translocation of P. gingivalis to placentas, remarkable changes in gestational tissues were found, followed by significantly enhanced expression of Toll-like receptor 2 (TLR2) as well as Fas and Fas ligand (FasL). These results support the concept that severe cases of periodontitis caused by P. gingivalis infection may be indicative of rats being more susceptible to PTB/LBW, probably through the activation of the TLR2 and Fas/FasL pathways within the placental tissues. This study gave us new insight into how maternal periodontopathogens might be linked to placental damage and premature pathogenesis.

Paracrine and endocrine actions of bone-the functions of secretory proteins from osteoblasts, osteocytes, and osteoclasts

The skeleton is a dynamic organ that is constantly remodeled. Proteins secreted from bone cells, namely osteoblasts, osteocytes, and osteoclasts exert regulation on osteoblastogenesis, osteclastogenesis, and angiogenesis in a paracrine manner. Osteoblasts secrete a range of different molecules including RANKL/OPG, M-CSF, SEMA3A, WNT5A, and WNT16 that regulate osteoclastogenesis. Osteoblasts also produce VEGFA that stimulates osteoblastogenesis and angiogenesis. Osteocytes produce sclerostin (SOST) that inhibits osteoblast differentiation and promotes osteoclast differentiation. Osteoclasts secrete factors including BMP6, CTHRC1, EFNB2, S1P, WNT10B, SEMA4D, and CT-1 that act on osteoblasts and osteocytes, and thereby influenceaA osteogenesis. Osteoclast precursors produce the angiogenic factor PDGF-BB to promote the formation of Type H vessels, which then stimulate osteoblastogenesis. Besides, the evidences over the past decades show that at least three hormones or "osteokines" from bone cells have endocrine functions. FGF23 is produced by osteoblasts and osteocytes and can regulate phosphate metabolism. Osteocalcin (OCN) secreted by osteoblasts regulates systemic glucose and energy metabolism, reproduction, and cognition. Lipocalin-2 (LCN2) is secreted by osteoblasts and can influence energy metabolism by suppressing appetite in the brain. We review the recent progresses in the paracrine and endocrine functions of the secretory proteins of osteoblasts, osteocytes, and osteoclasts, revealing connections of the skeleton with other tissues and providing added insights into the pathogenesis of degenerative diseases affecting multiple organs and the drug discovery process.

[Feasibility of using blood oxygen level-dependent magnetic resonance imaging to evaluate renal fibrosis of ureteral obstruction of rabbits]

Objective: To investigate the feasibility of blood oxygen level-dependent MR (MR-BOLD) in assessing renal fibrosis of ureteral obstruction of rabbits. Methods: Forty healthy New Zealand rabbits were randomly divided into control group (n=8) and unilateral ureteral obstruction (UUO) group (n=32). The rabbits in the UUO group were subjected to unilateral ureteral obstruction of the left kidney.Coronal T(2) weighted imaging (T(2)WI) and axial BOLD examinations were performed before operation, 2, 4, 6 and 8 W after operation (each subgroup n=8). After the examinations, nephrectomy was performed for histologic evaluation.The T(2)(*) relaxation rate of the renal cortex (CR(2)(*)) , medulla (MR(2)(*)) and the same level of muscle(R(2)(*)(muscle)) were measured separately.The normalization of the cortex and medulla (sR(2)(*)), and the difference of sR(2)(*) between renal cortex and medulla before and post UUO (ΔsR(2)(*)) were calculated.The differences of sR(2)(*) (sCR(2)(*), sMR(2)(*), sCR(2)(*)(control), sMR(2)(*)(control)) at each time point between control and UUO group were compared by using independent sample t test.The LSD test was used to compare the sR(2)(*) in the control with that in the UUO group.The ΔsCR(2)(*) and ΔsMR(2)(*) values of the subgroups at UUO 2, 4, 6 and 8 W were compared by independent sample t test. Results: The sR(2)(*) values of UUO group were all lower than those of control group (all P<0.05), while sR(2)(*)(control) and sR(2)(*) in UUO group before operation were not significant different (P>0.05). The sCR(2)(*) values of UUO 2, 4, 6 and 8 W were 0.32±0.01, 0.37±0.01, 0.47±0.02 and 0.50±0.03.The sMR(2)(*) values were 0.39±0.02, 0.48±0.02, 0.58±0.04 and 0.65±0.05.There were significant differences of sCR(2)(*)(between) UUO 2 W and 6 W, UUO 2 W and 8 W, UUO 4 W and 6 W, UUO 4 W and 8 W (all P<0.01). There were significant differences of sMR(2)(*)(between) UUO 2 W and 6 W, UUO 2 W and 8 W, UUO 4 W and 8 W (all P<0.01). No significant difference was founded between sCR(2)(*) and sMR(2)(*) at each time point in control group (P>0.05). The ΔsCR(2)(*) values of UUO 2, 4, 6 and 8 W were 0.31±0.02, 0.20±0.02, 0.14±0.20, 0.09±0.04; the ΔsMR(2)(*) values were 0.51±0.05, 0.36±0.04, 0.28±0.05, 0.19±0.05. The ΔsCR(2)(*) values of UUO 2, 4 and 6 W were less than ΔsMR(2)(*) (P<0.05). There was no significant difference between ΔsCR(2)(*) and ΔsMR(2)(*) in UUO 8 W (P>0.05). Conclusions: The R(2)(*) change in medulla resulted from renal fibrosis is more significant than cortex.MR-BOLD can reflect the process of renal fibrosis.It's feasible and of great value to use renal MR-BOLD for the assessment of renal fibrosis induced by unilateral ureteral obstruction.

[Feasibility of susceptibility weighted imaging in the evaluation of renal fibrosis induced by unilateral ureteral obstruction in white rabbits]

Objective: To evaluate the feasibility of susceptibility weighted imaging (SWI) to reflect the progression of renal fibrosis (RF) induced by unilateral ureteral obstruction (UUO) in rabbits. Methods: Total of 32 New Zealand white rabbits (aged 4-5 months) were used to establish the UUO model (RF group) and were divided into 4 subgroups (including RF-2W group, RF-4W group , RF-6W group and RF-8W group) according to the duration of obstruction(2, 4, 6 and 8 weeks). Another 8 rabbits, underwent sham operation, were used as the control group (Sham group). T2 weighted imaging (T(2)WI) and SWI sequences were performed on the rabbits at 2, 4, 6 and 8 weeks after UUO, respectively. Rabbits were sacrificed and specimens were obtained after MR examination, then hematoxylin&eosin staining and Masson staining were performed to explore the renal tubular injury and interstitial fibrosis. The relationship between the relative signal-to-noise ratio(rSNR)measured on the SWI fusion maps and the degree of renal fibrosis was analyzed. Analysis of variance was used to compare the difference of rSNRs among 5 groups. Results: The rSNRs of the inner medulla gradually decreased over time, the rSNRs of Sham group, RF-2W group, RF-4W group, RF-6W group and RF-8W group was 2.29±0.18, 1.73±0.30, 1.67±0.08, 1.42±0.28, 1.12±0.15, respectively (F=25.876, P<0.01). In RF-2W group, the rSNRs of the outer medulla and cortex increased when compared with those in the Sham group, then they decreased gradually over time in the other groups (F=5.230, 7.621, both P<0.05). The pathological results demonstrated that the morphology of renal tubules in Sham group was normal and no apoptosis or necrosis was seen in the tubular epithelial cells. In the RF-2W group, the main pathological manifestations included renal tubules dilatation, tubular epithelial cell degeneration or necrosis and the infiltration of inflammatory cell. In Masson staining, lots of blue areas were present in the interstitial. As the obstruction time prolonged, the renal tubular collapse, inflammatory cells and fibroblasts increased significantly. Meanwhile blue area also increased significantly in Masson staining. Conclusion: SWI can be used as a noninvasive method to evaluate the pathological progression of fibrosis in the rabbit UUO model.

Neoadjuvant versus adjuvant chemoradiation for stage II-III esophageal squamous cell carcinoma: a single institution experience

Esophageal cancer is the eighth most common cancer worldwide. It is the fourth most common cause of cancer death in China and esophageal squamous cell carcinoma (ESCC) is the most prevalent histologic type. Many clinical trials have explored the value of neoadjuvant or adjuvant chemoradiation therapy in potentially resectable ESCC; however, these studies have produced conflicting results. This retrospective study was performed to investigate whether patients with resectable stage II/III ESCC should receive neoadjuvant or adjuvant therapy in addition to surgery. A review of stage II/III thoracic ESCC patients who underwent esophagectomy and either neoadjuvant or adjuvant chemoradiation was performed. Chemotherapy regimen consisted of cisplatin 75 mg/m2 divided into 3 days and fluorouracil 500 mg/m2 on days 1 to 5. The patients who underwent neoadjuvant therapy were treated with one cycle of chemotherapy concurrently with radiotherapy (40 Gy in 20 fractions, 5 days/week), and those receiving adjuvant therapy were treated with two cycles of chemotherapy concurrently with radiotherapy (46-50 Gy in 23-25 fractions, 5 days/week). A total of 122 patients met inclusion criteria, of which 49 underwent neoadjuvant chemoradiation and 73 underwent adjuvant chemoradiation. Median follow up was 36.5 months. The median survival times and 3, 5-year overall survival (OS) rates for the neoadjuvant and adjuvant groups were 39.3 versus 31.5 months, and 53.0%, 45.7% versus 42.9%, 29.7%, respectively (P = 0.091). For the patients with stage III ESCC, the median survival times and 5-year OS rates for the neoadjuvant and adjuvant groups were 39.3 versus 21.3 months, and 43.4% versus 21.0%, respectively (P = 0.021). Among lymph node-positive patients, the median survival times and 5-year OS rates for the neoadjuvant and adjuvant groups were 55.6 versus 23.7 months, and 43.0% versus 25.7%, respectively (P = 0.085). The incidence of perioperative and postoperative complications was comparable between the two groups (P > 0.05). For patients with resectable stage II/III ESCC, neoadjuvant chemoradiation does not increase postoperative complications and is associated with a trend toward better OS when compared to adjuvant chemoradiation.

[Correlation between susceptibility weighted imaging manifestation and serum cystatin C for delayed graft function]

OBJECTIVE

To explore the correlation between susceptibility weighted imaging (SWI) manifestation and serum cystatin C level for delayed graft function (DGF).

METHODS

The conventional MRI, SWI and serum cystatin C of 27 cases with DGF in nephrotransplantation center in Third Affiliated Hospital of Suzhou University from September 2014 and August 2015 were retrospectively analyzed.By contrasting conventional MRI images of transplanted kidney in DGF, the imaging manifestations of benign tumors such as cysts and angiomyolipomas were excluded on SWI images, and then making the renal cortex as the reference, if the abnormal signal lesions were found in the transplanted kidney, the location and signal intensity would be analyzed. The differences in serum cystatin C level between DGF groups without and with abnormal signal lesions were compared by using independent-sample t-test.The correlation between SWI manifestation and serum cystatin C level for DGF was assessed with Spearman rank correlation analysis.

RESULTS

A total of 15 cases were found without abnormal signal lesions and the average value of their serum cystatin C level was (2.92±0.44) mg/L.A total of 12 cases were found with abnormal low signal lesions located at junctional zone between cortex and medulla, and the average value of their serum cystatin C level was (6.91±0.96) mg/L. The differences in serum cystatin C level between the two DGF groups were statistically significant (t=-4.040, P=0.000). There was a positive correlation between the abnormal low signal lesions on SWI and serum cystatin C level (r=0.660, P=0.000).

CONCLUSION

The status of renal function impairment could be reflected by being with or without abnormal signal lesions on SWI. A relatively big renal function impairment may be predicted by the appearance of abnormal low signal lesions at junctional zone between cortex and medulla on SWI.

Anisotropic Li intercalation in a Li(x)FePO4 nano-particle: a spectral smoothed boundary phase-field model

A spectral smoothed boundary phase-field model is implemented to study lithium (Li) intercalation in a LixFePO4 nano-particle immersed in a Li(+) rich electrolyte. It takes into account different physical processes on the particle surface, such as heterogeneous nucleation, Li flux and stress-free boundary conditions. We show the nucleation and growth of plate-like Li-rich crystallites along the (010) plane due to the high Li mobility along [001]. Since such plate-like crystallites, which are nucleated from (001) surfaces, align their phase boundaries along the (101) habit planes, a LixFePO4 nano-particle with prominent (010) and (001) surface facets and the longest axis length along [100] is proposed to exhibit great mechanical stability.