Tp53-dependent repression of cell cycle target genes: Global and mechanistic analysis

Most studies of p53 function have focused on genes transactivated by p53. It is less widely appreciated that p53 can repress target genes to affect a particular cellular response. There is evidence that repression is important for p53-induced apoptosis and cell cycle arrest. It is less clear if repression is important for other p53 functions. A comprehensive knowledge of the genes repressed by p53 and the cellular processes they affect is currently lacking. We used an expression profiling strategy to identify p53-responsive genes following adenoviral p53 gene transfer (Ad-p53) in PC3 prostate cancer cells. A total of 111 genes represented on the Affymetrix U133A microarray were repressed more than two fold (p ≤ 0.05) by p53. An objective assessment of array data quality was carried out using RT-PCR of 20 randomly selected genes. We estimate a confirmation rate of >95.5% for the complete data set. Functional over-representation analysis was used to identify cellular processes potentially affected by p53-mediated repression. Cell cycle regulatory genes exhibited significant enrichment (p ≤ 5E-28) within the repressed targets. Several of these genes are repressed in a p53-dependent manner following DNA damage, but preceding cell cycle arrest. These findings identify novel p53-repressed targets and indicate that p53-induced cell cycle arrest is a function of not only the transactivation of cell cycle inhibitors (e.g., p21), but also the repression of targets that act at each phase of the cell cycle. The mechanism of repression of this set of p53 targets was investigated. Most of the repressed genes identified here do not harbor consensus p53 DNA binding sites but do contain binding sites for E2F transcription factors. We demonstrate a role for E2F/RB repressor complexes in our system. Importantly, p53 is found at the promoter of CDC25A. CDC25A protein is rapidly degraded in response to DNA damage. Our group has demonstrated for the first time that CDC25A is also repressed at the transcript level by p53. This work has important implications for understanding the DNA damage cell cycle checkpoint response and the link between E2F/RB complexes and p53 in the repression of target genes. ^

Nonlinear dimension reduction in pathway-based genome-wide association analysis

Pathway based genome wide association study evolves from pathway analysis for microarray gene expression and is under rapid development as a complementary for single-SNP based genome wide association study. However, it faces new challenges, such as the summarization of SNP statistics to pathway statistics. The current study applies the ridge regularized Kernel Sliced Inverse Regression (KSIR) to achieve dimension reduction and compared this method to the other two widely used methods, the minimal-p-value (minP) approach of assigning the best test statistics of all SNPs in each pathway as the statistics of the pathway and the principal component analysis (PCA) method of utilizing PCA to calculate the principal components of each pathway. Comparison of the three methods using simulated datasets consisting of 500 cases, 500 controls and100 SNPs demonstrated that KSIR method outperformed the other two methods in terms of causal pathway ranking and the statistical power. PCA method showed similar performance as the minP method. KSIR method also showed a better performance over the other two methods in analyzing a real dataset, the WTCCC Ulcerative Colitis dataset consisting of 1762 cases, 3773 controls as the discovery cohort and 591 cases, 1639 controls as the replication cohort. Several immune and non-immune pathways relevant to ulcerative colitis were identified by these methods. Results from the current study provided a reference for further methodology development and identified novel pathways that may be of importance to the development of ulcerative colitis.^

Patenting genes and gene sequences: Analysis of public health implications and the law

