Free and integrated forms of hepatitis B virus DNA in human hepatocellular carcinoma cells (PLC/342) propagated in nude mice.

Primary hepatocellular carcinoma cells (PLC/342) propagated in nude mice produce hepatitis B surface antigen of subtype adr, as well as core particles containing viral DNA and DNA polymerase. Free and integrated forms of hepatitis B virus (HBV) DNA in the tumor were isolated by molecular cloning, and their nucleotide sequences were determined. Both of the two representative clones of free HBV DNA had the same genomic length (3,158 base pairs) and had two stop codons as well as two deletions in the envelope gene. None of the seven distinct clones of integrated HBV DNA possessed the entire viral genome. The integrated clone sequences had deletions and rearrangements, and only two clones possessed the envelope gene including the promoter and enhancer sequences. The C gene, which codes for core protein, was preserved in the two free clones and one of the integrated clones. The P gene, which codes for DNA polymerase, had deletions at two positions of 21 and 36 base pairs in both free clones, but was carried in toto by one of the integrated clones. The nucleotide sequences of the S genes of two free and four integrated clones, as well as their two inverted repeats, were compared. All of the eight sequences of the S gene possessed two nucleotide substitutions in common that were not displayed by any of the reported HBV genomes. The sequences differed from one another by only 1.2%. They differed, however, from 11 reported HBV genomes of subtype adr by 2.4%, from an ayr genome by 1.9%, from 2 adw genomes by 6.9%, and from 2 ayw genomes by 5.9%. These results indicate that all free and integrated HBV DNA species in the PLC/342 tumor cell evolved from a common progenitor. The free HBV DNA underwent nucleotide substitutions during several integration events, resulting in integrated HBV DNA copies that were similar in sequence but distinct from the reported HBV genomes.     

Repetitive satellite-like sequences are present within or upstream from 3 avian protein-coding genes.

Peculiar DNA sequences made up by the tandem repetition of a 5 bp unit have been identified within or upstream from three avian protein-coding genes. One sequence is located within an intron of the chicken "ovalbumin-X" gene with 5'-TCTCC-3' as basic repeat unit (36 repeats). Another sequence made of 27 repeats of a 5'-GGAAG-3' basic unit is found 2500 base pairs upstream from the promoter of the chicken ovotransferrin (conalbumin) gene. A related but different sequence is present in the corresponding region of the ovotransferrin gene in the pheasant, with 5'-GGAAA-3' as the basic unit (55 repeats). These three satellite-like elements are thus characterized by a total assymetry in base distribution, with purines restricted to one strand, and pyrimidines to the other. Two of the basic repeat units can be derived from the third one (GGAAA) by a single base pair change. These related sequences are found repeated in three avian genomes, at degrees which vary both with the sequence type and the genome type. Evolution of tandemly repeated sequences (including satellites) is in general studied by analysing randomly picked elements. The presence of conserved protein-coding regions neighbouring satellite-like sequences allow to follow their evolution at a single locus, as exemplified by the striking comparison of the pheasant and chicken sequences upstream from the ovotransferrin gene.     

Segmented Structure of Separate and Transposable DNA and RNA Elements as Suggested by Their Size Distributions

Abstract

A collection of about 1000 different eukaryotic and prokaryotic DNA mobile and separate elements is compiled from literature—;transposons, plasmids, extrachromosomal circular DNA, insertion sequences, as well as viral genomes and separate genome segments. Only small elements are collected, upto 2000 base pairs. Analysis of the sequence length distributions of the elements reveals that certain sizes are clearly preferred, namely those which correspond to multiples of about 345 bp in eukaryotes and multiples of about 210 bp in prokaryotes. This provides additional evidence in support of the theory (1) that segmented structure is characteristic of not only protein-coding sequences (2) but rather of genomes in general. In particular, it confirms the prediction (1) that mobile and separate elements would also be segmented.     

DNA, chromosomes, and in situ hybridization.

