Importance of replication in microarray gene expression studies: Statistical methods and evidence from repetitive cDNA hybridizations

We present statistical methods for analyzing replicated cDNA microarray expression data and report the results of a controlled experiment. The study was conducted to investigate inherent variability in gene expression data and the extent to which replication in an experiment produces more consistent and reliable findings. We introduce a statistical model to describe the probability that mRNA is contained in the target sample tissue, converted to probe, and ultimately detected on the slide. We also introduce a method to analyze the combined data from all replicates. Of the 288 genes considered in this controlled experiment, 32 would be expected to produce strong hybridization signals because of the known presence of repetitive sequences within them. Results based on individual replicates, however, show that there are 55, 36, and 58 highly expressed genes in replicates 1, 2, and 3, respectively. On the other hand, an analysis by using the combined data from all 3 replicates reveals that only 2 of the 288 genes are incorrectly classified as expressed. Our experiment shows that any single microarray output is subject to substantial variability. By pooling data from replicates, we can provide a more reliable analysis of gene expression data. Therefore, we conclude that designing experiments with replications will greatly reduce misclassification rates. We recommend that at least three replicates be used in designing experiments by using cDNA microarrays, particularly when gene expression data from single specimens are being analyzed.

Repetitive transcranial magnetic stimulation activates specific regions in rat brain

Repetitive transcranial magnetic stimulation (rTMS) is a noninvasive technique to induce electric currents in the brain. Although rTMS is being evaluated as a possible alternative to electroconvulsive therapy for the treatment of refractory depression, little is known about the pattern of activation induced in the brain by rTMS. We have compared immediate early gene expression in rat brain after rTMS and electroconvulsive stimulation, a well-established animal model for electroconvulsive therapy. Our result shows that rTMS applied in conditions effective in animal models of depression induces different patterns of immediate-early gene expression than does electroconvulsive stimulation. In particular, rTMS evokes strong neural responses in the paraventricular nucleus of the thalamus (PVT) and in other regions involved in the regulation of circadian rhythms. The response in PVT is independent of the orientation of the stimulation probe relative to the head. Part of this response is likely because of direct activation, as repetitive magnetic stimulation also activates PVT neurons in brain slices.

Information about movement direction obtained from synchronous activity of motor cortical neurons

Although neuronal synchronization has been shown to exist in primary motor cortex (MI), very little is known about its possible contribution to coding of movement. By using cross-correlation techniques from multi-neuron recordings in MI, we observed that activity of neurons commonly synchronized around the time of movement initiation. For some cell pairs, synchrony varied with direction in a manner not readily predicted by the firing of either neuron. Information theoretic analysis demonstrated quantitatively that synchrony provides information about movement direction beyond that expected by simple rate changes. Thus, MI neurons are not simply independent encoders of movement parameters but rather engage in mutual interactions that could potentially provide an additional coding dimension in cortex.

Systemic spread of an RNA insect virus in plants expressing plant viral movement protein genes

Flock house virus (FHV), a single-stranded RNA insect virus, has previously been reported to cross the kingdom barrier and replicate in barley protoplasts and in inoculated leaves of several plant species [Selling, B. H., Allison, R. F. & Kaesberg, P. (1990) Proc. Natl. Acad. Sci. USA 87, 434–438]. There was no systemic movement of FHV in plants. We tested the ability of movement proteins (MPs) of plant viruses to provide movement functions and cause systemic spread of FHV in plants. We compared the growth of FHV in leaves of nontransgenic and transgenic plants expressing the MP of tobacco mosaic virus or red clover necrotic mosaic virus (RCNMV). Both MPs mobilized cell-to-cell and systemic movement of FHV in Nicotiana benthamiana plants. The yield of FHV was more than 100-fold higher in the inoculated leaves of transgenic plants than in the inoculated leaves of nontransgenic plants. In addition, FHV accumulated in the noninoculated upper leaves of both MP-transgenic plants. RCNMV MP was more efficient in mobilizing FHV to noninoculated upper leaves. We also report here that FHV replicates in inoculated leaves of six additional plant species: alfalfa, Arabidopsis, Brassica, cucumber, maize, and rice. Our results demonstrate that plant viral MPs cause cell-to-cell and long-distance movement of an animal virus in plants and offer approaches to the study of the evolution of viruses and mechanisms governing mRNA trafficking in plants as well as to the development of promising vectors for transient expression of foreign genes in plants.