At issue is whether or not isolated DNA is patent eligible under the U.S. Patent Law and the implications of that determination on public health. The U.S. Patent and Trademark Office has issued patents on DNA since the 1980s, and scientists and researchers have proceeded under that milieu since that time. Today, genetic research and testing related to the human breast cancer genes BRCA1 and BRCA2 is conducted within the framework of seven patents that were issued to Myriad Genetics and the University of Utah Research Foundation between 1997 and 2000. In 2009, suit was filed on behalf of multiple researchers, professional associations and others to invalidate fifteen of the claims underlying those patents. The Court of Appeals for the Federal Circuit, which hears patent cases, has invalidated claims for analyzing and comparing isolated DNA but has upheld claims to isolated DNA. The specific issue of whether isolated DNA is patent eligible is now before the Supreme Court, which is expected to decide the case by year's end. In this work, a systematic review was performed to determine the effects of DNA patents on various stakeholders and, ultimately, on public health; and to provide a legal analysis of the patent eligibility of isolated DNA and the likely outcome of the Supreme Court's decision. ^ A literature review was conducted to: first, identify principle stakeholders with an interest in patent eligibility of the isolated DNA sequences BRCA1 and BRCA2; and second, determine the effect of the case on those stakeholders. Published reports that addressed gene patents, the Myriad litigation, and implications of gene patents on stakeholders were included. Next, an in-depth legal analysis of the patent eligibility of isolated DNA and methods for analyzing it was performed pursuant to accepted methods of legal research and analysis based on legal briefs, federal law and jurisprudence, scholarly works and standard practice legal analysis. ^ Biotechnology, biomedical and clinical research, access to health care, and personalized medicine were identified as the principle stakeholders and interests herein. Many experts believe that the patent eligibility of isolated DNA will not greatly affect the biotechnology industry insofar as genetic testing is concerned; unlike for therapeutics, genetic testing does not require tremendous resources or lead time. The actual impact on biomedical researchers is uncertain, with greater impact expected for researchers whose work is intended for commercial purposes (versus basic science). The impact on access to health care has been surprisingly difficult to assess; while invalidating gene patents might be expected to decrease the cost of genetic testing and improve access to more laboratories and physicians' offices that provide the test, a 2010 study on the actual impact was inconclusive. As for personalized medicine, many experts believe that the availability of personalized medicine is ultimately a public policy issue for Congress, not the courts. ^ Based on the legal analysis performed in this work, this writer believes the Supreme Court is likely to invalidate patents on isolated DNA whose sequences are found in nature, because these gene sequences are a basic tool of scientific and technologic work and patents on isolated DNA would unduly inhibit their future use. Patents on complementary DNA (cDNA) are expected to stand, however, based on the human intervention required to craft cDNA and the product's distinction from the DNA found in nature. ^ In the end, the solution as to how to address gene patents may lie not in jurisprudence but in a fundamental change in business practices to provide expanded licenses to better address the interests of the several stakeholders. ^

Cloning and characterization of the full length cDNA and promoter of mouse tissue transglutaminase

Transglutaminases are a family of enzymes that catalyze the covalent cross-linking of proteins through the formation of $\varepsilon$-($\gamma$-glutaminyl)-lysyl isopeptide bonds. Tissue transglutaminase (Tgase) is an intracellular enzyme which is expressed in terminally differentiated and senescent cells and also in cells undergoing apoptotic cell death. To characterize this enzyme and examine its relationship with other members of the transglutaminase family, cDNAs, the first two exons of the gene and 2 kb of the 5$\sp\prime$ flanking region, including the promoter, were isolated. The full length Tgase transcript consists of 66 bp of 5$\sp\prime$-UTR (untranslated) sequence, an open reading frame which encodes 686 amino acids and 1400 bp of 3$\sp\prime$-UTR sequence. Alignment of the deduced Tgase protein sequence with that of other transglutaminases revealed regions of strong homology, particularly in the active site region.^ The Tgase cDNA was used to isolate and characterize a genomic clone encompassing the 5$\sp\prime$ end of the mouse Tgase gene. The transcription start site was defined using genomic and cDNA clones coupled with S1 protection analysis and anchored PCR. This clone includes 2.3 kb upstream of the transcription start site and two exons that contain the first 256 nucleotides of the mouse Tgase cDNA sequence. The exon intron boundaries have been mapped and compared with the exon intron boundaries of three members of the transglutaminase family: human factor XIIIa, the human keratinocyte transglutaminase and human erythrocyte band 4.1. Tissue Tgase exon II is similar to comparable exons of these genes. However, exon I bears no resemblance with any of the other transglutaminase amino terminus exons.^ Previous work in our laboratory has shown that the transcription of the Tgase gene is directly controlled by retinoic acid and retinoic acid receptors. To identify the region of the Tgase gene responsible for regulating its expression, fragments of the Tgase promoter and 5$\sp\prime$-flanking region were cloned into the chloramphenicol actetyl transferase (CAT) reporter constructs. Transient transfection experiments with these constructs demonstrated that the upstream region of Tgase is a functional promoter which contains a retinoid response element within a 1573 nucleotide region spanning nucleotides $-$252 to $-$1825. ^

Functional analysis of the sea urchin {\it orthodenticle\/} -related gene, {\it SPOTX\/}, in aboral ectoderm-specific gene activation and aboral ectoderm cell fate specification

