Reproducibility is a fundamental requirement in scientific experiments. Some recent publications have claimed that microarrays are unreliable because lists of differentially expressed genes (DEGs) are not reproducible in similar experiments. Meanwhile, new statistical methods for identifying DEGs continue to appear in the scientific literature. The resultant variety of existing and emerging methods exacerbates confusion and continuing debate in the microarray community on the appropriate choice of methods for identifying reliable DEG lists.

Background: Transcriptome sequencing using next-generation sequencing platforms will soon be competing with DNA microarray technologies for global gene expression analysis. As a preliminary evaluation of these promising technologies, we performed deep sequencing of cDNA synthesized from the Microarray Quality Control (MAQC) reference RNA samples using Roche's 454 Genome Sequencer FLX. Results: We generated more that 3.6 million sequence reads of average length 250 bp for the MAQC A and B samples and introduced a data analysis pipeline for translating cDNA read counts into gene expression levels. Using BLAST, 90% of the reads mapped to the human genome and 64% of the reads mapped to the RefSeq database of well annotated genes with e-values ≤ 10^-20. We measured gene expression levels in the A and B samples by counting the numbers of reads that mapped to individual RefSeq genes in multiple sequencing runs to evaluate the MAQC quality metrics for reproducibility, sensitivity, specificity, and accuracy and compared the results with DNA microarrays and Quantitative RT-PCR (QRTPCR) from the MAQC studies. In addition, 88% of the reads were successfully aligned directly to the human genome using the AceView alignment programs with an average 90% sequence similarity to identify 137,899 unique exon junctions, including 22,193 new exon junctions not yet contained in the RefSeq database. Conclusion: Using the MAQC metrics for evaluating the performance of gene expression platforms, the ExpressSeq results for gene expression levels showed excellent reproducibility, sensitivity, and specificity that improved systematically with increasing shotgun sequencing depth, and quantitative accuracy that was comparable to DNA microarrays and QRTPCR. In addition, a careful mapping of the reads to the genome using the AceView alignment programs shed new light on the complexity of the human transcriptome including the discovery of thousands of new splice variants.

Background: Alternative splicing and isoform level expression profiling is an emerging field of interest within genomics. Splicing sensitive microarrays, with probes targeted to individual exons or exon-junctions, are becoming increasingly popular as a tool capable of both expression profiling and finer scale isoform detection. Despite their intuitive appeal, relatively little is known about the performance of such tools, particularly in comparison with more traditional 3' targeted microarrays. Here, we use the well studied Microarray Quality Control (MAQC) dataset to benchmark the Affymetrix Exon Array, and compare it to two other popular platforms: Illumina, and Affymetrix U133. Results: We show that at the gene expression level, the Exon Array performs comparably with the two 3' targeted platforms. However, the interplatform correlation of the results is slightly lower than between the two 3' arrays. We show that some of the discrepancies stem from the RNA amplification protocols, e.g. the Exon Array is able to detect expression of non-polyadenylated transcripts. More importantly, we show that many other differences result from the ability of the Exon Array to monitor more detailed isoform-level changes; several examples illustrate that changes detected by the 3' platforms are actually isoform variations, and that the nature of these variations can be resolved using Exon Array data. Finally, we show how the Exon Array can be used to detect alternative isoform differences, such as alternative splicing, transcript termination, and alternative promoter usage. We discuss the possible pitfalls and false positives resulting from isoform-level analysis. Conclusion: The Exon Array is a valuable tool that can be used to profile gene expression while providing important additional information regarding the types of gene isoforms that are expressed and variable. However, analysis of alternative splicing requires much more hands on effort and visualization of results in order to correctly interpret the data, and generally results in considerably higher false positive rates than expression analysis. One of the main sources of error in the MAQC dataset is variation in amplification efficiency across transcripts, most likely caused by joint effects of elevated GC content in the 5' ends of genes and reduced likelihood of random-primed first strand synthesis in the 3' ends of genes. These effects are currently not adequately corrected using existing statistical methods. We outline approaches to reduce such errors by filtering out potentially problematic data.

Toxin–antitoxin (TA) systems are key regulators of bacterial persistence, a multidrug-tolerant state found in bacterial species that is a major contributing factor to the growing human health crisis of antibiotic resistance. Type II TA systems consist of two proteins, a toxin and an antitoxin; the toxin is neutralized when they form a complex. The ratio of antitoxin to toxin is significantly greater than 1.0 in the susceptible population (non-persister state), but this ratio is expected to become smaller during persistence. Analysis of multiple datasets (RNA-seq, ribosome profiling) and results from translation initiation rate calculators reveal multiple mechanisms that ensure a high antitoxin-to-toxin ratio in the non-persister state. The regulation mechanisms include both translational and transcriptional regulation. We classified E. coli type II TA systems into four distinct classes based on the mechanism of differential protein production between toxin and antitoxin. We find that the most common regulation mechanism is translational regulation. This classification scheme further refines our understanding of one of the fundamental mechanisms underlying bacterial persistence, especially regarding maintenance of the antitoxin-to-toxin ratio.

BRCA1‐associated protein‐1 (BAP1) expression is commonly lost in several tumors including malignant pleural mesothelioma (MPM). Presence or absence of immunohistochemical BAP1 nuclear staining in tumor cells is currently used for differential diagnosis of MPM. In this study, a large cohort of 596 MPM tumors with available clinical data was analyzed to examine associations of BAP1 staining pattern with clinical and molecular features that may reflect the impact of BAP1 mutation on MPM biology. Cases were classified according to the BAP1 staining pattern of tumor cells. Exome and RNA‐sequencing data were available for subsets of cases. Levels of mRNA encoding claudin 15 (CLDN15) and vimentin (VIM) were determined using RT‐qPCR on 483 cases to estimate the relative proportions of epithelial‐like and mesenchymal‐like components in each tumor. Four BAP1 staining patterns were observed: single‐pattern nuclear staining (36%), single‐pattern cytoplasmic staining (25%), single‐pattern absent staining (12%), and combinations of these staining patterns (27%). This study confirmed prior reports that nuclear BAP1 is more frequently associated with wild‐type BAP1 and sarcomatoid histology. However, no associations between BAP1 staining pattern(s) and mutations in specific protein domains and/or mutation type were observed. BAP1 staining patterns were significantly associated (p < 0.001) with BAP1 gene expression, MPM histologic subtypes, molecular clusters, and markers of epithelial‐to‐mesenchymal transition. Frequent observation of combinations of BAP1 staining patterns in MPM tumors indicated intra‐tumoral heterogeneity of BAP1 status. Cytoplasmic BAP1 staining was identified as a putative indicator of favorable prognosis in non‐epithelioid MPM. In conclusion, novel significant associations among different BAP1 staining patterns and subgroups of MPM tumors were observed, suggesting that the role of BAP1 in tumor progression may be more complex than its presumed tumor suppressor function. Cytoplasmic staining was identified as a putative indicator of favorable prognosis in non‐epithelioid MPM, potentially addressing a critical need in clinical decision‐making in this disease. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.