The Reproducibility of Lists of Differentially Expressed Genes in Microarray Studies

Reproducibility is a fundamental requirement in scientific experiments and clinical contexts. Recent publications raise concerns about the reliability of microarray technology because of the apparent lack of agreement between lists of differentially expressed genes (DEGs). In this study we demonstrate that (1) such discordance may stem from ranking and selecting DEGs solely by statistical significance (P) derived from widely used simple t-tests; (2) when fold change (FC) is used as the ranking criterion, the lists become much more reproducible, especially when fewer genes are selected; and (3) the instability of short DEG lists based on P cutoffs is an expected mathematical consequence of the high variability of the t-values. We recommend the use of FC ranking plus a non-stringent P cutoff as a baseline practice in order to generate more reproducible DEG lists. The FC criterion enhances reproducibility while the P criterion balances sensitivity and specificity

The ‘heritagisation’ of the British seaside resort: The rise of the ‘old penny’ arcade.

Amusement arcades have long been a key component of the British seaside resort. For almost a century, they enjoyed popularity and success and became established as a quintessential feature of the British seaside holiday. However, the advent of home-based video games along with recent gambling legislation has led to a decline of the seaside amusement arcade sector. Arcades gained a reputation as unsavoury places and their appearance and fortunes often mirrored those of the resorts in which they were located. However, over the past decade, a new variant of the seaside amusement arcade has appeared, featuring mechanical machines working on pre-decimal currency. Such ‘old penny arcades’ frequently describe themselves as museums or heritage centres and they offer an experience based on a nostalgic affection for the ‘traditional’ seaside holiday. They have appeared in the context of an increasing interest in the heritage of the British seaside resort and constitute one element of the ‘heritagisation’ of such resorts. This paper argues that such arcades can be important elements of strategies to reposition and rebrand resorts for the heritage tourism market

Genetic analysis of the molecular regulation of electric fields-guided glia migration

© 2020, The Author(s). Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.In a developing nervous system, endogenous electric field (EF) influence embryonic growth. We reported the EF-directed migration of both rat Schwann cells (SCs) and oligodendrocyte precursor cells (OPCs) and explored the molecular mechanism using RNA-sequencing assay. However, previous studies revealed the differentially expressed genes (DEGs) associated with EF-guided migration of SCs or OPCs alone. In this study, we performed joint differential expression analysis on the RNA-sequencing data from both cell types. We report a number of significantly enriched gene ontology (GO) terms that are related to the cytoskeleton, cell adhesion, and cell migration. Of the DEGs associated with these terms, nine up-regulated DEGs and 32 down-regulated DEGs showed the same direction of effect in both SCs and OPCs stimulated with EFs, while the remaining DEGs responded differently. Thus, our study reveals the similarities and differences in gene expression and cell migration regulation of different glial cell types in response to EF stimulation.National Institute of General Medical Sciences (P20 GM103418) from the National Institutes of Health

Analysis of the Tribolium homeotic complex: insights into mechanisms constraining insect Hox clusters

The remarkable conservation of Hox clusters is an accepted but little understood principle of biology. Some organizational constraints have been identified for vertebrate Hox clusters, but most of these are thought to be recent innovations that may not apply to other organisms. Ironically, many model organisms have disrupted Hox clusters and may not be well-suited for studies of structural constraints. In contrast, the red flour beetle, Tribolium castaneum, which has a long history in Hox gene research, is thought to have a more ancestral-type Hox cluster organization. Here, we demonstrate that the Tribolium homeotic complex (HOMC) is indeed intact, with the individual Hox genes in the expected colinear arrangement and transcribed from the same strand. There is no evidence that the cluster has been invaded by non-Hox protein-coding genes, although expressed sequence tag and genome tiling data suggest that noncoding transcripts are prevalent. Finally, our analysis of several mutations affecting the Tribolium HOMC suggests that intermingling of enhancer elements with neighboring transcription units may constrain the structure of at least one region of the Tribolium cluster. This work lays a foundation for future studies of the Tribolium HOMC that may provide insights into the reasons for Hox cluster conservation

Large-scale insertional mutagenesis of a coleopteran stored grain pest, the red flour beetle Tribolium castaneum, identifies embryonic lethal mutations and enhancer traps

<p>Abstract</p> <p>Background</p> <p>Given its sequenced genome and efficient systemic RNA interference response, the red flour beetle <it>Tribolium castaneum </it>is a model organism well suited for reverse genetics. Even so, there is a pressing need for forward genetic analysis to escape the bias inherent in candidate gene approaches.</p> <p>Results</p> <p>To produce easy-to-maintain insertional mutations and to obtain fluorescent marker lines to aid phenotypic analysis, we undertook a large-scale transposon mutagenesis screen. In this screen, we produced more than 6,500 new <it>piggyBac </it>insertions. Of these, 421 proved to be recessive lethal, 75 were semi-lethal, and eight indicated recessive sterility, while 505 showed new enhancer-trap patterns. Insertion junctions were determined for 403 lines and often appeared to be located within transcription units. Insertion sites appeared to be randomly distributed throughout the genome, with the exception of a preference for reinsertion near the donor site.</p> <p>Conclusion</p> <p>A large collection of enhancer-trap and embryonic lethal beetle lines has been made available to the research community and will foster investigations into diverse fields of insect biology, pest control, and evolution. Because the genetic elements used in this screen are species-nonspecific, and because the crossing scheme does not depend on balancer chromosomes, the methods presented herein should be broadly applicable for many insect species.</p

