Improving Short DNA Sequence Alignment with Parallel Computing

Variations in different types of genomes have been found to be responsible for a large degree of physical diversity such as appearance and susceptibility to disease. Identification of genomic variations is difficult and can be facilitated through computational analysis of DNA sequences. Newly available technologies are able to sequence billions of DNA base pairs relatively quickly. These sequences can be used to identify variations within their specific genome but must be mapped to a reference sequence first. In order to align these sequences to a reference sequence, we require mapping algorithms that make use of approximate string matching and string indexing methods. To date, few mapping algorithms have been tailored to handle the massive amounts of output generated by newly available sequencing technologies. In otrder to handle this large amount of data, we modified the popular mapping software BWA to run in parallel using OpenMPI. Parallel BWA matches the efficiency of multithreaded BWA functions while providing efficient parallelism for BWA functions that do not currently support multithreading. Parallel BWA shows significant wall time speedup in comparison to multithreaded BWA on high-performance computing clusters, and will thus facilitate the analysis of genome sequencing data.

Factors affecting DNA uptake by mammalian cells

The ability to introduce DNA and express custom DNA sequences in bacteria opened the door for improvements in a large number of fields including agriculture, pharmacology, medicine, nutrition, etc. The ability to introduce foreign DNA sequences into mammalian cells in an efficient manner would have a large impact on therapeutic applications especially gene therapy. The methods in use today suffer from low efficiencies and sometimes toxicity. In this work a number of factors were evaluated for their effect onONA uptake efficiency. The factors studied included exposure to sublethal concentration of hydrogen peroxide which have been show to lead to destabilisation ofthe lysosomes. These exposures have proven to be very toxic to cells when combined with either the calcium phosphate or the lipofectAMINE® transfection methods. Another factor evaluated was exposure to Electro-Magnetic Fields (EMF). This was fuelled by the fact that EMF have been shown to mediate a number of effects on cell structure and/or physiology. EMF exposure by itself was not sufficient to induce the cells to pick up the DNA, therefore its effect on calcium phosphate and lipofectAMINE® was tested. Although some positive results were obtained, the variability of these results exceeded by far any observed enhancements which discouraged any further work on EMF. Also tested was the possible effect the presence of the cytomegalovirus (CMV) sequence might have on DNA uptake (based on previous results in this lab). It was found that the presence ofCMV in the DNA sequence does not enhance uptake or slow down degradation of the internalised DNA. The final factor tested was the effect of basic amino acids on transfection efficiency. It was found that arginine can enhance DNA uptake by about 170% v/ith calcium phosphate and about 200% with LipofectAMINE®. A model was proposed to explain the effect of arginine as well as the lack of effect from other amino acids.

Isolation and analysis of a corprinus cinereus DNA fragment showing homology to a fimbrial cDNA from ustilago violacea

Surface proteinaceous fibrils, termed fimbriae, were first identified on gram negative bacteria in the 1940s. Fungal fimbriae, discovered some 25 years later, are found on members of all fungal classes. In the present study, polyclonal antiserum raised against the fimbrial proteins of U. vio/acea were used in order to identify antigenically related proteins from Coprinus cinereus and Schizophy//um commune. Two polypeptides with molecular masses of 37 and 39 kDa from C. cinereus were observed and confirm earlier results. A single previously unidentified 50 kDa polypeptide in S. commune crossreacted with the antiserum. The 50 kDa protein was found to consist of 3 isoforms with isoelectric points ranging from 5.6 to 5.8. A fimbrial cDNA derived from U. vio/acea was used to identify DNA restriction fragments from C. cinereus and S. commune showing homology to the fimbrial transcript of U. vio/acea. Heterologous hybridization with this cDNA was used in order to screen a C. cinereus genomic DNA library. A single clone, A2-3A, with a 14 kbp insert showed strong homology to the pfim3-1 cDNA. The region of homology, a 700 bp Xba I fragment, was subcloned into pUG19. This plasmid was refered to as pXX8. DNA sequence determinations of pXX8 and adjacent fragments from A2-3A suggested that the cloned DNA was a portion of the rONA repeat encoding the small subunit rRNA. DNA sequence analysis of pfim3-1 yielded an incomplete open reading frame. The predicted amino acid sequence codes for a 206 amino acid, 22 kDa polypeptide which contains a domain similar to a transmembrane domain from rat leukocyte antigen, GDS3. As well, an untranslated 576 nucleotide domain showed 81 % homology to pXX8 and 830/0 homology to the 188 rRNA sequence of Ustilago maydis. This sequence was found adjacent to a region of adenine-thymine base pairs presumed to represent the polyadenylation sequence of the fimbrial transcript. The size and extent of homology is sufficient to account for the hybridization of pfim3-1 to rDNA. It is suggested that this domain represents a completely novel regulatory domain within eukaryotes that may enable the observed rapid regeneration of fimbriae in U. violacea.

