Phylogenetic structural modeling of molecular evolution

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Mitochondrial Inheritance, Mito-nuclear Coevolution, and Sex-associated Genes in Bivalve Molluscs: Transcriptomics and Molecular Evolution

Bivalvia represents an ancient taxon including around 25,000 living species that have adapted to a wide range of environmental conditions, and show a great diversity in body size, shell shapes, and anatomic structure. Bivalves are characterized by highly variable genome sizes and extremely high levels of heterozygosity, which obstacle complete and accurate genome assemblies and hinder further genomic studies. Moreover, some bivalve species presented a stable evolutionary exception to the strictly maternal inheritance of mitochondria, namely doubly uniparental inheritance (DUI), making these species a precious model to study mitochondrial biology. During my PhD, I focused on a DUI species, the Manila clam Ruditapes philippinarum, and my work was two-folded. First, taking advantage of a newly assembled draft genome and a large RNA-seq dataset from different tissues of both sexes, I investigated 1) the role of gene expression and alternative splicing in tissue differentiation; 2) the relationship across tissue specificity, regulatory network connectivity, and sequence evolution; 3) sexual contrasting genetic markers potentially associated with sexual differentiation. The detailed information for this part is in Chapter 2. Second, using the same RNA-seq data, I investigated how nuclear oxidative phosphorylation (OXPHOS) genes coordinate with two divergent mitochondrial genomes in DUI species (mito-nuclear coordination and coevolution). To address this question, I compared transcription, polymorphism, and synonymous codon usage in the mitochondrial and nuclear OXPHOS genes of R. philippinarum in Chapter 3. To my knowledge, this thesis represents the first study exploring the role of alternative splicing in tissue differentiation, and the first study analyzing both transcriptional regulation and sequence evolution to investigate the coordination of OXPHOS genes in bivalves.

Evolutionary history of Ramphastos toucans: Molecular phylogenetics, temporal diversification, and biogeography

The toucan genus Ramphastos (Piciformes: Ramphastidae) has been a model in the formulation of Neotropical paleobiogeographic hypotheses. Weckstein (2005) reported on the phylogenetic history of this genus based on three mitochondrial genes, but some relationships were weakly supported and one of the subspecies of R. vitellinus (citreolaemus) was unsampled. This study expands on Weckstein (2005) by adding more DNA sequence data (including a nuclear marker) and more samples, including R v. citreolaemus. Maximum parsimony, maximum likelihood, and Bayesian methods recovered similar trees, with nodes showing high support. A monophyletic R. vitellinus complex was strongly supported as the sister-group to R. brevis. The results also confirmed that the southeastern and northern populations of R. vitellinus ariel are paraphyletic. X v. citreolaemus is sister to the Amazonian subspecies of the vitellinus complex. Using three protein-coding genes (COI, cytochrome-b and ND2) and interval-calibrated nodes under a Bayesian relaxed-clock framework, we infer that ramphastid genera originated in the middle Miocene to early Pliocene, Ramphastos species originated between late Miocene and early Pleistocene, and intra-specific divergences took place throughout the Pleistocene. Parsimony-based reconstruction of ancestral areas indicated that evolution of the four trans-Andean Ramphastos taxa (R. v. citreolaemus, R. a. swainsonii, R. brevis and R. sulfuratus) was associated with four independent dispersals from the cis-Andean region. The last pulse of Andean uplift may have been important for the evolution of R. sulfuratus, whereas the origin of the other trans-Andean Ramphastos taxa is consistent with vicariance due to drying events in the lowland forests north of the Andes. Estimated rates of molecular evolution were higher than the ""standard"" bird rate of 2% substitutions/site/million years for two of the three genes analyzed (cytochrome-b and ND2). (C) 2009 Elsevier Inc. All rights reserved.

Linking life-history traits, ecology, and niche breadth evolution in north american eriogonoids (polygonaceae).

