Description: Genome-wide association studies have identified susceptibility loci for esophageal squamous cell carcinoma (ESCC). We conducted a meta-analysis of all single-nucleotide polymorphisms (SNPs) that showed nominally significant P -values in two previously published genome-wide scans that included a total of 2961 ESCC cases and 3400 controls. The meta-analysis revealed five SNPs at 2q33 with P 〈 5 x 10 –8 , and the strongest signal was rs13016963, with a combined odds ratio (95% confidence interval) of 1.29 (1.19–1.40) and P = 7.63 x 10 –10 . An imputation analysis of 4304 SNPs at 2q33 suggested a single association signal, and the strongest imputed SNP associations were similar to those from the genotyped SNPs. We conducted an ancestral recombination graph analysis with 53 SNPs to identify one or more haplotypes that harbor the variants directly responsible for the detected association signal. This showed that the five SNPs exist in a single haplotype along with 45 imputed SNPs in strong linkage disequilibrium, and the strongest candidate was rs10201587, one of the genotyped SNPs. Our meta-analysis found genome-wide significant SNPs at 2q33 that map to the CASP8 / ALS2CR12/TRAK2 gene region. Variants in CASP8 have been extensively studied across a spectrum of cancers with mixed results. The locus we identified appears to be distinct from the widely studied rs3834129 and rs1045485 SNPs in CASP8 . Future studies of esophageal and other cancers should focus on comprehensive sequencing of this 2q33 locus and functional analysis of rs13016963 and rs10201587 and other strongly correlated variants.

Description: Motivation: The accurate prediction of disease status is a central challenge in clinical cancer research. Microarray-based gene biomarkers have been identified to predict outcome and outperform traditional clinical parameters. However, the robustness of the individual gene biomarkers is questioned because of their little reproducibility between different cohorts of patients. Substantial progress in treatment requires advances in methods to identify robust biomarkers. Several methods incorporating pathway information have been proposed to identify robust pathway markers and build classifiers at the level of functional categories rather than of individual genes. However, current methods consider the pathways as simple gene sets but ignore the pathway topological information, which is essential to infer a more robust pathway activity. Results: Here, we propose a directed random walk (DRW)-based method to infer the pathway activity. DRW evaluates the topological importance of each gene by capturing the structure information embedded in the directed pathway network. The strategy of weighting genes by their topological importance greatly improved the reproducibility of pathway activities. Experiments on 18 cancer datasets showed that the proposed method yielded a more accurate and robust overall performance compared with several existing gene-based and pathway-based classification methods. The resulting risk-active pathways are more reliable in guiding therapeutic selection and the development of pathway-specific therapeutic strategies. Availability: DRW is freely available at http://210.46.85.180:8080/DRWPClass/ Contact: lixia@hrbmu.edu.cn or dm42298@126.com Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Despite growing consensus that long intergenic non-coding ribonucleic acids (lincRNAs) are modulators of cancer, the knowledge about the deoxyribonucleic acid (DNA) methylation patterns of lincRNAs in cancers remains limited. In this study, we constructed DNA methylation profiles for 4629 tumors and 705 normal tissue samples from 20 different types of human cancer by reannotating data of DNA methylation arrays. We found that lincRNAs had different promoter methylation patterns in cancers. We classified 2461 lincRNAs into two categories and three subcategories, according to their promoter methylation patterns in tumors. LincRNAs with resistant methylation patterns in tumors had conserved transcriptional regulation regions and were ubiquitously expressed across normal tissues. By integrating cancer subtype data and patient clinical information, we identified lincRNAs with promoter methylation patterns that were associated with cancer status, subtype or prognosis for several cancers. Network analysis of aberrantly methylated lincRNAs in cancers showed that lincRNAs with aberrant methylation patterns might be involved in cancer development and progression. The methylated and demethylated lincRNAs identified in this study provide novel insights for developing cancer biomarkers and potential therapeutic targets.

