Contig Assembly Research Articles

Random Amplification and Pyrosequencing for Identification of Novel Viral Genome Sequences

ssRNA viruses have high levels of genomic divergence, which can lead to difficulty in genomic characterization of new viruses using traditional PCR amplification and sequencing methods. In this study, random reverse transcription, anchored random PCR amplification, and high-throughput pyrosequencing were used to identify orthobunyavirus sequences from total RNA extracted from viral cultures of acute febrile illness specimens. Draft genome sequence for the orthobunyavirus L segment was assembled and sequentially extended using de novo assembly contigs from pyrosequencing reads and orthobunyavirus sequences in GenBank as guidance. Accuracy and continuous coverage were achieved by mapping all reads to the L segment draft sequence. Subsequently, RT-PCR and Sanger sequencing were used to complete the genome sequence. The complete L segment was found to be 6936 bases in length, encoding a 2248-aa putative RNA polymerase. The identified L segment was distinct from previously published South American orthobunyaviruses, sharing 63% and 54% identity at the nucleotide and amino acid level, respectively, with the complete Oropouche virus L segment and 73% and 81% identity at the nucleotide and amino acid level, respectively, with a partial Caraparu virus L segment. The result demonstrated the effectiveness of a sequence-independent amplification and next-generation sequencing approach for obtaining complete viral genomes from total nucleic acid extracts and its use in pathogen discovery.

Draft Genome of the Pearl Oyster Pinctada fucata: A Platform for Understanding Bivalve Biology

The study of the pearl oyster Pinctada fucata is key to increasing our understanding of the molecular mechanisms involved in pearl biosynthesis and biology of bivalve molluscs. We sequenced ∼1150-Mb genome at ∼40-fold coverage using the Roche 454 GS-FLX and Illumina GAIIx sequencers. The sequences were assembled into contigs with N50 = 1.6 kb (total contig assembly reached to 1024 Mb) and scaffolds with N50 = 14.5 kb. The pearl oyster genome is AT-rich, with a GC content of 34%. DNA transposons, retrotransposons, and tandem repeat elements occupied 0.4, 1.5, and 7.9% of the genome, respectively (a total of 9.8%). Version 1.0 of the P. fucata draft genome contains 23 257 complete gene models, 70% of which are supported by the corresponding expressed sequence tags. The genes include those reported to have an association with bio-mineralization. Genes encoding transcription factors and signal transduction molecules are present in numbers comparable with genomes of other metazoans. Genome-wide molecular phylogeny suggests that the lophotrochozoan represents a distinct clade from ecdysozoans. Our draft genome of the pearl oyster thus provides a platform for the identification of selection markers and genes for calcification, knowledge of which will be important in the pearl industry.

Next Generation Transcriptome Sequencing and Quantitative Real-Time PCR Technologies for Characterisation of the Bemisia tabaci Asia 1 mtCOI Phylogenetic Clade

A programme of functional genomics research is underway at the University of Greenwich, UK, to develop and apply genomics technologies to characterise an economically-important but under-researched Bemisia tabaci (Hemiptera: Aleyrodidae), the Asia 1 mtCOI phylogenetic group. A comparison of this putative species from India with other important B. tabaci populations and insect species may provide targets for the development of more effective whitefly control strategies. As a first step, next-generation sequencing (NGS) has been used to survey the transcriptome of adult female whitefly, with high quality RNA preparations being used to generate cDNA libraries for NGS using the Roche 454 Titanium DNA sequencing platform. Contig assemblies constructed from the resultant sequences (301 094 reads) using the software program CLC Genomics Workbench generated 3 821 core contigs. Comparison of a selection of these contigs with related sequences from other B. tabaci genetic groups has revealed good alignment for some genes (e.g., HSP90) but misassemblies in other datasets (e.g., the vitellogenin gene family), highlighting the need for manual curation as well as collaborative international efforts to obtain accurate assemblies from the existing next generation sequence datasets. Nevertheless, data emerging from the NGS has facilitated the development of accurate and reliable methods for analysing gene expression based on quantitative real-time RT-PCR, illustrating the power of this approach to enable rapid expression analyses in an organism for which a complete genome sequence is currently lacking.

Transcriptome Analysis of Silver Carp (Hypophthalmichthys molitrix) by Paired-End RNA Sequencing

