Paired-End Sequencing of Long-Range DNA Fragments for De Novo Assembly of Large,Complex Mammalian Genomes by Direct Intra-Molecule Ligation

Background

The relatively short read lengths from next generation sequencing (NGS) technologies still pose a challenge for de novo assembly of complex mammal genomes. One important solution is to use paired-end (PE) sequence information experimentally obtained from long-range DNA fragments (>1 kb). Here, we characterize and extend a long-range PE library construction method based on direct intra-molecule ligation (or molecular linker-free circularization) for NGS.

Results

We found that the method performs stably for PE sequencing of 2- to 5- kb DNA fragments, and can be extended to 10–20 kb (and even in extremes, up to ∼35 kb). We also characterized the impact of low quality input DNA on the method, and develop a whole-genome amplification (WGA) based protocol using limited input DNA (<1 µg). Using this PE dataset, we accurately assembled the YanHuang (YH) genome, the first sequenced Asian genome, into a scaffold N50 size of >2 Mb, which is over100-times greater than the initial size produced with only small insert PE reads(17 kb). In addition, we mapped two 7- to 8- kb insertions in the YH genome using the larger insert sizes of the long-range PE data.

Conclusions

In conclusion, we demonstrate here the effectiveness of this long-range PE sequencing method and its use for the de novo assembly of a large, complex genome using NGS short reads.     

Amplification Biases and Consistent Recovery of Loci in a Double-Digest RAD-seq Protocol

A growing variety of “genotype-by-sequencing” (GBS) methods use restriction enzymes and high throughput DNA sequencing to generate data for a subset of genomic loci, allowing the simultaneous discovery and genotyping of thousands of polymorphisms in a set of multiplexed samples. We evaluated a “double-digest” restriction-site associated DNA sequencing (ddRAD-seq) protocol by 1) comparing results for a zebra finch (Taeniopygia guttata) sample with in silico predictions from the zebra finch reference genome; 2) assessing data quality for a population sample of indigobirds (Vidua spp.); and 3) testing for consistent recovery of loci across multiple samples and sequencing runs. Comparison with in silico predictions revealed that 1) over 90% of predicted, single-copy loci in our targeted size range (178–328 bp) were recovered; 2) short restriction fragments (38–178 bp) were carried through the size selection step and sequenced at appreciable depth, generating unexpected but nonetheless useful data; 3) amplification bias favored shorter, GC-rich fragments, contributing to among locus variation in sequencing depth that was strongly correlated across samples; 4) our use of restriction enzymes with a GC-rich recognition sequence resulted in an up to four-fold overrepresentation of GC-rich portions of the genome; and 5) star activity (i.e., non-specific cutting) resulted in thousands of “extra” loci sequenced at low depth. Results for three species of indigobirds show that a common set of thousands of loci can be consistently recovered across both individual samples and sequencing runs. In a run with 46 samples, we genotyped 5,996 loci in all individuals and 9,833 loci in 42 or more individuals, resulting in <1% missing data for the larger data set. We compare our approach to similar methods and discuss the range of factors (fragment library preparation, natural genetic variation, bioinformatics) influencing the recovery of a consistent set of loci among samples.     

Identification of acquired DNA in Neisseria lactamica

Anomalous DNA (aDNA) in prokaryotic genomes, identified by its aberrant nucleotide composition, generally represents horizontally acquired DNA. Previous studies showed that frequent DNA transfer occurs between commensal Neisseriae and Neisseria meningitidis. Currently, it is unknown whether aDNA regions are also transferred between these species. The genome of Neisseria lactamica strain 892586 was assessed by a strategy that enables the selective isolation of aDNA, using endonucleases with recognition sites that are overrepresented in aDNA. Of eight regions with aDNA, five displayed similarity to virulence-associated meningococcal sequences. Of three aDNA fragments with limited or no similarity to neisserial sequences, one encodes a novel putative autotransporter/adhesin. The remaining two fragments are adjacent in the N. lactamica genome, and encode a novel putative ATPase/subtilisin-like protease operon. A similar operon is present in the genomes of different respiratory tract pathogens. The identification of aDNA from N. lactamica with similarity to meningococcal aDNA shows that genetic exchange between the Neisseriae is not limited to the neisserial core genome. The discovery of aDNA in N. lactamica similar to a locus in other pathogens substantially expands the neisserial gene pool.     

