HUMAN GENOME VARIATIONS AND EVOLUTION WITH A FOCUS ON THE ANALYSIS OF TRANSPOSABLE ELEMENTS

Genome sequence varies in numerous ways among individuals although the gross architecture is fixed for all humans. Retrotransposons create one of the most abundant structural variants in the human genome and are divided in many families, with certain members in some families, e.g., L1, Alu, SVA, and HERV-K, remaining active for transposition. Along with other types of genomic variants, retrotransponson-derived variants contribute to the whole spectrum of genome variants in humans. With the advancement of sequencing techniques, many human genomes are being sequenced at the individual level, fueling the comparative research on these variants among individuals. In this thesis, the evolution and functional impact of structural variations is examined primarily focusing on retrotransposons in the context of human evolution. The thesis comprises of three different studies on the topics that are presented in three data chapters. First, the recent evolution of all human specific AluYb members, representing the second most active subfamily of Alus, was tracked to identify their source/master copy using a novel approach. All human-specific AluYb elements from the reference genome were extracted, aligned with one another to construct clusters of similar copies and each cluster was analyzed to generate the evolutionary relationship between the members of the cluster. The approach resulted in identification of one major driver copy of all human specific Yb8 and the source copy of the Yb9 lineage. Three new subfamilies within the AluYb family – Yb8a1, Yb10 and Yb11 were also identified, with Yb11 being the youngest and most polymorphic. Second, an attempt to construct a relation between transposable elements (TEs) and tandem repeats (TRs) was made at a genome-wide scale for the first time. Upon sequence comparison, positional cross-checking and other relevant analyses, it was observed that over 20% of all TRs are derived from TEs. This result established the first connection between these two types of repetitive elements, and extends our appreciation for the impact of TEs on genomes. Furthermore, only 6% of these TE-derived TRs follow the already postulated initiation and expansion mechanisms, suggesting that the others are likely to follow a yet-unidentified mechanism. Third, by taking a combination of multiple computational approaches involving all types of genetic variations published so far including transposable elements, the first whole genome sequence of the most recent common ancestor of all modern human populations that diverged into different populations around 125,000-100,000 years ago was constructed. The study shows that the current reference genome sequence is 8.89 million base pairs larger than our common ancestor’s genome, contributed by a whole spectrum of genetic mechanisms. The use of this ancestral reference genome to facilitate the analysis of personal genomes was demonstrated using an example genome and more insightful recent evolutionary analyses involving the Neanderthal genome. The three data chapters presented in this thesis conclude that the tandem repeats and transposable elements are not two entirely distinctly isolated elements as over 20% TRs are actually derived from TEs. Certain subfamilies of TEs themselves are still evolving with the generation of newer subfamilies. The evolutionary analyses of all TEs along with other genomic variants helped to construct the genome sequence of the most recent common ancestor to all modern human populations which provides a better alternative to human reference genome and can be a useful resource for the study of personal genomics, population genetics, human and primate evolution.

A Study of Ordered Gene Problems Featuring DNA Error Correction and DNA Fragment Assembly with a Variety of Heuristics, Genetic Algorithm Variations, and Dynamic Representations

Ordered gene problems are a very common classification of optimization problems. Because of their popularity countless algorithms have been developed in an attempt to find high quality solutions to the problems. It is also common to see many different types of problems reduced to ordered gene style problems as there are many popular heuristics and metaheuristics for them due to their popularity. Multiple ordered gene problems are studied, namely, the travelling salesman problem, bin packing problem, and graph colouring problem. In addition, two bioinformatics problems not traditionally seen as ordered gene problems are studied: DNA error correction and DNA fragment assembly. These problems are studied with multiple variations and combinations of heuristics and metaheuristics with two distinct types or representations. The majority of the algorithms are built around the Recentering- Restarting Genetic Algorithm. The algorithm variations were successful on all problems studied, and particularly for the two bioinformatics problems. For DNA Error Correction multiple cases were found with 100% of the codes being corrected. The algorithm variations were also able to beat all other state-of-the-art DNA Fragment Assemblers on 13 out of 16 benchmark problem instances.

Regional variations of deformation mechanisms in the Lorrain Quartzite near Whitefish Falls, Ontario, Canada /

