Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.

The Accumulation of Reproductive Incompatibilities in African Cichlid Fish

The rate at which different components of reproductive isolation accumulate with divergence time between species has only been studied in a limited, but growing, number of species. We measured premating isolation and hybrid inviability at four different ontogenetic stages from zygotes to adults in interspecific hybrids of 26 pairs of African cichlid species, spanning the entire East African haplochromine radiation. We then used multiple relaxed molecular clock calibrations to translate genetic distances into absolute ages to compare evolutionary rates of different components of reproductive isolation. We find that premating isolation accumulates fast initially but then changes little with increasing genetic distance between species. In contrast, postmating isolation between closely related species is negligible but then accumulates rapidly, resulting in complete hybrid inviability after 4.4/8.5/18.4 million years (my). Thus, the rate at which complete intrinsic incompatibilities arise in this system is orders of magnitude lower than rates of speciation within individual lake radiations. Together these results suggest divergent ecological adaptations may prevent populations from interbreeding and help maintain cichlid species diversity, which may be vulnerable to environmental degradation. By quantifying the capacity to produce viable hybrids between allopatric, distantly related lineages our results also provide an upper divergence time limit for the "hybrid swarm origin" model of adaptive radiation.

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P.a.axillaris, P.a.parodii and P.a.subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P.axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100000years, which is compatible with a divergence time between 35000 and 107000years ago between P.a.axillaris and P.a.parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P.a.axillaris and P.a.parodii in their contact zone.

Hybrid Breakdown in Cichlid Fish

Studies from a wide diversity of taxa have shown a negative relationship between genetic compatibility and the divergence time of hybridizing genomes. Theory predicts the main breakdown of fitness to happen after the F1 hybrid generation, when heterosis subsides and recessive allelic (Dobzhansky-Muller) incompatibilities are increasingly unmasked. We measured the fitness of F2 hybrids of African haplochromine cichlid fish bred from species pairs spanning several thousand to several million years divergence time. F2 hybrids consistently showed the lowest viability compared to F1 hybrids and non-hybrid crosses (crosses within the grandparental species), in agreement with hybrid breakdown. Especially the short- and long-term survival (2 weeks to 6 months) of F2 hybrids was significantly reduced. Overall, F2 hybrids showed a fitness reduction of 21% compared to F1 hybrids, and a reduction of 43% compared to the grandparental, non-hybrid crosses. We further observed a decrease of F2 hybrid viability with the genetic distance between grandparental lineages, suggesting an important role for negative epistatic interactions in cichlid fish postzygotic isolation. The estimated time window for successful production of F2 hybrids resulting from our data is consistent with the estimated divergence time between the multiple ancestral lineages that presumably hybridized in three major adaptive radiations of African cichlids.

Hox cluster genomics in the horn shark, Heterodontus francisci

Reconstructing the evolutionary history of Hox cluster origins will lead to insights into the developmental and evolutionary significance of Hox gene clusters in vertebrate phylogeny and to their role in the origins of various vertebrate body plans. We have isolated two Hox clusters from the horn shark, Heterodontus francisci. These have been sequenced and compared with one another and with other chordate Hox clusters. The results show that one of the horn shark clusters (HoxM) is orthologous to the mammalian HoxA cluster and shows a structural similarity to the amphioxus cluster, whereas the other shark cluster (HoxN) is orthologous to the mammalian HoxD cluster based on cluster organization and a comparison with noncoding and Hox gene-coding sequences. The persistence of an identifiable HoxA cluster over an 800-million-year divergence time demonstrates that the Hox gene clusters are highly integrated and structured genetic entities. The data presented herein identify many noncoding sequence motifs conserved over 800 million years that may function as genetic control motifs essential to the developmental process.

Molecular evidence for an African origin of the Hawaiian endemic Hesperomannia (Asteraceae)