In situ hybridization is a powerful and unique technique that correlates molecular information of a DNA sequence with its physical location along chromosomes and genomes. It thus provides valuable information about physical map position of sequences and often is the only means to determine abundance and distribution of repetitive sequences making up the majority of most genomes. Repeated DNA sequences, composed of units of a few to a thousand base pairs in size, occur in blocks (tandem or satellite repeats) or are dispersed (including transposable elements) throughout the genome. They are often the most variable components of a genome, often being species and, occasionally, chromosome specific. Their variability arises through amplification, diversification and dispersion, as well as homogenization and loss; there is a remarkable correlation of molecular sequence features with chromosomal organization including the length of repeat units, their higher order structures, chromosomal locations, and dispersion mechanisms. Our understanding of the structure, function, organization, and evolution of genomes and their evolving repetitive components enabled many new cytogenetic applications to both medicine and agriculture, particularly in diagnosis and plant breeding.     

Base pair opening in three DNA-unwinding elements

DNA-unwinding elements are specific base sequences that are located in the origin of DNA replication where they provide the start point for strand separation and unwinding of the DNA double helix. In the present work we have obtained the first characterization of the opening of individual base pairs in DNA-unwinding elements. The three DNA molecules investigated reproduce the 13-mer DNA-unwinding elements present in the Escherichia coli chromosome. The base sequences of the three 13-mers are conserved in the origins of replication of enteric bacterial chromosomes. The exchange of imino protons with solvent protons was measured for each DNA as a function of the concentration of exchange catalyst using nuclear magnetic resonance spectroscopy. The exchange rates provided the rates and the equilibrium constants for opening of individual base pairs in each DNA at 20 degrees C. The results reveal that the kinetics and energetics of the opening reactions for AT/TA base pairs are different in the three DNA-unwinding elements due to long range effects of the base sequence. These differences encompass the AT/TA base pairs that are conserved in various bacterial genomes. Furthermore, a qualitative correlation is observed between the kinetics and energetics of opening of AT/TA base pairs and the location of the corresponding DNA-unwinding element in the origin of DNA replication.     

Tuning in to the signals: noncoding sequence conservation in vertebrate genomes

Aligning and comparing genomic sequences enables the identification of conserved sequence signatures and can enrich for coding and noncoding functional regions. In vertebrates, the comparison of human and rodent genomes and the comparison of evolutionarily distant genomes, such as human and pufferfish, have identified specific sets of 'ultraconserved' sequence elements associated with the control of early development. However, is this just the tip of a 'conservation iceberg' or do these sequences represent a specific class of regulatory element? Studies on the zebrafish phox2b gene region and the ENCODE project suggest that many regulatory elements are not highly conserved, posing intriguing questions about the relationship between noncoding sequence conservation and function and the evolution of regulatory sequences.     

Identification of a conserved sequence in the non-coding regions of many human genes.

We have analyzed a sequence of approximately 70 base pairs (bp) that shows a high degree of similarity to sequences present in the non-coding regions of a number of human and other mammalian genes. The sequence was discovered in a fragment of human genomic DNA adjacent to an integrated hepatitis B virus genome in cells derived from human hepatocellular carcinoma tissue. When one of the viral flanking sequences was compared to nucleotide sequences in GenBank, more than thirty human genes were identified that contained a similar sequence in their non-coding regions. The sequence element was usually found once or twice in a gene, either in an intron or in the 5' or 3' flanking regions. It did not share any similarities with known short interspersed nucleotide elements (SINEs) or presently known gene regulatory elements. This element was highly conserved at the same position within the corresponding human and mouse genes for myoglobin and N-myc, indicating evolutionary conservation and possible functional importance. Preliminary DNase I footprinting data suggested that the element or its adjacent sequences may bind nuclear factors to generate specific DNase I hypersensitive sites. The size, structure, and evolutionary conservation of this sequence indicates that it is distinct from other types of short interspersed repetitive elements. It is possible that the element may have a cis-acting functional role in the genome.     

MetaPlatanus: a metagenome assembler that combines long-range sequence links and species-specific features

