Polyploidy, or whole genome duplication (WGD), is a fundamental evolutionary mechanism for diverse organisms including many plants and some animals. Duplicate genes from WGD are a major source of expression and functional diversity. However, the biological and evolutionary mechanisms for gene expression changes within and between species following WGD are poorly understood. Using genome-wide gene expression microarrays and high-throughput sequencing technology, I studied the genetic and evolutionary mechanisms for gene expression changes in synthetic and natural allopolyploids that are derived from hybridization between closely related species. To investigate evolutionary fate of duplicate genes, I tested how duplicate genes respond to developmental and environmental changes within species and how ancient duplicate genes contribute to gene expression diversity in resynthesized allopolyploids. We found that expression divergence between gene duplicates was significantly higher in response to environmental stress than to developmental process. Furthermore, duplicate genes related to external stresses showed higher expression divergence between two closely related species and in resynthesized and natural allotetraploids than single-copy genes. A slow rate of expression divergence of duplicate genes during development may offer dosage-dependent selective advantage, whereas a high rate of expression divergence between gene duplicates in response to external changes may enhance adaptation. To investigate molecular mechanisms of expression diversity among allopolyploids, I analyzed high-throughput sequencing data of small RNAs in allopolyploids and their progenitors. Small interfering RNAs (siRNAs) induce epigenetic modification and gene silencing of repeats, while microRNAs (miRNAs) and trans-acting siRNAs (ta-siRNAs) induce expression modulation of protein coding genes. Our data showed that siRNA populations in progenitors were highly maintained in allopolyploids, and alteration of miRNA abundance in allopolyploids was significantly correlated with expression changes of miRNA target genes. These results suggest that stable inheritance of parental siRNAs in allopolyploids helps maintain genome stability in response to genome duplication, whereas expression diversity of miRNAs leads to interspecies variation in gene expression, growth, and development. Results from these research objectives show that genome-wide analysis of high throughput gene expression and small RNAs provides new insights into molecular and evolutionary mechanisms for gene expression diversity and phenotypic variation between closely related species and in the new allopolyploids. / text

Mycobacterium tuberculosis (Mtb) poses a global health catastrophe that has been compounded by the emergence of highly drug resistant Mtb strains. We used whole genome sequencing (WGS) to directly compare the accumulation of mutations in Mtb isolated from cynomolgus macaques with active, latent and early reactivation disease. Based on the distribution of single nucleotide polymorphisms (SNPs) observed, we calculated the mutation rates for these disease states. Our data suggest that during latency, Mtb acquires a similar number of chromosomal mutations as would be expected to emerge in a logarithmically growing culture over the same period of time despite reduced bacterial replication during latent infection. The pattern of polymorphisms suggests that the mutational burden in vivo is due to oxidative DNA damage. We next sought to determine why some strains of Mtb are preferentially associated with high-level drug resistance. We demonstrate that Mtb strains from the East Asian lineage acquire drug resistances in vitro more quickly than Mtb strains from the Euro-American lineage. Their higher drug resistance rate in vitro reflects a higher basal mutation. Moreover, the in vitro mutation rate correlates well with the bacterial mutation rate in humans as determined by whole genome sequencing of clinical isolates. Finally, using an agent-based model, we show that the observed differences in mutation rate predict a significantly higher probability of multi-drug resistance in patients infected with East Asian lineage strains of Mtb. Lastly, we sought to determine the mechanisms Mtb uses to proofread nascently polymerized DNA. Through fluctuation analysis of deletion mutants of two potential \(polIII\epsilon\) homologs, we demonstrate that neither is responsible for the maintenance of DNA replication fidelity. To explore the possibility that one of these homologs, Rv3711c, participates in an unknown redundant pathway, we used transposon capture and sequence (TraCS) to identify genes conditionally essential in an Rv3711c deletion mutant. Our analysis suggests that while Rv3711c does not participate in proofreading, it may act in an alternative novel DNA repair pathway. Taken together, our fluctuation analysis and TraCS data suggest that mycobacteria do not use canonical methods of proofreading to maintain genomic fidelity.