The silver carp (Hypophthalmichthys molitrix) is among the most intensively pond-cultured fish species and is used in the wild to counteract water bloom in China. However, little genomic information is available for this species, especially regarding its ability to grow rapidly in water, even water contaminated with high concentrations of poisonous microcystin. In this study, we performed de novo transcriptome assembly and analysis of the 17.10 million short-read sequences produced by the Illumina paired-end sequencing technology. Using an improved multiple k-mer contig assembly method coupled with further scaffolding, 85 759 sequences were obtained. There were 23 044 sequences annotated with 3423 gene ontology terms for 104 196 term occurrences and the three corresponding organizing principles. A total of 38 200 assembled sequences were involved in 218 predicted Kyoto Encyclopedia of Genes and Genomes metabolic pathways. We also recovered 41 of 44 genes involved in the biosynthesis of glutathione. Of these, five genes were identified as experienced positive selection between silver carp and zebrafish, as determined by the likelihood ratio test. This report is the first annotated review of the silver carp transcriptome. These data will be of interest to researchers investigating the evolution and biological processes of the silver carp. This work also provides an archive for future studies of recent speciation and evolution of Cyprinidae fishes and can be used in comparative studies of other fishes.

Discovery of DNA viruses in wild-caught mosquitoes using small RNA high throughput sequencing.

BackgroundMosquito-borne infectious diseases pose a severe threat to public health in many areas of the world. Current methods for pathogen detection and surveillance are usually dependent on prior knowledge of the etiologic agents involved. Hence, efficient approaches are required for screening wild mosquito populations to detect known and unknown pathogens.Methodology/principal findingsIn this study, we explored the use of Next Generation Sequencing to identify viral agents in wild-caught mosquitoes. We extracted total RNA from different mosquito species from South China. Small 18–30 bp length RNA molecules were purified, reverse-transcribed into cDNA and sequenced using Illumina GAIIx instrumentation. Bioinformatic analyses to identify putative viral agents were conducted and the results confirmed by PCR. We identified a non-enveloped single-stranded DNA densovirus in the wild-caught Culex pipiens molestus mosquitoes. The majority of the viral transcripts (.>80% of the region) were covered by the small viral RNAs, with a few peaks of very high coverage obtained. The +/− strand sequence ratio of the small RNAs was approximately 7∶1, indicating that the molecules were mainly derived from the viral RNA transcripts. The small viral RNAs overlapped, enabling contig assembly of the viral genome sequence. We identified some small RNAs in the reverse repeat regions of the viral 5′- and 3′ -untranslated regions where no transcripts were expected.Conclusions/significanceOur results demonstrate for the first time that high throughput sequencing of small RNA is feasible for identifying viral agents in wild-caught mosquitoes. Our results show that it is possible to detect DNA viruses by sequencing the small RNAs obtained from insects, although the underlying mechanism of small viral RNA biogenesis is unclear. Our data and those of other researchers show that high throughput small RNA sequencing can be used for pathogen surveillance in wild mosquito vectors.

Rapid whole genome optical mapping of Plasmodium falciparum

BackgroundImmune evasion and drug resistance in malaria have been linked to chromosomal recombination and gene copy number variation (CNV). These events are ideally studied using comparative genomic analyses; however in malaria these analyses are not as common or thorough as in other infectious diseases, partly due to the difficulty in sequencing and assembling complete genome drafts. Recently, whole genome optical mapping has gained wide use in support of genomic sequence assembly and comparison. Here, a rapid technique for producing whole genome optical maps of Plasmodium falciparum is described and the results of mapping four genomes are presented.MethodsFour laboratory strains of P. falciparum were analysed using the Argus™ optical mapping system to produce ordered restriction fragment maps of all 14 chromosomes in each genome. Plasmodium falciparum DNA was isolated directly from blood culture, visualized using the Argus™ system and assembled in a manner analogous to next generation sequence assembly into maps (AssemblyViewer™, OpGen Inc.®). Full coverage maps were generated for P. falciparum strains 3D7, FVO, D6 and C235. A reference P. falciparum in silico map was created by the digestion of the genomic sequence of P. falciparum with the restriction enzyme AflII, for comparisons to genomic optical maps. Maps were then compared using the MapSolver™ software.ResultsGenomic variation was observed among the mapped strains, as well as between the map of the reference strain and the map derived from the putative sequence of that same strain. Duplications, deletions, insertions, inversions and misassemblies of sizes ranging from 3,500 base pairs up to 78,000 base pairs were observed. Many genomic events occurred in areas of known repetitive sequence or high copy number genes, including var gene clusters and rifin complexes.ConclusionsThis technique for optical mapping of multiple malaria genomes allows for whole genome comparison of multiple strains and can assist in identifying genetic variation and sequence contig assembly. New protocols and technology allowed us to produce high quality contigs spanning four P. falciparum genomes in six weeks for less than $1,000.00 per genome. This relatively low cost and quick turnaround makes the technique valuable compared to other genomic sequencing technologies for studying genetic variation in malaria.