Identification of the progenitors of plants endemic to oceanic islands often is complicated by extreme morphological divergence between island and continental taxa. This is especially true for the Hawaiian Islands, which are 3,900 km from any continental source. We examine the origin of Hesperomannia, a genus of three species endemic to Hawaii that always have been placed in the tribe Mutisieae of the sunflower family. Phylogenetic analyses of representatives from all tribes in this family using the chloroplast gene ndhF (where ndhF is the ND5 protein of chloroplast NADH dehydrogenase) indicate that Hesperomannia belongs to the tribe Vernonieae. Phylogenetic comparisons within the Vernonieae using sequences of both ndhF and the internal transcribed spacer regions of nuclear ribosomal DNA reveal that Hesperomannia is sister to African species of Vernonia. Long-distance dispersal northeastward from Africa to southeast Asia and across the many Pacific Ocean island chains is the most likely explanation for this unusual biogeographic connection. The 17- to 26-million-year divergence time between African Vernonia and Hesperomannia estimated by the DNA sequences predates the age of the eight existing Hawaiian Islands. These estimates are consistent with an hypothesis that the progenitor of Hesperomannia arrived at one of the low islands of the Hawaiian-Emperor chain between the late Oligocene and mid-Miocene when these islands were above sea level. Subsequent to its arrival the southeast Pacific island chains served as steppingstones for dispersal to the existing Hawaiian Islands.

Diversificação de espécies e da morfologia em serpentes da família Viperidae: padrões e processos

A diversidade de espécies e fenotípica pode variar consideravelmente entre grupos taxonômicos e ao longo do tempo em uma mesma linhagem. O estudo de tais variações tornou-se um dos principais objetivos da biologia evolutiva fornecendo informações importantes a respeito dos possíveis mecanismos que regulam a biodiversidade. Dessa forma, o objetivo geral da presente tese foi investigar os padrões da diversificação de espécies e da morfologia em um grupo cosmopolita de serpentes, a família Viperidae, e os potenciais processos subjacentes. Primeiramente, (1) reconstruímos as relações filogenéticas e estimamos os tempos de divergência entre as linhagens da família Viperidae utilizando uma abordagem Bayesiana. (2) Aplicando um método recentemente desenvolvido (BAMM), exploramos como as taxas de especiação e extinção variaram ao longo da radiação do grupo inferindo os possíveis processos reguladores. Por fim, (3) analisamos se a evolução do tamanho do corpo e as taxas de especiação variam nos diferentes habitats ocupados pelos viperídeos (terrestres vs arborícola). Nesta tese geramos a filogenia molecular de viperídeos mais completa até o momento utilizando sequências para 11 genes mitocondriais e nucleares abrangendo 79% das espécies viventes (264 terminais) e todos com exceção de um gênero. De maneira geral, foi possível obter relações filogenéticas robustas para o grupo com a maioria dos gêneros sendo monofilética. Os tempos de divergência obtidos indicam que os viperídeos começaram a diversificar em meados do Paleoceno tardio/meio do Eoceno inferindo idades um pouco mais tardias que o encontrado em estudos anteriores. Durante a radiação do grupo, um aumento nas taxas de especiação parece ter ocorrido durante a diversificação dos crotalíneos (pit vipers) em decorrência não só da evolução das fossetas loreais mas também como resultado de mudanças geológicas e climáticas na Ásia e da invasão do novo mundo. Após este rápido aumento inicial, as taxas de especiação desaceleraram em direção ao presente. Por fim, os resultados aqui apresentados indicam que apesar dos habitats arborícolas limitarem a evolução morfológica nos viperídeos, a evolução da arborealidade parece não afetar as taxas de especiação que permanecem similares entre linhagens arborícolas e terrestres. Isto sugere dois cenários: (1) a especiação acontece de forma independente das mudanças morfológicas nos viperídeos; ou (2) o isolamento geográfico seria um mecanismo importante na diversificação de linhagens arborícolas contrabalançando decréscimos nas oportunidades de especiação possivelmente relacionados às pressões seletivas impostas pelo ambiente arborícola. A presente tese contribui para entendermos mais sobre como evoluíram os viperídeos ao longo dos seus ∼50 milhões de anos. Além de propor cenários e hipóteses a serem futuramente explorados com os viperídeos, elaboramos uma discussão ampla e conceitual a respeito dos possíveis mecanismos por trás da diversificação de espécies e da morfologia que poderiam também ser contemplados para outros grupos de organismos. Portanto, a presente tese contribui não só para entendermos os mecanismos que geram e mantém a diversidade de serpentes, mas também para enriquecer a discussão dos mecanismos que geram e mantém a biodiversidade como um todo

A comprehensive and integrative reconstruction of evolutionary history for Anomura (Crustacea: Decapoda)