Description: Endometriosis is a complex and enigmatic disease that arises from the interplay among multiple genetic and environmental factors. The defining feature of endometriosis is the deposition and growth of endometrial tissues at sites outside of the uterine cavity. Studies to date have established that endometriosis is heritable but have not addressed the causal genetic variants for this disease. Here, we conducted whole-exome sequencing to comprehensively search for somatic mutations in both eutopic and ectopic endometrium from 16 endometriosis patients and five normal control patients using laser capture microdissection. We compared the mutational landscape of ectopic endometrium with the corresponding eutopic sample from endometriosis patients compared with endometrium from normal women and identified previously unreported mutated genes and pathway alternations. Statistical analysis of exome data identified that most genes were specifically mutated in both eutopic and ectopic endometrium cells. In particular, genes that are involved in biological adhesion, cell–cell junctions, and chromatin-remodeling complex(es) were identified, which partially supports the retrograde menstruation theory that proposes that endometrial cells are refluxed through the fallopian tubes during menstruation and implanted onto the peritoneum or pelvic organs. Conspicuously, when we compared exomic mutation data for paired eutopic and ectopic endometrium, we identified a mutational signature in both endometrial types for which no overlap in somatic single nucleotide variants were observed. These mutations occurred in a mutually exclusive manner, likely because of the discrepancy in endometriosis pathology and physiology, as eutopic endometrium rapidly regrows, and ectopic endometrial growth is inert. Our findings provide, to our knowledge, an unbiased view of the landscape of genetic alterations in endometriosis and vital information for indicating that genetic alterations in cytoskeletal and chromatin-remodeling proteins could be involved in the pathogenesis of endometriosis, thus implicating a novel therapeutic possibility for endometriosis.

Description: Previous studies have implicated cancer stem cells in tumor recurrence and revealed that the stem cell gene SOX2 plays an important role in the tumor cell resistance to apoptosis. Nonetheless, the mechanism by which SOX2 regulates apoptosis signals remained undefined. Here, we demonstrated the surprising finding that silencing of the SOX2 gene effectively induces apoptosis via the activation of death receptor and mitochondrial signaling pathways in human non-small cell lung cancer cells. Unexpectedly, reverse transcription–PCR analysis suggested that downregulation of SOX2 leads to activation of MAP4K4, previously implicated in cell survival. Evaluation of the apoptotic pathways revealed an increased expression of key inducers of apoptosis, including tumor necrosis factor-α and p53, with concurrent attenuation of Survivin. Although p53 appeared dispensable for this pathway, the loss of Survivin in SOX2-deficient cells appeared critical for the observed MAP4K4 induced cell death. Rescue experiments revealed that SOX2 -silencing-mediated killing was blocked by ectopic expression of Survivin, or by reduction of MAP4K4 expression. Clinically, expressions of Survivin and SOX2 were highly correlated with each other. The results reveal a key target of SOX2 expression and highlight the unexpected context-dependent role for MAP4K4, a pluripotent activator of several mitogen-activated protein kinase pathways, in regulating tumor cell survival.

Description: Motivation: Alzheimer’s disease (AD) is a severe neurodegenerative disease of the central nervous system that may be caused by perturbation of regulatory pathways rather than the dysfunction of a single gene. However, the pathology of AD has yet to be fully elucidated. Results: In this study, we systematically analyzed AD-related mRNA and miRNA expression profiles as well as curated transcription factor (TF) and miRNA regulation to identify active TF and miRNA regulatory pathways in AD. By mapping differentially expressed genes and miRNAs to the curated TF and miRNA regulatory network as active seed nodes, we obtained a potential active subnetwork in AD. Next, by using the breadth-first-search technique, potential active regulatory pathways, which are the regulatory cascade of TFs, miRNAs and their target genes, were identified. Finally, based on the known AD-related genes and miRNAs, the hypergeometric test was used to identify active pathways in AD. As a result, nine pathways were found to be significantly activated in AD. A comprehensive literature review revealed that eight out of nine genes and miRNAs in these active pathways were associated with AD. In addition, we inferred that the pathway hsa-miR-146a-〉STAT1-〉MYC, which is the source of all nine significantly active pathways, may play an important role in AD progression, which should be further validated by biological experiments. Thus, this study provides an effective approach to finding active TF and miRNA regulatory pathways in AD and can be easily applied to other complex diseases. Contact: lixia@hrbmu.edu.cn or lw2247@gmail.com . Supplementary information: Supplementary data are available at Bioinformatics online.

