Generalization of DNA microarray dispersion properties: microarray equivalent of t-distribution

Background

DNA microarrays are a powerful technology that can provide a wealth of gene expression data for disease studies, drug development, and a wide scope of other investigations. Because of the large volume and inherent variability of DNA microarray data, many new statistical methods have been developed for evaluating the significance of the observed differences in gene expression. However, until now little attention has been given to the characterization of dispersion of DNA microarray data.

Results

Here we examine the expression data obtained from 682 Affymetrix GeneChips^® with 22 different types and we demonstrate that the Gaussian (normal) frequency distribution is characteristic for the variability of gene expression values. However, typically 5 to 15% of the samples deviate from normality. Furthermore, it is shown that the frequency distributions of the difference of expression in subsets of ordered, consecutive pairs of genes (consecutive samples) in pair-wise comparisons of replicate experiments are also normal. We describe a consecutive sampling method, which is employed to calculate the characteristic function approximating standard deviation and show that the standard deviation derived from the consecutive samples is equivalent to the standard deviation obtained from individual genes. Finally, we determine the boundaries of probability intervals and demonstrate that the coefficients defining the intervals are independent of sample characteristics, variability of data, laboratory conditions and type of chips. These coefficients are very closely correlated with Student's t-distribution.

Conclusion

In this study we ascertained that the non-systematic variations possess Gaussian distribution, determined the probability intervals and demonstrated that the K _α coefficients defining these intervals are invariant; these coefficients offer a convenient universal measure of dispersion of data. The fact that the K _α distributions are so close to t-distribution and independent of conditions and type of arrays suggests that the quantitative data provided by Affymetrix technology give "true" representation of physical processes, involved in measurement of RNA abundance.

Reviewers

This article was reviewed by Yoav Gilad (nominated by Doron Lancet), Sach Mukherjee (nominated by Sandrine Dudoit) and Amir Niknejad and Shmuel Friedland (nominated by Neil Smalheiser).     

Chk1 and p21 cooperate to prevent apoptosis during DNA replication fork stress

Cells respond to DNA replication stress by triggering cell cycle checkpoints, repair, or death. To understand the role of the DNA damage response pathways in determining whether cells survive replication stress or become committed to death, we examined the effect of loss of these pathways on cellular response to agents that slow or arrest DNA synthesis. We show that replication inhibitors such as excess thymidine, hydroxyurea, and camptothecin are normally poor inducers of apoptosis. However, these agents become potent inducers of death in S-phase cells upon small interfering RNA-mediated depletion of the checkpoint kinase Chk1. This death response is independent of p53 and Chk2. p21-deficient cells, on the other hand, produce a more robust apoptotic response upon Chk1 depletion. p21 is normally induced only late after thymidine treatment. In Chk1-depleted cells p21 induction occurs earlier and does not require p53. Thus, Chk1 plays a primary role in the protection of cells from death induced by replication fork stress, whereas p21 mediates through its role in regulating entry into S phase. These findings are of potential importance to cancer therapy because we demonstrate that the efficacy of clinically relevant agents can be enhanced by manipulation of these signaling pathways.     

Qatari DNA variation at a crossroad of human migrations

Genomic diversity of the Qatari population was investigated by screening 15 autosomal short tandem repeats (STRs). Significant departures from genetic equilibrium were detected at the D13S317, D19S433 and VWA loci, which persisted after applying Bonferroni-type corrections. Gene diversity (GD) values ranged from 0.6851 (TPOX) to 0.8813 (D2S1338), while observed heterozygosity (Ho) oscillated between 0.3388 (D19S433) and 0.8397 (D2S1338). Interestingly, Ho was lower than expected (He) for 14 of the loci analyzed. The information provided by these microsatellite markers was analyzed by means of genetic distances, multidimensional scaling, hierarchical analyses of the molecular variance (AMOVA) and admixture estimations to assess the genetic relationships of Qatar with European, Asian, African and other Middle Eastern populations. The main findings of the study were the genetic uniqueness of the Qatari population, its strong similarity to the United Arab Emirates (UAE) group, a slight genetic differentiation with respect to other Arab populations (Syria and Egypt) and Turkey, and a certain genetic affinity with sub-Saharan African populations. These results are discussed in light of two major issues: the high consanguinity rates characterizing the Qatari population and its strategic geographic position in the Arabian Peninsula close to major migratory routes, an important pivotal contact zone for bidirectional dispersals between Eurasia and Africa.     