Abstract Macroevolutionary and microevolutionary studies provide complementary explanations of the processes shaping the evolution of niche breadth. Macroevolutionary approaches scrutinize factors such as the temporal and spatial environmental heterogeneities that drive differentiation among species. Microevolutionary studies, in contrast, focus on the processes that affect intraspecific variability. We combine these perspectives by using macroevolutionary models in a comparative study of intraspecific variability. We address potential differences in rates of evolution of niche breadth and position in annual and perennial plants of the Eriogonoideae subfamily of the Polygonaceae. We anticipated higher rates of evolution in annuals than in perennials owing to differences in generation time that are paralleled by rates of molecular evolution. Instead, we found that perennial eriogonoid species present greater environmental tolerance (wider climate niche) than annual species. Niche breadth of perennial species has evolved two to four times faster than in annuals, while niche optimum has diversified more rapidly among annual species than among perennials. Niche breadth and average elevation of species are correlated. Moreover, niche breadth increases more rapidly with mean species elevation in perennials than in annuals. Our results suggest that both environmental gradients and life-history strategy influence rates and patterns of niche breadth evolution.

A phylogenetic study of the phosphoenolpyruvate carboxylase multigene family in Poaceae: understanding the molecular changes linked to C4 photosynthesis evolution

Phylogenetic reconstructions have supported several independent appearances of C₄ photosynthesis within grasses (Poaceae). These recurrent appearances appear to contradict the large number of biochemical and morphological changes required to change from C₃ to C₄, a paradox that leads to questions about the genetic changes underlying C₄ evolution. In this study, we analysed sequences encoding phosphoenolpyruvate carboxylases (PEPCs) in grasses in order to gain insights into the origin of the ppc-C₄ gene, which encodes a key enzyme in the C₄ pathway. We screened databanks for PEPC genes or cDNAs in grasses. A coding sequence of 1130 base pairs was used to build phylogenetic trees that supported the existence of four distinct PEPC gene lineages. Ppc-C₄ present in all C₄ grasses was also found in two C₃ species. The ppc-C₄ clade was congruent with the species tree, suggesting orthologous evolution. This result would imply that ppc-C₄ appeared without any duplication event. Nevertheless, caution is needed since the sampling of our study is still far from comprehensive. Further investigation with an increased sampling is recommended to elucidate the evolutionary changes underlying ppc-C₄ gene evolution in grasses.

Temporal evolution of sulfated zirconia sol-gel during dip coating: Analysis of mass drainage with time-varying refractive index

The formation of sulfated zirconia films from a sol-gel derived aqueous suspension is subjected to double-optical monitoring during batch dip coating. Interpretation of interferometric patterns, previously obscured by a variable refractive index, is now made possible by addition of its direct measurement by a polarimetric technique in real time. Significant sensitivity of the resulting physical thickness and refractive index curves (uncertainties of ±7 nm and ±0.005, respectively) to temporal film evolution is shown under different withdrawal speeds. As a first contribution to quantitative understanding of temporal film formation with varying nanostructure during dip coating, detailed analysis is directed to the stage of the process dominated by mass drainage, whose simple modeling with temporal t-1/2 dependence is verified experimentally. © 2006 Elsevier B.V. All rights reserved.

Molecular phylogeny analysis of fiddler crabs: test of the hypothesis of increasing behavioral complexity in evolution.

The current phylogenetic hypothesis for the evolution and biogeography of fiddler crabs relies on the assumption that complex behavioral traits are assumed to also be evolutionary derived. Indo-west Pacific fiddler crabs have simpler reproductive social behavior and are more marine and were thought to be ancestral to the more behaviorally complex and more terrestrial American species. It was also hypothesized that the evolution of more complex social and reproductive behavior was associated with the colonization of the higher intertidal zones. Our phylogenetic analysis, based upon a set of independent molecular characters, however, demonstrates how widely entrenched ideas about evolution and biogeography led to a reasonable, but apparently incorrect, conclusion about the evolutionary trends within this pantropical group of crustaceans. Species bearing the set of "derived traits" are phylogenetically ancestral, suggesting an alternative evolutionary scenario: the evolution of reproductive behavioral complexity in fiddler crabs may have arisen multiple times during their evolution. The evolution of behavioral complexity may have arisen by coopting of a series of other adaptations for high intertidal living and antipredator escape. A calibration of rates of molecular evolution from populations on either side of the Isthmus of Panama suggest a sequence divergence rate for 16S rRNA of 0.9% per million years. The divergence between the ancestral clade and derived forms is estimated to be approximately 22 million years ago, whereas the divergence between the American and Indo-west Pacific is estimated to be approximately 17 million years ago.

