Trans-driven variation in expression is common among detoxification genes in the extreme generalist herbivore Tetranychus urticae.

The extreme adaptation potential of the generalist herbivore Tetranychus urticae (the two-spotted spider mite) to pesticides as well as diverse host plants has been associated with clade-specific gene expansions in known detoxifying enzyme families, and with extensive and rapid transcriptional responses. However, how this broad transcriptional potential is regulated remains largely unknown. Using a parental/F1 design in which four inbred strains were crossed to a common inbred strain, we assessed the genetic basis and inheritance of gene expression variation in T. urticae. Mirroring known phenotypic variation in the progenitor strains of the inbreds, we confirmed that the inbred strains we created were genetically distinct, varied markedly in pesticide resistance, and also captured variation in host plant fitness as is commonly observed in this species. By examining differences in gene expression between parents and allele-specific expression in F1s, we found that variation in RNA abundance was more often explained in trans as compared to cis, with the former associated with dominance in inheritance. Strikingly, in a gene ontology analysis, detoxification genes of the cytochrome P450 monooxygenase (CYP) family, as well as dioxygenases (DOGs) acquired from horizontal gene transfer from fungi, were specifically enriched at the extremes of trans-driven up- and downregulation. In particular, multiple CYPs and DOGs with broad substrate-specificities for pesticides or plant specialized compounds were exceptionally highly upregulated as a result of trans-regulatory variation, or in some cases synergism of cis and trans, in the most multi-pesticide resistant strains. Collectively, our findings highlight the potential importance of trans-driven expression variation in genes associated with xenobiotic metabolism and host plant use for rapid adaptation in T. urticae, and also suggests modular control of these genes, a regulatory architecture that might ameliorate negative pleiotropic effects.

Dicer-like 3a mediates intergenerational resistance against root-knot nematodes in rice via hormone responses.

In a continuously changing and challenging environment, passing down the memory of encountered stress factors to offspring could provide an evolutionary advantage. In this study, we demonstrate the existence of "intergenerational acquired resistance" in the progeny of rice (Oryza sativa) plants attacked by the belowground parasitic nematode Meloidogyne graminicola. Transcriptome analyses revealed that genes involved in defense pathways are generally downregulated in progeny of nematode-infected plants under uninfected conditions but show a stronger induction upon nematode infection. This phenomenon was termed "spring loading" and depends on initial downregulation by the 24-nucleotide (nt) siRNA biogenesis gene dicer-like 3a (dcl3a) involved in the RNA-directed DNA methylation pathway. Knockdown of dcl3a led to increased nematode susceptibility and abolished intergenerational acquired resistance, as well as jasmonic acid/ethylene spring loading in the offspring of infected plants. The importance of ethylene signaling in intergenerational resistance was confirmed by experiments on a knockdown line of ethylene insensitive 2 (ein2b), which lacks intergenerational acquired resistance. Taken together, these data indicate a role for DCL3a in regulating plant defense pathways during both within-generation and intergenerational resistance against nematodes in rice.

Liver transcriptomic and methylomic analyses identify transcriptional mitogen-activated protein kinase regulation in facultative hibernation of Syrian hamster.

Hibernation consists of alternating torpor-arousal phases, during which animals cope with repetitive hypothermia and ischaemia-reperfusion. Due to limited transcriptomic and methylomic information for facultative hibernators, we here conducted RNA and whole-genome bisulfide sequencing in liver of hibernating Syrian hamster (Mesocricetus auratus). Gene ontology analysis was performed on 844 differentially expressed genes and confirmed the shift in metabolic fuel utilization, inhibition of RNA transcription and cell cycle regulation as found in seasonal hibernators. Additionally, we showed a so far unreported suppression of mitogen-activated protein kinase (MAPK) and protein phosphatase 1 pathways during torpor. Notably, hibernating hamsters showed upregulation of MAPK inhibitors (dual-specificity phosphatases and sproutys) and reduced levels of MAPK-induced transcription factors (TFs). Promoter methylation was found to modulate the expression of genes targeted by these TFs. In conclusion, we document gene regulation between hibernation phases, which may aid the identification of pathways and targets to prevent organ damage in transplantation or ischaemia-reperfusion.

Concurrent invasions of European starlings in Australia and North America reveal population-specific differentiation in shared genomic regions.