Important role for the transmembrane domain of severe acute respiratory syndrome coronavirus spike protein during entry

The spike protein (S) of severe acute respiratory syndrome coronavirus (SARS-CoV) is responsible for receptor binding and membrane fusion. It contains a highly conserved transmembrane domain that consists of three parts: an N-terminal tryptophan-rich domain, a central domain, and a cysteine-rich C-terminal domain. The cytoplasmic tail of S has previously been shown to be required for assembly. Here, the roles of the transmembrane and cytoplasmic domains of S in the infectivity and membrane fusion activity of SARS-CoV have been studied. SARS-CoV S-pseudotyped retrovirus (SARSpp) was used to measure S-mediated infectivity. In addition, the cell-cell fusion activity of S was monitored by a Renilla luciferase-based cell-cell fusion assay. S(VSV-Cyt), an S chimera with a cytoplasmic tail derived from vesicular stomatitis virus G protein (VSV-G), and S(MHV-TMDCyt), an S chimera with the cytoplasmic and transmembrane domains of mouse hepatitis virus, displayed wild-type-like activity in both assays. S(VSV-TMDCyt), a chimera with the cytoplasmic and transmembrane domains of VSV-G, was impaired in the SARSpp and cell-cell fusion assays, showing 3 to 25% activity compared to the wild type, depending on the assay and the cells used. Examination of the oligomeric state of the chimeric S proteins in SARSpp revealed that S(VSV-TMDCyt) trimers were less stable than wild-type S trimers, possibly explaining the lowered fusogenicity and infectivity.     

Mass spectrometric investigation of the DNA-binding properties of an anthracycline with two trisaccharide chains

Cosmomycin D (CosD) is an anthracycline that has two trisaccharide chains linked to its ring system. Gel electrophoresis showed that CosD formed stable complexes with plasmid DNA under conditions where daunorubicin (Dn) and doxorubicin (Dx) dissociated to some extent during the experiments. The footprint and stability of CosD complexed with 10- and 16 mer DNA was investigated using several applications of electrospray ionisation mass spectrometry (ESI-MS). ESI-MS binding profiles showed that fewer CosD molecules bound to the sequences than Dn or Dx. In agreement with this, ESI-MS analysis of nuclease digestion products of the complexes showed that CosD protected the DNA to a greater extent than Dn or Dx. In tandem MS experiments, all CosD–DNA complexes were more stable than Dn- and Dx-DNA complexes. These results support that CosD binds more tightly to DNA and exerts a larger footprint than Dn or Dx. ESI-MS investigations of the binding properties of CosD could be carried out rapidly and using only small amounts of sample.     

Co-regulation of Gremlin and Notch signalling in diabetic nephropathy

Diabetic nephropathy is currently the leading cause of end-stage renal disease worldwide, and occurs in approximately one third of all diabetic patients. The molecular pathogenesis of diabetic nephropathy has not been fully characterized and novel mediators and drivers of the disease are still being described. Previous data from our laboratory has identified the developmentally regulated gene Gremlin as a novel target implicated in diabetic nephropathy in vitro and in vivo. We used bioinformatic analysis to examine whether Gremlin gene sequence and structure could be used to identify other genes implicated in diabetic nephropathy. The Notch ligand Jagged1 and its downstream effector, hairy enhancer of split-1 (Hes1), were identified as genes with significant similarity to Gremlin in terms of promoter structure and predicted microRNA binding elements. This led us to discover that transforming growth factor-beta (TGFbeta1), a primary driver of cellular changes in the kidney during nephropathy, increased Gremlin, Jagged1 and Hes1 expression in human kidney epithelial cells. Elevated levels of Gremlin, Jagged1 and Hes1 were also detected in extracts from renal biopsies from diabetic nephropathy patients, but not in control living donors. In situ hybridization identified specific upregulation and co-expression of Gremlin, Jagged1 and Hes1 in the same tubuli of kidneys from diabetic nephropathy patients, but not controls. Finally, Notch pathway gene clustering showed that samples from diabetic nephropathy patients grouped together, distinct from both control living donors and patients with minimal change disease. Together, these data suggest that Notch pathway gene expression is elevated in diabetic nephropathy, co-incident with Gremlin, and may contribute to the pathogenesis of this disease.     

