Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

Alignment of a 1.2-Mb chromosomal region from three strains of Rhodobacter capsulatus reveals a significantly mosaic structure.

High-resolution physical maps of the genomes of three Rhodobacter capsulatus strains, derived from ordered cosmid libraries, were aligned. The 1.2-Mb segment of the SB1003 genome studied here is adjacent to a 1-Mb region analyzed previously [Fonstein, M., Nikolskaya, T. & Haselkorn, H. (1995) J. Bacteriol. 177, 2368-2372]. Probes derived from the ordered cosmid set of R. capsulatus SB1003 were used to link cosmids from the St. Louis and 2.3.1 strain libraries. Cosmids selected this way did not merge into a single contig but formed several unlinked groups. EcoRV restriction maps of the ordered cosmids were then constructed using lambda terminase and fused to derive fragments of the chromosomal map. In order to link these fragments, their ends were transcribed to produce secondary probes for hybridization to gridded cosmid libraries of the same strains. This linking reduced the number of subcontigs to three for the St. Louis strain and one for the 2.3.1 strain. Hybridization of the same probes back to the ordered cosmid set of SB1003 positioned the subcontigs on the high-resolution physical map of SB1003. The final alignment of the restriction maps shows numerous large and small translocations in this 1.2-Mb chromosomal region of the three Rhodobacter strains. In addition, the chromosomes of the three strains, whose fine-structure maps can now be compared over 2.2 Mb, are seen to contain regions of 15-80 kb in which restriction sites are highly polymorphic, interspersed among regions in which the positions of restriction sites are highly conserved.

Molecular identification of a novel retrovirus repeatedly isolated from patients with multiple sclerosis

The partial molecular characterization of multiple sclerosis (MS)-associated retrovirus (MSRV), a novel retrovirus previously called LM7, is reported. MSRV has been isolated repeatedly from leptomeningeal, choroid plexus and from Epstein–Barr virus-immortalized B cells of MS patients. A strategy based on reverse transcriptase PCR with RNA-purified extracellular virions yielded an initial pol fragment from which other regions of the retroviral genome were subsequently obtained by sequence extension. MSRV-specific PCR primers amplified a pol region from RNA present at the peak of reverse transcriptase activity, coinciding with extracellular viral particles in sucrose density gradients. The same sequence was detected in noncellular RNA from MS patient plasma and in cerebrospinal fluid from untreated MS patients. MSRV is related to, but distinct from, the endogenous retroviral sequence ERV9. Whether MSRV represents an exogenous retrovirus with closely related endogenous elements or a replication-competent, virion-producing, endogenous provirus is as yet unknown. Further molecular epidemiological studies are required to determine precisely the apparent association of virions containing MSRV RNA with MS.

Identification of a phospholipase C beta2 region that interacts with Gbeta-gamma.

To delineate the phospholipase C (PLC; EC 3.1.4.3) beta2 sequences involved in interactions with the beta-gamma subunits of G proteins, we prepared a number of mammalian expression plasmids encoding a series of PLC beta2 segments that span the region from the beginning of the X box to the end of the Y box. We found the sequence extending from residue Glu-435 to residue Val-641 inhibited Gbeta-gamma-mediated activation of PLC beta2 in transfected COS-7 cells. This PLC beta2 sequence also inhibited ligand-induced activation of PLC in COS-7 cells cotransfected with cDNAs encoding the complement component C5a receptor and PLC beta2 but not in cells transfected with the alpha1B-adrenergic receptor, suggesting that the PLC beta2 residues (Glu-435 to Val-641) inhibit the Gbeta-gamma-mediated but not the Galpha-mediated effect. The inhibitory effect on Gbeta-gamma-mediated activation of PLC beta2 may be the result of the interaction between Gbeta-gamma and the PLC beta2 fragment. This idea was confirmed by the observation that a fusion protein comprising these residues (Glu-435 to Val-641) of PLC beta2 and glutathione S-transferase (GST) bound to Gbeta-gamma in an in vitro binding assay. The Gbeta-gamma-binding region was further narrowed down to 62 amino acids (residues Leu-580 to Val-641) by testing fusion proteins comprising various PLC beta2 sequences and GST in the in vitro binding assay.

Characterization of the promoter region of the mouse Xist gene.

