Mutant huntingtin's effects on striatal gene expression in mice recapitulate changes observed in human Huntington's disease brain and do not differ with mutant huntingtin length or wild-type huntingtin dosage.

To test the hypotheses that mutant huntingtin protein length and wild-type huntingtin dosage have important effects on disease-related transcriptional dysfunction, we compared the changes in mRNA in seven genetic mouse models of Huntington's disease (HD) and postmortem human HD caudate. Transgenic models expressing short N-terminal fragments of mutant huntingtin (R6/1 and R6/2 mice) exhibited the most rapid effects on gene expression, consistent with previous studies. Although changes in the brains of knock-in and full-length transgenic models of HD took longer to appear, 15- and 22-month CHL2(Q150/Q150), 18-month Hdh(Q92/Q92) and 2-year-old YAC128 animals also exhibited significant HD-like mRNA signatures. Whereas it was expected that the expression of full-length huntingtin transprotein might result in unique gene expression changes compared with those caused by the expression of an N-terminal huntingtin fragment, no discernable differences between full-length and fragment models were detected. In addition, very high correlations between the signatures of mice expressing normal levels of wild-type huntingtin and mice in which the wild-type protein is absent suggest a limited effect of the wild-type protein to change basal gene expression or to influence the qualitative disease-related effect of mutant huntingtin. The combined analysis of mouse and human HD transcriptomes provides important temporal and mechanistic insights into the process by which mutant huntingtin kills striatal neurons. In addition, the discovery that several available lines of HD mice faithfully recapitulate the gene expression signature of the human disorder provides a novel aspect of validation with respect to their use in preclinical therapeutic trials.

Integrin-mediated adhesion and soluble ligand binding stabilize COX-2 protein levels in endothelial cells by inducing expression and preventing degradation.

Cyclooxygenase-2 (COX-2), a key enzyme in prostaglandin synthesis, is highly expressed during inflammation and cellular transformation and promotes tumor progression and angiogenesis. We have previously demonstrated that endothelial cell COX-2 is required for integrin alphaVbeta3-dependent activation of Rac-1 and Cdc-42 and for endothelial cell spreading, migration, and angiogenesis (Dormond, O., Foletti, A., Paroz, C., and Ruegg, C. (2001) Nat. Med. 7, 1041-1047; Dormond, O., Bezzi, M., Mariotti, A., and Ruegg, C. (2002) J. Biol. Chem. 277, 45838-45846). In this study, we addressed the question of whether integrin-mediated cell adhesion may regulate COX-2 expression in endothelial cells. We report that cell detachment from the substrate caused rapid degradation of COX-2 protein in human umbilical vein endothelial cells (HUVEC) independent of serum stimulation. This effect was prevented by broad inhibition of cellular proteinases and by neutralizing lysosomal activity but not by inhibiting the proteasome. HUVEC adhesion to laminin, collagen I, fibronectin, or vitronectin induced rapid COX-2 protein expression with peak levels reached within 2 h and increased COX-2-dependent prostaglandin E2 production. In contrast, nonspecific adhesion to poly-L-lysine was ineffective in inducing COX-2 expression. Furthermore, the addition of matrix proteins in solution promoted COX-2 protein expression in suspended or poly-L-lysine-attached HUVEC. Adhesion-induced COX-2 expression was strongly suppressed by pharmacological inhibition of c-Src, phosphatidylinositol 3-kinase, p38, extracellular-regulated kinase 1/2, and, to a lesser extent, protein kinase C and by the inhibition of mRNA or protein synthesis. In conclusion, this work demonstrates that integrin-mediated cell adhesion and soluble integrin ligands contribute to maintaining COX-2 steady-state levels in endothelial cells by the combined prevention of lysosomal-dependent degradation and the stimulation of mRNA synthesis involving multiple signaling pathways.

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.

Transgenic expression of Ly49A on T cells impairs a specific antitumor response.

Inhibitory MHC receptors determine the reactivity and specificity of NK cells. These receptors can also regulate T cells by modulating TCR-induced effector functions such as cytotoxicity, cytokine production, and proliferation. Here we have assessed the capacity of mouse T cells expressing the inhibitory MHC class I receptor Ly49A to respond to a well-defined tumor Ag in vivo using Ly49A transgenic mice. We find that the presence of Ly49A on the vast majority of lymphocytes prevents the development of a significant Ag-specific CD8+ T cell response and, consequently, the rejection of the tumor. Despite minor alterations in the TCR repertoire of CD8+ T cells in the transgenic lines, precursors of functional tumor-specific CD8+ T cells exist but could not be activated most likely due to a lack of appropriate CD4+ T cell help. Surprisingly, all of these effects are observed in the absence of a known ligand for the Ly49A receptor as defined by its ability to regulate NK cell function. Indeed, we found that the above effects on T cells may be based on a weak interaction of Ly49A with Kb or Db class I molecules. Thus, our data demonstrate that enforced expression of a Ly49A receptor on conventional T cells prevents a specific immune response in vivo and suggest that the functions of T and NK cells are differentially sensitive to the presence of inhibitory MHC class I receptors.

