Filarial and Wolbachia genomics.

Filarial nematode parasites, the causative agents for a spectrum of acute and chronic diseases including lymphatic filariasis and river blindness, threaten the well-being and livelihood of hundreds of millions of people in the developing regions of the world. The 2007 publication on a draft assembly of the 95-Mb genome of the human filarial parasite Brugia malayi- representing the first helminth parasite genome to be sequenced - has been followed in rapid succession by projects that have resulted in the genome sequencing of six additional filarial species, seven nonfilarial nematode parasites of animals and nearly 30 plant parasitic and free-living species. Parallel to the genomic sequencing, transcriptomic and proteomic projects have facilitated genome annotation, expanded our understanding of stage-associated gene expression and provided a first look at the role of epigenetic regulation of filarial genomes through microRNAs. The expansion in filarial genomics will also provide a significant enrichment in our knowledge of the diversity and variability in the genomes of the endosymbiotic bacterium Wolbachia leading to a better understanding of the genetic principles that govern filarial-Wolbachia mutualism. The goal here is to provide an overview of the trends and advances in filarial and Wolbachia genomics.

Determination of Wolbachia genome size by pulsed-field gel electrophoresis.

Genome sizes of six different Wolbachia strains from insect and nematode hosts have been determined by pulsed-field gel electrophoresis of purified DNA both before and after digestion with rare-cutting restriction endonucleases. Enzymes SmaI, ApaI, AscI, and FseI cleaved the studied Wolbachia strains at a small number of sites and were used for the determination of the genome sizes of wMelPop, wMel, and wMelCS (each 1.36 Mb), wRi (1.66 Mb), wBma (1.1 Mb), and wDim (0.95 Mb). The Wolbachia genomes studied were all much smaller than the genomes of free-living bacteria such as Escherichia coli (4.7 Mb), as is typical for obligate intracellular bacteria. There was considerable genome size variability among Wolbachia strains, especially between the more parasitic A group Wolbachia infections of insects and the mutualistic C and D group infections of nematodes. The studies described here found no evidence for extrachromosomal plasmid DNA in any of the strains examined. They also indicated that the Wolbachia genome is circular.

Denaturing polyacrylamide gel electrophoresis.

Polyacrylamide gels that contain a high concentration of urea as a denaturant are capable of resolving short (<500 nucleotides) single-stranded fragments of DNA or RNA that differ in length by as little as one nucleotide. Such gels are commonly used for DNA sequence analysis, as well as in PCR amplification of SSLPs (simple sequence length polymorphisms) for genotyping, genotyping by the ligase chain reaction (LCR), analysis of mutations by RNase A cleavage, chemical cleavage of heteroduplex DNA to identify mutations, and molecular analysis of fragile X syndrome by PCR. This appendix presents a protocol for the pouring, running, and processing of a typical gel which is 40-cm long with a uniform thickness of 0.4 mm, containing 7 M urea and 4% to 8% acrylamide.

Denaturing polyacrylamide gel electrophoresis.

Thin polyacrylamide gels that contain a high concentration of urea as a denaturant are capable of resolving short (<500 nucleotides) single-stranded fragments of DNA or RNA that differ in length by as little as one nucleotide. Such gels are uniquely suited for nucleic acid sequence analysis, which is required, for instance, for all footprinting protocols. Thicker gels are often used to purify oligonucleotides. This appendix describes the pouring, running, and processing of a typical sequencing gel, which is 40 cm long with a uniform thickness of 0.4 mm, containing 7 M urea and 4% to 8% acrylamide.

Constructing nested deletions for use in DNA sequencing.

Nested deletions useful for dideoxy DNA sequencing are a set of deletions originating at one end of a target DNA fragment and extending various lengths along the target DNA. Each successively longer deletion brings "new" regions of the target DNA into sequencing range (about 300 bp for normal sequencing gels) of the primer site for a general discussion of nested deletions in DNA sequencing). Two protocols for generating nested subclones via enzymatic digestion are included in this unit. In the first, a set of nested deletions is generated by exonuclease III. The primary advantage of this method is that the deletion products generated from the original clone can be recircularized to generate functional plasmids and thus do not require subcloning into another vector. An alternate method utilizes Bal 31 nuclease to generate the deletions. This method requires subcloning of the deletion fragments into a separate vector for subsequent use. Both methods require the presence of unique restriction sites in the vector that are not present in the insert DNA. Bal 31 can also be used to generate nested deletions for chemical sequencing in conjunction with specialized chemical sequencing vectors.

