Barnacle cement: a polymerization model based on evolutionary concepts.

Enzymes and biochemical mechanisms essential to survival are under extreme selective pressure and are highly conserved through evolutionary time. We applied this evolutionary concept to barnacle cement polymerization, a process critical to barnacle fitness that involves aggregation and cross-linking of proteins. The biochemical mechanisms of cement polymerization remain largely unknown. We hypothesized that this process is biochemically similar to blood clotting, a critical physiological response that is also based on aggregation and cross-linking of proteins. Like key elements of vertebrate and invertebrate blood clotting, barnacle cement polymerization was shown to involve proteolytic activation of enzymes and structural precursors, transglutaminase cross-linking and assembly of fibrous proteins. Proteolytic activation of structural proteins maximizes the potential for bonding interactions with other proteins and with the surface. Transglutaminase cross-linking reinforces cement integrity. Remarkably, epitopes and sequences homologous to bovine trypsin and human transglutaminase were identified in barnacle cement with tandem mass spectrometry and/or western blotting. Akin to blood clotting, the peptides generated during proteolytic activation functioned as signal molecules, linking a molecular level event (protein aggregation) to a behavioral response (barnacle larval settlement). Our results draw attention to a highly conserved protein polymerization mechanism and shed light on a long-standing biochemical puzzle. We suggest that barnacle cement polymerization is a specialized form of wound healing. The polymerization mechanism common between barnacle cement and blood may be a theme for many marine animal glues.

Ferrochelatase is a conserved downstream target of the blue light-sensing White collar complex in fungi.

Light is a universal signal perceived by organisms, including fungi, in which light regulates common and unique biological processes depending on the species. Previous research has established that conserved proteins, originally called White collar 1 and 2 from the ascomycete Neurospora crassa, regulate UV/blue light sensing. Homologous proteins function in distant relatives of N. crassa, including the basidiomycetes and zygomycetes, which diverged as long as a billion years ago. Here we conducted microarray experiments on the basidiomycete fungus Cryptococcus neoformans to identify light-regulated genes. Surprisingly, only a single gene was induced by light above the commonly used twofold threshold. This gene, HEM15, is predicted to encode a ferrochelatase that catalyses the final step in haem biosynthesis from highly photoreactive porphyrins. The C. neoformans gene complements a Saccharomyces cerevisiae hem15Delta strain and is essential for viability, and the Hem15 protein localizes to mitochondria, three lines of evidence that the gene encodes ferrochelatase. Regulation of HEM15 by light suggests a mechanism by which bwc1/bwc2 mutants are photosensitive and exhibit reduced virulence. We show that ferrochelatase is also light-regulated in a white collar-dependent fashion in N. crassa and the zygomycete Phycomyces blakesleeanus, indicating that ferrochelatase is an ancient target of photoregulation in the fungal kingdom.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

The receptor kinase family: primary structure of rhodopsin kinase reveals similarities to the beta-adrenergic receptor kinase.

Light-dependent deactivation of rhodopsin as well as homologous desensitization of beta-adrenergic receptors involves receptor phosphorylation that is mediated by the highly specific protein kinases rhodopsin kinase (RK) and beta-adrenergic receptor kinase (beta ARK), respectively. We report here the cloning of a complementary DNA for RK. The deduced amino acid sequence shows a high degree of homology to beta ARK. In a phylogenetic tree constructed by comparing the catalytic domains of several protein kinases, RK and beta ARK are located on a branch close to, but separate from the cyclic nucleotide-dependent protein kinase and protein kinase C subfamilies. From the common structural features we conclude that both RK and beta ARK are members of a newly delineated gene family of guanine nucleotide-binding protein (G protein)-coupled receptor kinases that may function in diverse pathways to regulate the function of such receptors.

Expansion of the alpha 2-adrenergic receptor family: cloning and characterization of a human alpha 2-adrenergic receptor subtype, the gene for which is located on chromosome 2.

