Études de réseaux d’expression génique : utilité pour l’élucidation des déterminants génétiques des traits complexes

Les traits quantitatifs complexes sont des caractéristiques mesurables d’organismes vivants qui résultent de l’interaction entre plusieurs gènes et facteurs environnementaux. Les locus génétiques liés à un caractère complexe sont appelés «locus de traits quantitatifs » (QTL). Récemment, en considérant les niveaux d’expression tissulaire de milliers de gènes comme des traits quantitatifs, il est devenu possible de détecter des «QTLs d’expression» (eQTL). Alors que ces derniers ont été considérés comme des phénotypes intermédiaires permettant de mieux comprendre l’architecture biologique des traits complexes, la majorité des études visent encore à identifier une mutation causale dans un seul gène. Cette approche ne peut remporter du succès que dans les situations où le gène incriminé a un effet majeur sur le trait complexe, et ne permet donc pas d’élucider les situations où les traits complexes résultent d’interactions entre divers gènes. Cette thèse propose une approche plus globale pour : 1) tenir compte des multiples interactions possibles entre gènes pour la détection de eQTLs et 2) considérer comment des polymorphismes affectant l’expression de plusieurs gènes au sein de groupes de co-expression pourraient contribuer à des caractères quantitatifs complexes. Nos contributions sont les suivantes : Nous avons développé un outil informatique utilisant des méthodes d’analyse multivariées pour détecter des eQTLs et avons montré que cet outil augmente la sensibilité de détection d’une classe particulière de eQTLs. Sur la base d’analyses de données d’expression de gènes dans des tissus de souris recombinantes consanguines, nous avons montré que certains polymorphismes peuvent affecter l’expression de plusieurs gènes au sein de domaines géniques de co-expression. En combinant des études de détection de eQTLs avec des techniques d’analyse de réseaux de co-expression de gènes dans des souches de souris recombinantes consanguines, nous avons montré qu’un locus génétique pouvait être lié à la fois à l’expression de plusieurs gènes au niveau d’un domaine génique de co-expression et à un trait complexe particulier (c.-à-d. la masse du ventricule cardiaque gauche). Au total, nos études nous ont permis de détecter plusieurs mécanismes par lesquels des polymorphismes génétiques peuvent être liés à l’expression de plusieurs gènes, ces derniers pouvant eux-mêmes être liés à des traits quantitatifs complexes.

Heterologous expression of equine CYP3A94 and investigation of a tunable system to regulate co-expressed NADPH P450 oxidoreductase levels.

The activity of cytochrome P450 enzymes depends on the enzyme NADPH P450 oxidoreductase (POR). The aim of this study was to investigate the activity of the equine CYP3A94 using a system that allows to regulate the POR protein levels in mammalian cells. CYP3A94 and the equine POR were heterologously expressed in V79 cells. In the system used, the POR protein regulation is based on a destabilizing domain (DD) that transfers its instability to a fused protein. The resulting fusion protein is therefore degraded by the ubiquitin-proteasome system (UPS). Addition of "Shield-1" prevents the DD fusion protein from degradation. The change of POR levels at different Shield-1 concentrations was demonstrated by cytochrome c reduction, Western immunoblot analysis, and immunocytochemistry. The alteration of CYP3A94 activity was investigated using a substrate (BFC) known to detect CYP3A4 activity. Equine CYP3A94 was demonstrated to be metabolically active and its activity could be significantly elevated by co-expression of POR. Cytochrome c reduction was significantly increased in V79-CYP3A94/DD-POR cells compared to V79-CYP3A94 cells. Surprisingly, incubation with different Shield-1 concentrations resulted in a decrease in POR protein shown by Western immunoblot analysis. Cytochrome c reduction did not change significantly, but the CYP3A94 activity decreased more than 4-fold after incubation with 500 nM and 1 µM Shield-1 for 24 hours. No differences were obtained when V79-CYP3A94 POR cells with and without Shield-1 were compared. The basal activity levels of V79-CYP3A94/DD-POR cells were unexpectedly high, indicating that DD/POR is not degraded without Shield-1. Shield-1 decreased POR protein levels and CYP3A94 activity suggesting that Shield-1 might impair POR activity by an unknown mechanism. Although regulation of POR with the pPTuner system could not be obtained, the cell line V79-CYP3A94/DD-POR system can be used for further experiments to characterize the equine CYP3A94 since the CYP activity was significantly enhanced with co-expressed POR.

