Integrative analysis of genomic and epigenomic regulation reveals miRNA mediated tumor heterogeneity and immune evasion in lower grade glioma.

The expression dysregulation of microRNAs (miRNA) has been widely reported during cancer development, however, the underling mechanism remains largely unanswered. In the present work, we performed a systematic integrative study for genome-wide DNA methylation, copy number variation and miRNA expression data to identify mechanisms underlying miRNA dysregulation in lower grade glioma. We identify 719 miRNAs whose expression was associated with alterations of copy number variation or promoter methylation. Integrative multi-omics analysis revealed four subtypes with differing prognoses. These glioma subtypes exhibited distinct immune-related characteristics as well as clinical and genetic features. By construction of a miRNA regulatory network, we identified candidate miRNAs associated with immune evasion and response to immunotherapy. Finally, eight prognosis related miRNAs were validated to promote cell migration, invasion and proliferation through in vitro experiments. Our study reveals the crosstalk among DNA methylation, copy number variation and miRNA expression for immune regulation in glioma, and could have important implications for patient stratification and development of biomarkers for immunotherapy approaches.

The defense of Shangri-La: Protecting isolated communities by periodic infection screening in the worst future pandemic.

In the worst future pandemic, effective vaccines and medicines could be unavailable for a long time. In such circumstances, it is necessary to evaluate whether a periodic screening can protect isolated communities and critical facilities and avoid a complete shutdown. In this study, we introduced an epidemiological model that included the essential parameters of infection transmission and screening. With varying parameters, we studied the dynamics of viral infection in the isolated communities. In the scenario with a periodic infection screening once per 3 days and a viral basic reproduction number 3.0, >85% of the infection waves have a duration <7 days and the infection size in each of the waves is generally <4 individuals when the efficiency of infection discovery is 0.9 in the screening. When the period of screening was elongated to once per 7 days, the cases of infection dramatically increased to 5 folds of that mentioned previously. Further, with a weak discovery efficiency of 0.7 and the aforementioned low screening frequency, the spread of infection would be out of control. Our study suggests that frequent periodic screening is capable of controlling a future epidemic in isolated communities without other measures.

Insights into the role of long non-coding RNAs in DNA methylation mediated transcriptional regulation.

DNA methylation is one of the most important epigenetic mechanisms that governing regulation of gene expression, aberrant DNA methylation patterns are strongly associated with human malignancies. Long non-coding RNAs (lncRNAs) have being discovered as a significant regulator on gene expression at the epigenetic level. Emerging evidences have indicated the intricate regulatory effects between lncRNAs and DNA methylation. On one hand, transcription of lncRNAs are controlled by the promoter methylation, which is similar to protein coding genes, on the other hand, lncRNA could interact with enzymes involved in DNA methylation to affect the methylation pattern of downstream genes, thus regulating their expression. In addition, circular RNAs (circRNAs) being an important class of noncoding RNA are also found to participate in this complex regulatory network. In this review, we summarize recent research progress on this crosstalk between lncRNA, circRNA, and DNA methylation as well as their potential functions in complex diseases including cancer. This work reveals a hidden layer for gene transcriptional regulation and enhances our understanding for epigenetics regarding detailed mechanisms on lncRNA regulatory function in human cancers.

Tracing Bai-Yue Ancestry in Aboriginal Li People on Hainan Island.

As the most prevalent aboriginal group on Hainan Island located between South China and the mainland of Southeast Asia, the Li people are believed to preserve some unique genetic information due to their isolated circumstances, although this has been largely uninvestigated. We performed the first whole-genome sequencing of 55 Hainan Li (HNL) individuals with high coverage (∼30-50×) to gain insight into their genetic history and potential adaptations. We identified the ancestry enriched in HNL (∼85%) is well preserved in present-day Tai-Kadai speakers residing in South China and North Vietnam, that is, Bai-Yue populations. A lack of admixture signature due to the geographical restriction exacerbated the bottleneck in the present-day HNL. The genetic divergence among Bai-Yue populations began ∼4,000-3,000 years ago when the proto-HNL underwent migration and the settling of Hainan Island. Finally, we identified signatures of positive selection in the HNL, some outstanding examples included FADS1 and FADS2 related to a diet rich in polyunsaturated fatty acids. In addition, we observed that malaria-driven selection had occurred in the HNL, with population-specific variants of malaria-related genes (e.g., CR1) present. Interestingly, HNL harbors a high prevalence of malaria leveraged gene variants related to hematopoietic function (e.g., CD3G) that may explain the high incidence of blood disorders such as B-cell lymphomas in the present-day HNL. The results have advanced our understanding of the genetic history of the Bai-Yue populations and have provided new insights into the adaptive scenarios of the Li people.