A timescale for diatom evolution based on four molecular markers: reassessment of ghost lineages and major steps defining diatom evolution.

ABSTRACT. – Phylogenies and molecular clocks of the diatoms have largely been inferred from SSU rDNA sequences. A new phylogeny of diatoms was estimated using four gene markers SSU and LSU rDNA rbcL and psbA (total 4352 bp) with 42 diatom species. The four gene trees analysed with a maximum likelihood (ML) and Baysian (BI) analysis recovered a monophyletic origin of the new diatom classes with high bootstrap support, which has been controversial with single gene markers using single outgroups and alignments that do not take secondary structure of the SSU gene into account. The divergence time of the classes were calculated from a ML tree in the MultliDiv Time program using a Bayesian estimation allowing for simultaneous constraints from the fossil record and varying rates of molecular evolution of different branches in the phylogenetic tree. These divergence times are generally in agreement with those proposed by other clocks using single genes with the exception that the pennates appear much earlier and suggest a longer Cretaceous fossil record that has yet to be sampled. Ghost lineages (i.e. the discrepancy between first appearance (FA) and molecular clock age of origin from an extant taxon) were revealed in the pennate lineage, whereas those ghost lineages in the centric lineages previously reported by others are reviewed and referred to earlier literature.

A timescale for diatom evolution based on four molecular markers: reassessment of ghost lineages and major steps defining diatom evolution.

ABSTRACT. – Phylogenies and molecular clocks of the diatoms have largely been inferred from SSU rDNA sequences. A new phylogeny of diatoms was estimated using four gene markers SSU and LSU rDNA rbcL and psbA (total 4352 bp) with 42 diatom species. The four gene trees analysed with a maximum likelihood (ML) and Baysian (BI) analysis recovered a monophyletic origin of the new diatom classes with high bootstrap support, which has been controversial with single gene markers using single outgroups and alignments that do not take secondary structure of the SSU gene into account. The divergence time of the classes were calculated from a ML tree in the MultliDiv Time program using a Bayesian estimation allowing for simultaneous constraints from the fossil record and varying rates of molecular evolution of different branches in the phylogenetic tree. These divergence times are generally in agreement with those proposed by other clocks using single genes with the exception that the pennates appear much earlier and suggest a longer Cretaceous fossil record that has yet to be sampled. Ghost lineages (i.e. the discrepancy between first appearance (FA) and molecular clock age of origin from an extant taxon) were revealed in the pennate lineage, whereas those ghost lineages in the centric lineages previously reported by others are reviewed and referred to earlier literature.

A total-evidence phylogeny of ticks provides insights into the evolution of life cycles and biogeography

We inferred the phylogeny of 33 species of ticks from the subfamilies Rhipicephalinae and Hyalomminae from analyses of nuclear and mitochondrial DNA and morphology. We used nucleotide sequences from 12S rRNA, cytochrome c oxidase I, internal transcribed spacer 2 of the nuclear rRNA, and 18S rRNA. Nucleotide sequences and morphology were analyzed separately and together in a total-evidence analysis. Analyses of the five partitions together (3303 characters) gave the best-resolved and the best-supported hypothesis so far for the phylogeny of ticks in the Rhipicephalinae and Hyalomminae, despite the fact that some partitions did not have data for some taxa. However, most of the hidden conflict (lower support in the total-evidence analyses compared to that in the individual analyses) was found in those partitions that had taxa without data. The partitions with complete taxonomic sampling had more hidden support (higher support in the total-evidence analyses compared to that in the separate-partition analyses) than hidden conflict. Mapping of geographic origins of ticks onto our phylogeny indicates an African origin for the Rhipicephalinae sensu lato (i.e., including Hyalomma spp.), the Rhipicephalus-Boophilus lineage, the Dermacentor-Anocentor lineage, and the Rhipicephalus-Booophilus-Nosomma-Hyalomma-Rhipicentor lineage. The Nosomma-Hyalomma lineage appears to have evolved in Asia. Our total-evidence phylogeny indicates that (i) the genus Rhipicephalus is paraphyletic with respect to the genus Boophilus, (ii) the genus Dermacentor is paraphyletic with respect to the genus Anocentor, and (iii) some subgenera of the genera Hyalomma and Rhipicephalus are paraphyletic with respect to other subgenera in these genera. Study of the Rhipicephalinae and Hyalomminae over the last 7 years has shown that analyses of individual datasets (e.g., one gene or morphology) seldom resolve many phylogenetic relationships, but analyses of more than one dataset can generate well-resolved phylogenies for these ticks. (C) 2001 Academic Press.

