A simple implementation of a normal mixture approach to differential gene expression in multiclass microarrays

Motivation: An important problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. We provide a straightforward and easily implemented method for estimating the posterior probability that an individual gene is null. The problem can be expressed in a two-component mixture framework, using an empirical Bayes approach. Current methods of implementing this approach either have some limitations due to the minimal assumptions made or with more specific assumptions are computationally intensive. Results: By converting to a z-score the value of the test statistic used to test the significance of each gene, we propose a simple two-component normal mixture that models adequately the distribution of this score. The usefulness of our approach is demonstrated on three real datasets.

Differential gene expression in the developing mouse ureter

In many instances, kidney dysgenesis results as a secondary consequence to defects in the development of the ureter. Through the use of mouse genetics a number of genes associated with such malformations have been identified, however, the cause of many other abnormalities remain unknown. In order to identify novel genes involved in ureter development we compared gene expression in embryonic day (E) 12.5, E15.5 and postnatal day (P) 75 ureters using the Compugen mouse long oligo microarrays. A total of 248 genes were dynamically upregulated and 208 downregulated between E12.5 and P75. At E12.5, when the mouse ureter is comprised of a simple cuboidal epithelium surrounded by ureteric mesenchyme, genes previously reported to be expressed in the ureteric mesenchyme, foxC1 and foxC2 were upregulated. By E15.5 the epithelial layer develops into urothelium, impermeable to urine, and smooth muscle develops for the peristaltic movement of urine towards the bladder. The development of these two cell types coincided with the upregulation of UPIIIa, RAB27b and PPAR gamma reported to be expressed in the urothelium, and several muscle genes, Acta1, Tnnt2, Myocd, and Tpm2. In situ hybridization identified several novel genes with spatial expression within the smooth muscle, Acta1; ureteric mesenchyme and smooth muscle, Thbs2 and Co15a2; and urothelium, Kcnj8 and Adh1. This study marks the first known report defining global gene expression of the developing mouse ureter and will provide insight into the molecular mechanisms underlying kidney and lower urinary tract malformations. (c) 2005 Elsevier B.V. All rights reserved.

Detection of differentially expressed genes in synovial fibroblasts by restriction fragment differential display

Objective. To identify differentially expressed genes in synovial fibroblasts and examine the effect on gene expression of exposure to TNF-alpha and IL-1beta. Methods. Restriction fragment differential display was used to isolate genes using degenerate primers complementary to the lysophosphatidic acid acyl transferase gene family. Differential gene expression was confirmed by reverse transcription-polymerase chain reaction and immunohistochemistry using a variety of synovial fibroblasts, including cells from patients with osteoarthritis and self-limiting parvovirus arthritis. Results. Irrespective of disease process, synovial fibroblasts constitutively produced higher levels of IL-6 and monocyte chemoattractant protein 1 (MCP-1) (CCL2) than skin fibroblasts. Seven genes were differentially expressed in synovial fibroblasts compared with skin fibroblasts. Of these genes, four [tissue factor pathway inhibitor 2 (TFPI2), growth regulatory oncogene beta (GRObeta), manganese superoxide dismutase (MnSOD) and granulocyte chemotactic protein 2 (GCP-2)] were all found to be constitutively overexpressed in synoviocytes derived from patients with osteoarthritis. These four genes were only weakly expressed in other synovial fibroblasts (rheumatoid and self-limiting parvovirus infection). However, expression in all types of fibroblasts was increased after stimulation with TNF-alpha and IL-1beta. Three other genes (aggrecan, biglycan and caldesmon) were expressed at higher levels in all types of synovial fibroblasts compared with skin fibroblasts even after stimulation with TNF-alpha and IL-1. Conclusions. Seven genes have been identified with differential expression patterns in terms of disease process (osteoarthritis vs rheumatoid arthritis), state of activation (resting vs cytokine activation) and anatomical location (synovium vs skin). Four of these genes, TFPI2, GRObeta (CXCL2), MnSOD and GCP-2 (CXCL6), were selectively overexpressed in osteoarthritis fibroblasts rather than rheumatoid fibroblasts. While these differences may represent differential behaviour of synovial fibroblasts in in vitro culture, these observations suggest that TFPI2, GRObeta (CXCL2), MnSOD and GCP-2 (CXCL6) may represent new targets for treatments specifically tailored to osteoarthritis.

