Restriction-modification gene complexes as selfish gene entities: Roles of a regulatory system in their establishment, maintenance, and apoptotic mutual exclusion

We have reported some type II restriction-modification (RM) gene complexes on plasmids resist displacement by an incompatible plasmid through postsegregational host killing. Such selfish behavior may have contributed to the spread and maintenance of RM systems. Here we analyze the role of regulatory genes (C), often found linked to RM gene complexes, in their interaction with the host and the other RM gene complexes. We identified the C gene of EcoRV as a positive regulator of restriction. A C mutation eliminated postsegregational killing by EcoRV. The C system has been proposed to allow establishment of RM systems in new hosts by delaying the appearance of restriction activity. Consistent with this proposal, bacteria preexpressing ecoRVC were transformed at a reduced efficiency by plasmids carrying the EcoRV RM gene complex. Cells carrying the BamHI RM gene complex were transformed at a reduced efficiency by a plasmid carrying a PvuII RM gene complex, which shares the same C specificity. The reduction most likely was caused by chromosome cleavage at unmodified PvuII sites by prematurely expressed PvuII restriction enzyme. Therefore, association of the C genes of the same specificity with RM gene complexes of different sequence specificities can confer on a resident RM gene complex the capacity to abort establishment of a second, incoming RM gene complex. This phenomenon, termed “apoptotic mutual exclusion,” is reminiscent of suicidal defense against virus infection programmed by other selfish elements. pvuIIC and bamHIC genes define one incompatibility group of exclusion whereas ecoRVC gene defines another.

Deletion of long-range regulatory elements upstream of SOX9 causes campomelic dysplasia

Campomelic dysplasia (CD) is a rare, neonatal human chondrodysplasia characterized by bowing of the long bones and often associated with male-to-female sex-reversal. Patients present with either heterozygous mutations in the SOX9 gene or chromosome rearrangements mapping at least 50 kb upstream of SOX9. Whereas mutations in SOX9 ORF cause haploinsufficiency, the effects of translocations 5′ to SOX9 are unclear. To test whether these rearrangements also cause haploinsufficiency by altering spatial and temporal expression of SOX9, we generated mice transgenic for human SOX9-lacZ yeast artificial chromosomes containing variable amounts of DNA sequences upstream of SOX9. We show that elements necessary for SOX9 expression during skeletal development are highly conserved between mouse and human and reveal that a rearrangement upstream of SOX9, similar to those observed in CD patients, leads to a substantial reduction of SOX9 expression, particularly in chondrogenic tissues. These data demonstrate that important regulatory elements are scattered over a large region upstream of SOX9 and explain how particular aspects of the CD phenotype are caused by chromosomal rearrangements 5′ to SOX9.

Nonsynonymous substitution in abalone sperm fertilization genes exceeds substitution in introns and mitochondrial DNA

Strong positive Darwinian selection acts on two sperm fertilization proteins, lysin and 18-kDa protein, from abalone (Haliotis). To understand the phylogenetic context for this dramatic molecular evolution, we obtained sequences of mitochondrial cytochrome c oxidase subunit I (mtCOI), and genomic sequences of lysin, 18-kDa, and a G protein subunit. Based on mtDNA differentiation, four north Pacific abalone species diverged within the past 2 million years (Myr), and remaining north Pacific species diverged over a period of 4–20 Myr. Between-species nonsynonymous differences in lysin and 18-kDa exons exceed nucleotide differences in introns by 3.5- to 24-fold. Remarkably, in some comparisons nonsynonymous substitutions in lysin and 18-kDa genes exceed synonymous substitutions in mtCOI. Lysin and 18-kDa intron/exon segments were sequenced from multiple red abalone individuals collected over a 1,200-km range. Only two nucleotide changes and two sites of slippage variation were detected in a total of >29,000 nucleotides surveyed. However, polymorphism in mtCOI and a G protein intron was found in this species. This finding suggests that positive selection swept one lysin allele and one 18-kDa allele to fixation. Similarities between mtCOI and lysin gene trees indicate that rapid adaptive evolution of lysin has occurred consistently through the history of the group. Comparisons with mtCOI molecular clock calibrations suggest that nonsynonymous substitutions accumulate 2–50 times faster in lysin and 18-kDa genes than in rapidly evolving mammalian genes.

Oocytes are a source of catecholamines in the primate ovary: Evidence for a cell–cell regulatory loop

