Tumor suppressor genes (TSGs) are commonly inactivated by somatic mutation and/or promoter methylation; yet, recent high-throughput genomic studies have not identified key TSGs inactivated by both mechanisms. We pursued an integrated molecular analysis based on methylation binding domain sequencing (MBD-seq), 450K Methylation arrays, whole exome sequencing, and whole genome gene expression arrays in primary head and neck squamous cell carcinoma (HNSCC) tumors and matched uvulopalatopharyngoplasty tissue samples (UPPPs). We uncovered 186 downregulated genes harboring cancer specific promoter methylation including PAX1 and PAX5 and we identified 10 key tumor suppressor genes (GABRB3, HOXC12, PARP15, SLCO4C1, CDKN2A, PAX1, PIK3AP1, HOXC6, PLCB1, and ZIC4) inactivated by both promoter methylation and/or somatic mutation. Among the novel tumor suppressor genes discovered with dual mechanisms of inactivation, we found a high frequency of genomic and epigenomic alterations in the PAX gene family of transcription factors, which selectively impact canonical NOTCH and TP53 pathways to determine cell fate, cell survival, and genome maintenance. Our results highlight the importance of assessing TSGs at the genomic and epigenomic level to identify key pathways in HNSCC, deregulated by simultaneous promoter methylation and somatic mutations.

Using high-throughput analyses and the TRANSFAC database, we characterized TF signatures of head and neck squamous cell carcinoma (HNSCC) subgroups by inferential analysis of target gene expression, correcting for the effects of DNA methylation and copy number. Using this discovery pipeline, we determined that human papillomavirus-related (HPV+) and HPV− HNSCC differed significantly based on the activity levels of key TFs including AP1, STATs, NF-κB, and p53. Immunohistochemical analysis confirmed that HPV− HNSCC is characterized by co-activated STAT3 and NF-κB pathways, and functional studies demonstrate that this phenotype can be effectively targeted with combined anti-NF-κB and anti-STAT therapies. These discoveries correlate strongly with previous findings connecting STATs, NF-κB, and AP1 in HNSCC. We identified 5 top-scoring pair biomarkers from STATs, NF-κB and AP1 pathways that distinguish HPV+ from HPV− HNSCC based on TF activity, and validated these biomarkers on TCGA and on independent validation cohorts. We conclude that a novel approach to TF pathway analysis can provide insight into therapeutic targeting of patient subgroup for heterogeneous disease such as HNSCC.

Over a half million new cases of Head and Neck Squamous Cell Carcinoma (HNSCC) are diagnosed annually worldwide, however, 5 year overall survival is only 50% for HNSCC patients. Recently, high throughput technologies have accelerated the genome-wide characterization of HNSCC. However, comprehensive pipelines with statistical algorithms that account for HNSCC biology and perform independent confirmatory and functional validation of candidates are needed to identify the most biologically relevant genes. We applied outlier statistics to high throughput gene expression data, and identified 76 top-scoring candidates with significant differential expression in tumors compared to normal tissues. We identified 15 epigenetically regulated candidates by focusing on a subset of the genes with a negative correlation between gene expression and promoter methylation. Differential expression and methylation of 3 selected candidates (BANK1, BIN2, and DTX1) were confirmed in an independent HNSCC cohorts from Johns Hopkins and TCGA (The Cancer Genome Atlas). We further performed functional evaluation of NOTCH regulator, DTX1, which was downregulated by promoter hypermethylation in tumors, and demonstrated that decreased expression of DTX1 in HNSCC tumors maybe associated with NOTCH pathway activation and increased migration potential.

Recently, miR-23b has emerged as a promising new cancer biomarker but its role in lung cancer has not been established yet. Patients still do not respond well to available treatments, probably due to expression of multidrug resistance (MDR) proteins, such as P-gp, MRP and LRP/MVP. The aim of this study was to determine the role of miR-23b in non-small cell lung cancer (NSCLC) and its relationship to the patient outcome together with MDR transporter proteins. We immunohistochemically evaluated expression of P-gp, MRP and LRP/MVP and quantified the relative levels of miR-23b in 62 NSCLC patients´ samples. The prognostic significance of miR-23b and MDR proteins was tested by Kaplan-Meier and Cox-regression analysis. Our results showed that miR-23b is mostly downregulated in NSCLC samples (57/62) and that its upregulation in tumors is connected with longer progression-free survival (PFS; P = 0.065) and overall survival (OS; P = 0.048). The Cox proportional hazard model revealed that the risk of death or relapse in NSCLC patients with miR-23b downregulation increases together with LRP/MVP expression and both risks decrease with miR-23b upregulation (HRPFS = 4.342, PPFS = 0.022; HROS = 4.408, POS = 0.015). Our findings indicate that miR-23b, especially in combination with LRP/MVP expression, might serve as a suitable prognostic biomarker for NSCLC patients.

