Chunhua Yan

Genetic Alterations Activating Kinase and Cytokine Receptor Signaling in High-Risk Acute Lymphoblastic Leukemia

Genomic profiling has identified a subtype of high-risk B-progenitor acute lymphoblastic leukemia (B-ALL) with alteration of IKZF1, a gene expression profile similar to BCR-ABL1-positive ALL and poor outcome (Ph-like ALL). The genetic alterations that activate kinase signaling in Ph-like ALL are poorly understood. We performed transcriptome and whole genome sequencing on 15 cases of Ph-like ALL and identified rearrangements involving ABL1, JAK2, PDGFRB, CRLF2, and EPOR, activating mutations of IL7R and FLT3, and deletion of SH2B3, which encodes the JAK2-negative regulator LNK. Importantly, several of these alterations induce transformation that is attenuated with tyrosine kinase inhibitors, suggesting the treatment outcome of these patients may be improved with targeted therapy.

Genetic variations at loci involved in the immune response are risk factors for hepatocellular carcinoma.

Primary liver cancer is the third most common cause of cancer‐related death worldwide, with a rising incidence in Western countries. Little is known about the genetic etiology of this disease. To identify genetic factors associated with hepatocellular carcinoma (HCC) and liver cirrhosis (LC), we conducted a comprehensive, genome‐wide variation analysis in a population of unrelated Asian individuals. Copy number variation (CNV) and single nucleotide polymorphisms (SNPs) were assayed in peripheral blood with the high‐density Affymetrix SNP6.0 microarray platform. We used a two‐stage discovery and replication design to control for overfitting and to validate observed results. We identified a strong association with CNV at the T‐cell receptor gamma and alpha loci (P < 1 × 10−15) in HCC cases when contrasted with controls. This variation appears to be somatic in origin, reflecting differences between T‐cell receptor processing in lymphocytes from individuals with liver disease and healthy individuals that is not attributable to chronic hepatitis virus infection. Analysis of constitutional variation identified three susceptibility loci including the class II MHC complex, whose protein products present antigen to T‐cell receptors and mediate immune surveillance. Statistical analysis of biologic networks identified variation in the “antigen presentation and processing” pathway as being highly significantly associated with HCC (P = 1 × 10−11). SNP analysis identified two variants whose allele frequencies differ significantly between HCC and LC. One of these (P = 1.74 × 10−12) lies in the PTEN homolog TPTE2. Conclusion: Combined analysis of CNV, individual SNPs, and pathways suggest that HCC susceptibility is mediated by germline factors affecting the immune response and differences in T‐cell receptor processing. (HEPATOLOGY 2010)

Recurrent DGCR8, DROSHA, and SIX Homeodomain Mutations in Favorable Histology Wilms Tumors

We report the most common single-nucleotide substitution/deletion mutations in favorable histology Wilms tumors (FHWTs) to occur within SIX1/2 (7% of 534 tumors) and microRNA processing genes (miRNAPGs) DGCR8 and DROSHA (15% of 534 tumors). Comprehensive analysis of 77 FHWTs indicates that tumors with SIX1/2 and/or miRNAPG mutations show a pre-induction metanephric mesenchyme gene expression pattern and are significantly associated with both perilobar nephrogenic rests and 11p15 imprinting aberrations. Significantly decreased expression of mature Let-7a and the miR-200 family (responsible for mesenchymal-to-epithelial transition) in miRNAPG mutant tumors is associated with an undifferentiated blastemal histology. The combination of SIX and miRNAPG mutations in the same tumor is associated with evidence of RAS activation and a higher rate of relapse and death.

OmicCircos: A Simple-to-Use R Package for the Circular Visualization of Multidimensional Omics Data

Summary OmicCircos is an R software package used to generate high-quality circular plots for visualizing genomic variations, including mutation patterns, copy number variations (CNVs), expression patterns, and methylation patterns. Such variations can be displayed as scatterplot, line, or text-label figures. Relationships among genomic features in different chromosome positions can be represented in the forms of polygons or curves. Utilizing the statistical and graphic functions in an R/Bioconductor environment, OmicCircos performs statistical analyses and displays results using cluster, boxplot, histogram, and heatmap formats. In addition, OmicCircos offers a number of unique capabilities, including independent track drawing for easy modification and integration, zoom functions, link-polygons, and position-independent heatmaps supporting detailed visualization. Availability and Implementation OmicCircos is available through Bioconductor at http://www.bioconductor.org/packages/devel/bioc/html/OmicCircos.html. An extensive vignette in the package describes installation, data formatting, and workflow procedures. The software is open source under the Artistic—2.0 license.

A Children's Oncology Group and TARGET initiative exploring the genetic landscape of Wilms tumor

