Genome-wide association study and trans-ethnic meta-analysis identify novel susceptibility loci for type 2 diabetes mellitus.

BACKGROUND: The genetic basis of type 2 diabetes (T2D) is under-investigated in the Middle East, despite the rapidly growing disease prevalence. We aimed to define the genetic determinants of T2D in Qatar. METHODS: Using whole genome sequencing of 11,436 participants (2765 T2D cases and 8671 controls) from the population-based Qatar Biobank (QBB), we conducted a genome-wide association study (GWAS) of T2D with and without body mass index (BMI) adjustment. RESULTS: We replicated 93 known T2D-associated loci in a BMI-unadjusted model, while 96 known loci were replicated in a BMI-adjusted model. The effect sizes and allele frequencies of replicated SNPs in the Qatari population generally concurred with those from European populations. We identified a locus specific to our cohort located between the APOBEC3H and CBX7 genes in the BMI-unadjusted model. Also, we performed a transethnic meta-analysis of our cohort with a previous GWAS on T2D in multi-ancestry individuals (180,834 T2D cases and 1,159,055 controls). One locus in DYNC2H1 gene reached genome-wide significance in the meta-analysis. Assessing polygenic risk scores derived from European- and multi-ancestries in the Qatari population showed higher predictive performance of the multi-ancestry panel compared to the European panel. CONCLUSION: Our study provides new insights into the genetic architecture of T2D in a Middle Eastern population and identifies genes that may be explored further for their involvement in T2D pathogenesis.

Identification of distinct circulating microRNAs in acute ischemic stroke patients with type 2 diabetes mellitus.

Stroke is the second leading cause of global mortality and continued efforts aim to identify predictive, diagnostic, or prognostic biomarkers to reduce the disease burden. Circulating microRNAs (miRNAs) have emerged as potential biomarkers in stroke. We performed comprehensive circulating miRNA profiling of ischemic stroke patients with or without type 2 diabetes mellitus (T2DM), an important risk factor associated with worse clinical outcomes in stroke. Serum samples were collected within 24 h of acute stroke diagnosis and circulating miRNAs profiled using RNA-Seq were compared between stroke patients with T2DM (SWDM; n = 92) and those without T2DM (SWoDM; n = 98). Our analysis workflow involved random allocation of study cohorts into discovery (n = 96) and validation (n = 94) datasets. Five miRNAs were found to be differentially regulated in SWDM compared to SWoDM patients. Hsa-miR-361-3p and -664a-5p were downregulated, whereas miR-423-3p, -140-5p, and -17-3p were upregulated. We also explored the gene targets of these miRNAs and investigated the downstream pathways associated with them to decipher the potential pathways impacted in stroke with diabetes as comorbidity. Overall, our novel findings provide important insights into the differentially regulated miRNAs, their associated pathways and potential utilization for clinical benefits in ischemic stroke patients with diabetes.

The Prevalence and Genetic Spectrum of Familial Hypercholesterolemia in Qatar Based on Whole Genome Sequencing of 14,000 Subjects.

Familial hypercholesterolemia (FH) is an inherited disease characterized by reduced efficiency of low-density lipoprotein-cholesterol (LDL-C) removal from the blood and, consequently, an increased risk of life-threatening early cardiovascular complications. In Qatar, the prevalence of FH has not been determined and the disease, as in many countries, is largely underdiagnosed. In this study, we combined whole-genome sequencing data from the Qatar Genome Program with deep phenotype data from Qatar Biobank for 14,056 subjects to determine the genetic spectrum and estimate the prevalence of FH in Qatar. We used the Dutch Lipid Clinic Network (DLCN) as a diagnostic tool and scrutinized 11 FH-related genes for known pathogenic and possibly pathogenic mutations. Results revealed an estimated prevalence of 0.8% (1:125) for definite/probable cases of FH in the Qatari population. We detected 16 known pathogenic/likely pathogenic mutations in LDLR and one in PCSK9; all in a heterozygous state with high penetrance. The most common mutation was rs1064793799 (c.313+3A >C) followed by rs771019366 (p.Asp90Gly); both in LDLR. In addition, we identified 18 highly penetrant possibly pathogenic variants, of which 5 were Qatari-specific, in LDLR, APOB, PCSK9 and APOE, which are predicted to be among the top 1% most deleterious mutations in the human genome but further validations are required to confirm their pathogenicity. We did not detect any homozygous FH or autosomal recessive mutations in our study cohort. This pioneering study provides a reliable estimate of FH prevalence in Qatar based on a significantly large population-based cohort, whilst uncovering the spectrum of genetic variants associated with FH.