Expression-based strategies for discovery of genes involved in testis and ovary development

In recent years, strategies for gene identification based on differential gene expression have become increasingly popular, due in part to the development of microarray technology. These strategies are particularly well suited to the identification of genes involved in sex determination and gonadal development, which unlike the development of other organ systems, proceeds along two very different alternative courses, depending on the sex of the embryo. We have used a high-throughput, array-based expression screen to identify several genes expressed sex-specifically in developing mouse gonads. One of these, vanin 1, appears to play a role in mediating migration of mesonephric cells into the male genital ridge. Progress in characterizing other genes arising from the screen is discussed.

A microarray study of post-mortem mRNA degradation in mouse brain tissue

Background: Although there is evidence that post-mortem interval (PMI) is not a major contributor to reduced overall RNA integrity, it may differentially affect a subgroup of gene transcripts that are susceptible to PMI-related degradation. This would particularly have ramifications for microarray studies that include a broad spectrum of genes. Method: Brain tissue was removed from adult mice at 0, 6, 12, 18, 24,36 and 48 h post-mortem. RNA transcript abundance was measured by hybridising RNA from the zero time point with test RNA from each PMI time point, and differential gene expression was assessed using cDNA microarrays. Sequence and ontological analyses were performed on the group of RNA transcripts showing greater than two-fold reduction. Results: Increasing PMI was associated with decreased tissue pH and increased RNA degradation as indexed by 28S/18S ribosomal RNA ratio. Approximately 12% of mRNAs detected on the arrays displayed more than a two-fold decrease in abundance by 48 It post-mortem. An analysis of nucleotide composition provided evidence that transcripts with the AUUUA motif in the 3' untranslated region (3'UTR) were more susceptible to PMI-related RNA degradation, compared to transcripts not carrying the 3'UTR AUUUA motif. Consistent with this finding, ontological analysis showed transcription factors and elements to be over-represented in the group of transcripts susceptible to degradation. Conclusion: A subgroup of mammalian mRNA transcripts are particularly susceptible to PMI-related degradation, and as a group, they are more likely to carry the YUTR AUUUA motif. PMI should be controlled for in human and animal model post-mortem brain studies, particularly those including a broad spectrum of mRNA transcripts. (c) 2005 Elsevier B.V. All rights reserved.

Temporal and spatial transcriptional programs in murine kidney development

We have performed a systematic temporal and spatial expression profiling of the developing mouse kidney using Compugen long-oligonucleotide microarrays. The activity of 18,000 genes was monitored at 24-h intervals from 10.5-day-postcoitum (dpc) metanephric mesenchyme (MM) through to neonatal kidney, and a cohort of 3,600 dynamically expressed genes was identified. Early metanephric development was further surveyed by directly comparing RNA from 10.5 vs. 11.5 vs. 13.5dpc kidneys. These data showed high concordance with the previously published dynamic profile of rat kidney development (Stuart RO, Bush KT, and Nigam SK. Proc Natl Acad Sci USA 98: 5649-5654, 2001) and our own temporal data. Cluster analyses were used to identify gene ontological terms, functional annotations, and pathways associated with temporal expression profiles. Genetic network analysis was also used to identify biological networks that have maximal transcriptional activity during early metanephric development, highlighting the involvement of proliferation and differentiation. Differential gene expression was validated using whole mount and section in situ hybridization of staged embryonic kidneys. Two spatial profiling experiments were also undertaken. MM (10.5dpc) was compared with adjacent intermediate mesenchyme to further define metanephric commitment. To define the genes involved in branching and in the induction of nephrogenesis, expression profiling was performed on ureteric bud (GFP+) FACS sorted from HoxB7-GFP transgenic mice at 15.5dpc vs. the GFP- mesenchymal derivatives. Comparisons between temporal and spatial data enhanced the ability to predict function for genes and networks. This study provides the most comprehensive temporal and spatial survey of kidney development to date, and the compilation of these transcriptional surveys provides important insights into metanephric development that can now be functionally tested.

Identification of three gene candidates for multicellular resistance in colon carcinoma

