Weiyin Zhou

Epidermal growth factor receptor polymorphisms and clinical outcomes in non-small-cell lung cancer patients treated with gefitinib

The−216G/T, −191C/A, intron 1 and Arg497Lys epidermal growth factor receptor (EGFR) polymorphisms were evaluated in 92 advanced non-small-cell lung cancer patients treated with gefitinib, an EGFR tyrosine-kinase inhibitor. Improved progression free survival (PFS) was found in patients homozygous for the shorter lengths of intron 1 polymorphism (S/S; S=16 or fewer CA repeats; log-rank test (LRT) P=0.03) and for patients carrying any T allele of the −216G/T polymorphism (LRT, P=0.005). When considered together, patients with intron 1 S/S genotype and at least one T allele of −216G/T had improved PFS (LRT P=0.0006; adjusted hazard ratio (AHR), 0.60 (95% confidence interval, 0.36–0.98)) and overall survival (LRT P=0.02; AHR, 0.60 (0.36–1.00)) when compared with all others. The T allele of −216G/T was also associated with significantly higher rates of stable disease/partial response (P=0.01) and a significantly higher risk of treatment-related rash/diarrhea (P=0.004, multivariate model). EGFR intron 1 and –216G/T polymorphisms influence clinical outcomes in gefitinib-treated non-small-cell lung cancer patients.

Characterizing the genetic basis of methylome diversity in histologically normal human lung tissue

The genetic regulation of the human epigenome is not fully appreciated. Here we describe the effects of genetic variants on the DNA methylome in human lung based on methylation-quantitative trait loci (meQTL) analyses. We report 34,304 cis- and 585 trans-meQTLs, a genetic-epigenetic interaction of surprising magnitude, including a regulatory hotspot. These findings are replicated in both breast and kidney tissues and show distinct patterns: cis-meQTLs mostly localize to CpG sites outside of genes, promoters, and CpG islands (CGIs), while trans-meQTLs are over-represented in promoter CGIs. meQTL SNPs are enriched in CTCF binding sites, DNaseI hypersensitivity regions and histone marks. Importantly, 4 of the 5 established lung cancer risk loci in European ancestry are cis-meQTLs and, in aggregate, cis-meQTLs are enriched for lung cancer risk in a genome-wide analysis of 11,587 subjects. Thus, inherited genetic variation may affect lung carcinogenesis by regulating the human methylome.

Mosaic loss of chromosome Y is associated with common variation near TCL1A

Mosaic loss of chromosome Y (mLOY) leading to gonosomal XY/XO commonly occurs during aging, particularly in smokers. We investigated whether mLOY was associated with non-hematological cancer in three prospective cohorts (8,679 cancer cases and 5,110 cancer-free controls) and genetic susceptibility to mLOY. Overall, mLOY was observed in 7% of men, and its prevalence increased with age (per-year odds ratio (OR) = 1.13, 95% confidence interval (CI) = 1.12–1.15; P < 2 × 10−16), reaching 18.7% among men over 80 years old. mLOY was associated with current smoking (OR = 2.35, 95% CI = 1.82–3.03; P = 5.55 × 10−11), but the association weakened with years after cessation. mLOY was not consistently associated with overall or specific cancer risk (for example, bladder, lung or prostate cancer) nor with cancer survival after diagnosis (multivariate-adjusted hazard ratio = 0.87, 95% CI = 0.73–1.04; P = 0.12). In a genome-wide association study, we observed the first example of a common susceptibility locus for genetic mosaicism, specifically mLOY, which maps to TCL1A at 14q32.13, marked by rs2887399 (OR = 1.55, 95% CI = 1.36–1.78; P = 1.37 × 10−10).

An integrated transcriptome and epigenome analysis identifies a novel candidate gene for pancreatic cancer

BackgroundPancreatic cancer is a highly lethal cancer with limited diagnostic and therapeutic modalities.MethodsTo begin to explore the genomic landscape of pancreatic cancer, we used massively parallel sequencing to catalog and compare transcribed regions and potential regulatory elements in two human cell lines derived from normal and cancerous pancreas.ResultsBy RNA-sequencing, we identified 2,146 differentially expressed genes in these cell lines that were enriched in cancer related pathways and biological processes that include cell adhesion, growth factor and receptor activity, signaling, transcription and differentiation. Our high throughput Chromatin immunoprecipitation (ChIP) sequence analysis furthermore identified over 100,000 regions enriched in epigenetic marks, showing either positive (H3K4me1, H3K4me3, RNA Pol II) or negative (H3K27me3) correlation with gene expression. Notably, an overall enrichment of RNA Pol II binding and depletion of H3K27me3 binding were seen in the cancer derived cell line as compared to the normal derived cell line. By selecting genes for further assessment based on this difference, we confirmed enhanced expression of aldehyde dehydrogenase 1A3 (ALDH1A3) in two larger sets of pancreatic cancer cell lines and in tumor tissues as compared to normal derived tissues.ConclusionsAs aldehyde dehydrogenase (ALDH) activity is a key feature of cancer stem cells, our results indicate that a member of the ALDH superfamily, ALDH1A3, may be upregulated in pancreatic cancer, where it could mark pancreatic cancer stem cells.