Circulating and Tumor-Infiltrating Immune Checkpoint-Expressing CD8+ Treg/T Cell Subsets and Their Associations with Disease-Free Survival in Colorectal Cancer Patients.

T cells in the tumor microenvironment (TME) have diverse roles in anti-tumor immunity, including orchestration of immune responses and anti-tumor cytotoxic attack. However, different T cell subsets may have opposing roles in tumor progression, especially in inflammation-related cancers such as colorectal cancer (CRC). In this study, we phenotypically characterized CD3+CD4- (CD8+) T cells in colorectal tumor tissues (TT), normal colon tissues (NT) and in circulation of CRC patients. We investigated the expression levels of key immune checkpoints (ICs) and Treg-related markers in CD8+ T cells. Importantly, we investigated associations between different tumor-infiltrating CD8+ T cell subpopulations and disease-free survival (DFS) in CRC patients. We found that FoxP3 expression and ICs including PD-1, CTLA-4, TIM-3, and LAG-3 were significantly increased in tumor-infiltrating CD8+ T cells compared with NT and peripheral blood. In the TME, we found that TIM-3 expression was significantly increased in patients with early stages and absent lymphovascular invasion (LVI) compared to patients with advanced stages and LVI. Importantly, we report that high levels of certain circulating CD8+ T cell subsets (TIM-3-expressing, FoxP3-Helios-TIM-3+ and FoxP3-Helios+TIM-3+ cells) in CRC patients were associated with better DFS. Moreover, in the TME, we report that elevated levels of CD25+ and TIM-3+ T cells, and FoxP3+Helios-TIM-3+ Tregs were associated with better DFS.

Identification of Novel Circulating miRNAs in Patients with Acute Ischemic Stroke.

Ischemic strokes are associated with significant morbidity and mortality, but currently there are no reliable prognostic or diagnostic blood biomarkers. MicroRNAs (miRNAs) regulate various molecular pathways and may be used as biomarkers. Using RNA-Seq, we conducted comprehensive circulating miRNA profiling in patients with ischemic stroke compared with healthy controls. Samples were collected within 24 h of clinical diagnosis. Stringent analysis criteria of discovery (46 cases and 95 controls) and validation (47 cases and 96 controls) cohorts led to the identification of 10 differentially regulated miRNAs, including 5 novel miRNAs, with potential diagnostic significance. Hsa-miR-451a was the most significantly upregulated miRNA (FC; 4.8, FDR; 3.78 × 10-85), while downregulated miRNAs included hsa-miR-574-5p and hsa-miR-142-3p, among others. Importantly, we computed a multivariate classifier based on the identified miRNA panel to differentiate between ischemic stroke patients and healthy controls, which showed remarkably high sensitivity (0.94) and specificity (0.99). The area under the ROC curve was 0.97 and it is superior to other current available biomarkers. Moreover, in samples collected one month following stroke, we found sustained upregulation of hsa-miR-451a and downregulation of another 5 miRNAs. Lastly, we report 3 miRNAs that were significantly associated with poor clinical outcomes of stroke, as defined by the modified Rankin scores. The clinical translation of the identified miRNA panel may be explored further.

Circulating MicroRNA Profiling Identifies Distinct MicroRNA Signatures in Acute Ischemic Stroke and Transient Ischemic Attack Patients.

