NKG2A is a late immune checkpoint on CD8 T cells and marks repeated stimulation and cell division.

The surface inhibitory receptor NKG2A forms heterodimers with the invariant CD94 chain and is expressed on a subset of activated CD8 T cells. As antibodies to block NKG2A are currently tested in several efficacy trials for different tumor indications, it is important to characterize the NKG2A+ CD8 T cell population in the context of other inhibitory receptors. Here we used a well-controlled culture system to study the kinetics of inhibitory receptor expression. Naïve mouse CD8 T cells were synchronously and repeatedly activated by artificial antigen presenting cells in the presence of the homeostatic cytokine IL-7. The results revealed NKG2A as a late inhibitory receptor, expressed after repeated cognate antigen stimulations. In contrast, the expression of PD-1, TIGIT and LAG-3 was rapidly induced, hours after first contact and subsequently down regulated during each resting phase. This late, but stable expression kinetics of NKG2A was most similar to that of TIM-3 and CD39. Importantly, single-cell transcriptomics of human tumor-infiltrating lymphocytes (TILs) showed indeed that these receptors were often coexpressed by the same CD8 T cell cluster. Furthermore, NKG2A expression was associated with cell division and was promoted by TGF-ß in vitro, although TGF-ß signaling was not necessary in a mouse tumor model in vivo. In summary, our data show that PD-1 reflects recent TCR triggering, but that NKG2A is induced after repeated antigen stimulations and represents a late inhibitory receptor. Together with TIM-3 and CD39, NKG2A might thus mark actively dividing tumor-specific TILs.

Hippocampal glucocorticoid target genes associated with enhancement of memory consolidation.

Glucocorticoids enhance memory consolidation of emotionally arousing events via largely unknown molecular mechanisms. This glucocorticoid effect on the consolidation process also requires central noradrenergic neurotransmission. The intracellular pathways of these two stress mediators converge on two transcription factors: the glucocorticoid receptor (GR) and phosphorylated cAMP response element-binding protein (pCREB). We therefore investigated, in male rats, whether glucocorticoid effects on memory are associated with genomic interactions between the GR and pCREB in the hippocampus. In a two-by-two design, object exploration training or no training was combined with post-training administration of a memory-enhancing dose of corticosterone or vehicle. Genomic effects were studied by chromatin immunoprecipitation followed by sequencing (ChIP-seq) of GR and pCREB 45 min after training and transcriptome analysis after 3 hr. Corticosterone administration induced differential GR DNA-binding and regulation of target genes within the hippocampus, largely independent of training. Training alone did not result in long-term memory nor did it affect GR or pCREB DNA-binding and gene expression. No strong evidence was found for an interaction between GR and pCREB. Combination of the GR DNA-binding and transcriptome data identified a set of novel, likely direct, GR target genes that are candidate mediators of corticosterone effects on memory consolidation. Cell-specific expression of the identified target genes using single-cell expression data suggests that the effects of corticosterone reflect in part non-neuronal cells. Together, our data identified new GR targets associated with memory consolidation that reflect effects in both neuronal and non-neuronal cells.

Templated insertions at VD and DJ junctions create unique B-cell receptors in the healthy B-cell repertoire.

RAG complexes recognise (cryptic) RSS sites both in and outside immunoglobulin sites. Excision circles may be reinserted into V(D)J rearrangements as long templated insertions to diversify the adaptive immune repertoire. We show that such VDJ with templated insertions are incidentally found in the repertoire of healthy donors.

Statistical method for modeling sequencing data from different technologies in longitudinal studies with application to Huntington disease.

