Tumor-associated T cell receptor repertoires in low- and high-grade gliomas

Glioblastoma (GBM) remains prognostically dismal, with care centered on resection, motivating research into novel therapies. Although inducing anti-tumor immunity remains an attractive target for therapeutic and preventative intervention, the interplay between evolving dysregulation of the glioma microenvironment and T cell inefficacy remains poorly understood. In our murine model of proneural glioma, retroviral delivery of PDGF and cre-mediated knockout of PTEN in glial progenitors of adult C57BL/6 gives rise to slow-growing tumors, which were harvested at early- mid- and late-stage progression timepoints following induction, along with peripheral blood. From human patients, tissue from low- and high-grade glioma resections and corresponding peripheral lymphocytes were cryofrozen during surgery at New York Presbyterian-CUMC. For both species, we employed a commercially available primer set (iRepertoire) for nested PCR of the complementarity-determining region 3 (CDR3) of the TCR-alpha and TCR-beta chains from the T cell RNA, followed by next-generation sequencing on an Illumina MiSeq. We developed a computational pipeline for mapping TCR cassettes, in silico translation, pairwise analysis of tissue/periphery per subject, and error analysis. In the murine model, we observe that at late-stage, the intratumoral TCR repertoire diverges significantly from the peripheral, including dramatic expansion of single tumor-associated CDR3s, while the peripheral repertoire itself diverges from those of healthy mice. In both human patients and mice, we observed tumor-associated CDR3s, disproportionately abundant in tumor tissue compared to the corresponding peripheral blood, at both the amino acid and nucleotide level. In human samples we observed tumor-specific TCR expansions that were associated with particular functional subsets (CD8+, CD4+, Treg, NKT). Sequence-level study of the TCR repertoire promises new insight into the scope of glioma immunosuppression, especially systemic effects which remain elusive and the origins of intratumoral suppressive populations, and holds the potential for immunotherapeutic interventions, non-invasive diagnostics, and direct assessment of global responses to immunotherapy

TCR repertoire divergence reflects micro-environmental immune phenotypes in glioma

Background & significance: Glioblastoma (GBM) remains prognostically dismal, with only modest gains in mean survival time with chemo- and radiotherapy motivating research into reversing its characteristic local and systemic immunosuppression with precision in this high-risk tissue. While whole-repertoire amplification of the TCR repertoire allows unprecedented depth regarding the potentiation of anti-tumor responses, most studies utilize TCRseq for monitoring reactivity to specific tumor antigens, or the identities of particular TCRs as biomarkers. In this study, we have utilized whole-repertoire analysis to describe the relationship between intra-tumoral T cells and peripheral circulation, and leverage mutual information between gene expression and the behavior of the T cell population to characterize glioma-reactive states, driven by the gene expression of the principal resident monocyte population, and perturbable by immunological interventions. Methods & results: From resected tumor tissue and peripheral lymphocytes of low- and high-grade human glioma patients, TCRseq libraries were generated using reverse transcription and nested PCR (iRepertoire 1) of the complementarity-determining region 3 (CDR3) of the TCR-alpha and TCR-beta chains, then sequenced on an Illumina MiSeq. We developed a computational pipeline for mapping TCR cassettes, in silico translation, and sequence error correction from these libraries, enabling sensitive calculation of tumor-infiltrating lymphocyte (TIL) and peripheral TCR diversity (Shannon entropy) 2, as well as the divergence (Jensen-Shannon divergence metric) between the two T cell populations. By integrating amino acid identity and V-J cassette combination, we observed varying levels of divergence between the TIL and peripheral lymphocytes of glioma patients, and changes in this divergence over tumor progression in a PDGF-driven murine model. Correlation of these properties with tumor tissue RNA profiling, by differential gene expression and mutual-information gene ontology, revealed an association between tumor growth and high blood-brain TCR divergence - particularly in amino-acid sequence, suggesting antigen-driven selection - while high expression of inflammatory and certain immune pathway markers computationally attributed to microglia 3 were anti-correlated with divergence. Preliminary murine experiments suggest that TCR divergence can be altered by induction and blockade of cytokine-mediated activation of these pathways. Conclusion: The expression of a subset of microglia-associated genes appears to describe micro-environmental states which are strongly tied to the tumor-specificity of the intra-tumoral TCR repertoire, complementary to the tumor-centric classifications of TCGA. TCRseq-based profiling not only promises to inform tailoring of local and systemic immunotherapy to target the most relevant immunosuppressive mechanisms, but may also provide non-invasive assessment of the intra-tumoral environment for refined diagnosis and monitoring during clinical trials

MRI-localized biopsies reveal subtype-specific differences in molecular and cellular composition at the margins of glioblastoma

Glioblastomas (GBMs) diffusely infiltrate the brain, making complete removal by surgical resection impossible. The mixture of neoplastic and nonneoplastic cells that remain after surgery form the biological context for adjuvant therapeutic intervention and recurrence. We performed RNA-sequencing (RNA-seq) and histological analysis on radiographically guided biopsies taken from different regions of GBM and showed that the tissue contained within the contrast-enhancing (CE) core of tumors have different cellular and molecular compositions compared with tissue from the nonenhancing (NE) margins of tumors. Comparisons with the The Cancer Genome Atlas dataset showed that the samples from CE regions resembled the proneural, classical, or mesenchymal subtypes of GBM, whereas the samples from the NE regions predominantly resembled the neural subtype. Computational deconvolution of the RNA-seq data revealed that contributions from nonneoplastic brain cells significantly influence the expression pattern in the NE samples. Gene ontology analysis showed that the cell type-specific expression patterns were functionally distinct and highly enriched in genes associated with the corresponding cell phenotypes. Comparing the RNA-seq data from the GBM samples to that of nonneoplastic brain revealed that the differentially expressed genes are distributed across multiple cell types. Notably, the patterns of cell type-specific alterations varied between the different GBM subtypes: the NE regions of proneural tumors were enriched in oligodendrocyte progenitor genes, whereas the NE regions of mesenchymal GBM were enriched in astrocytic and microglial genes. These subtypespecific patterns provide new insights into molecular and cellular composition of the infiltrative margins of GBM

Developmental Regulation of Genes Encoding Universal Stress Proteins in Schistosoma mansoni

The draft nuclear genome sequence of the snail-transmitted, dimorphic, parasitic, platyhelminth Schistosoma mansoni revealed eight genes encoding proteins that contain the Universal Stress Protein (USP) domain. Schistosoma mansoni is a causative agent of human schistosomiasis, a severe and debilitating Neglected Tropical Disease (NTD) of poverty, which is endemic in at least 76 countries. The availability of the genome sequences of Schistosoma species presents opportunities for bioinformatics and genomics analyses of associated gene families that could be targets for understanding schistosomiasis ecology, intervention, prevention and control. Proteins with the USP domain are known to provide bacteria, archaea, fungi, protists and plants with the ability to respond to diverse environmental stresses. In this research investigation, the functional annotations of the USP genes and predicted nucleotide and protein sequences were initially verified. Subsequently, sequence clusters and distinctive features of the sequences were determined. A total of twelve ligand binding sites were predicted based on alignment to the ATP-binding universal stress protein from Methanocaldococcus jannaschii. In addition, six USP sequences showed the presence of ATP-binding motif residues indicating that they may be regulated by ATP. Public domain gene expression data and RT-PCR assays confirmed that all the S. mansoni USP genes were transcribed in at least one of the developmental life cycle stages of the helminth. Six of these genes were up-regulated in the miracidium, a free-swimming stage that is critical for transmission to the snail intermediate host. It is possible that during the intra-snail stages, S. mansoni gene transcripts for universal stress proteins are low abundant and are induced to perform specialized functions triggered by environmental stressors such as oxidative stress due to hydrogen peroxide that is present in the snail hemocytes. This report serves to catalyze the formation of a network of researchers to understand the function and regulation of the universal stress proteins encoded in genomes of schistosomes and their snail intermediate hosts

Genetic Analysis of Completely Sequenced Disease-Associated MHC Haplotypes Identifies Shuffling of Segments in Recent Human History

The major histocompatibility complex (MHC) is recognised as one of the most important genetic regions in relation to common human disease. Advancement in identification of MHC genes that confer susceptibility to disease requires greater knowledge of sequence variation across the complex. Highly duplicated and polymorphic regions of the human genome such as the MHC are, however, somewhat refractory to some whole-genome analysis methods. To address this issue, we are employing a bacterial artificial chromosome (BAC) cloning strategy to sequence entire MHC haplotypes from consanguineous cell lines as part of the MHC Haplotype Project. Here we present 4.25 Mb of the human haplotype QBL (HLA-A26-B18-Cw5-DR3-DQ2) and compare it with the MHC reference haplotype and with a second haplotype, COX (HLA-A1-B8-Cw7-DR3-DQ2), that shares the same HLA-DRB1, -DQA1, and -DQB1 alleles. We have defined the complete gene, splice variant, and sequence variation contents of all three haplotypes, comprising over 259 annotated loci and over 20,000 single nucleotide polymorphisms (SNPs). Certain coding sequences vary significantly between different haplotypes, making them candidates for functional and disease-association studies. Analysis of the two DR3 haplotypes allowed delineation of the shared sequence between two HLA class II–related haplotypes differing in disease associations and the identification of at least one of the sites that mediated the original recombination event. The levels of variation across the MHC were similar to those seen for other HLA-disparate haplotypes, except for a 158-kb segment that contained the HLA-DRB1, -DQA1, and -DQB1 genes and showed very limited polymorphism compatible with identity-by-descent and relatively recent common ancestry (<3,400 generations). These results indicate that the differential disease associations of these two DR3 haplotypes are due to sequence variation outside this central 158-kb segment, and that shuffling of ancestral blocks via recombination is a potential mechanism whereby certain DR–DQ allelic combinations, which presumably have favoured immunological functions, can spread across haplotypes and populations

Psychometric Properties of the Chinese Version of the Perceived Stress Scale in Policewomen

BACKGROUND: The 10-item Perceived Stress Scale (PSS-10) is one of most widely used instruments to measure a global level of perceived stress in a range of clinical and research settings. This study was conducted to examine the psychometric properties of the Simplified Chinese version of the PSS-10 in policewomen. METHODOLOGY: A total of 240 policewomen were recruited in this study. The Simplified Chinese versions of the PSS-10, the Beck Depression Inventory Revised (BDI-II), and the Beck Anxiety Inventory (BAI) were administered to all participants, and 36 of the participants were re-tested two weeks after the initial testing. PRINCIPAL FINDINGS: The overall Cronbach's alpha was 0.86, and the test-retest reliability coefficient was 0.68. Exploratory Factor Analysis (EFA) yielded 2 factors with eigenvalues of 4.76 and 1.48, accounting for 62.41% of variance. Factor 1 consisted of 6 items representing "negative feelings"; whereas Factor 2 consisted of 4 items representing "positive feelings". The item loadings ranged from 0.72 to 0.83. The Confirmatory factor analysis (CFA) indicated a very good fit of this two-factor model to this sample. The PSS-10 significantly correlated with both BDI-II and BAI, indicating an acceptable concurrent validity. CONCLUSIONS: The Simplified Chinese version of the PSS-10 demonstrated adequate psychometric properties for evaluating stress levels. The results support its use among the Chinese population

Effects of antiplatelet therapy on stroke risk by brain imaging features of intracerebral haemorrhage and cerebral small vessel diseases: subgroup analyses of the RESTART randomised, open-label trial

Background Findings from the RESTART trial suggest that starting antiplatelet therapy might reduce the risk of recurrent symptomatic intracerebral haemorrhage compared with avoiding antiplatelet therapy. Brain imaging features of intracerebral haemorrhage and cerebral small vessel diseases (such as cerebral microbleeds) are associated with greater risks of recurrent intracerebral haemorrhage. We did subgroup analyses of the RESTART trial to explore whether these brain imaging features modify the effects of antiplatelet therapy

