Gap Junctional Communication Modulates Gene Expression in Osteoblastic Cells

Bone-forming cells are organized in a multicellular network interconnected by gap junctions. In these cells, gap junctions are formed by connexin43 (Cx43) and connexin45 (Cx45). Cx43 gap junctions form pores that are more permeable to negatively charged dyes such as Lucifer yellow and calcein than are Cx45 pores. We studied whether altering gap junctional communication by manipulating the relative expression of Cx43 and Cx45 affects the osteoblast phenotype. Transfection of Cx45 in cells that express primarily Cx43 (ROS 17/2.8 and MC3T3-E1) decreased both dye transfer and expression of osteocalcin (OC) and bone sialoprotein (BSP), genes pivotal to bone matrix formation and calcification. Conversely, transfection of Cx43 into cells that express predominantly Cx45 (UMR 106–01) increased both cell coupling and expression of OC and BSP. Transient cotransfection of promoter–luciferase constructs and connexin expression vectors demonstrated that OC and BSP gene transcription was down-regulated by Cx45 cotransfection in ROS 17/2.8 and MC3T3-E1 cells, in association with a decrease in dye coupling. Conversely, cotransfection of Cx43 in UMR 106–01 cells up-regulated OC and BSP gene transcription. Activity of other less specific osteoblast promoters, such as osteopontin and osteonectin, was less sensitive to changes in gap junctional communication. Thus, altering gap junctional permeability by manipulating the expression of Cx43 and Cx45 in osteoblastic cells alters transcriptional activity of osteoblast-specific promoters, presumably via modulation of signals that can diffuse from cell to cell. A communicating intercellular network is required for the full elaboration of a differentiated osteoblastic phenotype.

Receptor-mediated gene transfer into macrophages.

Gene transfer systems targeting various receptors have been developed to introduce functional genes into cells in culture and into intact animals. A synthetic molecular conjugate, consisting of mannosylated polylysine that exploits endocytosis via the macrophage mannose receptor, was constructed and complexed to expression plasmids containing either the Photinus pyralis luciferase or Escherichia coli beta-galactosidase (lacZ) reporter genes. The DNA complexes were used to transfect murine macrophages isolated from peritoneal exudates in vitro. Luciferase and beta-galactosidase activity was found in transfected cells in culture, whereas complexes consisting of an irrelevant plasmid bound to mannosylated polylysine or the expression plasmid bound to galactosylated polylysine resulted in no detectable transgene expression. Gene transfer was inhibited by the addition of excess mannosylated bovine serum albumin to the culture medium before transfection. Reporter genes were also transferred into macrophages residing in the spleen and liver of adult animals using this system. Luciferase activity was maximal at 4 days after transfection and decreased to lower levels by 16 days. Transgene expression conformed to the distribution of cells that had nonspecific esterase, a cytochemical marker for macrophages. Thus, this system can be used to introduce functional genes into macrophages and may be an approach to the treatment of storage diseases that affect the reticuloendothelial system.

From seed germination to flowering, light controls plant development via the pigment phytochrome.

Plant growth and development are regulated by interactions between the environment and endogenous developmental programs. Of the various environmental factors controlling plant development, light plays an especially important role, in photosynthesis, in seasonal and diurnal time sensing, and as a cue for altering developmental pattern. Recently, several laboratories have devised a variety of genetic screens using Arabidopsis thaliana to dissect the signal transduction pathways of the various photoreceptor systems. Genetic analysis demonstrates that light responses are not simply endpoints of linear signal transduction pathways but are the result of the integration of information from a variety of photoreceptors through a complex network of interacting signaling components. These signaling components include the red/far-red light receptors, phytochromes, at least one blue light receptor, and negative regulatory genes (DET, COP, and FUS) that act downstream from the photoreceptors in the nucleus. In addition, a steroid hormone, brassinolide, also plays a role in light-regulated development and gene expression in Arabidopsis. These molecular and genetic data are allowing us to construct models of the mechanisms by which light controls development and gene expression in Arabidopsis. In the future, this knowledge can be used as a framework for understanding how all land plants respond to changes in their environment.