Genomics and genetics of Sulfolobus islandicus LAL14/1, a model hyperthermophilic archaeon

The 2 465 177 bp genome of Sulfolobus islandicus LAL14/1, host of the model rudivirus SIRV2, was sequenced. Exhaustive comparative genomic analysis of S. islandicus LAL14/1 and the nine other completely sequenced S. islandicus strains isolated from Iceland, Russia and USA revealed a highly syntenic common core genome of approximately 2 Mb and a long hyperplastic region containing most of the strain-specific genes. In LAL14/1, the latter region is enriched in insertion sequences, CRISPR (clustered regularly interspaced short palindromic repeats), glycosyl transferase genes, toxin–antitoxin genes and MITE (miniature inverted-repeat transposable elements). The tRNA genes of LAL14/1 are preferential targets for the integration of mobile elements but clusters of atypical genes (CAG) are also integrated elsewhere in the genome. LAL14/1 carries five CRISPR loci with 10 per cent of spacers matching perfectly or imperfectly the genomes of archaeal viruses and plasmids found in the Icelandic hot springs. Strikingly, the CRISPR_2 region of LAL14/1 carries an unusually long 1.9 kb spacer interspersed between two repeat regions and displays a high similarity to pING1-like conjugative plasmids. Finally, we have developed a genetic system for S. islandicus LAL14/1 and created ΔpyrEF and ΔCRISPR_1 mutants using double cross-over and pop-in/pop-out approaches, respectively. Thus, LAL14/1 is a promising model to study virus–host interactions and the CRISPR/Cas defence mechanism in Archaea.     

Pulling out the 1%: Whole-Genome Capture for the Targeted Enrichment of Ancient DNA Sequencing Libraries

Most ancient specimens contain very low levels of endogenous DNA, precluding the shotgun sequencing of many interesting samples because of cost. Ancient DNA (aDNA) libraries often contain <1% endogenous DNA, with the majority of sequencing capacity taken up by environmental DNA. Here we present a capture-based method for enriching the endogenous component of aDNA sequencing libraries. By using biotinylated RNA baits transcribed from genomic DNA libraries, we are able to capture DNA fragments from across the human genome. We demonstrate this method on libraries created from four Iron Age and Bronze Age human teeth from Bulgaria, as well as bone samples from seven Peruvian mummies and a Bronze Age hair sample from Denmark. Prior to capture, shotgun sequencing of these libraries yielded an average of 1.2% of reads mapping to the human genome (including duplicates). After capture, this fraction increased substantially, with up to 59% of reads mapped to human and enrichment ranging from 6- to 159-fold. Furthermore, we maintained coverage of the majority of regions sequenced in the precapture library. Intersection with the 1000 Genomes Project reference panel yielded an average of 50,723 SNPs (range 3,062–147,243) for the postcapture libraries sequenced with 1 million reads, compared with 13,280 SNPs (range 217–73,266) for the precapture libraries, increasing resolution in population genetic analyses. Our whole-genome capture approach makes it less costly to sequence aDNA from specimens containing very low levels of endogenous DNA, enabling the analysis of larger numbers of samples.     

How Hot Are Drosophila Hotspots? Examining Recombination Rate Variation and Associations with Nucleotide Diversity,Divergence, and Maternal Age in Drosophila pseudoobscura