Background The infraorder Anomura has long captivated the attention of evolutionary biologists due to its impressive morphological diversity and ecological adaptations. To date, 2500 extant species have been described but phylogenetic relationships at high taxonomic levels remain unresolved. Here, we reconstruct the evolutionary history—phylogeny, divergence times, character evolution and diversification—of this speciose clade. For this purpose, we sequenced two mitochondrial (16S and 12S) and three nuclear (H3, 18S and 28S) markers for 19 of the 20 extant families, using traditional Sanger and next-generation 454 sequencing methods. Molecular data were combined with 156 morphological characters in order to estimate the largest anomuran phylogeny to date. The anomuran fossil record allowed us to incorporate 31 fossils for divergence time analyses. Results Our best phylogenetic hypothesis (morphological + molecular data) supports most anomuran superfamilies and families as monophyletic. However, three families and eleven genera are recovered as para- and polyphyletic. Divergence time analysis dates the origin of Anomura to the Late Permian ~259 (224–296) MYA with many of the present day families radiating during the Jurassic and Early Cretaceous. Ancestral state reconstruction suggests that carcinization occurred independently 3 times within the group. The invasion of freshwater and terrestrial environments both occurred between the Late Cretaceous and Tertiary. Diversification analyses found the speciation rate to be low across Anomura, and we identify 2 major changes in the tempo of diversification; the most significant at the base of a clade that includes the squat-lobster family Chirostylidae. Conclusions Our findings are compared against current classifications and previous hypotheses of anomuran relationships. Many families and genera appear to be poly- or paraphyletic suggesting a need for further taxonomic revisions at these levels. A divergence time analysis provides key insights into the origins of major lineages and events and the timing of morphological (body form) and ecological (habitat) transitions. Living anomuran biodiversity is the product of 2 major changes in the tempo of diversification; our initial insights suggest that the acquisition of a crab-like form did not act as a key innovation.

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.

CHARACTERIZATION OF EXPONENTIAL DIVERGENCE OF THE KALMAN FILTER FOR TIME-VARYING SYSTEMS

This paper studies semistability of the recursive Kalman filter in the context of linear time-varying (LTV), possibly nondetectable systems with incorrect noise information. Semistability is a key property, as it ensures that the actual estimation error does not diverge exponentially. We explore structural properties of the filter to obtain a necessary and sufficient condition for the filter to be semistable. The condition does not involve limiting gains nor the solution of Riccati equations, as they can be difficult to obtain numerically and may not exist. We also compare semistability with the notions of stability and stability w.r.t. the initial error covariance, and we show that semistability in a sense makes no distinction between persistent and nonpersistent incorrect noise models, as opposed to stability. In the linear time invariant scenario we obtain algebraic, easy to test conditions for semistability and stability, which complement results available in the context of detectable systems. Illustrative examples are included.

Measuring graph similarity through continuous-time quantum walks and the quantum Jensen-Shannon divergence

In this paper we propose a quantum algorithm to measure the similarity between a pair of unattributed graphs. We design an experiment where the two graphs are merged by establishing a complete set of connections between their nodes and the resulting structure is probed through the evolution of continuous-time quantum walks. In order to analyze the behavior of the walks without causing wave function collapse, we base our analysis on the recently introduced quantum Jensen-Shannon divergence. In particular, we show that the divergence between the evolution of two suitably initialized quantum walks over this structure is maximum when the original pair of graphs is isomorphic. We also prove that under special conditions the divergence is minimum when the sets of eigenvalues of the Hamiltonians associated with the two original graphs have an empty intersection.

Lineage divergence detected in the malaria vector Anopheles marajoara (Diptera: Culicidae) in Amazonian Brazil

