Tissue microarrays: one size does not fit all

<p>Abstract</p> <p>Background</p> <p>Although tissue microarrays (TMAs) are commonly employed in clinical and basic-science research, there are no guidelines for evaluating the appropriateness of a TMA for a given biomarker and tumor type. Furthermore, TMA performance across multiple biomarkers has not been systematically explored.</p> <p>Methods</p> <p>A simulated TMA with between 1 and 10 cores was designed to study tumor expression of 6 biomarkers with varied expression patterns (B7-H1, B7-H3, survivin, Ki-67, CAIX, and IMP3) using 100 patients with clear cell renal cell carcinoma (RCC). We evaluated agreement between whole tissue section and TMA immunohistochemical biomarker quantification to assess how many TMA cores are necessary to adequately represent RCC whole tissue section expression. Additionally, we evaluated associations of whole tissue section and TMA expression with RCC-specific death.</p> <p>Results</p> <p>The number of simulated TMA cores necessary to adequately represent whole tissue section quantification is biomarker specific. Although 2-3 cores appeared adequate for B7-H3, Ki-67, CAIX, and IMP3, even as many as 10 cores resulted in poor agreement for B7-H1 and survivin compared to RCC whole tissue sections. While whole tissue section B7-H1 was significantly associated with RCC-specific death, no significant associations were detected using as many as 10 TMA cores, suggesting that TMAs can result in false-negative findings if the TMA is not optimally designed.</p> <p>Conclusions</p> <p>Prior to TMA analysis, the number of TMA cores necessary to accurately represent biomarker expression on whole tissue sections should be established as there is not a one-size-fits-all TMA. We illustrate the use of a simulated TMA as a cost-effective tool for this purpose.</p

Protocols for the assurance of microarray data quality and process control

Microarrays represent a powerful technology that provides the ability to simultaneously measure the expression of thousands of genes. However, it is a multi-step process with numerous potential sources of variation that can compromise data analysis and interpretation if left uncontrolled, necessitating the development of quality control protocols to ensure assay consistency and high-quality data. In response to emerging standards, such as the minimum information about a microarray experiment standard, tools are required to ascertain the quality and reproducibility of results within and across studies. To this end, an intralaboratory quality control protocol for two color, spotted microarrays was developed using cDNA microarrays from in vivo and in vitro dose-response and time-course studies. The protocol combines: (i) diagnostic plots monitoring the degree of feature saturation, global feature and background intensities, and feature misalignments with (ii) plots monitoring the intensity distributions within arrays with (iii) a support vector machine (SVM) model. The protocol is applicable to any laboratory with sufficient datasets to establish historical high- and low-quality data

Software comparison for evaluating genomic copy number variation for Affymetrix 6.0 SNP array platform

<p>Abstract</p> <p>Background</p> <p>Copy number data are routinely being extracted from genome-wide association study chips using a variety of software. We empirically evaluated and compared four freely-available software packages designed for Affymetrix SNP chips to estimate copy number: Affymetrix Power Tools (APT), Aroma.Affymetrix, PennCNV and CRLMM. Our evaluation used 1,418 GENOA samples that were genotyped on the Affymetrix Genome-Wide Human SNP Array 6.0. We compared bias and variance in the locus-level copy number data, the concordance amongst regions of copy number gains/deletions and the false-positive rate amongst deleted segments.</p> <p>Results</p> <p>APT had median locus-level copy numbers closest to a value of two, whereas PennCNV and Aroma.Affymetrix had the smallest variability associated with the median copy number. Of those evaluated, only PennCNV provides copy number specific quality-control metrics and identified 136 poor CNV samples. Regions of copy number variation (CNV) were detected using the hidden Markov models provided within PennCNV and CRLMM/VanillaIce. PennCNV detected more CNVs than CRLMM/VanillaIce; the median number of CNVs detected per sample was 39 and 30, respectively. PennCNV detected most of the regions that CRLMM/VanillaIce did as well as additional CNV regions. The median concordance between PennCNV and CRLMM/VanillaIce was 47.9% for duplications and 51.5% for deletions. The estimated false-positive rate associated with deletions was similar for PennCNV and CRLMM/VanillaIce.</p> <p>Conclusions</p> <p>If the objective is to perform statistical tests on the locus-level copy number data, our empirical results suggest that PennCNV or Aroma.Affymetrix is optimal. If the objective is to perform statistical tests on the summarized segmented data then PennCNV would be preferred over CRLMM/VanillaIce. Specifically, PennCNV allows the analyst to estimate locus-level copy number, perform segmentation and evaluate CNV-specific quality-control metrics within a single software package. PennCNV has relatively small bias, small variability and detects more regions while maintaining a similar estimated false-positive rate as CRLMM/VanillaIce. More generally, we advocate that software developers need to provide guidance with respect to evaluating and choosing optimal settings in order to obtain optimal results for an individual dataset. Until such guidance exists, we recommend trying multiple algorithms, evaluating concordance/discordance and subsequently consider the union of regions for downstream association tests.</p

Adult infiltrating gliomas with WHO 2016 integrated diagnosis: additional prognostic roles of ATRX and TERT