Duplicated genes are present in the teleost fish lineage owing to a whole-genome duplication (WGD) event that occured ~ 230-400 million years ago. In the duplication-degeneration-complementation (DDC) model, partitioning of ancestral functions (subfunctionalization) and acquisition of novel functions (neofunctionalization) have been proposed as principal processes for the retention of duplicated genes in the genome. The DDC model was tested by analyzing the differential tissue-specific distribution of transcripts for the duplicated fatty acid-binding protein 10 (fabp10) genes in embryos, larvae and adult zebrafish (Danio rerio). The distribution of zebrafish fabp10a and fabp10b transcripts show a strikingly different tissue-specific pattern leading us to suggest that the zebrafish fabp10 duplicates had been retained in the genome owing to neofunctionalization. In another experiment to test the DDC model, transcriptional regulation of duplicated fabp genes was analyzed in zebrafish fed clofibrate, a peroxisome proliferator-activated receptor (PPAR) agonist. Clofibrate increased the steady-state level of both the duplicated copies of fabp1a/fabp1b.1, and fabp7a/fabp7b mRNA and heteronuclear RNA (hnRNA), but in different tissues of zebrafish. The steady-state level of fabp10a and fabp11a mRNA and hnRNA was elevated in liver of zebrafish, but not for fabp10b and fabp11b. We also investigated the effect of dietary fatty acids (FAs) and clofibrate on the transcriptional regulation of single copy fabp genes, fabp2, fabp3 and fabp6 in zebrafish. The steady-state level of fabp2 transcripts increased in intestine, while fabp3 mRNA increased in liver of zebrafish fed diets differing in FA content. In zebrafish fed clofibrate, fabp3 mRNA in intestine, and fabp6 mRNA in intestine and heart, was elevated. Whether the regulation of fabp gene transcription by clofibrate is controlled either directly or indirectly, the regulatory elements in the zebrafish fabp genes have diverged markedly since the WGD event, thereby supporting the DDC model.

No description available.