dbDEMC 3.0: Functional Exploration of Differentially Expressed miRNAs in Cancers of Human and Model Organisms.

MicroRNAs (miRNAs) are important regulators in gene expression. The dysregulation of miRNA expression is widely reported in the transformation from physiological to pathological states of cells. A large number of differentially expressed miRNAs (DEMs) have been identified in various human cancers by using high-throughput technologies, such as microarray and miRNA-seq. Through mining of published studies with high-throughput experiment information, the database of DEMs in human cancers (dbDEMC) was constructed with the aim of providing a systematic resource for the storage and query of the DEMs. Here we report an update of the dbDEMC to version 3.0, which contains two-fold more data entries than the second version and now includes also data from mice and rats. The dbDEMC 3.0 contains 3268 unique DEMs in 40 different cancer types. The current datasets for differential expression analysis have expanded to 9 generalized categories. Moreover, the current release integrates functional annotations of DEMs obtained by using experimentally validated targets. The annotations can be of great benefit to the intensive analysis of the roles of DEMs in cancer. In summary, dbDEMC 3.0 provides a valuable resource for characterizing molecular functions and regulatory mechanisms of DEMs in human cancers. The dbDEMC 3.0 is freely accessible at https://www.biosino.org/dbDEMC.

IFNα Potentiates Anti-PD-1 Efficacy by Remodeling Glucose Metabolism in the Hepatocellular Carcinoma Microenvironment.

The overall response rate for anti-PD-1 therapy remains modest in hepatocellular carcinoma (HCC). We found that a combination of IFNα and anti-PD-1-based immunotherapy resulted in enhanced antitumor activity in patients with unresectable HCC. In both immunocompetent orthotopic and spontaneous HCC models, IFNα therapy synergized with anti-PD-1 and the combination treatment led to significant enrichment of cytotoxic CD27+CD8+ T cells. Mechanistically, IFNα suppressed HIF1α signaling by inhibiting FosB transcription in HCC cells, resulting in reduced glucose consumption capacity and consequentially establishing a high-glucose microenvironment that fostered transcription of the T-cell costimulatory molecule Cd27 via mTOR-FOXM1 signaling in infiltrating CD8+ T cells. Together, these data reveal that IFNα reprograms glucose metabolism within the HCC tumor microenvironment, thereby liberating T-cell cytotoxic capacities and potentiating the PD-1 blockade-induced immune response. Our findings suggest that IFNα and anti-PD-1 cotreatment is an effective novel combination strategy for patients with HCC. SIGNIFICANCE: Our study supports a role of tumor glucose metabolism in IFNα-mediated antitumor immunity in HCC, and tumor-infiltrating CD27+CD8+ T cells may be a promising biomarker for stratifying patients for anti-PD-1 therapy. See related commentary by Kao et al., p. 1615. This article is highlighted in the In This Issue feature, p. 1599.

Matrisome provides a supportive microenvironment for oral squamous cell carcinoma progression.