An important question in developmental biology is how embryonic cell types are derived from a fertilized egg. To address this question, this thesis investigates the mechanisms by which the aboral ectoderm-specific Spec2a gene is spatially and temporally regulated during sea urchin embryogenesis. The Spec2a gene of the sea urchin Strongylocentratus purpuratus has served as a valuable maker to understand the basis of lineage-specific gene activation and the role of transcription factors in cell fate specification. The hypothesis is that transcription factors responsible for cell type-specific gene activation are key components in the initial cell specification step. The Spec2a gene, which encodes a small cytosolic calcium-binding protein, is expressed exclusively in aboral ectoderm cell lineages. The 1516-bp control region of the Spec2a gene contains a 188-bp enhancer element required for temporal activation and aboral ectoderm/mesenchyme cell expression, while an unidentified element upstream of the enhancer represses expression in mesenchyme cells. Using an enhancer activation assay, combined with site-directed mutagenesis, I showed that three TAATCC/T sites within the enhancer are responsible for enhancer activity. Mutagenizing these sites and a fourth one just upstream abolished all activity from the Spec2a control region. A 77-bp DNA fragment from the Spec2a enhancer containing two of the TAATCC/T sites is sufficient for aboral ectoderm/mesenchyme cell expression. A cDNA encoding SpOtx, an orthodenticle-related protein, was cloned from S. purpuratus and shown to bind with high affinity to the TAATCC/T sequences within the Spec2a control region. SpOtx transcripts were found initially in all cells of the cleaving embryo, but they gradually became restricted to oral ectoderm and endoderm cells, suggesting that SpOtx might play a role in the initial temporal activation of the Spec2a gene and most likely has additional functions in the developing embryo. To reveal the broader biological functions of SpOtx, I injected SpOtx mRNA into living sea urchin eggs to determine what effects overexpressing the SpOtx protein might have on embryo development. SpOtx mRNA-injected embryos displayed dramatic alterations in development. Instead of developing into pluteus larvae with 15 different cell types, uniform epithelia balls were formed. These balls consisted of a thin layer of squamous cells with short cilia highly reminiscent of aboral ectoderm. Immunohistochemical staining and RT-PCR demonstrated that the SpOtx-injected embryoids expressed aboral ectoderm markers uniformly, but showed very weak or no expression of markers for non-aboral ectoderm cell types. These data strongly suggested that overexpression of SpOtx redirected the normal fate of non-aboral ectoderm cells to that of aboral ectoderm. These results show that SpOtx is involved in aboral ectoderm differentiation by activating aboral ectoderm-specific genes and that modulating its expression can lead to changes in cell fate. ^

Two way clustering of Microarray Data using a Hybrid Approach

The Microarray technique is rather powerful, as it allows to test up thousands of genes at a time, but this produces an overwhelming set of data files containing huge amounts of data, which is quite difficult to pre-process, separate, classify and correlate for interesting conclusions to be extracted. Modern machine learning, data mining and clustering techniques based on information theory, are needed to read and interpret the information contents buried in those large data sets. Independent Component Analysis method can be used to correct the data affected by corruption processes or to filter the uncorrectable one and then clustering methods can group similar genes or classify samples. In this paper a hybrid approach is used to obtain a two way unsupervised clustering for a corrected microarray data.

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.

Transcripts from a single full-length cDNA clone of hepatitis C virus are infectious when directly transfected into the liver of a chimpanzee

We have succeeded in constructing a stable full-length cDNA clone of strain H77 (genotype 1a) of hepatitis C virus (HCV). We devised a cassette vector with fixed 5′ and 3′ termini and constructed multiple full-length cDNA clones of H77 in a single step by cloning of the entire ORF, which was amplified by long reverse transcriptase–PCR, directly into this vector. The infectivity of two complete full-length cDNA clones was tested by the direct intrahepatic injection of a chimpanzee with RNA transcripts. However, we found no evidence for HCV replication. Sequence analysis of these and 16 additional full-length clones revealed that seven clones were defective for polyprotein synthesis, and the remaining nine clones had 6–28 amino acid mutations in the predicted polyprotein compared with the consensus sequence of H77. Next, we constructed a consensus chimera from four of the full-length cDNA clones with just two ligation steps. Injection of RNA transcripts from this consensus clone into the liver of a chimpanzee resulted in viral replication. The sequence of the virus recovered from the chimpanzee was identical to that of the injected RNA transcripts. This stable infectious molecular clone should be an important tool for developing a better understanding of the molecular biology and pathogenesis of HCV.

Lactone-ring-cleaving enzyme: Genetic analysis, novel RNA editing, and evolutionary implications