Fine scale meiotic recombination maps have uncovered a large amount of variation in crossover rate across the genomes of many species, and such variation in mammalian and yeast genomes is concentrated to <5kb regions of highly elevated recombination rates (10–100x the background rate) called “hotspots.” Drosophila exhibit substantial recombination rate heterogeneity across their genome, but evidence for these highly-localized hotspots is lacking. We assayed recombination across a 40Kb region of Drosophila pseudoobscura chromosome 2, with one 20kb interval assayed every 5Kb and the adjacent 20kb interval bisected into 10kb pieces. We found that recombination events across the 40kb stretch were relatively evenly distributed across each of the 5kb and 10kb intervals, rather than concentrated in a single 5kb region. This, in combination with other recent work, indicates that the recombination landscape of Drosophila may differ from the punctate recombination pattern observed in many mammals and yeast. Additionally, we found no correlation of average pairwise nucleotide diversity and divergence with recombination rate across the 20kb intervals, nor any effect of maternal age in weeks on recombination rate in our sample.     

Differential detection of transposable elements between Saccharum species

Cultivars of sugarcane (Saccharum) are hybrids between species S. officinarum (x = 10, 2n = 8x = 80) and S. spontaneum (x = 8, 2n = 5 – 16x = 40 – 128). These accessions have 100 to 130 chromosomes, 80–85% of which are derived from S. officinarum, 10–15% from S. spontaneum, and 5–10% are possible recombinants between the two genomes. The aim of this study was to analyze the repetition of DNA sequences in S. officinarum and S. spontaneum. For this purpose, genomic DNA from S. officinarum was digested with restriction enzymes and the fragments cloned. Sixty-eight fragments, approximately 500 bp, were cloned, sequenced and had their identity analyzed in NCBI, and in the rice, maize, and sorghum genome databases using BLAST. Twelve clones containing partial transposable elements, one single-copy control, one DNA repetitive clone control and two genome controls were analyzed by DNA hybridization on membrane, using genomic probes from S. officinarum and S. spontaneum. The hybridization experiment revealed that six TEs had a similar repetitive DNA pattern in the genomes of S. officinarum and S. spontaneum, while six TEs were more abundant in the genome of S. officinarum. We concluded that the species S. officinarum and S. spontaneum have differential accumulation LTR retrotransposon families, suggesting distinct insertion or modification patterns.     

Fragmented mitochondrial genomes are present in both major clades of the blood-sucking lice (suborder Anoplura): evidence from two Hoplopleura rodent lice (family Hoplopleuridae)

Background

The suborder Anoplura contains 540 species of blood-sucking lice that parasitize over 840 species of eutherian mammals. Fragmented mitochondrial (mt) genomes have been found in the lice of humans, pigs, horses and rats from four families: Pediculidae, Pthiridae, Haematopinidae and Polyplacidae. These lice, eight species in total, are from the same major clade of the Anoplura. The mt genomes of these lice consist of 9–20 minichromosomes; each minichromosome is 1.5–4 kb in size and has 1–8 genes. To understand mt genome fragmentation in the other major clade of the Anoplura, we sequenced the mt genomes of two species of rodent lice in the genus Hoplopleura (family Hoplopleuridae).

Results

We identified 28 mt genes on 10 minichromosomes in the mouse louse, Ho. akanezumi; each minichromosome is 1.7–2.7 kb long and has 1–6 genes. We identified 34 mt genes on 11 minichromosomes in the rat louse, Ho. kitti; each minichromosome is 1.8–2.8 kb long and has 1–5 genes. Ho. akanezumi also has a chimeric minichromosome with parts of two rRNA genes and a full-length tRNA gene for tyrosine. These two rodent lice share the same pattern for the distribution of all of the protein-coding and rRNA genes but differ in tRNA gene content and gene arrangement in four minichromosomes. Like the four genera of blood-sucking lice that have been investigated in previous studies, the Hoplopleura species have four minichromosomes that are only found in this genus.

Conclusions