Emergence of secretion-defective sublines of Pseudomonas aeruginosa PAO1 resulting from spontaneous mutations in the vfr global regulatory gene.

Pseudomonas aeruginosa undergoes spontaneous mutation that impairs secretion of several extracellular enzymes during extended cultivation in vitro in rich media, as well as during long-term colonization of the cystic fibrosis lung. A frequent type of strong secretion deficiency is caused by inactivation of the quorum-sensing regulatory gene lasR. Here we analyzed a spontaneously emerging subline of strain PAO1 that exhibited moderate secretion deficiency and partial loss of quorum-sensing control. Using generalized transduction, we mapped the secretion defect to the vfr gene, which is known to control positively the expression of the lasR gene and type II secretion of several proteases. We confirmed this secretion defect by sequencing and complementation of the vfr mutation. In a reconstruction experiment conducted with a 1:1 mixture of wild-type strain PAO1 and a vfr mutant of PAO1, we observed that the vfr mutant had a selective advantage over the wild type after growth in static culture for 4 days. Under these conditions, spontaneous vfr emerged in a strain PAO1 population after four growth cycles, and these mutants accounted for more than 40% of the population after seven cycles. These results suggest that partial or complete loss of quorum sensing and secretion can be beneficial to P. aeruginosa under certain environmental conditions.

Validated prediction of clinical outcome in sarcomas and multiple types of cancer on the basis of a gene expression signature related to genome complexity.

Sarcomas are heterogeneous and aggressive mesenchymal tumors. Histological grading has so far been the best predictor for metastasis-free survival, but it has several limitations, such as moderate reproducibility and poor prognostic value for some histological types. To improve patient grading, we performed genomic and expression profiling in a training set of 183 sarcomas and established a prognostic gene expression signature, complexity index in sarcomas (CINSARC), composed of 67 genes related to mitosis and chromosome management. In a multivariate analysis, CINSARC predicts metastasis outcome in the training set and in an independent 127 sarcomas validation set. It is superior to the Fédération Francaise des Centres de Lutte Contre le Cancer grading system in determining metastatic outcome for sarcoma patients. Furthermore, it also predicts outcome for gastrointestinal stromal tumors (GISTs), breast carcinomas and lymphomas. Application of the signature will permit more selective use of adjuvant therapies for people with sarcomas, leading to decreased iatrogenic morbidity and improved outcomes for such individuals.

Pro-inflammatory gene expression and neurotoxic effects of activated microglia are attenuated by absence of CCAAT/enhancer binding protein beta

Background. Microglia and astrocytes respond to homeostatic disturbances with profound changes of gene expression. This response, known as glial activation or neuroinflammation, can be detrimental to the surrounding tissue. The transcription factor CCAAT/enhancer binding protein ß (C/EBPß) is an important regulator of gene expression in inflammation but little is known about its involvement in glial activation. To explore the functional role of C/EBPß in glial activation we have analyzed pro-inflammatory gene expression and neurotoxicity in murine wild type and C/EBPß-null glial cultures. Methods. Due to fertility and mortality problems associated with the C/EBPß-null genotype we developed a protocol to prepare mixed glial cultures from cerebral cortex of a single mouse embryo with high yield. Wild-type and C/EBPß-null glial cultures were compared in terms of total cell density by Hoechst-33258 staining; microglial content by CD11b immunocytochemistry; astroglial content by GFAP western blot; gene expression by quantitative real-time PCR, western blot, immunocytochemistry and Griess reaction; and microglial neurotoxicity by estimating MAP2 content in neuronal/microglial cocultures. C/EBPß DNA binding activity was evaluated by electrophoretic mobility shift assay and quantitative chromatin immunoprecipitation. Results. C/EBPß mRNA and protein levels, as well as DNA binding, were increased in glial cultures by treatment with lipopolysaccharide (LPS) or LPS + interferon ¿ (IFN¿). Quantitative chromatin immunoprecipitation showed binding of C/EBPß to pro-inflammatory gene promoters in glial activation in a stimulus- and gene-dependent manner. In agreement with these results, LPS and LPS+IFN¿ induced different transcriptional patterns between pro-inflammatory cytokines and NO synthase-2 genes. Furthermore, the expressions of IL-1ß and NO synthase-2, and consequent NO production, were reduced in the absence of C/EBPß. In addition, neurotoxicity elicited by LPS+IFN¿-treated microglia co-cultured with neurons was completely abolished by the absence of C/EBPß in microglia.