We performed genome-wide sequencing and analyzed mRNA and miRNA expression, DNA copy number, and DNA methylation in 117 Wilms tumors, followed by targeted sequencing of 651 Wilms tumors. In addition to genes previously implicated in Wilms tumors (WT1, CTNNB1, AMER1, DROSHA, DGCR8, XPO5, DICER1, SIX1, SIX2, MLLT1, MYCN, and TP53), we identified mutations in genes not previously recognized as recurrently involved in Wilms tumors, the most frequent being BCOR, BCORL1, NONO, MAX, COL6A3, ASXL1, MAP3K4, and ARID1A. DNA copy number changes resulted in recurrent 1q gain, MYCN amplification, LIN28B gain, and MIRLET7A loss. Unexpected germline variants involved PALB2 and CHEK2. Integrated analyses support two major classes of genetic changes that preserve the progenitor state and/or interrupt normal development.

Discovery Analysis of TCGA Data Reveals Association between Germline Genotype and Survival in Ovarian Cancer Patients

BACKGROUND Ovarian cancer remains a significant public health burden, with the highest mortality rate of all the gynecological cancers. This is attributable to the late stage at which the majority of ovarian cancers are diagnosed, coupled with the low and variable response of advanced tumors to standard chemotherapies. To date, clinically useful predictors of treatment response remain lacking. Identifying the genetic determinants of ovarian cancer survival and treatment response is crucial to the development of prognostic biomarkers and personalized therapies that may improve outcomes for the late-stage patients who comprise the majority of cases. METHODS To identify constitutional genetic variations contributing to ovarian cancer mortality, we systematically investigated associations between germline polymorphisms and ovarian cancer survival using data from The Cancer Genome Atlas Project (TCGA). Using stage-stratified Cox proportional hazards regression, we examined >650,000 SNP loci for association with survival. We additionally examined whether the association of significant SNPs with survival was modified by somatic alterations. RESULTS Germline polymorphisms at rs4934282 (AGAP11/C10orf116) and rs1857623 (DNAH14) were associated with stage-adjusted survival (p= 1.12e-07 and 1.80e-07, FDR q= 1.2e-04 and 2.4e-04, respectively). A third SNP, rs4869 (C10orf116), was additionally identified as significant in the exome sequencing data; it is in near-perfect LD with rs4934282. The associations with survival remained significant when somatic alterations. CONCLUSIONS Discovery analysis of TCGA data reveals germline genetic variations that may play a role in ovarian cancer survival even among late-stage cases. The significant loci are located near genes previously reported as having a possible relationship to platinum and taxol response. Because the variant alleles at the significant loci are common (frequencies for rs4934282 A/C alleles = 0.54/0.46, respectively; rs1857623 A/G alleles = 0.55/0.45, respectively) and germline variants can be assayed noninvasively, our findings provide potential targets for further exploration as prognostic biomarkers and individualized therapies.

Systematic genetic analysis identifies Cis-eQTL target genes associated with glioblastoma patient survival

Prior expression quantitative trait locus (eQTL) studies have demonstrated heritable variation determining differences in gene expression. The majority of eQTL studies were based on cell lines and normal tissues. We performed cis-eQTL analysis using glioblastoma multiforme (GBM) data sets obtained from The Cancer Genome Atlas (TCGA) to systematically investigate germline variation’s contribution to tumor gene expression levels. We identified 985 significant cis-eQTL associations (FDR<0.05) mapped to 978 SNP loci and 159 unique genes. Approximately 57% of these eQTLs have been previously linked to the gene expression in cell lines and normal tissues; 43% of these share cis associations known to be associated with functional annotations. About 25% of these cis-eQTL associations are also common to those identified in Breast Cancer from a recent study. Further investigation of the relationship between gene expression and patient clinical information identified 13 eQTL genes whose expression level significantly correlates with GBM patient survival (p<0.05). Most of these genes are also differentially expressed in tumor samples and organ-specific controls (p<0.05). Our results demonstrated a significant relationship of germline variation with gene expression levels in GBM. The identification of eQTLs-based expression associated survival might be important to the understanding of genetic contribution to GBM cancer prognosis.

Genomic Landscape of Somatic Alterations in Esophageal Squamous Cell Carcinoma and Gastric Cancer