Advancement of gene expression measurements in longitudinal studies enables the identification of genes associated with disease severity over time. However, problems arise when the technology used to measure gene expression differs between time points. Observed differences between the results obtained at different time points can be caused by technical differences. Modeling the two measurements jointly over time might provide insight into the causes of these different results. Our work is motivated by a study of gene expression data of blood samples from Huntington disease patients, which were obtained using two different sequencing technologies. At time point 1, DeepSAGE technology was used to measure the gene expression, with a subsample also measured using RNA-Seq technology. At time point 2, all samples were measured using RNA-Seq technology. Significant associations between gene expression measured by DeepSAGE and disease severity using data from the first time point could not be replicated by the RNA-Seq data from the second time point. We modeled the relationship between the two sequencing technologies using the data from the overlapping samples. We used linear mixed models with either DeepSAGE or RNA-Seq measurements as the dependent variable and disease severity as the independent variable. In conclusion, (1) for one out of 14 genes, the initial significant result could be replicated with both technologies using data from both time points; (2) statistical efficiency is lost due to disagreement between the two technologies, measurement error when predicting gene expressions, and the need to include additional parameters to account for possible differences.

Exacerbated inflammatory signaling underlies aberrant response to BMP9 in pulmonary arterial hypertension lung endothelial cells.

Imbalanced transforming growth factor beta (TGFß) and bone morphogenetic protein (BMP) signaling are postulated to favor a pathological pulmonary endothelial cell (EC) phenotype in pulmonary arterial hypertension (PAH). BMP9 is shown to reinstate BMP receptor type-II (BMPR2) levels and thereby mitigate hemodynamic and vascular abnormalities in several animal models of pulmonary hypertension (PH). Yet, responses of the pulmonary endothelium of PAH patients to BMP9 are unknown. Therefore, we treated primary PAH patient-derived and healthy pulmonary ECs with BMP9 and observed that stimulation induces transient transcriptional signaling associated with the process of endothelial-to-mesenchymal transition (EndMT). However, solely PAH pulmonary ECs showed signs of a mesenchymal trans-differentiation characterized by a loss of VE-cadherin, induction of transgelin (SM22α), and reorganization of the cytoskeleton. In the PAH cells, a prolonged EndMT signaling was found accompanied by sustained elevation of pro-inflammatory, pro-hypoxic, and pro-apoptotic signaling. Herein we identified interleukin-6 (IL6)-dependent signaling to be the central mediator required for the BMP9-induced phenotypic change in PAH pulmonary ECs. Furthermore, we were able to target the BMP9-induced EndMT process by an IL6 capturing antibody that normalized autocrine IL6 levels, prevented mesenchymal transformation, and maintained a functional EC phenotype in PAH pulmonary ECs. In conclusion, our results show that the BMP9-induced aberrant EndMT in PAH pulmonary ECs is dependent on exacerbated pro-inflammatory signaling mediated through IL6.

RNA-Seq in 296 phased trios provides a high-resolution map of genomic imprinting.

BACKGROUND: Identification of imprinted genes, demonstrating a consistent preference towards the paternal or maternal allelic expression, is important for the understanding of gene expression regulation during embryonic development and of the molecular basis of developmental disorders with a parent-of-origin effect. Combining allelic analysis of RNA-Seq data with phased genotypes in family trios provides a powerful method to detect parent-of-origin biases in gene expression. RESULTS: We report findings in 296 family trios from two large studies: 165 lymphoblastoid cell lines from the 1000 Genomes Project and 131 blood samples from the Genome of the Netherlands (GoNL) participants. Based on parental haplotypes, we identified > 2.8 million transcribed heterozygous SNVs phased for parental origin and developed a robust statistical framework for measuring allelic expression. We identified a total of 45 imprinted genes and one imprinted unannotated transcript, including multiple imprinted transcripts showing incomplete parental expression bias that was located adjacent to strongly imprinted genes. For example, PXDC1, a gene which lies adjacent to the paternally expressed gene FAM50B, shows a 2:1 paternal expression bias. Other imprinted genes had promoter regions that coincide with sites of parentally biased DNA methylation identified in the blood from uniparental disomy (UPD) samples, thus providing independent validation of our results. Using the stranded nature of the RNA-Seq data in lymphoblastoid cell lines, we identified multiple loci with overlapping sense/antisense transcripts, of which one is expressed paternally and the other maternally. Using a sliding window approach, we searched for imprinted expression across the entire genome, identifying a novel imprinted putative lncRNA in 13q21.2. Overall, we identified 7 transcripts showing parental bias in gene expression which were not reported in 4 other recent RNA-Seq studies of imprinting. CONCLUSIONS: Our methods and data provide a robust and high-resolution map of imprinted gene expression in the human genome.