Solid tumours display elevated resistance to chemo- and radiotherapies compared to individual tumour derived cells. This so-called multicellular resistance (MCR) phenomenon can only be partly explained by reduced diffusion and altered cell cycle status; even fast growing cells on the surface of solid tumours display MCR. Multicellular spheroids (MCS) recapture this phenomenon ex vivo and here we compare gene expression in exponentially growing MCS with gene expression in monolayer culture. Using an 18,664 gene microarray, we identified 42 differentially expressed genes and three of these genes can be linked to potential mechanisms of MCR. A group of interferon response genes were also up-regulated in MCS, as were a number of genes that that are indicative of greater differentiation in three-dimensional cultures.

Differential expression of N-methyl-D-aspartate receptor NR2 isoforms in Alzheimer's disease

We have previously shown that the expression of NMDA receptor NR1 subunit mRNA splice variants in Alzheimer's disease (AD) brain varies according to regional susceptibility to pathological damage. Here we investigated the expression of the modulatory NR2 subunits of the NMDA receptor using quantitative RT-PCR to assay all NR2 isoforms. Significantly lower expression of NR2A and NR2B transcripts was found in susceptible regions of AD brain, whereas expression of NR2C and NR2D transcripts did not differ from that in controls. Western blot analysis confirmed a lower expression of the NR2A and NR2B isoforms at the protein level. The results suggest that NR2 subunit composition may modulate NMDA receptor-mediated excitotoxicity. NMDA receptor dysfunction might give rise to the regionally selective pattern of neuronal loss that is characteristic of AD.

Differential expression of calcium channels in sympathetic and parasympathetic preganglionic inputs to neurons in paracervical ganglia of guinea-pigs

Neurons in pelvic ganglia receive nicotinic excitatory post-synaptic potentials (EPSPs) from sacral preganglionic neurons via the pelvic nerve, lumbar preganglionic neurons via the hypogastric nerve or both. We tested the effect of a range of calcium channel antagonists on EPSPs evoked in paracervical ganglia of female guinea-pigs after pelvic or hypogastric nerve stimulation. omega-Conotoxin GVIA (CTX GVIA, 100 nM) or the novel N-type calcium channel antagonist, CTX CVID (100 nM) reduced the amplitude of EPSPs evoked after pelvic nerve stimulation by 50-75% but had no effect on EPSPs evoked by hypogastric nerve stimulation. Combined addition of CTX GVIA and CTX CVID was no more effective than either antagonist alone. EPSPs evoked by stimulating either nerve trunk were not inhibited by the P/Q calcium channel antagonist, omega-agatoxin IVA (100 nM), nor the L-type calcium channel antagonist, nifedipine (30 muM). SNX 482 (300 nM), an antagonist at some R-type calcium channels, inhibited EPSPs after hypogastric nerve stimulation by 20% but had little effect on EPSPs after pelvic nerve stimulation. Amiloride (100 muM) inhibited EPSPs after stimulation of either trunk by 40%, while nickel (100 muM) was ineffective. CTX GVIA or CTX CVID (100 nM) also slowed the rate of action potential repolarization and reduced afterhyperpolarization amplitude in paracervical neurons. Thus, release of transmitter from the terminals of sacral preganglionic neurons is largely dependent on calcium influx through N-type calcium channels, although an unknown calcium channel which is resistant to selective antagonists also contributes to release. Release of transmitter from lumbar preganglionic neurons does not require calcium entry through either conventional N-type calcium channels or the variant CTX CVID-sensitive N-type calcium channel and seems to be mediated largely by a novel calcium channel. (C) 2004 IBRO. Published by Elsevier Ltd. All rights reserved.

Differential expression of muscle damage in humans following acute fast and slow velocity eccentric exercise

We sought to determine if the velocity of an acute bout of eccentric contractions influenced the duration and severity of several common indirect markers of muscle damage. Subjects performed 36 maximal fast (FST, n=8: 3.14 rad center dot s(-1)) or slow (SLW, n=7: 0.52 rad center dot s(-1)) velocity isokinetic eccentric contractions with the elbow flexors of the non-dominant arm. Muscle soreness, limb girth, plasma creatine kinase (CK) activity, isometric torque and concentric and eccentric torque at 0.52 and 3.14 rad center dot s(-1) were assessed prior to and for several days following the eccentric bout. Peak plasma CK activity was similar in SLW (4030 +/- 1029 U center dot l(-1)) and FST (5864 +/- 2664 U center dot l(-1)) groups, (p > 0.05). Both groups experienced similar decrement in all strength variables during the 48 hr following the eccentric bout. However, recovery occurred more rapidly in the FST group during eccentric (0.52 and 3.14 rad center dot s(-1)) and concentric (3.14 rad center dot s(-1)) post-testing. The severity of muscle soreness was similar in both groups. However, the FST group experienced peak muscle soreness 48 hr later than the SLW group (24 hr vs. 72 hr). The SLW group experienced a greater increase in upper arm girth than the FST group 20 min, 24 hr and 96 hr following the eccentric exercise bout. The contraction velocity of an acute bout of eccentric exercise differentially influences the magnitude and time course of several indirect markers of muscle damage.