Gastric cancer and esophageal cancer are the second and sixth leading causes of cancer-related death worldwide. Multiple genomic alterations underlying gastric cancer and esophageal squamous cell carcinoma (ESCC) have been identified, but the full spectrum of genomic structural variations and mutations have yet to be uncovered. Here, we report the results of whole-genome sequencing of 30 samples comprising tumor and blood from 15 patients, four of whom presented with ESCC, seven with gastric cardia adenocarcinoma (GCA), and four with gastric noncardia adenocarcinoma. Analyses revealed that an A>C mutation was common in GCA, and in addition to the preferential nucleotide sequence of A located 5 prime to the mutation as noted in previous studies, we found enrichment of T in the 5 prime base. The A>C mutations in GCA suggested that oxidation of guanine may be a potential mechanism underlying cancer mutagenesis. Furthermore, we identified genes with mutations in gastric cancer and ESCC, including well-known cancer genes, TP53, JAK3, BRCA2, FGF2, FBXW7, MSH3, PTCH, NF1, ERBB2, and CHEK2, and potentially novel cancer-associated genes, KISS1R, AMH, MNX1, WNK2, and PRKRIR Finally, we identified recurrent chromosome alterations in at least 30% of tumors in genes, including MACROD2, FHIT, and PARK2 that were often intragenic deletions. These structural alterations were validated using the The Cancer Genome Atlas dataset. Our studies provide new insights into understanding the genomic landscape, genome instability, and mutation profile underlying gastric cancer and ESCC development. Cancer Res; 76(7); 1714-23. ©2016 AACR.

Multi-SNP analysis of GWAS data identifies pathways associated with nonalcoholic fatty liver disease.

Non-alcoholic fatty liver disease (NAFLD) is a common liver disease; the histological spectrum of which ranges from steatosis to steatohepatitis. Nonalcoholic steatohepatitis (NASH) often leads to cirrhosis and development of hepatocellular carcinoma. To better understand pathogenesis of NAFLD, we performed the pathway of distinction analysis (PoDA) on a genome-wide association study dataset of 250 non-Hispanic white female adult patients with NAFLD, who were enrolled in the NASH Clinical Research Network (CRN) Database Study, to investigate whether biologic process variation measured through genomic variation of genes within these pathways was related to the development of steatohepatitis or cirrhosis. Pathways such as Recycling of eIF2:GDP, biosynthesis of steroids, Terpenoid biosynthesis and Cholesterol biosynthesis were found to be significantly associated with NASH. SNP variants in Terpenoid synthesis, Cholesterol biosynthesis and biosynthesis of steroids were associated with lobular inflammation and cytologic ballooning while those in Terpenoid synthesis were also associated with fibrosis and cirrhosis. These were also related to the NAFLD activity score (NAS) which is derived from the histological severity of steatosis, inflammation and ballooning degeneration. Eukaryotic protein translation and recycling of eIF2:GDP related SNP variants were associated with ballooning, steatohepatitis and cirrhosis. Il2 signaling events mediated by PI3K, Mitotic metaphase/anaphase transition, and Prostanoid ligand receptors were also significantly associated with cirrhosis. Taken together, the results provide evidence for additional ways, beyond the effects of single SNPs, by which genetic factors might contribute to the susceptibility to develop a particular phenotype of NAFLD and then progress to cirrhosis. Further studies are warranted to explain potential important genetic roles of these biological processes in NAFLD.

Significance of TP53 mutation in Wilms tumors with diffuse anaplasia : A report from the Children's Oncology Group

Purpose: To investigate the role and significance of TP53 mutation in diffusely anaplastic Wilms tumors (DAWTs). Experimental Design: All DAWTs registered on National Wilms Tumor Study-5 (n = 118) with available samples were analyzed for TP53 mutations and copy loss. Integrative genomic analysis was performed on 39 selected DAWTs. Results: Following analysis of a single random sample, 57 DAWTs (48%) demonstrated TP53 mutations, 13 (11%) copy loss without mutation, and 48 (41%) lacked both [defined as TP53-wild-type (wt)]. Patients with stage III/IV TP53-wt DAWTs (but not those with stage I/II disease) had significantly lower relapse and death rates than those with TP53 abnormalities. In-depth analysis of a subset of 39 DAWTs showed seven (18%) to be TP53-wt: These demonstrated gene expression evidence of an active p53 pathway. Retrospective pathology review of TP53-wt DAWT revealed no or very low volume of anaplasia in six of seven tumors. When samples from TP53-wt tumors known to contain anaplasia histologically were available, abnormal p53 protein accumulation was observed by immunohistochemistry. Conclusions: These data support the key role of TP53 loss in the development of anaplasia in WT, and support its significant clinical impact in patients with residual anaplastic tumor following surgery. These data also suggest that most DAWTs will show evidence of TP53 mutation when samples selected for the presence of anaplasia are analyzed. This suggests that modifications of the current criteria to also consider volume of anaplasia and documentation of TP53 aberrations may better reflect the risk of relapse and death and enable optimization of therapeutic stratification. Clin Cancer Res; 22(22); 5582–91. ©2016 AACR.