RNA sequencing data integration reveals an miRNA interactome of osteoarthritis cartilage.

OBJECTIVE: To uncover the microRNA (miRNA) interactome of the osteoarthritis (OA) pathophysiological process in the cartilage. METHODS: We performed RNA sequencing in 130 samples (n=35 and n=30 pairs for messenger RNA (mRNA) and miRNA, respectively) on macroscopically preserved and lesioned OA cartilage from the same patient and performed differential expression (DE) analysis of miRNA and mRNAs. To build an OA-specific miRNA interactome, a prioritisation scheme was applied based on inverse Pearson's correlations and inverse DE of miRNAs and mRNAs. Subsequently, these were filtered by those present in predicted (TargetScan/microT-CDS) and/or experimentally validated (miRTarBase/TarBase) public databases. Pathway enrichment analysis was applied to elucidate OA-related pathways likely mediated by miRNA regulatory mechanisms. RESULTS: We found 142 miRNAs and 2387 mRNAs to be differentially expressed between lesioned and preserved OA articular cartilage. After applying prioritisation towards likely miRNA-mRNA targets, a regulatory network of 62 miRNAs targeting 238 mRNAs was created. Subsequent pathway enrichment analysis of these mRNAs (or genes) elucidated that genes within the 'nervous system development' are likely mediated by miRNA regulatory mechanisms (familywise error=8.4×10-5). Herein NTF3 encodes neurotrophin-3, which controls survival and differentiation of neurons and which is closely related to the nerve growth factor. CONCLUSIONS: By an integrated approach of miRNA and mRNA sequencing data of OA cartilage, an OA miRNA interactome and related pathways were elucidated. Our functional data demonstrated interacting levels at which miRNA affects expression of genes in the cartilage and exemplified the complexity of functionally validating a network of genes that may be targeted by multiple miRNAs.

The distinct transcriptomes of slow and fast adult muscles are delineated by noncoding RNAs.

Adult muscles have a vast adaptation capacity, enabling function switches in response to altered conditions. During intensive physical activity, disease, or aging, adult skeletal muscles change and adjust their functions. The competence to adjust varies among muscles. Muscle-specific molecular mechanisms in healthy and normal conditions could designate changes in physiologic and pathologic conditions. We generated deep mRNA-sequencing data in adult fast and slow mouse muscles, and applying paired analysis, we identified that the muscle-specific signatures are composed of half of the muscle transcriptome. The fast muscles showed a more compact gene network that is concordant with homogenous myofiber typing, compared with the pattern in the slow muscle. The muscle-specific mRNA landscape did not correlate with alternative spicing, alternative polyadenylation, or the expression of muscle transcription factor gene networks. However, we found significant correlation between the differentially expressed noncoding RNAs, microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) and their target genes. More than 25% of the genes expressed in a muscle-specific fashion were found to be targets of muscle-specific miRNAs and lncRNAs. We suggest that muscle-specific miRNAs and lncRNAs contribute to the establishment of muscle-specific transcriptomes in adult muscles.-Raz, V., Riaz, M., Tatum, Z., Kielbasa, S. M., 't Hoen, P. A. C. The distinct transcriptomes of slow and fast adult muscles are delineated by noncoding RNAs.

Integrating omics datasets with the OmicsPLS package.

BACKGROUND: With the exponential growth in available biomedical data, there is a need for data integration methods that can extract information about relationships between the data sets. However, these data sets might have very different characteristics. For interpretable results, data-specific variation needs to be quantified. For this task, Two-way Orthogonal Partial Least Squares (O2PLS) has been proposed. To facilitate application and development of the methodology, free and open-source software is required. However, this is not the case with O2PLS. RESULTS: We introduce OmicsPLS, an open-source implementation of the O2PLS method in R. It can handle both low- and high-dimensional datasets efficiently. Generic methods for inspecting and visualizing results are implemented. Both a standard and faster alternative cross-validation methods are available to determine the number of components. A simulation study shows good performance of OmicsPLS compared to alternatives, in terms of accuracy and CPU runtime. We demonstrate OmicsPLS by integrating genetic and glycomic data. CONCLUSIONS: We propose the OmicsPLS R package: a free and open-source implementation of O2PLS for statistical data integration. OmicsPLS is available at https://cran.r-project.org/package=OmicsPLS and can be installed in R via install.packages("OmicsPLS").