Transcription and Expression of Plasmodium falciparum Histidine-Rich Proteins in Different Stages and Strains: Implications for Rapid Diagnostic Tests

Background: Although rapid diagnostic tests (RDTs) for Plasmodium falciparum infection that target histidine rich protein 2 (PfHRP2) are generally sensitive, their performance has been reported to be variable. One possible explanation for variable test performance is differences in expression level of PfHRP in different parasite isolates. Methods: Total RNA and protein were extracted from synchronised cultures of 7 P. falciparum lines over 5 time points of the life cycle, and from synchronised ring stages of 10 falciparum lines. Using quantitative real-time polymerase chain reaction, Western blot analysis and ELISA we investigated variations in the transcription and protein levels of pfhrp2, pfhrp3 and PfHRP respectively in the different parasite lines, over the parasite intraerythrocytic life cycle. Results: Transcription of pfhrp2 and pfhrp3 in different parasite lines over the parasite life cycle was observed to vary relative to the control parasite K1. In some parasite lines very low transcription of these genes was observed. The peak transcription was observed in ring-stage parasites. Pfhrp2 transcription was observed to be consistently higher than pfhrp3 transcription within parasite lines. The intraerythrocytic lifecycle stage at which the peak level of protein was present varied across strains. Total protein levels were more constant relative to total mRNA transcription, however a maximum 24 fold difference in expression at ring-stage parasites relative to the K1 strain was observed. Conclusions: The levels of transcription of pfhrp2 and pfhrp3, and protein expression of PfHRP varied between different P. falciparum strains. This variation may impact on the detection sensitivity of PfHRP2-detecting RDTs

Mitochondrial Structure, Function and Dynamics Are Temporally Controlled by c-Myc

Although the c-Myc (Myc) oncoprotein controls mitochondrial biogenesis and multiple enzymes involved in oxidative phosphorylation (OXPHOS), the coordination of these events and the mechanistic underpinnings of their regulation remain largely unexplored. We show here that re-expression of Myc in myc−/− fibroblasts is accompanied by a gradual accumulation of mitochondrial biomass and by increases in membrane polarization and mitochondrial fusion. A correction of OXPHOS deficiency is also seen, although structural abnormalities in electron transport chain complexes (ETC) are not entirely normalized. Conversely, the down-regulation of Myc leads to a gradual decrease in mitochondrial mass and a more rapid loss of fusion and membrane potential. Increases in the levels of proteins specifically involved in mitochondrial fission and fusion support the idea that Myc affects mitochondrial mass by influencing both of these processes, albeit favoring the latter. The ETC defects that persist following Myc restoration may represent metabolic adaptations, as mitochondrial function is re-directed away from producing ATP to providing a source of metabolic precursors demanded by the transformed cell

Resolving taxonomic uncertainty in vulnerable elasmobranchs : are the Madeira skate (Raja maderensis) and the thornback ray (Raja clavata) distinct species?

Skates and rays constitute the most speciose group of chondrichthyan fishes, yet are characterised by remarkable levels of morphological and ecological conservatism. They can be challenging to identify, which makes monitoring species compositions for fisheries management purposes problematic. Owing to their slow growth and low fecundity, skates are vulnerable to exploitation and species exhibiting endemism or limited ranges are considered to be the most at risk. The Madeira skate Raja maderensis is endemic and classified as ‘Data Deficient’ by the IUCN, yet its taxonomic distinctiveness from the morphologically similar and more wide-ranging thornback ray Raja clavate is unresolved. This study evaluated the sequence divergence of both the variable control region and cytochrome oxidase I ‘DNA barcode’ gene of the mitochondrial genome to elucidate the genetic differentiation of specimens identified as R. maderensis and R. clavate collected across much of their geographic ranges. Genetic evidence was insufficient to support the different species designations. However regardless of putative species identification, individuals occupying waters around the Azores and North African Seamounts represent an evolutionarily significant unit worthy of special consideration for conservation management