Transcriptome profiling of human thymic CD4+ and CD8+ T cells compared to primary peripheral T cells.

BACKGROUND: The thymus is a highly specialized organ of the immune system where T cell precursors develop and differentiate into self-tolerant CD4+ or CD8+ T cells. No studies to date have investigated how the human transcriptome profiles differ, between T cells still residing in the thymus and T cells in the periphery. RESULTS: We have performed high-throughput RNA sequencing to characterize the transcriptomes of primary single positive (SP) CD4+ and CD8+ T cells from infant thymic tissue, as well as primary CD4+ and CD8+ T cells from infant and adult peripheral blood, to enable the comparisons across tissues and ages. In addition, we have assessed the expression of candidate genes related to autoimmune diseases in thymic CD4+ and CD8+ T cells. The thymic T cells showed the largest number of uniquely expressed genes, suggesting a more diverse transcription in thymic T cells. Comparing T cells of thymic and blood origin, revealed more differentially expressed genes, than between infant and adult blood. Functional enrichment analysis revealed an over-representation of genes involved in cell cycle and replication in thymic T cells, whereas infant blood T cells were dominated by immune related terms. Comparing adult and infant blood T cells, the former was enriched for inflammatory response, cytokine production and biological adhesion, while upregulated genes in infant blood T cells were associated with cell cycle, cell death and gene expression. CONCLUSION: This study provides valuable insight into the transcriptomes of the human primary SP T cells still residing within the thymus, and offers a unique comparison to primary blood derived T cells. Interestingly, the majority of autoimmune disease associated genes were expressed in one or more T cell subset, however ~ 11% of these were not expressed in frequently studied adult peripheral blood.

Genetic risk variants for autoimmune diseases that influence gene expression in thymus.

Genome-wide association studies (GWAS) have boosted our knowledge of genetic risk variants in autoimmune diseases (AIDs). Most risk variants are located within or near genes with immunological functions, and the majority is found to be non-coding, pointing towards a regulatory role. In this study, we performed a cis expression quantitative trait locus (eQTL) screen restricted to 353 AID associated risk variants selected from the GWAS catalog to investigate whether these single nucleotide polymorphisms (SNPs) influence gene expression in thymus. Genotypes were obtained by Immunochip (Ichip) and tested against expression of surrounding genes (±1 Mb) in human thymic tissue (n = 42). We identified eight significant eQTLs located within seven genetic regions (FCRL3, RNASET2, C2orf74, NPIPB8, SIRPG, SYS1 and AJ006998.2) where the expression was associated with AID risk SNPs at a study-wide level of significance (P < 2.7 × 10-5). In NPIPB8 and AJ006998.2, the eQTL signals appeared to be thymus-specific. Furthermore, many AID risk SNPs from GWAS have been subsequently fine-mapped in recent Ichip projects, and fine-mapped AID SNPs overlapped with the thymic eQTLs within RNASET2 and SIRPG Finally, in all the eQTL regions, except C2orf74, SNPs underlying the thymic eQTLs were predicted to interfere with transcription factors important in T cell development. Our study therefore reveals autoimmune risk variants that act as eQTLs in thymus, and suggest that thymic gene regulation may play a functional role at some AID risk loci.

Prostate cancer radiogenomics reveals proliferative gene expression programs associated with distinct MRI-based hypoxia levels.

BACKGROUND AND PURPOSE: The biology behind individual hypoxia levels in patient tumors is poorly understood. Here, we used radiogenomics to identify associations between magnetic resonance imaging (MRI)-based hypoxia levels and biological processes derived from gene expression data in prostate cancer. MATERIALS AND METHODS: For 85 prostate cancer patients, MRI-based hypoxia images were constructed by combining diffusion-weighted images reflecting oxygen consumption and supply. The ability to differentiate hypoxia levels in these images was verified by comparison with matched biopsy sections stained for the hypoxia marker pimonidazole. For MRI-defined hypoxia levels, corresponding hypoxic fractions were calculated and correlated with biopsy gene expression profiles. Biological processes were predicted by gene set enrichment analysis (GSEA) and validated by immunohistochemistry (Ki67 proliferation marker, reactive stroma grade) and RT-PCR (MYC). RESULTS: Genes with correlation between expression level and hypoxic fraction were identified for 56 MRI-based hypoxia levels. At all levels, GSEA identified proliferation as the predominant biological process enriched among the correlating genes. Two independent proliferative gene signatures were developed. The Peak1 signature, upregulated at moderate/severe hypoxia, reflected MYC upregulation and high Ki67-proliferation index of cancer cells in pimonidazole-positive regions. The Peak2 signature, upregulated at mild to non-hypoxic levels, was associated with fibroblast gene signature and reactive stroma grade. High scores of both Peak1 and Peak2 indicated elevated risk of biochemical recurrence in multiple cohorts. CONCLUSION: Radiogenomics identified two gene expression programs activated at different hypoxia levels, reflecting proliferation of cancer cells and stroma cells. Genes involved in these programs could be candidate targets for intervention.

Characterisation of a novel paralog of scavenger receptor class B member I (SCARB1) in Atlantic salmon (Salmo salar).