A species' success during the invasion of new areas hinges on an interplay between the demographic processes common to invasions and the specific ecological context of the novel environment. Evolutionary genetic studies of invasive species can investigate how genetic bottlenecks and ecological conditions shape genetic variation in invasions, and our study pairs two invasive populations that are hypothesized to be from the same source population to compare how each population evolved during and after introduction. Invasive European starlings (Sturnus vulgaris) established populations in both Australia and North America in the 19th century. Here, we compare whole-genome sequences among native and independently introduced European starling populations to determine how demographic processes interact with rapid evolution to generate similar genetic patterns in these recent and replicated invasions. Demographic models indicate that both invasive populations experienced genetic bottlenecks as expected based on invasion history, and we find that specific genomic regions have differentiated even on this short evolutionary timescale. Despite genetic bottlenecks, we suggest that genetic drift alone cannot explain differentiation in at least two of these regions. The demographic boom intrinsic to many invasions as well as potential inversions may have led to high population-specific differentiation, although the patterns of genetic variation are also consistent with the hypothesis that this infamous and highly mobile invader adapted to novel selection (e.g., extrinsic factors). We use targeted sampling of replicated invasions to identify and evaluate support for multiple, interacting evolutionary mechanisms that lead to differentiation during the invasion process.

Spatiotemporal expression profile of novel and known small RNAs throughout rice plant development focussing on seed tissues.

BACKGROUND: Small RNAs (sRNAs) regulate numerous plant processes directly related to yield, such as disease resistance and plant growth. To exploit this yield-regulating potential of sRNAs, the sRNA profile of one of the world's most important staple crops - rice - was investigated throughout plant development using next-generation sequencing. RESULTS: Root and leaves were investigated at both the vegetative and generative phase, and early-life sRNA expression was characterized in the embryo and endosperm. This led to the identification of 49,505 novel sRNAs and 5581 tRNA-derived sRNAs (tsRNAs). In all tissues, 24 nt small interfering RNAs (siRNAs) were highly expressed and associated with euchromatic, but not heterochromatic transposable elements. Twenty-one nt siRNAs deriving from genic regions in the endosperm were exceptionally highly expressed, mimicking previously reported expression levels of 24 nt siRNAs in younger endosperm samples. In rice embryos, sRNA content was highly diverse while tsRNAs were underrepresented, possibly due to snoRNA activity. Publicly available mRNA expression and DNA methylation profiles were used to identify putative siRNA targets in embryo and endosperm. These include multiple genes related to the plant hormones gibberellic acid and ethylene, and to seed phytoalexin and iron content. CONCLUSIONS: This work introduces multiple sRNAs as potential regulators of rice yield and quality, identifying them as possible targets for the continuous search to optimize rice production.

Non-coding RNAs in the interaction between rice and Meloidogyne graminicola.

BACKGROUND: Root knot nematodes (RKN) are plant parasitic nematodes causing major yield losses of widely consumed food crops such as rice (Oryza sativa). Because non-coding RNAs, including small interfering RNAs (siRNA), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are key regulators of various plant processes, elucidating their regulation during this interaction may lead to new strategies to improve crop protection. In this study, we aimed to identify and characterize rice siRNAs, miRNAs and lncRNAs responsive to early infection with RKN Meloidogyne graminicola (Mg), based on sequencing of small RNA, degradome and total RNA libraries from rice gall tissues compared with uninfected root tissues. RESULTS: We found 425 lncRNAs, 3739 siRNAs and 16 miRNAs to be differentially expressed between both tissues, of which a subset was independently validated with RT-qPCR. Functional prediction of the lncRNAs indicates that a large part of their potential target genes code for serine/threonine protein kinases and transcription factors. Differentially expressed siRNAs have a predominant size of 24 nts, suggesting a role in DNA methylation. Differentially expressed miRNAs are generally downregulated and target transcription factors, which show reduced degradation according to the degradome data. CONCLUSIONS: To our knowledge, this work is the first to focus on small and long non-coding RNAs in the interaction between rice and Mg, and provides an overview of rice non-coding RNAs with the potential to be used as a resource for the development of new crop protection strategies.

The genome of the extremophile Artemia provides insight into strategies to cope with extreme environments.