Novel and known ribosomal causes of Diamond-Blackfan anaemia identified through comprehensive genomic characterisation

Background Diamond-Blackfan anaemia (DBA) is an inherited bone marrow failure syndrome (IBMFS) characterised by erythroid hypoplasia. It is associated with congenital anomalies and a high risk of developing specific cancers. DBA is caused predominantly by autosomal dominant pathogenic variants in at least 15 genes affecting ribosomal biogenesis and function. Two X-linked recessive genes have been identified. Objectives We aim to identify the genetic aetiology of DBA. Methods Of 87 families with DBA enrolled in an institutional review board-approved cohort study (ClinicalTrials.gov Identifier:NCT00027274), 61 had genetic testing information available. Thirty-five families did not have a known genetic cause and thus underwent comprehensive genomic evaluation with whole exome sequencing, deletion and CNV analyses to identify their disease-associated pathogenic variant. Controls for functional studies were healthy mutation-negative individuals enrolled in the same study. Results Our analyses uncovered heterozygous pathogenic variants in two previously undescribed genes in two families. One family had a non-synonymous variant (p.K77N) in RPL35; the second family had a non-synonymous variant (p. L51S) in RPL18. Both of these variants result in pre-rRNA processing defects. We identified heterozygous pathogenic variants in previously known DBA genes in 16 of 35 families. Seventeen families who underwent genetic analyses are yet to have a genetic cause of disease identified. Conclusions Overall, heterozygous pathogenic variants in ribosomal genes were identified in 44 of the 61 families (72%). De novo pathogenic variants were observed in 57% of patients with DBA. Ongoing studies of DBA genomics will be important to understand this complex disorder.

Genome Analysis of Latin American Cervical Cancer: Frequent Activation of the PIK3CA Pathway

Purpose: Cervical cancer is one of the most common causes of cancer mortality for women living in poverty, causing more than 28,000 deaths annually in Latin America and 266,000 worldwide. To better understand the molecular basis of the disease, we ascertained blood and tumor samples from Guatemala and Venezuela and performed genomic characterization. Experimental Design: We performed human papillomavirus (HPV) typing and identified somatically mutated genes using exome and ultra-deep targeted sequencing with confirmation in samples from Mexico. Copy number changes were also assessed in the exome sequence. Results: Cervical cancer cases in Guatemala and Venezuela have an average age of diagnosis of 50 years and 5.6 children. Analysis of 675 tumors revealed activation of PIK3CA and other PI3K/AKT pathway genes in 31% of squamous carcinomas and 24% of adeno- and adenosquamous tumors, predominantly at two sites (E542K, E545K) in the helical domain of the PIK3CA gene. This distribution of PIK3CA mutations is distinct from most other cancer types and does not result in the in vitro phosphorylation of AKT. Somatic mutations were more frequent in squamous carcinomas diagnosed after the age of 50 years. Frequent gain of chromosome 3q was found, and low PIK3CA mutation fractions in many tumors suggest that PI3K mutation can be a late event in tumor progression. Conclusions: PI3K pathway mutation is important to cervical carcinogenesis in Latin America. Therapeutic agents that directly target PI3K could play a role in the therapy of this common malignancy. Clin Cancer Res; 21(23); 5360–70. ©2015 AACR.

Genome-wide association study confirms lung cancer susceptibility loci on chromosomes 5p15 and 15q25 in an African-American population

OBJECTIVES Genome-wide association studies (GWAS) of lung cancer have identified regions of common genetic variation with lung cancer risk in Europeans who smoke and never-smoking Asian women. This study aimed to conduct a GWAS in African Americans, who have higher rates of lung cancer despite smoking fewer cigarettes per day when compared with Caucasians. This population provides a different genetic architecture based on underlying African ancestry allowing the identification of new regions and exploration of known regions for finer mapping. MATERIALS AND METHODS We genotyped 1,024,001 SNPs in 1737 cases and 3602 controls in stage 1, followed by a replication phase of 20 SNPs (p<1.51×10(-5)) in an independent set of 866 cases and 796 controls in stage 2. RESULTS AND CONCLUSION In the combined analysis, we confirmed two loci to be associated with lung cancer that achieved the threshold of genome-wide significance: 15q25.1 marked by rs2036527 (p=1.3×10(-9); OR=1.32; 95% CI=1.20-1.44) near CHRNA5, and 5p15.33 marked by rs2853677 (p=2.8×10(-9); OR=1.28; 95% CI=1.18-1.39) near TERT. The association with rs2853677 is driven by the adenocarcinoma subtype of lung cancer (p=1.3×10(-8); OR=1.37; 95% CI=1.23-1.54). No SNPs reached genome-wide significance for either of the main effect models examining smoking - cigarettes per day and current or former smoker. Our study was powered to identify strong risk loci for lung cancer in African Americans; we confirmed results previously reported in African Americans and other populations for two loci near plausible candidate genes, CHRNA5 and TERT, on 15q25.1 and 5p15.33 respectively, are associated with lung cancer. Additional work is required to map and understand the biological underpinnings of the strong association of these loci with lung cancer risk in African Americans.