In eukaryotic genomes ranging from plants to mammals, DNA methylation is a major epigenetic modification of DNA by adding a methyl group exclusively to cytosine residuals. In mammalian genomes such as humans, these cytosine bases are usually followed by guanine. Although it does not change the primary DNA sequence, this covalent modification plays critical roles in several regulatory processes and can impact gene activity in a heritable fashion. What is more important, DNA methylation is essential for mammalian embryonic development and aberrant DNA methylation is implicated in several human diseases, in particular in neuro-developmental syndromes (such as the fragile X and Rett syndromes) and cancer. These biological significances disclose the importance of understanding genomic patterns and function role of DNA methylation in human, as a initial step to get to know the epigenotype and its manner in connecting the phenotype and genotype. Two key papers back in 1975 independently suggested that methylation of CpG dinucleotides in vertebrates could be established de novo and inherited through somatic cell divisions by protein machineries of DNA methyltransferases that recognizes hemi-methylated CpG palindromes. They also indicated that the methyl group could be recognized by DNA-binding proteins and that DNA methylation directly silences gene expression. After almost four decades, several key points in these foundation papers are proved to be true. Take the mammalian genome for example, there are several findings indicating the epigenetic repression of gene expression by DNA methylation. These include X-chromosome inactivation, gene imprinting and suppressing the proliferation of transposable elements and repeat elements of viral or retroviral origin. In addition to these, many novel roles of DNA methylation have also been revealed. For example, DNA methylation can regulate alternative splicing by preventing CTCF, an evolutionarily conserved zinc-finger protein, binding to DNA. By using the technique of fluorescence resonance energy transfer (FRET) and fluorescence polarization, DNA methylation has also been shown to increase nucleosome compaction through DNA-histone contacts. What is more important, DNA methylation is essential for mammalian embryonic development and aberrant change of DNA methylation has been related to disease such as cancer. However, it is also notable there are several lines of evidence contradicting the relationship between DNA methylation and gene silencing. For example, comparison of DNA methylation levels in human genome on the active and inactive X chromosomes showed reduced methylation specifically over gene bodies on inactive X chromosomes. Not only in human, DNA methylation is found to be usually targeted to the transcription units of actively transcribed genes in invertebrate species. These results prove that the function of DNA methylation is challenging to be unravel. Besides, due to the development of sequencing technique, whole genome DNA methylation profiles have been detected in diverse species. Comparing genomic patterns of DNA methylation shows considerable variation among taxa, especially between vertebrates and invertebrates. However, even though extensive studies reveal the patterns and functions of DNA methylation in different species, in the mean time, they also highlight the limits to our understanding of this complex epigenetic system. During my Ph.D., in order to perform in-depth studies of DNA methylation in diverse animals as a way to understand the complexity of DNA methylation and its functions, I dedicated my efforts in investigating and analyzing the DNA methylation profiles in diverse species, ranging from insects to primates, including both model and non-model organisms. This dissertation, which constitutes an important part of my research, mainly focuses on the DNA methylation profile in primates including human and chimpanzee. In general, I will use three chapters to elucidate my work in generating and interpreting the whole genome DNA methylation data. Firstly, we generated nucleotide-resolution whole-genome methylation maps of the prefrontal cortex of multiple humans and chimpanzees, then comprehensive comparative studies for these DNA methylation maps have been performed, by integrating data on gene expression as well. This work demonstrates that differential DNA methylation might be an important molecular mechanism driving gene-expression divergence between human and chimpanzee brains and also potentially contribute to the human-specific traits, such as evolution of disease vulnerabilities. Secondly , we performed global analyses of CpG islands (CGIs) methylation across multiple methylomes of distinctive cellular origins in human. The results from this work show that the human CpG islands can be distinctly classified into different clusters solely based upon the DNA methylation profiles, and these CpG islands clusters reflect their distinctive nature at many biological levels, including both genomic characteristics and evolutionary features. Moreover, these CpG islands clusters are non-randomly associated with several important biological phenomena and processes such as diseases, aging, and gene imprinting. These new findings shed lights in deciphering the regulatory mechanisms of CpG islands in human health and diseases. At last, by utilizing the DNA methylome from human sperm and genetic map generated from the International HapMap Consortium project, we investigated the hypothesis suggesting a potential role of germ line DNA methylation in affecting meiotic recombination, which is essential for successful meiosis and various evolutionary processes. Even thought the results imply that DNA methylation is a important factor affecting regional recombination rate, the strength of correlation between these two is not as strong as the previous report. Besides, high-throughput analyses indicate that other epigenetic modifications, tri-methylation of histone 3 lysine 4 and histone 3 lysine 27 are also global features at the recombination hotspots, and may interact with methylation to affect the recombination pattern simultaneously. This work suggests epigenetic mechanisms as additional factors affecting recombination, which cannot be fully explained by the DNA sequence itself. In summary, I hope the results from these work can expand our knowledge regarding the common and variable patterns of DNA methylation in different taxa, and shed light about the function role and its major change during animal evolution.

Degenerative suspensory ligament desmitis (DSLD) is a debilitating disease of connective tissues seen in many breeds but has become prevalent in the Peruvian Pasohorse. DSLD is believed to be a genetic disorder caused by one primary founder and most likely has a recessive mode of inheritance although a dominant or co-dominant mode of inheritance has not been ruled out. A genome scan using 259 microsatellite markers was used to test for linkage disequilibrium between one or more markers and DSLD. Two groups of Peruvian Pasohorses were selected from one population including the US and Canada. The only difference between the two groups of horses besides the size of the two groups was the presence of DSLD in the affected group and the absence of DSLD in the unaffected group. It was assumed that differences seen between the two groups in homozygosity and or common allele frequency could be an indication of linkage to DSLD. As a connective tissue disorder, there were a large number of candidate genes forDSLD to consider, yet no identical human or animal model exists. The genome scan identified five chromosomal regions where statistically significant differences were seen between affected and unaffected sample populations that could be indications of linkage to DSLD. Those chromosomes were: ECA 6, 7, 11, 14, and 26. Sequencing of a portion of the G domain in the Chondroitin Sulfate Proteoglycan2 (CSPG2) gene has mostly ruled out that segment of chromosome 14 as having linkage to DSLD. Further research needs to be conducted in the regions of ECA 6,7,11 and 26 where statistically significant differences were seen between the affected and unaffected groups, especially on ECA 6 and 11 since possible candidate genes are located in those regions based on the human comparative map.