Preparation of templates for DNA sequencing.

This unit contains protocols for preparing DNA suitable for use as dideoxy sequencing templates and as material for end labeling and chemical sequencing. In all protocols, the starting material contains the recombinant molecule to be sequenced. DNA from M13mp-derived phage is easily prepared and is currently the most reliable source of template for large-scale dideoxy sequencing projects. Because it is occasionally necessary or convenient to use a lambda-derived phage as a source of DNA, a protocol for preparing lambda phage DNA from plate lysates is provided. Two protocols for minipreps of plasmid DNA are provided, one intended for dideoxy sequencing, the other for end labeling and chemical sequencing; they differ primarily in the way in which cellular RNA is removed. Alkali denaturation of double-stranded DNA (necessary prior to annealing) is described, and a final protocol describes the preparation of template for thermal cycle sequencing from a single phage plaque or bacterial colony.

DNA sequencing by the dideoxy method.

In the basic dideoxy sequencing reaction, an oligonucleotide primer is annealed to a single-stranded DNA template and extended by DNA polymerase in the presence of four deoxyribonucleoside triphosphates (dNTPs), one of which is 35S-labeled. The reaction also contains one of four dideoxyribonucleoside triphosphates (ddNTPs), which terminate elongation when incorporated into the growing DNA chain. After completion of the sequencing reactions, the products are subjected to electrophoresis on a high-resolution denaturing polyacrylamide gel and then autoradiographed to visualize the DNA sequence. Three variations of the dideoxy sequencing procedure are currently in use and are presented in this unit. In the "labeling/termination" procedure, primer chains are initially extended and labeled in the absence of terminating ddNTPs, whereas in the traditional "Sanger" procedure, labeling and termination of primer chains occur in a single step. A recent variation of the dideoxy sequencing method is thermal cycle sequencing in which the reaction mixture, containing template DNA, primer, thermostable DNA polymerase, dNTPs, and ddNTPs, is subjected to repeated rounds of denaturation, annealing, and elongation steps. The resulting linear amplification of the sequencing products allows much less template DNA to be used and eliminates independent primer annealing and template denaturation steps, which are required for the labeling/termination or Sanger procedures. The use of automated fluorescent sequencers for four-color dideoxy DNA sequencing is also described in detail.

Denaturing gel electrophoresis for sequencing.

The accuracy of DNA sequence determination depends largely upon resolution of the sequencing products in denaturing polyacrylamide gels. This unit provides a detailed description of the setup, electrophoresis, and processing of such gels. In general, the gels required for DNA sequencing are 40-cm long, of uniform thickness, and contain 4% to 8% acrylamide and 7 M urea. Modifications of this protocol increase the length of readable sequence information which can be obtained from a single gel (i.e., forming the gel with wedge-shaped spacers to create a field gradient, or incorporating a buffer gradient, an electrolyte gradient, or an acrylamide step gradient into the gel). A modification to the Basic Protocol--inclusion of formamide in the sequencing gel--is designed to overcome gel compressions arising from secondary structure in the sequencing products during gel electrophoresis. A discussion of acrylamide concentrations and electrophoresis conditions is included in the Commentary.

The filarial genome project: analysis of the nuclear, mitochondrial and endosymbiont genomes of Brugia malayi.

The Filarial Genome Project (FGP) was initiated in 1994 under the auspices of the World Health Organisation. Brugia malayi was chosen as the model organism due to the availability of all life cycle stages for the construction of cDNA libraries. To date, over 20000 cDNA clones have been partially sequenced and submitted to the EST database (dbEST). These ESTs define approximately 7000 new Brugia genes. Analysis of the EST dataset provides useful information on the expression pattern of the most abundantly expressed Brugia genes. Some highly expressed genes have been identified that are expressed in all stages of the parasite's life cycle, while other highly expressed genes appear to be stage-specific. To elucidate the structure of the Brugia genome and to provide a basis for comparison to the Caenorhabditis elegans genome, the FGP is also constructing a physical map of the Brugia chromosomes and is sequencing genomic BAC clones. In addition to the nuclear genome, B. malayi possesses two other genomes: the mitochondrial genome and the genome of a bacterial endosymbiont. Eighty percent of the mitochondrial genome of B. malayi has been sequenced and is being compared to mitochondrial sequences of other nematodes. The bacterial endosymbiont genome found in B. malayi is closely related to the Wolbachia group of rickettsia-like bacteria that infects many insect species. A set of overlapping BAC clones is being assembled to cover the entire bacterial genome. Currently, half of the bacterial genome has been assembled into four contigs. A consortium has been established to sequence the entire genome of the Brugia endosymbiont. The sequence and mapping data provided by the FGP is being utilised by the nematode research community to develop a better understanding of the biology of filarial parasites and to identify new vaccine candidates and drug targets to aid the elimination of human filariasis.

