Molecular risk assessment of BIG 1-98 participants by expression profiling using RNA from archival tissue.

Background: The purpose of the work reported here is to test reliable molecular profiles using routinely processed formalin-fixed paraffin-embedded (FFPE) tissues from participants of the clinical trial BIG 1-98 with a median follow-up of 60 months. Methods: RNA from fresh frozen (FF) and FFPE tumor samples of 82 patients were used for quality control, and independent FFPE tissues of 342 postmenopausal participants of BIG 1-98 with ER-positive cancer were analyzed by measuring prospectively selected genes and computing scores representing the functions of the estrogen receptor (eight genes, ER_8), the progesterone receptor (five genes, PGR_5), Her2 (two genes, HER2_2), and proliferation (ten genes, PRO_10) by quantitative reverse transcription PCR (qRT-PCR) on TaqMan Low Density Arrays. Molecular scores were computed for each category and ER_8, PGR_5, HER2_2, and PRO_10 scores were combined into a RISK_25 score. Results: Pearson correlation coefficients between FF- and FFPE-derived scores were at least 0.94 and high concordance was observed between molecular scores and immunohistochemical data. The HER2_2, PGR_ 5, PRO_10 and RISK_25 scores were significant predictors of disease free-survival (DFS) in univariate Cox proportional hazard regression. PRO_10 and RISK_25 scores predicted DFS in patients with histological grade II breast cancer and in lymph node positive disease. The PRO_10 and PGR_ 5 scores were independent predictors of DFS in multivariate Cox regression models incorporating clinical risk indicators; PRO_10 outperformed Ki-67 labeling index in multivariate Cox proportional hazard analyses. Conclusions: Scores representing the endocrine responsiveness and proliferation status of breast cancers were developed from gene expression analyses based on RNA derived from FFPE tissues. The validation of the molecular scores with tumor samples of participants of the BIG 1-98 trial demonstrates that such scores can serve as independent prognostic factors to estimate disease free survival (DFS) in postmenopausal patients with estrogen receptor positive breast cancer.

Alteration in neuromuscular function after a 5 km running time trial

The aim of this study was to characterize the effect of a 5 km running time trial on the neuromuscular properties of the plantar flexors. Eleven well-trained triathletes performed a series of neuromuscular tests before and immediately after the run on a 200 m indoor track. Muscle activation (twitch interpolation) and normalized EMG activity were assessed during maximal voluntary contraction (MVC) of plantar flexors. Maximal soleus H-reflexes and M-waves were evoked at rest (i.e. H (MAX) and M (MAX), respectively) and during MVC (i.e. H (SUP) and M (SUP), respectively). MVC significantly declined (-27%; P < 0.001) after the run, due to decrease in muscle activation (-8%; P < 0.05) and M (MAX)-normalized EMG activity (-13%; P < 0.05). Significant reductions in M-wave amplitudes (M (MAX): -13% and M (SUP): -16%; P < 0.05) as well as H (MAX)/M (MAX) (-37%; P < 0.01) and H (SUP)/M (SUP) (-25%; P < 0.05) ratios occurred with fatigue. Following exercise, the single twitch was characterized by lower peak torque (-16%; P < 0.001) as well as shorter contraction (-19%; P < 0.001) and half-relaxation (-24%; P < 0.001) times. In conclusion, the reduction in plantar flexors strength induced by a 5 km running time trial is caused by peripheral adjustments, which are attributable to a failure of the neuromuscular transmission and excitation-contraction coupling. Fatigue also decreased the magnitude of efferent motor outflow from spinal motor neurons to the plantar flexors and part of this suboptimal neural drive is the result of an inhibition of soleus motoneuron pool reflex excitability.

Rapid profiling of intact glucosinolates in Arabidopsis leaves by UHPLC-QTOFMS using a charged surface hybrid column.

