Gene expression profiling in the human pathogenic dermatophyte Trichophyton rubrum during growth on proteins.

Dermatophytes are highly specialized filamentous fungi which cause the majority of superficial mycoses in humans and animals. The high secreted proteolytic activity of these microorganisms during growth on proteins is assumed to be linked to their particular ability to exclusively infect keratinized host structures such as the skin stratum corneum, hair, and nails. Individual secreted dermatophyte proteases were recently described and linked with the in vitro digestion of keratin. However, the overall adaptation and transcriptional response of dermatophytes during protein degradation are largely unknown. To address this question, we constructed a cDNA microarray for the human pathogenic dermatophyte Trichophyton rubrum that was based on transcripts of the fungus grown on proteins. Profiles of gene expression during the growth of T. rubrum on soy and keratin protein displayed the activation of a large set of genes that encode secreted endo- and exoproteases. In addition, other specifically induced factors potentially implicated in protein utilization were identified, including heat shock proteins, transporters, metabolic enzymes, transcription factors, and hypothetical proteins with unknown functions. Of particular interest is the strong upregulation of key enzymes of the glyoxylate cycle in T. rubrum during growth on soy and keratin, namely, isocitrate lyase and malate synthase. This broad-scale transcriptional analysis of dermatophytes during growth on proteins reveals new putative pathogenicity-related host adaptation mechanisms of these human pathogenic fungi.

Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Test of four colon cancer risk-scores in formalin fixed paraffin embedded microarray gene expression data.

BACKGROUND: Prognosis prediction for resected primary colon cancer is based on the T-stage Node Metastasis (TNM) staging system. We investigated if four well-documented gene expression risk scores can improve patient stratification. METHODS: Microarray-based versions of risk-scores were applied to a large independent cohort of 688 stage II/III tumors from the PETACC-3 trial. Prognostic value for relapse-free survival (RFS), survival after relapse (SAR), and overall survival (OS) was assessed by regression analysis. To assess improvement over a reference, prognostic model was assessed with the area under curve (AUC) of receiver operating characteristic (ROC) curves. All statistical tests were two-sided, except the AUC increase. RESULTS: All four risk scores (RSs) showed a statistically significant association (single-test, P < .0167) with OS or RFS in univariate models, but with HRs below 1.38 per interquartile range. Three scores were predictors of shorter RFS, one of shorter SAR. Each RS could only marginally improve an RFS or OS model with the known factors T-stage, N-stage, and microsatellite instability (MSI) status (AUC gains < 0.025 units). The pairwise interscore discordance was never high (maximal Spearman correlation = 0.563) A combined score showed a trend to higher prognostic value and higher AUC increase for OS (HR = 1.74, 95% confidence interval [CI] = 1.44 to 2.10, P < .001, AUC from 0.6918 to 0.7321) and RFS (HR = 1.56, 95% CI = 1.33 to 1.84, P < .001, AUC from 0.6723 to 0.6945) than any single score. CONCLUSIONS: The four tested gene expression-based risk scores provide prognostic information but contribute only marginally to improving models based on established risk factors. A combination of the risk scores might provide more robust information. Predictors of RFS and SAR might need to be different.

Hierarchal utilization of different T-cell receptor Vbeta gene segments in the CD8(+)-T-cell response to an immunodominant Moloney leukemia virus-encoded epitope in vivo.

The CD8(+)-T-cell response to Moloney murine leukemia virus (M-MuLV)-associated antigens in C57BL/6 mice is directed against an immunodominant gag-encoded epitope (CCLCLTVFL) presented in the context of H-2D(b) and is restricted primarily to cytotoxic T lymphocytes (CTL) expressing the Valpha3.2 and Vbeta5.2 gene segments. We decided to examine the M-MuLV response in congenic C57BL/6 Vbeta(a) mice which are unable to express the dominant Valpha3.2(+) Vbeta5.2(+) T-cell receptor (TCR) due to a large deletion at the TCR locus that includes the Vbeta5.2 gene segment. Interestingly, M-MuLV-immune C57BL/6 Vbeta(a) mice were still able to reject M-MuLV-infected tumor cells and direct ex vivo analysis of peripheral blood lymphocytes from these immune mice revealed a dramatic increase in CD8(+) cells utilizing the same Valpha3.2 gene segment in association with two different Vbeta segments (Vbeta3 and Vbeta17). Surprisingly, all these CTL recognized the same immunodominant M-MuLV gag epitope. Analysis of the TCR repertoire of individual M-MuLV-immune (C57BL/6 x C57BL/6 Vbeta(a))F(1) mice revealed a clear hierarchy in Vbeta utilization, with a preferential usage of the Vbeta17 gene segment, whereas Vbeta3 and especially Vbeta5.2 were used to much lesser extents. Sequencing of TCRalpha- and -beta-chain junctional regions of CTL clones specific for the M-MuLV gag epitope revealed a diverse repertoire of TCRbeta chains in Vbeta(a) mice and a highly restricted TCRbeta-chain repertoire in Vbeta(b) mice, whereas TCRalpha-chain sequences were highly conserved in both cases. Collectively, our data indicate that the H-2D(b)-restricted M-MuLV gag epitope can be recognized in a hierarchal fashion by different Vbeta domains and that the degree of beta-chain diversity varies according to Vbeta utilization.

