Genetic modulation during thymocyte development an integrative approach based on microarray meta-analysis

The thymus is a primary lymphoid organ where T cell development takes place. The alpha-beta T lymphocyte lineage has well defined developmental stages, wherein the T cell receptor (TCR) structure and interaction ability are tested so that to guarantee central tolerance in the immune system. We conducted a meta-analysis using microarray data of T cell precursors (thymocytes) but only experiments where the thymocytes were isolated based on the developmental stages defined by the expression of CD4 and CD8 markers. Our approach allowed us to isolate meaningfully different gene categories regardless of the microarray platform, and to provide a perspective on gene modulation during this process. After filtering the data, the remaining 991 genes were sorted according to their expression profile across 27 comparisons, resulting in 15 clusters. We coupled the cluster expression profiles with the evidence in the literature, providing a comprehensive view of thymocyte development. As expected, our results associate the gene expression with TCR signaling pathways. Additionally, they revealed new genes associated with this process, and provided more evidence to understand the role of genes already known to be important to T cell development. Lastly, the clusters were related with other surface receptor signaling pathways, as well as chemotaxis and cell migration process, highlighting the importance of external stimuli to thymocyte differentiation. In conclusion, we have identified, through data integration, pathways and molecules that should have relevant impact on T cell development. We believe that the meta-analysis reported herein improved the comprehension on how the gene expression is organized temporally, providing a global perspective on the thymocyte gene modulation waves.

[1]  M. Ciofani,et al.  The thymus as an inductive site for T lymphopoiesis. , 2007, Annual review of cell and developmental biology.

[2]  L. Klein,et al.  Positive and negative selection of the T cell repertoire: what thymocytes see (and don't see) , 2014, Nature Reviews Immunology.

[3]  R. Bosselut,et al.  Decision checkpoints in the thymus , 2010, Nature Immunology.

[4]  S. Ghosh,et al.  Differential role of the transcription factor NF-kappaB in selection and survival of CD4+ and CD8+ thymocytes. , 2008, Immunity.

[5]  Susan R. Wilson,et al.  Impairment of organ-specific T cell negative selection by diabetes susceptibility genes: genomic analysis by mRNA profiling , 2007, Genome Biology.

[6]  Russ B. Altman,et al.  Missing value estimation methods for DNA microarrays , 2001, Bioinform..

[7]  W. Savino Intrathymic T Cell Migration Is a Multivectorial Process under a Complex Neuroendocrine Control , 2010, Neuroimmunomodulation.

[8]  Renaud Gaujoux,et al.  A flexible R package for nonnegative matrix factorization , 2010, BMC Bioinformatics.

[9]  Corey M. Carlson,et al.  Kruppel-like factor 2 regulates thymocyte and T-cell migration , 2006, Nature.

[10]  C. Benoist,et al.  Thymic negative selection is functional in NOD mice , 2012, The Journal of experimental medicine.

[11]  D. Littman,et al.  Transcription factor AP4 modulates reversible and epigenetic silencing of the Cd4 gene , 2011, Proceedings of the National Academy of Sciences.

[12]  Gordon K. Smyth,et al.  limma: Linear Models for Microarray Data , 2005 .

[13]  A. Kretz-Rommel,et al.  Linkage of immune self-tolerance with the positive selection of T cells. , 1999, Critical reviews in immunology.

[14]  A. Singer,et al.  Coreceptor reversal in the thymus: signaled CD4+8+ thymocytes initially terminate CD8 transcription even when differentiating into CD8+ T cells. , 2000, Immunity.

[15]  Douglas G Altman,et al.  Key Issues in Conducting a Meta-Analysis of Gene Expression Microarray Datasets , 2008, PLoS medicine.

[16]  R. Treisman,et al.  Ternary Complex Factors SAP-1 and Elk-1, but Not Net, Are Functionally Equivalent in Thymocyte Development , 2010, The Journal of Immunology.

[17]  Marie Malissen,et al.  The chemokine TECK is expressed by thymic and intestinal epithelial cells and attracts double‐ and single‐positive thymocytes expressing the TECK receptor CCR9 , 2000, European journal of immunology.

[18]  Ann Nowé,et al.  Comparison of Merging and Meta-Analysis as Alternative Approaches for Integrative Gene Expression Analysis , 2014, ISRN bioinformatics.

[19]  J. Epstein,et al.  PlexinD1 glycoprotein controls migration of positively selected thymocytes into the medulla. , 2008, Immunity.

[20]  Sean R. Davis,et al.  GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor , 2007, Bioinform..

[21]  Pablo Tamayo,et al.  Metagenes and molecular pattern discovery using matrix factorization , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[22]  Grace X. Y. Zheng,et al.  Dynamic regulation of miRNA expression in ordered stages of cellular development. , 2007, Genes & development.

[23]  Alex E. Lash,et al.  Gene Expression Omnibus: NCBI gene expression and hybridization array data repository , 2002, Nucleic Acids Res..

[24]  D. Littman,et al.  Regulation of the TCRalpha repertoire by the survival window of CD4(+)CD8(+) thymocytes. , 2002, Nature immunology.

[25]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[26]  C. Riccardi,et al.  Modulation of Pro- and Antiapoptotic Molecules in Double-Positive (CD4+CD8+) Thymocytes following Dexamethasone Treatment , 2006, Journal of Pharmacology and Experimental Therapeutics.

[27]  Dorothy N. Kakoola,et al.  Molecular pathway alterations in CD4 T-cells of nonobese diabetic (NOD) mice in the preinsulitis phase of autoimmune diabetes , 2014, Results in immunology.

[28]  Dorothy N. Kakoola,et al.  Molecular Phenotyping of Immune Cells from Young NOD Mice Reveals Abnormal Metabolic Pathways in the Early Induction Phase of Autoimmune Diabetes , 2012, PloS one.

[29]  Jill P. Mesirov,et al.  Consensus Clustering: A Resampling-Based Method for Class Discovery and Visualization of Gene Expression Microarray Data , 2003, Machine Learning.

[30]  P. Herzmark,et al.  Elimination of Self-Reactive T Cells in the Thymus: A Timeline for Negative Selection , 2013, PLoS biology.

[31]  L. Aravind,et al.  Corrigendum: Themis is a member of a new metazoan gene family and is required for the completion of thymocyte positive selection , 2010, Nature Immunology.

[32]  Gary D Bader,et al.  GeneMANIA: Fast gene network construction and function prediction for Cytoscape , 2014, F1000Research.

[33]  W. Savino,et al.  Thymocyte migration: an affair of multiple cellular interactions? , 2003, Brazilian journal of medical and biological research = Revista brasileira de pesquisas medicas e biologicas.

[34]  Jin Lu,et al.  Maturation and Emigration of Single-Positive Thymocytes , 2013, Clinical & developmental immunology.

[35]  V. Engelhard,et al.  Regulation of T-cell Tolerance by Lymphatic Endothelial Cells , 2014, Journal of clinical & cellular immunology.

[36]  D. Oh,et al.  Calcium oscillations regulate thymocyte motility during positive selection in the three-dimensional thymic environment , 2005, Nature Immunology.

[37]  Sean R. Davis,et al.  NCBI GEO: archive for functional genomics data sets—update , 2012, Nucleic Acids Res..

[38]  D. Mottet,et al.  Histone Deacetylase 7 Regulates Cell Survival and TCR Signaling in CD4/CD8 Double-Positive Thymocytes , 2011, The Journal of Immunology.