An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia.

Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

Observations that suggest a contribution of altered dermal papilla mitochondrial function to androgenetic alopecia.

Androgenetic alopecia (AGA) is a prevalent hair loss condition in males that develops due to the influence of androgens and genetic predisposition. With the aim of elucidating genes involved in AGA pathogenesis, we modelled AGA with three-dimensional culture of keratinocyte-surrounded dermal papilla (DP) cells. We co-cultured immortalised balding and non-balding human DP cells (DPCs) derived from male AGA patients with epidermal keratinocyte (NHEK) using multi-interfacial polyelectrolyte complexation technique. We observed up-regulated mitochondria-related gene expression in balding compared with non-balding DP aggregates which indicated altered mitochondria metabolism. Further observation of significantly reduced electron transport chain complex activity (complexes I, IV and V), ATP levels and ability to uptake metabolites for ATP generation demonstrated compromised mitochondria function in balding DPC. Balding DP was also found to be under significantly higher oxidative stress than non-balding DP. Our experiments suggest that application of antioxidants lowers oxidative stress levels and improves metabolite uptake in balding DPC. We postulate that the observed up-regulation of mitochondria-related genes in balding DP aggregates resulted from an over-compensatory effort to rescue decreased mitochondrial function in balding DP through the attempted production of new functional mitochondria. In all, our three-dimensional co-culturing revealed mitochondrial dysfunction in balding DPC, suggesting a metabolic component in the aetiology of AGA.

Shortened Telomere Length in Sputum Cells of Bronchiectasis Patients is Associated with Dysfunctional Inflammatory Pathways.

Telomere attrition is an established ageing biomarker and shorter peripheral blood leukocyte telomere length has been associated with increased risks of respiratory diseases. However, whether telomere length in disease-relevant sputum immune cells of chronic respiratory disease patients is shortened and which pathways are dysfunctional are not clear. Here we measured telomere length from sputum samples of bronchiectasis and asthmatic subjects and determined that telomere length in sputum of bronchiectasis subjects was significantly shorter (Beta = - 1.167, PAdj = 2.75 × 10-4). We further performed global gene expression analysis and identified genes involved in processes such as NLRP3 inflammasome activation and regulation of adaptive immune cells when bronchiectasis sputum telomere length was shortened. Our study provides insights on dysfunctions related to shortened telomere length in sputum immune cells of bronchiectasis patients.

Patient-specific driver gene prediction and risk assessment through integrated network analysis of cancer omics profiles.

Extensive and multi-dimensional data sets generated from recent cancer omics profiling projects have presented new challenges and opportunities for unraveling the complexity of cancer genome landscapes. In particular, distinguishing the unique complement of genes that drive tumorigenesis in each patient from a sea of passenger mutations is necessary for translating the full benefit of cancer genome sequencing into the clinic. We address this need by presenting a data integration framework (OncoIMPACT) to nominate patient-specific driver genes based on their phenotypic impact. Extensive in silico and in vitro validation helped establish OncoIMPACT's robustness, improved precision over competing approaches and verifiable patient and cell line specific predictions (2/2 and 6/7 true positives and negatives, respectively). In particular, we computationally predicted and experimentally validated the gene TRIM24 as a putative novel amplified driver in a melanoma patient. Applying OncoIMPACT to more than 1000 tumor samples, we generated patient-specific driver gene lists in five different cancer types to identify modes of synergistic action. We also provide the first demonstration that computationally derived driver mutation signatures can be overall superior to single gene and gene expression based signatures in enabling patient stratification and prognostication. Source code and executables for OncoIMPACT are freely available from http://sourceforge.net/projects/oncoimpact.

Comprehensive long-span paired-end-tag mapping reveals characteristic patterns of structural variations in epithelial cancer genomes.

Somatic genome rearrangements are thought to play important roles in cancer development. We optimized a long-span paired-end-tag (PET) sequencing approach using 10-Kb genomic DNA inserts to study human genome structural variations (SVs). The use of a 10-Kb insert size allows the identification of breakpoints within repetitive or homology-containing regions of a few kilobases in size and results in a higher physical coverage compared with small insert libraries with the same sequencing effort. We have applied this approach to comprehensively characterize the SVs of 15 cancer and two noncancer genomes and used a filtering approach to strongly enrich for somatic SVs in the cancer genomes. Our analyses revealed that most inversions, deletions, and insertions are germ-line SVs, whereas tandem duplications, unpaired inversions, interchromosomal translocations, and complex rearrangements are over-represented among somatic rearrangements in cancer genomes. We demonstrate that the quantitative and connective nature of DNA-PET data is precise in delineating the genealogy of complex rearrangement events, we observe signatures that are compatible with breakage-fusion-bridge cycles, and we discover that large duplications are among the initial rearrangements that trigger genome instability for extensive amplification in epithelial cancers.