De novo metagenome assembly is effective in assembling multiple draft genomes, including those of uncultured organisms. However, heterogeneity in the metagenome hinders assembly and introduces interspecies misassembly deleterious for downstream analysis. For this purpose, we developed a hybrid metagenome assembler, MetaPlatanus. First, as a characteristic function, it assembles the basic contigs from accurate short reads and then iteratively utilizes long-range sequence links, species-specific sequence compositions, and coverage depth. The binning information was also used to improve contiguity. Benchmarking using mock datasets consisting of known bacteria with long reads or mate pairs revealed the high contiguity MetaPlatanus with a few interspecies misassemblies. For published human gut data with nanopore reads from potable sequencers, MetaPlatanus assembled many biologically important elements, such as coding genes, gene clusters, viral sequences, and over-half bacterial genomes. In the benchmark with published human saliva data with high-throughput nanopore reads, the superiority of MetaPlatanus was considerably more evident. We found that some high-abundance bacterial genomes were assembled only by MetaPlatanus as near-complete. Furthermore, MetaPlatanus can circumvent the limitations of highly fragmented assemblies and frequent interspecies misassembles obtained by the other tools. Overall, the study demonstrates that MetaPlatanus could be an effective approach for exploring large-scale structures in metagenomes.     

Analysis of regulatory elements involved in the induction of two tobacco genes by salicylate treatment and virus infection.

Tobacco genes encoding the PR-1a protein and a glycine-rich protein are expressed after treatment of plants with salicylate or infection with tobacco mosaic virus. Upstream sequences of these genes were fused to reporter genes, and these constructs were used to transform tobacco. Upstream sequences of the PR-1a gene of 689 base pairs or longer were sufficient for induction of the reporter gene in tobacco mosaic virus-inoculated leaves, systemically induced leaves from infected plants, and leaves treated with salicylate. No such induction was found with upstream sequences of 643 base pairs or shorter of the PR-1a gene. When the PR-1a upstream sequence from nucleotides -625 to -902 was fused to the cauliflower mosaic virus 35S core promoter, a construct was obtained that conferred tobacco mosaic virus and salicylate inducibility to the reporter gene in transgenic plants. This confirmed the localization of tobacco mosaic virus- and salicylate-responsive elements between positions -643 and -689 in the PR-1a promoter. With the glycine-rich protein gene, an upstream sequence of 645 base pairs was sufficient for tobacco mosaic virus and salicylate inducibility of the reporter gene, whereas constructs containing 400 base pairs or fewer of the glycine-rich protein promoter were largely inactive.     

Characterization of genomic poly(dT-dG).poly(dC-dA) sequences: structure, organization, and conformation.

Hybridization studies suggest the abundant presence of poly(dT-dG).poly(dC-dA) (TG-element), a potential Z-DNA sequence, in eucaryotic genomes. We have isolated and characterized TG-elements from different locations in the human genome: from randomly isolated clones, associated with the actin gene family, and linked to another repeated element. The results indicate that the following features are typical of these TG-elements: the elements consist of 20 to 60 base pairs of (dT-dG)n.(dC-dA)n, the sequences characterized in our study were not flanked by direct or inverted repeats, the sequences are interspersed rather than in satellite blocks, the elements are not usually associated with other repeated elements, and some of the elements are found near coding sequences or in introns. Studies on the conformation of a genomic TG-element in a supercoiled plasmid indicate several distinct properties of the TG-element: it is in the Z-form only at low ionic strength, S1 nuclease recognizes its Z-form with a marked preference for one of the B-Z junctions, and the sensitive region extends for 20 base pairs near the B-Z junction. In contrast to the result with the supercoiled plasmid, S1 nuclease failed to recognize the TG-element in minichromosomes.     

Molecular cloning of Snyder-Theilen feline leukemia and sarcoma viruses: comparative studies of feline sarcoma virus with its natural helper virus and with Moloney murine sarcoma virus

Extrachromosomal DNA obtained from mink cells acutely infected with the Snyder-Theilen (ST) strain of feline sarcoma virus (feline leukemia virus) [FeSV(FeLV)] was fractionated electrophoretically, and samples enriched for FeLV and FeSV linear intermediates were digested with EcoRI and cloned in lambda phage. Hybrid phages were isolated containing either FeSV or FeLV DNA "inserts" and were characterized by restriction enzyme analysis, R-looping with purified 26 to 32S viral RNA, and heteroduplex formation. The recombinant phages (designated lambda FeSV and lambda FeLV) contain all of the genetic information represented in FeSV and FeLV RNA genomes but lack one extended terminally redundant sequence of 750 bases which appears once at each end of parental linear DNA intermediates. Restriction enzyme and heteroduplex analyses confirmed that sequences unique to FeSV (src sequences) are located at the center of the FeSV genome and are approximately 1.5 kilobase pairs in length. With respect to the 5'-3' orientation of genes in viral RNA, the order of genes in the FeSV genome is 5'-gag-src-env-c region-3'; only 0.9 kilobase pairs of gag and 0.6 kilobase pairs of env-derived FeLV sequences are represented in ST FeSV. Heteroduplex analyses between lambda FeSV or lambda FeLV DNA and Moloney murine sarcoma virus DNA (strain m1) were performed under conditions of reduced stringency to demonstrate limited regions of base pair homology. Two such regions were identified: the first occurs at the extreme 5' end of the leukemia and both sarcoma viral genomes, whereas the second corresponds to a 5' segment of leukemia virus "env" sequences conserved in both sarcoma viruses. The latter sequences are localized at the 3' end of FeSV src and at the 5' end of murine sarcoma virus src and could possibly correspond to regions of helper virus genomes that are required for retroviral transforming functions.     

