Li Yan

IMA: An R package for high-throughput analysis of Illumina’s 450K Infinium methylation data

UNLABELLED The Illumina Infinium HumanMethylation450 BeadChip is a newly designed high-density microarray for quantifying the methylation level of over 450 000 CpG sites within human genome. Illumina Methylation Analyzer (IMA) is a computational package designed to automate the pipeline for exploratory analysis and summarization of site-level and region-level methylation changes in epigenetic studies utilizing the 450K DNA methylation microarray. The pipeline loads the data from Illumina platform and provides user-customized functions commonly required to perform exploratory methylation analysis for individual sites as well as annotated regions. AVAILABILITY IMA is implemented in the R language and is freely available from http://www.rforge.net/IMA.

Evaluation of microRNA expression profiles and their associations with risk alleles in lymphoblastoid cell lines of familial ovarian cancer

Interindividual variations of microRNA expression are likely to influence the expression of microRNA target genes and, therefore, contribute to phenotypic differences in humans, including cancer susceptibility. Whether microRNA expression variation has any role in ovarian cancer development is still unknown. Here, we evaluated microRNA expression profiles in lymphoblastoid cell lines from 74 women with familial ovarian cancer and 47 unrelated controls matched on gender and race. We found that the cases and unrelated controls can be clustered using 95 differentially expressed microRNAs with 91% accuracy. To assess the potential implications of microRNAs in ovarian cancer, we investigated the associations between microRNA expression and seven ovarian cancer risk variants discovered from genome-wide association studies (GWAS), namely, rs3814113 on 9p22.2, rs2072590 on 2q31, rs2665390 on 3q25, rs10088218, rs1516982, rs10098821 on 8q24.21 and rs2363956 on 19p13. We observed 130 significant associations at a permutation level of 0.01. Compared with other risk variants, rs3814113 and rs2072590 had the greatest number of significant associations (68 and 37, respectively). Interestingly, 14 microRNAs that were associated with ovarian cancer risk alleles belong to five microRNA clusters. The most notable cluster is the tumorigenic miR-17-92 cluster with five microRNAs, all of which are significantly associated with rs3814113. Using pathway analysis, several key biological pathways were significantly overrepresented, such as cellular response to stress (P = 2.87 × 10(-06)), etc. Further characterization of significant associations between microRNAs and risk alleles could facilitate the understanding of the functions of these GWAS discovered risk alleles in the genetic etiology of ovarian cancer.

Effects of Preanalytic Variables on Circulating MicroRNAs in Whole Blood

Research in the last decade suggests the clinical potential of circulating microRNAs in whole blood as biomarkers for cancer detection. However, before applying the identified circulating microRNAs clinically, biospecimen-focused research has to be performed to identify possible preanalytic variables that may significantly affect the levels of circulating microRNAs. In this study, using a unique resource of the Data Bank and BioRepository (DBBR) at Roswell Park Cancer Institute, we conducted a two-step analysis to identify internal control circulating microRNAs in whole blood and then to study how selected major preanalytic variables (namely, processing delay, storage condition, storage time, and freeze/thaw cycles) might affect the detection of circulating microRNAs. In the discovery phase of the first step, we identified three microRNAs, including miR346, miR134, and miR934, whose levels exhibited the smallest variation between the case–control groups, as well as within each group interindividually. In the further validation analysis, the consistency was validated for miR346 and miR134 but not for miR934. At the second step, using miR346 and miR134 as internal controls, we observed that as the numbers of freeze/thaw cycles increased, levels of both miR346 and miR134 were significantly decreased (Ptrend < 0.0001); varying other processing and storage conditions did not affect miRNA levels. In the paralleled analysis in plasma samples, levels of miR16 were significantly decreased by increasing processing delay and increasing numbers of freeze/thaw cycles but not affected by storage condition and duration. The results from this study highlight the necessity of biospecimen-focused research on circulating microRNAs before clinical utilization. See all the articles in this CEBP Focus section, “Biomarkers, Biospecimens, and New Technologies in Molecular Epidemiology.” Cancer Epidemiol Biomarkers Prev; 23(12); 2643–8. ©2014 AACR.

OSAT: a tool for sample-to-batch allocations in genomics experiments