INTRODUCTION: The analysis of glucosinolates (GS) is traditionally performed by reverse-phase liquid chromatography coupled to ultraviolet detection after a time-consuming desulphation step, which is required for increased retention. Simpler and more efficient alternative methods that can shorten both sample preparation and analysis are much needed. OBJECTIVE: To evaluate the feasibility of using ultrahigh-pressure liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (UHPLC-QTOFMS) for the rapid profiling of intact GS. METHODOLOGY: A simple and short extraction of GS from Arabidopsis thaliana leaves was developed. Four sub-2 µm reverse-phase columns were tested for the rapid separation of these polar compounds using formic acid as the chromatographic additive. High-resolution QTOFMS was used to detect and identify GS. RESULTS: A novel charged surface hybrid (CSH) column was found to provide excellent retention and separation of GS within a total running time of 11 min. Twenty-one GS could be identified based on their accurate mass as well as isotopic and fragmentation patterns. The method was applied to determine the changes in GS content that occur after herbivory in Arabidopsis. In addition, we evaluated its applicability to the profiling of other Brassicaceae species. CONCLUSION: The method developed can profile the full range of GS, including the most polar ones, in a shorter time than previous methods, and is highly compatible with mass spectrometric detection.

Profiling of 19-norandrosterone sulfate and glucuronide in human urine: implications in athlete's drug testing.

19-Norandrosterone (19-NA) as its glucuronide derivative is the target metabolite in anti-doping testing to reveal an abuse of nandrolone or nandrolone prohormone. To provide further evidence of a doping with these steroids, the sulfoconjugate form of 19-norandrosterone in human urine might be monitored as well. In the present study, the profiling of sulfate and glucuronide derivatives of 19-norandrosterone together with 19-noretiocholanolone (19-NE) were assessed in the spot urines of 8 male subjects, collected after administration of 19-nor-4-androstenedione (100mg). An LC/MS/MS assay was employed for the direct quantification of sulfoconjugates, whereas a standard GC/MS method was applied for the assessment of glucuroconjugates in urine specimens. Although the 19-NA glucuronide derivative was always the most prominent at the excretion peak, inter-individual variability of the excretion patterns was observed for both conjugate forms of 19-NA and 19-NE. The ratio between the glucuro- and sulfoconjugate derivatives of 19-NA and 19-NE could not discriminate the endogenous versus the exogenous origin of the parent compound. However, after ingestion of 100mg 19-nor-4-androstenedione, it was observed in the urine specimens that the sulfate conjugates of 19-NA was detectable over a longer period of time with respect to the other metabolites. These findings indicate that more interest shall be given to this type of conjugation to deter a potential doping with norsteroids.

False identity documents forensic profiling: study of known-sources and cross-border seizures

A forensic intelligence process was conducted over cross-border seizures of false identity documents whose sources were partly known to be the same. Visual features of 300 counterfeit Portuguese and French identity cards seized in France and Switzerland were observed and integrated in a structured database developed to detect and analyze forensic links. Based on a few batches of documents known to come from common sources, the forensic profiling method could be validated and its performance evaluated. The method also proved efficient and complementary to conventional means of detecting connections between cases. Cross-border links were detected, highlighting the need for more collaboration. Forensic intelligence could be produced, uncovering the structure of counterfeits' illegal trade, the concentration of their sources and the evolution of their quality over time. In addition, two case examples illustrated how forensic profiling may support specific investigations. The forensic intelligence process and its results will underline the need to develop such approaches to support the fight against fraudulent documents and organized crime.

Subject-based steroid profiling and the determination of novel biomarkers for DHT and DHEA misuse in sports.

Doping with natural steroids can be detected by evaluating the urinary concentrations and ratios of several endogenous steroids. Since these biomarkers of steroid doping are known to present large inter-individual variations, monitoring of individual steroid profiles over time allows switching from population-based towards subject-based reference ranges for improved detection. In an Athlete Biological Passport (ABP), biomarkers data are collated throughout the athlete's sporting career and individual thresholds defined adaptively. For now, this approach has been validated on a limited number of markers of steroid doping, such as the testosterone (T) over epitestosterone (E) ratio to detect T misuse in athletes. Additional markers are required for other endogenous steroids like dihydrotestosterone (DHT) and dehydroepiandrosterone (DHEA). By combining comprehensive steroid profiles composed of 24 steroid concentrations with Bayesian inference techniques for longitudinal profiling, a selection was made for the detection of DHT and DHEA misuse. The biomarkers found were rated according to relative response, parameter stability, discriminative power, and maximal detection time. This analysis revealed DHT/E, DHT/5β-androstane-3α,17β-diol and 5α-androstane-3α,17β-diol/5β-androstane-3α,17β-diol as best biomarkers for DHT administration and DHEA/E, 16α-hydroxydehydroepiandrosterone/E, 7β-hydroxydehydroepiandrosterone/E and 5β-androstane-3α,17β-diol/5α-androstane-3α,17β-diol for DHEA. The selected biomarkers were found suitable for individual referencing. A drastic overall increase in sensitivity was obtained. The use of multiple markers as formalized in an Athlete Steroidal Passport (ASP) can provide firm evidence of doping with endogenous steroids.