The pathogen Campylobacter fetus (CF) is classified into two subspecies, Campylobacter fetus subspecies fetus (CFF) and Campylobacter fetus subspecies venerealis (CFV). Even though CFF and CFV are genetically closely related, they exhibit differences in their host adaptation; CFF inhabits the gastrointestinal tract of both humans and several animal species, while classical CFV is specific to the bovine genital tract and is of particular concern with respect to international bovine trade regulation. Traditionally, differentiation between the two subspecies has been achieved using a limited number of biochemical tests but more rapid and definitive genetic methods of discrimination are desired. A recent study suggested that the presence of a genomic island only in CFV could discriminate between the two sub- species but this hypothesis could not be confirmed on a collection of isolates originating in Canada. To identify alternative gene targets that would support accurate subspecies discrimination, this study has applied several approaches including suppression subtractive hybridization and whole genome sequencing supplemented with optical mapping. A subtractive hybridization screen, using a well-characterized CFV isolate recovered during routine screening of bulls in an Artificial Insemination center in western Canada and that lacked much of the genomic island and a typical Canadian CFF isolate, yielded 50 clones; characterization of these clones by hybridization screening against selected CF isolates and by nucleotide sequence BLAST analysis identified three potentially CFV-specific clones that contained inserts originating from a second genomic island. Further screening using a larger CF sample set found that only Clone #35 was truly CFV-specific. Optical maps (NcoI digest) of the Canadian CFF and CFV isolates used for the subtractive hybridization showed that certain regions of these genomes were quite distinct from those of two reference strains. Whole genome sequencing of these two isolates identified two target genes (PICFV5_ORF548 and CFF_Feature #3) that appear to be selectively retained in the two subspecies. Screening of a collection of CF isolates by PCRs targeting these three loci (SSH_Clone #35, PICFV5_ORF548 and CFF_Feature #3) supported their use for subspecies discrimination. This work demonstrates the complex genomic diversity associated with these CF subtypes and the challenge posed by their discrimination using limited genetic loci.

Breast cancer is a heterogeneous, polygenic disease and is influenced by genetic, environmental and life-style factors. Many single nucleotide polymorphisms (SNPs) associated with breast cancer risk have been identified in genome-wide association studies (GWASs) by several research groups for different populations. However, the variants identified so far contribute to a small proportion of disease risk. The objectives of the work described in this thesis were (i) to seek relevance/replicability of reported risk alleles from SNP scans to our study population; and (ii) to perform an independent GWAS for identification of additional/novel polymorphisms in the Albertan population. We approached these two end points by using cases and controls recruited in Alberta (total sample size, n=3064) in a two-stage association study (discovery study followed by replication study). We reproduced 14 of the 28 variants reported by others and also identified seven novel variants associated with breast cancer risk in our study population.