BackgroundBatch effect is one type of variability that is not of primary interest but ubiquitous in sizable genomic experiments. To minimize the impact of batch effects, an ideal experiment design should ensure the even distribution of biological groups and confounding factors across batches. However, due to the practical complications, the availability of the final collection of samples in genomics study might be unbalanced and incomplete, which, without appropriate attention in sample-to-batch allocation, could lead to drastic batch effects. Therefore, it is necessary to develop effective and handy tool to assign collected samples across batches in an appropriate way in order to minimize the impact of batch effects.ResultsWe describe OSAT (Optimal Sample Assignment Tool), a bioconductor package designed for automated sample-to-batch allocations in genomics experiments.ConclusionsOSAT is developed to facilitate the allocation of collected samples to different batches in genomics study. Through optimizing the even distribution of samples in groups of biological interest into different batches, it can reduce the confounding or correlation between batches and the biological variables of interest. It can also optimize the homogeneous distribution of confounding factors across batches. It can handle challenging instances where incomplete and unbalanced sample collections are involved as well as ideally balanced designs.

NANOG reprograms prostate cancer cells to castration resistance via dynamically repressing and engaging the AR/FOXA1 signaling axis

The pluripotency transcription factor NANOG has been implicated in tumor development, and NANOG-expressing cancer cells manifest stem cell properties that sustain tumor homeostasis, mediate therapy resistance and fuel tumor progression. However, how NANOG converges on somatic circuitry to trigger oncogenic reprogramming remains obscure. We previously reported that inducible NANOG expression propels the emergence of aggressive castration-resistant prostate cancer phenotypes. Here we first show that endogenous NANOG is required for the growth of castration-resistant prostate cancer xenografts. Genome-wide chromatin immunoprecipitation sequencing coupled with biochemical assays unexpectedly reveals that NANOG co-occupies a distinctive proportion of androgen receptor/Forkhead box A1 genomic loci and physically interacts with androgen receptor and Forkhead box A1. Integrative analysis of chromatin immunoprecipitation sequencing and time-resolved RNA sequencing demonstrates that NANOG dynamically alters androgen receptor/Forkhead box A1 signaling leading to both repression of androgen receptor-regulated pro-differentiation genes and induction of genes associated with cell cycle, stem cells, cell motility and castration resistance. Our studies reveal global molecular mechanisms whereby NANOG reprograms prostate cancer cells to a clinically relevant castration-resistant stem cell-like state driven by distinct NANOG-regulated gene clusters that correlate with patient survival. Thus, reprogramming factors such as NANOG may converge on and alter lineage-specific master transcription factors broadly in somatic cancers, thereby facilitating malignant disease progression and providing a novel route for therapeutic resistance.

Overexpression of suppressive microRNAs, miR-30a and miR-200c are associated with improved survival of breast cancer patients

Some microRNAs (miRNAs) are known to suppress breast cancer. However, whether the expressions of these tumor suppressive miRNAs translate to patient survival were not investigated in large cohort. Nine miRNAs (miR-30a, miR-30c, miR-31, miR-126, miR-140, miR-146b, miR-200c, miR-206, and miR-335) known to be tumor suppressive miRNAs in breast cancer were investigated in Genomic Data Common data portal miRNA-Seq dataset and The Cancer Genome Atlas (TCGA) (n = 1052). Of the 9 miRNAs, miR-30a, miR-30c, miR-126, miR-140, miR-206, and miR-335 were found to have significantly lower expression in breast cancer tissues compared to paired normal breast tissue. High expression of miR-30a or miR-200c was associated with significantly better overall survival (OS). Gene Set Enrichment Analysis (GSEA) demonstrated that low expression levels of miR-30a had the tendency to associate with gene enrichment of EMT, while miR-200c did not, in TCGA cohort, and our findings support the need of validation using large cohort to use miRNA as prognostic biomarker for patients with breast cancer.

Replication and validation of genetic polymorphisms associated with survival after allogeneic blood or marrow transplant

Multiple candidate gene-association studies of non-HLA single-nucleotide polymorphisms (SNPs) and outcomes after blood or marrow transplant (BMT) have been conducted. We identified 70 publications reporting 45 SNPs in 36 genes significantly associated with disease-related mortality, progression-free survival, transplant-related mortality, and/or overall survival after BMT. Replication and validation of these SNP associations were performed using DISCOVeRY-BMT (Determining the Influence of Susceptibility COnveying Variants Related to one-Year mortality after BMT), a well-powered genome-wide association study consisting of 2 cohorts, totaling 2888 BMT recipients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome, and their HLA-matched unrelated donors, reported to the Center for International Blood and Marrow Transplant Research. Gene-based tests were used to assess the aggregate effect of SNPs on outcome. None of the previously reported significant SNPs replicated at P < .05 in DISCOVeRY-BMT. Validation analyses showed association with one previously reported donor SNP at P < .05 and survival; more associations would be anticipated by chance alone. No gene-based tests were significant at P < .05. Functional annotation with publicly available data shows these candidate SNPs most likely do not have biochemical function; only 13% of candidate SNPs correlate with gene expression or are predicted to impact transcription factor binding. Of these, half do not impact the candidate gene of interest; the other half correlate with expression of multiple genes. These findings emphasize the peril of pursing candidate approaches and the importance of adequately powered tests of unbiased genome-wide associations with BMT clinical outcomes given the ultimate goal of improving patient outcomes.