A lactonohydrolase from Fusarium oxysporum AKU 3702 is an enzyme catalyzing the hydrolysis of aldonate lactones to the corresponding aldonic acids. The amino acid sequences of the NH2 terminus and internal peptide fragments of the enzyme were determined to prepare synthetic oligonucleotides as primers for the PCR. An approximate 1,000-base genomic DNA fragment thus amplified was used as the probe to clone both genomic DNA and cDNA for the enzyme. The lactonohydrolase genomic gene consists of six exons separated by five short introns. A novel type of RNA editing, in which lactonohydrolase mRNA included the insertion of guanosine and cytidine residues, was observed. The predicted amino acid sequence of the cloned lactonohydrolase cDNA showed significant similarity to those of the gluconolactonase from Zymomonas mobilis, and paraoxonases from human and rabbit, forming a unique superfamily consisting of C-O cleaving enzymes and P-O cleaving enzymes. Lactonohydrolase was expressed under the control of the lac promoter in Escherichia coli.

Gene discovery in the wood-forming tissues of poplar: Analysis of 5,692 expressed sequence tags

A rapidly growing area of genome research is the generation of expressed sequence tags (ESTs) in which large numbers of randomly selected cDNA clones are partially sequenced. The collection of ESTs reflects the level and complexity of gene expression in the sampled tissue. To date, the majority of plant ESTs are from nonwoody plants such as Arabidopsis, Brassica, maize, and rice. Here, we present a large-scale production of ESTs from the wood-forming tissues of two poplars, Populus tremula L. × tremuloides Michx. and Populus trichocarpa ‘Trichobel.’ The 5,692 ESTs analyzed represented a total of 3,719 unique transcripts for the two cDNA libraries. Putative functions could be assigned to 2,245 of these transcripts that corresponded to 820 protein functions. Of specific interest to forest biotechnology are the 4% of ESTs involved in various processes of cell wall formation, such as lignin and cellulose synthesis, 5% similar to developmental regulators and members of known signal transduction pathways, and 2% involved in hormone biosynthesis. An additional 12% of the ESTs showed no significant similarity to any other DNA or protein sequences in existing databases. The absence of these sequences from public databases may indicate a specific role for these proteins in wood formation. The cDNA libraries and the accompanying database are valuable resources for forest research directed toward understanding the genetic control of wood formation and future endeavors to modify wood and fiber properties for industrial use.

Chromosomal aberrations in PARP−/− mice: Genome stabilization in immortalized cells by reintroduction of poly(ADP-ribose) polymerase cDNA

Depletion of poly(ADP-ribose) polymerase (PARP) increases the frequency of recombination, gene amplification, sister chromatid exchanges, and micronuclei formation in cells exposed to genotoxic agents, implicating PARP in the maintenance of genomic stability. Flow cytometric analysis now has revealed an unstable tetraploid population in immortalized fibroblasts derived from PARP−/− mice. Comparative genomic hybridization detected partial chromosomal gains in 4C5-ter, 5F-ter, and 14A1-C1 in PARP−/−mice and immortalized PARP−/−fibroblasts. Neither the chromosomal gains nor the tetraploid population were apparent in PARP−/− cells stably transfected with PARP cDNA [PARP−/−(+PARP)], indicating negative selection of cells with these genetic aberrations after reintroduction of PARP cDNA. Although the tumor suppressor p53 was not detectable in PARP−/− cells, p53 expression was partially restored in PARP−/− (+PARP) cells. Loss of 14D3-ter that encompasses the tumor suppressor gene Rb-1 in PARP−/− mice was associated with a reduction in retinoblastoma(Rb) expression; increased expression of the oncogene Jun was correlated with a gain in 4C5-ter that harbors this oncogene. These results further implicate PARP in the maintenance of genomic stability and suggest that altered expression of p53, Rb, and Jun, as well as undoubtedly many other proteins may be a result of genomic instability associated with PARP deficiency.

Cloning and expression of rat 25-hydroxyvitamin D3-1α-hydroxylase cDNA

A full-length cDNA for the rat kidney mitochondrial cytochrome P450 mixed function oxidase, 25-hydroxyvitamin D3-1α-hydroxylase (P4501α), was cloned from a vitamin D-deficient rat kidney cDNA library and subcloned into the mammalian expression vector pcDNA 3.1(+). When P4501α cDNA was transfected into COS-7 transformed monkey kidney cells, they expressed 25-hydroxyvitamin D3-1α-hydroxylase activity. The sequence analysis showed that P4501α was of 2,469 bp long and contained an ORF encoding 501 amino acids. The deduced amino acid sequence showed a 53% similarity and 44% identity to the vitamin D3-25-hydroxylase (CYP27), whereas it has 42.6% similarity and 34% identity with the 25-hydroxyvitamin D3-24-hydroxylase (CYP24). Thus, it composes a new subfamily of the CYP27 family. Further, it is more closely related to the CYP27 than to the CYP24. The expression of P4501α mRNA was greatly increased in the kidney of vitamin D-deficient rats. In rats with the enhanced renal production of 1α,25-dihydroxyvitamin D3 (rats fed a low Ca diet), P4501α mRNA was greatly increased in the renal proximal convoluted tubules.

