A ratiometric-based measure of gene co-expression

Current measures of global gene expression analyses, such as correlation and mutual information-based approaches, largely depend on the degree of association between mRNA levels and to a lesser extent on variability. I develop and implement a new approach, called the Ratiometric method, which is based on the coefficient of variation of the expression ratio of two genes, relying more on variation than previous methods. The advantage of such modus operandi is the ability to detect possible gene pair interactions regardless of the degree of expression dispersion across the sample group. Gene pairs with low expression dispersion, i.e., their absolute expressions remain constant across the sample group, are systematically missed by correlation and mutual information analyses. The superiority of the Ratiometric method in finding these gene pair interactions is demonstrated in a data set of RNA-seq B-cell samples from the 1000 Genomes Project Consortium. The Ratiometric method renders a more comprehensive recovery of KEGG pathways and GO-terms.

Surface markers of regionalization in the vertebrate nervous system

In order to identify new molecules that might play a role in regional specification of the nervous system, we generated and characterized monoclonal antibodies (mAbs) that have positionally-restricted labeling patterns.

The FORSE-1 mAb was generated using a strategy designed to produce mAbs against neuronal cell surface antigens that might be regulated by regionally-restricted transcription factors in the developing central nervous system (CNS). FORSE-1 staining is enriched in the forebrain as compared to the rest of the CNS until E18. Between E11.5-E13.5, only certain areas of the forebrain are labeled. There is also a dorsoventrally-restricted region of labeling in the hindbrain and spinal cord. The mAb labels a large proteoglycan-like cell-surface antigen (>200 kD). The labeling pattern of FORSE-1 is conserved in various mammals and in chick.

To determine whether the FORSE-1 labeling pattern is similar to that of known transcription factors, the expression of BF-1 and Dlx-2 was compared with FORSE-1. There is a striking overlap between BF-1 and FORSE-1 in the telencephalon. In contrast, FORSE-1 and Dlx-2 have very different patterns of expression in the forebrain, suggesting that regulation by Dlx-2 alone cannot explain the distribution of FORSE-1. They do, however, share some sharp boundaries in the diencephalon. In addition, FORSE-1 identifies some previously unknown boundaries in the developing forebrain. Thus, FORSE-1 is a new cell surface marker that can be used to subdivide the embryonic forebrain into regions smaller than previously described, providing further complexity necessary for developmental patterning.

I also studied the expression of the cell surface protein CD9 in the developing and adult rat nervous system. CD9 is implicated in intercellular signaling and cell adhesion in the hematopoetic system. In the nervous system, CD9 may perform similar functions in early sympathetic ganglia, chromaffin cells, and motor neurons, all of which express the protein. The presence of CD9 on the surfaces of Schwann cells and axons at the appropriate time may allow the protein to participate in the cellular interactions involved in myelination.

Characterization and developmental regulation of a gene expressed specifically in the skeletogenic lineage of the sea urchin embryo

The sea urchin embryonic skeleton, or spicule, is deposited by mesenchymal progeny of four precursor cells, the micromeres, which are determined to the skeletogenic pathway by a process known as cytoplasmic localization. A gene encoding one of the major products of the skeletogenic mesenchyme, a prominent 50 kD protein of the spicule matrix, has been characterized in detail. cDNA clones were first isolated by antibody screening of a phage expression library, followed by isolation of homologous genomic clones. The gene, known as SM50, is single copy in the sea urchin genome, is divided into two exons of 213 and 1682 bp, and is expressed only in skeletogenic cells. Transcripts are first detectable at the 120 cell stage, shortly after the segregation of the skeletogenic precursors from the rest of the embryo. The SM50 open reading frame begins within the first exon, is 450 amino acids in length, and contains a loosely repeated 13 amino acid motif rich in acidic residues which accounts for 45% of the protein and which is possibly involved in interaction with the mineral phase of the spicule.

The important cis-acting regions of the SM50 gene necessary for proper regulation of expression were identified by gene transfer experiments. A 562 bp promoter fragment, containing 438 bp of 5' promoter sequence and 124 bp of the SM50 first exon (including the SM50 initiation codon), was both necessary and sufficient to direct high levels of expression of the bacterial chloramphenicol acetyltransferase (CAT) reporter gene specifically in the skeletogenic cells. Removal of promoter sequences between positions -2200 and -438, and of transcribed regions downstream of +124 (including the SM50 intron), had no effect on the spatial or transcriptional activity of the transgenes.

Regulatory proteins that interact with the SM50 promoter were identified by the gel retardation assay, using bulk embryo mesenchyme blastula stage nuclear proteins. Five protein binding sites were identified and mapped to various degrees of resolution. Two sites are homologous, may be enhancer elements, and at least one is required for expression. Two additional sites are also present in the promoter of the aboral ectoderm specific cytoskeletal actin gene CyIIIa; one of these is a CCAA T element, the other a putative repressor element. The fifth site overlaps the binding site of the putative repressor and may function as a positive regulator by interfering with binding of the repressor. All of the proteins are detectable in nuclear extracts prepared from 64 cell stage embryos, a stage just before expression of SM50 is initiated, as well as from blastula and gastrula stage; the putative enhancer binding protein may be maternal as well.

A Detailed Analysis of Transcriptional Regulators Affecting the Saccharomyces cerevisiae Heat-shock Gene SSA1

The yeast Saccharomyces cerevisiae contains a family of hsp70 related genes. One member of this family, SSA1, encodes a 70kD heat-shock protein which in addition to its heat inducible expression has a significant basal level of expression. The first 500 bp upstream of the SSA1 start point of transcription was examined by DNAse I protection analysis. The results reveal the presence of at least 14 factor binding sites throughout the upstream promoter region. The function of these binding sites has been examined using a series of 5' promoter deletions fused to the recorder gene lacZ in a centromere-containing yeast shuttle vector. The following sites have been identified in the promoter and their activity in yeast determined individually with a centromere-based recorder plasmid containing a truncated CYC1 /lacZ fusion: a heat-shock element or HSE which is sufficient to convey heat-shock response on the recorder plasmid; a homology to the SV40 'core' sequence which can repress the GCN4 recognition element (GCRE) and the yAP1 recognition element (ARE), and has been designated a upstream repression element or URE; a 'G'-rich region named G-box which can also convey heatshock response on the recorder plasmid; and a purine-pyrimidine alternating sequence name GT-box which is an activator of transcription. A series of fusion constructs were made to identify a putative silencer-like element upstream of SSA1. This element is position dependent and has been localized to a region containing both an ABF1 binding site and a RAP1 binding site. Five site-specific DNA-binding factors are identified and their purification is presented: the heat-shock transcription factor or HSTF, which recognizes the HSE; the G-box binding factor or GBF; the URE recognition factor or URF; the GT-box binding factor; and the GC-box binding factor or yeast Sp1.