The p92 polymerase coding region contains an internal RNA element required at an early step in Tombusvirus genome replication

The replication of positive-strand RNA viral genomes involves various cis-acting RNA sequences. Generally, regulatory RNA sequences are present at or near genomic termini; however, internal replication elements (IREs) also exist. Here we report the structural and functional characterization of an IRE present in the readthrough portion of the p92 polymerase gene of Tomato bushy stunt virus. Analysis of this element in the context of a noncoding defective interfering RNA revealed a functional core structure composed of two noncontiguous segments of sequence that interact with each other to form an extended helical conformation. IRE activity required maintenance of several base-paired sections as well as two distinct structural features: (i) a short, highly conserved segment that can potentially form two different and mutually exclusive structures and (ii) an internal loop that contains a critical CC mismatch. The IRE was also shown to play an essential role within the context of the viral genome. In vivo analysis with novel RNA-based temperature-sensitive genomic mutants and translationally active subgenomic viral replicons revealed the following about the IRE: (i) it is active in the positive strand, (ii) it is dispensable late in the viral RNA replication process, and (iii) it is functionally inhibited by active translation over its sequence. Together, these results suggest that IRE activity is required in the cytosol at an early step in the viral replication process, such as template recruitment and/or replicase complex assembly.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

Nucleotide sequence analysis of avian retroviruses: structural similarities with transposable elements

Integrated retroviral genomes are flanked by direct repeats of sequences derived from the termini of the viral RNA genome. These sequences are designated long terminal repeats (LTRs). We have determined and analyzed the nucleotide sequence of the LTRs from several exogenous and endogenous avian retroviruses. These LTRs possess several structural similarities with eukaryotic and prokaryotic transposable elements: 1) inverted complementary repeats at the termini, 2) deletions of sequences adjacent to the LTR, 3) small duplications of host sequences flanking the integrated provirus, and 4) sequence homologies with transposable and other genetic elements. These observations suggest that LTRs function in the integration and perhaps transposition of retrovirus genomes. Evidence exists for the presence of a strong promoter sequence within the LTR. The retroviral LTR also contains a "Hogness box" up-stream of the capping site and a poly(A) signal. These features suggest an additional role for the LTR in the regulation of gene expression.     

Vertebrate DNAs contain nucleotide sequences related to the transforming gene of avian myeloblastosis virus.

Avian myeloblastosis virus contains a continuous sequence of approximately 1,000 nucleotides which may represent a gene (amv) responsible for acute myeloblastic leukemia in chickens. This sequence appears to have been acquired from chicken DNA and to be substituted for the envelope gene in the viral genome. We used hybridization probes enriched for the amv sequences and conditions that facilitate annealing of partially homologous nucleotide sequences to show that cellular sequences related to amv are present in the genomes of all vertebrates ranging from amphibians to humans but were not detected in fish, sea urchins, or Escherichia coli. In contrast to the preceding findings, nontransforming endogenous proviral nucleotide sequences closely related to the remainder of the avian myeloblastosis virus genome and to the entire myeloblastosis-associated helper virus are present only in chicken DNA. The amv-related cellular sequences appear to be highly conserved during evolution and to be contained at only one or a few locations in the genome of vertebrates. Within closely related species, they appear to share common evolutionary genetic loci. These findings and similar ones obtained with other highly oncogenic retroviruses containing a transforming gene suggest a general mechanism for acquisition of viral oncogenic sequences and an essential role for these sequences in the normal cellular state.