Anti-inflammatory activity and neutrophil reductions mediated by the JAK1/JAK3 inhibitor, CP-690,550, in rat adjuvant-induced arthritis

<p>Abstract</p> <p>Background</p> <p>The Janus kinase (JAK) family of tyrosine kinases includes JAK1, JAK2, JAK3 and TYK2, and is required for signaling through Type I and Type II cytokine receptors. CP-690,550 is a potent and selective JAK inhibitor currently in clinical trials for rheumatoid arthritis (RA) and other autoimmune disease indications. In RA trials, dose-dependent decreases in neutrophil counts (PBNC) were observed with CP-690,550 treatment. These studies were undertaken to better understand the relationship between JAK selectivity and PBNC decreases observed with CP-690,550 treatment.</p> <p>Methods</p> <p>Potency and selectivity of CP-690,550 for mouse, rat and human JAKs was evaluated in a panel of <it>in vitro </it>assays. The effect of CP-690,550 on granulopoiesis from progenitor cells was also assessed <it>in vitro </it>using colony forming assays. <it>In vivo </it>the potency of orally administered CP-690,550 on arthritis (paw edema), plasma cytokines, PBNC and bone marrow differentials were evaluated in the rat adjuvant-induced arthritis (AIA) model.</p> <p>Results</p> <p>CP-690,550 potently inhibited signaling through JAK1 and JAK3 with 5-100 fold selectivity over JAK2 in cellular assays, despite inhibiting all four JAK isoforms with nM potency in <it>in vitro </it>enzyme assays. Dose-dependent inhibition of paw edema was observed <it>in vivo </it>with CP-690,550 treatment. Plasma cytokines (IL-6 and IL-17), PBNC, and bone marrow myeloid progenitor cells were elevated in the context of AIA disease. At efficacious exposures, CP-690,550 returned all of these parameters to pre-disease levels. The plasma concentration of CP-690,550 at efficacious doses was above the <it>in vitro </it>whole blood IC50 of JAK1 and JAK3 inhibition, but not that of JAK2.</p> <p>Conclusion</p> <p>Results from this investigation suggest that CP-690,550 is a potent inhibitor of JAK1 and JAK3 with potentially reduced cellular potency for JAK2. In rat AIA, as in the case of human RA, PBNC were decreased at efficacious exposures of CP-690,550. Inflammatory end points were similarly reduced, as judged by attenuation of paw edema and cytokines IL-6 and IL-17. Plasma concentration at these exposures was consistent with inhibition of JAK1 and JAK3 but not JAK2. Decreases in PBNC following CP-690,550 treatment may thus be related to attenuation of inflammation and are likely not due to suppression of granulopoiesis through JAK2 inhibition.</p

Effects of Maternal Swine Appeasing Substance on Growth Performance of Nursery Pigs

Two experiments using a total of 640 weanling pigs were conducted to determine the effects of a maternal swine appeasing substance (FerAppease, FERA Diagnostics and Biologicals, College Station, TX) on growth performance, salivary cortisol concentration, and fecal microbiome of nursery pigs. In exp. 1, a total of 360 weanling pigs (DNA 600 × 241; initially 12.6 lb) were used in a 24-d study with five pigs per pen. Due to the potential contact of pigs in adjacent pens transferring the test product, pens were grouped by treatment within the barn and groups were distributed in the barn to minimize this risk. In exp. 1, there were two blocks of pens for each treatment (four total blocks of pens), with each block having 18 pens. In exp. 2, a total of 280 weanling pigs (DNA 600 × 241; initially 12.1 lb) were used in a 45-d study with five pigs per pen in 56 pens split in two rooms. Treatment groups were located in blocks of seven pens to avoid contact between treatment groups, and each treatment was assigned to four blocks (eight blocks total). Block of like treatments was considered the experimental unit for analysis of growth performance data. In both experiments, pens were assigned either a control (placebo of mineral oil) or the appeasing substance (FerAppease, FERA diagnostics and Biologicals, College Station TX). Applications of 3 mL/pig were done at d 0 (weaning) and d 10 by spraying the substance on the forehead of each individual pig. To limit any cross-contamination, either empty pens or the walkway separated the groups to physically segregate treatments from each other. All pigs were fed the same two-phase diet regimens: from d 0 to 10 and d 11 to 24, with phase 3 added in exp. 2 from d 25 to 45. In phase 1 (d 0 to 10, Table 2), control and FerAppease pigs had similar (P \u3e 0.05) ADG, ADFI, and F/G, as well as BW at d 10. Similar results were observed in phase 2 with no evidence of a difference observed between treatments (P \u3e 0.05). For the overall period, when combining phases 1 and 2 (both experiments) and 1, 2 and 3 (only in exp. 2), no differences (P \u3e 0.05) between the two treatments were found for ADG, ADFI, F/G, and BW at d 24 and d 45. Salivary cortisol was measured on d 1 post-weaning and again on d 11 (both 24 h after product administration). There was no evidence (P \u3e 0.10) of a treatment × day interaction or main effect of treatment. Salivary cortisol was greater (P \u3c 0.0001) on d 1 post-weaning compared to d 11 post-weaning. Numerically, on d 1 post-weaning pigs applied with the FerAppease had reduced salivary cortisol concentrations; however, there were no statistically significant differences observed. Fecal samples were collected from two pigs per pen on d 0 and 10, and fecal microbiome analysis did not result in any differences between treatments for alpha or beta diversity, and no meaningful differences were observed in evaluation of class and genus relative abundance data. In summary, the post-weaning application of maternal pheromone to the forehead of nursery pigs post-weaning did not result in any improvements in feed intake, growth rate, feed efficiency, or fecal microbial composition