TaNF-YC11, one of the light-upregulated NF-YC members in Triticum aestivum, is co-regulated with photosynthesis-related genes

NF-Y is a heterotrimeric transcription factor complex. Each of the NF-Y subunits (NF-YA, NF-YB and NF-YC) in plants is encoded by multiple genes. Quantitative RT-PCR analysis revealed that five wheat NF-YC members (TaNF-YC5, 8, 9, 11 & 12) were upregulated by light in both the leaf and seedling shoot. Co-expression analysis of Affymetrix wheat genome array datasets revealed that transcript levels of a large number of genes were consistently correlated with those of the TaNF-YC11 and TaNF-YC8 genes in 3-4 separate Affymetrix array datasets. TaNF-YC11-correlated transcripts were significantly enriched with the Gene Ontology term photosynthesis. Sequence analysis in the promoters of TaNF-YC11-correlated genes revealed the presence of putative NF-Y complex binding sites (CCAAT motifs). Quantitative RT-PCR analysis of a subset of potential TaNF-YC11 target genes showed that ten out of the thirteen genes were also light-upregulated in both the leaf and seedling shoot and had significantly correlated expression profiles with TaNF-YC11. The potential target genes for TaNF-YC11 include subunit members from all four thylakoid membrane bound complexes required for the conversion of solar energy into chemical energy and rate limiting enzymes in the Calvin cycle. These data indicate that TaNF-YC11 is potentially involved in regulation of photosynthesis-related genes.

Inverse expression states of the BRN2 and MITF transcription factors in melanoma spheres and tumour xenografts regulate the NOTCH pathway

The use of adherent monolayer cultures have produced many insights into melanoma cell growth and differentiation, but often novel therapeutics demonstrated to act on these cells are not active in vivo. It is imperative that new methods of growing melanoma cells that reflect growth in vivo are investigated. To this end, a range of human melanoma cell lines passaged as adherent cultures or induced to form melanoma spheres (melanospheres) in stem cell media have been studied to compare cellular characteristics and protein expression. Melanoma spheres and tumours grown from cell lines as mouse xenografts had increased heterogeneity when compared with adherent cells and 3D-spheroids in agar (aggregates). Furthermore, cells within the melanoma spheres and mouse xenografts each displayed a high level of reciprocal BRN2 or MITF expression, which matched more closely the pattern seen in human melanoma tumours in situ, rather than the propensity for co-expression of these important melanocytic transcription factors seen in adherent cells and 3D-spheroids. Notably, when the levels of the BRN2 and MITF proteins were each independently repressed using siRNA treatment of adherent melanoma cells, members of the NOTCH pathway responded by decreasing or increasing expression, respectively. This links BRN2 as an activator, and conversely, MITF as a repressor of the NOTCH pathway in melanoma cells. Loss of the BRN2-MITF axis in antisense-ablated cell lines decreased the melanoma sphere-forming capability, cell adhesion during 3D-spheroid formation and invasion through a collagen matrix. Combined, this evidence suggests that the melanoma sphere-culture system induces subpopulations of cells that may more accurately portray the in vivo disease, than the growth as adherent melanoma cells.

Expression of PTRF in PC3 cells modulates cholesterol dynamics and the actin cytoskeleton impacting secretion pathways

Expression of caveolin-1 is up-regulated in prostate cancer metastasis and is associated with aggressive recurrence of the disease. Intriguingly, caveolin-1 is also secreted from prostate cancer cell lines and has been identified in secreted prostasomes. Caveolin-1 is the major structural component of the plasma membrane invaginations called caveolae. Co-expression of the coat protein Polymerase I and transcript release factor (PTRF) is required for caveolae formation. We recently found that expression of caveolin-1 in the aggressive prostate cancer cell line PC-3 is not accompanied by PTRF, leading to noncaveolar caveolin-1 lipid rafts. Moreover, ectopic expression of PTRF in PC-3 cells sequesters caveolin-1 into caveolae. Here we quantitatively analyzed the effect of PTRF expression on the PC-3 proteome using stable isotope labeling by amino acids in culture and subcellular proteomics. We show that PTRF reduced the secretion of a subset of proteins including secreted proteases, cytokines, and growth regulatory proteins, partly via a reduction in prostasome secretion. To determine the cellular mechanism accounting for the observed reduction in secreted proteins we analyzed total membrane and the detergent-resistant membrane fractions. Our data show that PTRF expression selectively impaired the recruitment of actin cytoskeletal proteins to the detergent-resistant membrane, which correlated with altered cholesterol distribution in PC-3 cells expressing PTRF. Consistent with this, modulating cellular cholesterol altered the actin cytoskeleton and protein secretion in PC-3 cells. Intriguingly, several proteins that function in ER to Golgi trafficking were reduced by PTRF expression. Taken together, these results suggest that the noncaveolar caveolin-1 found in prostate cancer cells generates a lipid raft microenvironment that accentuates secretion pathways, possibly at the step of ER sorting/exit. Importantly, these effects could be modulated by PTRF expression.

A transcription factor map as revealed by a genome-wide gene expression analysis of whole-blood mRNA transcriptome in multiple sclerosis

Background Several lines of evidence suggests that transcription factors are involved in the pathogenesis of Multiple Sclerosis (MS) but a complete mapping the whole network has been elusive. One of the reasons is that there are several clinical subtypes of MS and transcription factors which may be involved in one subtype may not be in others. We investigated the possibility that this network could be mapped using microarray technologies and modern bioinformatics methods on a dataset from whole blood in 99 untreated MS patients (36 Relapse Remitting MS, 43 Primary Progressive MS, and 20 Secondary Progressive MS) and 45 age-matched healthy controls, Methodology/Principal Findings We have used two different analytical methodologies: a differential expression analysis and a differential co-expression analysis, which have converged on a significant number of regulatory motifs that seem to be statistically overrepresented in genes which are either differentially expressed (or differentially co-expressed) in cases and controls (e.g. V$KROX_Q6, p-value < 3.31E-6; V$CREBP1_Q2, p-value < 9.93E-6, V$YY1_02, p-value < 1.65E-5). Conclusions/significance: Our analysis uncovered a network of transcription factors that potentially dysregulate several genes in MS or one or more of its disease subtypes. Analysing the published literature we have found that these transcription factors are involved in the early T-lymphocyte specification and commitment as well as in oligodendrocytes dedifferentiation and development. The most significant transcription factors motifs were for the Early Growth response EGR/KROX family, ATF2, YY1 (Yin and Yang 1), E2F-1/DP-1 and E2F-4/DP-2 heterodimers, SOX5, and CREB and ATF families.

Expression and prognostic influence of MMP-2, MMP-9 and TIMP-2 in non-small cell lung cancer

