Developmental and Environmental Regulation of Aquaporin Gene Expression across Populus Species: Divergence or Redundancy?

Aquaporins (AQPs) are membrane channels belonging to the major intrinsic proteins family and are known for their ability to facilitate water movement. While in Populus trichocarpa, AQP proteins form a large family encompassing fifty-five genes, most of the experimental work focused on a few genes or subfamilies. The current work was undertaken to develop a comprehensive picture of the whole AQP gene family in Populus species by delineating gene expression domain and distinguishing responsiveness to developmental and environmental cues. Since duplication events amplified the poplar AQP family, we addressed the question of expression redundancy between gene duplicates. On these purposes, we carried a meta-analysis of all publicly available Affymetrix experiments. Our in-silico strategy controlled for previously identified biases in cross-species transcriptomics, a necessary step for any comparative transcriptomics based on multispecies design chips. Three poplar AQPs were not supported by any expression data, even in a large collection of situations (abiotic and biotic constraints, temporal oscillations and mutants). The expression of 11 AQPs was never or poorly regulated whatever the wideness of their expression domain and their expression level. Our work highlighted that PtTIP1;4 was the most responsive gene of the AQP family. A high functional divergence between gene duplicates was detected across species and in response to tested cues, except for the root-expressed PtTIP2;3/PtTIP2;4 pair exhibiting 80% convergent responses. Our meta-analysis assessed key features of aquaporin expression which had remained hidden in single experiments, such as expression wideness, response specificity and genotype and environment interactions. By consolidating expression profiles using independent experimental series, we showed that the large expansion of AQP family in poplar was accompanied with a strong divergence of gene expression, even if some cases of functional redundancy could be suspected.

Comparative genomics of epidemic versus sporadic Staphylococcus aureus strains does not reveal molecular markers for epidemicity.

Staphylococcus aureus, especially when it is methicillin resistant, has been recognised as a major cause of nosocomial and community-acquired infections. It has also been shown that certain strains were able to cause clonal epidemics whereas others showed a more incidental occurrence. On the basis of this behavioural distinction, a genetic feature underlying this difference in epidemicity can be assumed. Understanding the difference will not only contribute to the development of markers for the identification of epidemic strains but will also shed light on the evolution of clones. Genomes of strains from two independent collections (n=18 and n=10 strains) were analysed. Both collections were composed of carefully selected, genetically diverse strains with clinically well-defined epidemic and sporadic behaviour. Comparative genome hybridisation (CGH) was performed using an Agilent array for one collection (up to 11 probes per open reading frame - ORF), and an Affymetrix array for the other (up to 30 probes per ORF). Presence and absence information of probe homologues and ORFs was taken for analysis of molecular variance (AMOVA) at the strain and behaviour levels. Not a single probe showed 100% concordant differences between epidemic and sporadic strains. Moreover, probe differences between groups were always smaller than those within groups. This was also true, when the analysis was focussed on presence versus absence of ORF's or when probe information was transformed into allelic profiles. These findings present strong evidence against the presence or absence of a single common specific genetic factor differentiating epidemic from sporadic S. aureus clones.

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.

Identification of prognostic molecular features in the reactive stroma of human breast and prostate cancer.

Primary tumor growth induces host tissue responses that are believed to support and promote tumor progression. Identification of the molecular characteristics of the tumor microenvironment and elucidation of its crosstalk with tumor cells may therefore be crucial for improving our understanding of the processes implicated in cancer progression, identifying potential therapeutic targets, and uncovering stromal gene expression signatures that may predict clinical outcome. A key issue to resolve, therefore, is whether the stromal response to tumor growth is largely a generic phenomenon, irrespective of the tumor type or whether the response reflects tumor-specific properties. To address similarity or distinction of stromal gene expression changes during cancer progression, oligonucleotide-based Affymetrix microarray technology was used to compare the transcriptomes of laser-microdissected stromal cells derived from invasive human breast and prostate carcinoma. Invasive breast and prostate cancer-associated stroma was observed to display distinct transcriptomes, with a limited number of shared genes. Interestingly, both breast and prostate tumor-specific dysregulated stromal genes were observed to cluster breast and prostate cancer patients, respectively, into two distinct groups with statistically different clinical outcomes. By contrast, a gene signature that was common to the reactive stroma of both tumor types did not have survival predictive value. Univariate Cox analysis identified genes whose expression level was most strongly associated with patient survival. Taken together, these observations suggest that the tumor microenvironment displays distinct features according to the tumor type that provides survival-predictive value.