Background: Cryptic species complexes are common among anophelines. Previous phylogenetic analysis based on the complete mtDNA COI gene sequences detected paraphyly in the Neotropical malaria vector Anopheles marajoara. The ""Folmer region"" detects a single taxon using a 3% divergence threshold. Methods: To test the paraphyletic hypothesis and examine the utility of the Folmer region, genealogical trees based on a concatenated (white + 3' COI sequences) dataset and pairwise differentiation of COI fragments were examined. The population structure and demographic history were based on partial COI sequences for 294 individuals from 14 localities in Amazonian Brazil. 109 individuals from 12 localities were sequenced for the nDNA white gene, and 57 individuals from 11 localities were sequenced for the ribosomal DNA (rDNA) internal transcribed spacer 2 (ITS2). Results: Distinct A. marajoara lineages were detected by combined genealogical analysis and were also supported among COI haplotypes using a median joining network and AMOVA, with time since divergence during the Pleistocene (< 100,000 ya). COI sequences at the 3' end were more variable, demonstrating significant pairwise differentiation (3.82%) compared to the more moderate 2.92% detected by the Folmer region. Lineage 1 was present in all localities, whereas lineage 2 was restricted mainly to the west. Mismatch distributions for both lineages were bimodal, likely due to multiple colonization events and spatial expansion (similar to 798 - 81,045 ya). There appears to be gene flow within, not between lineages, and a partial barrier was detected near Rio Jari in Amapa state, separating western and eastern populations. In contrast, both nDNA data sets (white gene sequences with or without the retention of the 4th intron, and ITS2 sequences and length) detected a single A. marajoara lineage. Conclusions: Strong support for combined data with significant differentiation detected in the COI and absent in the nDNA suggest that the divergence is recent, and detectable only by the faster evolving mtDNA. A within subgenus threshold of >2% may be more appropriate among sister taxa in cryptic anopheline complexes than the standard 3%. Differences in demographic history and climatic changes may have contributed to mtDNA lineage divergence in A. marajoara.

Phylogeography of the cactophilic species Drosophila gouveai: demographic events and divergence timing in dry vegetation enclaves in eastern Brazil

Aim The aim of this study was to assess the causal mechanisms underlying populational subdivision in Drosophila gouveai, a cactophilic species associated with xeric vegetation enclaves in eastern Brazil. A secondary aim was to investigate the genetic effects of Pleistocene climatic fluctuations on these environments. Location Dry vegetation enclaves within the limits of the Cerrado domain in eastern Brazil. Methods We determined the mitochondrial DNA haplotypes of 55 individuals (representing 12 populations) based on sequence data of a 483-bp fragment from the cytochrome c oxidase subunit II (COII) gene. Phylogenetic and coalescent analyses were used to test for the occurrence of demographic events and to infer the time of divergence amongst genetically independent groups. Results Our analyses revealed the existence of two divergent subclades (G1 and G2) plus an introgressed clade restricted to the southernmost range of D. gouveai. Subclades G1 and G2 displayed genetic footprints of range expansion and segregated geographical distributions in south-eastern and some central highland regions, east and west of the Parana River valley. Molecular dating indicated that the main demographic and diversification events occurred in the late to middle Pleistocene. Main conclusions The phylogeographical and genetic patterns observed for D. gouveai in this study are consistent with changes in the distribution of dry vegetation in eastern Brazil. All of the estimates obtained by molecular dating indicate that range expansion and isolation pre-dated the Last Glacial Maximum, occurring during the late to middle Pleistocene, and were probably triggered by climatic changes during the Pleistocene. The current patchy geographical distribution and population subdivision in D. gouveai is apparently closely linked to these past events.

Microevolution in island forms: the roles of drift and directional selection in morphological divergence of a passerine bird

Theory predicts that in small isolated populations random genetic drift can lead to phenotypic divergence; however this prediction has rarely been tested quantitatively in natural populations. Here we utilize natural repeated island colonization events by members of the avian species complex, Zosterops lateralis, to assess whether or not genetic drift alone is an adequate explanation for the observed patterns of microevolutionary divergence in morphology. Morphological and molecular genetic characteristics of island and mainland populations are compared to test three predictions of drift theory: (1) that the pattern of morphological change is idiosyncratic to each island; (2) that there is concordance between morphological and neutral genetic shifts across island populations; and (3) for populations whose time of colonization is known, that the rate of morphological change is sufficiently slow to be accounted for solely by genetic drift. Our results are not consistent with these predictions. First, the direction of size shifts was consistently towards larger size, suggesting the action of a nonrandom process. Second, patterns of morphological divergence among recently colonized populations showed little concordance with divergence in neutral genetic characters. Third, rate tests of morphological change showed that effective population sizes were not small enough for random processes alone to account for the magnitude of microevolutionary change. Altogether, these three lines of evidence suggest that drift alone is not an adequate explanation of morphological differentiation in recently colonized island Zosterops and therefore we suggest that the observed microevolutionary changes are largely a result of directional natural selection.