Oral squamous cell carcinoma (OSCC) is a common pernicious tumor in the head and neck regions. However, the function of tumor extracellular matrix (ECM) has not been elucidated. A tissue engineering method was applied for remodeling ECM through decellularization. The cellular components were removed, and the biological composition was mostly preserved. Proteomics was performed to analyze the characterization between normal and tumor ECM. According to LC-MS/MS results, 26 proteins just showed in tumor ECM, and 14 proteins only showed in late-stage tumor ECM. KEGG pathway analysis showed that most variant proteins were linked to metabolic regulation and tumor immunity (such as SCC-Ag1, LOX). To affirm the influence of tumor ECM on the progression of OSCC, tumor cells and macrophages were co-cultured with ECM scaffold. Marked differences in proliferation, apoptosis, and migration of OSCC cells were observed between tumor and normal ECM. Tumor ECM polarized macrophages towards an anti-inflammatory phenotype (higher IL-10 and CD68, and relatively lower CD86 and IL1-ß). Collectively, these findings suggest that tumor ECM served as a permissive role in OSCC progression. SIGNIFICANCE: The variation between OSCC ECM and normal ECM confirm tumor ECM plays a significant role in OSCC deterioration, which is conducive to exploring the occurrence and progression mechanisms of OSCC, and further improving the curative effect of this disease.

Bone-targeted erythrocyte-cancer hybrid membrane-camouflaged nanoparticles for enhancing photothermal and hypoxia-activated chemotherapy of bone invasion by OSCC.

BACKGROUND: Jaw bones are the most common organs to be invaded by oral malignancies, such as oral squamous cell carcinoma (OSCC), because of their special anatomical relationship. Various serious complications, such as pathological fractures and bone pain can significantly decrease the quality of life or even survival outcomes for a patient. Although chemotherapy is a promising strategy for bone invasion treatment, its clinical applications are limited by the lack of tumor-specific targeting and poor permeability in bone tissue. Therefore, it is necessary to develop a smart bone and cancer dual targeting drug delivery platform. RESULTS: We designed a dual targeting nano-biomimetic drug delivery vehicle Asp8[H40-TPZ/IR780@(RBC-H)] that has excellent bone and cancer targeting as well as immune escape abilities to treat malignancies in jaw bones. These nanoparticles were camouflaged with a head and neck squamous cell carcinoma WSU-HN6 cell (H) and red blood cell (RBC) hybrid membrane, which were modified by an oligopeptide of eight aspartate acid (Asp8). The spherical morphology and typical core-shell structure of biomimetic nanoparticles were observed by transmission electron microscopy. These nanoparticles exhibited the same surface proteins as those of WSU-HN6 and RBC. Flow cytometry and confocal microscopy showed a greater uptake of the biomimetic nanoparticles when compared to bare H40-PEG nanoparticles. Biodistribution of the nanoparticles in vivo revealed that they were mainly localized in the area of bone invasion by WSU-HN6 cells. Moreover, the Asp8[H40-TPZ/IR780@(RBC-H)] nanoparticles exhibited effective cancer growth inhibition properties when compared to other TPZ or IR780 formulations. CONCLUSIONS: Asp8[H40-TPZ/IR780@(RBC-H)] has bone targeting, tumor-homing and immune escape abilities, therefore, it is an efficient multi-targeting drug delivery platform for achieving precise anti-cancer therapy during bone invasion.

Pan-cancer characterization of long non-coding RNA and DNA methylation mediated transcriptional dysregulation.

BACKGROUND: Disruption of DNA methylation (DNAm) is one of the key signatures of cancer, however, detailed mechanisms that alter the DNA methylome in cancer remain to be elucidated. METHODS: Here we present a novel integrative analysis framework, called MeLncTRN (Methylation mediated LncRNA Transcriptional Regulatory Network), that integrates genome-wide transcriptome, DNA methylome and copy number variation profiles, to systematically identify the epigenetically-driven lncRNA-gene regulation circuits across 18 cancer types. FINDING: We show that a significant fraction of the aberrant DNAm and gene expression landscape in cancer is associated with long noncoding RNAs (lncRNAs). We reveal distinct types of regulation between lncRNA modulators and target genes that are operative in either only specific cancers or across cancers. Functional studies identified a common theme of cancer hallmarks that lncRNA modulators may participate in. The coupled lncRNA gene interactions via DNAm also serve as markers for classifications of cancer subtypes with different prognoses. INTERPRETATION: Our study reveals a vital layer of DNAm and associated expression regulation for many cancer-related genes and we also provide a valuable database resource for interrogating epigenetically mediated lncRNA-gene interactions in cancer. FUNDING: National Natural Science Foundation of China [91959106, 31871255].

CMTTdb: the cancer molecular targeted therapy database.