ChIA-PET analysis of transcriptional chromatin interactions.

Long-range chromatin contacts between specific DNA regulatory elements play a pivotal role in gene expression regulation, and a global characterization of these interactions in the 3-dimensional (3D) chromatin structure is imperative in understanding signaling networks and cell states. Chromatin Interaction Analysis using Paired-End Tag sequencing (ChIA-PET) is a method which converts functional chromatin structure into millions of short tag sequences. Combining Chromatin Immunoprecipitation (ChIP), proximity ligation and high-throughput sequencing, ChIA-PET provides a global and unbiased interrogation of higher-order chromatin structures associated with specific protein factors. Here, we describe the detailed procedures of the ChIA-PET methodology, unraveling transcription-associated chromatin contacts in a model human cell line.

Experimental and bioinformatics considerations in cancer application of single cell genomics.

Single cell genomics offers an unprecedented resolution to interrogate genetic heterogeneity in a patient's tumour at the intercellular level. However, the DNA yield per cell is insufficient for today's sequencing library preparation protocols. This necessitates DNA amplification which is a key source of experimental noise. We provide an evaluation of two protocols using micro-fluidics based amplification for whole exome sequencing, which is an experimental scenario commonly used in single cell genomics. The results highlight their respective biases and relative strengths in identification of single nucleotide variations. Towards this end, we introduce a workflow SoVaTSiC, which allows for quality evaluation and somatic variant identification of single cell data. As proof of concept, the framework was applied to study a lung adenocarcinoma tumour. The analysis provides insights into tumour phylogeny by identifying key mutational events in lung adenocarcinoma evolution. The consequence of this inference is supported by the histology of the tumour and demonstrates usefulness of the approach.

Differential Expression between Human Dermal Papilla Cells from Balding and Non-Balding Scalps Reveals New Candidate Genes for Androgenetic Alopecia.

Androgenetic alopecia (AGA) is a common heritable and androgen-dependent hair loss condition in men. Twelve genetic risk loci are known to date, but it is unclear which genes at these loci are relevant for AGA. Dermal papilla cells (DPCs) located in the hair bulb are the main site of androgen activity in the hair follicle. Widely used monolayer-cultured primary DPCs in hair-related studies often lack dermal papilla characteristics. In contrast, immortalized DPCs have high resemblance to intact dermal papilla. We derived immortalized human DPC lines from balding (BAB) and non-balding (BAN) scalp. Both BAB and BAN retained high proportions of dermal papilla signature gene and versican protein expression. We performed expression analysis of BAB and BAN and annotated AGA risk loci with differentially expressed genes. We found evidence for AR but not EDA2R as the candidate gene at the AGA risk locus on chromosome X. Further, our data suggest TWIST1 (twist family basic helix-loop-helix transcription factor 1) and SSPN (sarcospan) to be the functionally relevant AGA genes at the 7p21.1 and 12p12.1 risk loci, respectively. Down-regulated genes in BAB compared to BAN were highly enriched for vasculature-related genes, suggesting that deficiency of DPC from balding scalps in fostering vascularization around the hair follicle may contribute to the development of AGA.

Recurrent Fusion Genes in Gastric Cancer: CLDN18-ARHGAP26 Induces Loss of Epithelial Integrity.

Genome rearrangements, a hallmark of cancer, can result in gene fusions with oncogenic properties. Using DNA paired-end-tag (DNA-PET) whole-genome sequencing, we analyzed 15 gastric cancers (GCs) from Southeast Asians. Rearrangements were enriched in open chromatin and shaped by chromatin structure. We identified seven rearrangement hot spots and 136 gene fusions. In three out of 100 GC cases, we found recurrent fusions between CLDN18, a tight junction gene, and ARHGAP26, a gene encoding a RHOA inhibitor. Epithelial cell lines expressing CLDN18-ARHGAP26 displayed a dramatic loss of epithelial phenotype and long protrusions indicative of epithelial-mesenchymal transition (EMT). Fusion-positive cell lines showed impaired barrier properties, reduced cell-cell and cell-extracellular matrix adhesion, retarded wound healing, and inhibition of RHOA. Gain of invasion was seen in cancer cell lines expressing the fusion. Thus, CLDN18-ARHGAP26 mediates epithelial disintegration, possibly leading to stomach H(+) leakage, and the fusion might contribute to invasiveness once a cell is transformed.