The quartzite microfabric found in the Lorrain Formation was studied across the La Cloche syncline, along a regional north-south transect along highway 6, near Whitefish Falls, Ontario. The complete stratigraphic sequence across the syncline is preserved, and is present on each fold limb. The lithostratigraphic units with the smallest grains size and lowest mica content are located close to the core of the fold, while coarser grained mica and feldspar rich units are situated at the northern and southern most extent of the transect. Deformation mechanisms vary with lithology and with position across the fold. Pressure solution appears to be the dominant deformation mechanism in the feldspathic, micaceous and ferruginous units. In the finer grained, mica poor white medium grained and cherty sandstone units, grain boundary migration (GBM) characteristics show dominance over those of pressure solution and show high amounts of fracturing which cut migrated boundaries and therefore post date GBM. All samples across the fold display a preferred orientation of quartz c-axes. The senses of asymmetry of fabrics are found to be similar across the syncline, with the exception of the ferruginous sandstone unit. Formation of these similar fabrics synmietries can not be the result of strain related to first order folding. The mica content appears to be related to the percentage of quartz lost due to pressure solution as a result of strain; the more mica present, the less quartz was lost. Calculations based on the shape of initial grains suggest that conservatively 30% of the quartz volume has been dissolved out of the Lorrain quartzite, and potentially migrated hundreds of meters to other members of the Huronian Supergroup as there was no meso or macroscopic evidence observed in outcrop.

Developmental variations in the peripheral erythrocytic system of the rainbow trout, Salmo gairdneri

The peripheral circulating erythrocytic system of the rainbow trout, l3 almo gairdner , was examined in vitro in relation differences in the morphology and multiple hemoglobin system organization of adult and juvenile red cells. Cells were separated by velocity sedimentation under unit gravity, a procedure requiring red cell exposure to an incubation medium for periods of at least three hours. Therefore , this must provide an environment in which red cells remain in a condition approximaing normalcy. Previous studies having demonstrated commonly employed media to be ineffective in this regard , a medium was developed through modification of Cortl and saline. One of the principal additions to this me dium , norepinephrine, altered cell regulation of intracellular calcium, magnesium and chloride concentrations. Catecholamine involvement was also suggeste d in the synthes is of hemoglobin . The procedure was found to separtate cells primarily by density and, to a lesser extent, by shape. Characterization of red cells revealed two subpopulations to exist . The first comprised the bulk of the cell population, and were of greater l ength, width, volume and major:minor axis ratio than the smaller population; these were adult cells. The later, juvenile cells were of smaller overall size and were more spherical in shape . Juvenile cells also possessed fewer electrophore tpically distinguishable isomorphs than did adults with only eight of eleven hemoglobin component s typically found With maturation,hemoglobin complement with the development of three more bands. The total complement of the adult cell contained 7 cathodal bands and four anodal hemoglobin isomorphs. Bands acquired with maturation comprised the smallest percentage of the cells hemoglobin. each averaging less than one-percent of the total. Whether these additional bands are derived through degradation and reaggregation of existing components or are the product of pe gQy2 synthesis is not yet known.

Thermoacclimatory variations in the activities of enzymes implicated in ion transport in the rainbow trout, salmo gairdneri

Two groups of rainbow trout were acclimated to 20 , 100 , and 18 o C. Plasma sodium, potassium, and chloride levels were determined for both. One group was employed in the estimation of branchial and renal (Na+-K+)-stimulated, (HC0 3-)-stimulated, and CMg++)-dependent ATPase activities, while the other was used in the measurement of carbonic anhydrase activity in the blood, gill and kidney. Assays were conducted using two incubation temperature schemes. One provided for incubation of all preparations at a common temperature of 2S oC, a value equivalent to the upper incipient lethal level for this species. In the other procedure the preparations were incubated at the appropriate acclimation temperature of the sampled fish. Trout were able to maintain plasma sodium and chloride levels essentially constant over the temperature range employed. The different incubation temperature protocols produced different levels of activity, and, in some cases, contrary trends with respect to acclimation temperature. This information was discussed in relation to previous work on gill and kidney. The standing-gradient flow hypothesis was discussed with reference to the structure of the chloride cell, known thermallyinduced changes in ion uptake, and the enzyme activities obtained in this study. Modifications of the model of gill lon uptake suggested by Maetz (1971) were proposed; high and low temperature models resulting. In short, ion transport at the gill at low temperatures appears to involve sodium and chloride 2 uptake by heteroionic exchange mechanisms working in association w.lth ca.rbonlc anhydrase. G.l ll ( Na + -K + ) -ATPase and erythrocyte carbonic anhydrase seem to provide the supplemental uptake required at higher temperatures. It appears that the kidney is prominent in ion transport at low temperatures while the gill is more important at high temperatures. 3 Linear regression analyses involving weight, plasma ion levels, and enzyme activities indicated several trends, the most significant being the interrelationship observed between plasma sodium and chloride. This, and other data obtained in the study was considered in light of the theory that a link exists between plasma sodium and chloride regulatory mechanisms.