BACKGROUND: Brine shrimp Artemia have an unequalled ability to endure extreme salinity and complete anoxia. This study aims to elucidate its strategies to cope with these stressors. RESULTS AND DISCUSSION: Here, we present the genome of an inbred A. franciscana Kellogg, 1906. We identified 21,828 genes of which, under high salinity, 674 genes and under anoxia, 900 genes were differentially expressed (42%, respectively 30% were annotated). Under high salinity, relevant stress genes and pathways included several Heat Shock Protein and Leaf Embryogenesis Abundant genes, as well as the trehalose metabolism. In addition, based on differential gene expression analysis, it can be hypothesized that a high oxidative stress response and endocytosis/exocytosis are potential salt management strategies, in addition to the expression of major facilitator superfamily genes responsible for transmembrane ion transport. Under anoxia, genes involved in mitochondrial function, mTOR signalling and autophagy were differentially expressed. Both high salt and anoxia enhanced degradation of erroneous proteins and protein chaperoning. Compared with other branchiopod genomes, Artemia had 0.03% contracted and 6% expanded orthogroups, in which 14% of the genes were differentially expressed under high salinity or anoxia. One phospholipase D gene family, shown to be important in plant stress response, was uniquely present in both extremophiles Artemia and the tardigrade Hypsibius dujardini, yet not differentially expressed under the described experimental conditions. CONCLUSIONS: A relatively complete genome of Artemia was assembled, annotated and analysed, facilitating research on its extremophile features, and providing a reference sequence for crustacean research.

Molecular correlates of hypothalamic development in songbird ontogeny in comparison with the telencephalon.

Development of the songbird brain provides an excellent experimental model for understanding the regulation of sex differences in ontogeny. Considering the regulatory role of the hypothalamus in endocrine, in particular reproductive, physiology, we measured the structural (volume) and molecular correlates of hypothalamic development during ontogeny of male and female zebra finches. We quantified by relative quantitative polymerase chain reaction (rqPCR) the expression of 14 genes related to thyroid and steroid hormones actions as well as 12 genes related to brain plasticity at four specific time points during ontogeny and compared these expression patterns with the expression of the same genes as detected by transcriptomics in the telencephalon. These two different methodological approaches detected specific changes with age and demonstrated that in a substantial number of cases changes observed in both brain regions are nearly identical. Other genes however had a tissue-specific developmental pattern. Sex differences or interactions of sex by age were detected in the expression of a subset of genes, more in hypothalamus than telencephalon. These results correlate with multiple known aspects of the developmental and reproductive physiology but also raise a number of new functional questions.

MEXPRESS update 2019.

