Structure and evolution of a proviral locus of Glyptapanteles indiensis bracovirus

Background

Bracoviruses (BVs), a group of double-stranded DNA viruses with segmented genomes, are mutualistic endosymbionts of parasitoid wasps. Virus particles are replication deficient and are produced only by female wasps from proviral sequences integrated into the wasp genome. Virus particles are injected along with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and therefore perpetuation of proviral DNA. Here we describe a 223 kbp region of Glyptapanteles indiensis genomic DNA which contains a part of the G. indiensis bracovirus (GiBV) proviral genome.

Results

Eighteen of ~24 GiBV viral segment sequences are encoded by 7 non-overlapping sets of BAC clones, revealing that some proviral segment sequences are separated by long stretches of intervening DNA. Two overlapping BACs, which contain a locus of 8 tandemly arrayed proviral segments flanked on either side by ~35 kbp of non-packaged DNA, were sequenced and annotated. Structural and compositional analyses of this cluster revealed it exhibits a G+C and nucleotide composition distinct from the flanking DNA. By analyzing sequence polymorphisms in the 8 GiBV viral segment sequences, we found evidence for widespread selection acting on both protein-coding and non-coding DNA. Comparative analysis of viral and proviral segment sequences revealed a sequence motif involved in the excision of proviral genome segments which is highly conserved in two other bracoviruses.

Conclusion

Contrary to current concepts of bracovirus proviral genome organization our results demonstrate that some but not all GiBV proviral segment sequences exist in a tandem array. Unexpectedly, non-coding DNA in the 8 proviral genome segments which typically occupies ~70% of BV viral genomes is under selection pressure suggesting it serves some function(s). We hypothesize that selection acting on GiBV proviral sequences maintains the genetic island-like nature of the cluster of proviral genome segments described herein. In contrast to large differences in the predicted gene composition of BV genomes, sequences that appear to mediate processes of viral segment formation, such as proviral segment excision and circularization, appear to be highly conserved, supporting the hypothesis of a single origin for BVs.     

Homologous sequences in the Campoletis sonorensis polydnavirus genome are implicated in replication and nesting of the W segment family.

Polydnaviruses (PDVs) are double-stranded DNA viruses with segmented genomes that replicate only in the oviducts of some species of parasitic wasps and are required for the successful parasitization of lepidopteran insects. PDV DNA segments are integrated in the genomes of their associated wasp hosts, and some are nested; i.e., smaller segments are produced from and largely colinear with larger segments. To determine the internal structure of nested viral segments, the first complete nucleotide sequence of a PDV genome segment and its integration locus was determined. By restriction mapping, Southern blot, and sequence analyses, we demonstrated that the Campoletis sonorensis PDV segment W is integrated into wasp genomic DNA. DNA sequence analysis revealed that proviral segment W terminates in two 1,185-bp direct long terminal repeats (LTRs) in the wasp chromosome, while only one LTR copy is present in the extrachromosomal (viral) W. The results suggest that terminal direct repeats are a general feature of PDV DNA segment integration but that the homology and size of the repeats can vary extensively. Segment W contains 12 imperfect direct repeats of six different types between 89 bp and 1.9 kbp with 65 to 90% homology. The orientation and structure of the repeats suggest that W itself may have arisen through sequence duplication and subsequent divergence. Mapping, hybridization, and sequence analyses of cloned R and M demonstrated that these segments are nested within segment W and that internal imperfect direct repeats of one type are implicated in the homologous intramolecular recombination events that generate segments R and M. Interestingly, segment nesting differentially increases the copy number of genes encoded by segment W, suggesting that the unusual genomic organization of PDVs may be directly linked to the unique functions of this virus in its obligate mutualistic association with parasitic wasps.     

A programmed site-specific DNA rearrangement in Tetrahymena thermophila requires flanking polypurine tracts