Role of Notch signaling in epidermis and appendages

Résumé Dans la peau, il a été montré que Notch1 induit l'arrêt de la prolifération et la différentiation des keratinocytes. L'inactivation de Notch1 cause une hyperplasie de l'épiderme et la formation de carcinomes basaux cellulaires. Notre groupe a principalement identifié deux voies de signalisations, la voie Shh et la voie Wnt, qui sont dérégulées en conséquence de l'inactivation de Notch1 dans la peau. Nous avons démontré l'habilité de Notch1 à réprimer la voie Wnt induite par ß-catenin dans les keratinocytes primaires ainsi que dans d'autres types de cellules épithéliales humaines. De plus, nous avons pu déterminer que Notch1 régule cette voie, probablement en favorisant la phosphorylation de ß-catenin par le complexe axin/APC/GSK-3ß. La protéine faisant partie de la voie Wnt, ou la protéine affectant la voie Wnt, qui est régulée par Notch1 est sujette à de plus amples investigations. Un autre but de cette étude a été l'identification de potentiels gènes cibles de Notch1 autres que ceux faisant partie des voies de signalisation Shh et Wnt précédemment évoquées. Ce projet fut abordé par l'analyse de puces à ADN (ISREC et Affymetrix) qui ont été utilisées pour des expériences de gain et de perte de fonction de Notch1 dans des keratinocytes prúmaires. En plus de l'hyperplasie épidermale, les souris Notch1 déficiente ont une perte importante de poils. Nous avons montré que Notch1 est nécessaire pour le développement et l'homéostasie des follicules pileux. En effet, l'inactivation du gène Notch1 mediée par l'activation des kératines 5 ou 14 dans l'épiderme, cause des défauts du cycle ainsi que de la structure des poils. De plus, d'autres appendices de la peau, comme les glandes sudoripares et de Meibomius, ont une structure anormale et sont non fonctionnelles dans les souris Notch1 déficiente. Finalement, nous avons observé que la déficience de Notch1 dans l'épithélium cornéen mène à la formation d'une plaque épidermale opaque sur la cornée. Basé sur l'hypothèse que le défaut des glandes de Meibomius des souris Notch1 déficientes cause des lésions de la surface oculaire, nous avons montré que Notch1 est essentiel pour la cicatrisation de la cornée. Lorsque Notch1 est absent, les cellules souches de l'épithélium cornéen ne sont plus capables de se différentier en cellules cornéennes, mais réparent la blessure en se différentiant en épiderme. Ce résultat indique que Notch1 est essentiel pour la différentiation de cellules souches de la cornée qui sont spécifiquement impliquées dans la réparation de la cornée. De plus, nous avons montré que l'expression de CRBP1 dans l'épithélium cornéen est diminuée en l'absence de Notch1, ceci étant possiblement à l'origine de la formation de la plaque épidermale. Abstract: In the skin, Notch1 has been shown to trigger cell growth arrest and differentiation of keratinocytes. Notch1 inactivation results in epidermal hyperplasia and subsequent formation of basal cell carcinoma-like (BCC-like) tumors. So far our group has identified two main pathways, the Shh and the Wnt pathway, that are deregulated as a consequence of Notch1 inactivation in the skin. We showed the ability of Notch1 to represses ß-catenin-mediated Wnt signaling in primary keratinocytes as well as in other types of human epithelial cells. In addition we were able to determine that Notch1 regulates this pathway possibly by enhancing ß-catenin phosphorylation by the axin/APC/GSK-3ß complex. The exact target protein of the Wnt pathway or target protein that affects the Wnt pathway, and that is regulated by Notch1, is subject of current investigation. Another aim of this study was the identification of possible Notch1 target genes in addition to those of the Shh and Wnt signaling pathways. This was addressed by gene chip analysis using ISREC as well as Affymetrix microarrays for gain and loss of function of Notch1 in mouse primary keratinocytes. In addition to epidermal hyperplasia, Notch1 deficient mice show an important hair loss. We showed that Notch1 is required for postnatal development and homeostasis of hair follicles. Indeed, keratin5 or keratinl4-driven Cre recombinase-mediated inactivation of the Notch1 gene in the epidermis causes perturbations of the hair cycle and structural defects of the hair follicle. Moreover, other skin appendages, like the sweat and Meibomian glands show abnormal morphology and are not functional in the Notch 1 deficient mice. Finally, we observed that Notch1 deficiency in the corneal epithelium leads to the formation of an epidermal corneal plaque. Based on the hypothesis that the Meiboinian gland defect in the Notch1 deficient mice results in lesions of the eye surface, we showed that Notch1 is essential for wound-healing of the cornea. In absence of Notch1 the stem cells of the corneal epithelium are no longer able to differentiate in the corneal fate but instead repair the wound by differentiating into skin-like epidermis. This result indicated that Notch1 is essential for the differentiation of corneal stem cells specifically implicated in corneal wound-healing. Moreover, we showed that CRBP1 expression in the corneal epithelium was lost in the absence of Notch1, possibly being at the origin of plaque formation.

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.