Transient ischemic attack (TIA) refers to a momentary neurologic deficit caused by focal cerebral, spinal or retinal ischemic insult. TIA is associated with a high risk of impending acute ischemic stroke (AIS), a neurologic dysfunction characterized by focal cerebral, spinal or retinal infarction. Understanding the differences in molecular pathways in AIS and TIA has merit for deciphering the underlying cause for neuronal deficits with long-term effects and high risks of morbidity and mortality. In this study, we performed comprehensive investigations into the circulating microRNA (miRNA) profiles of AIS (n = 191) and TIA (n = 61) patients. We performed RNA-Seq on serum samples collected within 24 hrs of clinical diagnosis and randomly divided the study populations into discovery and validation cohorts. We identified a panel of 11 differentially regulated miRNAs at FDR < 0.05. Hsa-miR-548c-5p, -20a-5p, -18a-5p, -484, -652-3p, -486-3p, -24-3p, -181a-5p and -222-3p were upregulated, while hsa-miR-500a-3p and -206 were downregulated in AIS patients compared to TIA patients. We also probed the previously validated gene targets of our identified miRNA panel to highlight the molecular pathways affected in AIS. Moreover, we developed a multivariate classifier with potential utilization as a discriminative biomarker for AIS and TIA patients. The underlying molecular pathways in AIS compared to TIA may be explored further in functional studies for therapeutic targeting in clinical translation.

The Genetic Spectrum of Maturity-Onset Diabetes of the Young (MODY) in Qatar, a Population-Based Study.

Maturity-onset diabetes of the young (MODY) is a rare monogenic form of diabetes mellitus. In this study, we estimated the prevalence and genetic spectrum of MODY in the Middle Eastern population of Qatar using whole-genome sequencing (WGS) of 14,364 subjects from the population-based Qatar biobank (QBB) cohort. We focused our investigations on 14 previously identified genes ascribed to the cause of MODY and two potentially novel MODY-causing genes, RFX6 and NKX6-1. Genetic variations within the 16 MODY-related genes were assessed for their pathogenicity to identify disease-causing mutations. Analysis of QBB phenotype data revealed 72 subjects (0.5%) with type 1 diabetes, 2915 subjects (20.3%) with type 2 diabetes and 11,377 (79.2%) without diabetes. We identified 22 mutations in 67 subjects that were previously reported in the Human Genetic Mutation Database (HGMD) as disease-causing (DM) or likely disease causing (DM?) for MODY. We also identified 28 potentially novel MODY-causing mutations, predicted to be among the top 1% most deleterious mutations in the human genome, which showed complete (100%) disease penetrance in 34 subjects. Overall, we estimated that MODY accounts for around 2.2-3.4% of diabetes patients in Qatar. This is the first population-based study to determine the genetic spectrum and estimate the prevalence of MODY in the Middle East. Further research to characterize the newly identified mutations is warranted.

Comprehensive Transcriptomic Profiling of Murine Osteoclast Differentiation Reveals Novel Differentially Expressed Genes and LncRNAs.