Our results indicate that fragmented mt genomes were present in the most recent common ancestor of the two major clades of the blood-sucking lice, which lived ~75 million years ago. Intra-genus variation in the pattern of mt genome fragmentation is common in the blood-sucking lice (suborder Anoplura) and genus-specific minichromosomes are potential synapomorphies. Future studies should expand into more species, genera and families of blood-sucking lice to explore further the phylogenetic utility of the novel features associated with fragmented mt genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-751) contains supplementary material, which is available to authorized users.     

Meningococcal genetic variation mechanisms viewed through comparative analysis of serogroup C strain FAM18

The bacterium Neisseria meningitidis is commonly found harmlessly colonising the mucosal surfaces of the human nasopharynx. Occasionally strains can invade host tissues causing septicaemia and meningitis, making the bacterium a major cause of morbidity and mortality in both the developed and developing world. The species is known to be diverse in many ways, as a product of its natural transformability and of a range of recombination and mutation-based systems. Previous work on pathogenic Neisseria has identified several mechanisms for the generation of diversity of surface structures, including phase variation based on slippage-like mechanisms and sequence conversion of expressed genes using information from silent loci. Comparison of the genome sequences of two N. meningitidis strains, serogroup B MC58 and serogroup A Z2491, suggested further mechanisms of variation, including C-terminal exchange in specific genes and enhanced localised recombination and variation related to repeat arrays. We have sequenced the genome of N. meningitidis strain FAM18, a representative of the ST-11/ET-37 complex, providing the first genome sequence for the disease-causing serogroup C meningococci; it has 1,976 predicted genes, of which 60 do not have orthologues in the previously sequenced serogroup A or B strains. Through genome comparison with Z2491 and MC58 we have further characterised specific mechanisms of genetic variation in N. meningitidis, describing specialised loci for generation of cell surface protein variants and measuring the association between noncoding repeat arrays and sequence variation in flanking genes. Here we provide a detailed view of novel genetic diversification mechanisms in N. meningitidis. Our analysis provides evidence for the hypothesis that the noncoding repeat arrays in neisserial genomes (neisserial intergenic mosaic elements) provide a crucial mechanism for the generation of surface antigen variants. Such variation will have an impact on the interaction with the host tissues, and understanding these mechanisms is important to aid our understanding of the intimate and complex relationship between the human nasopharynx and the meningococcus.     

Analysis of six prophages in Lactococcus lactis IL1403: different genetic structure of temperate and virulent phage populations

We report the genetic organisation of six prophages present in the genome of Lactococcus lactis IL1403. The three larger prophages (36–42 kb), belong to the already described P335 group of temperate phages, whereas the three smaller ones (13–15 kb) are most probably satellites relying on helper phage(s) for multiplication. These data give a new insight into the genetic structure of lactococcal phage populations. P335 temperate phages have variable genomes, sharing homology over only 10–33% of their length. In contrast, virulent phages have highly similar genomes sharing homology over >90% of their length. Further analysis of genetic structure in all known groups of phages active on other bacterial hosts such as Escherichia coli, Bacillus subtilis, Mycobacterium and Streptococcus thermophilus confirmed the existence of two types of genetic structure related to the phage way of life. This might reflect different intensities of horizontal DNA exchange: low among purely virulent phages and high among temperate phages and their lytic homologues. We suggest that the constraints on genetic exchange among purely virulent phages reflect their optimal genetic organisation, adapted to a more specialised and extreme form of parasitism than temperate/lytic phages.     

Comparative Genometrics (CG): a database dedicated to biometric comparisons of whole genomes