The mouse Xist gene is expressed exclusively from the inactive X chromosome and may be implicated in initiating X inactivation. To better understand the mechanisms underlying the control of Xist expression, we investigated the upstream regulatory region of the mouse Xist promoter. A 1.2-kb upstream region of the Xist gene was sequenced and promoter activity was studied by chloramphenicol acetyltransferase (CAT) assays after transfection in murine XX and XY cell lines. The region analyzed (-1157 to +917 showed no in vitro sex-specific promoter activity. However, a minimal constitutional promoter was assigned to a region from -81 to +1, and a cis element from -41 to -15 regulates promoter activity. We showed that a nuclear factor binds to an element located at -30 to -25 (TTAAAG). A second sequence at -41 to -15 does not act as an enhancer and is unable to confer transcriptional activity to the Xist gene on its own. A third region from -82 to -41 is needed for correct expression. Deletion of the segment -441 to -231 is associated with an increase in CAT activity and may represent a silencer element.

Probing the Saccharomyces cerevisiae centromeric DNA (CEN DNA)–binding factor 3 (CBF3) kinetochore complex by using atomic force microscopy

Yeast centromeric DNA (CEN DNA) binding factor 3 (CBF3) is a multisubunit protein complex that binds to the essential CDEIII element in CEN DNA. The four CBF3 proteins are required for accurate chromosome segregation and are considered to be core components of the yeast kinetochore. We have examined the structure of the CBF3–CEN DNA complex by atomic force microscopy. Assembly of CBF3–CEN DNA complexes was performed by combining purified CBF3 proteins with a DNA fragment that includes the CEN region from yeast chromosome III. Atomic force microscopy images showed DNA molecules with attached globular bodies. The contour length of the DNA containing the complex is ≈9% shorter than the DNA alone, suggesting some winding of DNA within the complex. The measured location of the single binding site indicates that the complex is located asymmetrically to the right of CDEIII extending away from CDEI and CDEII, which is consistent with previous data. The CEN DNA is bent ≈55° at the site of complex formation. A significant fraction of the complexes are linked in pairs, showing three to four DNA arms, with molecular volumes approximately three times the mean volumes of two-armed complexes. These multi-armed complexes indicate that CBF3 can bind two DNA molecules together in vitro and, thus, may be involved in holding together chromatid pairs during mitosis.

DsrA RNA regulates translation of RpoS message by an anti-antisense mechanism, independent of its action as an antisilencer of transcription

DsrA RNA regulates both transcription, by overcoming transcriptional silencing by the nucleoid-associated H-NS protein, and translation, by promoting efficient translation of the stress σ factor, RpoS. These two activities of DsrA can be separated by mutation: the first of three stem-loops of the 85 nucleotide RNA is necessary for RpoS translation but not for anti-H-NS action, while the second stem-loop is essential for antisilencing and less critical for RpoS translation. The third stem-loop, which behaves as a transcription terminator, can be substituted by the trp transcription terminator without loss of either DsrA function. The sequence of the first stem-loop of DsrA is complementary with the upstream leader portion of rpoS messenger RNA, suggesting that pairing of DsrA with the rpoS message might be important for translational regulation. Mutations in the Rpos leader and compensating mutations in DsrA confirm that this predicted pairing is necessary for DsrA stimulation of RpoS translation. We propose that DsrA pairing stimulates RpoS translation by acting as an anti-antisense RNA, freeing the translation initiation region from the cis-acting antisense RNA and allowing increased translation.

RAG2 is regulated differentially in B and T cells by elements 5′ of the promoter

To study RAG2 gene regulation in vivo, we developed a blastocyst complementation method in which RAG2-deficient embryonic stem cells were transfected with genomic clones containing RAG2 and then assessed for their ability to generate lymphocytes. A RAG2 genomic clone that contained only the RAG2 promoter sequences rescued V(D)J recombination in RAG2-deficient pro-B cell lines, but did not rescue development of RAG2-deficient lymphocytes in vivo. However, inclusion of varying lengths of sequences 5′ of the RAG2 promoter generated constructs capable of rescuing only in vivo B cell development, as well as other constructs that rescued both B and T cell development. In particular, the 2-kb 5′ region starting just upstream of the RAG2 promoter, as well as the region from 2–7 kb 5′, could independently drive B cell development, but not efficient T cell development. Deletion of the 2-kb 5′ region from the murine germ line demonstrated that this region was not required for RAG expression sufficient to generate normal B or T cell numbers, implying redundancy among 5′ elements. We conclude that RAG2 expression in vivo requires elements beyond the core promoter, that such elements contribute to differential regulation in the B vs. T lineages, and that sequences sufficient to direct B cell expression are located in the promoter-proximal 5′ region.

Structural mimicry of a native protein by a minimized binding domain

The affinity between molecules depends both on the nature and presentation of the contacts. Here, we observe coupling of functional and structural elements when a protein binding domain is evolved to a smaller functional mimic. Previously, a 38-residue form of the 59-residue B-domain of protein A, termed Z38, was selected by phage display. Z38 contains 13 mutations and binds IgG only 10-fold weaker than the native B-domain. We present the solution structure of Z38 and show that it adopts a tertiary structure remarkably similar to that observed for the first two helices of B-domain in the B-domain/Fc complex [Deisenhofer, J. (1981) Biochemistry 20, 2361–2370], although it is significantly less stable. Based on this structure, we have improved on Z38 by designing a 34-residue disulfide-bonded variant (Z34C) that has dramatically enhanced stability and binds IgG with 9-fold higher affinity. The improved stability of Z34C led to NMR spectra with much greater chemical shift dispersion, resulting in a more precisely determined structure. Z34C, like Z38, has a structure virtually identical to the equivalent region from native protein A domains. The well-defined hydrophobic core of Z34C reveals key structural features that have evolved in this small, functional domain. Thus, the stabilized two-helix peptide, about half the size and having one-third of the remaining residues altered, accurately mimics both the structure and function of the native domain.

Sequence Determinants for Regulated Degradation of Yeast 3-Hydroxy-3-Methylglutaryl-CoA Reductase, an Integral Endoplasmic Reticulum Membrane Protein

The degradation rate of 3-hydroxy-3-methylglutaryl CoA reductase (HMG-R), a key enzyme of the mevalonate pathway, is regulated through a feedback mechanism by the mevalonate pathway. To discover the intrinsic determinants involved in the regulated degradation of the yeast HMG-R isozyme Hmg2p, we replaced small regions of the Hmg2p transmembrane domain with the corresponding regions from the other, stable yeast HMG-R isozyme Hmg1p. When the first 26 amino acids of Hmg2p were replaced with the same region from Hmg1p, Hmg2p was stabilized. The stability of this mutant was not due to mislocalization, but rather to an inability to be recognized for degradation. When amino acid residues 27–54 of Hmg2p were replaced with those from Hmg1p, the mutant was still degraded, but its degradation rate was poorly regulated. The degradation of this mutant was still dependent on the first 26 amino acid residues and on the function of the HRD genes. These mutants showed altered ubiquitination levels that were well correlated with their degradative phenotypes. Neither determinant was sufficient to impart regulated degradation to Hmg1p. These studies provide evidence that there are sequence determinants in Hmg2p necessary for degradation and optimal regulation, and that independent processes may be involved in Hmg2p degradation and its regulation.

The photoisomerization of retinal in bacteriorhodopsin: Experimental evidence for a three-state model

The primary events in the all-trans to 13-cis photoisomerization of retinal in bacteriorhodopsin have been investigated with femtosecond time-resolved absorbance spectroscopy. Spectra measured over a broad range extending from 7000 to 22,400 cm−1 reveal features whose dynamics are inconsistent with a model proposed earlier to account for the highly efficient photoisomerization process. Emerging from this work is a new three-state model. Photoexcitation of retinal with visible light accesses a shallow well on the excited state potential energy surface. This well is bounded by a small barrier, arising from an avoided crossing that separates the Franck–Condon region from the nearby reactive region of the photoisomerization coordinate. At ambient temperatures, the reactive region is accessed with a time constant of ≈500 fs, whereupon the retinal rapidly twists and encounters a second avoided crossing region. The protein mediates the passage into the second avoided crossing region and thereby exerts control over the quantum yield for forming 13-cis retinal. The driving force for photoisomerization resides in the retinal, not in the surrounding protein. This view contrasts with an earlier model where photoexcitation was thought to access directly a reactive region of the excited-state potential and thereby drive the retinal to a twisted conformation within 100–200 fs.

A functional genetic screen identifies regions at the C-terminal tail and death-domain of death-associated protein kinase that are critical for its proapoptotic activity

Death-associated protein kinase (DAP-kinase) is a Ca+2/calmodulin-regulated serine/threonine kinase with a multidomain structure that participates in apoptosis induced by a variety of signals. To identify regions in this protein that are critical for its proapoptotic activity, we performed a genetic screen on the basis of functional selection of short DAP-kinase-derived fragments that could protect cells from apoptosis by acting in a dominant-negative manner. We expressed a library of randomly fragmented DAP-kinase cDNA in HeLa cells and treated these cells with IFN-γ to induce apoptosis. Functional cDNA fragments were recovered from cells that survived the selection, and those in the sense orientation were examined further in a secondary screen for their ability to protect cells from DAP-kinase-dependent tumor necrosis factor-α-induced apoptosis. We isolated four biologically active peptides that mapped to the ankyrin repeats, the “linker” region, the death domain, and the C-terminal tail of DAP-kinase. Molecular modeling of the complete death domain provided a structural basis for the function of the death-domain-derived fragment by suggesting that the protective fragment constitutes a distinct substructure. The last fragment, spanning the C-terminal serine-rich tail, defined a new regulatory region. Ectopic expression of the tail peptide (17 amino acids) inhibited the function of DAP-kinase, whereas removal of this region from the complete protein caused enhancement of the killing activity, indicating that the C-terminal tail normally plays a negative regulatory role. Altogether, this unbiased screen highlighted functionally important regions in the protein and revealed an additional level of regulation of DAP-kinase apoptotic function that does not affect the catalytic activity.

A human protein containing multiple types of protease-inhibitory modules

By using sensitive homology-search and gene-finding programs, we have found that a genomic region from the tip of the short arm of human chromosome 16 (16p13.3) encodes a putative secreted protein consisting of a domain related to the whey acidic protein (WAP) domain, a domain homologous with follistatin modules of the Kazal-domain family (FS module), an immunoglobulin-related domain (Ig domain), two tandem domains related to Kunitz-type protease inhibitor modules (KU domains), and a domain belonging to the recently defined NTR-module family (NTR domain). The gene encoding these WAP, FS, Ig, KU, and NTR modules (hereafter referred to as the WFIKKN gene) is intron-depleted—its single 1,157-bp intron splits the WAP module. The validity of our gene prediction was confirmed by sequencing a WFIKKN cDNA cloned from a lung cDNA library. Studies on the tissue-expression pattern of the WFIKKN gene have shown that the gene is expressed primarily in pancreas, kidney, liver, placenta, and lung. As to the function of the WFIKKN protein, it is noteworthy that it contains FS, WAP, and KU modules, i.e., three different module types homologous with domains frequently involved in inhibition of serine proteases. The protein also contains an NTR module, a domain type implicated in inhibition of zinc metalloproteinases of the metzincin family. On the basis of its intriguing homologies, we suggest that the WFIKKN protein is a multivalent protease inhibitor that may control the action of multiple types of serine proteases as well as metalloproteinase(s).

Multiple maternal origins and weak phylogeographic structure in domestic goats

Domestic animals have played a key role in human history. Despite their importance, however, the origins of most domestic species remain poorly understood. We assessed the phylogenetic history and population structure of domestic goats by sequencing a hypervariable segment (481 bp) of the mtDNA control region from 406 goats representing 88 breeds distributed across the Old World. Phylogeographic analysis revealed three highly divergent goat lineages (estimated divergence >200,000 years ago), with one lineage occurring only in eastern and southern Asia. A remarkably similar pattern exists in cattle, sheep, and pigs. These results, combined with recent archaeological findings, suggest that goats and other farm animals have multiple maternal origins with a possible center of origin in Asia, as well as in the Fertile Crescent. The pattern of goat mtDNA diversity suggests that all three lineages have undergone population expansions, but that the expansion was relatively recent for two of the lineages (including the Asian lineage). Goat populations are surprisingly less genetically structured than cattle populations. In goats only ≈10% of the mtDNA variation is partitioned among continents. In cattle the amount is ≥50%. This weak structuring suggests extensive intercontinental transportation of goats and has intriguing implications about the importance of goats in historical human migrations and commerce.

Analysis of Promoter Activity for the Gene Encoding Pyruvate Orthophosphate Dikinase in Stably Transformed C4 Flaveria Species1

The C4 enzyme pyruvate orthophosphate dikinase is encoded by a single gene, Pdk, in the C4 plant Flaveria trinervia. This gene also encodes enzyme isoforms located in the chloroplast and in the cytosol that do not have a function in C4 photosynthesis. Our goal is to identify cis-acting DNA sequences that regulate the expression of the gene that is active in the C4 cycle. We fused 1.5 kb of a 5′ flanking region from the Pdk gene, including the entire 5′ untranslated region, to the uidA reporter gene and stably transformed the closely related C4 species Flaveria bidentis. β-Glucuronidase (GUS) activity was detected at high levels in leaf mesophyll cells. GUS activity was detected at lower levels in bundle-sheath cells and stems and at very low levels in roots. This lower-level GUS expression was similar to the distribution of mRNA encoding the nonphotosynthetic form of the enzyme. We conclude that cis-acting DNA sequences controlling the expression of the C4 form in mesophyll cells and the chloroplast form in other cells and organs are co-located within the same 5′ region of the Pdk gene.