Testing the relationship between morphological and molecular rates of change along phylogenies

Molecular evolution has been considered to be essentially a stochastic process, little influenced by the pace of phenotypic change. This assumption was challenged by a study that demonstrated an association between rates of morphological and molecular change estimated for total-evidence phylogenies, a finding that led some researchers to challenge molecular date estimates of major evolutionary radiations. Here we show that Omland's (1997) result is probably due to methodological bias, particularly phylogenetic nonindependence, rather than being indicative of an underlying evolutionary phenomenon. We apply three new methods specifically designed to overcome phylogenetic bias to 13 published phylogenetic datasets for vertebrate taxa, each of which includes both morphological characters and DNA sequence data. We find no evidence of an association between rates of molecular and morphological rates of change.

Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda)

ITS2 sequences are used extensively in molecular taxonomy and population genetics of arthropods and other animals yet little is known about the molecular evolution of ITS2. We studied the secondary structure of ITS2 in species from each of the six main lineages of hard ticks (family Ixodidae). The ITS2 of these ticks varied in length from 679 bp in Ixodes scapularis to 1547 bp in Aponomma concolor. Nucleotide content varied also: the ITS2 of ticks from the Prostriata lineage (Ixodes spp.) had 46-49% GC whereas ITS2 sequences of ticks from the Metastriata lineage (all other hard ticks) had 61-62% GC. Despite variation in nucleotide sequence, the secondary structure of the ITS2 of all of these ticks apparently has five domains. Stems 1, 3, 4 and 5 of this secondary structure were obvious in all of the species studied. However, stem 2 was not always obvious despite the fact that it is flanked by highly conserved sequence motifs in the adjacent stems, stems 1 and 3. The ITS2 of hard ticks has apparently evolved mostly by increases and decreases in length of the nucleotide sequences, which caused increases, and decreases in the length of stems of the secondary structure. This is most obvious when stems of the secondary structures of the Prostriata (Ixodes spp.) are compared to those of the Metastriata (all other hard ticks). Increases in the size of the ITS2 may have been caused by replication slippage which generated large repeats, like those seen in Haemaphysalis humerosa and species from the Rhipicepalinae lineage, and the small repeats found in species from the other lineages of ticks.

Engineering Temporal and Spatial Aspects in OWL using Patterns

WWW is a huge, open, heterogeneous system, however its contents data is mainly human oriented. The Semantic Web needs to assure that data is readable and “understandable” to intelligent software agents, though the use of explicit and formal semantics. Ontologies constitute a privileged artifact for capturing the semantic of the WWW data. Temporal and spatial dimensions are transversal to the generality of knowledge domains and therefore are fundamental for the reasoning process of software agents. Representing temporal/spatial evolution of concepts and their relations in OWL (W3C standard for ontologies) it is not straightforward. Although proposed several strategies to tackle this problem but there is still no formal and standard approach. This work main goal consists of development of methods/tools to support the engineering of temporal and spatial aspects in intelligent systems through the use of OWL ontologies. An existing method for ontology engineering, Fonte was used as framework for the development of this work. As main contributions of this work Fonte was re-engineered in order to: i) support the spatial dimension; ii) work with OWL Ontologies; iii) and support the application of Ontology Design Patterns. Finally, the capabilities of the proposed approach were demonstrated by engineering time and space in a demo ontology about football.