Distinct roles of the NH2- and COOH-terminal domains of the protein inhibitor of activated signal transducer and activator of transcription (STAT) 1 (PIAS1) in cytokine-induced PIAS1–Stat1 interaction

STATs are activated by tyrosine phosphorylation on cytokine stimulation. A tyrosine-phosphorylated STAT forms a functional dimer through reciprocal Src homology 2 domain (SH2)–phosphotyrosyl peptide interactions. IFN treatment induces the association of PIAS1 and Stat1, which results in the inhibition of Stat1-mediated gene activation. The molecular basis of the cytokine-dependent PIAS1–Stat1 interaction has not been understood. We report here that a region near the COOH terminus of PIAS1 (amino acids 392–541) directly interacts with the NH2-terminal domain of Stat1 (amino acids 1–191). A mutant PIAS1 lacking the Stat1-interacting domain failed to inhibit Stat1-mediated gene activation. By using a modified yeast two-hybrid assay, we demonstrated that PIAS1 specifically interacts with the Stat1 dimer, but not tyrosine-phosphorylated or -unphosphorylated Stat1 monomer. In addition, whereas the NH2-terminal region of PIAS1 does not interact with Stat1, it serves as a modulatory domain by preventing the interaction of the COOH-terminal domain of PIAS1 with the Stat1 monomer. Thus, the cytokine-induced PIAS1–Stat1 interaction is mediated through the specific recognition of the dimeric form of Stat1 by PIAS1.

Identification of proacrosin binding protein sp32 precursor as a human cancer/testis antigen

Serological expression cloning of antigens eliciting a humoral immune response to a syngeneic mouse sarcoma identified pem (mouse placenta and embryonic expression gene) as a new member of the cancer/testis family. To identify the human homologue of pem, mouse pem sequences and pem-related expressed sequence tags from human testis were used as PCR primers for amplification using human testis cDNA. However, rather than pem, another gene, designated OY-TES-1, was isolated and found to be the human homologue of proacrosin binding protein sp32 precursor originally identified in mouse, guinea pig, and pig. OY-TES-1 maps to chromosome 12p12-p13 and contains 10 exons. Southern blot analysis suggests the presence of two OY-TES-1-related genes in the human genome. In normal tissues, OY-TES-1 mRNA was expressed only in testis, whereas in malignant tissues, a variable proportion of a wide array of cancers, including bladder, breast, lung, liver, and colon cancers, expressed OY-TES-1. Serological survey of 362 cancer patients with a range of different cancers showed antibody to OY-TES-1 in 25 patients. No OY-TES-1 sera reactivity was found in 20 normal individuals. These findings indicate that OY-TES-1 is an additional member of the cancer/testis family of antigens and that OY-TES-1 is immunogenic in humans.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.

Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization

Tuberculosis is a chronic infectious disease that is transmitted by cough-propelled droplets that carry the etiologic bacterium, Mycobacterium tuberculosis. Although currently available drugs kill most isolates of M. tuberculosis, strains resistant to each of these have emerged, and multiply resistant strains are increasingly widespread. The growing problem of drug resistance combined with a global incidence of seven million new cases per year underscore the urgent need for new antituberculosis therapies. The recent publication of the complete sequence of the M. tuberculosis genome has made possible, for the first time, a comprehensive genomic approach to the biology of this organism and to the drug discovery process. We used a DNA microarray containing 97% of the ORFs predicted from this sequence to monitor changes in M. tuberculosis gene expression in response to the antituberculous drug isoniazid. Here we show that isoniazid induced several genes that encode proteins physiologically relevant to the drug’s mode of action, including an operonic cluster of five genes encoding type II fatty acid synthase enzymes and fbpC, which encodes trehalose dimycolyl transferase. Other genes, not apparently within directly affected biosynthetic pathways, also were induced. These genes, efpA, fadE23, fadE24, and ahpC, likely mediate processes that are linked to the toxic consequences of the drug. Insights gained from this approach may define new drug targets and suggest new methods for identifying compounds that inhibit those targets.