The “integrated diagnosis” for infiltrating gliomas in the 2016 revised World Health Organization (WHO) classification of tumors of the central nervous system requires assessment of the tumor for IDH mutations and 1p/19q codeletion. Since TERT promoter mutations and ATRX alterations have been shown to be associated with prognosis, we analyzed whether these tumor markers provide additional prognostic information within each of the five WHO 2016 categories. We used data for 1206 patients from the UCSF Adult Glioma Study, the Mayo Clinic and The Cancer Genome Atlas (TCGA) with infiltrative glioma, grades II-IV for whom tumor status for IDH, 1p/19q codeletion, ATRX, and TERT had been determined. All cases were assigned to one of 5 groups following the WHO 2016 diagnostic criteria based on their morphologic features, and IDH and 1p/19q codeletion status. These groups are: (1) Oligodendroglioma, IDH-mutant and 1p/19q-codeleted; (2) Astrocytoma, IDH-mutant; (3) Glioblastoma, IDH-mutant; (4) Glioblastoma, IDH-wildtype; and (5) Astrocytoma, IDH-wildtype. Within each group, we used univariate and multivariate Cox proportional hazards models to assess associations of overall survival with patient age at diagnosis, grade, and ATRX alteration status and/or TERT promoter mutation status. Among Group 1 IDH-mutant 1p/19q-codeleted oligodendrogliomas, the TERT-WT group had significantly worse overall survival than the TERT-MUT group (HR: 2.72, 95% CI 1.05–7.04, p&nbsp;=&nbsp;0.04). In both Group 2, IDH-mutant astrocytomas and Group 3, IDH-mutant glioblastomas, neither TERT mutations nor ATRX alterations were significantly associated with survival. Among Group 4, IDH-wildtype glioblastomas, ATRX alterations were associated with favorable outcomes (HR: 0.36, 95% CI 0.17–0.81, p&nbsp;=&nbsp;0.01). Among Group 5, IDH-wildtype astrocytomas, the TERT-WT group had significantly better overall survival than the TERT-MUT group (HR: 0.48, 95% CI 0.27–0.87), p&nbsp;=&nbsp;0.02). Thus, we present evidence that in certain WHO 2016 diagnostic groups, testing for TERT promoter mutations or ATRX alterations may provide additional useful prognostic information

Integrated Analysis of Gene Expression, CpG Island Methylation, and Gene Copy Number in Breast Cancer Cells by Deep Sequencing

We used deep sequencing technology to profile the transcriptome, gene copy number, and CpG island methylation status simultaneously in eight commonly used breast cell lines to develop a model for how these genomic features are integrated in estrogen receptor positive (ER+) and negative breast cancer. Total mRNA sequence, gene copy number, and genomic CpG island methylation were carried out using the Illumina Genome Analyzer. Sequences were mapped to the human genome to obtain digitized gene expression data, DNA copy number in reference to the non-tumor cell line (MCF10A), and methylation status of 21,570 CpG islands to identify differentially expressed genes that were correlated with methylation or copy number changes. These were evaluated in a dataset from 129 primary breast tumors. Gene expression in cell lines was dominated by ER-associated genes. ER+ and ER− cell lines formed two distinct, stable clusters, and 1,873 genes were differentially expressed in the two groups. Part of chromosome 8 was deleted in all ER− cells and part of chromosome 17 amplified in all ER+ cells. These loci encoded 30 genes that were overexpressed in ER+ cells; 9 of these genes were overexpressed in ER+ tumors. We identified 149 differentially expressed genes that exhibited differential methylation of one or more CpG islands within 5 kb of the 5′ end of the gene and for which mRNA abundance was inversely correlated with CpG island methylation status. In primary tumors we identified 84 genes that appear to be robust components of the methylation signature that we identified in ER+ cell lines. Our analyses reveal a global pattern of differential CpG island methylation that contributes to the transcriptome landscape of ER+ and ER− breast cancer cells and tumors. The role of gene amplification/deletion appears to more modest, although several potentially significant genes appear to be regulated by copy number aberrations

Influence of obesity-related risk factors in the aetiology of glioma

BACKGROUND: Obesity and related factors have been implicated as possible aetiological factors for the development of glioma in epidemiological observation studies. We used genetic markers in a Mendelian randomisation framework to examine whether obesity-related traits influence glioma risk. This methodology reduces bias from confounding and is not affected by reverse causation. METHODS: Genetic instruments were identified for 10 key obesity-related risk factors, and their association with glioma risk was evaluated using data from a genome-wide association study of 12,488 glioma patients and 18,169 controls. The estimated odds ratio of glioma associated with each of the genetically defined obesity-related traits was used to infer evidence for a causal relationship. RESULTS: No convincing association with glioma risk was seen for genetic instruments for body mass index, waist-to-hip ratio, lipids, type-2 diabetes, hyperglycaemia or insulin resistance. Similarly, we found no evidence to support a relationship between obesity-related traits with subtypes of glioma-glioblastoma (GBM) or non-GBM tumours. CONCLUSIONS: This study provides no evidence to implicate obesity-related factors as causes of glioma

Genome-wide association study identifies multiple risk loci for renal cell carcinoma

Previous genome-wide association studies (GWAS) have identified six risk loci for renal cell carcinoma (RCC). We conducted a meta-analysis of two new scans of 5,198 cases and 7,331 controls together with four existing scans, totalling 10,784 cases and 20,406 controls of European ancestry. Twenty-four loci were tested in an additional 3,182 cases and 6,301 controls. We confirm the six known RCC risk loci and identify seven new loci at 1p32.3 (rs4381241, P=3.1 × 10−10), 3p22.1 (rs67311347, P=2.5 × 10−8), 3q26.2 (rs10936602, P=8.8 × 10−9), 8p21.3 (rs2241261, P=5.8 × 10−9), 10q24.33-q25.1 (rs11813268, P=3.9 × 10−8), 11q22.3 (rs74911261, P=2.1 × 10−10) and 14q24.2 (rs4903064, P=2.2 × 10−24). Expression quantitative trait analyses suggest plausible candidate genes at these regions that may contribute to RCC susceptibility

Genomic copy number variation in Mus musculus.

BACKGROUND: Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously. RESULTS: We found 9634 putative autosomal CNVs across the samples affecting 6.87% of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR). CONCLUSION: The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies