Multi-omic analyses in immune cell development with lessons learned from T cell development.

Traditionally, flow cytometry has been the preferred method to characterize immune cells at the single-cell level. Flow cytometry is used in immunology mostly to measure the expression of identifying markers on the cell surface, but-with good antibodies-can also be used to assess the expression of intracellular proteins. The advent of single-cell RNA-sequencing has paved the road to study immune development at an unprecedented resolution. Single-cell RNA-sequencing studies have not only allowed us to efficiently chart the make-up of heterogeneous tissues, including their most rare cell populations, it also increasingly contributes to our understanding how different omics modalities interplay at a single cell resolution. Particularly for investigating the immune system, this means that these single-cell techniques can be integrated to combine and correlate RNA and protein data at the single-cell level. While RNA data usually reveals a large heterogeneity of a given population identified solely by a combination of surface protein markers, the integration of different omics modalities at a single cell resolution is expected to greatly contribute to our understanding of the immune system.

Single-cell immune profiling reveals thymus-seeding populations, T cell commitment, and multilineage development in the human thymus.

T cell development in the mouse thymus has been studied extensively, but less is known regarding T cell development in the human thymus. We used a combination of single-cell techniques and functional assays to perform deep immune profiling of human T cell development, focusing on the initial stages of prelineage commitment. We identified three thymus-seeding progenitor populations that also have counterparts in the bone marrow. In addition, we found that the human thymus physiologically supports the development of monocytes, dendritic cells, and NK cells, as well as limited development of B cells. These results are an important step toward monitoring and guiding regenerative therapies in patients after hematopoietic stem cell transplantation.

Molecular pathways in post-colonoscopy versus detected colorectal cancers: results from a nested case-control study.

BACKGROUND: Post-colonoscopy colorectal cancers (PCCRCs) pose challenges in clinical practice. PCCRCs occur due to a combination of procedural and biological causes. In a nested case-control study, we compared clinical and molecular features of PCCRCs and detected CRCs (DCRCs). METHODS: Whole-genome chromosomal copy number changes and mutation status of genes commonly affected in CRC were examined by low-coverage WGS and targeted sequencing, respectively. MSI and CIMP status was also determined. RESULTS: In total, 122 PCCRCs and 98 DCRCs with high-quality DNA were examined. PCCRCs were more often located proximally (P < 0.001), non-polypoid appearing (P = 0.004), early stage (P = 0.009) and poorly differentiated (P = 0.006). PCCRCs showed significantly less 18q loss (FDR < 0.2), compared to DCRCs. No significant differences in mutations were observed. PCCRCs were more commonly CIMP high (P = 0.014) and MSI (P = 0.029). After correction for tumour location, only less 18q loss remained significant (P = 0.005). CONCLUSION: Molecular features associated with the sessile serrated lesions (SSLs) and non-polypoid colorectal neoplasms (CRNs) are more commonly seen in PCCRCs than in DCRCs. These together with the clinical features observed support the hypothesis that SSLs and non-polypoid CRNs are contributors to the development of PCCRCs. The future focus should be directed at improving the detection and endoscopic removal of these non-polypoid CRN and SSLs. CLINICAL TRIAL REGISTRATION: NTR3093 in the Dutch trial register ( www.trialregister.nl ).

Functional definition of a transcription factor hierarchy regulating T cell lineage commitment.

T cell factor 1 (Tcf1) is the first T cell-specific protein induced by Notch signaling in the thymus, leading to the activation of two major target genes, Gata3 and Bcl11b. Tcf1 deficiency results in partial arrests in T cell development, high apoptosis, and increased development of B and myeloid cells. Phenotypically, seemingly fully T cell-committed thymocytes with Tcf1 deficiency have promiscuous gene expression and an altered epigenetic profile and can dedifferentiate into more immature thymocytes and non-T cells. Restoring Bcl11b expression in Tcf1-deficient cells rescues T cell development but does not strongly suppress the development of non-T cells; in contrast, expressing Gata3 suppresses their development but does not rescue T cell development. Thus, T cell development is controlled by a minimal transcription factor network involving Notch signaling, Tcf1, and the subsequent division of labor between Bcl11b and Gata3, thereby ensuring a properly regulated T cell gene expression program.

Successful Preclinical Development of Gene Therapy for Recombinase-Activating Gene-1-Deficient SCID.

Recombinase-activating gene-1 (RAG1)-deficient severe combined immunodeficiency (SCID) patients lack B and T lymphocytes due to the inability to rearrange immunoglobulin and T cell receptor genes. Gene therapy is an alternative for those RAG1-SCID patients who lack a suitable bone marrow donor. We designed lentiviral vectors with different internal promoters driving codon-optimized RAG1 to ensure optimal expression. We used Rag1 -/- mice as a preclinical model for RAG1-SCID to assess the efficacy of the various vectors. We observed that B and T cell reconstitution directly correlated with RAG1 expression. Mice with low RAG1 expression showed poor immune reconstitution; however, higher expression resulted in phenotypic and functional lymphocyte reconstitution comparable to mice receiving wild-type stem cells. No signs of genotoxicity were found. Additionally, RAG1-SCID patient CD34+ cells transduced with our clinical RAG1 vector and transplanted into NSG mice led to improved human B and T cell development. Considering this efficacy outcome, together with favorable safety data, these results substantiate the need for a clinical trial for RAG1-SCID.

iPSC-Based Modeling of RAG2 Severe Combined Immunodeficiency Reveals Multiple T Cell Developmental Arrests.

RAG2 severe combined immune deficiency (RAG2-SCID) is a lethal disorder caused by the absence of functional T and B cells due to a differentiation block. Here, we generated induced pluripotent stem cells (iPSCs) from a RAG2-SCID patient to study the nature of the T cell developmental blockade. We observed a strongly reduced capacity to differentiate at every investigated stage of T cell development, from early CD7-CD5- to CD4+CD8+. The impaired differentiation was accompanied by an increase in CD7-CD56+CD33+ natural killer (NK) cell-like cells. T cell receptor D rearrangements were completely absent in RAG2SCID cells, whereas the rare T cell receptor B rearrangements were likely the result of illegitimate rearrangements. Repair of RAG2 restored the capacity to induce T cell receptor rearrangements, normalized T cell development, and corrected the NK cell-like phenotype. In conclusion, we succeeded in generating an iPSC-based RAG2-SCID model, which enabled the identification of previously unrecognized disorder-related T cell developmental roadblocks.

IBD-Associated Dysplastic Lesions Show More Chromosomal Instability Than Sporadic Adenomas.

