Fold change and p-value cutoffs significantly alter microarray interpretations

<p>Abstract</p> <p>Background</p> <p>As context is important to gene expression, so is the preprocessing of microarray to transcriptomics. Microarray data suffers from several normalization and significance problems. Arbitrary fold change (FC) cut-offs of >2 and significance p-values of <0.02 lead data collection to look only at genes which vary wildly amongst other genes. Therefore, questions arise as to whether the biology or the statistical cutoff are more important within the interpretation. In this paper, we reanalyzed a zebrafish (<it>D. rerio</it>) microarray data set using GeneSpring and different differential gene expression cut-offs and found the data interpretation was drastically different. Furthermore, despite the advances in microarray technology, the array captures a large portion of genes known but yet still leaving large voids in the number of genes assayed, such as leptin a pleiotropic hormone directly related to hypoxia-induced angiogenesis.</p> <p>Results</p> <p>The data strongly suggests that the number of differentially expressed genes is more up-regulated than down-regulated, with many genes indicating conserved signalling to previously known functions. Recapitulated data from Marques et al. (2008) was similar but surprisingly different with some genes showing unexpected signalling which may be a product of tissue (heart) or that the intended response was transient.</p> <p>Conclusions</p> <p>Our analyses suggest that based on the chosen statistical or fold change cut-off; microarray analysis can provide essentially more than one answer, implying data interpretation as more of an art than a science, with follow up gene expression studies a must. Furthermore, gene chip annotation and development needs to maintain pace with not only new genomes being sequenced but also novel genes that are crucial to the overall gene chips interpretation.</p

A Novel Role for Mc1r in the Parallel Evolution of Depigmentation in Independent Populations of the Cavefish Astyanax mexicanus

The evolution of degenerate characteristics remains a poorly understood phenomenon. Only recently has the identification of mutations underlying regressive phenotypes become accessible through the use of genetic analyses. Focusing on the Mexican cave tetra Astyanax mexicanus, we describe, here, an analysis of the brown mutation, which was first described in the literature nearly 40 years ago. This phenotype causes reduced melanin content, decreased melanophore number, and brownish eyes in convergent cave forms of A. mexicanus. Crosses demonstrate non-complementation of the brown phenotype in F2 individuals derived from two independent cave populations: Pachón and the linked Yerbaniz and Japonés caves, indicating the same locus is responsible for reduced pigmentation in these fish. While the brown mutant phenotype arose prior to the fixation of albinism in Pachón cave individuals, it is unclear whether the brown mutation arose before or after the fixation of albinism in the linked Yerbaniz/Japonés caves. Using a QTL approach combined with sequence and functional analyses, we have discovered that two distinct genetic alterations in the coding sequence of the gene Mc1r cause reduced pigmentation associated with the brown mutant phenotype in these caves. Our analysis identifies a novel role for Mc1r in the evolution of degenerative phenotypes in blind Mexican cavefish. Further, the brown phenotype has arisen independently in geographically separate caves, mediated through different mutations of the same gene. This example of parallelism indicates that certain genes are frequent targets of mutation in the repeated evolution of regressive phenotypes in cave-adapted species

An intestinal epithelial defect conferring ER stress results in inflammation involving both innate and adaptive immunity

We recently characterized Winnie mice carrying a missense mutation in Muc2, leading to severe endoplasmic reticulum stress in intestinal goblet cells and spontaneous colitis. In this study, we characterized the immune responses due to this intestinal epithelial dysfunction. In Winnie, there was a fourfold increase in activated dendritic cells (DCs; CD11c+ major histocompatibility complex (MHC) class IIhi) in the colonic lamina propria accompanied by decreased colonic secretion of an inhibitor of DC activation, thymic stromal lymphopoietin (TSLP). Winnie also displayed a significant increase in mRNA expression of the mucosal TH17 signature genes Il17a, IL17f, Tgfb, and Ccr6, particularly in the distal colon. Winnie mesenteric lymph node leukocytes secreted multiple TH1, TH2, and TH17 cytokines on activation, with a large increase in interleukin-17A (IL-17A) progressively with age. A major source of mucosal IL-17A in Winnie was CD4+ T lymphocytes. Loss of T and B lymphocytes in Rag1-/- × Winnie (RaW) crosses did not prevent spontaneous inflammation but did prevent progression with age in the colon but not the cecum. Adoptive transfer of naive T cells into RaW mice caused more rapid and severe colitis than in Rag1-/-, indicating that the epithelial defect results in an intestinal microenvironment conducive to T-cell activation. Thus, the Winnie primary epithelial defect results in complex multicytokine-mediated colitis involving both innate and adaptive immune components with a prominent IL-23/TH17 response, similar to that of human ulcerative colitis