The metal ion binding protein Cnnm4 is mutated in rod-cone dystrophy/amelogenesis imperfecta syndrome

Purpose:To identify the gene causing rod-cone dystrophy/amelogenesis imperfecta Methods:Homozygosity mapping was performed using the Affymetrix 50K XbaI array in one family and candidate genes in the linked interval were sequenced with ABI Dye Terminator, vers. 1 in the index patient of 3 families. The identified mutations were screened in normal control individuals. Expression analyses were performed on RNA extracted from the brain, various parts of the eye and teeth; immunostaining was done on mouse eyes and jaw and knock-down experiments were carried out in zebrafish embroys. Results:Sequencing the coding regions of ancient conserved domain protein 4 (CNNM4), a metal ions transporter, revealed a 1-base pair duplication (p.L438fs) in family A, a p.R236Q mutation in family B and a p.L324P in family C. All these mutations were homozygous and involved very conserved amino acids in paralogs and orthologs. Immunostaining and RT-PCR confirmed that CNNM4 was strongly expressed in various parts of the eye and in the teeth. Morpholino experiments in zebrafish showed a loss of ganglion cells at 5 days post fertilization. Conclusions:The rod-cone dystrophy/amelogenesis imperfecta syndrome is caused by mutation in CNNM4 and is due to aberrant metal ion homeostasis.

Rpe65 microarrays, from genes to phosphorylated proteins

Purpose:To describe a novel in silico method to gather and analyze data from high-throughput heterogeneous experimental procedures, i.e. gene and protein expression arrays. Methods:Each microarray is assigned to a database which handles common data (names, symbols, antibody codes, probe IDs, etc.). Links between informations are automatically generated from knowledge obtained in freely accessible databases (NCBI, Swissprot, etc). Requests can be made from any point of entry and the displayed result is fully customizable. Results:The initial database has been loaded with two sets of data: a first set of data originating from an Affymetrix-based retinal profiling performed in an RPE65 knock-out mouse model of Leber's congenital amaurosis. A second set of data generated from a Kinexus microarray experiment done on the retinas from the same mouse model has been added. Queries display wild type versus knock out expressions at several time points for both genes and proteins. Conclusions:This freely accessible database allows for easy consultation of data and facilitates data mining by integrating experimental data and biological pathways.

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.

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

Defining the genomic signature of the parous breast.