The ever increasing rate at which whole genome sequences are becoming accessible to the scientific community has created an urgent need for tools enabling comparison of chromosomes of different species. We have applied biometric methods to available chromosome sequences and posted the results on our Comparative Genometrics (CG) web site. By genometrics, a term coined by Elston and Wilson [Genet. Epidemiol. (1990), 7, 17–19], we understand a biometric analysis of chromosomes. During the initial phase, our web site displays, for all completely sequenced prokaryotic genomes, three genometric analyses: the DNA walk [Lobry (1999) Microbiology Today, 26, 164–165] and two complementary representations, i.e. the cumulative GC- and TA-skew analyses, capable of identifying, at the level of whole genomes, features inherent to chromosome organization and functioning. It appears that the latter features are taxon-specific. Although primarily focused on prokaryotic chromosomes, the CG web site contains genometric information on paradigm plasmids, phages, viruses and eukaryotic organelles. Relevant data and methods can be readily used by the scientific community for further analyses as well as for tutorial purposes. Our data posted at the CG web site are freely available on the World Wide Web at http://www.unil.ch/comparativegenometrics.     

The Genetic Properties of the Primary Endosymbionts of Mealybugs Differ from Those of Other Endosymbionts of Plant Sap-Sucking Insects

Mealybugs (Hemiptera, Coccoidea, Pseudococcidae), like aphids and psyllids, are plant sap-sucking insects that have an obligate association with prokaryotic endosymbionts that are acquired through vertical, maternal transmission. We sequenced two fragments of the genome of Tremblaya princeps, the endosymbiont of mealybugs, which is a member of the β subdivision of the Proteobacteria. Each of the fragments (35 and 30 kb) contains a copy of 16S-23S-5S rRNA genes. A total of 37 open reading frames were detected, which corresponded to putative rRNA proteins, chaperones, and enzymes of branched-chain amino acid biosynthesis, DNA replication, protein translation, and RNA synthesis. The genome of T. princeps has a number of properties that distinguish it from the genomes of Buchnera aphidicola and Carsonella ruddii, the endosymbionts of aphids and psyllids, respectively. Among these properties are a high G+C content (57.1 mol%), the same G+C content in intergenic spaces and structural genes, and similar G+C contents of the genes encoding highly and poorly conserved proteins. The high G+C content has a substantial effect on protein composition; about one-third of the residues consist of four amino acids with high-G+C-content codons. Sequence analysis of DNA fragments containing the rRNA operon and adjacent regions from endosymbionts of several mealybug species suggested that there was a single duplication of the rRNA operon and the adjacent genes in an ancestor of the present T. princeps. Subsequently, in one mealybug lineage rpS15, one of the duplicated genes, was retained, while in another lineage it decayed. These results extend the diversity of the types of endosymbiotic associations found in plant sap-sucking insects.     

High-throughput,cost-effective verification of structural DNA assembly

DNA ‘assembly’ from ‘building blocks’ remains a cornerstone in synthetic biology, whether it be for gene synthesis (∼1 kb), pathway engineering (∼10 kb) or synthetic genomes (>100 kb). Despite numerous advances in the techniques used for DNA assembly, verification of the assembly is still a necessity, which becomes cost-prohibitive and a logistical challenge with increasing scale. Here we describe for the first time a comprehensive, high-throughput solution for structural DNA assembly verification by restriction digest using exhaustive in silico enzyme screening, rolling circle amplification of plasmid DNA, capillary electrophoresis and automated digest pattern recognition. This low-cost and robust methodology has been successfully used to screen over 31 000 clones of DNA constructs at <$1 per sample.     

Genome Wide Search for Biomarkers to Diagnose Yersinia Infections

Bacterial identification on the basis of the highly conserved 16S rRNA (rrs) gene is limited by its presence in multiple copies and a very high level of similarity among them. The need is to look for other genes with unique characteristics to be used as biomarkers. Fifty-one sequenced genomes belonging to 10 different Yersinia species were used for searching genes common to all the genomes. Out of 304 common genes, 34 genes of sizes varying from 0.11 to 4.42 kb, were selected and subjected to in silico digestion with 10 different Restriction endonucleases (RE) (4–6 base cutters). Yersinia species have 6–7 copies of rrs per genome, which are difficult to distinguish by multiple sequence alignments or their RE digestion patterns. However, certain unique combinations of other common gene sequences—carB, fadJ, gluM, gltX, ileS, malE, nusA, ribD, and rlmL and their RE digestion patterns can be used as markers for identifying 21 strains belonging to 10 Yersinia species: Y. aldovae, Y. enterocolitica, Y. frederiksenii, Y. intermedia, Y. kristensenii, Y. pestis, Y. pseudotuberculosis, Y. rohdei, Y. ruckeri, and Y. similis. This approach can be applied for rapid diagnostic applications.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0552-6) contains supplementary material, which is available to authorized users.     