Partial characterization of the cloned dihydrofolate reductase gene of Saccharomyces cerevisiae

The cloned dihydrofolate reductase gene of Saccharomyces cerevisiae (DFR 1) is expressed in Escherichia coli. Bacterial strain JF1754 transformed with plasmids containing DFR 1 is at least 5X more resistant to inhibition by the folate antagonist trimethoprim. Expression of yeast DFR 1 in E. coli suggests it is likely that the gene lacks intervening sequences. The 1.8 kbp DNA fragment encoding yeast dhfr activity probably has its own promotor, as the gene is expressed in both orientations in E. coli. Expression of the yeast dhfr gene cloned into M13 viral vectors allowed positive selection of DFR 1 - M13 bacterial transfectants in medium supplemented with trimethoprim. A series of nested deletions generated by nuclease Bal 31 digestion and by restriction endonuclease cleavage of plasmids containing DFR 1 physically mapped the gene to a 930 bp region between the Pst 1 and Sal 1 cut sites. This is consistent with the 21,000 molecular weight attributed to yeast dhfr in previous reports. From preliminary DNA sequence analysis of the dhfr DNA fragment the 3' terminus of DFR 1 was assigned to a position 27 nucleotides from the Eco Rl cut site on the Bam Hi - Eco Rl DNA segment. Several putative yeast transcription termination consensus sequences were identified 3' to the opal stop codon. DFR 1 is expressed in yeast and it confers resistance to the antifolate methotrexate when the gene is present in 2 - 10 copies per cell. Plasmid-dependent resistance to methotrexate is also observed in a rad 6 background although the effect is somewhat less than that conferred to wild-type or rad 18 cells. Integration of DFR 1 into the yeast genome showed an intermediate sensitivity to folate antagonists. This may suggest a gene dosage effect. No change in petite induction in these yeast strains was observed in transformed cells containing yeast dhfr plasmids. The sensitivity of rad 6 , rad 18 and wild-type cell populations to trimethoprim were unaffected by the presence of DFR 1 in transformants. Moreover, trimethoprim did not induce petites in any strain tested, which normally results if dhfr is inhibited by other antifolates such as methotrexate. This may suggest that the dhfr enzyme is not the only possible target of trimethoprim in yeast. rad 6 mutants showed a very low level of spontaneous petite formation. Methotrexate failed to induce respiratory deficient mutants in this strain which suggested that rad 6 might be an obligate grande. However, ethidium bromide induced petites to a level approximately 50% of that exhibited by wild-type and rad 18 strains.

Cryptic species status of Anopheles (Diptera: Culicidae) mosquitoes in Canada using a multidisciplinary approach