Catecholamines, thought to derive from the extrinsic innervation of the ovary, participate in the regulation of ovarian development and mature gonadal function. Recently, intraovarian neurons containing tyrosine hydroxylase (TH), the rate-limiting enzyme in catecholamine biosynthesis, were described in the ovary of nonhuman primates. We now show that the primate ovary expresses both the genes encoding TH and dopamine β-hydroxylase (DBH), the key enzymes in norepinephrine (NE) biosynthesis. Ovarian neurons were identified as a site of TH and DBH gene expression, and surprisingly, oocytes were identified as an exclusive site of DBH synthesis. Oocytes contain neither TH mRNA nor protein, indicating that they are unable to synthesize dopamine (DA). They did, however, express a DA transporter gene identical to that found in human brain. The physiological relevance of this transporter system and DBH in oocytes was indicated by the ability of isolated oocytes to metabolize exogenous DA into NE. Isolated follicles containing oocytes—but not those from which the oocytes had been removed—responded to DA with an elevation in cAMP levels; this elevation was prevented by propranolol, a β-adrenoreceptor antagonist. The results suggest that oocytes and somatic cells are linked by a neuroendocrine loop consisting of NE synthesized in oocytes from actively transported DA and cAMP produced by somatic follicular cells in response to NE-induced β-adrenoreceptor activation.

Comparative studies on mammalian Hoxc8 early enhancer sequence reveal a baleen whale-specific deletion of a cis-acting element

Variations in regulatory regions of developmental control genes have been implicated in the divergence of axial morphologies. To find potentially significant changes in cis-regulatory regions, we compared nucleotide sequences and activities of mammalian Hoxc8 early enhancers. The nucleotide sequence of the early enhancer region is extremely conserved among mammalian clades, with five previously described cis-acting elements, A–E, being invariant. However, a 4-bp deletion within element C of the Hoxc8 early enhancer sequence is observed in baleen whales. When assayed in transgenic mouse embryos, a baleen whale enhancer (unlike other mammalian enhancers) directs expression of the reporter gene to more posterior regions of the neural tube but fails to direct expression to posterior mesoderm. We suggest that regulation of Hoxc8 in baleen whales differs from other mammalian species and may be associated with variation in axial morphology.

Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix–loop–helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

Backtracking leukemia to birth: Identification of clonotypic gene fusion sequences in neonatal blood spots

Epidemiological evidence has suggested that some pediatric leukemias may be initiated in utero and, for some pairs of identical twins with concordant leukemia, this possibility has been strongly endorsed by molecular studies of clonality. Direct evidence for a prenatal origin can only be derived by prospective or retrospective detection of leukemia-specific molecular abnormalities in fetal or newborn samples. We report a PCR-based method that has been developed to scrutinize neonatal blood spots (Guthrie cards) for the presence of numerically infrequent leukemic cells at birth in individuals who subsequently developed leukemia. We demonstrate that unique or clonotypic MLL-AF4 genomic fusion sequences are present and detectable in neonatal blood spots from individuals who were diagnosed with acute lymphoblastic leukemia at ages 5 months to 2 years and, therefore, have arisen during fetal hematopoiesis in utero. This result provides unequivocal evidence for a prenatal initiation of acute leukemia in young patients. The method should be applicable to other fusion genes in children with common subtypes of leukemia and will be of value in attempts to unravel the natural history and etiology of this major subtype of pediatric cancer.

The Ath5 proneural genes function upstream of Brn3 POU domain transcription factor genes to promote retinal ganglion cell development

During retinogenesis, the Xenopus basic helix–loop–helix transcription factor Xath5 has been shown to promote a ganglion cell fate. In the developing mouse and chicken retinas, gene targeting and overexpression studies have demonstrated critical roles for the Brn3 POU domain transcription factor genes in the promotion of ganglion cell differentiation. However, the genetic relationship between Ath5 and Brn3 genes is unknown. To understand the genetic regulatory network(s) that controls retinal ganglion cell development, we analyzed the relationship between Ath5 and Brn3 genes by using a gain-of-function approach in the chicken embryo. We found that during retinogenesis, the chicken Ath5 gene (Cath5) is expressed in retinal progenitors and in differentiating ganglion cells but is absent in terminally differentiated ganglion cells. Forced expression of both Cath5 and the mouse Ath5 gene (Math5) in retinal progenitors activates the expression of cBrn3c following central-to-peripheral and temporal-to-nasal gradients. As a result, similar to the Xath5 protein, both Cath5 and Math5 proteins have the ability to promote the development of ganglion cells. Moreover, we found that forced expression of all three Brn3 genes also can stimulate the expression of cBrn3c. We further found that Ath5 and Brn3 proteins are capable of transactivating a Brn3b promoter. Thus, these data suggest that the expression of cBrn3c in the chicken and Brn3b in the mouse is initially activated by Ath5 factors in newly generated ganglion cells and later maintained by a feedback loop of Brn3 factors in the differentiated ganglion cells.