Somatic Genomics and Clinical Features of Lung Adenocarcinoma: A Retrospective Study

Background Lung adenocarcinoma (LUAD) is the most common histologic subtype of lung cancer and has a high risk of distant metastasis at every disease stage. We aimed to characterize the genomic landscape of LUAD and identify mutation signatures associated with tumor progression. Methods and Findings We performed an integrative genomic analysis, incorporating whole exome sequencing (WES), determination of DNA copy number and DNA methylation, and transcriptome sequencing for 101 LUAD samples from the Environment And Genetics in Lung cancer Etiology (EAGLE) study. We detected driver genes by testing whether the nonsynonymous mutation rate was significantly higher than the background mutation rate and replicated our findings in public datasets with 724 samples. We performed subclonality analysis for mutations based on mutant allele data and copy number alteration data. We also tested the association between mutation signatures and clinical outcomes, including distant metastasis, survival, and tumor grade. We identified and replicated two novel candidate driver genes, POU class 4 homeobox 2 (POU4F2) (mutated in 9 [8.9%] samples) and ZKSCAN1 (mutated in 6 [5.9%] samples), and characterized their major deleterious mutations. ZKSCAN1 was part of a mutually exclusive gene set that included the RTK/RAS/RAF pathway genes BRAF, EGFR, KRAS, MET, and NF1, indicating an important driver role for this gene. Moreover, we observed strong associations between methylation in specific genomic regions and somatic mutation patterns. In the tumor evolution analysis, four driver genes had a significantly lower fraction of subclonal mutations (FSM), including TP53 (p = 0.007), KEAP1 (p = 0.012), STK11 (p = 0.0076), and EGFR (p = 0.0078), suggesting a tumor initiation role for these genes. Subclonal mutations were significantly enriched in APOBEC-related signatures (p < 2.5×10−50). The total number of somatic mutations (p = 0.0039) and the fraction of transitions (p = 5.5×10−4) were associated with increased risk of distant metastasis. Our study’s limitations include a small number of LUAD patients for subgroup analyses and a single-sample design for investigation of subclonality. Conclusions These data provide a genomic characterization of LUAD pathogenesis and progression. The distinct clonal and subclonal mutation signatures suggest possible diverse carcinogenesis pathways for endogenous and exogenous exposures, and may serve as a foundation for more effective treatments for this lethal disease. LUAD’s high heterogeneity emphasizes the need to further study this tumor type and to associate genomic findings with clinical outcomes.

Mosaic 13q14 deletions in peripheral leukocytes of non-hematologic cancer cases and healthy controls

Loss of 13q14.3 is a chromosomal event found in ~50% of B-cell chronic lymphocytic leukemia (CLL) and monoclonal B-cell lymphocytosis (MBL) cases. Surveys of somatic alterations in solid tumors have shown sporadic 13q14.3 loss in many different tumor types, but not at high frequency in any specific tumor type. In our recent survey of the single-nucleotide polymorphism (SNP) microarray data from 127 000 cancer-free or solid tumor cases, we observed mosaic 13q14.3 loss as common autosomal somatic large structural events (>2 Mb in size) in blood and buccal-derived DNA. Herein, we examined this region more closely investigating structural mosaic events <2 Mb using SNP microarray data in 46 254 non-hematologic cancer cases and 36 229 controls. We detected 60 individuals with 13q14.3 mosaic loss, 1 mosaic copy neutral uniparental disomy and 13 individuals with homozygosity. Although 13q14.3 loss size was variable, the minimally deleted region (MDR) (chr13: 49 590 000–49 983 100; GRCh36) was comparable to what is classically reported in MBL and CLL. Breakpoint analysis of the estimated boundaries reveals enrichment for genes and open chromatin. The frequency of 13q14.3 loss significantly increases with increasing age (P-value=0.028), but was not significantly different between non-hematological cancer cases and controls (0.084% versus 0.058%; P-value=0.19). These findings suggest that mosaic 13q14.3 losses accumulate with age. Individuals with detected mosaic 13q14.3 deletions may be early, undetected cases of MBL or CLL, but not necessarily all will develop MBL and CLL.

Evolution of multiple cell clones over a 29-year period of a CLL patient.

Chronic lymphocytic leukaemia (CLL) is a frequent B-cell malignancy, characterized by recurrent somatic chromosome alterations and a low level of point mutations. Here we present single-nucleotide polymorphism microarray analyses of a single CLL patient over 29 years of observation and treatment, and transcriptome and whole-genome sequencing at selected time points. We identify chromosome alterations 13q14−, 6q− and 12q+ in early cell clones, elimination of clonal populations following therapy, and subsequent appearance of a clone containing trisomy 12 and chromosome 10 copy-neutral loss of heterogeneity that marks a major population dominant at death. Serial single-cell RNA sequencing reveals an expression pattern with high FOS, JUN and KLF4 at disease acceleration, which resolves following therapy, but reoccurs following relapse and death. Transcriptome evolution indicates complex changes in expression occur over time. In conclusion, CLL can evolve gradually during indolent phases, and undergo rapid changes following therapy.