BACKGROUND: The cancer molecular targeted therapy has achieved unprecedented progress in the past decade and is thought to be the most promising direction for cancer treatment in future. As the fast growing of the clinical trials of targeted anticancer agents for different cancer types, it is critical to collect and integrate such information to guide clinical practice. METHODS: We constructed the Cancer Molecular Targeted Therapy database (CMTTdb) to store and retrieve molecular targeted therapy data about randomized clinical trials (RCTs) of targeted agents and also accompanied targets, biomarkers, targeted cancer subtypes, etc. RESULTS: Different with some existing resources, CMTTdb particularly focuses on clinical application of the trails. Design of the trails, such as treatment modalities (monotherapy or combination with other therapies), as well as results on clinical efficacy parameters, adverse events are also collected. In this current version, CMTTdb contains data for 1,088 clinical trials which cover 165 agents, 80 targets, 15 cancer types (95 molecular subtypes and 56 histological or cytological subtypes) from public literatures. This database is freely available at http://www.biosino.org/CMTTdb. A user-friendly web interface was designed so that these data can be easily retrieved. CONCLUSIONS: CMTTdb will be a valuable source for providing access to information of clinical trials on the rapidly growing number of novel targeted agent and be useful in guiding oncologists for the optimization of the therapy strategy for cancer treatment.

The DRD4 exon III VNTR, bupropion, and associations with prospective abstinence.

INTRODUCTION: DRD4 Exon III Variable Number of Tandem Repeat (VNTR) variation was found to interact with bupropion to influence prospective smoking abstinence, in a recently published longitudinal analyses of N = 331 individuals from a randomized double-blind placebo-controlled trial of bupropion and intensive cognitive-behavioral mood management therapy. METHODS: We used univariate, multivariate, and longitudinal logistic regression to evaluate gene, treatment, time, and interaction effects on point prevalence and continuous abstinence at end of treatment, 6 months, and 12 months, respectively, in N = 416 European ancestry participants in a double-blind pharmacogenetic efficacy trial randomizing participants to active or placebo bupropion. Participants received 10 weeks of pharmacotherapy and 7 sessions of behavioral therapy, with a target quit date 2 weeks after initiating both therapies. VNTR genotypes were coded with the long allele dominant resulting in 4 analysis categories. Covariates included demographics, dependence measures, depressive symptoms, and genetic ancestry. We also performed genotype-stratified secondary analyses. RESULTS: We observed significant effects of time in longitudinal analyses of both abstinence outcomes, of treatment in individuals with VNTR long allele genotypes for both abstinence outcomes, and of covariates in some analyses. We observed non-significantly larger differences in active versus placebo effect sizes in individuals with VNTR long allele genotypes than in individuals without the VNTR long allele, in the directions previously reported. CONCLUSIONS: VNTR by treatment interaction differences between these and previous analyses may be attributable to insufficient size of the replication sample. Analyses of multiple randomized clinical trials will enable identification and validation of factors mediating treatment response.

Intron 3 16 bp duplication polymorphism of TP53 contributes to cancer susceptibility: a meta-analysis.

A few genetic polymorphisms of TP53 are known to have a significant effect on cancer susceptibility. Intron 3 16 bp duplication polymorphism of TP53 has been reported to be associated with breast cancer, colorectal cancer, lung cancer and other cancers, but the reported results remain inconclusive. The present study, a meta-analysis including a total of 9801 cases and 10,391 controls from 26 studies, revealed that the 16 bp insertion (Ins) allele is significantly associated with an increased cancer risk in overall analysis [Ins/Ins + deletion (Del)/Ins versus Del/Del: odds ratio (OR) = 1.14, 95% confidence interval (CI) = 1.02-1.27, P = 0.02; Ins/Ins versus Del/Del: OR = 1.35, 95% CI = 1.11-1.63, P = 0.002; Del/Ins versus Del/Del: OR = 1.10, 95% CI = 0.98-1.23, P = 0.11.), particularly in breast cancer subgroup (Ins/Ins + Del/Ins versus Del/Del: OR = 1.16, 95% CI = 1.03-1.31, P = 0.02; Ins/Ins versus Del/Del: OR = 1.81, 95% CI = 1.30-2.52, P < 0.001; Del/Ins versus Del/Del: OR = 1.10, 95% CI = 0.97-1.25, P = 0.13). The relative risks to the colorectal and lung cancers increased but their association power was relatively weak, which may result from a limited number of studies of these two cancer types. These results suggest that intron 3 16 bp duplication polymorphism of TP53 is potentially an important and clinically relevant genetic marker contributing to cancer susceptibility.