ABSTRACT: BACKGROUND: It is accepted that a woman's lifetime risk of developing breast cancer after menopause is reduced by early full term pregnancy and multiparity. This phenomenon is thought to be associated with the development and differentiation of the breast during pregnancy. METHODS: In order to understand the underlying molecular mechanisms of pregnancy induced breast cancer protection, we profiled and compared the transcriptomes of normal breast tissue biopsies from 71 parous (P) and 42 nulliparous (NP) healthy postmenopausal women using Affymetrix Human Genome U133 Plus 2.0 arrays. To validate the results, we performed real time PCR and immunohistochemistry. RESULTS: We identified 305 differentially expressed probesets (208 distinct genes). Of these, 267 probesets were up- and 38 down-regulated in parous breast samples; bioinformatics analysis using gene ontology enrichment revealed that up-regulated genes in the parous breast represented biological processes involving differentiation and development, anchoring of epithelial cells to the basement membrane, hemidesmosome and cell-substrate junction assembly, mRNA and RNA metabolic processes and RNA splicing machinery. The down-regulated genes represented biological processes that comprised cell proliferation, regulation of IGF-like growth factor receptor signaling, somatic stem cell maintenance, muscle cell differentiation and apoptosis. CONCLUSIONS: This study suggests that the differentiation of the breast imprints a genomic signature that is centered in the mRNA processing reactome. These findings indicate that pregnancy may induce a safeguard mechanism at post-transcriptional level that maintains the fidelity of the transcriptional process.

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

Characterization of a genomic signature of pregnancy identified in the breast.

The objective of this study was to comprehensively compare the genomic profiles in the breast of parous and nulliparous postmenopausal women to identify genes that permanently change their expression following pregnancy. The study was designed as a two-phase approach. In the discovery phase, we compared breast genomic profiles of 37 parous with 18 nulliparous postmenopausal women. In the validation phase, confirmation of the genomic patterns observed in the discovery phase was sought in an independent set of 30 parous and 22 nulliparous postmenopausal women. RNA was hybridized to Affymetrix HG_U133 Plus 2.0 oligonucleotide arrays containing probes to 54,675 transcripts, scanned and the images analyzed using Affymetrix GCOS software. Surrogate variable analysis, logistic regression, and significance analysis of microarrays were used to identify statistically significant differences in expression of genes. The false discovery rate (FDR) approach was used to control for multiple comparisons. We found that 208 genes (305 probe sets) were differentially expressed between parous and nulliparous women in both discovery and validation phases of the study at an FDR of 10% and with at least a 1.25-fold change. These genes are involved in regulation of transcription, centrosome organization, RNA splicing, cell-cycle control, adhesion, and differentiation. The results provide initial evidence that full-term pregnancy induces long-term genomic changes in the breast. The genomic signature of pregnancy could be used as an intermediate marker to assess potential chemopreventive interventions with hormones mimicking the effects of pregnancy for prevention of breast cancer.

Concordance among digital gene expression, microarrays, and qPCR when measuring differential expression of microRNAs.

Profiling microRNA (miRNA) expression is of widespread interest given the critical role of miRNAs in many cellular functions. Profiling can be achieved via hybridization-based (microarrays), sequencing-based, or amplification-based (quantitative reverse transcription-PCR, qPCR) technologies. Among these, microarrays face the significant challenge of accurately distinguishing between mature and immature miRNA forms, and different vendors have developed different methods to meet this challenge. Here we measure differential miRNA expression using the Affymetrix, Agilent, and Illumina microarray platforms, as well as qPCR (Applied Biosystems) and ultra high-throughput sequencing (Illumina). We show that the differential expression measurements are more divergent when the three types of microarrays are compared than when the Agilent microarray, qPCR, and sequencing technology measurements are compared, which exhibit a good overall concordance.

Tumor-host interactions Part 1: Stromal signatures of breast and prostate cancer Part 2: GLG1 and the control of the secretory pathway in tumor cells