Osteoclasts are the sole bone resorbing cells, which undertake opposing roles to osteoblasts to affect skeletal mass and structure. However, unraveling the comprehensive molecular mechanisms behind osteoclast differentiation is necessitated to overcome limitations and scarcity of available data, particularly in relation with the emerging roles of long non-coding RNAs (LncRNAs) in gene expression. In this study, we performed comprehensive and progressive analyses of the dynamic transcriptomes of murine osteoclasts, generated in vitro. We compared the total RNA-based transcriptomes of murine bone marrow derived cells with differentiated osteoclasts, while focusing on potentially novel genes and LncRNAs, to uncover critical genes and their associated pathways, which are differentially regulated during osteoclast differentiation. We found 4,214 differentially regulated genes during osteoclast differentiation, which included various types of LncRNAs. Among the upregulated protein coding genes not previously associated with osteoclast are Pheta1, Hagh, Gfpt1 and Nol4, while downregulated genes included Plau, Ltf, Sell and Zfp831. Notably, we report Nol4 as a novel gene related to osteoclast activity since Nol4 knockout mice Nol4 em1(International Mouse Phenotyping Consortium)J exhibit increased bone mineral density. Moreover, the differentially expressed LncRNAs included antisense and long intergenic non-coding RNAs, among others. Overall, immune-related and metabolism-related genes were downregulated, while anatomical morphogenesis and remodeling-related genes were upregulated in early-differentiated osteoclasts with sustained downregulation of immune-related genes in mature osteoclasts. The gene signatures and the comprehensive transcriptome of osteoclast differentiation provided herein can serve as an invaluable resource for deciphering gene dysregulation in osteoclast-related pathologic conditions.

Integrated whole transcriptome and small RNA analysis revealed multiple regulatory networks in colorectal cancer.

Colorectal cancer (CRC) remains a global disease burden and a leading cause of cancer related deaths worldwide. The identification of aberrantly expressed messenger RNA (mRNA), long non-coding RNA (lncRNA), and microRNA (miRNA), and the resulting molecular interactions and signaling networks is essential for better understanding of CRC, identification of novel diagnostic biomarkers and potential development of therapeutic interventions. Herein, we performed microRNA (miRNA) sequencing on fifteen CRC and their non-tumor adjacent tissues and whole transcriptome RNA-Seq on six paired samples from the same cohort and identified alterations in miRNA, mRNA, and lncRNA expression. Computational analyses using Ingenuity Pathway Analysis (IPA) identified multiple activated signaling networks in CRC, including ERBB2, RABL6, FOXM1, and NFKB networks, while functional annotation highlighted activation of cell proliferation and migration as the hallmark of CRC. IPA in combination with in silico prediction algorithms and experimentally validated databases gave insight into the complex associations and interactions between downregulated miRNAs and upregulated mRNAs in CRC and vice versa. Additionally, potential interaction between differentially expressed lncRNAs such as H19, SNHG5, and GATA2-AS1 with multiple miRNAs has been revealed. Taken together, our data provides thorough analysis of dysregulated protein-coding and non-coding RNAs in CRC highlighting numerous associations and regulatory networks thus providing better understanding of CRC.

Transcriptome of Tumor-Infiltrating T Cells in Colorectal Cancer Patients Uncovered a Unique Gene Signature in CD4+ T Cells Associated with Poor Disease-Specific Survival.

Colorectal cancer (CRC) is influenced by infiltration of immune cell populations in the tumor microenvironment. While elevated levels of cytotoxic T cells are associated with improved prognosis, limited studies have reported associations between CD4+ T cells and disease outcomes. We recently performed transcriptomic profiling and comparative analyses of sorted CD4+ and CD8+ tumor-infiltrating lymphocytes (TILs) from bulk tumors of CRC patients with varying disease stages. In this study, we compared the transcriptomes of CD4+ with CD8+ TILs. Functional annotation pathway analyses revealed the downregulation of inflammatory response-related genes, while T cell activation and angiogenesis-related genes were upregulated in CD4+ TILs. The top 200 deregulated genes in CD4+ TILs were aligned with the cancer genome atlas (TCGA) CRC dataset to identify a unique gene signature associated with poor prognosis. Moreover, 69 upregulated and 20 downregulated genes showed similar trends of up/downregulation in the TCGA dataset and were used to calculate "poor prognosis score" (ppScore), which was significantly associated with disease-specific survival. High ppScore patients showed lower expression of Treg-, Th1-, and Th17-related genes, and higher expression of Th2-related genes. Our data highlight the significance of T cells within the TME and identify a unique candidate prognostic gene signature for CD4+ TILs in CRC patients.

Metabolic reprogramming of T regulatory cells in the hypoxic tumor microenvironment.

Metabolic dysregulation in the hypoxic tumor microenvironment (TME) is considered as a hallmark of solid tumors, leading to changes in biosynthetic pathways favoring onset, survival and proliferation of malignant cells. Within the TME, hypoxic milieu favors metabolic reprogramming of tumor cells, which subsequently affects biological properties of tumor-infiltrating immune cells. T regulatory cells (Tregs), including both circulating and tissue-resident cells, are particularly susceptible to hypoxic metabolic signaling that can reprogram their biological and physicochemical properties. Furthermore, metabolic reprogramming modifies Tregs to utilize alternative substrates and undergo a plethora of metabolic events to meet their energy demands. Major impact of this metabolic reprogramming can result in differentiation, survival, excessive secretion of immunosuppressive cytokines and proliferation of Tregs within the TME, which in turn dampen anti-tumor immune responses. Studies on fine-tuning of Treg metabolism are challenging due to heterogenicity of tissue-resident Tregs and their dynamic functions. In this review, we highlight tumor intrinsic and extrinsic factors, which can influence Treg metabolism in the hypoxic TME. Moreover, we focus on metabolic reprogramming of Tregs that could unveil potential regulatory networks favoring tumorigenesis/progression, and provide novel insights, including inhibitors against acetyl-coA carboxylase 1 and transforming growth factor beta into targeting Treg metabolism for therapeutic benefits.

Tumor-Infiltrating Lymphoid Cells in Colorectal Cancer Patients with Varying Disease Stages and Microsatellite Instability-High/Stable Tumors.

Immune checkpoint inhibition is an effective anti-cancer therapeutic approach but has shown limited efficacy in treating colorectal cancer (CRC) patients. Importantly, immune constituents of the tumor microenvironment (TME) can influence therapy response and cancer progression. We investigated the expression of immune checkpoints (ICs) on lymphoid populations within the CRC TME and compared with cells from normal colon tissues using samples from 50 patients with varying disease stages. We found that the levels of B cells, T cells, and NK cells were similar, IC-expressing CD4+ and CD4+CD8+ double positive T cells were higher, while CD8+ T cells and CD4-CD8- double negative T cells were significantly lower in CRC tumors. Notably, patients with mismatch-repair deficiency/microsatellite instability-high tumors had higher levels of IC-expressing CD4+ and CD8+ T cells than patients with proficient MMR and microsatellite stable tumors. Lastly, The Cancer Genome Atlas Colon Adenocarcinoma datasets showed associations between low expression of selective genes and poorer progression-free interval. Our findings highlight differential expression of ICs on lymphoid cells in CRC tumors in the era of cancer immunotherapy, which at present is solely approved for anti-PD-1 therapy in patients with dMMR/MSI-H tumors. Further investigations into their functionality have potentials for deciphering resistance mechanisms to IC inhibition.

Transcriptomic Profiling of Circulating HLA-DR- Myeloid Cells, Compared with HLA-DR+ Myeloid Antigen-presenting Cells.

Myeloid-derived suppressor cells (MDSCs) are a heterogeneous population of cells with potent immunosuppressive functions, which can inhibit the activation of immune responses under a steady-state condition and pathological conditions. We performed transcriptomic profiling of circulating CD33+HLA-DR+ myeloid antigen-presenting cells (APCs) and CD33+HLA-DR- myeloid cells (potentially MDSCs) in healthy individuals. We sorted both subpopulations from peripheral blood mononuclear cells (PBMCs) of 10 healthy donors and performed RNA sequencing (RNA-Seq). We found that several signaling pathways associated with the positive regulation of immune responses, such as antigen presentation/processing, FcγR-mediated phagocytosis and immune cell trafficking, phosphoinositide 3-kinase (PI3K)/Akt signaling, DC maturation, triggering receptor expressed on myeloid cells 1 (TREM1) signaling, nuclear factor of activated T cells (NFAT) and IL-8 signaling were downregulated in CD33+HLA-DR- myeloid cells. In contrast, pathways implicated in tumor suppression and anti-inflammation, including peroxisome proliferator-activated receptor (PPAR) and phosphatase and tensin homolog (PTEN), were upregulated in CD33+HLA-DR- myeloid cells. These data indicate that PPAR/PTEN axis could be upregulated in myeloid cells to keep the immune system in check in normal physiological conditions. Our data reveal some of the molecular and functional differences between CD33+HLA-DR+ APCs and CD33+HLA-DR- myeloid cells in a steady-state condition, reflecting the potential suppressive function of CD33+HLA-DR- myeloid cells to maintain immune tolerance. For future studies, the same methodological approach could be applied to perform transcriptomic profiling of myeloid subsets in pathological conditions.