Recent proceedings in biotechnology have enabled completely new avenues in life science research to be explored. By allowing increased parallelization an ever-increasing complexity of cell samples or experiments can be investigated in shorter time and at a lower cost. This facilitates for example large-scale efforts to study cell heterogeneity at the single cell level, by analyzing cells in parallel that also can include global genomic analyses. The work presented in this thesis focuses on massively parallel analysis of cells or nucleic acid samples, demonstrating technology developments in the field as well as use of the technology in life sciences. In stem cell research issues such as cell morphology, cell differentiation and effects of reprogramming factors are frequently studied, and to obtain information on cell heterogeneity these experiments are preferably carried out on single cells. In paper I we used a high-density microwell device in silicon and glass for culturing and screening of stem cells. Maintained pluripotency in stem cells from human and mouse was demonstrated in a screening assay by antibody staining and the chip was furthermore used for studying neural differentiation. The chip format allows for low sample volumes and rapid high-throughput analysis of single cells, and is compatible with Fluorescence Activated Cell Sorting (FACS) for precise cell selection. Massively parallel DNA sequencing is revolutionizing genomics research throughout the life sciences by constantly producing increasing amounts of data from one sequencing run. However, the reagent costs and labor requirements in current massively parallel sequencing protocols are still substantial. In paper II-IV we have focused on flow-sorting techniques for improved sample preparation in bead-based massive sequencing platforms, with the aim of increasing the amount of quality data output, as demonstrated on the Roche/454 platform. In paper II we demonstrate a rapid alternative to the existing shotgun sample titration protocol and also use flow-sorting to enrich for beads that carry amplified template DNA after emulsion PCR, thus obtaining pure samples and with no downstream sacrifice of DNA sequencing quality. This should be seen in comparison to the standard 454-enrichment protocol, which gives rise to varying degrees of sample purity, thus affecting the sequence data output of the sequencing run. Massively parallel sequencing is also useful for deep sequencing of specific PCR-amplified targets in parallel. However, unspecific product formation is a common problem in amplicon sequencing and since these shorter products may be difficult to fully remove by standard procedures such as gel purification, and their presence inevitably reduces the number of target sequence reads that can be obtained in each sequencing run. In paper III a gene-specific fluorescent probe was used for target-specific FACS enrichment to specifically enrich for beads with an amplified target gene on the surface. Through this procedure a nearly three-fold increase in fraction of informative sequences was obtained and with no sequence bias introduced. Barcode labeling of different DNA libraries prior to pooling and emulsion PCR is standard procedure to maximize the number of experiments that can be run in one sequencing lane, while also decreasing the impact of technical noise. However, variation between libraries in quality and GC content affects amplification efficiency, which may result in biased fractions of the different libraries in the sequencing data. In paper IV barcode specific labeling and flow-sorting for normalization of beads with different barcodes on the surface was used in order to weigh the proportion of data obtained from different samples, while also removing mixed beads, and beads with no or poorly amplified product on the surface, hence also resulting in an increased sequence quality. In paper V, cell heterogeneity within a human being is being investigated by low-coverage whole genome sequencing of single cell material. By focusing on the most variable portion of the human genome, polyguanine nucleotide repeat regions, variability between different cells is investigated and highly variable polyguanine repeat loci are identified. By selectively amplifying and sequencing polyguanine nucleotide repeats from single cells for which the phylogenetic relationship is known, we demonstrate that massively parallel sequencing can be used to study cell-cell variation in length of these repeats, based on which a phylogenetic tree can be drawn. / QC 20111031

Breast cancer is the most common cancer in women worldwide. However, the genetic alterations that lead to breast cancer are not fully understood. This thesis aims to identify novel genes of potential mechanistic, diagnostic or therapeutic interest in breast cancers by mutational analysis and whole-genome sequencing. In paper I, sequencing of 36 previously identified candidate genes in 96 breast tumors with patient-matched normal DNA determined the somatic mutation prevalence of these candidate genes and identified additional mutations in Notch, NF-κB, PI3K, and Hedgehog pathways as well as in processes mediating DNA methylation, RNA processing and calcium signaling. In paper II, comparison of massively parallel mate-pair sequencing results of a human genome before and after phi29-mediated multiple displacement amplification (MDA) revealed that MDA introduces structural alteration artifacts, with an emphasis on false positive inversions, and impairs the sensitivity to detect true inversions. Therefore, MDA has limited value in sample preparation for whole-genome sequencing for structural alteration detection. In paper III, massively parallel paired-end sequencing identified gene rearrangements in 15 hormone receptor negative breast cancers. Forty validated rearrangements were predicted to directly affect 30 genes, involved in epigenetic regulation, cell mitosis, signalling transduction and glycolytic flux. RNA interference-based assays revealed the potential roles in cell growth of some affected genes, among which DDX10 was implicated to be involved in apoptosis. In paper IV, a method for statistical evaluation of putative translocations detected by massively parallel paired-end sequencing was proposed. In an application of this method to analyse translocations detected by cancer genome deep paired-end sequencing, 76 putative translocations were classified into four categories, with the majority likely to be caused by mismapping due to repetitive regions. Taken together, this thesis provides insights into genes and pathways mutated in sporadic breast cancer genomes, which broaden our understanding of the genetic basis of breast cancer and may ultimately facilitate the diagnosis and treatment of this disease.