Novel markers of the human follicle-associated epithelium identified by genomic profiling and microdissection.

BACKGROUND & AIMS: Regulation of gene expression in the follicle-associated epithelium (FAE) over Peyer's patches is largely unknown. CCL20, a chemokine that recruits immature dendritic cells, is one of the few FAE-specific markers described so far. Lymphotoxin beta (LTalpha1beta2) expressed on the membrane of immune cells triggers CCL20 expression in enterocytes. In this study, we measured expression profiles of LTalpha1beta2-treated intestinal epithelial cells and selected CCL20 -coregulated genes to identify new FAE markers. METHODS: Genomic profiles of T84 and Caco-2 cell lines treated with either LTalpha1beta2, flagellin, or tumor necrosis factor alpha were measured using the Affymetrix GeneChip U133A. Clustering analysis was used to select CCL20 -coregulated genes, and laser dissection microscopy and real-time polymerase chain reaction on human biopsy specimens was used to assess the expression of the selected markers. RESULTS: Applying a 2-way analysis of variance, we identified regulated genes upon the different treatments. A subset of genes involved in inflammation and related to the nuclear factor kappaB pathway was coregulated with CCL20 . Among these genes, the antiapoptotic factor TNFAIP3 was highly expressed in the FAE. CCL23 , which was not coregulated in vitro with CCL20 , was also specifically expressed in the FAE. CONCLUSIONS: We have identified 2 novel human FAE specifically expressed genes. Most of the CCL20 -coregulated genes did not show FAE-specific expression, suggesting that other signaling pathways are critical to modulate FAE-specific gene expression.

An ultra performance liquid chromatography-tandem MS assay for tamoxifen metabolites profiling in plasma: first evidence of 4'-hydroxylated metabolites in breast cancer patients.

There is increasing evidence that the clinical efficacy of tamoxifen, the first and most widely used targeted therapy for estrogen-sensitive breast cancer, depends on the formation of the active metabolites 4-hydroxy-tamoxifen and 4-hydroxy-N-desmethyl-tamoxifen (endoxifen). Large inter-individual variability in endoxifen plasma concentrations has been observed and related both to genetic and environmental (i.e. drug-induced) factors altering CYP450s metabolizing enzymes activity. In this context, we have developed an ultra performance liquid chromatography-tandem mass spectrometry method (UPLC-MS/MS) requiring 100 μL of plasma for the quantification of tamoxifen and three of its major metabolites in breast cancer patients. Plasma is purified by a combination of protein precipitation, evaporation at room temperature under nitrogen, and reconstitution in methanol/20 mM ammonium formate 1:1 (v/v), adjusted to pH 2.9 with formic acid. Reverse-phase chromatographic separation of tamoxifen, N-desmethyl-tamoxifen, 4-hydroxy-tamoxifen and 4-hydroxy-N-desmethyl-tamoxifen is performed within 13 min using elution with a gradient of 10 mM ammonium formate and acetonitrile, both containing 0.1% formic acid. Analytes quantification, using matrix-matched calibration samples spiked with their respective deuterated internal standards, is performed by electrospray ionization-triple quadrupole mass spectrometry using selected reaction monitoring detection in the positive mode. The method was validated according to FDA recommendations, including assessment of relative matrix effects variability, as well as tamoxifen and metabolites short-term stability in plasma and whole blood. The method is precise (inter-day CV%: 2.5-7.8%), accurate (-1.4 to +5.8%) and sensitive (lower limits of quantification comprised between 0.4 and 2.0 ng/mL). Application of this method to patients' samples has made possible the identification of two further metabolites, 4'-hydroxy-tamoxifen and 4'-hydroxy-N-desmethyl-tamoxifen, described for the first time in breast cancer patients. This UPLC-MS/MS assay is currently applied for monitoring plasma levels of tamoxifen and its metabolites in breast cancer patients within the frame of a clinical trial aiming to assess the impact of dose increase on tamoxifen and endoxifen exposure.

Development of a system for 3D high-resolution seismic reflection profiling on lakes