Thermoaclimatory variations in the microenvironment of hemoglobin in the rainbow trout, Salmo gairdneri and the carp, Cyprinus carpio

Several inorganic substances (e.g., C£ , Mg , Ca , H ) are potent negative modulators of hemoglobin-oxygen affinity. To evaluate the possibility that potentially adaptive changes in the red cell ionic environment of hemoglobin may take place during acclimation of fishes to increased environmental temperature, hematological status (hemoglobin, hematocrit, red cell numbers, mean erythrocytic volume and hemoglobin content), plasma + + 2+ 2+ and packed red cell electrolyte levels (Na , K , Ca , Mg , C£ ) were evaluated in summer and winter populations of the stenothermal rainbow trout, Salmo gairdneri, following acclimation to 2°, 10°, 18°C, and in a spring population of eurythermal carp, Cyprinus carpio, held at 2°, 16° and 30°C. From these data cell ion concentrations and ion:hemoglobin ratios were estimated. In view of the role of red cell carbonic anhydrase in the reductions of blood C02 tensions and the recruitment of Na and C£~ lost by fishes, a preliminary investigation of thermoacclimatory changes in the activity of this system in rainbow trout erythrocytes was conducted. Few changes in hematological status were encountered following acclimation. There was, however, some evidence of weight-specific differential hematological response in carp. This lead to markedly greater increases in hemoglobin, hematocrit and red cell numbers in smaller rather than in larger specimens at higher temperatures; variations which were 2+ well correlated with changes in plasma Ca . Plasma composition in summer trout was not altered by acclimation. In winter trout plasma Na and K increased at higher temperatures. Carp were characterized by increases in plasma calcium, and reductions in sodium and magnesium under these conditions. Several significant seasonal differences in plasma ion levels were observed in the trout. (n) In trout, only erythrocytic K and K :Hb were altered by acclimation, rising at higher temperatures. In carp Na , Na :Hb, C£~ and C£~:Hb in- 2+ 2+ creased with temperature, while Mg and Mg :Hb declined. Changes in overall ionic composition in carp red cells were consistent with increases in H content. In both species significant reciprocal variations in C£~ 2+ - + and Mg were found. In mammalian systems increases in C£ and H reduce hemoglobin-oxygen affinity by interaction with hemoglobin. Reduction in 2+ 2+ Mg maximizes organophosphate modulator availability by decreasing ATP»Mg complex formation. Thus, the changes observed may be of adaptive value in reducing hemoglobin-oxygen affinity, and facilitating oxygen release to cells at higher temperatures. Trout appear to maintain a high chloridelow magnesium state over the entire thermal tolerance zone. Carp, however, achieved this state only at higher temperatures. In both species mean erythrocytic volume was decreased at higher temperatures and this may facilitate branchial oxygen loading. Since mean erythrocytic volume was inversely related to red cell ion content, it is hypothesized that reductions in cell volume are achieved by export of some unidentified solute or solutes. Variations in the carbonic anhydrase activity that could be attributed to the thermoacclimatory process were quite modest. On the other hand, assays performed at the temperature of acclimation showed a large temperature effect where under in vivo conditions of temperature fish acclimated to higher temperatures might be expected to have higher activities. Furthermore, since hematocrit increased with temperature in these fish, while carbonic anhydrase is present only in the erythrocyte, the whole blood levels of this enzyme are expected to increase and further augment the temperature effect. This, in turn, could aid in the reduction of C02 (111) tension and increase the production of H and HC0~~ used in the active uptake of Na and C£ at higher temperatures.

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.