Inhibition of Auxin Movement from the Shoot into the Root Inhibits Lateral Root Development in Arabidopsis1

In roots two distinct polar movements of auxin have been reported that may control different developmental and growth events. To test the hypothesis that auxin derived from the shoot and transported toward the root controls lateral root development, the two polarities of auxin transport were uncoupled in Arabidopsis. Local application of the auxin-transport inhibitor naphthylphthalamic acid (NPA) at the root-shoot junction decreased the number and density of lateral roots and reduced the free indoleacetic acid (IAA) levels in the root and [3H]IAA transport into the root. Application of NPA to the basal half of or at several positions along the root only reduced lateral root density in regions that were in contact with NPA or in regions apical to the site of application. Lateral root development was restored by application of IAA apical to NPA application. Lateral root development in Arabidopsis roots was also inhibited by excision of the shoot or dark growth and this inhibition was reversible by IAA. Together, these results are consistent with auxin transport from the shoot into the root controlling lateral root development.

Strontium-Induced Repetitive Calcium Spikes in a Unicellular Green Alga1

The divalent cation Sr2+ induced repetitive transient spikes of the cytosolic Ca2+ activity [Ca2+]cy and parallel repetitive transient hyperpolarizations of the plasma membrane in the unicellular green alga Eremosphaera viridis. [Ca2+]cy measurements, membrane potential measurements, and cation analysis of the cells were used to elucidate the mechanism of Sr2+-induced [Ca2+]cy oscillations. Sr2+ was effectively and rapidly compartmentalized within the cell, probably into the vacuole. The [Ca2+]cy oscillations cause membrane potential oscillations, and not the reverse. The endoplasmic reticulum (ER) Ca2+-ATPase blockers 2,5-di-tert-butylhydroquinone and cyclopiazonic acid inhibited Sr2+-induced repetitive [Ca2+]cy spikes, whereas the compartmentalization of Sr2+ was not influenced. A repetitive Ca2+ release and Ca2+ re-uptake by the ER probably generated repetitive [Ca2+]cy spikes in E. viridis in the presence of Sr2+. The inhibitory effect of ruthenium red and ryanodine indicated that the Sr2+-induced Ca2+ release from the ER was mediated by a ryanodine/cyclic ADP-ribose type of Ca2+ channel. The blockage of Sr2+-induced repetitive [Ca2+]cy spikes by La3+ or Gd3+ indicated the necessity of a certain influx of divalent cations for sustained [Ca2+]cy oscillations. Based on these data we present a mathematical model that describes the baseline spiking [Ca2+]cy oscillations in E. viridis.

Instability of repetitive DNA sequences: The role of replication in multiple mechanisms

Rearrangements between tandem sequence homologies of various lengths are a major source of genomic change and can be deleterious to the organism. These rearrangements can result in either deletion or duplication of genetic material flanked by direct sequence repeats. Molecular genetic analysis of repetitive sequence instability in Escherichia coli has provided several clues to the underlying mechanisms of these rearrangements. We present evidence for three mechanisms of RecA-independent sequence rearrangements: simple replication slippage, sister-chromosome exchange-associated slippage, and single-strand annealing. We discuss the constraints of these mechanisms and contrast their properties with RecA-dependent homologous recombination. Replication plays a critical role in the two slipped misalignment mechanisms, and difficulties in replication appear to trigger rearrangements via all these mechanisms.

Antagonistic forces generated by myosin II and cytoplasmic dynein regulate microtubule turnover, movement, and organization in interphase cells