Differential expression of the GABA transporters GAT-1 and GAT-3 in brains of rats, cats, monkeys and humans

The homeostasis of GABA is critical to normal brain function. Extracellular levels of GABA are regulated mainly by plasmalemmal gamma-aminobutyric acid (GABA) transporters. Whereas the expression of GABA transporters has been extensively studied in rodents, validation of this data in other species, including humans, has been limited. As this information is crucial for our understanding of therapeutic options in human diseases such as epilepsy, we have compared, by immunocytochemistry, the distributions of the GABA transporters GAT-1 and GAT-3 in rats, cats, monkeys and humans. We demonstrate subtle differences between the results reported in the literature and our results, such as the predominance of GAT-1 labelling in neurons rather than astrocytes in the rat cortex. We note that the optimal localisation of GAT-1 in cats, monkeys and humans requires the use of an antibody against the human sequence carboxyl terminal region of GAT-1 rather than against the slightly different rat sequence. We demonstrate that GAT-3 is localised mainly to astrocytes in hindbrain and midbrain regions of rat brains. However, in species such as cats, monkeys and humans, additional strong immunolabelling of oligodendrocytes has also been observed. We suggest that differences in GAT distribution, especially the expression of GAT-3 by oligodendrocytes in humans, must be accommodated in extrapolating rodent models of GABA homeostasis to humans.

Respiratory gene clusters of Metallosphaera sedula - differential expression and transcriptional organization

Metallosphaera sedula is a thermoacidophilic Crenarchaeon which is capable of leaching metals from sulfidic ores. The authors have investigated the presence and expression of genes encoding respiratory complexes in this organism when grown heterotrophically or chemolithotrophically on either sulfur or pyrite. The presence of three gene clusters, encoding two terminal oxidase complexes, the quinol oxidase SoxABCD and the SoxM oxidase supercomplex, and a gene cluster encoding a high-potential cytochrome b and components of a bc(1) complex analogue (cbsBA-soxL2N gene cluster) was established. Expression studies showed that the soxM gene was expressed to high levels during heterotrophic growth of M. sedula on yeast extract, while the soxABCD mRNA was most abundant in cells grown on sulfur. Reduced-minus-oxidized difference spectra of cell membranes showed cytochrome-related peaks that correspond to published spectra of Sulfolobus-type terminal oxidase complexes. In pyrite-grown cells, expression levels of the two monitored oxidase gene clusters were reduced by a factor of 10-12 relative to maximal expression levels, although spectra of membranes clearly contained oxidase-associated haems, suggesting the presence of additional gene clusters encoding terminal oxidases in M. sedula. Pyrite- and sulfur-grown cells contained high levels of the cbsA transcript, which encodes a membrane-bound cytochrome b with a possible role in iron oxidation or chemolithotrophy. The cbsA gene is not co-transcribed with the soxL2N genes, and therefore does not appear to be an integral part of this bc(1) complex analogue. The data show for the first time the differential expression of the Sulfolobus-type terminal oxidase gene clusters in a Crenarchaeon in response to changing growth modes.

Differential expression of the NMDA NR2B receptor subunit in motoneuron populations susceptible and resistant to amyotrophic lateral sclerosis

We have compared the expression pattern of NMDA receptor subunits (NR1 and NR2A-D)and NRI splice variants (NR1-1a/1b,-2a/2b,-3a/3b,4a/4b) in motor neuron populations from adult Wistar rats that are vulnerable (hypoglossal, XII) or resistant (oculomotor, III) to death in amyotrophic lateral sclerosis (ALS). The major finding was higher levels of expression of the NR2B subunit in the hypoglossal nucleus. Quantitative real-time PCR showed that NR1 was expressed at a greater level than any of the NR2 subunits (> 15 fold greater, P