Description: Dysregulation of microRNA (miRNA) biogenesis is implicated in cancer development and progression. Dicer and Drosha are established regulators of miRNA biogenesis. In this study, we used a miRNA array to evaluate the miRNA expression profiles in nasopharyngeal carcinoma (NPC) samples. The significance analysis of microarrays showed a global downregulation of miRNA expression in NPC samples compared with normal nasopharyngeal epithelial tissues. Notably, miR-18a, a member of the oncogenic miR-17–92 cluster, was upregulated in the NPC samples and ell lines. Clinical parameter studies showed that higher levels of miR-18a correlated with NPC advanced stage, lymph node metastasis, Epstein-Barr virus infection and a higher death rate from NPC, indicating oncogenic roles in NPC development. The expression levels of miR-18a and Dicer1 were inversely related in NPC tissues. Further studies demonstrated that miR-18a negatively regulated Dicer1 by binding to the 3' untranslated regions of Dicer1. In vitro and in vivo biological function assays showed that miR-18a promoted the growth, migration and invasion of NPC cells by regulating Dicer1 expression, which caused the global downregulation of miRNA expression levels including miR-200 family and miR-143. Furthermore, we found that the epithelial mesenchymal transition marker E-cadherin and the oncogene K-Ras were aberrantly expressed after miR-18a transduction, and these alterations were directly induced by downregulation of the miR-200 family and miR-143. Collectively, our findings indicate that miR-18a plays an oncogenic role in the development of NPC by widespread downregulation of the miRNome and could be a potential therapeutic target for NPC.

Description: Basaltic lavas from Hainan Island near the northern edge of the South China Sea have an age range of between late Miocene (about 13 Ma) and Holocene, with a peak age of late Pliocene to middle Pleistocene. The basaltic province is dominated by tholeiites with subordinate alkali basalts. Most analysed samples display light rare earth element (LREE) enriched REE patterns and ocean island basalt (OIB)-like incompatible element distributions. The basalts contain abundant undeformed high-Mg olivine phenocrysts (up to Fo 90·7 ) that are high in CaO and MnO, indicating high-magnesian parental magmas. Independent barometers indicate that clinopyroxenes in the basalts crystallized over a wide range of pressures of 2–25 kbar (dominantly at 10–15 kbar) and that the melt cooled from about 1350°C to 1100°C during their crystallization. The compositional characteristics of the basalts indicate that their generation most probably involved both low-silica and high-silica melts, as represented by the alkali basalts and tholeiites, respectively. Our results show that the source region for the Hainan basalts is highly heterogeneous. The source for the tholeiites is mainly composed of peridotite and recycled oceanic crust, whereas the source for the bulk of the low-Th alkali basalts consists predominantly of peridotite and low-silica eclogite (garnet pyroxenite). Some high-Th (≥ 4 ppm) alkali basalts may have been produced by partial melting of low-silica garnet pyroxenite (eclogite). We estimated the primary melt compositions for the Hainan basalts using the most forsteritic olivine (Fo 90·7 ) composition and the most primitive bulk-rock samples (MgO 〉 9·0 wt % and CaO 〉8·0 wt %), assuming a constant Fe–Mg exchange partition coefficient of K D = 0·31 and Fe 3+ /Fe T = 0·1. The effective melting pressure ( P f ) and melting temperature ( T ) of the primary melts are P f = 18–32 kbar (weighted average = 23·8 ± 1·8 kbar) and T = 1420–1520°C for the tholeiites, and P f = 25–32 (weighted average = 28·3 ± 1·4 kbar) and T = 1480–1530°C for the alkali basalts. The P f – T data form an array that plots systematically above the dry lherzolite solidus but below the base of the lithosphere (~55 km) and intersects the dry peridotite solidus at a pressure of about 50 kbar. The mantle potential temperature beneath Hainan Island, based on the estimate primary melt compositions, varies from about 1500 to 1580°C with a weighted average of 1541 ± 10°C. The high-magnesian olivine phenocrysts, high mantle potential temperature, and the presence of recycled oceanic crust in the source region provide independent support for the Hainan plume model that has previously been proposed largely based on geophysical observations. The Hainan plume thus provides a rare example of a young mantle plume associated with deep slab subduction.