Implementation of automation in a small-scale DNA sequencing core facility.

New England Biolabs (NEB) sequencing core facility provides automated sequencing services to support various company-wide projects in house, but on a very small scale of about 1000 to 1500 reactions per month. A procedure has been implemented at the NEB core sequencing facility to integrate simplified methods and robotics to provide a more efficient small-scale process. This has been done using a Beckman Biomek 2000 robot combined with an MJ DNA Engine, 96-well plate cycler (PTC-200), AB 373 and 377 sequencers, BMA Singel gels, and several other materials that help reduce the time required for otherwise lengthy procedures in a cost-efficient manner. Protocols have also been developed for efficient sequencing of a variety of templates submitted to the NEB core facility.

Alternative splicing creates sex-specific transcripts and truncated forms of the furin protease in the parasite Dirofilaria immitis.

Many extracellular proteins are activated by specific cleavage with an endoprotease. In nematodes, several proteins are cleaved after RX(K/R)R, the recognition site for the subtilisin-like proprotein convertases, furin and blisterase. To characterize furin in the parasitic nematode Dirofilaria immitis, we determined the sequence of the difur gene and its multiple transcripts. The gene spans 11 kb; encodes 16 exons and has a complex pattern of alternative splicing which generates at least 16 distinct mRNAs. The major transcript is a 4.4 kb mRNA which codes for a protein of 834 aa with an unusually long prodomain of 254 aa. Sex-specific splice variants of difur were observed by RT-PCR. The three female-specific and five male-specific transcripts are the first reported examples of sex-specific splicing in parasitic nematodes. This suggests that nematodes have sex-specific factors which regulate RNA splicing. Other splice variants are predicted to alter the phosphorylation and localization of the protease. Alternative splicing after the prodomain encodes a truncated protein that may be an inhibitor and/or substrate of Difurin.

Chemiluminescent detection of sequential DNA hybridizations to high-density, filter-arrayed cDNA libraries: a subtraction method for novel gene discovery.

A chemiluminescent approach for sequential DNA hybridizations to high-density filter arrays of cDNAs, using a biotin-based random priming method followed by a streptavidin/alkaline phosphatase/CDP-Star detection protocol, is presented. The method has been applied to the Brugia malayi genome project, wherein cDNA libraries, cosmid and bacterial artificial chromosome (BAC) libraries have been gridded at high density onto nylon filters for subsequent analysis by hybridization. Individual probes and pools of rRNA probes, ribosomal protein probes and expressed sequence tag probes show correct specificity and high signal-to-noise ratios even after ten rounds of hybridization, detection, stripping of the probes from the membranes and rehybridization with additional probe sets. This approach provides a subtraction method that leads to a reduction in redundant DNA sequencing, thus increasing the rate of novel gene discovery. The method is also applicable for detecting target sequences, which are present in one or only a few copies per cell; it has proven useful for physical mapping of BAC and cosmid high-density filter arrays, wherein multiple probes have been hybridized at one time (multiplexed) and subsequently "deplexed" into individual components for specific probe localizations.

Cloning and expression of the ApaLI, NspI, NspHI, SacI, ScaI, and SapI restriction-modification systems in Escherichia coli.

The genes encoding the ApaLI (5'-GTGCAC-3'), NspI (5'-RCATGY-3'), NspHI (5'-RCATGY-3'), SacI (5'-GAGCTC-3'), SapI (5'-GCTCTTCN1-3', 5'-N4GAAGAGC-3') and ScaI (5'-AGTACT-3') restriction-modification systems have been cloned in E. coli. Amino acid sequence comparison of M.ApaLI, M.NspI, M.NspHI, and M.SacI with known methylases indicated that they contain the ten conserved motifs characteristic of C5 cytosine methylases. NspI and NspHI restriction-modification systems are highly homologous in amino acid sequence. The C-termini of the NspI and NlaIII (5'-CATG-3') restriction endonucleases share significant similarity. 5mC modification of the internal C in a SacI site renders it resistant to SacI digestion. External 5mC modification of a SacI site has no effect on SacI digestion. N4mC modification of the second base in the sequence 5'-GCTCTTC-3' blocks SapI digestion. N4mC modification of the other cytosines in the SapI site does not affect SapI digestion. N4mC modification of ScaI site blocks ScaI digetion. A DNA invertase homolog was found adjacent to the ApaLI restriction-modification system. A DNA transposase subunit homolog was found upstream of the SapI restriction endonuclease gene.

