ALBERT

Description: This paper develops an optimizing model for the long-term exploitation of limited natural gas reserves in China. In addition to describing the life cycle characteristics of natural gas production and introducing the inter-temporal allocation theory, this paper builds the optimal exploitation model of natural gas resources within a gas field in the Ordos Basin as an example to analyze its exploitation scale and how influence factors, such as recovery rate, discount rate and the gas well exhausting cycle, affect the optimal exploration path of this gas field. We determine that an increase in the discount rate stimulates investors to invest more aggressively in natural gas exploitation in the early period due to the lower discounted value, thereby increasing the pace of the exploitation of natural gas and the exhaustion of gas fields. A higher recoverable factor implies more recoverable reserves and greater potential of increasing the output of gas fields. The exhaustion rate of gas wells affects the capability of converting capacity to output. When exhaustion occurs quickly in gas wells, the output will likely increase in the output rising period, and the output will likely decrease at a faster rate in the output reduction period. Price reform affects the economic recoverable reserves of gas fields.

Description: Brain–computer interfaces (BCI) have witnessed a rapid development in recent years. However, the active BCI paradigm is still underdeveloped with a lack of variety. It is imperative to adapt more voluntary mental activities for the active BCI control, which can induce separable electroencephalography (EEG) features. This study aims to demonstrate the brain function of timing prediction, i.e., the expectation of upcoming time intervals, is accessible for BCIs. Eighteen subjects were selected for this study. They were trained to have a precise idea of two sub-second time intervals, i.e., 400 ms and 600 ms, and were asked to measure a time interval of either 400 ms or 600 ms in mind after a cue onset. The EEG features induced by timing prediction were analyzed and classified using the combined discriminative canonical pattern matching and common spatial pattern. It was found that the ERPs in low-frequency (0~4 Hz) and energy in high-frequency (20~60 Hz) were separable for distinct timing predictions. The accuracy reached the highest of 93.75% with an average of 76.45% for the classification of 400 vs. 600 ms timing. This study first demonstrates that the cognitive EEG features induced by timing prediction are detectable and separable, which is feasible to be used in active BCIs controls and can broaden the category of BCIs.