Description: Inorganic arsenic is a well-documented human carcinogen associated with cancers of the skin, lung, liver, and bladder. However, the underlying mechanisms explaining the tumorigenic role of arsenic are not well understood. The present study explored a potential mechanism of cell transformation induced by arsenic exposure. Exposure to a low dose (0.5 μm) of arsenic trioxide (As2O3) caused transformation of BALB/c 3T3 cells. In addition, in a xenograft mouse model, tumor growth of the arsenic-induced transformed cells was dramatically increased. In arsenic-induced transformed cells, polycomb group (PcG) proteins, including BMI1 and SUZ12, were activated resulting in enhanced histone H3K27 tri-methylation levels. On the other hand, tumor suppressor p16INK4a and p19ARF mRNA and protein expression were dramatically suppressed. Introduction of small hairpin (sh) RNA-BMI1 or -SUZ12 into BALB/c 3T3 cells resulted in suppression of arsenic-induced transformation. Histone H3K27 tri-methylation returned to normal in BMI1- or SUZ12-knockdown BALB/c 3T3 cells compared with BMI1- or SUZ12-wildtype cells after arsenic exposure. As a consequence, the expression of p16INK4a and p19ARF was recovered in arsenic-treated BMI1- or SUZ12-knockdown cells. Thus, arsenic-induced cell transformation was blocked by inhibition of PcG function. Taken together, these results strongly suggest that the polycomb proteins, BMI1 and SUZ12 are required for cell transformation induced by organic arsenic exposure.

Description: BMP2 (bone morphogenetic protein 2) is known to activate unfolded protein response signaling molecules, including XBP1S and ATF6. However, the influence on XBP1S and ATF6 in BMP2-induced chondrocyte differentiation has not yet been elucidated. In this study, we demonstrate that BMP2 mediates mild endoplasmic reticulum stress-activated ATF6 and directly regulates XBP1S splicing in the course of chondrogenesis. XBP1S is differentially expressed during BMP2-stimulated chondrocyte differentiation and exhibits prominent expression in growth plate chondrocytes. This expression is probably due to the activation of the XBP1 gene by ATF6 and splicing by IRE1a. ATF6 directly binds to the 5′-flanking regulatory region of the XBP1 gene at its consensus binding elements. Overexpression of XBP1S accelerates chondrocyte hypertrophy, as revealed by enhanced expression of type II collagen, type X collagen, and RUNX2; however, knockdown of XBP1S via the RNAi approach abolishes hypertrophic chondrocyte differentiation. In addition, XBP1S associates with RUNX2 and enhances RUNX2-induced chondrocyte hypertrophy. Altered expression of XBP1S in chondrocyte hypertrophy was accompanied by altered levels of IHH (Indian hedgehog) and PTHrP (parathyroid hormone-related peptide). Collectively, XBP1S may be a novel regulator of hypertrophic chondrocyte differentiation by 1) acting as a cofactor of RUNX2 and 2) affecting IHH/PTHrP signaling.