A retinal code for motion along the gravitational and body axes

Self-motion triggers complementary visual and vestibular reflexes supporting image-stabilization and balance. Translation through space produces one global pattern of retinal image motion (optic flow), rotation another. We examined the direction preferences of direction-sensitive ganglion cells (DSGCs) in flattened mouse retinas in vitro. Here we show that for each subtype of DSGC, direction preference varies topographically so as to align with specific translatory optic flow fields, creating a neural ensemble tuned for a specific direction of motion through space. Four cardinal translatory directions are represented, aligned with two axes of high adaptive relevance: the body and gravitational axes. One subtype maximizes its output when the mouse advances, others when it retreats, rises or falls. Two classes of DSGCs, namely, ON-DSGCs and ON-OFF-DSGCs, share the same spatial geometry but weight the four channels differently. Each subtype ensemble is also tuned for rotation. The relative activation of DSGC channels uniquely encodes every translation and rotation. Although retinal and vestibular systems both encode translatory and rotatory self-motion, their coordinate systems differ

A voting approach to identify a small number of highly predictive genes using multiple classifiers

<p>Abstract</p> <p>Background</p> <p>Microarray gene expression profiling has provided extensive datasets that can describe characteristics of cancer patients. An important challenge for this type of data is the discovery of gene sets which can be used as the basis of developing a clinical predictor for cancer. It is desirable that such gene sets be compact, give accurate predictions across many classifiers, be biologically relevant and have good biological process coverage.</p> <p>Results</p> <p>By using a new type of multiple classifier voting approach, we have identified gene sets that can predict breast cancer prognosis accurately, for a range of classification algorithms. Unlike a wrapper approach, our method is not specialised towards a single classification technique. Experimental analysis demonstrates higher prediction accuracies for our sets of genes compared to previous work in the area. Moreover, our sets of genes are generally more compact than those previously proposed. Taking a biological viewpoint, from the literature, most of the genes in our sets are known to be strongly related to cancer.</p> <p>Conclusion</p> <p>We show that it is possible to obtain superior classification accuracy with our approach and obtain a compact gene set that is also biologically relevant and has good coverage of different biological processes.</p

The Leucine Zipper Domains of the Transcription Factors GCN4 and c-Jun Have Ribonuclease Activity

Basic-region leucine zipper (bZIP) proteins are one of the largest transcription factor families that regulate a wide range of cellular functions. Owing to the stability of their coiled coil structure leucine zipper (LZ) domains of bZIP factors are widely employed as dimerization motifs in protein engineering studies. In the course of one such study, the X-ray structure of the retro-version of the LZ moiety of yeast transcriptional activator GCN4 suggested that this retro-LZ may have ribonuclease activity. Here we show that not only the retro-LZ but also the authentic LZ of GCN4 has weak but distinct ribonuclease activity. The observed cleavage of RNA is unspecific, it is not suppressed by the ribonuclease A inhibitor RNasin and involves the breakage of 3′,5′-phosphodiester bonds with formation of 2′,3′-cyclic phosphates as the final products as demonstrated by HPLC/electrospray ionization mass spectrometry. Several mutants of the GCN4 leucine zipper are catalytically inactive, providing important negative controls and unequivocally associating the enzymatic activity with the peptide under study. The leucine zipper moiety of the human factor c-Jun as well as the entire c-Jun protein are also shown to catalyze degradation of RNA. The presented data, which was obtained in the test-tube experiments, adds GCN4 and c-Jun to the pool of proteins with multiple functions (also known as moonlighting proteins). If expressed in vivo, the endoribonuclease activity of these bZIP-containing factors may represent a direct coupling between transcription activation and controlled RNA turnover. As an additional result of this work, the retro-leucine zipper of GCN4 can be added to the list of functional retro-peptides