The recent growth in the number of publicly available cancer omics databases has been accompanied by the development of various tools that allow researchers to visually explore these data. In 2015, we built MEXPRESS, an online tool for the integration and visualization of gene expression, DNA methylation and clinical data from The Cancer Genome Atlas (TCGA), a large collection of publicly available multi-omics cancer data. MEXPRESS addresses the need for an easy-to-use, interactive application that allows researchers to identify dysregulated genes and their clinical relevance in cancer. Furthermore, while other tools typically do not support integrated visualization of expression and DNA methylation data in combination with the precise genomic location of the methylation, MEXPRESS is unique in how it depicts these diverse data types together. Motivated by the large number of users MEXPRESS has managed to attract over the past 3 years and the recent migration of all TCGA data to a new data portal, we developed a new version of MEXPRESS (https://mexpress.be). It contains the latest TCGA data, additional types of omics and clinical data and extra functionality, allowing users to explore mechanisms of gene dysregulation beyond expression and DNA methylation.

Molecular insights into the compatible and incompatible interactions between sugar beet and the beet cyst nematode.

BACKGROUND: Sugar beet (Beta vulgaris subsp. vulgaris) is an economically important crop that provides nearly one third of the global sugar production. The beet cyst nematode (BCN), Heterodera schachtii, causes major yield losses in sugar beet and other crops worldwide. The most effective and economic approach to control this nematode is growing tolerant or resistant cultivars. To identify candidate genes involved in susceptibility and resistance, the transcriptome of sugar beet and BCN in compatible and incompatible interactions at two time points was studied using mRNA-seq. RESULTS: In the susceptible cultivar, most defense-related genes were induced at 4 dai while suppressed at 10 dai but in the resistant cultivar Nemakill, induction of genes involved in the plant defense response was observed at both time points. In the compatible interaction, alterations in phytohormone-related genes were detected. The effect of exogenous application of Methyl Jasmonate and ET-generator ethephon on susceptible plants was therefore investigated and the results revealed significant reduction in plant susceptibility. Genes putatively involved in the resistance of Nemakill were identified, such as genes involved in phenylpropanoid pathway and genes encoding CYSTM domain-containing proteins, F-box proteins, chitinase, galactono-1,4-lactone dehydrogenase and CASP-like protein. Also, the transcriptome of the BCN was analyzed in infected root samples and several novel potential nematode effector genes were found. CONCLUSIONS: Our data provides detailed insights into the plant and nematode transcriptional changes occurring during compatible and incompatible interactions between sugar beet and BCN. Many important genes playing potential roles in susceptibility or resistance of sugar beet against BCN, as well as some BCN effectors with a potential role as avr proteins were identified. In addition, our findings indicate the effective role of jasmonate and ethylene in enhancing sugar beet defense response against BCN. This research provides new molecular insights into the plant-nematode interactions that can be used to design novel management strategies against BCN.

Genome-wide DNA hypomethylation shapes nematode pattern-triggered immunity in plants.

A role for DNA hypomethylation has recently been suggested in the interaction between bacteria and plants; it is unclear whether this phenomenon reflects a conserved response. Treatment of plants of monocot rice and dicot tomato with nematode-associated molecular patterns from different nematode species or bacterial pathogen-associated molecular pattern flg22 revealed global DNA hypomethylation. A similar hypomethylation response was observed during early gall induction by Meloidogyne graminicola in rice. Evidence for the causal impact of hypomethylation on immunity was revealed by a significantly reduced plant susceptibility upon treatment with DNA methylation inhibitor 5-azacytidine. Whole-genome bisulphite sequencing of young galls revealed massive hypomethylation in the CHH context, while not for CG or CHG nucleotide contexts. Further, CHH hypomethylated regions were predominantly associated with gene promoter regions, which was not correlated with activated gene expression at the same time point but, rather, was correlated with a delayed transcriptional gene activation. Finally, the relevance of CHH hypomethylation in plant defence was confirmed in rice mutants of the RNA-directed DNA methylation pathway and DECREASED DNA METHYLATION 1. We demonstrated that DNA hypomethylation is associated with reduced susceptibility in rice towards root-parasitic nematodes and is likely to be part of the basal pattern-triggered immunity response in plants.

Targeted RNA-seq successfully identifies normal and pathogenic splicing events in breast/ovarian cancer susceptibility and Lynch syndrome genes.

A subset of genetic variants found through screening of patients with hereditary breast and ovarian cancer syndrome (HBOC) and Lynch syndrome impact RNA splicing. Through target enrichment of the transcriptome, it is possible to perform deep-sequencing and to identify the different and even rare mRNA isoforms. A targeted RNA-seq approach was used to analyse the naturally-occurring splicing events for a panel of 8 breast and/or ovarian cancer susceptibility genes (BRCA1, BRCA2, RAD51C, RAD51D, PTEN, STK11, CDH1, TP53), 3 Lynch syndrome genes (MLH1, MSH2, MSH6) and the fanconi anaemia SLX4 gene, in which monoallelic mutations were found in non-BRCA families. For BRCA1, BRCA2, RAD51C and RAD51D the results were validated by capillary electrophoresis and were compared to a non-targeted RNA-seq approach. We also compared splicing events from lymphoblastoid cell-lines with those from breast and ovarian fimbriae tissues. The potential of targeted RNA-seq to detect pathogenic changes in RNA-splicing was validated by the inclusion of samples with previously well characterized BRCA1/2 genetic variants. In our study, we update the catalogue of normal splicing events for BRCA1/2, provide an extensive catalogue of normal RAD51C and RAD51D alternative splicing, and list splicing events found for eight other genes. Additionally, we show that our approach allowed the identification of aberrant splicing events due to the presence of BRCA1/2 genetic variants and distinguished between complete and partial splicing events. In conclusion, targeted-RNA-seq can be very useful to classify variants based on their putative pathogenic impact on splicing.

Selection of miRNA reference genes for plant defence studies in rice (Oryza sativa).

MAIN CONCLUSION: MicroRNAs miR390-5p, miR7694-3p miR1868 and miR1849 were found to be suitable miRNA reference genes for rice, under either infection with the root-knot nematode Meloidogyne graminicola or treatment with BABA. RT-qPCR is a widely used method to investigate the expression levels of genes under certain conditions. A key step, however, to have reliable results is the normalization of expression. For every experimental condition, suitable reference genes must be chosen. These reference genes must not be affected by differences in experimental conditions. MicroRNAs are regulatory RNA molecules, able to direct the expression levels of protein coding genes. In plants, their attributed functions range from roles in development to immunity. In this work, microRNAs (miRNAs) are evaluated for their suitability as reference genes in rice after infection with root-knot nematode Meloidogyne graminicola or after priming with beta-amino butyric acid. The evaluation was based on their amplification efficiency and their stability estimates according to geNorm, NormFinder and BestKeeper. All tested miRNAs, excluding one, were considered acceptable for normalization. Furthermore, miRNAs were validated using miRNA sequencing data. The set of microRNAs miR390-5p and miR7694-3p was found to be the most stable combination under the tested conditions. Another miRNA set consisting of miR7694-3p, miR1868 and miR1849 also shows potential to be used for miRNA expression normalization under experimental conditions beyond the scope of this study. This work is the first report on reference miRNAs in rice for the purpose of plant defence studies.

Below-Ground Attack by the Root Knot Nematode Meloidogyne graminicola Predisposes Rice to Blast Disease.

Magnaporthe oryzae (rice blast) and the root-knot nematode Meloidogyne graminicola are causing two of the most important pathogenic diseases jeopardizing rice production. Here, we show that root-knot nematode infestation on rice roots leads to important above-ground changes in plant immunity gene expression, which is correlated with significantly enhanced susceptibility to blast disease. A detailed metabolic analysis of oxidative stress responses and hormonal balances demonstrates that the above-ground tissues have a disturbed oxidative stress level, with accumulation of H2O2, as well as hormonal disturbances. Moreover, double infection experiments on an oxidative stress mutant and an auxin-deficient rice line indicate that the accumulation of auxin in the above-ground tissue is at least partly responsible for the blast-promoting effect of root-knot nematode infection.

Methylome analysis of extreme chemoresponsive patients identifies novel markers of platinum sensitivity in high-grade serous ovarian cancer.

BACKGROUND: Despite an early response to platinum-based chemotherapy in advanced stage high-grade serous ovarian cancer (HGSOC), the majority of patients will relapse with drug-resistant disease. Aberrant epigenetic alterations like DNA methylation are common in HGSOC. Differences in DNA methylation are associated with chemoresponse in these patients. The objective of this study was to identify and validate novel epigenetic markers of chemoresponse using genome-wide analysis of DNA methylation in extreme chemoresponsive HGSOC patients. METHODS: Genome-wide next-generation sequencing was performed on methylation-enriched tumor DNA of two HGSOC patient groups with residual disease, extreme responders (≥18 months progression-free survival (PFS), n = 8) and non-responders (≤6 months PFS, n = 10) to platinum-based chemotherapy. DNA methylation and expression data of the same patients were integrated to create a gene list. Genes were validated on an independent cohort of extreme responders (n = 21) and non-responders (n = 31) using pyrosequencing and qRT-PCR. In silico validation was performed using publicly available DNA methylation (n = 91) and expression (n = 208) datasets of unselected advanced stage HGSOC patients. Functional validation of FZD10 on chemosensitivity was carried out in ovarian cancer cell lines using siRNA-mediated silencing. RESULTS: Integrated genome-wide methylome and expression analysis identified 45 significantly differentially methylated and expressed genes between two chemoresponse groups. Four genes FZD10, FAM83A, MYO18B, and MKX were successfully validated in an external set of extreme chemoresponsive HGSOC patients. High FZD10 and MKX methylation were related with extreme responders and high FAM83A and MYO18B methylation with non-responders. In publicly available advanced stage HGSOC datasets, FZD10 and MKX methylation levels were associated with PFS. High FZD10 methylation was strongly associated with improved PFS in univariate analysis (hazard ratio (HR) = 0.43; 95% CI, 0.27-0.71; P = 0.001) and multivariate analysis (HR = 0.39; 95% CI, 0.23-0.65; P = 0.003). Consistently, low FZD10 expression was associated with improved PFS (HR = 1.36; 95% CI, 0.99-1.88; P = 0.058). FZD10 silencing caused significant sensitization towards cisplatin treatment in survival assays and apoptosis assays. CONCLUSIONS: By applying genome-wide integrated methylome analysis on extreme chemoresponsive HGSOC patients, we identified novel clinically relevant, epigenetically-regulated markers of platinum-sensitivity in HGSOC patients. The clinical potential of these markers in predictive and therapeutic approaches has to be further validated in prospective studies.

PROTEOFORMER: deep proteome coverage through ribosome profiling and MS integration.

An increasing amount of studies integrate mRNA sequencing data into MS-based proteomics to complement the translation product search space. However, several factors, including extensive regulation of mRNA translation and the need for three- or six-frame-translation, impede the use of mRNA-seq data for the construction of a protein sequence search database. With that in mind, we developed the PROTEOFORMER tool that automatically processes data of the recently developed ribosome profiling method (sequencing of ribosome-protected mRNA fragments), resulting in genome-wide visualization of ribosome occupancy. Our tool also includes a translation initiation site calling algorithm allowing the delineation of the open reading frames (ORFs) of all translation products. A complete protein synthesis-based sequence database can thus be compiled for mass spectrometry-based identification. This approach increases the overall protein identification rates with 3% and 11% (improved and new identifications) for human and mouse, respectively, and enables proteome-wide detection of 5'-extended proteoforms, upstream ORF translation and near-cognate translation start sites. The PROTEOFORMER tool is available as a stand-alone pipeline and has been implemented in the galaxy framework for ease of use.

Identification and validation of WISP1 as an epigenetic regulator of metastasis in oral squamous cell carcinoma.

Lymph node (LN) metastasis is the most important prognostic factor in oral squamous cell carcinoma (OSCC) patients. However, in approximately one third of OSCC patients nodal metastases remain undetected, and thus are not adequately treated. Therefore, clinical assessment of LN metastasis needs to be improved. The purpose of this study was to identify DNA methylation biomarkers to predict LN metastases in OSCC. Genome wide methylation assessment was performed on six OSCC with (N+) and six without LN metastases (N0). Differentially methylated sequences were selected based on the likelihood of differential methylation and validated using an independent OSCC cohort as well as OSCC from The Cancer Genome Atlas (TCGA). Expression of WISP1 using immunohistochemistry was analyzed on a large OSCC cohort (n = 204). MethylCap-Seq analysis revealed 268 differentially methylated markers. WISP1 was the highest ranking annotated gene that showed hypomethylation in the N+ group. Bisulfite pyrosequencing confirmed significant hypomethylation within the WISP1 promoter region in N+ OSCC (P = 0.03) and showed an association between WISP1 hypomethylation and high WISP1 expression (P = 0.01). Both these results were confirmed using 148 OSCC retrieved from the TCGA database. In a large OSCC cohort, high WISP1 expression was associated with LN metastasis (P = 0.05), disease-specific survival (P = 0.022), and regional disease-free survival (P = 0.027). These data suggest that WISP1 expression is regulated by methylation and WISP1 hypomethylation contributes to LN metastasis in OSCC. WISP1 is a potential biomarker to predict the presence of LN metastases.

Molecular and epigenetic features of melanomas and tumor immune microenvironment linked to durable remission to ipilimumab-based immunotherapy in metastatic patients.

BACKGROUND: Ipilimumab (Ipi) improves the survival of advanced melanoma patients with an incremental long-term benefit in 10-15 % of patients. A tumor signature that correlates with this survival benefit could help optimizing individualized treatment strategies. METHODS: Freshly frozen melanoma metastases were collected from patients treated with either Ipi alone (n: 7) or Ipi combined with a dendritic cell vaccine (TriMixDC-MEL) (n: 11). Samples were profiled by immunohistochemistry (IHC), whole transcriptome (RNA-seq) and methyl-DNA sequencing (MBD-seq). RESULTS: Patients were divided in two groups according to clinical evolution: durable benefit (DB; 5 patients) and no clinical benefit (NB; 13 patients). 20 metastases were profiled by IHC and 12 were profiled by RNA- and MBD-seq. 325 genes were identified as differentially expressed between DB and NB. Many of these genes reflected a humoral and cellular immune response. MBD-seq revealed differences between DB and NB patients in the methylation of genes linked to nervous system development and neuron differentiation. DB tumors were more infiltrated by CD8(+) and PD-L1(+) cells than NB tumors. B cells (CD20(+)) and macrophages (CD163(+)) co-localized with T cells. Focal loss of HLA class I and TAP-1 expression was observed in several NB samples. CONCLUSION: Combined analyses of melanoma metastases with IHC, gene expression and methylation profiling can potentially identify durable responders to Ipi-based immunotherapy.