Characterization of Chelonus inanitus polydnavirus segments: sequences and analysis, excision site and demonstration of clustering

Polydnaviruses (genera Ichnovirus and Bracovirus) have a segmented genome of circular double-stranded DNA molecules, replicate in the ovary of parasitic wasps and are essential for successful parasitism of the host. Here we show the first detailed analysis of various segments of a bracovirus, the Chelonus inanitus virus (CiV). Four segments were sequenced and two of them, CiV12 and CiV14, were found to be closely related while CiV14.5 and CiV16.8 were unrelated. CiV12, CiV14.5 and CiV16.8 are unique while CiV14 occurs also nested in another larger segment. All four segments are predicted to contain genes and predictions could be substantiated in most cases. Comparison with databases revealed no significant similarities at either the nucleotide or amino acid level. Inverted repeats with identities between 77% and 92% and lengths between 26 bp and 100 bp were found on all segments outside of predicted genes. Hybridization experiments indicate that CiV12 and CiV14 are both flanked by other virus segments, suggesting that proviral CiV segments are clustered in the genome of the wasp. The integration/excision site of CiV14 was analysed and compared to that of CiV12. On both termini of proviral CiV12 and CiV14 as well as in the excised circular molecule and the rejoined DNA a very similar repeat of 14 bp was found. A model to illustrate where the terminal repeats might recombine to yield the circular molecule is presented. Excision of CiV12 and CiV14 is restricted to the female and sets in at a very specific time-point in pupal-adult development.

Conserved sequences and cell-specific DNase I hypersensitive sites upstream from the co-ordinately expressed alpha I- and alpha II-globin genes of Xenopus laevis

The globin gene family of Xenopus laevis comprises pairs of closely related genes that are arranged in two clusters, each pair of genes being co-ordinately and stage-specifically expressed. To get information on putative regulatory elements, we compared the DNA sequences and the chromatin conformation 5' to the co-ordinately expressed adult alpha-globin genes. Sequence analysis revealed a relatively conserved region from the cap site up to position -289, and further upstream seven distinct boxes of homology, separated by more diverged sequences or deletions/insertions. The homology boxes comprise 22 to 194 base-pairs showing 78 to 95% homology. Analysis of chromatin conformation showed that DNase I preferentially cuts the upstream region of both genes at similar positions, 5' to the T-A-T-A and the C-C-A-A-T boxes, only in chromatin of adult erythroblasts and erythrocytes, where adult globin genes are expressed, but not in chromatin of adult liver cells or larval erythrocytes, where these genes are silent. This suggests that cell- and stage-specific activation of these genes coincides with specific changes in chromatin conformation within the proximal upstream region. No difference was found in the nucleotide sequence within the DNase I hypersensitive region proximal to the adult alpha 1-globin gene in DNA from embryonic cells, in which this gene is inactive, and adult erythrocytes, expressing this gene.

DNA incisions stimulate genetic instability of triplet repeat sequences related to human hereditary neurological diseases

The molecular mechanisms responsible for the expansion and deletion of trinucleotide repeat sequences (TRS) are the focus of our studies. Several hereditary neurological diseases including Huntington's disease, myotonic dystrophy, and fragile X syndrome are associated with the instability of TRS. Using the well defined and controllable model system of Escherichia coli, the influences of three types of DNA incisions on genetic instability of CTG•CAG repeats were studied: DNA double-strand breaks (DSB), single-strand nicks, and single-strand gaps. The DNA incisions were generated in pUC19 derivatives by in vitro cleavage with restriction endonucleases. The cleaved DNA was then transformed into E. coli parental and mutant strains. Double-strand breaks induced deletions throughout the TRS region in an orientation dependent manner relative to the origin of replication. The extent of instability was enhanced by the repeat length and sequence (CTG•CAG vs. CGG•CCG). Mutations in recA and recBC increased deletions, mutations in recF stabilized the TRS, whereas mutations in ruvA had no effect. DSB were repaired by intramolecular recombination, versus an intermolecular gene conversion or crossover mechanism. 30 nt gaps formed a distinct 30 nt deletion product, whereas single strand nicks and gaps of 15 nts did not induce expansions or deletions. Formation of this deletion product required the CTG•CAG repeats to be present in the single-stranded region and was stimulated by E. coli DNA ligase, but was not dependent upon the RecFOR pathway. Models are presented to explain the DSB induced instabilities and formation of the 30 nucleotide deletion product. In addition to the in vitro creation of DSBs, several attempts to generate this incision in vivo with the use of EcoR I restriction modification systems were conducted. ^

Brefeldin A-inhibited guanine nucleotide-exchange activity of Sec7 domain from yeast Sec7 with yeast and mammalian ADP ribosylation factors