Next-generation sequencing (NGS) is increasingly used in transplantation settings, but also as a method of choice for in-depth analysis of population-specific HLA genetic architecture and its linkage to various diseases. With respect to complex ethnic admixture characteristic for East Croatian population, we aimed to investigate class-I (HLA-A,-B,-C) and class-II (HLA-DRB1,-DQA1,-DQB1) HLA diversity at the highest, 4-field resolution level in 120 healthy, unrelated, blood donor volunteers. Genomic DNA was extracted and HLA genotypes of class I and DQA1 genes were defined in full-length,-DQB1 from intron 1 to 3′ UTR, and-DRB1 from intron 1 to intron 4 (Illumina MiSeq platform, Omixon Twin algorithms, IMGT/HLA release 3.30.0_5). Linkage disequilibrium statistics, Hardy-Weinberg departures, and haplotype frequencies were inferred by exact tests and iterative Expectation-Maximization algorithm using PyPop 0.7.0 and Arlequin v3.5.2.2 software. Our data provide first description of 4-field allele and haplotype frequencies in Croatian population, revealing 192 class-I and class-II alleles and extended haplotypic combinations not apparent from the existing 2-field HLA reports from Croatia. This established reference database complements current knowledge of HLA diversity and should prove useful in future population studies, transplantation settings, and disease-associated HLA screening. Croatia is a Mediterranean, crescent-shaped south European country bordering Slovenia in the northwest, Hungary in the northeast, Serbia in the east, Bosnia and Herzegovina and Montenegro in the southeast, and Italy along the maritime border. Croatia consists of three major geomorphologic areas, which can be further broken down into five traditional districts based on history, topography, and economy; Istria and Dalmatia in the northern and southern Croatian littoral, Gorski Kotar in country's mountainous area, central continental Croatia, and Slavonia in the Pannonian basin in the east (Fig. 1). Slavonia territory was originally populated by the southern branch of the Indo-European Slavic populations in the 7 th century 1 , and has been a witness of significant population admixture ever since, including the Hungarian migration to Slavonia in 10 th century, and the influx of Islamic and Orthodox Balkan and Asian populations during the Ottoman conquest in 16 th century, causing at the same time, the continuous shift of Catholics from Bosnia to Slavonia during several centuries 2. Under the auspices of Habsburg monarchy, the settlement of Germans and Austrians in Slavonian urban areas peaks between 18 th and 19 th century, while Orthodox Vlachs from Bosnia, immigrating Czechs, Slovaks, Ukrainians, Italians, and

No abstract

Idiopathic pulmonary fibrosis (IPF) affects lung parenchyma with progressing fibrosis. In this study, we aimed to replicate MUC5B rs35705950 variants and determine new plausible candidate variants for IPF among four different European populations. We genotyped 26 IPF candidate loci in 165 IPF patients from four European countries, such as Czech Republic (n = 41), Germany (n = 33), Greece (n = 40), France (n = 51), and performed association study comparing observed variant distribution with that obtained in a genetically similar Czech healthy control population (n = 96) described in our earlier data report. A highly significant association for a promoter variant (rs35705950) of mucin encoding MUC5B gene was observed in all IPF populations, individually and combined [odds ratio (95% confidence interval); p-value as 5.23 (8.94–3.06); 1.80 × 10−11]. Another non-coding variant, rs7934606 in MUC2 was significant among German patients [2.85 (5.05–1.60); 4.03 × 10−4] and combined European IPF cases [2.18 (3.16–1.50); 3.73 × 10−5]. The network analysis for these variants indicated gene–gene and gene–phenotype interactions in IPF and lung biology. With replication of MUC5B rs35705950 previously reported in U.S. populations of European descent and indicating other plausible polymorphic variants relevant for IPF, we provide additional reference information for future extended functional and population studies aimed, ideally with inclusion of clinical parameters, at identification of IPF genetic markers.

No abstract