Lamprey Dlx genes and early vertebrate evolution

Gnathostome vertebrates have multiple members of the Dlx family of transcription factors that are expressed during the development of several tissues considered to be vertebrate synapomorphies, including the forebrain, cranial neural crest, placodes, and pharyngeal arches. The Dlx gene family thus presents an ideal system in which to examine the relationship between gene duplication and morphological innovation during vertebrate evolution. Toward this end, we have cloned Dlx genes from the lamprey Petromyzon marinus, an agnathan vertebrate that occupies a critical phylogenetic position between cephalochordates and gnathostomes. We have identified four Dlx genes in P. marinus, whose orthology with gnathostome Dlx genes provides a model for how this gene family evolved in the vertebrate lineage. Differential expression of these lamprey Dlx genes in the forebrain, cranial neural crest, pharyngeal arches, and sensory placodes of lamprey embryos provides insight into the developmental evolution of these structures as well as a model of regulatory evolution after Dlx gene duplication events.

A single gene for lycopene cyclase, phytoene synthase, and regulation of carotene biosynthesis in Phycomyces

Previous complementation and mapping of mutations that change the usual yellow color of the Zygomycete Phycomyces blakesleeanus to white or red led to the definition of two structural genes for carotene biosynthesis. We have cloned one of these genes, carRA, by taking advantage of its close linkage to the other, carB, responsible for phytoene dehydrogenase. The sequences of the wild type and six mutants have been established, compared with sequences in other organisms, and correlated with the mutant phenotypes. The carRA and carB coding sequences are separated by 1,381 untranslated nucleotides and are divergently transcribed. Gene carRA contains separate domains for two enzymes, lycopene cyclase and phytoene synthase, and regulates the overall activity of the pathway and its response to physical and chemical stimuli from the environment. The lycopene cyclase domain of carRA derived from a duplication of a gene from a common ancestor of fungi and Brevibacterium linens; the phytoene synthase domain is similar to the phytoene and squalene synthases of many organisms; but the regulatory functions appear to be specific to Phycomyces.

Herpes simplex virus 1 α regulatory protein ICP0 functionally interacts with cellular transcription factor BMAL1

The infected cell protein no. 0 (ICP0) of herpes simplex virus 1 (HSV-1) is a promiscuous transactivator shown to enhance the expression of gene introduced into cells by infection or transfection. At the molecular level, ICP0 is a 775-aa ring finger protein localized initially in the nucleus and late in infection in the cytoplasm and mediates the degradation of several proteins and stabilization of others. None of the known functions at the molecular level account for the apparent activity of ICP0 as a transactivator. Here we report that ICP0 functionally interacts with cellular transcription factor BMAL1, a member of the basic helix–loop–helix PER-ARNT-SIM (PAS) super family of transcriptional regulators. Specifically, sequences mapped to the exon II of ICP0 interacted with BMAL1 in the yeast two-hybrid system and in reciprocal pull-down experiments in vitro. Moreover, the enhancement of transcription of a luciferase reporter construct whose promoter contained multiple BMAL1-binding sites by ICP0 and BMAL1 was significantly greater than that observed by ICP0 or BMAL1 alone. Although the level of BMAL1 present in nuclei of infected cells remained unchanged between 3 and 8 h after infection, the level of cytoplasmic BMAL1 was reduced at 8 h after infection. The reduction of cytoplasmic BMAL1 was significantly greater in cells infected with the ICP0-null mutant than in the wild-type virus-infected cells, suggesting that ICP0 mediates partial stabilization of the protein. These results indicate that ICP0 interacts physically and functionally with at least one cellular transcription-regulatory factor.

Arrayed transposase-binding sequences on the ends of transposon Tn5090/Tn402