A high-resolution three-dimensional (3D) seismic reflection system for small-scale targets in lacustrine settings has been developed. Its main characteristics include navigation and shot-triggering software that fires the seismic source at regular distance intervals (max. error of 0.25 m) with real-time control on navigation using differential GPS (Global Positioning System). Receiver positions are accurately calculated (error < 0.20 m) with the aid of GPS antennas attached to the end of each of three 24-channel streamers. Two telescopic booms hold the streamers at a distance of 7.5 m from each other. With a receiver spacing of 2.5 m, the bin dimension is 1.25 m in inline and 3.75 m in crossline direction. To test the system, we conducted a 3D survey of about 1 km(2) in Lake Geneva, Switzerland, over a complex fault zone. A 5-m shot spacing resulted in a nominal fold of 6. A double-chamber bubble-cancelling 15/15 in(3) air gun (40-650 Hz) operated at 80 bars and 1 m depth gave a signal penetration of 300 m below water bottom and a best vertical resolution of 1.1 m. Processing followed a conventional scheme, but had to be adapted to the high sampling rates, and our unconventional navigation data needed conversion to industry standards. The high-quality data enabled us to construct maps of seismic horizons and fault surfaces in three dimensions. The system proves to be well adapted to investigate complex structures by providing non-aliased images of reflectors with dips up to 30 degrees.

Transcript profiling suggests that differential metabolic adaptation of mice to a high fat diet is associated with changes in liver to muscle lipid fluxes.

Genetically homogenous C57Bl/6 mice display differential metabolic adaptation when fed a high fat diet for 9 months. Most become obese and diabetic, but a significant fraction remains lean and diabetic or lean and non-diabetic. Here, we performed microarray analysis of "metabolic" transcripts expressed in liver and hindlimb muscles to evaluate: (i) whether expressed transcript patterns could indicate changes in metabolic pathways associated with the different phenotypes, (ii) how these changes differed from the early metabolic adaptation to short term high fat feeding, and (iii) whether gene classifiers could be established that were characteristic of each metabolic phenotype. Our data indicate that obesity/diabetes was associated with preserved hepatic lipogenic gene expression and increased plasma levels of very low density lipoprotein and, in muscle, with an increase in lipoprotein lipase gene expression. This suggests increased muscle fatty acid uptake, which may favor insulin resistance. In contrast, the lean mice showed a strong reduction in the expression of hepatic lipogenic genes, in particular of Scd-1, a gene linked to sensitivity to diet-induced obesity; the lean and non-diabetic mice presented an additional increased expression of eNos in liver. After 1 week of high fat feeding the liver gene expression pattern was distinct from that seen at 9 months in any of the three mouse groups, thus indicating progressive establishment of the different phenotypes. Strikingly, development of the obese phenotype involved re-expression of Scd-1 and other lipogenic genes. Finally, gene classifiers could be established that were characteristic of each metabolic phenotype. Together, these data suggest that epigenetic mechanisms influence gene expression patterns and metabolic fates.

Strong time dependence of the 76-gene prognostic signature for node-negative breast cancer patients in the TRANSBIG multicenter independent validation series.

PURPOSE: Recently, a 76-gene prognostic signature able to predict distant metastases in lymph node-negative (N(-)) breast cancer patients was reported. The aims of this study conducted by TRANSBIG were to independently validate these results and to compare the outcome with clinical risk assessment. EXPERIMENTAL DESIGN: Gene expression profiling of frozen samples from 198 N(-) systemically untreated patients was done at the Bordet Institute, blinded to clinical data and independent of Veridex. Genomic risk was defined by Veridex, blinded to clinical data. Survival analyses, done by an independent statistician, were done with the genomic risk and adjusted for the clinical risk, defined by Adjuvant! Online. RESULTS: The actual 5- and 10-year time to distant metastasis were 98% (88-100%) and 94% (83-98%), respectively, for the good profile group and 76% (68-82%) and 73% (65-79%), respectively, for the poor profile group. The actual 5- and 10-year overall survival were 98% (88-100%) and 87% (73-94%), respectively, for the good profile group and 84% (77-89%) and 72% (63-78%), respectively, for the poor profile group. We observed a strong time dependence of this signature, leading to an adjusted hazard ratio of 13.58 (1.85-99.63) and 8.20 (1.10-60.90) at 5 years and 5.11 (1.57-16.67) and 2.55 (1.07-6.10) at 10 years for time to distant metastasis and overall survival, respectively. CONCLUSION: This independent validation confirmed the performance of the 76-gene signature and adds to the growing evidence that gene expression signatures are of clinical relevance, especially for identifying patients at high risk of early distant metastases.

Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

