Variations in gene sequence and expression underlie much of human variability. Despite the known biological roles of differential allelic gene expression resulting from X-chromosome inactivation and genomic imprinting, a large-scale analysis of allelic gene expression in human is lacking. We examined allele-specific gene expression of 1063 transcribed single-nucleotide polymorphisms (SNPs) by using Affymetrix HuSNP oligo arrays. Among the 602 genes that were heterozygous and expressed in kidney or liver tissues from seven individuals, 326 (54%) showed preferential expression of one allele in at least one individual, and 170 of those showed greater than fourfold difference between the two alleles. The allelic variation has been confirmed by real-time quantitative PCR experiments. Some of these 170 genes are known to be imprinted, such as SNRPN, IPW, HTR2A, and PEG3. Most of the differentially expressed genes are not in known imprinting domains but instead are distributed throughout the genome. Our studies demonstrate that variation of gene expression between alleles is common, and this variation may contribute to human variability.

Purpose: Esophageal squamous cell carcinoma (ESCC) is an aggressive tumor with poor prognosis. Understanding molecular changes in ESCC will enable identification of molecular subtypes and provide potential targets for early detection and therapy. Experimental Design: We followed up a previous array study with additional discovery and confirmatory studies in new ESCC cases by using alternative methods. We profiled global gene expression for discovery and confirmation, and validated selected dysregulated genes with additional RNA and protein studies. Results: A total of 159 genes showed differences with extreme statistical significance (P < E-15) and 2-fold differences or more in magnitude (tumor/normal RNA expression ratio, N = 53 cases), including 116 upregulated and 43 downregulated genes. Of 41 genes dysregulated in our prior array study, all but one showed the same fold change directional pattern in new array studies, including 29 with 2-fold changes or more. Alternative RNA expression methods validated array results: more than two thirds of 51 new cases examined by real-time PCR (RT-PCR) showed 2-fold differences or more for all seven genes assessed. Immunohistochemical protein expression results in 275 cases which were concordant with RNA for five of six genes. Conclusion: We identified an expanded panel of genes dysregulated in ESCC and confirmed previously identified differentially expressed genes. Microarray-based gene expression results were confirmed by RT-PCR and protein expression studies. These dysregulated genes will facilitate molecular categorization of tumor subtypes and identification of their risk factors, and serve as potential targets for early detection, outcome prediction, and therapy. Clin Cancer Res; 17(9); 2955–66. ©2011 AACR.

Background: Genomic instability plays an important role in human cancers. We previously characterized genomic instability in esophageal squamous cell carcinomas (ESCC) in terms of loss of heterozygosity (LOH) and copy number (CN) changes in tumors using the Affymetrix GeneChip Human Mapping 500K array in 30 cases from a high-risk region of China. In the current study we focused on copy number neutral (CN = 2) LOH (CNNLOH) and its relation to gene expression in ESCC.Results: Overall we found that 70% of all LOH observed was CNNLOH. Ninety percent of ESCCs showed CNNLOH (median frequency in cases = 60%) and this was the most common type of LOH in two-thirds of cases. CNNLOH occurred on all 39 autosomal chromosome arms, with highest frequencies on 19p (100%), 5p (96%), 2p (95%), and 20q (95%). In contrast, LOH with CN loss represented 19% of all LOH, occurred in just half of ESCCs (median frequency in cases = 0%), and was most frequent on 3p (56%), 5q (47%), and 21q (41%). LOH with CN gain was 11% of all LOH, occurred in 93% of ESCCs (median frequency in cases = 13%), and was most common on 20p (82%), 8q (74%), and 3q (42%). To examine the effect of genomic instability on gene expression, we evaluated RNA profiles from 17 pairs of matched normal and tumor samples (a subset of the 30 ESCCs) using Affymetrix U133A 2.0 arrays. In CN neutral regions, expression of 168 genes (containing 1976 SNPs) differed significantly in tumors with LOH versus tumors without LOH, including 101 genes that were up-regulated and 67 that were down-regulated.Conclusion: Our results indicate that CNNLOH has a profound impact on gene expression in ESCC, which in turn may affect tumor development.

To gain insight into the role of genomic alterations in breast cancer progression, we conducted a comprehensive genetic characterization of a series of four cell lines derived from MCF10A. MCF10A is an immortalized mammary epithelial cell line (MEC); MCF10AT is a premalignant cell line generated from MCF10A by transformation with an activated HRAS gene; MCF10CA1h and MCF10CA1a, both derived from MCF10AT xenografts, form well-differentiated and poorly-differentiated malignant tumors in the xenograft models, respectively. We analyzed DNA copy number variation using the Affymetrix 500 K SNP arrays with the goal of identifying gene-specific amplification and deletion events. In addition to a previously noted deletion in the CDKN2A locus, our studies identified MYC amplification in all four cell lines. Additionally, we found intragenic deletions in several genes, including LRP1B in MCF10CA1h and MCF10CA1a, FHIT and CDH13 in MCF10CA1h, and RUNX1 in MCF10CA1a. We confirmed the deletion of RUNX1 in MCF10CA1a by DNA and RNA analyses, as well as the absence of the RUNX1 protein in that cell line. Furthermore, we found that RUNX1 expression was reduced in high-grade primary breast tumors compared to low/mid-grade tumors. Mutational analysis identified an activating PIK3CA mutation, H1047R, in MCF10CA1h and MCF10CA1a, which correlates with an increase of AKT1 phosphorylation at Ser473 and Thr308. Furthermore, we showed increased expression levels for genes located in the genomic regions with copy number gain. Thus, our genetic analyses have uncovered sequential molecular events that delineate breast tumor progression. These events include CDKN2A deletion and MYC amplification in immortalization, HRAS activation in transformation, PIK3CA activation in the formation of malignant tumors, and RUNX1 deletion associated with poorly-differentiated malignant tumors.

To profile RNA expression in gastric cancer by anatomic subsites as an initial step in identifying molecular subtypes and providing targets for early detection and therapy. We performed transcriptome analysis using the Affymetrix GeneChip U133A in gastric cardia adenocarcinomas (n = 62) and gastric noncardia adenocarcinomas (n = 72) and their matched normal tissues from patients in Shanxi Province, and validated selected dysregulated genes with additional RNA studies. Expression of dysregulated genes was also related to survival of cases. Principal Component Analysis showed that samples clustered by tumor vs. normal, anatomic location, and histopathologic features. Paired t-tests of tumor/normal tissues identified 511 genes whose expression was dysregulated (P<4.7E-07 and at least two-fold difference in magnitude) in cardia or noncardia gastric cancers, including nearly one-half (n = 239, 47%) dysregulated in both cardia and noncardia, one-fourth dysregulated in cardia only (n = 128, 25%), and about one-fourth in noncardia only (n = 144, 28%). Additional RNA studies confirmed profiling results. Expression was associated with case survival for 20 genes in cardia and 36 genes in noncardia gastric cancers. The dysregulated genes identified here represent a comprehensive starting point for future efforts to understand etiologic heterogeneity, develop diagnostic biomarkers for early detection, and test molecularly-targeted therapies for gastric cancer.