BACKGROUND: Red flesh colour is a unique trait found in some salmonid genera. Carotenoid pigments are not synthesized de novo in the fish, but are provided by dietary uptake. A better understanding of the molecular mechanisms underlying the cellular uptake and deposition of carotenoids could potentially be used to improve the low muscle deposition rate that is typically found in farmed Atlantic salmon. In addition, from an evolutionary point of view, the establishment and maintenance of this trait is still poorly understood. It has been demonstrated in several species that scavenger receptor class B, member 1 (SCARB1) is involved in intestinal absorption of carotenoids, which makes this gene a possible source of genetic variation in salmonid flesh pigmentation. RESULTS: In this study, a novel paralog of SCARB1 (SCARB1-2) was detected through screening for genetic variation in Atlantic salmon SCARB1. Full length SCARB1-2 encodes a protein with 89% identity to Atlantic salmon SCARB1, except for the C-terminal cytoplasmic tail that shows only 12% identity. The most prominent site of SCARB1 mRNA expression was in the mid gut, while a five-fold lower level was detected in Atlantic salmon skeletal muscle and liver. The SCARB1-2 mRNA was equally expressed in liver, muscle and mid gut, and at a lower level than SCARB1 mRNA. A total of seven different SCARB1-2 alleles comprising repetitive enhancer of zeste motifs (EZH2) were identified in the founding parents of a resource Atlantic salmon population. We mapped the SCARB1-2 paralog to a region on Atlantic salmon chromosome 1, containing a putative QTL for flesh colour. Addition of the SCARB1-2 marker increased the significance of this QTL, however the large confidence interval surrounding the QTL precludes confirmation of SCARB1-2 as a causative gene underlying variation in this trait. CONCLUSION: We have characterised a novel paralog of SCARB1 (SCARB1-2), have mapped it to Atlantic salmon chromosome 1 and have described its expression in various tissues. Mapping with SCARB1-2 alleles added further evidence for a QTL affecting flesh colour on this chromosome, however further studies are needed to confirm a functional role for this gene in flesh colour pigmentation.

Genomic and functional gene studies suggest a key role of beta-carotene oxygenase 1 like (bco1l) gene in salmon flesh color.

Red coloration of muscle tissue (flesh) is a unique trait in several salmonid genera, including Atlantic salmon. The color results from dietary carotenoids deposited in the flesh, whereas the color intensity is affected both by diet and genetic components. Herein we report on a genome-wide association study (GWAS) to identify genetic variation underlying this trait. Two SNPs on ssa26 showed strong associations to the flesh color in salmon. Two genes known to be involved in carotenoid metabolism were located in this QTL- region: beta-carotene oxygenase 1 (bco1) and beta-carotene oxygenase 1 like (bco1l). To determine whether flesh color variation is caused by one, or both, of these genes, functional studies were carried out including mRNA and protein expression in fish with red and pale flesh color. The catalytic abilities of these two genes were also tested with different carotenoids. Our results suggest bco1l to be the most likely gene to explain the flesh color variation observed in this population.

Transcriptomes of antigen presenting cells in human thymus.

Antigen presenting cells (APCs) in the thymus play an essential role in the establishment of central tolerance, i.e. the generation of a repertoire of functional and self-tolerant T cells to prevent autoimmunity. In this study, we have compared the transcriptomes of four primary APCs from human thymus (mTECs, CD19+ B cells, CD141+ and CD123+ DCs). We investigated a set of genes including the HLA genes, genes encoding transcriptional regulators and finally, tissue-enriched genes, i.e, genes with a five-fold higher expression in a particular human tissue. We show that thymic CD141+ DCs express the highest levels of all classical HLA genes and 67% (14/21) of the HLA class I and II pathway genes investigated in this study. CD141+ DCs also expressed the highest levels of the transcriptional regulator DEAF1, whereas AIRE and FEZF2 expression were mainly found in primary human mTECs. We found expression of "tissue enriched genes" from the Human Protein Atlas (HPA) in all four APC types, but the mTECs were clearly dominating in the number of uniquely expressed tissue enriched genes (20% in mTECs, 7% in CD19+ B cells, 4% in CD123+ DCs and 2% in CD141+ DCs). The tissue enriched genes also overlapped with reported human autoantigens. This is, to our knowledge, the first study that performs RNA sequencing of mTECs, CD19+ B cells, CD141+ and CD123+ DCs isolated from the same individuals and provides insight into the transcriptomes of these human thymic APCs.

The evolution and functional divergence of the beta-carotene oxygenase gene family in teleost fish--exemplified by Atlantic salmon.

In mammals, two carotenoid cleaving oxygenases are known; beta-carotene 15,15'-monooxygenase (BCMO1) and beta-carotene 9',10'-oxygenase (BCO2). BCMO1 is a key enzyme in vitamin A synthesis by symmetrically cleaving beta-carotene into 2 molecules of all-trans-retinal, while BCO2 is responsible for asymmetric cleavage of a broader range of carotenoids. Here, we show that the Atlantic salmon beta-carotene oxygenase (bco) gene family contains 5 members, three bco2 and two bcmo1 paralogs. Using public sequence databases, multiple bco genes were also found in several additional teleost species. Phylogenetic analysis indicates that bco2a and bco2b originate from the teleost fish specific genome duplication (FSGD or 3R), while the third and more distant paralog, bco2 like, might stem from a prior duplication event in the teleost lineage. The two bcmo1 paralogs (bcmo1 and bcmo1 like) appear to be the result of an ancient duplication event that took place before the divergence of ray-finned (Actinopterygii) and lobe-finned fish (Sarcopterygii), with subsequent nonfunctionalization and loss of one Sarcopterygii paralog. Gene expression analysis of the bcmo1 and bco2 paralogs in Atlantic salmon reveals regulatory divergence with tissue specific expression profiles, suggesting that the beta-carotene oxygenase subtypes have evolved functional divergences. We suggest that teleost fish have evolved and maintained an extended repertoire of beta-carotene oxygenases compared to the investigated Sarcopterygii species, and hypothesize that the main driver behind this functional divergence is the exposure to a diverse set of carotenoids in the aquatic environment.