Many species of Anopheles mosquitoes (Diptera: Culicidae) are now recognized as species complexes whose members are often indistinguishable morphologically but identifiable based on ecological, genetic, or behavioural data. Because the members of species complexes often differ in their vector potential, accurate identification of vector species is essential for successful mosquito control. To investigate the cryptic species status of Anopheles mosquitoes in Canada, specimens were collected from across the country and examined using morphological, molecular, and ecological data. Six of the seven traditionally recognised species from Canada were collected from locations in British Columbia, Quebec, Newfoundland and Labrador, and throughout Ontario, including Anopheles barberi, An. earlei, An. freeborni, An. punctipennis, An. quadrimaculatus s.l., and An. walkeri. Variation in polymorphic traits within An. earlei, An. punctipennis, and An. quadrimaculatus s.l. were quantified and egg morphology examined using scanning electron microscopy. Morphological identification of adult and larval specimens suggested that two described cryptic species, An. perplexens and An. smaragdinus, were present in Canada. DNA sequence data were analysed for evidence of cryptic species using three molecular markers: COl, ITS2, and ITS!. Intraspecific COl variation was very low in most species «1 %), except for An. punctipennis with 2% sequence divergence between those from British Columbia (BC) and Ontario (ON), and An. walkeri with 7% sequence divergence between populations from Manitoulin Island (NO) and Long Point Provincial Park (LP). Similar patterns were also seen using ITS2 and ITS 1. Therefore, molecular data revealed the presence of two putative cryptic species within two species examined (i.e., An. walkeri and An. punctipennis), corresponding to collection location (i.e., NO vs. LP and BC vs. ON, respectively). Surprisingly, there was no molecular support for the presence of either An. perplexens or An. smaragdinus in Canada despite the morphological assessments. Ecological data from all collection sites were recorded and are available in an online database designed to manage all collection and identification data. Current bionomic information, including regional abundance, larval habitat, and species associations, was determined for each species. This multidisciplinary study of Anopheles mosquitoes is the first detailed investigation of these potential disease vectors in Canada and demonstrates the importance of an integrated approach to anopheline systematics that includes molecular data.

Construction of bovine adenovirus type 3 E1 and E3, substitution plasmids and the sequencing analysis of DNA pol and pTP /

The relative ease to concentrate and purify adenoviruses, their well characterized mid-sized genome, and the ability to delete non-essential regions from their genome to accommodate foreign gene, made adenoviruses a suitable candidate for the construction of vectors. The use of adenoviral vectors in gene therapy, vaccination, and as a general vector system for expressing foreign genes have been documented for some time. In this study, the objective was to rescue a BAV3 E1 or E3 recombinant vector carrying the kanamycin resistant gene, a dominant selectable marker with useful applications in studying vectored gene expression in mammalian cells. To accomplish the objective of this study, more information about BAV3 DNA sequences was required in order to make the manipulation of the virus genome accessible. Therefore, sequencing of the BAV3 genome from 1 1 .7% to 30.8% was carried out. Analysis of the determined sequences revealed the primary structure of important viral gene products coded by E2 including BAV3 DNA pol and precursor to terminal protein. Comparative analysis of these proteins with their counterparts from human and non human adenoviruses revealed important insights as to the evolutionary lineage of BAV3. In order to insert the kanamycin resistance gene in either E1 or E3, it was necessary to delete BAV3 sequences to accommodate the foreign gene so as not to exceed the limit of the packaging capacity of the virus. To construct a recombinant BAV3 in which a foreign gene was inserted in the deleted E1 region, an E1 shuttle vector was constructed. This involved the deletion from the viral sequences a region between 1.3% to 9% and inserting the kanamycin resistance gene to replace the deletion. The E1 shuttle vector contained the left (0%- 53.9%) segment of the genome and was expected to generate BAV3 recombinants that can be grown and propagated in cells that can complement the missing E1 functions. To construct a similar shuttle vector for E3 deletion, DNA sequences extending from 78.9% to 82.5% (1281 bp) were deleted from within the E3 region that had been cloned into a plasmid vector. The deleted region corresponds to those that have been shown to be non-essential for viral replication in cell culture. The resulting plasmid was used to construct another recombinant plasmid with BAV3 DNA sequences extending from 37.1% to 100% and with a deletion of E3 sequences that were replaced by kanamycin resistance gene. This shuttle plasmid was used in cotransfections with digested viral DNA in an attempt to rescue a recombinant BAV3 carrying the kanamycin resistance gene to replace the deleted E3. In spite of repeated attempts of transfection, El or E3 recombinant BAV3 were not isolated. It seems that other approaches should be applied to make a final conclusion on BAV3 infectivity.