T-cell responses and therapies against SARS-CoV-2 infection.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, a novel coronavirus strain. Some studies suggest that COVID-19 could be an immune-related disease, and failure of effective immune responses in initial stages of viral infection could contribute to systemic inflammation and tissue damage, leading to worse disease outcomes. T cells can act as a double-edge sword with both pro- and anti-roles in the progression of COVID-19. Thus, better understanding of their roles in immune responses to SARS-CoV-2 infection is crucial. T cells primarily react to the spike protein on the coronavirus to initiate antiviral immunity; however, T-cell responses can be suboptimal, impaired or excessive in severe COVID-19 patients. This review focuses on the multifaceted roles of T cells in COVID-19 pathogenesis and rationalizes their significance in eliciting appropriate antiviral immune responses in COVID-19 patients and unexposed individuals. In addition, we summarize the potential therapeutic approaches related to T cells to treat COVID-19 patients. These include adoptive T-cell therapies, vaccines activating T-cell responses, recombinant cytokines, Th1 activators and Th17 blockers, and potential utilization of immune checkpoint inhibitors alone or in combination with anti-inflammatory drugs to improve antiviral T-cell responses against SARS-CoV-2.

Differential gene expression of tumor-infiltrating CD33+ myeloid cells in advanced- versus early-stage colorectal cancer.

Colorectal cancer (CRC) has high mortality rates, especially in patients with advanced disease stages, who often do not respond to therapy. The cellular components of the tumor microenvironment are essentially responsible for dictating disease progression and response to therapy. Expansion of different myeloid cell subsets in CRC tumors has been reported previously. However, tumor-infiltrating myeloid cells have both pro- and anti-tumor roles in disease progression. In this study, we performed transcriptomic profiling of cells of myeloid lineage (CD33+) from bulk CRC tumors at varying disease stages. We identified differentially expressed genes and pathways between CRC patients with advanced stage and early stages. We found that pro-angiogenic and hypoxia-related genes were upregulated, while genes related to immune and inflammatory responses were downregulated in CD33+ myeloid cells from patients with advanced stages, implying that immune cell recruitment and activation could be compromised in advanced disease stages. Moreover, we identified a unique "poor prognosis CD33+ gene signature" by aligning top upregulated and downregulated genes in tumor-infiltrating myeloid cells from our analyses with data from The Cancer Genome Atlas. Our results showed that this gene signature is an independent prognostic indicator for disease-specific survival in CRC patients, potentially reflecting its clinical importance.

Epigenetic regulation of immune checkpoints and T cell exhaustion markers in tumor-infiltrating T cells of colorectal cancer patients.

