Sample size calculation for microarray experiments with blocked one-way design

<p>Abstract</p> <p>Background</p> <p>One of the main objectives of microarray analysis is to identify differentially expressed genes for different types of cells or treatments. Many statistical methods have been proposed to assess the treatment effects in microarray experiments.</p> <p>Results</p> <p>In this paper, we consider discovery of the genes that are differentially expressed among <it>K </it>(> 2) treatments when each set of <it>K </it>arrays consists of a block. In this case, the array data among <it>K </it>treatments tend to be correlated because of block effect. We propose to use the blocked one-way ANOVA <it>F</it>-statistic to test if each gene is differentially expressed among <it>K </it>treatments. The marginal p-values are calculated using a permutation method accounting for the block effect, adjusting for the multiplicity of the testing procedure by controlling the false discovery rate (FDR). We propose a sample size calculation method for microarray experiments with a blocked one-way design. With FDR level and effect sizes of genes specified, our formula provides a sample size for a given number of true discoveries.</p> <p>Conclusion</p> <p>The calculated sample size is shown via simulations to provide an accurate number of true discoveries while controlling the FDR at the desired level.</p

Alteration of gene expression in mammary gland tissue of dairy cows in response to dietary unsaturated fatty acids

The aim of this study was to determine the effects of supplementing unprotected dietary unsaturated fatty acids (UFAs) from different plant oils on gene expression in the mammary gland of grazing dairy cows. A total of 28 Holstein–Friesian dairy cows in mid-lactation were blocked according to parity, days in milk, milk yield and fat percentage. The cows were then randomly assigned to four UFA sources based on rapeseed, soybean, linseed or a mixture of the three oils for 23 days, after which, all 28 cows were switched to a control diet for an additional 28 days. On the last day of both periods, mammary gland biopsies were taken to study genome-wide differences in gene expression on Affymetrix GeneChip® Bovine Genome Arrays (no. 900493) by ServiceXS (Leiden, The Netherlands). Supplementation with UFAs resulted in increased milk yield but decreased milk fat and protein percentages. Furthermore, the proportion of de novo fatty acids (FAs) in the milk was reduced, whereas that of long-chain FAs increased. Applying a statistical cut-off of false discovery rate of q-value

Environmental stimuli shape microglial plasticity in glioma

In glioma, microglia and infiltrating macrophages are exposed to factors that force them to produce cytokines and chemokines, contributing to tumor growth and maintaining a pro-tumorigenic, immunosuppressed microenvironment. We demonstrate that housing glioma-bearing mice in enriched environment (EE) reverts the immunosuppressive phenotype of infiltrating myeloid cells, by modulating inflammatory gene expression. Under these conditions, branching and patrolling activity of myeloid cells is increased, and their phagocytic activity is promoted. Modulation of gene expression depends on interferon-(IFN) g produced by natural killer (NK) cells, disappearing in mice depleted of NK cells or lacking IFN-g, and was mimicked by exogenous interleukin-15 (IL-15). Further, we describe a key role for BDNF produced in the brain of mice housed in EE in mediating the expression of IL-15 in CD11b+ cells. These data define novel mechanisms linking environmental cues to the acquisition of a pro-inflammatory, anti-tumor microenvironment in mouse brain

Designing Toxicogenomics Studies that use DNA Array Technology

Background: Bioassays are routinely used to evaluate the toxicity of test agents. Experimental designs for bioassays are largely encompassed by fixed effects linear models. In toxicogenomics studies where DNA arrays measure mRNA levels, the tissue samples are typically generated in a bioassay. These measurements introduce additional sources of variation, which must be properly managed to obtain valid tests of treatment effects.Results: An analysis of covariance model is developed which combines a fixed-effects linear model for the bioassay with important variance components associated with DNA array measurements. These models can accommodate the dominant characteristics of measurements from DNA arrays, and they account for technical variation associated with normalization, spots, dyes, and batches as well as the biological variation associated with the bioassay. An example illustrates how the model is used to identify valid designs and to compare competing designs.Conclusions: Many toxicogenomics studies are bioassays which measure gene expression using DNA arrays. These studies can be designed and analyzed using standard methods with a few modifications to account for characteristics of array measurements, such as multiple endpoints and normalization. As much as possible, technical variation associated with probes, dyes, and batches are managed by blocking treatments within these sources of variation. An example shows how some practical constraints can be accommodated by this modelling and how it allows one to objectively compare competing designs