A multi-organ transcriptome resource for the Burmese Python (Python molurus bivittatus)

BackgroundSnakes provide a unique vertebrate system for studying a diversity of extreme adaptations, including those related to development, metabolism, physiology, and venom. Despite their importance as research models, genomic resources for snakes are few. Among snakes, the Burmese python is the premier model for studying extremes of metabolic fluctuation and physiological remodelling. In this species, the consumption of large infrequent meals can induce a 40-fold increase in metabolic rate and more than a doubling in size of some organs. To provide a foundation for research utilizing the python, our aim was to assemble and annotate a transcriptome reference from the heart and liver. To accomplish this aim, we used the 454-FLX sequencing platform to collect sequence data from multiple cDNA libraries.ResultsWe collected nearly 1 million 454 sequence reads, and assembled these into 37,245 contigs with a combined length of 13,409,006 bp. To identify known genes, these contigs were compared to chicken and lizard gene sets, and to all Genbank sequences. A total of 13,286 of these contigs were annotated based on similarity to known genes or Genbank sequences. We used gene ontology (GO) assignments to characterize the types of genes in this transcriptome resource. The raw data, transcript contig assembly, and transcript annotations are made available online for use by the broader research community.ConclusionThese data should facilitate future studies using pythons and snakes in general, helping to further contribute to the utilization of snakes as a model evolutionary and physiological system. This sequence collection represents a major genomic resource for the Burmese python, and the large number of transcript sequences characterized should contribute to future research in this and other snake species.

Transcriptome sequencing, and rapid development and application of SNP markers for the legume pod borer Maruca vitrata (Lepidoptera: Crambidae).

The legume pod borer, Maruca vitrata (Lepidoptera: Crambidae), is an insect pest species of crops grown by subsistence farmers in tropical regions of Africa. We present the de novo assembly of 3729 contigs from 454- and Sanger-derived sequencing reads for midgut, salivary, and whole adult tissues of this non-model species. Functional annotation predicted that 1320 M. vitrata protein coding genes are present, of which 631 have orthologs within the Bombyx mori gene model. A homology-based analysis assigned M. vitrata genes into a group of paralogs, but these were subsequently partitioned into putative orthologs following phylogenetic analyses. Following sequence quality filtering, a total of 1542 putative single nucleotide polymorphisms (SNPs) were predicted within M. vitrata contig assemblies. Seventy one of 1078 designed molecular genetic markers were used to screen M. vitrata samples from five collection sites in West Africa. Population substructure may be present with significant implications in the insect resistance management recommendations pertaining to the release of biological control agents or transgenic cowpea that express Bacillus thuringiensis crystal toxins. Mutation data derived from transcriptome sequencing is an expeditious and economical source for genetic markers that allow evaluation of ecological differentiation.

Evaluation of High-Throughput Sequencing for Identifying Known and Unknown Viruses in Biological Samples

High-throughput sequencing furnishes a large number of short sequence reads from uncloned DNA and has rapidly become a major tool for identifying viruses in biological samples, and in particular when the target sequence is undefined. In this study, we assessed the analytical sensitivity of a pipeline for detection of viruses in biological samples based on either the Roche-454 genome sequencer or Illumina genome analyzer platforms. We sequenced biological samples artificially spiked with a wide range of viruses with genomes composed of single or double-stranded DNA or RNA, including linear or circular single-stranded DNA. Viruses were added at a very low concentration most often corresponding to 3 or 0.8 times the validated level of detection of quantitative reverse transcriptase PCRs (RT-PCRs). For the viruses represented, or resembling those represented, in public nucleotide sequence databases, we show that the higher output of Illumina is associated with a much greater sensitivity, approaching that of optimized quantitative (RT-)PCRs. In this blind study, identification of viruses was achieved without incorrect identification. Nevertheless, at these low concentrations, the number of reads generated by the Illumina platform was too small to facilitate assembly of contigs without the use of a reference sequence, thus precluding detection of unknown viruses. When the virus load was sufficiently high, de novo assembly permitted the generation of long contigs corresponding to nearly full-length genomes and thus should facilitate the identification of novel viruses.