Aim: To elucidate the epigenetic alterations behind the upregulation of immune checkpoints and T cell exhaustion markers in colorectal cancer (CRC) patients. Materials & methods: mRNA expressions of different immune checkpoint/exhaustion markers were analyzed by quantitative real-time reverse transcriptase PCR and epigenetic investigations were performed using bisulfite sequencing and chromatin immunoprecipitation quantitative PCR. Results: mRNA expressions of PD-1, TIM-3, CTLA-4, PD-L1 and TOX2 were significantly upregulated in CD4+ and CD8+ tumor-infiltrating lymphocytes and bulk CRC tumor tissues. Histone 3 lysine 9 trimethylation was downregulated and histone 3 lysine 4 trimethylation was upregulated in PD-L1 and TOX2 promoters in tumor tissues, suggesting that PD-L1 and TOX2 upregulation in CRC tumors could be mediated by activating histone 3 lysine 4 trimethylation. Conclusion: Epigenetic modifications in promoters of immune checkpoint and T cell exhaustion genes could induce their upregulation, and potentially implicate the use of epigenetic modifiers to enhance antitumor immunity in CRC patients.

An evaluation of sorter induced cell stress (SICS) on peripheral blood mononuclear cells (PBMCs) after different sort conditions - Are your sorted cells getting SICS?

Flow cytometry and fluorescence-activated cell sorting have become invaluable tools to analyze and isolate specific cell populations in a wide range of biomedical research and clinical applications. In countless approaches worldwide, scientists are using single cell analyses to better understand the significance and variation within different cellular populations, and fluorescence-activated cell sorting has become a major technique for cell isolation in both basic and clinical research. However, majority of available cell sorters are pressurized, droplet-based systems, which apply significant environmental pressure and shear stress to cells during sorting. Recently, the flow cytometry community has become increasingly aware about the potential negative effects this could have on sorted cells and the term "sorter induced cell stress" (SICS) has been proposed. However, up to date only a limited number of studies have investigated the effects of cell sorting on cell viability and function. Therefore, solid data on the effects of sheath pressure and nozzle size on survival and function of sorted cells are surprisingly rare. With this in mind, we sorted "CD4+" T-cells and "live" cells from human peripheral blood mononuclear cells (PBMCs) at different sort conditions and analyzed their quality before and after sorting in a series of assays. Here we present our findings in reference to cell viability and cell proliferation following sorting on different instruments (BD FACSAria III SORP and BD FACSJazz), utilizing different nozzle sizes (70 to 100 µm) and sheath pressure settings (20 to 70 psi). The results show no significant differences in cell viability and proliferation after the different tested sort conditions, but rather differences between individual experiments. These findings are evaluated and their potential significance in cell sorting experiments is discussed.

Differential gene expression of tumor-infiltrating CD4+ T cells in advanced versus early stage colorectal cancer and identification of a gene signature of poor prognosis.

Tumor-infiltrating lymphocytes (TILs) play indispensable roles in the progression and response to treatment of solid tumors. However, the prognostic significance of CD4+ TILs is not fully disclosed in cancers generally and in CRC in particular, mainly due to the existence of different functional subsets of CD4+ T cells. We performed transcriptomic profiling of CD4+ TILs isolated from CRC patients in order to identify differentially expressed genes and their functional pathways in early versus advanced disease stages. We found that in advanced stages, genes related to immune and inflammatory responses, in particular Th1-mediated immune response and cytotoxicity-mediated genes, were downregulated; while epigenetic-mediated silencing genes were upregulated. Interestingly, we identified genes, which were steadily upregulated or downregulated in CD4+ TILs with CRC progression from stage I to IV. Additionally, of the top 200 deregulated genes, 43 upregulated and 64 downregulated genes showed similar deregulation trends in the cancer genome atlas CRC dataset. From these 97 deregulated genes, we identified a "poor prognosis CD4 gene signature (ppCD4sig)". Patients with high ppCD4sig score showed shorter disease-specific survival (DSS) and progression-free interval (PFI). The ppCD4sig was an independent prognostic indicator for DSS (HR = 1.73, 95% CI 1.32-2.27, P = 0.0001) and PFI (HR = 1.75, 95% CI 1.3-2.35, P = 0.0016). Additionally, patients at advanced stages and at a younger age (<55 years) were more likely to have a high ppCD4sig score. Altogether, our data provide novel insights and a unique prognostic gene signature of CD4+ TILs in the CRC microenvironment.