Enhanced identification and biological validation of differential gene expression via Illumina whole-genome expression arrays through the use of the model-based background correction methodology

Despite the tremendous growth of microarray usage in scientific studies, there is a lack of standards for background correction methodologies, especially in single-color microarray platforms. Traditional background subtraction methods often generate negative signals and thus cause large amounts of data loss. Hence, some researchers prefer to avoid background corrections, which typically result in the underestimation of differential expression. Here, by utilizing nonspecific negative control features integrated into Illumina whole genome expression arrays, we have developed a method of model-based background correction for BeadArrays (MBCB). We compared the MBCB with a method adapted from the Affymetrix robust multi-array analysis algorithm and with no background subtraction, using a mouse acute myeloid leukemia (AML) dataset. We demonstrated that differential expression ratios obtained by using the MBCB had the best correlation with quantitative RT–PCR. MBCB also achieved better sensitivity in detecting differentially expressed genes with biological significance. For example, we demonstrated that the differential regulation of Tnfr2, Ikk and NF-kappaB, the death receptor pathway, in the AML samples, could only be detected by using data after MBCB implementation. We conclude that MBCB is a robust background correction method that will lead to more precise determination of gene expression and better biological interpretation of Illumina BeadArray data

Impact of the spotted microarray preprocessing method on fold-change compression and variance stability

<p>Abstract</p> <p>Background</p> <p>The standard approach for preprocessing spotted microarray data is to subtract the local background intensity from the spot foreground intensity, to perform a log2 transformation and to normalize the data with a global median or a lowess normalization. Although well motivated, standard approaches for background correction and for transformation have been widely criticized because they produce high variance at low intensities. Whereas various alternatives to the standard background correction methods and to log2 transformation were proposed, impacts of both successive preprocessing steps were not compared in an objective way.</p> <p>Results</p> <p>In this study, we assessed the impact of eight preprocessing methods combining four background correction methods and two transformations (the log2 and the glog), by using data from the MAQC study. The current results indicate that most preprocessing methods produce fold-change compression at low intensities. Fold-change compression was minimized using the Standard and the Edwards background correction methods coupled with a log2 transformation. The drawback of both methods is a high variance at low intensities which consequently produced poor estimations of the p-values. On the other hand, effective stabilization of the variance as well as better estimations of the p-values were observed after the glog transformation.</p> <p>Conclusion</p> <p>As both fold-change magnitudes and p-values are important in the context of microarray class comparison studies, we therefore recommend to combine the Edwards correction with a hybrid transformation method that uses the log2 transformation to estimate fold-change magnitudes and the glog transformation to estimate p-values.</p

Evaluating methods for ranking differentially expressed genes applied to microArray quality control data

<p>Abstract</p> <p>Background</p> <p>Statistical methods for ranking differentially expressed genes (DEGs) from gene expression data should be evaluated with regard to high sensitivity, specificity, and reproducibility. In our previous studies, we evaluated eight gene ranking methods applied to only Affymetrix GeneChip data. A more general evaluation that also includes other microarray platforms, such as the Agilent or Illumina systems, is desirable for determining which methods are suitable for each platform and which method has better inter-platform reproducibility.</p> <p>Results</p> <p>We compared the eight gene ranking methods using the MicroArray Quality Control (MAQC) datasets produced by five manufacturers: Affymetrix, Applied Biosystems, Agilent, GE Healthcare, and Illumina. The area under the curve (AUC) was used as a measure for both sensitivity and specificity. Although the highest AUC values can vary with the definition of "true" DEGs, the best methods were, in most cases, either the weighted average difference (WAD), rank products (RP), or intensity-based moderated <it>t </it>statistic (ibmT). The percentages of overlapping genes (POGs) across different test sites were mainly evaluated as a measure for both intra- and inter-platform reproducibility. The POG values for WAD were the highest overall, irrespective of the choice of microarray platform. The high intra- and inter-platform reproducibility of WAD was also observed at a higher biological function level.</p> <p>Conclusion</p> <p>These results for the five microarray platforms were consistent with our previous ones based on 36 real experimental datasets measured using the Affymetrix platform. Thus, recommendations made using the MAQC benchmark data might be universally applicable.</p