RNA-seq analysis of Lgr6+ stem cells and identification of an Lgr6 isoform.

We studied Lgr6+ stem cells in experimental UV carcinogenesis in hairless mice. For further characterization through RNA-seq, these stem cells were isolated by FACS from transgenic hairless mice bearing an EGFP-Ires-CreERT2 reporter cassette inserted into exon 1 of the Lgr6 gene (purity confirmed by human ERT2 expression). Between Lgr6/EGFP+ and Lgr6/EGFP- basal cells (Tg/wt), 682 RNAs were differentially expressed, indicating stemness and expression of cancer-related pathways in Lgr6/EGFP+ cells. We discovered that suspected "Lgr6 null" mice (Tg/Tg) expressed RNA of an Lgr6 isoform (delta-Lgr6, lacking 74 N-terminal aa) which could be functional and explain the lack of a phenotype.

New function of the myostatin/activin type I receptor (ALK4) as a mediator of muscle atrophy and muscle regeneration.

Skeletal muscle fibrosis and impaired muscle regeneration are major contributors to muscle wasting in Duchenne muscular dystrophy (DMD). Muscle growth is negatively regulated by myostatin (MSTN) and activins. Blockage of these pathways may improve muscle quality and function in DMD. Antisense oligonucleotides (AONs) were designed specifically to block the function of ALK4, a key receptor for the MSTN/activin pathway in skeletal muscle. AON-induced exon skipping resulted in specific Alk4 down-regulation, inhibition of MSTN activity, and increased myoblast differentiation in vitro Unexpectedly, a marked decrease in muscle mass (10%) was found after Alk4 AON treatment in mdx mice. In line with in vitro results, muscle regeneration was stimulated, and muscle fiber size decreased markedly. Notably, when Alk4 was down-regulated in adult wild-type mice, muscle mass decreased even more. RNAseq analysis revealed dysregulated metabolic functions and signs of muscle atrophy. We conclude that ALK4 inhibition increases myogenesis but also regulates the tight balance of protein synthesis and degradation. Therefore, caution must be used when developing therapies that interfere with MSTN/activin pathways.-Pasteuning-Vuhman, S., Boertje-van der Meulen, J. W., van Putten, M., Overzier, M., ten Dijke, P., Kielbasa, S. M., Arindrarto, W., Wolterbeek, R., Lezhnina, K. V., Ozerov, I. V., Aliper, A. M., Hoogaars, W. M., Aartsma-Rus, A., Loomans, C. J. M. New function of the myostatin/activin type I receptor (ALK4) as a mediator of muscle atrophy and muscle regeneration.

Two Plasmodium 6-Cys family-related proteins have distinct and critical roles in liver-stage development.

The 10 Plasmodium 6-Cys proteins have critical roles throughout parasite development and are targets for antimalaria vaccination strategies. We analyzed the conserved 6-cysteine domain of this family and show that only the last 4 positionally conserved cysteine residues are diagnostic for this domain and identified 4 additional "6-Cys family-related" proteins. Two of these, sequestrin and B9, are critical to Plasmodium liver-stage development. RT-PCR and immunofluorescence assays show that B9 is translationally repressed in sporozoites and is expressed after hepatocyte invasion where it localizes to the parasite plasma membrane. Mutants lacking B9 expression in the rodent malaria parasites P. berghei and P. yoelii and the human parasite P. falciparum developmentally arrest in hepatocytes. P. berghei mutants arrest in the livers of BALB/c (100%) and C57BL6 mice (>99.9%), and in cultures of Huh7 human-hepatoma cell line. Similarly, P. falciparum mutants while fully infectious to primary human hepatocytes abort development 3 d after infection. This growth arrest is associated with a compromised parasitophorous vacuole membrane a phenotype similar to, but distinct from, mutants lacking the 6-Cys sporozoite proteins P52 and P36. Our results show that 6-Cys proteins have critical but distinct roles in establishment and maintenance of a parasitophorous vacuole and subsequent liver-stage development.

Adaptive seeds tame genomic sequence comparison.

The main way of analyzing biological sequences is by comparing and aligning them to each other. It remains difficult, however, to compare modern multi-billionbase DNA data sets. The difficulty is caused by the nonuniform (oligo)nucleotide composition of these sequences, rather than their size per se. To solve this problem, we modified the standard seed-and-extend approach (e.g., BLAST) to use adaptive seeds. Adaptive seeds are matches that are chosen based on their rareness, instead of using fixed-length matches. This method guarantees that the number of matches, and thus the running time, increases linearly, instead of quadratically, with sequence length. LAST, our open source implementation of adaptive seeds, enables fast and sensitive comparison of large sequences with arbitrarily nonuniform composition.

Identification of Y-box binding protein 1 as a core regulator of MEK/ERK pathway-dependent gene signatures in colorectal cancer cells.

Transcriptional signatures are an indispensible source of correlative information on disease-related molecular alterations on a genome-wide level. Numerous candidate genes involved in disease and in factors of predictive, as well as of prognostic, value have been deduced from such molecular portraits, e.g. in cancer. However, mechanistic insights into the regulatory principles governing global transcriptional changes are lagging behind extensive compilations of deregulated genes. To identify regulators of transcriptome alterations, we used an integrated approach combining transcriptional profiling of colorectal cancer cell lines treated with inhibitors targeting the receptor tyrosine kinase (RTK)/RAS/mitogen-activated protein kinase pathway, computational prediction of regulatory elements in promoters of co-regulated genes, chromatin-based and functional cellular assays. We identified commonly co-regulated, proliferation-associated target genes that respond to the MAPK pathway. We recognized E2F and NFY transcription factor binding sites as prevalent motifs in those pathway-responsive genes and confirmed the predicted regulatory role of Y-box binding protein 1 (YBX1) by reporter gene, gel shift, and chromatin immunoprecipitation assays. We also validated the MAPK-dependent gene signature in colorectal cancers and provided evidence for the association of YBX1 with poor prognosis in colorectal cancer patients. This suggests that MEK/ERK-dependent, YBX1-regulated target genes are involved in executing malignant properties.

Prostate-specific membrane antigen (PSMA) as a potential target for molecular imaging and treatment in bone and soft tissue sarcomas.

Bone and soft tissue sarcomas are a group of rare malignant tumours with major histological and anatomical varieties. In a metastatic setting, sarcomas have a poor prognosis due to limited response rates to chemotherapy. Radioligand therapy targeting prostate-specific membrane antigen (PSMA) may offer a new perspective. PSMA is a type II transmembrane glycoprotein which is present in all prostatic tissue and overexpressed in prostate cancer. Despite the name, PSMA is not prostate-specific. PSMA expression is also found in a multitude of non-prostatic diseases including a subgroup of sarcomas, mostly in its neovascular endothelial cells. On PET/CT imaging, multiple sarcomas have also shown intense PSMA-tracer accumulation. PSMA expression and PSMA-tracer uptake seem to be highest in patients with aggressive and advanced sarcomas, who are also in highest need of new therapeutic options. Although these results provide a good rationale for the future use of PSMA-targeted radioligand therapy in a selection of sarcoma patients, more research is needed to gain insight into optimal patient selection methods, PSMA-targeting antibodies and tracers, administered doses of radioligand therapy, and their efficacy and tolerability. In this review, mRNA expression of the FOLH1 gene which encodes PSMA, PSMA immunohistochemistry, PSMA-targeted imaging and PSMA-targeted therapy in sarcomas will be discussed.

Cancer-Associated Fibroblasts Are Key Determinants of Cancer Cell Invasion in the Earliest Stage of Colorectal Cancer.