A novel strategy for the identification of genomic islands by comparative analysis of the contents and contexts of tRNA sites in closely related bacteria

We devised software tools to systematically investigate the contents and contexts of bacterial tRNA and tmRNA genes, which are known insertion hotspots for genomic islands (GIs). The strategy, based on MAUVE-facilitated multigenome comparisons, was used to examine 87 Escherichia coli MG1655 tRNA and tmRNA genes and their orthologues in E.coli EDL933, E.coli CFT073 and Shigella flexneri Sf301. Our approach identified 49 GIs occupying ~1.7 Mb that mapped to 18 tRNA genes, missing 2 but identifying a further 30 GIs as compared with Islander [Y. Mantri and K. P. Williams (2004), Nucleic Acids Res., 32, D55–D58]. All these GIs had many strain-specific CDS, anomalous GC contents and/or significant dinucleotide biases, consistent with foreign origins. Our analysis demonstrated marked conservation of sequences flanking both empty tRNA sites and tRNA-associated GIs across all four genomes. Remarkably, there were only 2 upstream and 5 downstream deletions adjacent to the 328 loci investigated. In silico PCR analysis based on conserved flanking regions was also used to interrogate hotspots in another eight completely or partially sequenced E.coli and Shigella genomes. The tools developed are ideal for the analysis of other bacterial species and will lead to in silico and experimental discovery of new genomic islands.     

Genome Wide Analysis for Rapid Identification of Vibrio Species

The highly conserved 16S rRNA (rrs) gene is generally used for bacterial identification. In organisms possessing multiple copies of rrs, high intra-genomic heterogeneity does not allow easy distinction among different species. In order to identify Vibrio species, a wide range of genes have been employed. There is an urgent requirement of a consensus gene, which can be used as biomarker for rapid identification. Eight sequenced genomes of Vibrio species were screened for selecting genes which were common among all the genomes. Out of 108 common genes, 24 genes of sizes varying from 0.11 to 3.94 kb were subjected to in silico digestion with 10 type II restriction endonucleases (RE). A few unique genes—dapF, fadA, hisD, ilvH, lpxC, recF, recR, rph and ruvB in combination with certain REs provided unique digestion patterns, which can be used as biomarkers. This protocol can be exploited for rapid diagnosis of Vibrio species.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0553-5) contains supplementary material, which is available to authorized users.     

Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39

The genome sequences of Chlamydia trachomatis mouse pneumonitis (MoPn) strain Nigg (1 069 412 nt) and Chlamydia pneumoniae strain AR39 (1 229 853 nt) were determined using a random shotgun strategy. The MoPn genome exhibited a general conservation of gene order and content with the previously sequenced C.trachomatis serovar D. Differences between C.trachomatis strains were focused on an ~50 kb ‘plasticity zone’ near the termination origins. In this region MoPn contained three copies of a novel gene encoding a >3000 amino acid toxin homologous to a predicted toxin from Escherichia coli 0157:H7 but had apparently lost the tryptophan biosyntheis genes found in serovar D in this region. The C.pneumoniae AR39 chromosome was >99.9% identical to the previously sequenced C.pneumoniae CWL029 genome, however, comparative analysis identified an invertible DNA segment upstream of the uridine kinase gene which was in different orientations in the two genomes. AR39 also contained a novel 4524 nt circular single-stranded (ss)DNA bacteriophage, the first time a virus has been reported infecting C.pneumoniae. Although the chlamydial genomes were highly conserved, there were intriguing differences in key nucleotide salvage pathways: C.pneumoniae has a uridine kinase gene for dUTP production, MoPn has a uracil phosphororibosyl transferase, while C.trachomatis serovar D contains neither gene. Chromosomal comparison revealed that there had been multiple large inversion events since the species divergence of C.trachomatis and C.pneumoniae, apparently oriented around the axis of the origin of replication and the termination region. The striking synteny of the Chlamydia genomes and prevalence of tandemly duplicated genes are evidence of minimal chromosome rearrangement and foreign gene uptake, presumably owing to the ecological isolation of the obligate intracellular parasites. In the absence of genetic analysis, comparative genomics will continue to provide insight into the virulence mechanisms of these important human pathogens.     