De novo assembly and transcriptome analysis of five major tissues of Jatropha curcas L. using GS FLX titanium platform of 454 pyrosequencing

BackgroundJatropha curcas L. is an important non-edible oilseed crop with promising future in biodiesel production. However, factors like oil yield, oil composition, toxic compounds in oil cake, pests and diseases limit its commercial potential. Well established genetic engineering methods using cloned genes could be used to address these limitations. Earlier, 10,983 unigenes from Sanger sequencing of ESTs, and 3,484 unique assembled transcripts from 454 pyrosequencing of uncloned cDNAs were reported. In order to expedite the process of gene discovery, we have undertaken 454 pyrosequencing of normalized cDNAs prepared from roots, mature leaves, flowers, developing seeds, and embryos of J. curcas.ResultsFrom 383,918 raw reads, we obtained 381,957 quality-filtered and trimmed reads that are suitable for the assembly of transcript sequences. De novo contig assembly of these reads generated 17,457 assembled transcripts (contigs) and 54,002 singletons. Average length of the assembled transcripts was 916 bp. About 30% of the transcripts were longer than 1000 bases, and the size of the longest transcript was 7,173 bases. BLASTX analysis revealed that 2,589 of these transcripts are full-length. The assembled transcripts were validated by RT-PCR analysis of 28 transcripts. The results showed that the transcripts were correctly assembled and represent actively expressed genes. KEGG pathway mapping showed that 2,320 transcripts are related to major biochemical pathways including the oil biosynthesis pathway. Overall, the current study reports 14,327 new assembled transcripts which included 2589 full-length transcripts and 27 transcripts that are directly involved in oil biosynthesis.ConclusionThe large number of transcripts reported in the current study together with existing ESTs and transcript sequences will serve as an invaluable genetic resource for crop improvement in jatropha. Sequence information of those genes that are involved in oil biosynthesis could be used for metabolic engineering of jatropha to increase oil content, and to modify oil composition.

Local De Novo Assembly of RAD Paired-End Contigs Using Short Sequencing Reads

Despite the power of massively parallel sequencing platforms, a drawback is the short length of the sequence reads produced. We demonstrate that short reads can be locally assembled into longer contigs using paired-end sequencing of restriction-site associated DNA (RAD-PE) fragments. We use this RAD-PE contig approach to identify single nucleotide polymorphisms (SNPs) and determine haplotype structure in threespine stickleback and to sequence E. coli and stickleback genomic DNA with overlapping contigs of several hundred nucleotides. We also demonstrate that adding a circularization step allows the local assembly of contigs up to 5 kilobases (kb) in length. The ease of assembly and accuracy of the individual contigs produced from each RAD site sequence suggests RAD-PE sequencing is a useful way to convert genome-wide short reads into individually-assembled sequences hundreds or thousands of nucleotides long.

Cloning, sequencing and functional expression of cytosolic malate dehydrogenase from Taenia solium: Purification and characterization of the recombinant enzyme

We report herein the complete coding sequence of a Taenia solium cytosolic malate dehydrogenase (TscMDH). The cDNA fragment, identified from the T. solium genome project database, encodes a protein of 332 amino acid residues with an estimated molecular weight of 36517Da. For recombinant expression, the full length coding sequence was cloned into pET23a. After successful expression and enzyme purification, isoelectrofocusing gel electrophoresis allowed to confirm the calculated pI value at 8.1, as deduced from the amino acid sequence. The recombinant protein (r-TscMDH) showed MDH activity of 409U/mg in the reduction of oxaloacetate, with neither lactate dehydrogenase activity nor NADPH selectivity. Optimum pH for enzyme activity was 7.6 for oxaloacetate reduction and 9.6 for malate oxidation. Kcat values for oxaloacetate, malate, NAD, and NADH were 665, 47, 385, and 962s−1, respectively. Additionally, a partial characterization of TsMDH gene structure after analysis of a 1.56Kb genomic contig assembly is also reported.