Background Matrix metalloproteinases (MMPs) are a family of endopeptidases that digest the extracellular matrix (ECM). Overexpression of different MMPs has been shown to promote tumour cell invasion in vitro. Tissue inhibitors of matrix metalloproteinases (TIMPs) are specific inhibitors of MMPs that also possess growth-promoting properties. Aims To analyse the expression profile of MMP-2, MMP-9 and TIMP-2 in non-small cell lung cancer (NSCLC) and to assess the impact of expression on survival. Methods This is a retrospective study of patients who underwent resection for stage I-IIIa NSCLC with a post-operative survival >60 days. Patient follow up was a minimum of 2 years. Standard ABC immunohistochemistry was performed on 4μm paraffin-embedded sections from the tumour periphery using monoclonal antibodies to MMP-2, MMP-9 and TIMP-2. Results The results of the immunohistochemistry are set out below. marker tumour expression log-rank survival stromal expression log-rank survival MMP-2 9/72 (13%) p=0.10 34/72 (47%) p=0.34 MMP-9 79/152 (52%) p=0.04* 69/152 (45%) p=0.84 TIMP-2 28/90 (31%) p=0.04* 66/90 (73%) p=0.90 Two or more 16/59 (27%) p=0.007* There were no associations between expression and clinicopathological findings for any tumour marker. There was co-expression of MMP-2 and MMP-9 in tumour cells (p=0.01). Conclusions MMP-2, MMP-9 and TIMP-2 are expressed in NSCLC. MMP-9 and TIMP-2 tumour expression correlate with a poor outcome (both p=0.04) and are potential prognostic markers for NSCLC. Cumulative expression of two or more MMPs/TIMPs may also have increased prognostic significance. Proteases and their inhibitors are novel targets for therapeutic intervention and should be evaluated in NSCLC.

Transient expression vectors for functional genomics, quantification of promoter activity and RNA silencing in plants

Background We describe novel plasmid vectors for transient gene expression using Agrobacterium, infiltrated into Nicotiana benthamiana leaves. We have generated a series of pGreenII cloning vectors that are ideally suited to transient gene expression, by removing elements of conventional binary vectors necessary for stable transformation such as transformation selection genes. Results We give an example of expression of heme-thiolate P450 to demonstrate effectiveness of this system. We have also designed vectors that take advantage of a dual luciferase assay system to analyse promoter sequences or post-transcriptional regulation of gene expression. We have demonstrated their utility by co-expression of putative transcription factors and the promoter sequence of potential target genes and show how orthologous promoter sequences respond to these genes. Finally, we have constructed a vector that has allowed us to investigate design features of hairpin constructs related to their ability to initiate RNA silencing, and have used these tools to study cis-regulatory effect of intron-containing gene constructs. Conclusion In developing a series of vectors ideally suited to transient expression analysis we have provided a resource that further advances the application of this technology. These minimal vectors are ideally suited to conventional cloning methods and we have used them to demonstrate their flexibility to investigate enzyme activity, transcription regulation and post-transcriptional regulatory processes in transient assays.

Frizzled-7 receptor ectodomain expression in a colon cancer cell line induces morphological change and attenuates tumor growth

Frizzled (FZD) receptors have a conserved N-terminal extracellular cysteine-rich domain that interacts with Wnts and co-expression of the receptor ectodomain can antagonize FZD-mediated signalling. Using the ectodomain as an antagonist we have modulated endogenous FZD7 signalling in the moderately differentiated colon adenocarcinoma cell line, SK-CO-1. Unlike the parental cell line, which grows as tightly associated adherent cell clusters, the FZD7 ectodomain expressing cells display a spread out morphology and grow as a monolayer in tissue culture. This transition in morphology was associated with decreased levels of plasma membrane-associated E-cadherin and β-catenin, localized increased levels of vimentin and redistribution of α6 integrin to cellular processes in the FZD7 ectodomain expressing cells. The morphological and phenotype changes induced by FZD7 ectodomain expression in SK-CO-1 cells is thus consistent with the cells undergoing an epithelial-to-mesenchymal-like transition. Furthermore, initiation of tumor formation in a xenograft tumor growth assay was attenuated in the FZD7 ectodomain expressing cells. Our results indicate a pivotal role for endogenous FZD7 in morphology transitions that are associated with colon tumor initiation and progression.

Human mesenchymal stem cells retain multilineage differentiation capacity including neural marker expression after extended in vitro expansion