Hyperglycosylation and reduced GABA currents of mutated GABRB3 polypeptide in remitting childhood absence epilepsy

Childhood absence epilepsy (CAE) accounts for 10% to 12% of epilepsy in children under 16 years of age. We screened for mutations in the GABA(A) receptor (GABAR) beta 3 subunit gene (GABRB3) in 48 probands and families with remitting CAE. We found that four out of 48 families (8%) had mutations in GABRB3. One heterozygous missense mutation (P11S) in exon 1a segregated with four CAE-affected persons in one multiplex, two-generation Mexican family. P11S was also found in a singleton from Mexico. Another heterozygous missense mutation (S15F) was present in a singleton from Honduras. An exon 2 heterozygous missense mutation (G32R) was present in two CAE-affected persons and two persons affected with EEG-recorded spike and/or sharp wave in a two-generation Honduran family. All mutations were absent in 630 controls. We studied functions and possible pathogenicity by expressing mutations in HeLa cells with the use of Western blots and an in vitro translation and translocation system. Expression levels did not differ from those of controls, but all mutations showed hyperglycosylation in the in vitro translation and translocation system with canine microsomes. Functional analysis of human GABA(A) receptors (alpha 1 beta 3-v2 gamma 2S, alpha 1 beta 3-v2[P11S]gamma 2S, alpha 1 beta 3-v2[S15F]gamma 2S, and alpha 1 beta 3-v2[G32R]gamma 2S) transiently expressed in HEK293T cells with the use of rapid agonist application showed that each amino acid transversion in the beta 3-v2 subunit (P11S, S15F, and G32R) reduced GABA-evoked current density from whole cells. Mutated beta 3 subunit protein could thus cause absence seizures through a gain in glycosylation of mutated exon 1a and exon 2, affecting maturation and trafficking of GABAR from endoplasmic reticulum to cell surface and resulting in reduced GABA-evoked currents.     

Cursive word recognition based on interactive activation and early visual processing models

We present an off-line cursive word recognition system based completely on neural networks: reading models and models of early visual processing. The first stage (normalization) preprocesses the input image in order to reduce letter position uncertainty; the second stage (feature extraction) is based on the feedforward model of orientation selectivity; the third stage (letter pre-recognition) is based on a convolutional neural network, and the last stage (word recognition) is based on the interactive activation model.     

Determination of content of metallothionein and low molecular mass stress peptides in transgenic tobacco plants

Phytoremediation is a process that utilizes plants to remove, transfer, stabilize, or destroy pollutants in soil, sediment, and groundwater. Plants used for such purposes have several requirements. Genetic engineering these plants could be an effective tool used to acquire features needed for such purposes within a substantial amount of time. This paper aims to utilize electrochemical techniques to analyze transgenic tobacco and, thus, to reveal their heavy metals phytoremediation potential. Total thiol and metallothionein (MT) quantities were determined in the control and transgenic tobacco plants. The total content of thiols in transgenic plants varied within the range of 561 to 1,671 μg g⁻¹. Furthermore, the determination of MT was done on transgenic tobacco plants. The level of human MT in transgenic tobacco plants varied between 25 and 95 μg g⁻¹. However, a plant cell protects itself by synthesizing low molecular mass thiols such as reduced glutathione and phytochelatins to protect itself against heavy metals toxicity. The most important thiols, cysteine (Cys), glutathione (GSH), oxidised glutathione (GSSG) and phytochelatin 2 (PC2), were determined in the non-transgenic and transgenic tobacco plants by high performance liquid chromatography with electrochemical detection. Tobacco plants synthesizing the highest amount of metallothionein have the highest basal level of phytochelatin 2 as well as reduced glutathione and free cysteine. It clearly follows from the results obtained that the biosynthesis of particular thiols is mutually linked, which contributes to a better protection of a transgenic plant against heavy metals effects.