Characterization of Arabidopsis FPS isozymes and FPS gene expression analysis provide insight into the biosynthesis of isoprenoid precursors in seeds

Arabidopsis thaliana contains two genes encoding farnesyl diphosphate (FPP) synthase (FPS), the prenyl diphoshate synthase that catalyzes the synthesis of FPP from isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). In this study, we provide evidence that the two Arabidopsis short FPS isozymes FPS1S and FPS2 localize to the cytosol. Both enzymes were expressed in E. coli, purified and biochemically characterized. Despite FPS1S and FPS2 share more than 90% amino acid sequence identity, FPS2 was found to be more efficient as a catalyst, more sensitive to the inhibitory effect of NaCl, and more resistant to thermal inactivation than FPS1S. Homology modelling for FPS1S and FPS2 and analysis of the amino acid differences between the two enzymes revealed an increase in surface polarity and a greater capacity to form surface salt bridges of FPS2 compared to FPS1S. These factors most likely account for the enhanced thermostability of FPS2. Expression analysis of FPS::GUS genes in seeds showed that FPS1 and FPS2 display complementary patterns of expression particularly at late stages of seed development, which suggests that Arabidopsis seeds have two spatially segregated sources of FPP. Functional complementation studies of the Arabidopsis fps2 knockout mutant seed phenotypes demonstrated that under normal conditions FPS1S and FPS2 are functionally interchangeable. A putative role for FPS2 in maintaining seed germination capacity under adverse environmental conditions is discussed.

T cell receptor beta chain gene rearrangements and V beta gene usage in horse cytochrome c-specific T cell hybridomas.

Comparison of T cell receptor alpha and beta-chain genes in murine major histocompatibility complex (MHC) class I and class II-restricted T cell clones and hybridomas recognizing different antigens indicates that no simple correlation exists between the observed antigen/MHC specificity and the expression of certain alpha and beta-chain heterodimers. We have attempted to establish a possible correlation by analyzing T cell receptor beta chain gene rearrangements and V beta gene usage in five T cell hybridomas with identical antigen/MHC specificity and another hybridoma recognizing a different antigenic determinant in association with the same restriction molecule. We report here that in each of the five clones a uniquely rearranged beta chain gene is expressed in combination with at least two different V beta gene segments. The presence of the differently rearranged T cell receptor beta chain genes correlated with the finding of distinct fine specificity pattern of antigen recognition in each of the hybridomas. Interestingly, two hybridomas specific for different epitopes showed identical beta chain D-J rearrangements indicating that the differences might be encoded by the alpha chain gene or/and the V beta gene element.

Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling

Background Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life. Results The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR. Conclusions Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.

MFG-E8/lactadherin regulates cyclins D1/D3 expression and enhances the tumorigenic potential of mammary epithelial cells.

Milk fat globule-EGF factor 8 (MFG-E8) is a glycoprotein highly expressed in breast cancer that contributes to tumor progression through largely undefined mechanisms. By analyzing publicly available gene expression profiles of breast carcinomas, we found that MFG-E8 is highly expressed in primary and metastatic breast carcinomas, associated with absent estrogen receptor expression. Immunohistochemistry analysis of breast cancer biopsies revealed that MFG-E8 is expressed on the cell membrane as well as in the cytoplasm and nucleus. We also show that increased expression of MFG-E8 in mammary carcinoma cells increases their tumorigenicity in immunodeficient mice, and conversely, its downregulation reduces their in vivo growth. Moreover, expression of MFG-E8 in immortalized mammary epithelial cells promotes their growth and branching in three-dimensional collagen matrices and induces the expression of cyclins D1/D3 and N-cadherin. A mutant protein unable to bind integrins can in part exert these effects, indicating that MFG-E8 function is only partially dependent on integrin activation. We conclude that MFG-E8-dependent signaling stimulates cell proliferation and the acquisition of mesenchymal properties and contributes to mammary carcinoma development.

Human natural Treg microRNA signature: role of microRNA-31 and microRNA-21 in FOXP3 expression.