Identification of a novel splice variant form of the influenza a virus m2 ion channel with an antigenically distinct ectodomain

Segment 7 of influenza A virus produces up to four mRNAs. Unspliced transcripts encode M1, spliced mRNA2 encodes the M2 ion channel, while protein products from spliced mRNAs 3 and 4 have not previously been identified. The M2 protein plays important roles in virus entry and assembly, and is a target for antiviral drugs and vaccination. Surprisingly, M2 is not essential for virus replication in a laboratory setting, although its loss attenuates the virus. To better understand how IAV might replicate without M2, we studied the reversion mechanism of an M2-null virus. Serial passage of a virus lacking the mRNA2 splice donor site identified a single nucleotide pseudoreverting mutation, which restored growth in cell culture and virulence in mice by upregulating mRNA4 synthesis rather than by reinstating mRNA2 production. We show that mRNA4 encodes a novel M2-related protein (designated M42) with an antigenically distinct ectodomain that can functionally replace M2 despite showing clear differences in intracellular localisation, being largely retained in the Golgi compartment. We also show that the expression of two distinct ion channel proteins is not unique to laboratory-adapted viruses but, most notably, was also a feature of the 1983 North American outbreak of H5N2 highly pathogenic avian influenza virus. In identifying a 14th influenza A polypeptide, our data reinforce the unexpectedly high coding capacity of the viral genome and have implications for virus evolution, as well as for understanding the role of M2 in the virus life cycle

Role of educational level in the relationship between Body Mass Index (BMI) and health-related quality of life (HRQL) among rural Spanish women

BACKGROUND: The impact of obesity on health-related quality of life (HRQL) has been little explored in rural areas. The goal of this study is to ascertain the association between obesity and HRQL among Spanish women living in a rural area, and the influence of their educational level. METHODS: Cross-sectional study with personal interview of 1298 women (aged 18 to 60) randomly selected from the electoral rolls of 14 towns in Galicia, a region in the north-west of Spain. HRQL was assessed using the SF-36 questionnaire. The association between body mass index (BMI) and suboptimal scores in the different HRQL dimensions was summarised using odds ratios (ORs), obtained from multivariate logistic regression models. Separate analyses were conducted for women who had finished their education younger than 16 years old and women with secondary education to assess differences in the relationship between BMI and HRQL according to educational level. RESULTS: Among women with primary or lower education, obesity was associated with a higher prevalence of suboptimal values in the following dimensions: Physical functioning (OR: 1.97; 95%CI: 1.22-3.18); Role-physical (OR: 1.81; 95%CI: 1.04-3.14); General health (OR: 1.76; 95%CI: 1.10-2.81); and Role-emotional (OR: 2.52; 95%CI: 1.27-5.03). In women with higher education, physical functioning was the only dimension associated with obesity (OR: 2.02: 95%CI 0.83-4.97). CONCLUSION: The impact of obesity on women's HRQL is greater among those with a lower educational level. This group registered higher prevalence of obesity and poorer self-perceived health.This study was funded by the Instituto de Salud Carlos III (grant 001/05).S

Syndromics: A Bioinformatics Approach for Neurotrauma Research

Substantial scientific progress has been made in the past 50 years in delineating many of the biological mechanisms involved in the primary and secondary injuries following trauma to the spinal cord and brain. These advances have highlighted numerous potential therapeutic approaches that may help restore function after injury. Despite these advances, bench-to-bedside translation has remained elusive. Translational testing of novel therapies requires standardized measures of function for comparison across different laboratories, paradigms, and species. Although numerous functional assessments have been developed in animal models, it remains unclear how to best integrate this information to describe the complete translational “syndrome” produced by neurotrauma. The present paper describes a multivariate statistical framework for integrating diverse neurotrauma data and reviews the few papers to date that have taken an information-intensive approach for basic neurotrauma research. We argue that these papers can be described as the seminal works of a new field that we call “syndromics”, which aim to apply informatics tools to disease models to characterize the full set of mechanistic inter-relationships from multi-scale data. In the future, centralized databases of raw neurotrauma data will enable better syndromic approaches and aid future translational research, leading to more efficient testing regimens and more clinically relevant findings