Photoactivation of caged fluorescent tubulin was used mark the microtubule (MT) lattice and monitor MT behavior in interphase cells. A broadening of the photoactivated region occurred as MTs moved bidirectionally. MT movement was not inhibited when MT assembly was suppressed with nocodazole or Taxol; MT movement was suppressed by inhibition of myosin light chain kinase with ML7 or by a peptide inhibitor. Conversely, MT movement was increased after inhibition of cytoplasmic dynein with the antibody 70.1. In addition, the half-time for MT turnover was decreased in cells treated with ML7. These results demonstrate that myosin II and cytoplasmic dynein contribute to a balance of forces that regulates MT organization, movement, and turnover in interphase cells.

Targeted deletion of the mouse POU domain gene Brn-3a causes selective loss of neurons in the brainstem and trigeminal ganglion, uncoordinated limb movement, and impaired suckling.

The Brn-3 subfamily of POU domain genes are expressed in sensory neurons and in select brainstem nuclei. Earlier work has shown that targeted deletion of the Brn-3b and Brn-3c genes produce, respectively, defects in the retina and in the inner ear. We show herein that targeted deletion of the Brn-3a gene results in defective suckling and in uncoordinated limb and trunk movements, leading to early postnatal death. Brn-3a (-/-) mice show a loss of neurons in the trigeminal ganglia, the medial habenula, the red nucleus, and the caudal region of the inferior olivary nucleus but not in the retina and dorsal root ganglia. In the trigeminal and dorsal root ganglia, but not in the retina, there is a marked decrease in the frequency of neurons expressing Brn-3b and Brn-3c, suggesting that Brn-3a positively regulates Brn-3b and Brn-3c expression in somatosensory neurons. Thus, Brn-3a exerts its major developmental effects in somatosensory neurons and in brainstem nuclei involved in motor control. The pheno-types of Brn-3a, Brn-3b, and Brn-3c mutant mice indicate that individual Brn-3 genes have evolved to control development in the auditory, visual, or somatosensory systems and that despite differences between these systems in transduction mechanisms, sensory organ structures, and central information processing, there may be fundamental homologies in the genetic regulatory events that control their development.

Invasion of minor veins of tobacco leaves inoculated with tobacco mosaic virus mutants defective in phloem-dependent movement.

To fully understand vascular transport of plant viruses, the viral and host proteins, their structures and functions, and the specific vascular cells in which these factors function must be determined. We report here on the ability of various cDNA-derived coat protein (CP) mutants of tobacco mosaic virus (TMV) to invade vascular cells in minor veins of Nicotiana tabacum L. cv. Xanthi nn. The mutant viruses we studied, TMV CP-O, U1mCP15-17, and SNC015, respectively, encode a CP from a different tobamovirus (i.e., from odontoglossum ringspot virus) resulting in the formation of non-native capsids, a mutant CP that accumulates in aggregates but does not encapsidate the viral RNA, or no CP. TMV CP-O is impaired in phloem-dependent movement, whereas U1mCP15-17 and SNC015 do not accumulate by phloem-dependent movement. In developmentally-defined studies using immunocytochemical analyses we determined that all of these mutants invaded vascular parenchyma cells within minor veins in inoculated leaves. In addition, we determined that the CPs of TMV CP-O and U1mCP15-17 were present in companion (C) cells of minor veins in inoculated leaves, although more rarely than CP of wild-type virus. These results indicate that the movement of TMV into minor veins does not require the CP, and an encapsidation-competent CP is not required for, but may increase the efficiency of, movement into the conducting complex of the phloem (i.e., the C cell/sieve element complex). Also, a host factor(s) functions at or beyond the C cell/sieve element interface with other cells to allow efficient phloem-dependent accumulation of TMV CP-O.

Ultraviolet-induced movement of the human DNA repair protein, Xeroderma pigmentosum type G, in the nucleus.