Knowledge-based analysis of microarray gene expression data by using support vector machines

We introduce a method of functionally classifying genes by using gene expression data from DNA microarray hybridization experiments. The method is based on the theory of support vector machines (SVMs). SVMs are considered a supervised computer learning method because they exploit prior knowledge of gene function to identify unknown genes of similar function from expression data. SVMs avoid several problems associated with unsupervised clustering methods, such as hierarchical clustering and self-organizing maps. SVMs have many mathematical features that make them attractive for gene expression analysis, including their flexibility in choosing a similarity function, sparseness of solution when dealing with large data sets, the ability to handle large feature spaces, and the ability to identify outliers. We test several SVMs that use different similarity metrics, as well as some other supervised learning methods, and find that the SVMs best identify sets of genes with a common function using expression data. Finally, we use SVMs to predict functional roles for uncharacterized yeast ORFs based on their expression data.

Differential gene expression in p53-mediated apoptosis-resistant vs. apoptosis-sensitive tumor cell lines

Induction of wild-type p53 in the ECV-304 bladder carcinoma cell line by infection with a p53 recombinant adenovirus (Ad5CMV-p53) resulted in extensive apoptosis and eventual death of nearly all of the cells. As a strategy to determine the molecular events important to p53-mediated apoptosis in these transformed cells, ECV-304 cells were selected for resistance to p53 by repeated infections with Ad5CMV-p53. We compared the expression of 5,730 genes in p53-resistant (DECV) and p53-sensitive ECV-304 cells by reverse transcription–PCR, Northern blotting, and DNA microarray analysis. The expression of 480 genes differed by 2-fold or more between the two p53-infected cell lines. A number of potential targets for p53 were identified that play roles in cell cycle regulation, DNA repair, redox control, cell adhesion, apoptosis, and differentiation. Proline oxidase, a mitochondrial enzyme involved in the proline/pyrroline-5-carboxylate redox cycle, was up-regulated by p53 in ECV but not in DECV cells. Pyrroline-5-carboxylate (P5C), a proline-derived metabolite generated by proline oxidase, inhibited the proliferation and survival of ECV-304 and DECV cells and induced apoptosis in both cell lines. A recombinant proline oxidase protein tagged with a green fluorescent protein at the amino terminus localized to mitochondria and induced apoptosis in p53-null H1299 non-small cell lung carcinoma cells. The results directly implicate proline oxidase and the proline/P5C pathway in p53-induced growth suppression and apoptosis.

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.

Comprehensive copy number and gene expression profiling of the 17q23 amplicon in human breast cancer

The biological significance of DNA amplification in cancer is thought to be due to the selection of increased expression of a single or few important genes. However, systematic surveys of the copy number and expression of all genes within an amplified region of the genome have not been performed. Here we have used a combination of molecular, genomic, and microarray technologies to identify target genes for 17q23, a common region of amplification in breast cancers with poor prognosis. Construction of a 4-Mb genomic contig made it possible to define two common regions of amplification in breast cancer cell lines. Analysis of 184 primary breast tumors by fluorescence in situ hybridization on tissue microarrays validated these results with the highest amplification frequency (12.5%) observed for the distal region. Based on GeneMap'99 information, 17 known genes and 26 expressed sequence tags were localized to the contig. Analysis of genomic sequence identified 77 additional transcripts. A comprehensive analysis of expression levels of these transcripts in six breast cancer cell lines was carried out by using complementary DNA microarrays. The expression patterns varied from one cell line to another, and several overexpressed genes were identified. Of these, RPS6KB1, MUL, APPBP2, and TRAP240 as well as one uncharacterized expressed sequence tag were located in the two common amplified regions. In summary, comprehensive analysis of the 17q23 amplicon revealed a limited number of highly expressed genes that may contribute to the more aggressive clinical course observed in breast cancer patients with 17q23-amplified tumors.