Large-scale functional organization of long-range chromatin interaction networks.

Chromatin interactions play important roles in transcription regulation. To better understand the underlying evolutionary and functional constraints of these interactions, we implemented a systems approach to examine RNA polymerase-II-associated chromatin interactions in human cells. We found that 40% of the total genomic elements involved in chromatin interactions converged to a giant, scale-free-like, hierarchical network organized into chromatin communities. The communities were enriched in specific functions and were syntenic through evolution. Disease-associated SNPs from genome-wide association studies were enriched among the nodes with fewer interactions, implying their selection against deleterious interactions by limiting the total number of interactions, a model that we further reconciled using somatic and germline cancer mutation data. The hubs lacked disease-associated SNPs, constituted a nonrandomly interconnected core of key cellular functions, and exhibited lethality in mouse mutants, supporting an evolutionary selection that favored the nonrandom spatial clustering of the least-evolving key genomic domains against random genetic or transcriptional errors in the genome. Altogether, our analyses reveal a systems-level evolutionary framework that shapes functionally compartmentalized and error-tolerant transcriptional regulation of human genome in three dimensions.

Whole-genome reconstruction and mutational signatures in gastric cancer.

BACKGROUND: Gastric cancer is the second highest cause of global cancer mortality. To explore the complete repertoire of somatic alterations in gastric cancer, we combined massively parallel short read and DNA paired-end tag sequencing to present the first whole-genome analysis of two gastric adenocarcinomas, one with chromosomal instability and the other with microsatellite instability. RESULTS: Integrative analysis and de novo assemblies revealed the architecture of a wild-type KRAS amplification, a common driver event in gastric cancer. We discovered three distinct mutational signatures in gastric cancer--against a genome-wide backdrop of oxidative and microsatellite instability-related mutational signatures, we identified the first exome-specific mutational signature. Further characterization of the impact of these signatures by combining sequencing data from 40 complete gastric cancer exomes and targeted screening of an additional 94 independent gastric tumors uncovered ACVR2A, RPL22 and LMAN1 as recurrently mutated genes in microsatellite instability-positive gastric cancer and PAPPA as a recurrently mutated gene in TP53 wild-type gastric cancer. CONCLUSIONS: These results highlight how whole-genome cancer sequencing can uncover information relevant to tissue-specific carcinogenesis that would otherwise be missed from exome-sequencing data.

Effects of the demethylating agent, 5-azacytidine, on expression of the kallikrein-kinin genes in carcinoma cells of the lung and pleura.

Tissue kallikrein (KLK1) and plasma kallikrein (KLKB1) may regulate the growth and proliferation of tumours of the lung and pleura, through the generation of kinin peptides that signal through the kinin B(1) (BDKRB1) and B(2) (BDKRB2) receptors. The development and progression of cancer results from genetic mutations, as well as epigenetic changes that include methylation of DNA at CpG islands. The aim of this study was to assess whether expression of the kallikrein-kinin genes in lung cancer and mesothelioma cells is regulated by DNA methylation. Quantitative reverse transcriptase-PCR and immunocytochemistry showed differences in the basal expression of the kallikrein-kinin genes and proteins in lung carcinoma and mesothelioma cells, compared with non-malignant lung epithelial and mesothelial cells, respectively. Following treatment with the demethylating agent, 5-azacytidine (5-AZA), KLKB1 mRNA expression was consistently increased in both lung carcinoma and mesothelioma cells, whereas KLK1, BDKRB1 and BDKRB2 mRNA expression was decreased or unchanged. Increased expression of KLKB1 after 5-AZA treatment suggests it may function as a tumour suppressor gene in cancers of the lung and pleura. Studies on DNA methylation of the kallikrein-kinin genes will enhance understanding of their role in carcinogenesis and provide insights into the importance of kallikreins as tumour biomarkers.