Differences in somatic mutation landscape of hepatocellular carcinoma in Asian American and European American populations

The incidence rate of hepatocellular carcinoma (HCC) is higher in populations of Asian ancestry than European ancestry (EA). We sought to investigate HCC mutational differences between the two populations, which may reflect differences in the prevalence of etiological factors. We compared HCC somatic mutations in patients of self-reported Asian American and EA from The Cancer Genome Atlas (TCGA), and assessed associations of tumor mutations with established HCC risk factors. Although the average mutation burden was similar, TP53 and RB1 were mutated at a much higher frequency in Asian Americans than in EAs (TP53: 43% vs. 21%; RB1: 19% vs. 2%). Three putative oncogenic genes, including TRPM3, SAGE1, and ADAMTS7, were mutated exclusively in Asians. In addition, VEGF binding pathway, a druggable target by tyrosine kinase inhibitors such as sorafenib, was mutated at a higher frequency among Asians (13% vs. 2%); while the negative regulation of IL17 production, involved in inflammation and autoimmunity, was mutated only in EAs (12% vs. 0). Accounting for HCC risk factors had little impact on any of the mutational differences. In conclusion, we demonstrated here mutational differences in important cancer genes and pathways between Asian and European ancestries. These differences may have implications for the prevention and treatment of HCC.

Integrity of p53 Associated Pathways Determines Induction of Apoptosis of Tumor Cells Resistant to Aurora-A Kinase Inhibitors

We have previously shown that mammary tumorigenesis in MMTV-Aurora-A mice is further enhanced when p53 is inactivated, demonstrating that integrity of p53 pathway determines phenotypes induced by this oncogenic kinase. In this study, we investigated the roles of genes involved in p53 pathway (p53, Puma, p21, Bax, and Chk2) in response to Aurora-A inhibitors, VX680 and MK-8745, and explored whether chemoresistant tumor cells would further undergo apoptosis with other therapeutic agents. Isogenic HCT116 cell lines were treated with VX680 or MK-8745. Cell cycle analysis, apoptosis, and tumorigenesity were studied. Chemoresistant cells were recovered from xenograft, and further induction of apoptosis was studied. Induction of apoptosis and aneuploidy with VX680 is much stronger than MK-8745. Xenograft assay indicates that tumor growth of HCT116 and HCT116 p53(-) cells are strongly inhibited by VX680, while that of other cell types are similarly inhibited by two compounds. Among the established cell lines recovered from xenografts, MK-8745-resistant clones contain elevated phosphorylation of mTOR and Akt. When further treated with inhibitors of both mTOR and Akt, those cells undergo apoptosis. These results indicate that p53-associated pathway plays a crucial role in regulating growth inhibition of tumor cells when treated with Aurora-A inhibitors. Combined treatment with Akt/mTOR inhibitors can further induce apoptosis of Aurora-A tumors.

Ceramide species are elevated in human breast cancer and are associated with less aggressiveness

Sphingolipids have emerged as key regulatory molecules in cancer cell survival and death. Although important roles of sphingolipids in breast cancer progression have been reported in experimental models, their roles in human patients are yet to be revealed. The aim of this study was to investigate the ceramide levels and its biosynthesis pathways in human breast cancer patients. Breast cancer, peri-tumor and normal breast tissue samples were collected from surgical specimens from a series of 44 patients with breast cancer. The amount of sphingolipid metabolites in the tissue were determined by mass spectrometry. The Cancer Genome Atlas was used to analyze gene expression related to the sphingolipid metabolism. Ceramide levels were higher in breast cancer tissue compared to both normal and peri-tumor breast tissue. Substrates and enzymes that generate ceramide were significantly increased in all three ceramide biosynthesis pathways in cancer. Further, higher levels of ceramide in breast cancer were associated with less aggressive cancer biology presented by Ki-67 index and nuclear grade of the cancer. Interestingly, patients with higher gene expressions of enzymes in the three major ceramide synthesis pathways showed significantly worse prognosis. This is the first study to reveal the clinical relevance of ceramide metabolism in breast cancer patients. We demonstrated that ceramide levels in breast cancer tissue were significantly higher than those in normal tissue, with activation of the three ceramide biosynthesis pathways. We also identified that ceramide levels have a significant association with aggressive phenotype and its enzymes have prognostic impact on breast cancer patients.