BACKGROUND: Patients with longstanding inflammatory bowel disease (IBD; ie, ulcerative colitis and Crohn's disease) have an increased risk of colorectal cancer (CRC). Due to ongoing inflammation, IBD-associated dysplastic lesions can develop. These lesions have an increased risk to progress to cancer compared with sporadic adenomas, which are also found in these patients. Differentiating between these 2 types of dysplasia remains challenging, both clinically and histologically, while treatment strategies may differ. Therefore, the aim of this study was to investigate molecular alterations associated with colorectal dysplasia to cancer progression in IBD and evaluate to what extent these alterations differ from sporadic adenomas. METHODS: DNA copy number aberrations and mutation analyses of 48 genes were performed by next-generation sequencing in 43 IBD-associated dysplastic lesions, 30 of which were dysplastic and 13 of which were cancers. Results were compared with existing DNA copy number and mutation data from 118 sporadic adenomas and 24 sporadic cancers. RESULTS: Inflammatory bowel disease-associated dysplastic lesions harbor patterns of DNA copy number aberrations comparable to carcinomas, which are rare in sporadic adenomas. TP53 mutation was the most frequent mutation observed in IBD-associated dysplastic lesions and in cancers. FBXW7 was mutated significantly more often in IBD-associated dysplastic lesions than in sporadic adenomas. CONCLUSIONS: Inflammatory bowel disease-associated dysplastic lesions show more DNA copy number aberrations than sporadic adenomas. TP53 and FBXW7 mutations appear to be involved in the development of IBD-associated dysplastic lesions and cancer. These findings indicate that IBD-associated dysplastic lesions are more genomically unstable, possibly reflecting a faster progression toward cancer.

ImSpectR - R package to quantify immune repertoire diversity in spectratype and repertoire sequencing data.

SUMMARY: An effective immune system is characterized by a diverse immune repertoire. There is a strong demand for accurate and quantitative methods to assess the diversity of the immune repertoire for various (pre-)clinical applications, including the diagnosis and prognosis of primary immune deficiencies, or to assess the response to therapy. Current strategies for immune diversity assessment generally comprise the visual inspection of the length distribution of rearranged T- and B-cell receptors. Visual inspections, however, are prone to subjective assessments and thus lead to biases. Here, we introduce ImSpectR, a unified approach to quantify immunodiversity using either spectratype, repertoire sequencing or single cell RNA sequencing data. ImSpectR scores various types of deviations from the expected length distribution and integrates these into one measure, allowing for robust quantitative comparisons of immune diversity across individuals or conditions. AVAILABILITY: R-package is available for download on GitHub at https://github.com/martijn-cordes/ImSpectR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Loss of Chromosome 18q11.2-q12.1 Is Predictive for Survival in Patients With Metastatic Colorectal Cancer Treated With Bevacizumab.

Purpose Patients with metastatic colorectal cancer (mCRC) have limited benefit from the addition of bevacizumab to standard chemotherapy. However, a subset probably benefits substantially, highlighting an unmet clinical need for a biomarker of response to bevacizumab. Previously, we demonstrated that losses of chromosomes 5q34, 17q12, and 18q11.2-q12.1 had a significant correlation with progression-free survival (PFS) in patients with mCRC treated with bevacizumab in the CAIRO2 clinical trial but not in patients who did not receive bevacizumab in the CAIRO trial. This study was designed to validate these findings. Materials and Methods Primary mCRC samples were analyzed from two cohorts of patients who received bevacizumab as first-line treatment; 96 samples from the European multicenter study Angiopredict (APD) and 81 samples from the Italian multicenter study, MOMA. A third cohort of 90 samples from patients with mCRC who did not receive bevacizumab was analyzed. Copy number aberrations of tumor biopsy specimens were measured by shallow whole-genome sequencing and were correlated with PFS, overall survival (OS), and response. Results Loss of chromosome 18q11.2-q12.1 was associated with prolonged PFS most significantly in both the cohorts that received bevacizumab (APD: hazard ratio, 0.54; P = .01; PFS difference, 65 days; MOMA: hazard ratio, 0.55; P = .019; PFS difference, 49 days). A similar association was found for OS and overall response rate in these two cohorts, which became significant when combined with the CAIRO2 cohort. Median PFS in the cohort of patients with mCRC who did not receive bevacizumab and in the CAIRO cohort was similar to that of the APD, MOMA, and CAIRO2 patients without an 18q11.2-q12.1 loss. Conclusion We conclude that the loss of chromosome 18q11.2-q12.1 is consistently predictive for prolonged PFS in patients receiving bevacizumab. The predictive value of this loss is substantiated by a significant gain in OS and overall response rate.

Apelin: A putative novel predictive biomarker for bevacizumab response in colorectal cancer.

Bevacizumab (bvz) is currently employed as an anti-angiogenic therapy across several cancer indications. Bvz response heterogeneity has been well documented, with only 10-15% of colorectal cancer (CRC) patients benefitting in general. For other patients, clinical efficacy is limited and side effects are significant. This reinforces the need for a robust predictive biomarker of response. To identify such a biomarker, we performed a DNA microarray-based transcriptional profiling screen with primary endothelial cells (ECs) isolated from normal and tumour colon tissues. Thirteen separate populations of tumour-associated ECs and 10 of normal ECs were isolated using fluorescence-activated cell sorting. We hypothesised that VEGF-induced genes were overexpressed in tumour ECs; these genes could relate to bvz response and serve as potential predictive biomarkers. Transcriptional profiling revealed a total of 2,610 differentially expressed genes when tumour and normal ECs were compared. To explore their relation to bvz response, the mRNA expression levels of top-ranked genes were examined using quantitative PCR in 30 independent tumour tissues from CRC patients that received bvz in the adjuvant setting. These analyses revealed that the expression of MMP12 and APLN mRNA was significantly higher in bvz non-responders compared to responders. At the protein level, high APLN expression was correlated with poor progression-free survival in bvz-treated patients. Thus, high APLN expression may represent a novel predictive biomarker for bvz unresponsiveness.

Outcome of Colorectal Cancer Patients Treated with Combination Bevacizumab Therapy: A Pooled Retrospective Analysis of Three European Cohorts from the Angiopredict Initiative.

BACKGROUND/AIMS: This study is aimed at analyzing the survival rates and prognostic factors of stage IV colorectal cancer patients from 3 European cohorts undergoing combination chemotherapy with bevacizumab. METHODS: Progression free-survival (PFS) and overall survival (OS) were analyzed in 172 patients using the Kaplan-Meier method and uni- and multivariable Cox proportional hazards regression models. RESULTS: The median PFS was 9.7 and the median OS 27.4 months. Patients treated at centers in Germany (n = 97), Ireland (n = 32), and The Netherlands (n = 43) showed a median PFS of 9.9, 9.2, and 9.7 months, OS of 34.0, 20.5, and 25.1 months, respectively. Patients >65 years had a significantly shorter PFS (9.5 vs. 9.8 months) but not OS (27.4 vs. 27.5 months) than younger patients. High tumor grade (G3/4) was associated with a shorter PFS, T4 classification with both shorter PFS and OS. Fluoropyrimidine (FP) chemotherapy backbones (doublets and single) had comparable outcomes, while patients not receiving FP backbones had a shorter PFS. In multivariable analysis, age and non-FP backbone were associated with inferior PFS, T4 classification and therapy line >2nd were significantly associated with poor PFS and OS. CONCLUSION: The observed survival rates confirm previous studies and demonstrate reproducible benefits of combination bevacizumab regimens. Classification T4, non-FP chemotherapy backbone, and age >65 were associated with inferior outcome.

Overexpression of LMO2 causes aberrant human T-Cell development in vivo by three potentially distinct cellular mechanisms.

Overexpression of LMO2 is known to be one of the causes of T-cell acute lymphoblastic leukemia (T-ALL) development; however, the mechanisms behind its oncogenic activity are incompletely understood. LMO2-overexpressing transgenic mouse models suggest an accumulation of immature T-cell progenitors in the thymus as the main preleukemic event. The effects of LMO2 overexpression on human T-cell development in vivo are unknown. Here, we report studies of a humanized mouse model transplanted with LMO2-transduced human hematopoietic stem/progenitor cells. The effects of LMO2 overexpression were confined to the T-cell lineage; however, initially, multipotent cells were transduced. Three effects of LMO2 on human T-cell development were observed: (1) a block at the double-negative/immature single-positive stage, (2) an accumulation of CD4(+)CD8(+) double-positive CD3(-) cells, and (3) an altered CD8/CD4 ratio with enhanced peripheral T lymphocytes. Microarray analysis of sorted double-positive cells overexpressing LMO2 led to the identification of an LMO2 gene set that clustered with human T-ALL patient samples of the described "proliferative" cluster. In this article, we demonstrate previously unrecognized mechanisms by which LMO2 alters human T-cell development in vivo; these mechanisms correlate with human T-ALL leukemogenesis.

Spatial and temporal evolution of distal 10q deletion, a prognostically unfavorable event in diffuse low-grade gliomas.

BACKGROUND: The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behavior and can thereby complement the prognostically favorable 1p/19q co-deletion. RESULTS: Genome-wide, 50 base pair single-end sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analyzed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors. CONCLUSIONS: CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.

PePPER: a webserver for prediction of prokaryote promoter elements and regulons.

BACKGROUND: Accurate prediction of DNA motifs that are targets of RNA polymerases, sigma factors and transcription factors (TFs) in prokaryotes is a difficult mission mainly due to as yet undiscovered features in DNA sequences or structures in promoter regions. Improved prediction and comparison algorithms are currently available for identifying transcription factor binding sites (TFBSs) and their accompanying TFs and regulon members. RESULTS: We here extend the current databases of TFs, TFBSs and regulons with our knowledge on Lactococcus lactis and developed a webserver for prediction, mining and visualization of prokaryote promoter elements and regulons via a novel concept. This new approach includes an all-in-one method of data mining for TFs, TFBSs, promoters, and regulons for any bacterial genome via a user-friendly webserver. We demonstrate the power of this method by mining WalRK regulons in Lactococci and Streptococci and, vice versa, use L. lactis regulon data (CodY) to mine closely related species. CONCLUSIONS: The PePPER webserver offers, besides the all-in-one analysis method, a toolbox for mining for regulons, promoters and TFBSs and accommodates a new L. lactis regulon database in addition to already existing regulon data. Identification of putative regulons and full annotation of intergenic regions in any bacterial genome on the basis of existing knowledge on a related organism can now be performed by biologists and it can be done for a wide range of regulons. On the basis of the PePPER output, biologist can design experiments to further verify the existence and extent of the proposed regulons. The PePPER webserver is freely accessible at http://pepper.molgenrug.nl.