Additivity in protein-DNA interactions: how good an approximation is it?

Man and Stormo and Bulyk et al. recently presented their results on the study of the DNA binding affinity of proteins. In both of these studies the main conclusion is that the additivity assumption, usually applied in methods to search for binding sites, is not true. In the first study, the analysis of binding affinity data from the Mnt repressor protein bound to all possible DNA (sub)targets at positions 16 and 17 of the binding site, showed that those positions are not independent. In the second study, the authors analysed DNA binding affinity data of the wild-type mouse EGR1 protein and four variants differing on the middle finger. The binding affinity of these proteins was measured to all 64 possible trinucleotide (sub)targets of the middle finger using microarray technology. The analysis of the measurements also showed interdependence among the positions in the DNA target. In the present report, we review the data of both studies and we re- analyse them using various statistical methods, including a comparison with a multiple regression approach. We conclude that despite the fact that the additivity assumption does not fit the data perfectly, in most cases it provides a very good approximation of the true nature of the specific protein–DNA interactions. Therefore, additive models can be very useful for the discovery and prediction of binding sites in genomic DNA.     

Evolution of the 28S ribosomal RNA gene in anurans: regions of variability and their phylogenetic implications

Fifteen restriction sites were mapped to the 28S ribosomal RNA gene of individuals representing 54 species of frogs, two species of salamanders, a caecilian, and a lungfish. Eight of these sites were present in all species examined, and two were found in all but one species. Alignment of these conserved restriction sites revealed, among anuran 28S rRNA genes, five regions of major length variation that correspond to four of 12 previously identified divergent domains of this gene. One of the divergent domains (DD8) consists of two regions of length variation separated by a short segment that is conserved at least throughout tetrapods. Most of the insertions, deletions, and restriction-site variations identified in the 28S gene will require sequence-level analysis for a detailed reconstruction of their history. However, an insertion in DD9 that is coextensive with frogs in the suborder Neobatrachia, a BstEII site that is limited to representatives of two leptodactylid subfamilies, and a deletion in DD10 that is found only in three ranoid genera are probably synapomorphies.      

Analytical validation of quantitative immunohistochemical assays of tumor infiltrating lymphocyte biomarkers

Tumor infiltrating lymphocytes (TIL), especially T-cells, have both prognostic and therapeutic applications. The presence of CD8+ effector T-cells and the ratio of CD8+ cells to FOXP3+ regulatory T-cells have been used as biomarkers of disease prognosis to predict response to various immunotherapies. Blocking the interaction between inhibitory receptors on T-cells and their ligands with therapeutic antibodies including atezolizumab, nivolumab, pembrolizumab and tremelimumab increases the immune response against cancer cells and has shown significant improvement in clinical benefits and survival in several different tumor types. The improved clinical outcome is presumed to be associated with a higher tumor infiltration; therefore, it is thought that more accurate methods for measuring the amount of TIL could assist prognosis and predict treatment response. We have developed and validated quantitative immunohistochemistry (IHC) assays for CD3, CD8 and FOXP3 for immunophenotyping T-lymphocytes in tumor tissue. Various types of formalin fixed, paraffin embedded (FFPE) tumor tissues were immunolabeled with anti-CD3, anti-CD8 and anti-FOXP3 antibodies using an IHC autostainer. The tumor area of stained tissues, including the invasive margin of the tumor, was scored by a pathologist (visual scoring) and by computer-based quantitative image analysis. Two image analysis scores were obtained for the staining of each biomarker: the percent positive cells in the tumor area and positive cells/mm² tumor area. Comparison of visual vs. image analysis scoring methods using regression analysis showed high correlation and indicated that quantitative image analysis can be used to score the number of positive cells in IHC stained slides. To demonstrate that the IHC assays produce consistent results in normal daily testing, we evaluated the specificity, sensitivity and reproducibility of the IHC assays using both visual and image analysis scoring methods. We found that CD3, CD8 and FOXP3 IHC assays met the fit-for-purpose analytical acceptance validation criteria and that they can be used to support clinical studies.     

A global approach to identify differentially expressed genes in cDNA (two-color) microarray experiments

MOTIVATION: Currently most of the methods for identifying differentially expressed genes fall into the category of so called single-gene-analysis, performing hypothesis testing on a gene-by-gene basis. In a single-gene-analysis approach, estimating the variability of each gene is required to determine whether a gene is differentially expressed or not. Poor accuracy of variability estimation makes it difficult to identify genes with small fold-changes unless a very large number of replicate experiments are performed. RESULTS: We propose a method that can avoid the difficult task of estimating variability for each gene, while reliably identifying a group of differentially expressed genes with low false discovery rates, even when the fold-changes are very small. In this article, a new characterization of differentially expressed genes is established based on a theorem about the distribution of ranks of genes sorted by (log) ratios within each array. This characterization of differentially expressed genes based on rank is an example of all-gene-analysis instead of single gene analysis. We apply the method to a cDNA microarray dataset and many low fold-changed genes (as low as 1.3 fold-changes) are reliably identified without carrying out hypothesis testing on a gene-by-gene basis. The false discovery rate is estimated in two different ways reflecting the variability from all the genes without the complications related to multiple hypothesis testing. We also provide some comparisons between our approach and single-gene-analysis based methods. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Selection of optimal DNA oligos for gene expression arrays.

MOTIVATION: High density DNA oligo microarrays are widely used in biomedical research. Selection of optimal DNA oligos that are deposited on the microarrays is critical. Based on sequence information and hybridization free energy, we developed a new algorithm to select optimal short (20-25 bases) or long (50 or 70 bases) oligos from genes or open reading frames (ORFs) and predict their hybridization behavior. Having optimized probes for each gene is valuable for two reasons. By minimizing background hybridization they provide more accurate determinations of true expression levels. Having optimum probes minimizes the number of probes needed per gene, thereby decreasing the cost of each microarray, raising the number of genes on each chip and increasing its usage. RESULTS: In this paper we describe algorithms to optimize the selection of specific probes for each gene in an entire genome. The criteria for truly optimum probes are easily stated but they are not computable at all levels currently. We have developed an heuristic approach that is efficiently computable at all levels and should provide a good approximation to the true optimum set. We have run the program on the complete genomes for several model organisms and deposited the results in a database that is available on-line (http://ural.wustl.edu/~lif/probe.pl). AVAILABILITY: The program is available upon request.     

The induction of allophanate lyase during the vegetative cell cycle in light-synchronized cultures of Chlamydomonas reinhardi.

Allophanate lyase can be induced by urea or acetamide 20-40-fold within 4 h in NH4 + -deprived cultures of Chlamydomonas reinhardi. In light-synchronized cultures, allophanate lyase induction appeared to be limited to the light phase of the cell cycle, provided that culture samples were induced under ongoing illumination conditions (i.e. light induction of light phase cells and dark induction of dark phase cells). However, when culture samples were induced under constant light conditions this cell cycle pattern was abolished. Light was found to be required for allophanate lyase induction and this was shown to be due, in part, to the light requirement for inducer uptake. The relationship between allophanate lyase induction and gametogenesis is discussed.     

Preparation of Encapsulated Microbial Cells for Environmental Applications

An improved method for the encapsulation of bacteria into microspheres of alginate, agarose, or polyurethane is described. Cell suspensions were passed through a low-pressure nozzle into an aqueous phase where matrix polymerization or gelation yielded beads 2 to 50 μm in diameter. Trials with a chlorophenol-degrading Flavobacterium species showed that cells entrapped by these procedures were as catabolically active as free cells. These types of beads should have numerous applications in the fields of environmental science and engineering.