The Elucidation of the involvement of endonuclease DNase y in reducing DNA transfection efficiency in mammalian cells

Gene therapy is predicated upon efficient gene transfer. While viral vectors are the method of choice for transformation efficiency, the immunogenicity and safety concerns remain problematic. Non-viral vectors, on the other hand, have shown high degrees of safety and are mostly non-immunogenic in nature. However, non-viral vectors usually suffer from low levels oftransformation efficiency and transgene expression. Thus, increasing transformation efficiency ofnon-viral vectors, in particular by calcium phosphate co-precipitation technique, is a way of generating a suitable vector for gene therapy and is the aim of this study. It is a long known fact that different cell lines have different transfection efficiencies regardless oftransfection methodology (Lin et a!., 1994). Using commonly available cell lines Madine-Darby Bovine Kidney (MDBK), HeLa and Human Embryonic Kidney (HEK-293), we have shown a decreasing trend ofDNase activity based on a plasmid digestion assay. From densitometry studies, as much as a 40% reduction in DNase activity was observed when comparing HEK-293 (least active) to MDBK (most active). Using various biochemical assays, it was determined that DNase y, in particular, was expressed more highly in MDBK cells than both HeLa and HEK-293. Upon cloning of the bovine DNase y gene, we utilized the sequence information to construct antisense expressing plasmids via both traditional antisense RNA (pASDGneoM) and siRNA (psiRNA-S4, psiRNA-S11 and psiRNA-S16). For the construction ofpASDGneoM, the 3' end of the DNase y was inserted in opposite orientation under a cytomegalovirus (CMV) promoter such that the expression ofRNA complementary to the DNase 2 ymRNA occurred. For siRNA plasmids, the sequence was screened to yield optimal short sequences for siRNA inhibition. The silencing ofbovine DNase y led to an increase in transfection efficiency based on traditional calcium phosphate co-precipitation technique; stable clones of siRNA-producing MDBK cell lines (psiRNA-S4 Bland psiRNA-S4 B4) both demol).strated 4-fold increases in transfection efficiency. Furthermore, serial transfection of antisense DNase y plasmid pASDGneoM and reporter pCMV-~ showed a maximum of 8-fold increase in transfection efficiency when the two separate transfections were carried out 4 hours apart (i.e. transfection ofpASDGneoM, separated by four hours, then transfection ofpCMV-~). Together, these results demonstrate the involvement ofDNase y in reducing transfection efficiency, at least by traditional calcium phosphate technique.

Isolation of DNA sequences with homology to a degenerate FaRP oligonucleotide from the crayfish Procambarus clarkii