Editor's choice: AnABlast: a new in silico strategy for the genome-wide search of novel genes and fossil regions

Genome annotation, assisted by computer programs, is one of the great advances in modern biology. Nevertheless, the in silico identification of small and complex coding sequences is still challenging. We observed that amino acid sequences inferred from coding—but rarely from non-coding—DNA sequences accumulated alignments in low-stringency BLAST searches, suggesting that this alignments accumulation could be used to highlight coding regions in sequenced DNA. To investigate this possibility, we developed a computer program (AnABlast) that generates profiles of accumulated alignments in query amino acid sequences using a low-stringency BLAST strategy. To validate this approach, all six-frame translations of DNA sequences between every two annotated exons of the fission yeast genome were analysed with AnABlast. AnABlast-generated profiles identified three new copies of known genes, and four new genes supported by experimental evidence. New pseudogenes, ancestral carboxyl- and amino-terminal subtractions, complex gene rearrangements, and ancient fragments of mitDNA and of bacterial origin, were also inferred. Thus, this novel in silico approach provides a powerful tool to uncover new genes, as well as fossil-coding sequences, thus providing insight into the evolutionary history of annotated genomes.     

Molecular Cloning of the Fujinami Sarcoma Virus Genome and Its Comparison with Sequences of Other Related Transforming Viruses

Full-length proviral DNA of Fujinami sarcoma virus (FSV) of chickens was molecularly cloned and characterized. An analysis of FSV DNA integrated in mammalian cells showed that restriction endonuclease SacI has a single cleavage site on FSV DNA. Unintegrated closed circular FSV DNA obtained from newly infected cells was linearized by digestion with SacI and cloned into λgtWES·λB. The following three different molecules were isolated: FSV-1 (4.4 kilobases [kb]) and FSV-2 (4.7 kb), which appeared to be full-length FSV DNA molecules containing either one or two copies of the long terminal repeat structure, and FSV-3 (6 kb), which consisted of part FSV DNA and part DNA of unknown origin. An analysis of the structure of cloned FSV-1 and FSV-2 DNA molecules by restriction endonuclease mapping and hybridization with appropriate probes showed that about 2.6 kb of the FSV-unique sequence called FSV-fps is located in the middle of the FSV genome and is flanked by helper virus-derived sequences of about 1.3 kb at the 5′ end and 0.5 kb at the 3′ end. The long terminal repeats of FSV were found to have no cleavage site for either EcoRI or PvuI. Upon transfection, both FSV-1 DNA and FSV-2 DNA were able to transform mammalian fibroblasts. Four ³²P-labeled DNA fragments derived from different portions of the FSV-fps sequence were used for hybridization to viral RNAs. We found that sequences within the 3′ half of the FSV-fps gene are homologous to RNAs of PRCII avian sarcoma virus and the Snyder-Theilen strain of feline sarcoma virus, both of which were previously shown to contain transforming genes related to FSV-fps. These results suggest that the 3′ portion of the FSV-fps sequence may be crucial for the transforming activity of fps-related oncogenic sequences.