Tumor-host interaction is a key determinant during cancer progression, from primary tumor growth to metastatic dissemination. At each step, tumor cells have to adapt to and subvert different types of microenvironment, leading to major phenotypic and genotypic alterations that affect both tumor and surrounding stromal compartments. Understanding the molecular mechanisms that govern tumor-host interplay may be essential for better comprehension of tumorigenesis in an effort to improve current anti-cancer therapies. The present work is composed of two projects that address tumor-host interactions from two different perspectives, the first focusing on the characterization of tumor-associated stroma and the second on membrane trafficking in tumor cells. Part 1. To selectively address stromal gene expression changes during cancer progression, oligonucleotide-based Affymetrix microarray technology was used to analyze the transcriptomes of laser-microdissected stromal cells derived from invasive human breast and prostate carcinoma. Comparison showed that invasive breast and prostate cancer elicit distinct, tumor-specific stromal responses, with a limited panel of shared induced and/or repressed genes. Both breast and prostate tumor-specific deregulated stromal gene sets displayed statistically significant survival-predictive ability for their respective tumor type. By contrast, a stromal gene signature common to both tumor types did not display prognostic value, although expression of two individual genes within this common signature was found to be associated with patient survival. Part 2. GLG1 is known as an E-selectin ligand and an intracellular FGF receptor, depending on cell type and context. Immunohistochemical and immunofluorescence analyses showed that GLG1 is primarily localized in the Golgi of human tumor cells, a central location in the biosynthetic/secretory pathways. GLG1 has been shown to interact with and to recruit the ARF GEF BIGI to the Golgi membrane. Depletion of GLG1 or BIGI markedly reduced ARF3 membrane localization and activation, and altered the Golgi structure. Interestingly, these perturbations did not impair constitutive secretion in general, but rather seemed to impair secretion of a specific subset of proteins that includes MMP-9. Thus, GLG1 coordinates ARF3 activation by recruiting BIGI to the Golgi membrane, thereby affecting secretion of specific molecules. - Les interactions tumeur-hôte constituent un élément essentiel à la progression tumorale, de la croissance de la tumeur primaire à la dissémination des métastases. A chaque étape, les cellules tumorales doivent s'adapter à différents types de microenvironnement et les détourner à leur propre avantage, donnant lieu à des altérations phénotypiques et génotypiques majeures qui affectent aussi bien la tumeur elle-même que le compartiment stromal environnant. L'étude des mécanismes moléculaires qui régissent les interactions tumeur-hôte constitue une étape essentielle pour une meilleure compréhension du processus de tumorigenèse dans le but d'améliorer les thérapies anti cancer existantes. Le travail présenté ici est composé de deux projets qui abordent la problématique des interactions tumeur-hôte selon différentes perspectives, le premier se concentrant sur la caractérisation du stroma tumoral et le second sur le trafic intracellulaire des cellules tumorales. Partie 1. Pour examiner les changements d'expression des gènes dans le stroma en réponse à la progression du cancer, des puces à ADN Affymetrix ont été utilisées afin d'analyser les transcriptomes des cellules stromales issues de carcinomes invasifs du sein et de la prostate et collectées par microdissection au laser. L'analyse comparative a montré que les cancers invasifs du sein et de la prostate provoquent des réponses stromales spécifiques à chaque type de tumeur, et présentent peu de gènes induits ou réprimés de façon similaire. L'ensemble des gènes dérégulés dans le stroma associé au cancer du sein, ou à celui de la prostate, présente une valeur pronostique pour les patients atteints d'un cancer du sein, respectivement de la prostate. En revanche, la signature stromale commune aux deux types de cancer n'a aucune valeur prédictive, malgré le fait que l'expression de deux gènes présents dans cette liste soit liée à la survie des patients. Partie 2. GLG1 est connu comme un ligand des sélectines E ainsi que comme récepteur intracellulaire pour des facteurs de croissances FGFs selon le type de cellule dans lequel il est exprimé. Des analyses immunohistochimiques et d'immunofluorescence ont montré que dans les cellules tumorales, GLG1 est principalement localisé au niveau de l'appareil de Golgi, une place centrale dans la voie biosynthétique et sécrétoire. Nous avons montré que GLG1 interagit avec la protéine BIGI et participe à son recrutement à la membrane du Golgi. L'absence de GLG1 ou de BIGI réduit drastiquement le pool d'ARF3 associé aux membranes ainsi que la quantité d'ARF3 activés, et modifie la structure de l'appareil de Golgi. Il est particulièrement intéressant de constater que ces perturbations n'ont pas d'effet sur la sécrétion constitutive en général, mais semblent plutôt affecter la sécrétion spécifique d'un sous-groupe défini de protéines comprenant MMP-9. GLG1 coordonne donc l'activation de ARF3 en recrutant BIGI à la membrane du Golgi, agissant par ce moyen sur la sécrétion de molécules spécifiques.