The suitability of human mesenchymal stem cells (hMSCs) in regenerative medicine relies on retention of their proliferative expansion potential in conjunction with the ability to differentiate toward multiple lineages. Successful utilisation of these cells in clinical applications linked to tissue regeneration requires consideration of biomarker expression, time in culture and donor age, as well as their ability to differentiate towards mesenchymal (bone, cartilage, fat) or non-mesenchymal (e.g., neural) lineages. To identify potential therapeutic suitability we examined hMSCs after extended expansion including morphological changes, potency (stemness) and multilineage potential. Commercially available hMSC populations were expanded in vitro for > 20 passages, equating to > 60 days and > 50 population doublings. Distinct growth phases (A-C) were observed during serial passaging and cells were characterised for stemness and lineage markers at representative stages (Phase A: P+5, approximately 13 days in culture; Phase B: P+7, approximately 20 days in culture; and Phase C: P+13, approximately 43 days in culture). Cell surface markers, stem cell markers and lineage-specific markers were characterised by FACS, ICC and Q-PCR revealing MSCs maintained their multilineage potential, including neural lineages throughout expansion. Co-expression of multiple lineage markers along with continued CD45 expression in MSCs did not affect completion of osteogenic and adipogenic specification or the formation of neurospheres. Improved standardised isolation and characterisation of MSCs may facilitate the identification of biomarkers to improve therapeutic efficacy to ensure increased reproducibility and routine production of MSCs for therapeutic applications including neural repair.

Distress call-induced gene expression in the brain of the Indian short-nosed fruit bat, Cynopterus sphinx

Individuals in distress emit audible vocalizations to either warn or inform conspecifics. The Indian short-nosed fruit bat, Cynopterus sphinx, emits distress calls soon after becoming entangled in mist nets, which appear to attract conspecifics. Phase I of these distress calls is longer and louder, and includes a secondary peak, compared to phase II. Activity-dependent expression of egr-1 was examined in free-ranging C. sphinx following the emissions and responses to a distress call. We found that the level of expression of egr-1 was higher in bats that emitted a distress call, in adults that responded, and in pups than in silent bats. Up-regulated cDNA was amplified to identify the target gene (TOE1) of the protein Egr-1. The observed expression pattern Toe1 was similar to that of egr-1. These findings suggest that the neuronal activity related to recognition of a distress call and an auditory feedback mechanism induces the expression of Egr-1. Co-expression of egr-1 with Toe1 may play a role in initial triggering of the genetic mechanism that could be involved in the consolidation or stabilization of distress call memories.

Over-expression and purification strategies for recombinant multi-protein oligomers: A case study of Mycobacterium tuberculosis sigma/anti-sigma factor protein complexes

The function of a protein in a cell often involves coordinated interactions with one or several regulatory partners. It is thus imperative to characterize a protein both in isolation as well as in the context of its complex with an interacting partner. High resolution structural information determined by X-ray crystallography and Nuclear Magnetic Resonance offer the best route to characterize protein complexes. These techniques, however, require highly purified and homogenous protein samples at high concentration. This requirement often presents a major hurdle for structural studies. Here we present a strategy based on co-expression and co-purification to obtain recombinant multi-protein complexes in the quantity and concentration range that can enable hitherto intractable structural projects. The feasibility of this strategy was examined using the sigma factor/anti-sigma factor protein complexes from Mycobacterium tuberculosis. The approach was successful across a wide range of sigma factors and their cognate interacting partners. It thus appears likely that the analysis of these complexes based on variations in expression constructs and procedures for the purification and characterization of these recombinant protein samples would be widely applicable for other multi-protein systems. (C) 2010 Elsevier Inc. All rights reserved.

Endogenous expression and localization of myostatin and its relation to MHC distribution in C2C12 skeletal muscle cells