Sequence-based physical mapping of complex genomes by whole genome profiling

We present whole genome profiling (WGP), a novel next-generation sequencing-based physical mapping technology for construction of bacterial artificial chromosome (BAC) contigs of complex genomes, using Arabidopsis thaliana as an example. WGP leverages short read sequences derived from restriction fragments of two-dimensionally pooled BAC clones to generate sequence tags. These sequence tags are assigned to individual BAC clones, followed by assembly of BAC contigs based on shared regions containing identical sequence tags. Following in silico analysis of WGP sequence tags and simulation of a map of Arabidopsis chromosome 4 and maize, a WGP map of Arabidopsis thaliana ecotype Columbia was constructed de novo using a six-genome equivalent BAC library. Validation of the WGP map using the Columbia reference sequence confirmed that 350 BAC contigs (98%) were assembled correctly, spanning 97% of the 102-Mb calculated genome coverage. We demonstrate that WGP maps can also be generated for more complex plant genomes and will serve as excellent scaffolds to anchor genetic linkage maps and integrate whole genome sequence data.

LTC: a novel algorithm to improve the efficiency of contig assembly for physical mapping in complex genomes

BackgroundPhysical maps are the substrate of genome sequencing and map-based cloning and their construction relies on the accurate assembly of BAC clones into large contigs that are then anchored to genetic maps with molecular markers. High Information Content Fingerprinting has become the method of choice for large and repetitive genomes such as those of maize, barley, and wheat. However, the high level of repeated DNA present in these genomes requires the application of very stringent criteria to ensure a reliable assembly with the FingerPrinted Contig (FPC) software, which often results in short contig lengths (of 3-5 clones before merging) as well as an unreliable assembly in some difficult regions. Difficulties can originate from a non-linear topological structure of clone overlaps, low power of clone ordering algorithms, and the absence of tools to identify sources of gaps in Minimal Tiling Paths (MTPs).ResultsTo address these problems, we propose a novel approach that: (i) reduces the rate of false connections and Q-clones by using a new cutoff calculation method; (ii) obtains reliable clusters robust to the exclusion of single clone or clone overlap; (iii) explores the topological contig structure by considering contigs as networks of clones connected by significant overlaps; (iv) performs iterative clone clustering combined with ordering and order verification using re-sampling methods; and (v) uses global optimization methods for clone ordering and Band Map construction. The elements of this new analytical framework called Linear Topological Contig (LTC) were applied on datasets used previously for the construction of the physical map of wheat chromosome 3B with FPC. The performance of LTC vs. FPC was compared also on the simulated BAC libraries based on the known genome sequences for chromosome 1 of rice and chromosome 1 of maize.ConclusionsThe results show that compared to other methods, LTC enables the construction of highly reliable and longer contigs (5-12 clones before merging), the detection of "weak" connections in contigs and their "repair", and the elongation of contigs obtained by other assembly methods.

Integration of Genetic, Physical, and Cytogenetic Maps for Brassica rapa Chromosome A7

Bacterial artificial chromosome (BAC) contigs have been genetically mapped to the 10 linkage groups of Brassica rapa by BAC end sequences (BES). To integrate the genetic, physical, and cytogenetic maps, fluorescence in situ hybridization (FISH) was used to anchor the assembly of BAC contigs onto Brassica chromosomes using representative BACs. This BAC-FISH approach can be used to identify chromosome arms on separate mitotic metaphase chromosomes or to map multiple BACs to single long pachytene chromosomes. As part of an international consortium that is sequencing the B. rapa genome, we integrated the linkage and physical maps with the B. rapa cytogenetic map for chromosome A7 by hybridizing BACs to mitotic chromosomes and along the length of pachytene chromosome spreads. A total of 31 BACs that were putatively located on A7 were used as probes for FISH analyses; however, only 19 BACs mapped unambiguously to A7 while the remaining BACs either mapped to other chromosomes or hybridized to multiple locations. We then created a multicolor FISH cocktail of 16 BAC probes to simultaneously hybridize the entire length of the A7 chromosome. We successfully applied the 16 A7 BAC probe mix to B. rapa, B. oleracea, and domesticated and resynthesized genotypes of B. napus to demonstrate that this approach can facilitate studies of genome evolution by integrating the genetic, physical, and cytogenetic maps among closely related species of Brassica.