The Saccharomyces cerevisiae Sec7 protein (ySec7p), which is an important component of the yeast secretory pathway, contains a sequence of ≈200 amino acids referred to as a Sec7 domain. Similar Sec7 domain sequences have been recognized in several guanine nucleotide-exchange proteins (GEPs) for ADP ribosylation factors (ARFs). ARFs are ≈20-kDa GTPases that regulate intracellular vesicular membrane trafficking and activate phospholipase D. GEPs activate ARFs by catalyzing the replacement of bound GDP with GTP. We, therefore, undertook to determine whether a Sec7 domain itself could catalyze nucleotide exchange on ARF and found that it exhibited brefeldin A (BFA)-inhibitable ARF GEP activity. BFA is known to inhibit ARF GEP activity in Golgi membranes, thereby causing reversible apparent dissolution of the Golgi complex in many cells. The His6-tagged Sec7 domain from ySec7p (rySec7d) synthesized in Escherichia coli enhanced binding of guanosine 5′-[γ-[35S]thio]triphosphate by recombinant yeast ARF1 (ryARF1) and ryARF2 but not by ryARF3. The effects of rySec7d on ryARF2 were inhibited by BFA in a concentration-dependent manner but not by inactive analogues of BFA (B-17, B-27, and B-36). rySec7d also promoted BFA-sensitive guanosine 5′-[γ-thio]triphosphate binding by nonmyristoylated recombinant human ARF1 (rhARF1), rhARF5, and rhARF6, although the effect on rhARF6 was very small. These results are consistent with the conclusion that the yeast Sec7 domain itself contains the elements necessary for ARF GEP activity and its inhibition by BFA.

Isolation of a brefeldin A-inhibited guanine nucleotide-exchange protein for ADP ribosylation factor (ARF) 1 and ARF3 that contains a Sec7-like domain

Brefeldin A (BFA) inhibited the exchange of ADP ribosylation factor (ARF)-bound GDP for GTP by a Golgi-associated guanine nucleotide-exchange protein (GEP) [Helms, J. B. & Rothman, J. E. (1992) Nature (London) 360, 352–354; Donaldson, J. G., Finazzi, D. & Klausner, R. D. (1992) Nature (London) 360, 350–352]. Cytosolic ARF GEP was also inhibited by BFA, but after purification from bovine brain and rat spleen, it was no longer BFA-sensitive [Tsai, S.-C., Adamik, R., Moss, J. & Vaughan, M. (1996) Proc. Natl. Acad. Sci. USA 93, 305–309]. We describe here purification from bovine brain cytosol of a BFA-inhibited GEP. After chromatography on DEAE–Sephacel, hydroxylapatite, and Mono Q and precipitation at pH 5.8, GEP was eluted from Superose 6 as a large molecular weight complex at the position of thyroglobulin (≈670 kDa). After SDS/PAGE of samples from column fractions, silver-stained protein bands of ≈190 and 200 kDa correlated with activity. BFA-inhibited GEP activity of the 200-kDa protein was demonstrated following electroelution from the gel and renaturation by dialysis. Four tryptic peptides from the 200-kDa protein had amino acid sequences that were 47% identical to sequences in Sec7 from Saccharomyces cerevisiae (total of 51 amino acids), consistent with the view that the BFA-sensitive 200-kDa protein may be a mammalian counterpart of Sec7 that plays a similar role in cellular vesicular transport and Sec7 may be a GEP for one or more yeast ARFs.

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.

Interactive cloning with the SH3 domain of N-src identifies a new brain specific ion channel protein, with homology to Eag and cyclic nucleotide-gated channels

We have isolated a novel cDNA, that appears to represent a new class of ion channels, by using the yeast two-hybrid system and the SH3 domain of the neural form of Src (N-src) as a bait. The encoded polypeptide, BCNG-1, is distantly related to cyclic nucleotide-gated channels and the voltage-gated channels, Eag and H-erg. BCNG-1 is expressed exclusively in the brain, as a glycosylated protein of ≈132 kDa. Immunohistochemical analysis indicates that BCNG-1 is preferentially expressed in specific subsets of neurons in the neocortex, hippocampus, and cerebellum, in particular pyramidal neurons and basket cells. Within individual neurons, the BCNG-1 protein is localized to either the dendrites or the axon terminals depending on the cell type. Southern blot analysis shows that several other BCNG-related sequences are present in the mouse genome, indicating the emergence of an entire subfamily of ion channel coding genes. These findings suggest the existence of a new type of ion channel, which is potentially able to modulate membrane excitability in the brain and could respond to regulation by cyclic nucleotides.

Nucleotide sequence of the Kaposi sarcoma-associated herpesvirus (HHV8)