Nicotine withdrawal sensitivity, linkage to chr6q26, and association of OPRM1 SNPs in the SMOking in FAMilies (SMOFAM) sample.

BACKGROUND: Nicotine withdrawal symptoms are related to smoking cessation. A Rasch model has been used to develop a unidimensional sensitivity score representing multiple correlated measures of nicotine withdrawal. A previous autosome-wide screen identified a nonparametric linkage (NPL) log-likelihood ratio (LOD) score of 2.7 on chromosome 6q26 for the sum of nine withdrawal symptoms. METHODS: The objectives of these analyses were to (a) assess the influence of nicotine withdrawal sensitivity on relapse, (b) conduct autosome-wide NPL analysis of nicotine withdrawal sensitivity among 158 pedigrees with 432 individuals with microsatellite genotypes and nicotine withdrawal scores, and (c) explore family-based association of single nucleotide polymorphism (SNP) at the mu opioid receptor candidate gene (OPRM1) with nicotine withdrawal sensitivity in 172 nuclear pedigrees with 419 individuals with both SNP genotypes and nicotine withdrawal scores. RESULTS: An increased risk for relapse was associated with nicotine withdrawal sensitivity score (odds ratio, 1.25; 95% confidence interval, 1.10-1.42). A maximal NPL LOD score of 3.15, suggestive of significant linkage, was identified at chr6q26 for nicotine withdrawal sensitivity. Evaluation of 18 OPRM1 SNPs via the family-based association test with the nicotine withdrawal sensitivity score identified eight tagging SNPs with global P values <0.05 and false discovery rate Q values <0.06. CONCLUSION: An increased risk of relapse, suggestive linkage at chr6q26, and nominally significant association with multiple OPRM1 SNPs were found with Rasch-modeled nicotine withdrawal sensitivity scores in a multiplex smoking pedigree sample. Future studies should attempt to replicate these findings and investigate the relationship between nicotine withdrawal symptoms and variation at OPRM1.

Dopamine genes and nicotine dependence in treatment-seeking and community smokers.

We utilized a cohort of 828 treatment-seeking self-identified white cigarette smokers (50% female) to rank candidate gene single nucleotide polymorphisms (SNPs) associated with the Fagerström Test for Nicotine Dependence (FTND), a measure of nicotine dependence which assesses quantity of cigarettes smoked and time- and place-dependent characteristics of the respondent's smoking behavior. A total of 1123 SNPs at 55 autosomal candidate genes, nicotinic acetylcholine receptors and genes involved in dopaminergic function, were tested for association to baseline FTND scores adjusted for age, depression, education, sex, and study site. SNP P-values were adjusted for the number of transmission models, the number of SNPs tested per candidate gene, and their intragenic correlation. DRD2, SLC6A3, and NR4A2 SNPs with adjusted P-values <0.10 were considered sufficiently noteworthy to justify further genetic, bioinformatic, and literature analyses. Each independent signal among the top-ranked SNPs accounted for approximately 1% of the FTND variance in this sample. The DRD2 SNP appears to represent a novel association with nicotine dependence. The SLC6A3 SNPs have previously been shown to be associated with SLC6A3 transcription or dopamine transporter density in vitro, in vivo, and ex vivo. Analysis of SLC6A3 and NR4A2 SNPs identified a statistically significant gene-gene interaction (P=0.001), consistent with in vitro evidence that the NR4A2 protein product (NURR1) regulates SLC6A3 transcription. A community cohort of N=175 multiplex ever-smoking pedigrees (N=423 ever smokers) provided nominal evidence for association with the FTND at these top ranked SNPs, uncorrected for multiple comparisons.

Genetic and environmental influences on the ratio of 3'hydroxycotinine to cotinine in plasma and urine.