Pharmacologic, biochemical, and genetic analyses have demonstrated the existence of multiple alpha 2-adrenergic receptor (alpha 2AR) subtypes. We have cloned a human alpha 2AR by using the polymerase chain reaction with oligonucleotide primers homologous to conserved regions of the previously cloned alpha 2ARs, the genes for which are located on human chromosomes 4 (C4) and 10 (C10). The deduced amino acid sequence encodes a protein of 450 amino acids whose putative topology is similar to that of the family of guanine nucleotide-binding protein-coupled receptors, but whose structure most closely resembles that of the alpha 2ARs. Competition curve analysis of the binding properties of the receptor expressed in COS-7 cells with a variety of adrenergic ligands demonstrates a unique alpha 2AR pharmacology. Hybridization with somatic cell hybrids shows that the gene for this receptor is located on chromosome 2. Northern blot analysis of various rat tissues shows expression in liver and kidney. The unique pharmacology and tissue localization of this receptor suggest that this is an alpha 2AR subtype not previously identified by classical pharmacological or ligand binding approaches.

Inhibition of beta-adrenergic receptor kinase prevents rapid homologous desensitization of beta 2-adrenergic receptors.

Homologous (agonist-specific) desensitization of beta-adrenergic receptors (beta ARs) is accompanied by and appears to require phosphorylation of the receptors. We have recently described a novel protein kinase, beta AR kinase, which phosphorylates beta ARs in vitro in an agonist-dependent manner. This kinase is inhibited by two classes of compounds, polyanions and synthetic peptides derived from the beta 2-adrenergic receptor (beta 2AR). In this report we describe the effects of these inhibitors on the process of homologous desensitization induced by the beta-adrenergic agonist isoproterenol. Permeabilization of human epidermoid carcinoma A431 cells with digitonin was used to permit access of the charged inhibitors to the cytosol; this procedure did not interfere with the pattern of isoproterenol-induced homologous desensitization of beta 2AR-stimulated adenylyl cyclase. Inhibitors of beta AR kinase markedly inhibited homologous desensitization of beta 2ARs in the permeabilized cells. Inhibition of desensitization by heparin, the most potent of the polyanion inhibitors of beta AR kinase, occurred over the same concentration range (5-50 nM) as inhibition of purified beta AR kinase assessed in a reconstituted system. Inhibition of desensitization by heparin was accompanied by a marked reduction of receptor phosphorylation in the permeabilized cells. Whereas inhibitors of beta AR kinase inhibited homologous desensitization, inhibitors of protein kinase C and of cyclic-nucleotide-dependent protein kinases were ineffective. These data establish that phosphorylation of beta ARs by beta AR kinase is an essential step in homologous desensitization of the receptors. They further suggest a potential therapeutic value of inhibitors of beta AR kinase in inhibiting agonist-induced desensitization.

Cloning of the cDNA for the human beta 1-adrenergic receptor.

Screening of a human placenta lambda gt11 library has led to the isolation of the cDNA for the human beta 1-adrenergic receptor (beta 1AR). Used as the probe was the human genomic clone termed G-21. This clone, which contains an intronless gene for a putative receptor, was previously isolated by virtue of its cross hybridization with the human beta 2-adrenergic receptor (beta 2AR). The 2.4-kilobase cDNA for the human beta 1AR encodes a protein of 477 amino acid residues that is 69% homologous with the avian beta AR but only 54% homologous with the human beta 2AR. This suggests that the avian gene encoding beta AR and the human gene encoding beta 1AR evolved from a common ancestral gene. RNA blot analysis indicates a message of 2.5 kilobases in rat tissues, with a pattern of tissue distribution consistent with beta 1AR binding. This pattern is quite distinct from the pattern obtained when the beta 2AR cDNA is used as a probe. Expression of receptor protein in Xenopus laevis oocytes conveys adenylate cyclase responsiveness to catecholamines with a typical beta 1AR specificity. This contrasts with the typical beta 2 subtype specificity observed when the human beta 2AR cDNA is expressed in this system. Mammalian beta 1AR and beta 2AR are thus products of distinct genes, both of which are apparently related to the putative G-21 receptor.

Interactions of chromatin context, binding site sequence content, and sequence evolution in stress-induced p53 occupancy and transactivation.