BACKGROUND & AIMS: Improving clinical management of early stage colorectal cancers (T1CRCs) requires a better understanding of their underlying biology. Accumulating evidence shows that cancer-associated fibroblasts (CAFs) are important determinants of tumor progression in advanced colorectal cancer (CRC), but their role in the initial stages of CRC tumorigenesis is unknown. Therefore, we investigated the contribution of T1CAFs to early CRC progression. METHODS: Primary T1CAFs and patient-matched normal fibroblasts (NFs) were isolated from endoscopic biopsy specimens of histologically confirmed T1CRCs and normal mucosa, respectively. The impact of T1CAFs and NFs on tumor behavior was studied using 3-dimensional co-culture systems with primary T1CRC organoids and extracellular matrix (ECM) remodeling assays. Whole-transcriptome sequencing and gene silencing were used to pinpoint mediators of T1CAF functions. RESULTS: In 3-dimensional multicellular cultures, matrix invasion of T1CRC organoids was induced by T1CAFs, but not by matched NFs. Enhanced T1CRC invasion was accompanied by T1CAF-induced ECM remodeling and up-regulation of CD44 in epithelial cells. RNA sequencing of 10 NF-T1CAF pairs revealed 404 differentially expressed genes, with significant enrichment for ECM-related pathways in T1CAFs. Cathepsin H, a cysteine-type protease that was specifically up-regulated in T1CAFs but not in fibroblasts from premalignant lesions or advanced CRCs, was identified as a key factor driving matrix remodeling by T1CAFs. Finally, we showed high abundance of cathepsin H-expressing T1CAFs at the invasive front of primary T1CRC sections. CONCLUSIONS: Already in the earliest stage of CRC, cancer cell invasion is promoted by CAFs via direct interactions with epithelial cancer cells and stage-specific, cathepsin H-dependent ECM remodeling. RNA sequencing data of the 10 NF-T1CAF pairs can be found under GEO accession number GSE200660.

Inflammatory and tolerogenic myeloid cells determine outcome following human allergen challenge.

Innate mononuclear phagocytic system (MPS) cells preserve mucosal immune homeostasis. We investigated their role at nasal mucosa following allergen challenge with house dust mite. We combined single-cell proteome and transcriptome profiling on nasal immune cells from nasal biopsies cells from 30 allergic rhinitis and 27 non-allergic subjects before and after repeated nasal allergen challenge. Biopsies of patients showed infiltrating inflammatory HLA-DRhi/CD14+ and CD16+ monocytes and proallergic transcriptional changes in resident CD1C+/CD1A+ conventional dendritic cells (cDC)2 following challenge. In contrast, non-allergic individuals displayed distinct innate MPS responses to allergen challenge: predominant infiltration of myeloid-derived suppressor cells (MDSC: HLA-DRlow/CD14+ monocytes) and cDC2 expressing inhibitory/tolerogenic transcripts. These divergent patterns were confirmed in ex vivo stimulated MPS nasal biopsy cells. Thus, we identified not only MPS cell clusters involved in airway allergic inflammation but also highlight novel roles for non-inflammatory innate MPS responses by MDSC to allergens in non-allergic individuals. Future therapies should address MDSC activity as treatment for inflammatory airway diseases.

Methylation and deamination of CpGs generate p53-binding sites on a genomic scale.

The formation of transcription-factor-binding sites is an important evolutionary process. Here, we show that methylation and deamination of CpG dinucleotides generate in vivo p53-binding sites in numerous Alu elements and in non-repetitive DNA in a species-specific manner. In light of this, we propose that the deamination of methylated CpGs constitutes a universal mechanism for de novo generation of various transcription-factor-binding sites in Alus.

Targetfinder.org: a resource for systematic discovery of transcription factor target genes.

Targetfinder.org (http://targetfinder.org/) provides a web-based resource for finding genes that show a similar expression pattern to a group of user-selected genes. It is based on a large-scale gene expression compendium (>1200 experiments, >13,000 genes). The primary application of Targetfinder.org is to expand a list of known transcription factor targets by new candidate target genes. The user submits a group of genes (the 'seed'), and as a result the web site provides a list of other genes ranked by similarity of their expression to the expression of the seed genes. Additionally, the web site provides information on a recovery/cross-validation test to check for consistency of the provided seed and the quality of the ranking. Furthermore, the web site allows to analyse affinities of a selected transcription factor to the promoter regions of the top-ranked genes in order to select the best new candidate target genes for further experimental analysis.