The evolution of gene expression levels in mammalian organs.

Changes in gene expression are thought to underlie many of the phenotypic differences between species. However, large-scale analyses of gene expression evolution were until recently prevented by technological limitations. Here we report the sequencing of polyadenylated RNA from six organs across ten species that represent all major mammalian lineages (placentals, marsupials and monotremes) and birds (the evolutionary outgroup), with the goal of understanding the dynamics of mammalian transcriptome evolution. We show that the rate of gene expression evolution varies among organs, lineages and chromosomes, owing to differences in selective pressures: transcriptome change was slow in nervous tissues and rapid in testes, slower in rodents than in apes and monotremes, and rapid for the X chromosome right after its formation. Although gene expression evolution in mammals was strongly shaped by purifying selection, we identify numerous potentially selectively driven expression switches, which occurred at different rates across lineages and tissues and which probably contributed to the specific organ biology of various mammals.

In vivo fate mapping with SCL regulatory elements identifies progenitors for primitive and definitive hematopoiesis in mice.

One of the principal issues facing biomedical research is to elucidate developmental pathways and to establish the fate of stem and progenitor cells in vivo. Hematopoiesis, the process of blood cell formation, provides a powerful experimental system for investigating this process. Here, we employ transcriptional regulatory elements from the stem cell leukemia (SCL) gene to selectively label primitive and definitive hematopoiesis. We report that SCL-labelled cells arising in the mid to late streak embryo give rise to primitive red blood cells but fail to contribute to the vascular system of the developing embryo. Restricting SCL-marking to different stages of foetal development, we identify a second population of multilineage progenitors, proficient in contributing to adult erythroid, myeloid and lymphoid cells. The distinct lineage-restricted potential of SCL-labelled early progenitors demonstrates that primitive erythroid cell fate specification is initiated during mid gastrulation. Our data also suggest that the transition from a hemangioblastic precursors with endothelial and blood forming potential to a committed hematopoietic progenitor must have occurred prior to SCL-marking of definitive multilineage blood precursors.

Dramatic influence of V beta gene polymorphism on an antigen-specific CD8+ T cell response in vivo.

According to recent crystallographic studies, the TCR-alpha beta contacts MHC class I-bound antigenic peptides via the polymorphic V gene-encoded complementarity-determining region 1 beta (CDR1 beta) and the hypervariable (D)J-encoded CDR3 beta and CDR3 alpha domains. To evaluate directly the relative importance of CDR1 beta polymorphism on the fine specificity of T cell responses in vivo, we have taken advantage of congenic V beta a and V beta b mouse strains that differ by a CDR1 polymorphism in the V beta 10 gene segment. The V beta 10-restricted CD8+ T cell response to a defined immunodominant epitope was dramatically reduced in V beta a compared with V beta b mice, as measured either by the expansion of V beta 10+ cells or by the binding of MHC-peptide tetramers. These data indicate that V beta polymorphism has an important impact on TCR-ligand binding in vivo, presumably by modifying the affinity of CDR1 beta-peptide interactions.

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.

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.

Nutritional and insulin regulation of fatty acid synthetase and leptin gene expression through ADD1/SREBP1.