Cellular stresses activate the tumor suppressor p53 protein leading to selective binding to DNA response elements (REs) and gene transactivation from a large pool of potential p53 REs (p53REs). To elucidate how p53RE sequences and local chromatin context interact to affect p53 binding and gene transactivation, we mapped genome-wide binding localizations of p53 and H3K4me3 in untreated and doxorubicin (DXR)-treated human lymphoblastoid cells. We examined the relationships among p53 occupancy, gene expression, H3K4me3, chromatin accessibility (DNase 1 hypersensitivity, DHS), ENCODE chromatin states, p53RE sequence, and evolutionary conservation. We observed that the inducible expression of p53-regulated genes was associated with the steady-state chromatin status of the cell. Most highly inducible p53-regulated genes were suppressed at baseline and marked by repressive histone modifications or displayed CTCF binding. Comparison of p53RE sequences residing in different chromatin contexts demonstrated that weaker p53REs resided in open promoters, while stronger p53REs were located within enhancers and repressed chromatin. p53 occupancy was strongly correlated with similarity of the target DNA sequences to the p53RE consensus, but surprisingly, inversely correlated with pre-existing nucleosome accessibility (DHS) and evolutionary conservation at the p53RE. Occupancy by p53 of REs that overlapped transposable element (TE) repeats was significantly higher (p<10-7) and correlated with stronger p53RE sequences (p<10-110) relative to nonTE-associated p53REs, particularly for MLT1H, LTR10B, and Mer61 TEs. However, binding at these elements was generally not associated with transactivation of adjacent genes. Occupied p53REs located in L2-like TEs were unique in displaying highly negative PhyloP scores (predicted fast-evolving) and being associated with altered H3K4me3 and DHS levels. These results underscore the systematic interaction between chromatin status and p53RE context in the induced transactivation response. This p53 regulated response appears to have been tuned via evolutionary processes that may have led to repression and/or utilization of p53REs originating from primate-specific transposon elements.

Explicit DNase sequence bias modeling enables high-resolution transcription factor footprint detection.

DNaseI footprinting is an established assay for identifying transcription factor (TF)-DNA interactions with single base pair resolution. High-throughput DNase-seq assays have recently been used to detect in vivo DNase footprints across the genome. Multiple computational approaches have been developed to identify DNase-seq footprints as predictors of TF binding. However, recent studies have pointed to a substantial cleavage bias of DNase and its negative impact on predictive performance of footprinting. To assess the potential for using DNase-seq to identify individual binding sites, we performed DNase-seq on deproteinized genomic DNA and determined sequence cleavage bias. This allowed us to build bias corrected and TF-specific footprint models. The predictive performance of these models demonstrated that predicted footprints corresponded to high-confidence TF-DNA interactions. DNase-seq footprints were absent under a fraction of ChIP-seq peaks, which we show to be indicative of weaker binding, indirect TF-DNA interactions or possible ChIP artifacts. The modeling approach was also able to detect variation in the consensus motifs that TFs bind to. Finally, cell type specific footprints were detected within DNase hypersensitive sites that are present in multiple cell types, further supporting that footprints can identify changes in TF binding that are not detectable using other strategies.

Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.

Punctuated evolution and transitional hybrid network in an ancestral cell cycle of fungi.

Although cell cycle control is an ancient, conserved, and essential process, some core animal and fungal cell cycle regulators share no more sequence identity than non-homologous proteins. Here, we show that evolution along the fungal lineage was punctuated by the early acquisition and entrainment of the SBF transcription factor through horizontal gene transfer. Cell cycle evolution in the fungal ancestor then proceeded through a hybrid network containing both SBF and its ancestral animal counterpart E2F, which is still maintained in many basal fungi. We hypothesize that a virally-derived SBF may have initially hijacked cell cycle control by activating transcription via the cis-regulatory elements targeted by the ancestral cell cycle regulator E2F, much like extant viral oncogenes. Consistent with this hypothesis, we show that SBF can regulate promoters with E2F binding sites in budding yeast.