Recursive partitioning for tumor classification with gene expression microarray data

Precise classification of tumors is critically important for cancer diagnosis and treatment. It is also a scientifically challenging task. Recently, efforts have been made to use gene expression profiles to improve the precision of classification, with limited success. Using a published data set for purposes of comparison, we introduce a methodology based on classification trees and demonstrate that it is significantly more accurate for discriminating among distinct colon cancer tissues than other statistical approaches used heretofore. In addition, competing classification trees are displayed, which suggest that different genes may coregulate colon cancers.

Genome-wide gene expression profiles of the developing mouse hippocampus

We have analyzed the developmental molecular programs of the mouse hippocampus, a cortical structure critical for learning and memory, by means of large-scale DNA microarray techniques. Of 11,000 genes and expressed sequence tags examined, 1,926 showed dynamic changes during hippocampal development from embryonic day 16 to postnatal day 30. Gene-cluster analysis was used to group these genes into 16 distinct clusters with striking patterns that appear to correlate with major developmental hallmarks and cellular events. These include genes involved in neuronal proliferation, differentiation, and synapse formation. A complete list of the transcriptional changes has been compiled into a comprehensive gene profile database (http://BrainGenomics.Princeton.edu), which should prove valuable in advancing our understanding of the molecular and genetic programs underlying both the development and the functions of the mammalian brain.

Human carbonic anhydrase XII: cDNA cloning, expression, and chromosomal localization of a carbonic anhydrase gene that is overexpressed in some renal cell cancers

We report the cloning and characterization of a tumor-associated carbonic anhydrase (CA) that was identified in a human renal cell carcinoma (RCC) by serological expression screening with autologous antibodies. The cDNA sequence predicts a 354-amino acid polypeptide with a molecular mass of 39,448 Da that has features of a type I membrane protein. The predicted sequence includes a 29-amino acid signal sequence, a 261-amino acid CA domain, an additional short extracellular segment, a 26-amino acid hydrophobic transmembrane domain, and a hydrophilic C-terminal cytoplasmic tail of 29 amino acids that contains two potential phosphorylation sites. The extracellular CA domain shows 30–42% homology with known human CAs, contains all three Zn-binding histidine residues found in active CAs, and contains two potential sites for asparagine glycosylation. When expressed in COS cells, the cDNA produced a 43- to 44-kDa protein in membranes that had around one-sixth the CA activity of membranes from COS cells transfected with the same vector expressing bovine CA IV. We have designated this human protein CA XII. Northern blot analysis of normal tissues demonstrated a 4.5-kb transcript only in kidney and intestine. However, in 10% of patients with RCC, the CA XII transcript was expressed at much higher levels in the RCC than in surrounding normal kidney tissue. The CA XII gene was mapped by using fluorescence in situ hybridization to 15q22. CA XII is the second catalytically active membrane CA reported to be overexpressed in certain cancers. Its relationship to oncogenesis and its potential as a clinically useful tumor marker clearly merit further investigation.

Profound misregulation of muscle-specific gene expression in facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is a neuromuscular disorder characterized by an insidious onset and progressive course. The disease has a frequency of about 1 in 20,000 and is transmitted in an autosomal dominant fashion with almost complete penetrance. Deletion of an integral number of tandemly arrayed 3.3-kb repeat units (D4Z4) on chromosome 4q35 is associated with FSHD but otherwise the molecular basis of the disease and its pathophysiology remain obscure. Comparison of mRNA populations between appropriate cell types can facilitate identification of genes relevant to a particular biological or pathological process. In this report, we have compared mRNA populations of FSHD and normal muscle. Unexpectedly, the dystrophic muscle displayed profound alterations in gene expression characterized by severe underexpression or overexpression of specific mRNAs. Intriguingly, many of the deregulated mRNAs are muscle specific. Our results suggest that a global misregulation of gene expression is the underlying basis for FSHD, distinguishing it from other forms of muscular dystrophy. The experimental approach used here is applicable to any genetic disorder whose pathogenic mechanism is incompletely understood.