OBJECTIVES: The ratio of trans-3'hydroxycotinine/cotinine (3HC/COT) is a marker of CYP2A6 activity, an important determinant of nicotine metabolism. This analysis sought to conduct a combined genetic epidemiologic and pharmacogenetic investigation of the 3HC/COT ratio in plasma and urine. METHODS: One hundred and thirty-nine twin pairs [110 monozygotic and 29 dizygotic] underwent a 30-min infusion of stable isotope-labelled nicotine and its major metabolite, cotinine, followed by an 8-h in-hospital stay. Blood and urine samples were taken at regular intervals for analysis of nicotine, cotinine, and metabolites. DNA was genotyped to confirm zygosity and for variation in the gene for the primary nicotine metabolic enzyme, CYP2A6 (variants genotyped: *1B, *1 x 2, *2, *4, *9, *12). Univariate biometric analyses quantified genetic and environmental influences on each measure in the presence and absence of covariates, including measured CYP2A6 genotype. RESULTS: There was a substantial amount of variation in the free 3HC/COT ratio in plasma (6 h postinfusion) attributable to additive genetic influences (67.4%, 95% confidence interval=55.9-76.2%). The heritability estimate was reduced to 61.0 and 49.4%, respectively, after taking into account the effect of covariates and CYP2A6 genotype. In urine (collected over 8 h), the estimated amount of variation in the 3HC/COT ratio attributable to additive genetic influences was smaller (47.2%, 95% confidence interval=0-67.2%) and decreased to 44.6 and 42.0% after accounting for covariates and genotype. CONCLUSION: Additive genetic factors are prominent in determining variation in plasma 3HC/COT but less so in determining variation in urine 3HC/COT.

Genome-wide linkage of cotinine pharmacokinetics suggests candidate regions on chromosomes 9 and 11.

Characterizing cotinine pharmacokinetics is a useful way to study nicotine metabolism because the same liver enzyme is primarily responsible for the metabolism of both, and the clearances of nicotine and cotinine are highly correlated. We conducted a whole-genome linkage analysis to search for candidate regions influencing quantitative variation in cotinine pharmacokinetics in a large-scale pharmacokinetic study with 61 families containing 224 healthy adult participants. The strongest linkage signal was identified at 135 cM of chromosome 9 with LOD = 2.81 and P = 0.0002; two other suggestive linkage peaks appear at 31.4 and 73.5 cM of chromosome 11 with LOD = 1.96 (P = 0.0013) and 1.94 (P = 0.0014). The confidence level of the linkage between the three genome regions and cotinine pharmacokinetics is statistically significant with a genome-wide empirical probability of P = 0.029.

Cloning and functional analysis of human mTERFL encoding a novel mitochondrial transcription termination factor-like protein.

Serum plays an important role in the regulation of cell cycle and cell growth. To identify novel serum-inhibitory factors and study their roles in cell cycle regulation, we performed mRNA differential display analysis of U251 cells in the presence or absence of serum and cloned a novel gene encoding the human mitochondrial transcription termination factor-like protein (mTERFL). The full-length mTERFL cDNA has been isolated and the genomic structure determined. The mTERFL gene consists of three exons and encodes 385 amino acids with 52% sequence similarity to the human mitochondrial transcription termination factor (mTERF). However, mTERFL and mTERF have an opposite expression pattern in response to serum. The expression of mTERFL is dramatically inhibited by the addition of serum in serum-starved cells while the mTERF is rather induced. Northern blot analysis detected three mTERFL transcripts of 1.7, 3.2, and 3.5kb. Besides the 3.2kb transcript that is unique to skeletal muscle, other two transcripts express predominant in heart, liver, pancreas, and skeletal muscle. Expression of the GFP-mTERFL fusion protein in HeLa cells localized it to the mitochondria. Furthermore, ectopic expression of mTERFL suppresses cell growth and arrests cells in the G1 stage demonstrated by MTT and flow cytometry analysis. Collectively, our data suggest that mTERFL is a novel mTERF family member and a serum-inhibitory factor probably participating in the regulation of cell growth through the modulation of mitochondrial transcription.