During macronuclear development in ciliates, precise deletion events eliminate thousands of specific DNA segments. Each segment is bounded by a unique pair of short direct repeats, but no other common feature has been reported. To determine the critical cis-acting sequences, we developed an in vivo system for analyzing this process in Tetrahymena. We show that sequences essential for recognition and excision of one such region are located within the 70 bp of DNA flanking either side of it. Three authentic splice sites and one cryptic site are each adjacent to an unusual polypurine tract (5'-A5G5) situated 40-50 bp distal to each terminal repeat. Removal of this tract or substitution of 3 bp within it abolishes splicing to the adjacent site. The normal chromosomal environment and the integrity of the eliminated sequence are not required for its removal. We believe the polypurine tract is a signal essential for excision of this sequence.     

Characterization of an unusual DNA length polymorphism 5'' to the human antithrombin III gene

Nucleotide sequence analysis revealed that a DNA length polymorphism 5' to the human antithrombin III gene is due to the presence of 32bp or 108bp nonhomologous nucleotide sequences (variable segments) 345bp upstream from the translation initiation codon. Sequences at the 3' borders of both variable segments can form intrastrand inverted repeat structures with sequences further downstream. An inverted repeat is also found immediately 5' to the site where the variable segments are located. Thus, cruciform structures may form flanking the variable segments of both alleles of this DNA length polymorphism. DNA secondary structure may be detected with single strand specific nucleases. S1 nuclease sensitive sites were mapped in recombinant plasmids containing the cloned alleles of the ATIII length polymorphism. The site most sensitive to S1 is located upstream from the variable segments in an AT-rich segment flanked by 6bp direct repeats. A region of lesser nuclease sensitivity was also observed in the AT-rich loops formed between the inverted repeats 5' to the variable segments.     

DNA sequence of the murine gamma 1 switch segment reveals novel structural elements

Immunoglobulin heavy chain switch regions are segments of DNA considered to be important in mediating class switching in B lymphocytes. Whereas these segments vary in length among the different murine isotypes, their structural organization schemes are all based on the tandem repetition of unit sequences. We previously showed that the S gamma 1 segment unexpectedly contains sequence elements that differ significantly from its prevalent unit repeat (49mer). Here we extend this preliminary characterization by determining the complete nucleotide sequence of the cloned S gamma 1 segment from BALB/c DNA. We find that S gamma 1 consists of more than 120 tandemly repeated 49mers. In addition, we show that the previously identified non-49mer sequences are part of a direct repeat element about 350 bp in length (DR II), which exists in two copies at the 5' end of S gamma 1. We also show that another unrelated direct repeat element about 500 bp long (DR I) exists near the 5' and 3' ends of S gamma 1. Thus, the structure of the S gamma 1 segment might be may be abbreviated as 5'-DRII-(49mer)15-DRI-DRII-(49mer)n-DRI , where n is between 40 and 160. Our results of Southern hybridization experiments suggest that this basic structural scheme is maintained in eight different Igh haplotypes, although S gamma 1 segments in different Igh haplotypes include different numbers of 49mer elements. Other murine S gamma segments differ in size among various Igh loci, but to a lesser extent than S gamma 1. At the level of tandemly repeated sequences, S gamma 1, S gamma 3, and S gamma 2b represent three distinct, nonoverlapping sets of sequences.     

Internal eliminated sequences interrupting the Oxytricha 81 locus: allelic divergence, conservation, conversions, and possible transposon origins

Internal eliminated sequences (IESs) often interrupt ciliate genes in the silent germline nucleus but are exactly excised and eliminated from the developing somatic nucleus from which genes are then expressed. Some long IESs are transposons, supporting the hypothesis that short IESs are ancient transposon relics. In light of that hypothesis and to explore the evolutionary history of a collection of IESs, we have compared various alleles of a particular locus (the 81 locus) of the ciliated protozoa Oxytricha trifallax and O. fallax. Three short IESs that interrupt two genes of the locus are found in alleles from both species, and thus must be relatively ancient, consistent with the hypothesis that short IESs are transposon relics. In contrast, TBE1 transposon interruptions of the locus are allele-specific and probably the results of recent transpositions. These IESs (and the TBE1s) are precisely excised from the DNA of the developing somatic macronucleus. Each IES interrupts a highly conserved sequence. A few nucleotides at the ends of each IES are also conserved, suggesting that they interact critically with IES excision machinery. However, most IES nucleotide positions have evolved at high rates, showing little or no selective constraint for function. Nonetheless, the length of each IES has been maintained (+/- 3 bp). While one IES is approximately 33 bp long, three other IESs have very similar sizes, approximately 70 bp long. Two IESs are surrounded by direct repeats of the sequence TTCTT. No other sequence similarities were found between any of the four IESs. However, the ends of one IES do match the inverted terminal repeat consensus sequence of the "TA" IESs of Paramecium. Three O. trifallax alleles appear to have been recipients in recent conversion events that could have been provoked by double-strand breaks associated with IES ends subsequent to IES transposition. Our findings support the hypothesis that short IESs evolved from ancient transposons that have lost most of their sequences, except those necessary for precise excision during macronuclear development.      

Characterization of proviruses cloned from mink cell focus-forming virus-infected cellular DNA

Two proviruses were cloned from EcoRI-digested DNA extracted from mink cells chronically infected with AKR mink cell focus-forming (MCF) 247 murine leukemia virus (MuLV), using a lambda phage host vector system. One cloned MuLV DNA fragment (designated MCF 1) contained sequences extending 6.8 kilobases from an EcoRI restriction site in the 5' long terminal repeat (LTR) to an EcoRI site located in the envelope (env) region and was indistinguishable by restriction endonuclease mapping for 5.1 kilobases (except for the EcoRI site in the LTR) from the 5' end of AKR ecotropic proviral DNA. The DNA segment extending from 5.1 to 6.8 kilobases contained several restriction sites that were not present in the AKR ecotropic provirus. A 0.5-kilobase DNA segment located at the 3' end of MCF 1 DNA contained sequences which hybridized to a xenotropic env-specific DNA probe but not to labeled ecotropic env-specific DNA. This dual character of MCF 1 proviral DNA was also confirmed by analyzing heteroduplex molecules by electron microscopy. The second cloned proviral DNA (designated MCF 2) was a 6.9-kilobase EcoRI DNA fragment which contained LTR sequences at each end and a 2.0-kilobase deletion encompassing most of the env region. The MCF 2 proviral DNA proved to be a useful reagent for detecting LTRs electron microscopically due to the presence of nonoverlapping, terminally located LTR sequences which effected its circularization with DNAs containing homologous LTR sequences. Nucleotide sequence analysis demonstrated the presence of a 104-base-pair direct repeat in the LTR of MCF 2 DNA. In contrast, only a single copy of the reiterated component of the direct repeat was present in MCF 1 DNA.     

In vitro HIV-1 selective integration into the target sequence and decoy-effect of the modified sequence

Although there have been a few reports that the HIV-1 genome can be selectively integrated into the genomic DNA of cultured host cell, the biochemistry of integration selectivity has not been fully understood. We modified the in vitro integration reaction protocol and developed a reaction system with higher efficiency. We used a substrate repeat, 5'-(GTCCCTTCCCAGT)(n)(ACTGGGAAGGGAC)(n)-3', and a modified sequence DNA ligated into a circular plasmid. CAGT and ACTG (shown in italics in the above sequence) in the repeat units originated from the HIV-1 proviral genome ends. Following the incubation of the HIV-1 genome end cDNA and recombinant integrase for the formation of the pre-integration (PI) complex, substrate DNA was reacted with this complex. It was confirmed that the integration selectively occurred in the middle segment of the repeat sequence. In addition, integration frequency and selectivity were positively correlated with repeat number n. On the other hand, both frequency and selectivity decreased markedly when using sequences with deletion of CAGT in the middle position of the original target sequence. Moreover, on incubation with the deleted DNAs and original sequence, the integration efficiency and selectivity for the original target sequence were significantly reduced, which indicated interference effects by the deleted sequence DNAs. Efficiency and selectivity were also found to vary discontinuously with changes in manganese dichloride concentration in the reaction buffer, probably due to its influence on the secondary structure of substrate DNA. Finally, integrase was found to form oligomers on the binding site and substrate DNA formed a loop-like structure. In conclusion, there is a considerable selectivity in HIV-integration into the specified sequence; however, similar DNA sequences can interfere with the integration process, and it is therefore difficult for in vivo integration to occur selectively in the actual host genome DNA.     

Hypotheses for testing deviations from random integration: evidence for nonrandom retroviral integration

The use of a recently developed in vitro model for retroviral integration provides a means of statistically testing hypotheses concerning the distribution of integration sites and hypotheses about the sequence of proviral orientations. In this study, three null hypotheses are formulated and applied to previously published data. Statistical analyses of these data suggest that the distribution of integration sites may not be uniform, and the sequence of proviral orientations is not random. On the basis of these results and the observed clustering of orientations, it was postulated that if a DNA sequence was involved in nonrandom proviral integration, that sequence would be found in the regions where the orientations change direction with respect to the target DNA ("I" regions). Computer analyses for homologous and complementary DNA sequences were performed on all possible pairs of identifiable "I" regions. A common sequence (at least 8 bp in size) was found in three out of four regions and that sequence was absent elsewhere in the target DNA. A model, with features of recombination reminiscent of chi sequences in bacteria, is proposed that may account for these results.     

Comparative genomics of mutualistic viruses of Glyptapanteles parasitic wasps

Background

Polydnaviruses, double-stranded DNA viruses with segmented genomes, have evolved as obligate endosymbionts of parasitoid wasps. Virus particles are replication deficient and produced by female wasps from proviral sequences integrated into the wasp genome. These particles are co-injected with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and, thereby, survival of proviral DNA. Here we characterize and compare the encapsidated viral genome sequences of bracoviruses in the family Polydnaviridae associated with Glyptapanteles gypsy moth parasitoids, along with near complete proviral sequences from which both viral genomes are derived.

Results

The encapsidated Glyptapanteles indiensis and Glyptapanteles flavicoxis bracoviral genomes, each composed of 29 different size segments, total approximately 517 and 594 kbp, respectively. They are generated from a minimum of seven distinct loci in the wasp genome. Annotation of these sequences revealed numerous novel features for polydnaviruses, including insect-like sugar transporter genes and transposable elements. Evolutionary analyses suggest that positive selection is widespread among bracoviral genes.

Conclusions

The structure and organization of G. indiensis and G. flavicoxis bracovirus proviral segments as multiple loci containing one to many viral segments, flanked and separated by wasp gene-encoding DNA, is confirmed. Rapid evolution of bracovirus genes supports the hypothesis of bracovirus genes in an 'arms race' between bracovirus and caterpillar. Phylogenetic analyses of the bracoviral genes encoding sugar transporters provides the first robust evidence of a wasp origin for some polydnavirus genes. We hypothesize transposable elements, such as those described here, could facilitate transfer of genes between proviral segments and host DNA.     

PCR-based gene disruption and recombinatory marker excision to produce modified industrial Saccharomyces cerevisiae without added sequences

The dominant selectable Kanr marker, which confers geneticin resistance in yeast, is extensively used for PCR based disruption of genes in functional analysis studies in laboratory strains of Saccharomyces cerevisiae. We have developed a gene disruption cassette, which incorporates the Kanr marker, and direct repeat sequences designed from the target gene to enable the deletion of the gene without the introduction of added DNA sequences. We report on the disruption of the HO gene as a test case, using the hodr-Kanr-hodr cassette. The cassette was shown to integrate at the HO locus and the Kanr marker excised by recombination between the two direct repeat sequences. The disruption/excision event resulted in the removal of one direct repeat and the coding sequence of the gene, and hence in this case loss of HO function, with the introduction of no foreign or additional sequences, including the Kanr marker. Having been derived from the target site, the remaining direct repeat sequence is native sequence in its native location. This design template has the potential to be adapted to other genes, and as such will be of advantage in instances such as the optimization of strains by recombinant DNA technology where the retention of minimal or no foreign sequences is desired.     

Characterization of a porcine variable number tandem repeat sequence specific for the glucosephosphate isomerase locus

A variable number of tandem repeat from a porcine glucosephosphate isomerase intron has been isolated and sequenced. The repeat has a unit size of 39 bp, is highly conserved and is present in at least 14 copies. Flanking sequences show a sequence periodicity of 53-54 bp and some sequence homology to the 39 bp repeat. A considerable part of the genomic DNA has been lost during subcloning and is considered to be deletion prone or refractory to propagation in E. coli. The tandem repeat is locus specific and detects at least six alleles in BamHI digested porcine DNA. No homology to other tandem repeat sequences has been found.     

Site-specific integration and excision of pMEA100 in Nocardia mediterranei

Summary Nocardia mediterranei strain LBG A3136 contains the 23.7 kb element pMEA100 in a chromosomally integrated form as well as in the free state (Moretti et al. 1985). The integrated form of this element can be excised precisely from the Nocardia chromosome without any accompanying rearrangements in flanking chromosomal DNA. After transfer into plasmid-free mutant strains, pMEA100 reintegrates site specifically into its original chromosomal locus. The exact mapping of the pMEA100 integration site was accomplished by restriction analysis and DNA sequencing. The attachment site of pMEA100, the junctions of its integrated form and plasmid-free chromosomal DNA of N. mediterranei contain an identical 47 bp long sequence which is probably required for site-specific recombination connected with integration and excision of pMEA100. Only one such sequence was found in the chromosome of pMEA100-free N. mediterranei derivatives as suggested by the single integration locus.     

Characterization of a molecularly cloned retroviral sequence associated with Fv-4 resistance.

The murine leukemia virus (MuLV) sequence associated with the resistance allele of the Fv-4 gene (Fv-4r) was molecularly cloned from genomic DNA of uninfected mice carrying this allele. The 5.2-kilobase cloned EcoRI DNA fragment (pFv4) was shown by nucleotide sequencing to contain 3.4 kilobases of a colinear MuLV-related proviral sequence which began in the C-terminal end of the pol region and extended through the env region and the 3' long terminal repeat. Cellular sequences flanked the 3' as well as the 5' ends of the truncated MuLV sequence. Alignment of the N-terminal half of the pFv4 env sequence with ecotropic, mink cell focus-forming, and xenotropic MuLV env sequences established the relatedness of pFv4 and ecotropic MuLV env sequences. A subcloned 700-base pair segment (pFv4env) from the 5' env region of pFv4 was used as an Fv-4-specific probe; it hybridized specifically to the Fv-4r-associated proviral sequence but not to endogenous ecotropic MuLV proviral DNA under high stringency. All Fv-4-resistant mice contained the same retroviral segment associated with the same flanking cellular DNA. Expression of Fv-4r-specific mRNA was demonstrated in the spleens of Fv-4r mice but not Fv-4s mice, supporting the previously proposed resistance model based on interference.     

Serial Alu sequence transposition interrupting a human B creatine kinase pseudogene.

We have isolated, sequenced, and characterized a single-copy B creatine kinase pseudogene. The chromosomal assignment of this gene is 16p13 and a unique sequence probe from this locus detects EcoRI restriction fragment length polymorphisms of 7.8 and 5.4 kb. In 26 unrelated individuals, the frequencies for the 7.8- and 5.4-kb B creatine kinase pseudogene alleles were calculated to be 17.3 and 82.7%, respectively. The B creatine kinase pseudogene is interrupted by a 904-bp DNA insertion composed of three Alu repeat sequences in tandem flanked by an 18-bp direct repeat, derived from the pseudogene sequence. Nucleotide sequence analysis of the Alu elements suggests that the Alu sequences were incorporated into this locus in three separate integration events. Several complex clustered Alu repeat sequences without defined integration borders have been previously identified at different genomic loci. This is the first evidence that complex tandem Alu elements can integrate in an apparently serial manner in the human genome and supports the contention that Alu repeats integrate nonrandomly into the human genome.     

A member of a new repeated sequence family which is conserved throughout eucaryotic evolution is found between the human delta and beta globin genes

A new class of human interspersed repeated sequences distinct from the AluI family was found by screening a human gene library with a mouse ribosomal gene non-transcribed spacer probe (rDNA NTS). A member of this sequence family was localized to a 251 bp segment between the human delta and beta globin genes: a region previously judged to be devoid of repeated DNA. The complete nucleotide sequence of this segment revealed a tandem block of 17 TG dinucleotides, a feature hypothesized by others to be a recombination hot spot responsible for gene conversion in the gamma globin locus region. When the genomes of Xenopus, pigeon, slime mold and yeast were examined, reiterated sequences homologous to both the mouse rDNA NTS and human globin repeat were found in every case. The discovery of this extraordinarily conserved repeated sequence family appears to have depended upon not using salmon sperm DNA during hybridization. The use of eucaryotic carrier DNA may bias the search for repeated sequences against any which may be highly conserved during eucaryotic evolution.