Generation of a BAC-based physical map of the melon genome.

BackgroundCucumis melo (melon) belongs to the Cucurbitaceae family, whose economic importance among horticulture crops is second only to Solanaceae. Melon has high intra-specific genetic variation, morphologic diversity and a small genome size (450 Mb), which make this species suitable for a great variety of molecular and genetic studies that can lead to the development of tools for breeding varieties of the species. A number of genetic and genomic resources have already been developed, such as several genetic maps and BAC genomic libraries. These tools are essential for the construction of a physical map, a valuable resource for map-based cloning, comparative genomics and assembly of whole genome sequencing data. However, no physical map of any Cucurbitaceae has yet been developed. A project has recently been started to sequence the complete melon genome following a whole-genome shotgun strategy, which makes use of massive sequencing data. A BAC-based melon physical map will be a useful tool to help assemble and refine the draft genome data that is being produced.ResultsA melon physical map was constructed using a 5.7 × BAC library and a genetic map previously developed in our laboratories. High-information-content fingerprinting (HICF) was carried out on 23,040 BAC clones, digesting with five restriction enzymes and SNaPshot labeling, followed by contig assembly with FPC software. The physical map has 1,355 contigs and 441 singletons, with an estimated physical length of 407 Mb (0.9 × coverage of the genome) and the longest contig being 3.2 Mb. The anchoring of 845 BAC clones to 178 genetic markers (100 RFLPs, 76 SNPs and 2 SSRs) also allowed the genetic positioning of 183 physical map contigs/singletons, representing 55 Mb (12%) of the melon genome, to individual chromosomal loci. The melon FPC database is available for download at http://melonomics.upv.es/static/files/public/physical_map/.ConclusionsHere we report the construction of the first physical map of a Cucurbitaceae species described so far. The physical map was integrated with the genetic map so that a number of physical contigs, representing 12% of the melon genome, could be anchored to known genetic positions. The data presented is already helping to improve the quality of the melon genomic sequence available as a result of a project currently being carried out in Spain, adopting a whole genome shotgun approach based on 454 sequencing data.

Feasibility of physical map construction from fingerprinted bacterial artificial chromosome libraries of polyploid plant species

BackgroundThe presence of closely related genomes in polyploid species makes the assembly of total genomic sequence from shotgun sequence reads produced by the current sequencing platforms exceedingly difficult, if not impossible. Genomes of polyploid species could be sequenced following the ordered-clone sequencing approach employing contigs of bacterial artificial chromosome (BAC) clones and BAC-based physical maps. Although BAC contigs can currently be constructed for virtually any diploid organism with the SNaPshot high-information-content-fingerprinting (HICF) technology, it is currently unknown if this is also true for polyploid species. It is possible that BAC clones from orthologous regions of homoeologous chromosomes would share numerous restriction fragments and be therefore included into common contigs. Because of this and other concerns, physical mapping utilizing the SNaPshot HICF of BAC libraries of polyploid species has not been pursued and the possibility of doing so has not been assessed. The sole exception has been in common wheat, an allohexaploid in which it is possible to construct single-chromosome or single-chromosome-arm BAC libraries from DNA of flow-sorted chromosomes and bypass the obstacles created by polyploidy.ResultsThe potential of the SNaPshot HICF technology for physical mapping of polyploid plants utilizing global BAC libraries was evaluated by assembling contigs of fingerprinted clones in an in silico merged BAC library composed of single-chromosome libraries of two wheat homoeologous chromosome arms, 3AS and 3DS, and complete chromosome 3B. Because the chromosome arm origin of each clone was known, it was possible to estimate the fidelity of contig assembly. On average 97.78% or more clones, depending on the library, were from a single chromosome arm. A large portion of the remaining clones was shown to be library contamination from other chromosomes, a feature that is unavoidable during the construction of single-chromosome BAC libraries.ConclusionsThe negligibly low level of incorporation of clones from homoeologous chromosome arms into a contig during contig assembly suggested that it is feasible to construct contigs and physical maps using global BAC libraries of wheat and almost certainly also of other plant polyploid species with genome sizes comparable to that of wheat. Because of the high purity of the resulting assembled contigs, they can be directly used for genome sequencing. It is currently unknown but possible that equally good BAC contigs can be also constructed for polyploid species containing smaller, more gene-rich genomes.