The transposon Tn5090/Tn402 encodes a 559 amino acid transposase, TniA, with a DDE motif. Gel mobility shifting and cleavage protection analysis with DNase I and hydroxyl radical probes revealed that TniA binds to multiple repeat sequences on either terminus of Tn5090/Tn402. Four of these TniA-binding 19mers occurred on the left-hand (t) end and two on the right-hand (i) end. Hydroxyl radical cleavage protection demonstrated the presence of 3–6 bp contact sequences on one face of the DNA helix. The binding pattern and organisation of repeats suggested parallels between Tn5090/Tn402 and Mu, which controls its transpositional activity in the assembly step of a higher order transpososome complex. The complex terminal structure and genes of transposase and nucleotide-binding proteins in tandem are hallmarks of the handful of Mu-like elements that are known to date.

Identification of sample-specific sequences in mammalian cDNA and genomic DNA by the novel ligation-mediated subtraction (Limes)

The representational difference analysis (RDA) and other subtraction techniques are used to enrich sample-specific sequences by elimination of ubiquitous sequences existing in both the sample of interest (tester) and the subtraction partner (driver). While applying the RDA to genomic DNA of cutaneous lymphoma cells in order to identify tumor relevant alterations, we predominantly isolated repetitive sequences and artificial repeat-mediated fusion products of otherwise independent PCR fragments (PCR hybrids). Since these products severely interfered with the isolation of tester-specific fragments, we developed a considerably more robust and efficient approach, termed ligation-mediated subtraction (Limes). In first applications of Limes, genomic sequences and/or transcripts of genes involved in the regulation of transcription, such as transforming growth factor β stimulated clone 22 related gene (TSC-22R), cell death and cytokine production (caspase-1) or antigen presentation (HLA class II sequences), were found to be completely absent in a cutaneous lymphoma line. On the assumption that mutations in tumor-relevant genes can affect their transcription pattern, a protocol was developed and successfully applied that allows the identification of such sequences. Due to these results, Limes may substitute/supplement other subtraction/comparison techniques such as RDA or DNA microarray techniques in a variety of different research fields.

Conserved regulatory elements of the promoter sequence of the gene rpoH of enteric bacteria

The rpoH regulatory region of different members of the enteric bacteria family was sequenced or downloaded from GenBank and compared. In addition, the transcriptional start sites of rpoH of Yersinia frederiksenii and Proteus mirabilis, two distant members of this family, were determined. Sequences similar to the σ70 promoters P1, P4 and P5, to the σE promoter P3 and to boxes DnaA1, DnaA2, cAMP receptor protein (CRP) boxes CRP1, CRP2 and box CytR present in Escherichia coli K12, were identified in sequences of closely related bacteria such as: E.coli, Shigella flexneri, Salmonella enterica serovar Typhimurium, Citrobacter freundii, Enterobacter cloacae and Klebsiella pneumoniae. In more distant bacteria, Y.frederiksenii and P.mirabilis, the rpoH regulatory region has a distal P1-like σ70 promoter and two proximal promoters: a heat-induced σE-like promoter and a σ70 promoter. Sequences similar to the regulatory boxes were not identified in these bacteria. This study suggests that the general pattern of transcription of the rpoH gene in enteric bacteria includes a distal σ70 promoter, >200 nt upstream of the initiation codon, and two proximal promoters: a heat-induced σE-like promoter and a σ70 promoter. A second proximal σ70 promoter under catabolite-regulation is probably present only in bacteria closely related to E.coli.

Isolation of segments of homologous genes with only one conserved amino acid region via PCR

We present a method which allows the isolation of fragments from genes coding for homologous proteins via PCR when only one block of conserved amino acids is available. Sets of degenerated primers are defined by reverse translation of the conserved amino acids such that each set contains not more than 128 different sequences. The second primer binding site is provided by a special cassette that is designed such that it does not allow binding of the second primer prior to being copied by DNA synthesis. The cassette is ligated to partially-digested chromosomal DNA. The second primer is biotinylated to allow elimination of PCR products carrying degenerated primers on both sides via streptavidin binding. Fragments obtained after amplification and enrichment are cloned and sequenced. The feasibility of this method was demonstrated in a model experiment, where degenerated primers were deduced from six conserved amino acids within the family of homologs to the Escherichia coli Vsr protein.