The ability to regulate specific genes of energy metabolism in response to fasting and feeding is an important adaptation allowing survival of intermittent food supplies. However, little is known about transcription factors involved in such responses in higher organisms. We show here that gene expression in adipose tissue for adipocyte determination differentiation dependent factor (ADD) 1/sterol regulatory element binding protein (SREBP) 1, a basic-helix-loop-helix protein that has a dual DNA-binding specificity, is reduced dramatically upon fasting and elevated upon refeeding; this parallels closely the regulation of two adipose cell genes that are crucial in energy homeostasis, fatty acid synthetase (FAS) and leptin. This elevation of ADD1/SREBP1, leptin, and FAS that is induced by feeding in vivo is mimicked by exposure of cultured adipocytes to insulin, the classic hormone of the fed state. We also show that the promoters for both leptin and FAS are transactivated by ADD1/SREBP1. A mutation in the basic domain of ADD1/SREBP1 that allows E-box binding but destroys sterol regulatory element-1 binding prevents leptin gene transactivation but has no effect on the increase in FAS promoter function. Molecular dissection of the FAS promoter shows that most if not all of this action of ADD1/SREBP1 is through an E-box motif at -64 to -59, contained with a sequence identified previously as the major insulin response element of this gene. These results indicate that ADD1/SREBP1 is a key transcription factor linking changes in nutritional status and insulin levels to the expression of certain genes that regulate systemic energy metabolism.

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.

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.

Cardiac raptor ablation impairs adaptive hypertrophy, alters metabolic gene expression, and causes heart failure in mice.

Background- Cardiac hypertrophy involves growth responses to a variety of stimuli triggered by increased workload. It is an independent risk factor for heart failure and sudden death. Mammalian target of rapamycin (mTOR) plays a key role in cellular growth responses by integrating growth factor and energy status signals. It is found in 2 structurally and functionally distinct multiprotein complexes called mTOR complex (mTORC) 1 and mTORC2. The role of each of these branches of mTOR signaling in the adult heart is currently unknown. Methods and Results- We generated mice with deficient myocardial mTORC1 activity by targeted ablation of raptor, which encodes an essential component of mTORC1, during adulthood. At 3 weeks after the deletion, atrial and brain natriuretic peptides and β-myosin heavy chain were strongly induced, multiple genes involved in the regulation of energy metabolism were altered, but cardiac function was normal. Function deteriorated rapidly afterward, resulting in dilated cardiomyopathy and high mortality within 6 weeks. Aortic banding-induced pathological overload resulted in severe dilated cardiomyopathy already at 1 week without a prior phase of adaptive hypertrophy. The mechanism involved a lack of adaptive cardiomyocyte growth via blunted protein synthesis capacity, as supported by reduced phosphorylation of ribosomal S6 kinase 1 and 4E-binding protein 1. In addition, reduced mitochondrial content, a shift in metabolic substrate use, and increased apoptosis and autophagy were observed. Conclusions- Our results demonstrate an essential function for mTORC1 in the heart under physiological and pathological conditions and are relevant for the understanding of disease states in which the insulin/insulin-like growth factor signaling axis is affected such as diabetes mellitus and heart failure or after cancer therapy.

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.

Encapsulation of packaging cell line results in successful retroviral-mediated transfer of a suicide gene in vivo in an experimental model of glioblastoma.

AIMS: Retroviral-mediated gene therapy has been proposed as a primary or adjuvant treatment for advanced cancer, because retroviruses selectively infect dividing cells. Efficacy of retroviral-mediated gene transfer, however, is limited in vivo. Although packaging cell lines can produce viral vectors continuously, such allo- or xenogeneic cells are normally rejected when used in vivo. Encapsulation using microporous membranes can protect the packaging cells from rejection. In this study, we used an encapsulated murine packaging cell line to test the effects of in situ delivery of a retrovirus bearing the herpes simplex virus thymidine kinase suicide gene in a rat model of orthotopic glioblastoma. MATERIALS AND METHODS: To test gene transfer in vitro, encapsulated murine psi2-VIK packaging cells were co-cultured with baby hamster kidney (BHK) cells, and the percentage of transfected BHK cells was determined. For in vivo experiments, orthotopic C6 glioblastomas were established in Wistar rats. Capsules containing psi2-VIK cells were stereotaxically implanted into these tumours and the animals were treated with ganciclovir (GCV). Tumours were harvested 14 days after initiation of GCV therapy for morphometric analysis. RESULTS: Encapsulation of psi2-VIK cells increased transfection rates of BHK target cells significantly in vitro compared to psi2-VIK conditioned medium (3 x 10(6) vs 2.3 x 10(4) cells; P<0.001). In vivo treatment with encapsulated packaging cells resulted in 3% to 5% of C6 tumour cells transduced and 45% of tumour volume replaced by necrosis after GCV (P<0.01 compared to controls). CONCLUSION: In this experimental model of glioblastoma, encapsulation of a xenogeneic packaging cell line increased half-life and transduction efficacy of retrovirus-mediated gene transfer and caused significant tumour necrosis.