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.