Isolation, Sequence Analysis, and Linkage Mapping of Nucleotide Binding Site–Leucine-rich Repeat Disease Resistance Gene Analogs in Watermelon

Sixty-six watermelon (Citrullus lanatus var. lanatus) disease resistance gene analogs were cloned from ‘Calhoun Gray’, PI 296341, and PI 595203 using degenerate primers to select for the nucleotide binding sites (NBS) from the NBS–leucine-rich repeat (LRR) resistance gene family. After contig assembly, watermelon resistance gene analogs (WRGA) were identified and amino acid sequence alignment revealed that these groups contained motifs characteristic of NBS-LRR resistance genes. Using cluster analysis, eight groups of WRGA were identified and further characterized as having homology to Drosophila Toll and mammalian interleukin-1 receptor (TIR) and non-TIR domains. Three of these WRGA as well as three disease-related watermelon expressed sequence tag homologs were placed on a test-cross map. Linkage mapping placed the WRGA on linkage group XIII, an area on the watermelon map where resistance gene analogs cluster. In addition, these WRGA sequence-tagged sites (STS) were amplified from various genera of the Cucurbitaceae indicating that conservation of resistance gene analogs exists among cucurbits. These WRGA-STS markers may be useful in marker-assisted selection for the improvement for disease resistance in watermelon.

Cloning of TaCYP707A1 Gene that Encodes ABA 8′-Hydroxylase in Common Wheat ( Triticum aestivum L.)

The plant hormone abscisic acid (ABA) regulates many important physiological and developmental processes in plants. The objective of this study was to clone the ABA 8′-hydroxylase gene in common wheat. In the present study, we used the cDNA sequence of barley HvCYP707A1 gene (GenBank accession no. AB239299) as a probe for BLAST search against the common wheat ( Triticum aestivum L.) EST database in GenBank. All wheat ESTs sharing high similarity with the reference gene were subjected to contig assembly. Primers were designed based on the constructed contigs to clone the wheat CYP707A1 gene, designated as TaCYP707A1. The genomic DNA sequence of TaCYP707A1 gene comprised five exons and four introns, with a size of 2 225 bp. The corresponding cDNA sequence of TaCYP707A1 was 1 737 bp, containing an open reading frame (ORF) of 1 431 bp, a 42-bp 5′-untranslated region (UTR) and a 264-bp 3′UTR, with 94.9% of identical sequences to HvCYP707A1 gene ( AB239299). The neighbor joining tree indicated that the deduced amino acid sequences of TaCYP707A1 gene was highly similar to those of barley and rice. The TaCYP707A1 gene was located on chromosome 6BL using a set of Chinese Spring nullisomic-tetrasomic lines and ditelosomic line 6BS. These results will be of high importance in understanding of molecular mechanism of ABA catabolism.

Role of Fluorescence in situ Hybridization in Sequencing the Tomato Genome

The tomato (Solanum lycopersicum L.) genome is being sequenced by a consortium of laboratories in 10 countries. Seventy-seven percent of the tomato genome (DNA) is located in repeat-rich, gene-poor, pericentric heterochromatin, while 23% of the genome is located in repeat-poor, gene-rich, distal euchromatin. It is estimated that approximately 90% of tomato’s nuclear genes can be characterized by limiting the sequencing effort to euchromatin while avoiding the problems involved in sequencing the repetitive DNA in heterochromatin. Sequencing is being performed on tomato nuclear DNA cloned into bacterial artificial chromosome (BAC) vectors. Fluorescence in situ hybridization (FISH) is used to help direct the sequencing effort by cytologically demonstrating the location of selected BACs on tomato chromosomes. While mitotic metaphase chromosomes are too short and compact for this purpose, long pachytene chromosomes are ideal. BACs localized in euchromatin can be used confidently as anchors for the assembly of BAC contigs that extend through the euchromatic length of each chromosome arm. Another important role for FISH is identification of BACs near telomeres and near borders with pericentric heterochromatin to indicate that sequencing should not extend much further. This role of FISH is enhanced by our ability to estimate base pair distances between localized BACs and these chromosomal features. Finally, it is noteworthy that when BAC-FISH is combined with chromosomal in situ suppression (CISS) hybridization to block repeats and localize single/low copy sequences, the great majority of BACs localize to single sites. This observation is consistent with tomato being an ancient diploid.