Myostatin is a negative regulator of skeletal muscle growth. We have previously reported that recombinant myostatin protein inhibits DNA and protein synthesis in C2C12 cells. Our objective was to assess if C2C12 cells express myostatin, determine its sub-cellular localization and the developmental stage of C2C12 cells in which myostatin mRNA and protein are expressed. To study the endogenous expression of myostatin, C2C12 myoblasts were allowed to progress to myotubes, and changes in the levels of endogenous myostatin mRNA expression were determined by RT-PCR. The myostatin protein and the two major myosin heavy chain (MHC) isoforms (MHC-I and -II) were determined by Western blot. Confirmation of the relative MHC expression patterns was obtained by a modified polyacrylamide gel electropheretic (PAGE) procedure. Imunofluorescence staining was employed to localize the site of myostatin expression and the relative distribution of the MHC isoforms. Co-expression of these proteins was studied using a dual staining approach. Expression of myostatin mRNA was found in myotubes but not in myoblasts. Myostatin protein was seen in most but not all, of the nuclei of polynucleated fibers expressing MHC-II, and myostatin was detected in the cytoplasm of myotube. The localization of myostatin protein in myotube nuclei was confirmed by Western blot of isolated nuclear and cytoplasmic fractions. Incubation of C2C12 myotubes with graded doses of dexamethasone dose-dependently increased the intensity of nuclear myostatin immunostaining and also resulted in the appearance of cytoplasmic expression. In conclusion, myostatin was expressed mostly in C2C12 myotubes nuclei expressing MHC-II. Its predominant

Cyclooxygenase-2 expression correlates with phaeochromocytoma malignancy.

Abstract Aims: Phaeochromocytomas are rare but potentially life-threatening neuroendocrine tumours of the adrenal medulla or sympathetic nervous system ganglia. There are no histological features which reliably differentiate benign from malignant phaeochromocytomas. The current study evaluated cyclooxygenase-2 (Cox-2) and Bcl-2 as tissue-based biomarkers of phaeochromocytoma prognosis. Methods and Results: Cox-2 and Bcl-2 expression were examined immunohistochemically in tissue from forty-one sporadic phaeochromocytoma patients followed up for a minimum of five years after diagnosis. There was a statistically significant association between Cox-2 histoscore (intensity x porportion) and the development of tumour recurrence or metastases (p=0.006). A significant relationship between the co-expression of Cox-2 and Bcl-2 in the primary tumour and the presence of recurrent disease was observed (p=0.034). A highly significant association was observed between, (i) the tumour-associated expression of these two oncoproteins (p=0.001) and, (ii) Cox-2 histoscore and the presence of Bcl-2 expression (p=0.002). Cox regression analysis demonstrated no significant relationship between, (i) the presence or absence of either Cox-2 or Bcl-2 and patient survival or, (ii) between Cox-2 histoscore and patient survival. Conclusions: These results suggest that Cox-2 and Bcl-2 may promote phaeochromocytoma malignancy and that these oncoproteins may be valuable surrogate markers of an aggressive tumour phenotype.

Co-regulation of Gremlin and Notch signalling in diabetic nephropathy

Diabetic nephropathy is currently the leading cause of end-stage renal disease worldwide, and occurs in approximately one third of all diabetic patients. The molecular pathogenesis of diabetic nephropathy has not been fully characterized and novel mediators and drivers of the disease are still being described. Previous data from our laboratory has identified the developmentally regulated gene Gremlin as a novel target implicated in diabetic nephropathy in vitro and in vivo. We used bioinformatic analysis to examine whether Gremlin gene sequence and structure could be used to identify other genes implicated in diabetic nephropathy. The Notch ligand Jagged1 and its downstream effector, hairy enhancer of split-1 (Hes1), were identified as genes with significant similarity to Gremlin in terms of promoter structure and predicted microRNA binding elements. This led us to discover that transforming growth factor-beta (TGFß1), a primary driver of cellular changes in the kidney during nephropathy, increased Gremlin, Jagged1 and Hes1 expression in human kidney epithelial cells. Elevated levels of Gremlin, Jagged1 and Hes1 were also detected in extracts from renal biopsies from diabetic nephropathy patients, but not in control living donors. In situ hybridization identified specific upregulation and co-expression of Gremlin, Jagged1 and Hes1 in the same tubuli of kidneys from diabetic nephropathy patients, but not controls. Finally, Notch pathway gene clustering showed that samples from diabetic nephropathy patients grouped together, distinct from both control living donors and patients with minimal change disease. Together, these data suggest that Notch pathway gene expression is elevated in diabetic nephropathy, co-incident with Gremlin, and may contribute to the pathogenesis of this disease.