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.

Intracellular antioxidant and DNA repair enzymes as correlates of stress resistance and longevity in vertebrates

In animals, both stress resistance and longevity appear to be influenced by the insulin/insulin-like growth factor-l signaling (lIS) pathway, the basic organization of which is highly conserved from invertebrates to vertebrates. Reduced lIS or genetic disruption of the lIS pathway leads to the activation of forkhead box transcription factors, which is thought to upregulate the expression of genes involved in enhancing stress resistance, including perhaps key antioxidant enzymes as well as DNA repair enzymes. Enhanced antioxidant and DNA repair capacities may underlie the enhanced cellular stress resistance observed in long-lived animals, however little data is available that directly supports this idea. I used three. experimental approaches to test the association of intracellular antioxidant and DNA base excision repair (BER) capacities with stress resistance and longevity: (1) a comparison of multiple vertebrate endotherm species of varying body masses and longevities; (2) a comparison of long-lived Snell dwarf mice and their normallittermates; and (3) a comparison of hypometabolic animals undergoing hibernation or estivation with their active counterparts. The activities of the five major intracellular antioxidant enzymes as well as the two rate-limiting enzymes in the BER pathway, apurininc/apyrimidinic (AP) endonuclease and polymerase ~, were measured. These measurements were performed in one or more of the following: (1) cultured dermal fibroblasts; (2) brain tissue; (3) heart tissue; (4) liver tissue. My results indicate that antioxidant enzymes are not universally upregulated in association with enhanced stress resistance and longevity. I also did not find that BER enzyme activity was positively correlated with longevity, in an inter-species context, though there was evidence for enhanced BER in long-lived Snell dwarf mice. Thus, while there were instances in which enhanced antioxidant and BER enzyme activities were associated with increased stress resistance and/or longevity, this was not universally the case, indicating that other mechanisms must be involved. These results suggest the need to re-examine existing 'oxidative stress' hypotheses of longevity and probe further into the molecular physiology of longevity to discover its mechanistic basis.

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.

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.

New Contig Creation Algorithm for the de novo DNA Assembly Problem

DNA assembly is among the most fundamental and difficult problems in bioinformatics. Near optimal assembly solutions are available for bacterial and small genomes, however assembling large and complex genomes especially the human genome using Next-Generation-Sequencing (NGS) technologies is shown to be very difficult because of the highly repetitive and complex nature of the human genome, short read lengths, uneven data coverage and tools that are not specifically built for human genomes. Moreover, many algorithms are not even scalable to human genome datasets containing hundreds of millions of short reads. The DNA assembly problem is usually divided into several subproblems including DNA data error detection and correction, contig creation, scaffolding and contigs orientation; each can be seen as a distinct research area. This thesis specifically focuses on creating contigs from the short reads and combining them with outputs from other tools in order to obtain better results. Three different assemblers including SOAPdenovo [Li09], Velvet [ZB08] and Meraculous [CHS+11] are selected for comparative purposes in this thesis. Obtained results show that this thesis’ work produces comparable results to other assemblers and combining our contigs to outputs from other tools, produces the best results outperforming all other investigated assemblers.