Treg are the main mediators of dominant tolerance. Their mechanisms of action and applications are subjects of considerable debate currently. However, a human microRNA (miR) Treg signature has not been described yet. We investigated human natural Treg and identified a signature composed of five miR (21, 31, 125a, 181c and 374). Among those, two were considerably under-expressed (miR-31 and miR-125a). We identified a functional target sequence for miR-31 in the 3' untranslated region (3' UTR) of FOXP3 mRNA. Using lentiviral transduction of fresh cord blood T cells, we demonstrated that miR-31 and miR-21 had an effect on FOXP3 expression levels. We showed that miR-31 negatively regulates FOXP3 expression by binding directly to its potential target site in the 3' UTR of FOXP3 mRNA. We next demonstrated that miR-21 acted as a positive, though indirect, regulator of FOXP3 expression. Transduction of the remaining three miR had no direct effect on FOXP3 expression or on the phenotype and will remain the subject of future investigations. In conclusion, not only have we identified and validated a miR signature for human natural Treg, but also unveiled some of the mechanisms by which this signature was related to the control of FOXP3 expression in these cells.

Inhibition of estrogen-responsive gene activation by the retinoid X receptor beta: evidence for multiple inhibitory pathways.

The retinoid X receptor beta (RXR beta; H-2RIIBP) forms heterodimers with various nuclear hormone receptors and binds multiple hormone response elements, including the estrogen response element (ERE). In this report, we show that endogenous RXR beta contributes to ERE binding activity in nuclear extracts of the human breast cancer cell line MCF-7. To define a possible regulatory role of RXR beta regarding estrogen-responsive transcription in breast cancer cells, RXR beta and a reporter gene driven by the vitellogenin A2 ERE were transfected into estrogen-treated MCF-7 cells. RXR beta inhibited ERE-driven reporter activity in a dose-dependent and element-specific fashion. This inhibition occurred in the absence of the RXR ligand 9-cis retinoic acid. The RXR beta-induced inhibition was specific for estrogen receptor (ER)-mediated ERE activation because inhibition was observed in ER-negative MDA-MB-231 cells only following transfection of the estrogen-activated ER. No inhibition of the basal reporter activity was observed. The inhibition was not caused by simple competition of RXR beta with the ER for ERE binding, since deletion mutants retaining DNA binding activity but lacking the N-terminal or C-terminal domain failed to inhibit reporter activity. In addition, cross-linking studies indicated the presence of an auxiliary nuclear factor present in MCF-7 cells that contributed to RXR beta binding of the ERE. Studies using known heterodimerization partners of RXR beta confirmed that RXR beta/triiodothyronine receptor alpha heterodimers avidly bind the ERE but revealed the existence of another triiodothyronine-independent pathway of ERE inhibition. These results indicate that estrogen-responsive genes may be negatively regulated by RXR beta through two distinct pathways.

Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

Characterization of the PHO1 gene family and the responses to phosphate deficiency of Physcomitrella patens.

PHO1 was previously identified in Arabidopsis (Arabidopsis thaliana) as a protein involved in loading inorganic phosphate (Pi) into the xylem of roots and its expression was associated with the vascular cylinder. Seven genes homologous to AtPHO1 (PpPHO1;1-PpPHO1;7) have been identified in the moss Physcomitrella patens. The corresponding proteins harbor an SPX tripartite domain in the N-terminal hydrophilic portion and an EXS domain in the conserved C-terminal hydrophobic portion, both common features of the plant PHO1 family. Northern-blot analysis showed distinct expression patterns for the PpPHO1 genes, both at the tissue level and in response to phosphate deficiency. Transgenic P. patens expressing the beta-glucuronidase reporter gene under three different PpPHO1 promoters revealed distinct expression profiles in various tissues. Expression of PpPHO1;1 and PpPHO1;7 was specifically induced by Pi starvation. P. patens homologs to the Arabidopsis PHT1, DGD2, SQD1, and APS1 genes also responded to Pi deficiency by increased mRNA levels. Morphological changes associated with Pi deficiency included elongation of caulonemata with inhibition of the formation of side branches, resulting in colonies with greater diameter, but reduced mass compared to Pi-sufficient plants. Under Pi-deficient conditions, P. patens also increased the synthesis of ribonucleases and of an acid phosphatase, and increased the ratio of sulfolipids over phospholipids. These results indicate that P. patens and higher plants share some common strategies to adapt to Pi deficiency, although morphological changes are distinct, and that the PHO1 proteins are well conserved in bryophyte despite the lack of a developed vascular system.