Mutations in CHMP2B in lower motor neuron predominant amyotrophic lateral sclerosis (ALS)

Background: Amyotrophic lateral sclerosis (ALS), a common late-onset neurodegenerative disease, is associated with fronto-temporal dementia (FTD) in 3-10% of patients. A mutation in CHMP2B was recently identified in a Danish pedigree with autosomal dominant FTD. Subsequently, two unrelated patients with familial ALS, one of whom also showed features of FTD, were shown to carry missense mutations in CHMP2B. The initial aim of this study was to determine whether mutations in CHMP2B contribute more broadly to ALS pathogenesis. Methodology/Principal Findings: Sequencing of CHMP2B in 433 ALS cases from the North of England identified 4 cases carrying 3 missense mutations, including one novel mutation, p. Thr104Asn, none of which were present in 500 neurologically normal controls. Analysis of clinical and neuropathological data of these 4 cases showed a phenotype consistent with the lower motor neuron predominant (progressive muscular atrophy (PMA)) variant of ALS. Only one had a recognised family history of ALS and none had clinically apparent dementia. Microarray analysis of motor neurons from CHMP2B cases, compared to controls, showed a distinct gene expression signature with significant differential expression predicting disassembly of cell structure; increased calcium concentration in the ER lumen; decrease in the availability of ATP; down-regulation of the classical and p38 MAPK signalling pathways, reduction in autophagy initiation and a global repression of translation. Transfection of mutant CHMP2B into HEK-293 and COS-7 cells resulted in the formation of large cytoplasmic vacuoles, aberrant lysosomal localisation demonstrated by CD63 staining and impairment of autophagy indicated by increased levels of LC3-II protein. These changes were absent in control cells transfected with wild-type CHMP2B. Conclusions/Significance: We conclude that in a population drawn from North of England pathogenic CHMP2B mutations are found in approximately 1% of cases of ALS and 10% of those with lower motor neuron predominant ALS. We provide a body of evidence indicating the likely pathogenicity of the reported gene alterations. However, absolute confirmation of pathogenicity requires further evidence, including documentation of familial transmission in ALS pedigrees which might be most fruitfully explored in cases with a LMN predominant phenotype

Thermodynamically optimal whole-genome tiling microarray design and validation

Additional file 2. Assembled genome sequences for E. coli MG1655 and Agrobacterium tumefaciens C58. Two FASTA files containing the assembled genomic sequences are provided in this compressed ZIP archive

A microarray configuration to quantify expression levels and relative abundance of splice variants

Over the past decade, alternative RNA splicing has raised a great interest appearing to be of high importance in the generation of expression diversity. This regulatory process plays a critical role in the normal development and its impact on the initiation and development of human disorders as well as on the pharmacological properties of drugs is increasingly being recognized. Only few studies describe specific alternative splicing expression profiling. Microarray strategies have been conceived to address alternative splicing events but with very few experimental data related to their abilities to provide true quantification values. We have developed a specific microarray configuration relying on a few, well optimized probes per splice event. Basically, five probes of 24mer are used to fully characterize a splice event. These probes are of two types, exon probes and junction probes, and are either specific to a splice event or not. The performances of such a ‘splice array’ were validated on synthetic model systems and on complex biological materials. The results indicate that DNA chips based on this design combining exon and junction derived probes enable the detection and, absolute and relative quantification of splice variants. In addition, this strategy is compatible with all the microarrays that use oligonucleotide probes

Single-Cell Phosphoproteomics Resolves Adaptive Signaling Dynamics and Informs Targeted Combination Therapy in Glioblastoma

Intratumoral heterogeneity of signaling networks may contribute to targeted cancer therapy resistance, including in the highly lethal brain cancer glioblastoma (GBM). We performed single-cell phosphoproteomics on a patient-derived in vivo GBM model of mTOR kinase inhibitor resistance and coupled it to an analytical approach for detecting changes in signaling coordination. Alterations in the protein signaling coordination were resolved as early as 2.5 days after treatment, anticipating drug resistance long before it was clinically manifest. Combination therapies were identified that resulted in complete and sustained tumor suppression in vivo. This approach may identify actionable alterations in signal coordination that underlie adaptive resistance, which can be suppressed through combination drug therapy, including non-obvious drug combinations