Changes in leg-spring behavior during a 5000 m self-paced run in differently trained athletes

Summary Aims.-To explore whether fatigue-induced changes in spring-mass behavior during a 5000m self-paced run varied according to the runner's training status. Methods and results.-Six highly- and six well-trained triathletes completed a 5000m time trial. Running velocity and vertical stiffness decreased significantly (P < 0.05) with fatigue, whereas leg stiffness remained constant. None of these parameters displayed a significant interaction between fatigue and training status, despite vertical stiffness being higher (P < 0.05) in highly-trained triathletes. Conclusions.-During a 5000m self-paced run, impairments in leg-spring behavior that occur with fatigue are not affected by athletes' training status. © 2009 Elsevier Masson SAS. All rights reserved. Objectifs.-Étudier, chez des athlètes de niveaux différents, les modifications de raideur mécanique liées à l'apparition de la fatigue lors d'une course de 5000 m. Synthèse des faits.-Six triathlètes très entraînés et six autres bien entraînés ont réalisé une course de 5000 m. La vitesse de course et la raideur verticale diminuaient significativement (p < 0,05) avec la fatigue, alors que la raideur de la jambe demeurait inchangée. Aucune interaction entre la fatigue et le niveau d'entraînement n'a été détectée, malgré des niveaux de raideur verticale plus élevés (p < 0,05) chez les sujets les mieux entraînés.

Changes in leg spring behaviour, plantar loading and foot mobility magnitude induced by an exhaustive treadmill run in adolescent middle-distance runners.

OBJECTIVES: This study aimed to determine adjustments in spring-mass model characteristics, plantar loading and foot mobility induced by an exhaustive run. DESIGN: Within-participants repeated measures. METHODS: Eleven highly-trained adolescent middle-distance runners ran to exhaustion on a treadmill at a constant velocity corresponding to 95% of velocity associated with VO₂max (17.8 ± 1.4 kmh(-1), time to exhaustion=8.8 ± 3.4 min). Contact time obtained from plantar pressure sensors was used to estimate spring-mass model characteristics, which were recorded (during 30 s) 1 min after the start and prior to exhaustion using pressure insoles. Foot mobility magnitude (a composite measure of vertical and medial-lateral mobility of the midfoot) was measured before and after the run. RESULTS: Mean contact area (foot to ground), contact time, peak vertical ground reaction force, centre of mass vertical displacement and leg compression increased significantly with fatigue, while flight time, leg stiffness and mean pressure decreased. Leg stiffness decreased because leg compression increased to a larger extent than peak vertical ground reaction forces. Step length, step frequency and foot mobility magnitude did not change at exhaustion. CONCLUSIONS: The stride pattern of adolescents when running on a treadmill at high constant velocity deteriorates near exhaustion, as evidenced by impaired leg-spring behaviour (leg stiffness) and altered plantar loading.

Spatio-temporal transcript profiling of rice roots and shoots in response to phosphate starvation and recovery.

Using rice (Oryza sativa) as a model crop species, we performed an in-depth temporal transcriptome analysis, covering the early and late stages of Pi deprivation as well as Pi recovery in roots and shoots, using next-generation sequencing. Analyses of 126 paired-end RNA sequencing libraries, spanning nine time points, provided a comprehensive overview of the dynamic responses of rice to Pi stress. Differentially expressed genes were grouped into eight sets based on their responses to Pi starvation and recovery, enabling the complex signaling pathways involved in Pi homeostasis to be untangled. A reference annotation-based transcript assembly was also generated, identifying 438 unannotated loci that were differentially expressed under Pi starvation. Several genes also showed induction of unannotated splice isoforms under Pi starvation. Among these, PHOSPHATE2 (PHO2), a key regulator of Pi homeostasis, displayed a Pi starvation-induced isoform, which was associated with increased translation activity. In addition, microRNA (miRNA) expression profiles after long-term Pi starvation in roots and shoots were assessed, identifying 20 miRNA families that were not previously associated with Pi starvation, such as miR6250. In this article, we present a comprehensive spatio-temporal transcriptome analysis of plant responses to Pi stress, revealing a large number of potential key regulators of Pi homeostasis in plants.