The neuropeptide Th1RFamide with the sequence Phe-Met-Arg-Phe-amide was originally isolated in the clam Macrocallista nimbosa (price and Greenberg, 1977). Since its discovery, a large family ofFl\1RFamide-related peptides termed FaRPs have been found to be present in all major animal phyla with functions ranging from modulation of neuronal activity to alteration of muscular contractions. However, little is known about the genetics encoding these peptides, especially in invertebrates. As FaRP-encoding genes have yet to be investigated in the invertebrate Malacostracean subphylum, the isolation and characterization ofFaRP-encoding DNA and mRNA was pursued in this project. The immediate aims of this thesis were: (1) to amplify mRNA sequences of Procambarus clarkii using a degenerate oligonucleotide primer deduced from the common amino acid sequence ofisolated Procambarus FaRPS, (2) to determine if these amplification products encode FaRP gene sequences, and (3) to create a selective cDNA library of sequences recognized by the degenerate oligonucleotide primer. The polymerase chain reaction - rapid amplification of cDNA ends (PCR-RACE) is a procedure in which a single gene-specific primer is used in conjunction with a generalized 3' or 5' primer to amplify copies ofthe region between a single point in the transcript and the 3' or 5' end of cDNA of interest (Frohman et aI., 1988). PCRRACE reactions were optimized with respect to primers used, buffer composition, cycle number, nature ofgenetic substrate to be amplified, annealing, extension and denaturation temperatures and times, and use of reamplification procedures. Amplification products were cloned into plasmid vectors and recombinant products were isolated, as were the recombinant plaques formed in the selective cDNA library. Labeled amplification products were hybridized to recombinant bacteriophage to determine ligated amplification product presence. When sequenced, the five isolated PCR-RACE amplification products were determined not to possess FaRP-encoding sequences. The 200bp, 450bp, and 1500bp sequences showed homology to the Caenorhabditis elegans cosmid K09A11, which encodes for cytochrome P450; transfer-RNA; transposase; and tRNA-Tyr, while the 500bp and 750bp sequences showed homology with the complete genome of the Vaccinia virus. Under the employed amplification conditions the degenerate oligonucleotide primer was observed to bind to and to amplify sequences with either 9 or 10bp of 17bp identity. The selective cDNA library was obselVed to be of extremely low titre. When library titre was increased, white. plaques were isolated. Amplification analysis of eight isolated Agt11 sequences from these plaques indicated an absence of an insertion sequence. The degenerate 17 base oligonucleotide primer synthesized from the common amino acid sequence ofisolated Procambarus FaRPs was thus determined to be non-specific in its binding under the conditions required for its use, and to be insufficient for the isolation and identification ofFaRP-encoding sequences. A more specific primer oflonger sequence, lower degeneracy, and higher melting temperature (TJ is recommended for further investigation into the FaRP-encoding genes of Procambarlls clarkii.

Isolation and characterization of genomic DNA sequences that enhance the stability of plasmid DNA in mammalian cells /

The ease of production and manipulation has made plasmid DNA a prime target for its use in gene transfer technologies such as gene therapy and DNA vaccines. The major drawback of plasmid however is its stability within mammalian cells. Plasmid DNA is usually lost by cellular mechanisms or as a result of mitosis by simple dilution. This study set out to search for mammalian genomic DNA sequences that would enhance the stability of plasmid DNA in mammalian cells.Creating a plasmid based genomic DNA library, we were able to screen the human genome by transfecting the library into Human Embryonic Kidney (HEK 293) Cells. Cells that contained plasmid DNA were selected, using G418 for 14 days. The resulting population was then screened for the presence of biologically active plasmid DNA using the process of transformation as a detector.A commercially available plasmid DNA isolation kit was modified to extract plasmid DNA from mammalian cells. The standardized protocol had a detection limit of -0.6 plasmids per cell in one million cells. This allowed for the detection of 45 plasmids that were maintained for 32 days in the HEK 293 cells. Sequencing of selected inserts revealed a significantly higher thymine content in comparison to the human genome. Sequences with high A/T content have been associated with Scaffold/Matrix Attachment Region (S/MAR) sequences in mammalian cells. Therefore, association with the nuclear matrix might be required for the stability of plasmids in mammalian cells.

Alternative splicing of DNA polymerase beta in vertebrates

Alternative splicing (AS) is the predominant mechanism responsible for increasing eukaryotic transcriptome and proteome complexity. In this phenomenon, numerous mRNA transcripts are produced from a single pre-mRNA sequence. AS is reported to occur in 95% of human multi-exon genes; one specific gene that undergoes AS is DNA polymerase beta (POLB). POLB is the main DNA repair gene which performs short patch base excision repair (BER). In primate untransformed primary fibroblast cell lines, it was determined that the splice variant (SV) frequency of POLB correlates positively with species lifespan. To date, AS patterns of POLB have only been examined in mammals primarily through the use of cell lines. However, little attention has been devoted to investigating if such a relationship exists in non-mammals and whether cell lines reflect what is observed in vertebrate tissues. This idea was explored through cloning and characterization of 1,214 POLB transcripts from four non-mammalian species (Gallus gallus domesticus, Larus glaucescens, Xenopus laevis, and Pogona vitticeps) and two mammalian species (Sylvilagus floridanus and Homo sapiens) in two tissue types, liver and brain. POLB SV frequency occurred at low frequencies, < 3.2%, in non-mammalian tissues relative to mammalian (>20%). The highest POLB SV frequency was found in H. sapiens liver and brain tissues, occurring at 65.4% and 91.7%, respectively. Tissue specific AS of POLB was observed in L. glaucescens, P. vitticeps, and H. sapiens, but not G. gallus domesticus, X. laevis and S. floridanus.The AS patterns of a second gene, transient receptor potential cation channel subfamily V member 1 (TRPV1), were compared to those of POLB in liver and brain tissues of G. gallus domesticus, X. laevis and H. sapiens. This comparison was performed to investigate if any changes (either increase or decrease) observed in the AS of POLB were gene specific or if they were tissue specific, in which case similar changes in AS would be seen in POLB and TRPV1. Analysis did not reveal an increase or decrease in both the AS of POLB and TRPV1 in either the liver or brain tissues of G. gallus domesticus and H. sapiens. This result suggested that the AS patterns of POLB were not influenced by tissue specific rates of AS. Interestingly, an increase in the AS of both genes was only observed in X. laevis brain tissue. This result suggests that AS in general may be increased in the X. laevis brain as compared to liver tissue. No positive correlation between POLB SV frequency and species lifespan was found in non-mammalian tissues. The AS patterns of POLB in human primary untransformed fibroblast cell lines were representative of those seen in human liver tissue but not in brain tissue. Altogether, the AS patterns of POLB from vertebrate tissues and primate cell lines revealed a positive correlation between POLB SV frequency and lifespan in mammals, but not in non-mammals. It appears that this positive correlation does not exist in vertebrate species as a whole.

DNA-Directed DNA Polymerases Evolve From Reverse Transcriptase

Scientists have been debating for decades the origin of life on earth. A number of hypotheses were proposed as to what emerged first RNA or DNA; with most scientists are in favour of the "RNA World" hypothesis. Assuming RNA emerged first, it fellow that the RNA polymerases would've appeared before DNA polymerases. Using recombinant DNA technology and bioinformatics we undertook this study to explore the relationship between RNA polymerases, reverse transcriptase and DNA polymerases. The working hypothesis is that DNA polymerases evolved from reverse transcriptase and the latter evolved from RNA polymerases. If this hypothesis is correct then one would expect to find various ancient DNA polymerases with varying level of reverse transcriptase activity. In the first phase of this research project multiple sequence alignments were made on the protein sequence of 32 prokaryotic DNA-directed DNA polymerases originating from 11 prokaryotic families against 3 viral reverse transcriptase. The data from such alignments was not very conclusive. DNA polymerases with higher level of reverse transcriptase activity were non-confined to ancient organisms, as one would've expected. The second phase of this project was focused on conditions that may alter the DNA polymerase activity. Various reaction conditions, such as temperature, using various ions (Ni2+, Mn2+, Mg2+) were tested. Interestingly, it was found that the DNA polymerase from the Thermos aquatics family can be made to copy RNA into DNA (i.e. reverse transcriptase activity). Thus it was shown that under appropriate conditions (ions and reactions temperatures) reverse transcriptase activity can be induced in DNA polymerase. In the third phase of this study recombinant DNA technology was used to generate a chimeric DNA polymerase; in attempts to identify the region(s) of the polymerase responsible for RNA-directed DNA polymerase activity. The two DNA polymerases employed were the Thermus aquatic us and Thermus thermophiles. As in the second phase various reaction conditions were investigated. Data indicated that the newly engineered chimeric DNA polymerase can be induced to copy RNA into DNA. Thus the intrinsic reverse transcriptase activity found in ancient DNA polymerases was localized into a domain and can be induced via appropriate reaction conditions.