Description: BACKGROUND: The biological basis for the adverse prognosis of chr1q gain/amplification (1q+) present in ~30% of patients with multiple myeloma (MM) remains ill-defined. The transcription factor (TF) Pre-B-cell leukaemia homeobox 1 (PBX1), encoded on chr1q, acts as master regulator of early hematopoiesis and as an oncogene in leukemia and other malignancies. Herein, we hypothesized that PBX1 orchestrates proliferative regulatory networks that underpin the poor prognosis associated with 1q+ in MM. METHODS: We employed qPCR for mRNA quantification, western blotting and immunohistochemistry for protein analysis, lentiviral shRNA-mediated knock-down, Hoechst/Annexin V staining and flow-cytometry for apoptosis and cell cycle analysis, ChIPseq for cistrome and RNAseq after knock-down for transcriptome analysis. Additional data were obtained from MMRF/CoMMpass, Blueprint Consortium and Arkansas datasets. Computational analysis of clinical and "-omics" data was performed using standard bioinformatic work-flows; pathway enrichment analysis using EnrichR and GSEA. RESULTS: Combined genomic (WGS/WES) and transcriptomic (RNAseq) analysis of the CoMMpass dataset identified a subgroup of 1q+ MM patients (60%) characterized by aberrant PBX1 overexpression and amplification (PBX1amp); survival analysis revealed significantly worse outcome of this subgroup compared to 1q+ non-PBX1amp (p

Description: Overexpression of the transcription factor MAF, as a result of its juxtaposition to the IgH enhancer [MAF-translocated t(14;16)], is a myeloma-initiating event in 3-5% of patients with multiple myeloma (MM) and confers a poor prognosis. MAF is also overexpressed in another 40% of cases, often in co-operation with the oncogene MMSET. The mechanisms by which MAF overexpression impacts on the regulatory genome to generate the MAF-driven oncogenic transcriptome and its direct targets are not known. To address this, we employed a multi-layer -omics approach using primary myeloma plasma cells (PC) as well as myeloma cell lines (MMCL). First, we determined the chromatin accessibility and transcriptome profiles of MAF-translocated myeloma by performing ATAC-seq and RNA-seq, respectively, in purified bone marrow CD138+ PC from two patients with t(14;16) and three healthy donors. We identified 6,640 differentially accessible regions, 87% of which displayed enhanced chromatin accessibility in MAF samples compared to normal PC. Secondary analysis comparing this with ATAC-seq data from a set of 28 other MM samples, including hyperdiploid, MMSET and CCND1-translocated MM, revealed 33% of those regions to be MAF subgroup specific (1,949 regions), with the rest shared between MAF and other cytogenetic groups. Gene annotation and pathway enrichment analysis using GREAT confirmed overrepresentation of the MF myeloma patient signature, as previously identified in microarray datasets. RNA-seq detected significant upregulation of approximately 900 genes in MAF samples compared to normal counterparts, including MAF itself (top 4th hit) as well as its presumed targets (CCND2, ITGB7 and NUAK1). Next, we obtained the MAF cistrome using ChIP-seq in the MAF-translocated MMCL MM1.S and integrated it with the primary PC ATAC-seq data. This revealed that 31% (618/1,949) of the differentially accessible regions in MAF-translocated MM PC are also MAF-bound. Additional overlay with ENCODE ChromHMM epigenome map showed that 47% of MAF binding sites are on active enhancers and 42% on active promoters signifying potential direct regulation of the corresponding genes. Next, we superimposed the accessible and MAF-bound loci on the epigenomic landscapes of normal PC and other B-cell types using their corresponding ChromHMM maps (Blueprint consortium data). Interestingly, 56% (345/618) of the MAF-specific regions were not active in any stage of B cell development. This suggests that aberrant MAF overexpression and chromatin binding in PC is associated with de novo activation of these chromatin regions, over half of which (200/345; 58%) are enhancers; we termed these 'neo-enhancers'. Upon de novo motif analysis of MAF ChIP-seq in MAF-translocated JJN3 and MM1.S MMCL, we confirmed MAF as the first and, interestingly, IRF4 as the second top hit, suggesting a possible MAF-IRF4 functional interaction in myelomagenesis. Indeed, overlay of the accessible MAF-bound loci with IRF4 ChIP-seq data in MM1.S revealed 63% co-occupancy (including 62% of "neo-enhancers"), proposing a novel and extensive co-operative chromatin-based network between the two transcription factors. Final integration of the accessible MAF-bound regions with the paired transcriptomes of primary myeloma PC revealed a set 206 candidate enhancer-gene pairs. Strikingly, we identified two IRF4-cobound "neo-enhancers" linked to overexpression of TLR4 and CCR1, two genes known for their roles in myeloma cell proliferation and migration. We confirmed significant downregulation of both genes upon shRNA-mediated knockdown of MAF in the two MAF-translocated MMCL, MM1.S and JJN3, as well as the lethality of MAF depletion. Further, MAF overexpression in MAF-negative myeloma backgrounds led to transcriptional upregulation of these genes, further validating them as MAF targets. While CRISPR/Cas9i experiments targeting TLR4 are ongoing, preliminary results validated the functional role of the "neo-enhancer" in CCR1 gene expression. In conclusion, we demonstrate for the first time an extensive re-organisation of the PC chromatin conferred by oncogenic MAF in MM; we reveal its extensive co-operation with IRF4 in this process; we validate the directly MAF-regulated genes and functionally characterise neo-enhancers of key MAF-dependent genes that in addition to MAF itself are also critical for myeloma biology. Disclosures Hatjiharissi: Janssen: Honoraria. Caputo:GSK: Research Funding. Karadimitris:GSK: Research Funding.

Description: In multiple myeloma (MM), a malignancy of the bone marrow plasma cells (BMPC), hyperdiploidy (HY) and oncogene over-expression via chromosomal translocation [including CCND1- t(11;14), MAF- t(14;16), MMSET-t(4;14)] are the primary myeloma initiating events (MIE) that drive distinct transcriptional programs. These are further shaped by secondary SNV and CNV events. This genetic heterogeneity converges, in most cases, to a functionally dichotomous state of CCND1 or CCND2 overexpression. The molecular mechanisms underlying each of the distinct myelomagenic transcriptomes and the CCND1 vs CCND2 dichotomy have not been defined. To address these questions, we obtained highly purified BMPC from 3 healthy donors and 30 MM patients (HY: 15; CCND1: 4; MMSET: 5; MAF: 2; other: 4), either at diagnosis or relapse, and mapped their chromatin accessibility and transcriptome profiles by ATAC-seq and RNA-seq, respectively. In total, we obtained ~300K regions with accessible chromatin in either MM or normal PC. Overall chromatin accessibility increased in myeloma compared to normal PC, particularly in MAF- and MMSET-translocated subtypes. Analysis of combined ATAC-seq/RNA-seq by Multi-Omics Factor Analysis (MOFA) resulted in a clearer samples distinction than either ATAC-seq or RNA-seq alone, with altered chromatin accessibility accounting for more of the variance than expression. Of the top five identified factors, the top two (one transcriptome driven, one accessibility driven) distinguished normal from MM samples, whilst two more separated MMSET, MAF and CCND1 subgroups. Ninety seven, 157, 256 and 348 overexpressed genes in the CCND1, HY, MMSET and MAF subgroups, respectively, were predicted to be regulated by differentially accessible enhancers. Twenty percent (165/858) of these genes were overexpressed in 〉1 subgroup suggesting a process of chromatin accessibility-based convergence evolution. Enrichment analysis suggested direct or indirect involvement of Polycomb and chromatin remodellers; significant enrichment was also found for genes involved in neurogenesis. ATAC-seq footprinting predicted binding sites for 250 expressed transcription factors (TFs), 116 of which displayed higher binding frequency in myeloma than in normal PC and included both known (e.g., XBP1, RELA, IRF4, PRDM1) and potentially novel regulators of myeloma biology (e.g., CXXC1 and NFE2L1). The remaining 134 TF were predicted to be present in at least one MM subgroup, but absent in normal PC. Amongst them, as expected, MAF was active in the MMSET- and more so in the MAF-translocated subgroups. DepMap database analysis suggested myeloma cell dependency on 181/250 TF (CRISPR/Cas9 CERES score 〈 -0.1 in 〉3/14 MMCL analysed). In dissecting the regulatory basis of CCND2 vs CCND1 dichotomy, one MOFA factor completely separated MAF from CCND1 samples, placing extreme opposite weights on the expression of CCND2 and CCND1 respectively. Interestingly, the same factor identified open-chromatin clusters upstream of CCND2 and linked them to its over-expression. These clusters were also open in the MMSET group and in CCND2-expressing HY samples. Conversely, no accessibility was detected in the CCND1 group, the CCND1-expressing HY samples or in normal PCs. Further, super-enhancer calling using the H3K27ac histone mark in MAF-translocated JJN3 cells identified the region of interest as a bona fide super-enhancer. Chromatin long range interactions, as assessed by Capture-HiC, demonstrated high frequency interactions of the CCND2 promoter with the constituent elements of the putative super-enhancer. Experimental validation using a CRISPR/Cas9i system confirmed the functional role of all 4 super-enhancer constituents tested in the regulation of CCND2 expression, while TF footprinting predicted MAF binding to the super-enhancer in MAF-translocated PC. In conclusion, we show that distinct oncogenic transcriptomes in MM are underpinned by extensive chromatin changes, accompanied by TF activity 're-wiring' that does not necessarily require transcriptional deregulation of the TF themselves. We identify novel, non-oncogene TF dependencies that suggest therapeutic opportunities in MM and we discover and characterise the critical super-enhancer that drives overexpression of the CCND2 oncogene in MM. Disclosures Auner: Amgen: Other: Consultancy and Research Funding; Takeda: Consultancy; Karyopharm: Consultancy. Hatjiharissi:Janssen: Honoraria. Caputo:GSK: Research Funding. Karadimitris:GSK: Research Funding.

Description: BACKGROUND: Bone disease, a common source of morbidity in multiple myeloma (MM), is caused by RANKL-induced aberrant activation of osteoclasts (OC). RANKL-induced OC lineage commitment requires repression of an Irf-8 dependent macrophage inflammatory transcriptional programme commensurate with activation of an OC lineage-specific programme. Functional data have shown the requirement for the histone acetylation readers Brd2-4 BET proteins and of cMyc for OC lineage development. However, how Brd2-4 and Myc co-operate genome-wide to regulate transcriptome changes that underpin the very early stages of RANKL-induced OC lineage commitment has not been defined. METHODS: The OC progenitor-like murine RAW264.7 cell line was used for osteoclastogenesis. OC were assayed by TRAP staining. We performed RNA-seq for transcriptome analysis and ChIP-seq against Brd2-4, cMyc, and H3K27Ac mark for epigenomic profiling. The pan-Bet inhibitor IBET151 was used alone or in combination with RANKL. ChIP-seq/RNA-seq data were processed using standard bioinformatics pipelines; downstream analyses (pathway and motif enrichment, factor differential binding) were performed by various tools including EnrichR, R packages ChIPpeakAnno/DiffBind, Rose. RESULTS: Transcriptomic profiling of OC progenitors at 0, 4, 14 and 24h post-RANKL treatment identified 12 distinct clusters of expression trends. The 4h activated cluster includes OC master transcription factors (TFs; cMyc, Nfatc1, Fosl), and is enriched in OC-defining pathways. Notably, by 14h the majority of the genes required for mature OC formation and activation are already highly expressed (e.g. Ctsk, Mmp9). The downregulated clusters include monocyte defining TFs (e.g. Irf8, Mafb and Bcl6). These RANKL-dependent transcriptome changes are completely abrogated by iBET151, highlighting the critical role of Brd2-4 in osteoclastogenesis. Differential chromatin binding analysis upon RANKL induction revealed an overall enhanced Brd2-4 binding at already existing or de novo gained sites. This was more pronounced for Brd2&4 and much less for Brd3, with differentially binding sites (DBS) comprising 50% and 20% respectively of all binding sites in RANKL-treated cells. For Brd2&3, DBS were primarily distributed at promoters and for Brd4 at intergenic, candidate enhancers regions. Notably, nearly all gained DBS were sensitive to and abrogated by iBET151. Combinatorial profiling of Brd2 and Brd4 showed that almost half of Brd2 DBS peaks overlap with Brd4 (47%; 897/1896), while only 24% (766/3234) of Brd4 DBS peaks are co-occupied by Brd2. Transcriptome and Brd2&4 DBS integration in combination with motif enrichment analysis, identified genes that are predicted to be regulated by Brd2 and/or Brd4. EnrichR analysis suggests that enhanced binding of Brd2&4, singly or in combination, is required for activation of the critical OC lineage-specific and repression of the macrophage-defining transcriptional programs highlighting the non-redundant roles of Brd2&4 in OC development. Cell lineage commitment often requires 'commissioning' of cell-specific super-enhancers (SE). Combined analysis of genome-wide Brd4/H3K27ac profiles identified 678 RANKL-induced SE and their respective target genes. Further, 110 of these SE showed enhanced Brd4 binding in 2 peaks: 20/110 were linked to significantly up- and 90/100 to down-regulated genes. The repressed genes were significantly enriched to previously described Irf8, MafB and RunX1 targets, suggesting a critical role of SE in the repression of the monocyte/macrophage inflammatory programme during OC lineage commitment. Strikingly, among top hits, we detected a SE linked to the regulation of cMyc. To further investigate its role in OC development, we obtained the cistrome of cMyc after RANKL induction. We identified 560 binding sites which were highly enriched in cMyc, Max, Fli1, Fosl2 and Irf8 motifs. Cistrome-transcriptome integration suggested direct activation of 141 and repression of 52 genes by cMyc in response to RANKL; these are enriched in ribosome biogenesis pathways and Irf8-dependent targets respectively. CONCLUSIONS: Myc and Brd4 mark SE that repress an Irf8-dependent transcriptional programme, a requirement for OC lineage commitment. The non-redundant roles of Brd2&4 suggest that selective targeting of either could inhibit aberrant OC activation associated with MM. Disclosures Caputo: GSK: Research Funding. Auner:Amgen: Other: Consultancy and Research Funding; Takeda: Consultancy; Karyopharm: Consultancy. Karadimitris:GSK: Research Funding.

Description: ZBP1 is an inducible nucleic acid (NA) sensor that is activated when pathogen NA bind to its Zα and Zβ domains. ZBP1 is required for TBK1-dependent phosphorylation of the transcription factor IRF3 (pIRF3) followed by its direct activation of type I interferon genes. However, the role, if any, of ZBP1 in tumour biology is not known. By searching for genes selectively expressed in multiple myeloma (MM) we identified ZBP1 mRNA expressed in 29 MM cell lines (MMCL) but not in 〉1000 other cancer cell lines (CCLE dataset); ZBP1 was expressed in all 766 patient myeloma PC (CoMMpass dataset) but not in normal blood cells (Blueprint) or 53 healthy tissues (GTex). We confirmed expression of ZBP1 mRNA and/or protein in MMCL, primary human and murine germinal centre B (GCB) and plasma cells (PC) as well as in myeloma PC. By inducing T cell-dependent humoral immune responses after ip alum-NP-KLH immunisation, we explored the role of selective and constitutive expression of Zbp1 in GCB to PC transition. We found no differences in the frequency of splenic GCB cells and PC between control WT and Zbp1-/- mice and in GCB cell frequency between immunised WT and Zbp1-/- mice. However, compared to WT, the increase in PC frequency in immunised Zbp1-/- mice was 50% lower (n=10/group, p