The genome of the Kaposi sarcoma-associated herpesvirus (KSHV or HHV8) was mapped with cosmid and phage genomic libraries from the BC-1 cell line. Its nucleotide sequence was determined except for a 3-kb region at the right end of the genome that was refractory to cloning. The BC-1 KSHV genome consists of a 140.5-kb-long unique coding region flanked by multiple G+C-rich 801-bp terminal repeat sequences. A genomic duplication that apparently arose in the parental tumor is present in this cell culture-derived strain. At least 81 ORFs, including 66 with homology to herpesvirus saimiri ORFs, and 5 internal repeat regions are present in the long unique region. The virus encodes homologs to complement-binding proteins, three cytokines (two macrophage inflammatory proteins and interleukin 6), dihydrofolate reductase, bcl-2, interferon regulatory factors, interleukin 8 receptor, neural cell adhesion molecule-like adhesin, and a D-type cyclin, as well as viral structural and metabolic proteins. Terminal repeat analysis of virus DNA from a KS lesion suggests a monoclonal expansion of KSHV in the KS tumor.

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.

The EMBL nucleotide sequence database

The EMBL Nucleotide Sequence Database (http://www.ebi.ac.uk/embl/) is maintained at the European Bioinformatics Institute (EBI) in an international collaboration with the DNA Data Bank of Japan (DDBJ) and GenBank at the NCBI (USA). Data is exchanged amongst the collaborating databases on a daily basis. The major contributors to the EMBL database are individual authors and genome project groups. Webin is the preferred web-based submission system for individual submitters, whilst automatic procedures allow incorporation of sequence data from large-scale genome sequencing centres and from the European Patent Office (EPO). Database releases are produced quarterly. Network services allow free access to the most up-to-date data collection via ftp, email and World Wide Web interfaces. EBI’s Sequence Retrieval System (SRS), a network browser for databanks in molecular biology, integrates and links the main nucleotide and protein databases plus many specialized databases. For sequence similarity searching a variety of tools (e.g. Blitz, Fasta, BLAST) are available which allow external users to compare their own sequences against the latest data in the EMBL Nucleotide Sequence Database and SWISS-PROT.

trEST, trGEN and Hits: access to databases of predicted protein sequences

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits.isb-sib.ch).

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Isolation of a superfamily of candidate disease-resistance genes in soybean based on a conserved nucleotide-binding site.

The tobacco N and Arabidopsis RPS2 genes, among several recently cloned disease-resistance genes, share highly conserved structure, a nucleotide-binding site (NBS). Using degenerate oligonucleotide primers for the NBS region of N and RPS2, we have amplified and cloned the NBS sequences from soybean. Each of these PCR-derived NBS clones detected low-or moderate-copy soybean DNA sequences and belongs to 1 of 11 different classes. Sequence analysis showed that all PCR clones encode three motifs (P-loop, kinase-2, and kinase-3a) of NBS nearly identical to those in N and RPS2. The intervening region between P-loop and kinase-3a of the 11 classes has high (26% average) amino acid sequence similarity to the N gene although not as high (19% average) to RPS2. These 11 classes represent a superfamily of NBS-containing soybean genes that are homologous to N and RPS2. Each class or subfamily was assessed for its positional association with known soybean disease-resistance genes through near-isogenic line assays, followed by linkage analysis in F2 populations using restriction fragment length polymorphisms. Five of the 11 subfamilies have thus far been mapped to the vicinity of known soybean genes for resistance to potyviruses (Rsv1 and Rpv), Phytophthora root rot (Rps1, Rps2, and Rps3), and powdery mildew (rmd). The conserved N- or RPS2-homologous NBS sequences and their positional associations with mapped soybean-resistance genes suggest that a number of the soybean disease-resistance genes may belong to this superfamily. The candidate subfamilies of NBS-containing genes identified by genetic mapping should greatly facilitate the molecular cloning of disease-resistance genes.

A general additive distance with time-reversibility and rate variation among nucleotide sites.

As additivity is a very useful property for a distance measure, a general additive distance is proposed under the stationary time-reversible (SR) model of nucleotide substitution or, more generally, under the stationary, time-reversible, and rate variable (SRV) model, which allows rate variation among nucleotide sites. A method for estimating the mean distance and the sampling variance is developed. In addition, a method is developed for estimating the variance-covariance matrix of distances, which is useful for the statistical test of phylogenies and molecular clocks. Computer simulation shows (i) if the sequences are longer than, say, 1000 bp, the SR method is preferable to simpler methods; (ii) the SR method is robust against deviations from time-reversibility; (iii) when the rate varies among sites, the SRV method is much better than the SR method because the distance is seriously underestimated by the SR method; and (iv) our method for estimating the sampling variance is accurate for sequences longer than 500 bp. Finally, a test is constructed for testing whether DNA evolution follows a general Markovian model.