Xeroderma pigmentosum type G (XPG) is a human genetic disease exhibiting extreme sensitivity to sunlight. XPG patients are defective XPG endonuclease, which is an enzyme essential for DNA repair of the major kinds of solar ultraviolet (UV)-induced DNA damages. Here we describe a novel dynamics of this protein within the cell nucleus after UV irradiation of human cells. Using confocal microscopy, we have localized the immunofluorescent, antigenic signal of XPG protein to foci throughout the cell nucleus. Our biochemical studies also established that XPG protein forms a tight association with nuclear structure(s). In human skin fibroblast cells, the number of XPG foci decreased within 2 h after UV irradiation, whereas total nuclear XPG fluorescence intensity remained constant, suggesting redistribution of XPG from a limited number of nuclear foci to the nucleus overall. Within 8 h after UV, most XPG antigenic signal was found as foci. Using beta-galactosidase-XPG fusion constructs (beta-gal-XPG) transfected into HeLa cells, we have identified a single region of XPG that is evidently responsible both for foci formation and for the UV dynamic response. The fusion protein carrying the C terminus of XPG (amino acids 1146-1185) localized beta-gal specific antigenic signal to foci and to the nucleolus regions. After UV irradiation, antigenic beta-gal translocated reversibly from the subnuclear structures to the whole nucleus with kinetics very similar to the movements of XPG protein. These findings lead us to propose a model in which distribution of XPG protein may regulate the rate of DNA repair within transcriptionally active and inactive compartments of the cell nucleus.

Assembly and movement of a plant virus carrying a green fluorescent protein overcoat.

Potato virus X (PVX) is a filamentous plant virus infecting many members of the family Solanaceae. A modified form of PVX, PVX.GFP-CP which expressed a chimeric gene encoding a fusion between the 27-kDa Aequorea victoria green fluorescent protein and the amino terminus of the 25-kDa PVX coat protein, assembled into virions and moved both locally and systemically. The PVX.GFP-CP virions were over twice the diameter of wild-type PVX virions. Assembly of PVX.GFP-CP virions required the presence of free coat protein subunits in addition to the fusion protein subunits. PVX.GFP-CP virions accumulated as paracrystalline arrays in infected cells similar to those seen in cells infected with wild-type PVX The formation of virions carrying large superficial fusions illustrates a novel approach for production of high levels of foreign proteins in plants. Aggregates of PVX.GFP-CP particles were fluorescent, emitting green light when excited with ultraviolet light and could be imaged using confocal laser scanning microscopy. The detection of virus particles in infected tissue demonstrates the potential of fusions between the green fluorescent protein and virus coat protein for the non-invasive study of virus multiplication and spread.

Movement of yeast cortical actin cytoskeleton visualized in vivo.

Fusion proteins between the green fluorescent protein (GFP) and the cytoskeleton proteins Act1p (actin), Sac6p (yeast fimbrin homolog), and Abp1p in budding yeast (Saccharomyces cerevisiae) localize to the cortical actin patches. The actin fusions could not function as the sole actin source in yeast, but fusions between the actin-binding proteins Abp1p and Sac6p complement fully the phenotypes associated with their gene deletions. Direct observation in vivo reveals that the actin cortical patches move. Movement of actin patches is constrained to the asymmetric distribution of the patches in growing cells, and this movement is greatly reduced when metabolic inhibitors such as sodium azide are added. Fusion protein-labeled patches are normally distributed during the yeast cell cycle and during mating. In vivo observation made possible the visualization of actin patches during sporulation as well.

A methylated human 9-kb repetitive sequence on acrocentric chromosomes is homologous to a subtelomeric repeat in chimpanzees.

We have implemented an approach for the detection of DNA alterations in cancer by means of computerized analysis of end-labeled genomic fragments, separated in two dimensions. Analysis of two-dimensional patterns of neuroblastoma tumors, prepared by first digesting DNA with the methylation-sensitive restriction enzyme Not I, yielded a multicopy fragment which was detected in some tumor patterns but not in normal controls. Cloning and sequencing of the fragment, isolated from two-dimensional gels, yielded a sequence with a strong homology to a subtelomeric sequence in chimpanzees and which was previously reported to be undetectable in humans. Fluorescence in situ hybridization indicated the occurrence of this sequence in normal tissue, for the most part in the satellite regions of acrocentric chromosomes. A product containing this sequence was obtained by telomere-anchored PCR using as a primer an oligonucleotide sequence from the cloned fragment. Our data suggest demethylation of cytosines at the cloned Not I site and in neighboring DNA in some tumors, compared with normal tissue, and suggest a greater similarity between human and chimpanzee subtelomeric sequences than was previously reported.