Cloning and characterization of the Bg/II restriction-modification system reveals a possible evolutionary footprint.

Bg/II, a type II restriction-modification (R-M) system from Bacillus globigii, recognizes the sequence 5'-AGATCT-3'. The system has been cloned into E. coli in multiple steps: first the methyltransferase (MTase) gene, bglIIM, was cloned from B. globigii RUB561, a variant containing an inactivated endonuclease (ENase) gene (bglIIR). Next the ENase protein (R.BglII) was purified to homogeneity from RUB562, a strain expressing the complete R-M system. Oligonucleotide probes specific for the 5' end of the gene were then synthesized and used to locate bglIIR, and the gene was isolated and cloned in a subsequent step. The nucleotide sequence of the system has been determined, and several interesting features have been found. The genes are tandemly arranged, with bglIIR preceding bglIIM. The amino acid sequence of M.BglII is compared to those of other known MTases. A third gene encoding a protein with sequence similarity to known C elements of other R-M systems is found upstream of bglIIR. This is the first instance of a C gene being associated with an R-M system where the R and M genes are collinear. In addition, open reading frames (ORFs) resembling genes involved with DNA mobility are found in close association with BglII. These may shed light on the evolution of the R-M system.

Thermal cycle dideoxy DNA sequencing.

Thermal cycle dideoxy DNA sequencing eliminates the requirements for independent primer annealing and double-stranded DNA denaturation steps. The method enables sequencing from nanogram amounts of DNA from double-stranded and single-stranded PCR products, and plasmid or phage DNA templates. Thermal cycle sequencing also enables direct sequencing from bacterial colonies or phage plaques. Protocols using the Vent exo- DNA polymerase, helpful suggestions, and a troubleshooting guide are also presented.

Cloning, expression and sequence analysis of the SphI restriction-modification system.

SphI, a type II restriction-modification (R-M) system from the bacterium Streptomyces phaeochromogenes, recognizes the sequence 5'-GCATGC. The SphI methyltransferase (MTase)-encoding gene, sphIM, was cloned into Escherichia coli using MTase selection to isolate the clone. However, none of these clones contained the restriction endonuclease (ENase) gene. Repeated attempts to clone the complete ENase gene along with sphIM in one step failed, presumably due to expression of SphI ENase gene, sphIR, in the presence of inadequate expression of sphIM. The complete sphIR was finally cloned using a two-step process. PCR was used to isolate the 3' end of sphIR from a library. The intact sphIR, reconstructed under control of an inducible promoter, was introduced into an E. coli strain containing a plasmid with the NlaIII MTase-encoding gene (nlaIIIM). The nucleotide sequence of the SphI system was determined, analyzed and compared to previously sequenced R-M systems. The sequence was also examined for features which would help explain why sphIR unlike other actinomycete ENase genes seemed to be expressed in E. coli.

Cloning and sequence comparison of AvaI and BsoBI restriction-modification systems.

AvaI and BsoBI restriction endonucleases are isoschizomers which recognize the symmetric sequence 5'CYCGRG3' and cleave between the first C and second Y to generate a four-base 5' extension. The AvaI restriction endonuclease gene (avaIR) and methylase gene (avaIM) were cloned into Escherichia coli by the methylase selection method. The BsoBI restriction endonuclease gene (bsoBIR) and part of the BsoBI methylase gene (bsoBIM) were cloned by the "endo-blue" method (SOS induction assay), and the remainder of bsoBIM was cloned by inverse PCR. The nucleotide sequences of the two restriction-modification (RM) systems were determined. Comparisons of the predicted amino acid sequences indicated that AvaI and BsoBI endonucleases share 55% identity, whereas the two methylases share 41% identity. Although the two systems show similarity in protein sequence, their gene organization differs. The avaIM gene precedes avaIR in the AvaI RM system, while the bsoBI R gene is located upstream of bsoBI M in the BsoBI RM system. Both AvaI and BsoBI methylases contain motifs conserved among the N4 cytosine methylases.