Cloning and expression of a cDNA encoding a bovine brain brefeldin A-sensitive guanine nucleotide-exchange protein for ADP-ribosylation factor

A 200-kDa guanine nucleotide-exchange protein (p200 or GEP) for ADP-ribosylation factors 1 and 3 (ARF1 and ARF3) that was inhibited by brefeldin A (BFA) was purified earlier from cytosol of bovine brain cortex. Amino acid sequences of four tryptic peptides were 47% identical to that of Sec7 from Saccharomyces cerevisiae, which is involved in vesicular trafficking in the Golgi. By using a PCR-based procedure with two degenerate primers representing sequences of these peptides, a product similar in size to Sec7 that contained the peptide sequences was generated. Two oligonucleotides based on this product were used to screen a bovine brain library, which yielded one clone that was a partial cDNA for p200. The remainder of the cDNA was obtained by 5′ and 3′ rapid amplification of cDNA ends (RACE). The ORF of the cDNA encodes a protein of 1,849 amino acids (≈208 kDa) that is 33% identical to yeast Sec7 and 50% identical in the Sec7 domain region. On Northern blot analysis of bovine tissues, a ≈7.4-kb mRNA was identified that hybridized with a p200 probe; it was abundant in kidney, somewhat less abundant in lung, spleen, and brain, and still less abundant in heart. A six-His-tagged fusion protein synthesized in baculovirus-infected Sf9 cells demonstrated BFA-inhibited GEP activity, confirming that BFA sensitivity is an intrinsic property of this ARF GEP and not conferred by another protein component of the complex from which p200 was originally purified.

Cluster analysis and display of genome-wide expression patterns

A system of cluster analysis for genome-wide expression data from DNA microarray hybridization is described that uses standard statistical algorithms to arrange genes according to similarity in pattern of gene expression. The output is displayed graphically, conveying the clustering and the underlying expression data simultaneously in a form intuitive for biologists. We have found in the budding yeast Saccharomyces cerevisiae that clustering gene expression data groups together efficiently genes of known similar function, and we find a similar tendency in human data. Thus patterns seen in genome-wide expression experiments can be interpreted as indications of the status of cellular processes. Also, coexpression of genes of known function with poorly characterized or novel genes may provide a simple means of gaining leads to the functions of many genes for which information is not available currently.

Cloning and functional expression of a cDNA encoding a pheromone gland-specific acyl-CoA Δ11-desaturase of the cabbage looper moth, Trichoplusia ni

Desaturation of coenzyme-A esters of saturated fatty acids is a common feature of sex pheromone biosynthetic pathways in the Lepidoptera. The enzymes that catalyze this step share several biochemical properties with the ubiquitous acyl-CoA Δ9-desaturases of animals and fungi, suggesting a common ancestral origin. Unlike metabolic acyl-CoA Δ9-desaturases, pheromone desaturases have evolved unusual regio- and stereoselective activities that contribute to the remarkable diversity of chemical structures used as pheromones in this large taxonomic group. In this report, we describe the isolation of a cDNA encoding a pheromone gland desaturase from the cabbage looper moth, Trichoplusia ni, a species in which all unsaturated pheromone products are produced via a Δ11Z-desaturation mechanism. The largest ORF of the ≈1,250-bp cDNA encodes a 349-aa apoprotein (PDesat-Tn Δ11Z) with a predicted molecular mass of 40,240 Da. Its hydrophobicity profile is similar overall to those of rat and yeast Δ9-desaturases, suggesting conserved transmembrane topology. A 182-aa core domain delimited by conserved histidine-rich motifs implicated in iron-binding and catalysis has 72 and 58% similarity (including conservative substitutions) to acyl-CoA Δ9Z-desaturases of rat and yeast, respectively. Northern blot analysis revealed an ≈1,250-nt PDesat-Tn Δ11Z mRNA that is consistent with the spatial and temporal distribution of Δ11-desaturase enzyme activity. Genetic transformation of a desaturase-deficient strain of the yeast Saccharomyces cerevisiae with an expression plasmid encoding PDesat-Tn Δ11Z resulted in complementation of the strain’s fatty acid auxotrophy and the production of Δ11Z-unsaturated fatty acids.