The etiology of depression is still poorly understood, but two major causative hypotheses have been put forth: the monoamine deficiency and the stress hypotheses of depression. We evaluate these ...hypotheses using animal models of endogenous depression and chronic stress. The endogenously depressed rat and its control strain were developed by bidirectional selective breeding from the Wistar-Kyoto (WKY) rat, an accepted model of major depressive disorder (MDD). The WKY More Immobile (WMI) substrain shows high immobility/despair-like behavior in the forced swim test (FST), while the control substrain, WKY Less Immobile (WLI), shows no depressive behavior in the FST. Chronic stress responses were investigated by using Brown Norway, Fischer 344, Lewis and WKY, genetically and behaviorally distinct strains of rats. Animals were either not stressed (NS) or exposed to chronic restraint stress (CRS). Genome-wide microarray analyses identified differentially expressed genes in hippocampi and amygdalae of the endogenous depression and the chronic stress models. No significant difference was observed in the expression of monoaminergic transmission-related genes in either model. Furthermore, very few genes showed overlapping changes in the WMI vs WLI and CRS vs NS comparisons, strongly suggesting divergence between endogenous depressive behavior- and chronic stress-related molecular mechanisms. Taken together, these results posit that although chronic stress may induce depressive behavior, its molecular underpinnings differ from those of endogenous depression in animals and possibly in humans, suggesting the need for different treatments. The identification of novel endogenous depression-related and chronic stress response genes suggests that unexplored molecular mechanisms could be targeted for the development of novel therapeutic agents.
Microarray technology is now widely available and is being applied to address increasingly complex scientific questions. Consequently, there is a greater demand for statistical assessment of the ...conclusions drawn from microarray experiments. This review discusses fundamental issues of how to design an experiment to ensure that the resulting data are amenable to statistical analysis. The discussion focuses on two-color spotted cDNA microarrays, but many of the same issues apply to single-color gene-expression assays as well.
Implementation of wildfire- and climate-adaptation strategies in seasonally dry forests of western North America is impeded by numerous constraints and uncertainties. After more than a century of ...resource and land use change, some question the need for proactive management, particularly given novel social, ecological, and climatic conditions. To address this question, we first provide a framework for assessing changes in landscape conditions and fire regimes. Using this framework, we then evaluate evidence of change in contemporary conditions relative to those maintained by active fire regimes, i.e., those uninterrupted by a century or more of human-induced fire exclusion. The cumulative results of more than a century of research document a persistent and substantial fire deficit and widespread alterations to ecological structures and functions. These changes are not necessarily apparent at all spatial scales or in all dimensions of fire regimes and forest and nonforest conditions. Nonetheless, loss of the once abundant influence of low- and moderate-severity fires suggests that even the least fire-prone ecosystems may be affected by alteration of the surrounding landscape and, consequently, ecosystem functions. Vegetation spatial patterns in fire-excluded forested landscapes no longer reflect the heterogeneity maintained by interacting fires of active fire regimes. Live and dead vegetation (surface and canopy fuels) is generally more abundant and continuous than before European colonization. As a result, current conditions are more vulnerable to the direct and indirect effects of seasonal and episodic increases in drought and fire, especially under a rapidly warming climate. Long-term fire exclusion and contemporaneous social-ecological influences continue to extensively modify seasonally dry forested landscapes. Management that realigns or adapts fire-excluded conditions to seasonal and episodic increases in drought and fire can moderate ecosystem transitions as forests and human communities adapt to changing climatic and disturbance regimes. As adaptation strategies are developed, evaluated, and implemented, objective scientific evaluation of ongoing research and monitoring can aid differentiation of warranted and unwarranted uncertainties.
R/qtl is an extensible, interactive environment for mapping quantitative trait loci (QTLs) in experimental populations derived from inbred lines. It is implemented as an add-on package for the ...freely-available statistical software, R, and includes functions for estimating genetic maps, identifying genotyping errors, and performing single-QTL and two-dimensional, two-QTL genome scans by multiple methods, with the possible inclusion of covariates. Availability: The package is freely available at http://www.biostat.jhsph.edu/~kbroman/qtl. Contact: kbroman@jhsph.edu * To whom correspondence should be addressed. † Present address: Department of Epidemiology and Biostatistics, University of California, San Francisco, CA 94143, USA.
We introduce a general technique for making statistical inference from clustering tools applied to gene expression microarray data. The approach utilizes an analysis of variance model to achieve ...normalization and estimate differential expression of genes across multiple conditions. Statistical inference is based on the application of a randomization technique, bootstrapping. Bootstrapping has previously been used to obtain confidence intervals for estimates of differential expression for individual genes. Here we apply bootstrapping to assess the stability of results from a cluster analysis. We illustrate the technique with a publicly available data set and draw conclusions about the reliability of clustering results in light of variation in the data. The bootstrapping procedure relies on experimental replication. We discuss the implications of replication and good design in microarray experiments.
We describe a general statistical framework for the genetic analysis of quantitative trait data in inbred line crosses. Our main result is based on the observation that, by conditioning on the ...unobserved QTL genotypes, the problem can be split into two statistically independent and manageable parts. The first part involves only the relationship between the QTL and the phenotype. The second part involves only the location of the QTL in the genome. We developed a simple Monte Carlo algorithm to implement Bayesian QTL analysis. This algorithm simulates multiple versions of complete genotype information on a genomewide grid of locations using information in the marker genotype data. Weights are assigned to the simulated genotypes to capture information in the phenotype data. The weighted complete genotypes are used to approximate quantities needed for statistical inference of QTL locations and effect sizes. One advantage of this approach is that only the weights are recomputed as the analyst considers different candidate models. This device allows the analyst to focus on modeling and model comparisons. The proposed framework can accommodate multiple interacting QTL, nonnormal and multivariate phenotypes, covariates, missing genotype data, and genotyping errors in any type of inbred line cross. A software tool implementing this procedure is available. We demonstrate our approach to QTL analysis using data from a mouse backcross population that is segregating multiple interacting QTL associated with salt-induced hypertension.
Combining information across genes in the statistical analysis of microarray data is desirable because of the relatively small number of data points obtained for each individual gene. Here we develop ...an estimator of the error variance that can borrow information across genes using the James-Stein shrinkage concept. A new test statistic (FS) is constructed using this estimator. The new statistic is compared with other statistics used to test for differential expression: the gene-specific F test (F1), the pooled-variance F statistic (F3), a hybrid statistic (F2) that uses the average of the individual and pooled variances, the regularized t-statistic, the posterior odds statistic B, and the SAM t-test. The FS-test shows best or nearly best power for detecting differentially expressed genes over a wide range of simulated data in which the variance components associated with individual genes are either homogeneous or heterogeneous. Thus FS provides a powerful and robust approach to test differential expression of genes that utilizes information not available in individual gene testing approaches and does not suffer from biases of the pooled variance approach.
Evolution may depend more strongly on variation in gene expression than on differences between variant forms of proteins. Regions of DNA that affect gene expression are highly variable, containing ...0.6% polymorphic sites. These naturally occurring polymorphic nucleotides can alter in vivo transcription rates. Thus, one might expect substantial variation in gene expression between individuals. But the natural variation in mRNA expression for a large number of genes has not been measured. Here we report microarray studies addressing the variation in gene expression within and between natural populations of teleost fish of the genus Fundulus. We observed statistically significant differences in expression between individuals within the same population for approximately 18% of 907 genes. Expression typically differed by a factor of 1.5, and often more than 2.0. Differences between populations increased the variation. Much of the variation between populations was a positive function of the variation within populations and thus is most parsimoniously described as random. Some genes showed unexpected patterns of expression-changes unrelated to evolutionary distance. These data suggest that substantial natural variation exists in gene expression and that this quantitative variation is important in evolution.
Spotted cDNA microarrays are emerging as a powerful and cost-effective tool for large-scale analysis of gene expression. Microarrays can be used to measure the relative quantities of specific mRNAs ...in two or more tissue samples for thousands of genes simultaneously. While the power of this technology has been recognized, many open questions remain about appropriate analysis of microarray data. One question is how to make valid estimates of the relative expression for genes that are not biased by ancillary sources of variation. Recognizing that there is inherent "noise" in microarray data, how does one estimate the error variation associated with an estimated change in expression, i.e., how does one construct the error bars? We demonstrate that ANOVA methods can be used to normalize microarray data and provide estimates of changes in gene expression that are corrected for potential confounding effects. This approach establishes a framework for the general analysis and interpretation of microarray data.
The problem of detecting minor quantitative trait loci (QTL) responsible for genetic variation not explained by major QTL is of importance in the complete dissection of quantitative characters. Two ...extensions of the permutation-based method for estimating empirical threshold values are presented. These methods, the conditional empirical threshold (CET) and the residual empirical threshold (RET), yield critical values that can be used to construct tests for the presence of minor QTL effects while accounting for effects of known major QTL. The CET provides a completely nonparametric test through conditioning on markers linked to major QTL. It allows for general nonadditive interactions among QTL but its practical application is restricted to regions of the genome that are unlinked to the major QTL. The RET assumes a structural model for the effect of major QTL, and a threshold is constructed using residuals from this structural model. The search space for minor QTL is unrestricted, and RET-based tests may be more powerful than the CET-based test when the structural model is approximately true