Evolutionary models for insertions and deletions in a probabilistic modeling framework

Background  

Probabilistic models for sequence comparison (such as hidden Markov models and pair hidden Markov models for proteins and mRNAs, or their context-free grammar counterparts for structural RNAs) often assume a fixed degree of divergence. Ideally we would like these models to be conditional on evolutionary divergence time.     

Trypanosoma (Schizotrypanum) cruzi: repetitive DNA sequence evolution in three geographically distinct isolates

1. Middle-repetitive DNA sequences were analyzed by molecular hybridization to determine both the extent of complementarity and time of evolutionary divergence between isolates of Trypanosoma cruzi from Argentina, Mexico, or Venezuela. 2. Although molecular hybridizations showed no significant difference between the middle-repetitive DNAs of the Mexican and Venezuelan isolates, there was a 2.7% base pair mismatch in hybrid molecules formed by association of strands from both the Mexican and Argentine isolates. 3. Using the rates for divergence of middle-repetitive DNA in sea urchins, the Mexican and Argentine isolates were estimated to have diverged approximately 20-25 million years ago. 4. Analysis of the kinetics for the DNA hybridizations indicates that only minor amplifications of specific gene sequences or changes in the complexity of the genomes could have occurred during the divergence of the three isolates studied.     

Phylogeny of Organisms Investigated by the Base-Pair Changes in the Stem Regions of Small and Large Ribosomal Subunit RNAs

In order to obtain the evolutionary distance data that are as purely additive as possible, we have developed a novel method for evaluating the evolutionary distances from the base-pair changes in stem regions of ribosomal RNAs (rRNAs). The application of this method to small-subunit (SSU) and large-subunit (LSU) rRNAs provides the distance data, with which both the unweighted pair group method of analysis and the neighbor-joining method give almost the same tree topology of most organisms except for some Protoctista, thermophilic bacteria, parasitic organisms, and endosymbionts. Although the evolutionary distances calculated with LSU rRNAs are somewhat longer than those with SSU rRNAs, the difference, probably due to a slight difference in functional constraint, is substantially decreased when the distances are converted into the divergence times of organisms by the measure of the time scale estimated in each type of rRNAs. The divergence times of main branches agree fairly well with the geological record of organisms, at least after the appearance of oxygen-releasing photosynthesis, although the divergence times of Eukaryota, Archaebacteria, and Eubacteria are somewhat overestimated in comparison with the geological record of Earth formation. This result is explained by considering that the mutation rate is determined by the accumulation of misrepairs for DNA damage caused by radiation and that the effect of radiation had been stronger before the oxygen molecules became abundant in the atmosphere of the Earth. Received: 23 October 1997 / Accepted: 12 August 1998     

Early diversification of sperm size in the evolutionary history of the old world leaf warblers (Phylloscopidae)

Sperm morphological traits are highly variable among species and are commonly thought to evolve by post‐copulatory sexual selection. However, little is known about the evolutionary dynamics of sperm morphology, and whether rates of evolutionary change are variable over time and among taxonomic groups. Here, we examine sperm morphology from 21 species of Old World leaf warblers (Phylloscopidae), a group of generally dull, sexually monochromatic birds, which are known to have high levels of extra‐pair paternity. We found that sperm length differs markedly across species, spanning about 40% of the range observed across a larger selection of passerine birds. Furthermore, we found strong support for an ‘early‐burst’ model of trait evolution, implying that the majority of divergence in sperm length has occurred early in the evolutionary history of this clade with subsequent evolutionary stasis. This large early divergence matches the early divergence reported in ecological traits (i.e. body size and feeding behaviour). Our findings demonstrate that rates of evolution in sperm morphology can change over time in passerine taxa, and that evolutionary stasis in sperm traits can occur even in species exhibiting characteristics consistent with moderate‐to‐high levels of sperm competition. It remains a major challenge to identify the selection mechanisms and possible constraints responsible for these variable rates of sperm evolution.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

Robust Demographic Inference from Genomic and SNP Data

We introduce a flexible and robust simulation-based framework to infer demographic parameters from the site frequency spectrum (SFS) computed on large genomic datasets. We show that our composite-likelihood approach allows one to study evolutionary models of arbitrary complexity, which cannot be tackled by other current likelihood-based methods. For simple scenarios, our approach compares favorably in terms of accuracy and speed with , the current reference in the field, while showing better convergence properties for complex models. We first apply our methodology to non-coding genomic SNP data from four human populations. To infer their demographic history, we compare neutral evolutionary models of increasing complexity, including unsampled populations. We further show the versatility of our framework by extending it to the inference of demographic parameters from SNP chips with known ascertainment, such as that recently released by Affymetrix to study human origins. Whereas previous ways of handling ascertained SNPs were either restricted to a single population or only allowed the inference of divergence time between a pair of populations, our framework can correctly infer parameters of more complex models including the divergence of several populations, bottlenecks and migration. We apply this approach to the reconstruction of African demography using two distinct ascertained human SNP panels studied under two evolutionary models. The two SNP panels lead to globally very similar estimates and confidence intervals, and suggest an ancient divergence (>110 Ky) between Yoruba and San populations. Our methodology appears well suited to the study of complex scenarios from large genomic data sets.     

芦鹀种群中的文化、遗传和形态进化的同时比较

我们比较了芦 (Emberizaschoeniclus)两个亚种组 ,即北部的薄喙亚种组和南部的厚喙亚种组的 10个种群中的文化、遗传和形态变异。使用了四个不同的变异标记物 ,其中两个用来测量文化分化 ,一个用来测量遗传分化 ,即微卫星等位基因的频次 ,一个用来测量种群的形态分化 ,即喙的高度。将遗传分化作为进化时间的尺度 ,我们计算了亚种组间和组内分化指标与所估计的进化率之间的相关性 ,发现只有文化定量指标和遗传分化与种群的形态分化相关 ,而两个文化分化指标之间没有关系 ,文化分化与遗传分化之间也没有关系。使用文化 -定量分化指标 ,发现亚种组间的文化进化率高于亚种内的文化进化率 ,提示鸣唱在防止杂交方面只有微弱的、也许是次要的作用。鸣唱定量特征的变异与微卫星频次相同 ,实际上在自然界中更可能是遗传决定的 ,这可以解释由于分析两个文化变异指标所得出的结果的不一致性。鸣唱的声学特性可能由于栖息地的差异或形态上的限制而发生了演变 ,而文化传播单位 (Meme)的特性可能由于学习鸣唱和文化传播而受到了影响     

The origin and evolution of porcine reproductive and respiratory syndrome viruses

Porcine reproductive and respiratory syndrome viruses (PRRSV) are divided into North American and European types, which show about 40% difference in their amino acid sequences. The divergence time of these two types has been estimated to be about 1980 from epidemiological data. This suggested that PRRSV have evolved at a higher evolutionary rate (order of 10(-2)/site/year) compared with other RNA viruses of 10(-3) to 10(-5)/site/year. Here, to test the evolutionary history of PRRSV speculated by the epidemiological background, we estimated the divergence time and evolutionary rate of PRRSV with molecular evolutionary analysis. Estimated divergence time (1972-1988) corresponded well to that estimated by the epidemiological data, and the evolutionary rate (4.71-9.8) x 10(-2) of PRRSV was indeed the highest among RNA viruses so far reported. Furthermore, we inferred important sites for the adaptation in order to examine how PRRSV have adapted to swine since they emerged. The adaptive sites were located not only in the epitopes related to immunity but also in the transmembrane regions including a signal peptide. In particular, the adaptive sites in the transmembrane regions were considered to affect compatibility to the host cell membrane. We conclude that PRRSV were transmitted from another host species to swine in about 1980 and have adapted to swine by altering the transmembrane regions.     

Calibrating the molecular clock beyond cytochrome b: assessing the evolutionary rate of COI in birds

Estimating the age of species or their component lineages based on sequence data is crucial for many studies in avian evolutionary biology. Although calibrations of the molecular clock in birds have been performed almost exclusively using cytochrome b (cyt b), they are commonly extrapolated to other mitochondrial genes. The existence of a large, standardized cytochrome c oxidase subunit I (COI) library generated as a result of the DNA barcoding initiative provides the opportunity to obtain a calibration for this mitochondrial gene in birds. In this study we compare the evolutionary rate of COI relative to cyt b across ten different avian orders. We obtained divergence estimates for both genes from nearly 300 phylogenetically independent pairs of species through the analysis of almost 5000 public sequences. For each pair of species we calculated the difference in divergence between COI and cyt b. Our results indicate that COI evolves on average 14% slower than cyt b, but also reveal considerable variation both among and within avian orders, precluding the use of this value as a standard adjustment for the COI molecular clock for birds. Our findings suggest that this variation is partially explained by a clear negative relationship between the difference in divergence in these genes and the age of species. Distances for cyt b are higher than those for COI for closely related species, but the values become similar as the divergence between the species increases. This appears to be the result of a stronger pattern of negative time‐dependency in the rate of cyt b than in that of COI, a difference that could be related to lower functional constraints on a small number of sites in cyt b that allow it to initially accumulate mutations more rapidly than COI.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Intragenomic and interspecific 5S rDNA sequence variation in five Asian pines

Patterns of intragenomic and interspecific variation of 5S rDNA in Pinus (Pinaceae) were studied by cloning and sequencing multiple 5S rDNA repeats from individual trees. Five pines, from both subgenera, Pinus and Strobus, were selected. The 5S rDNA repeat in pines has a conserved 120-base pair (bp) transcribed region and an intergenic spacer region of variable length (382-608 bp). The evolutionary rate in the spacer region is three- to sevenfold higher than in the genic region. We found substantial sequence divergence between the two subgenera. Intragenomic sequence heterogeneity was high for all species, and more than 86% of the clones within each individual were unique. The 5S gene tree revealed that different 5S repeats within individuals are polyphyletic, indicating that their ancestral divergence preceded the speciation events. The degrees of interspecific and intragenomic divergence among diploxylon pines are similar. The observed sequence patterns suggest that concerted evolution has been acting after the diversification of the two subgenera but very weak after the speciation of the four diploxylon pines. Sequence patterns in P. densata are consistent with hybrid origin. It had higher intragenomic diversity and maintained polymorphic copies of the parental types in addition to new and recombinant types unique to the hybrid.     

Evolutionary Mechanisms Affecting the Multivariate Divergence in Some Myotis Species (Chiroptera, Vespertilionidae)

Phenotypic evolutionary rates were measured for 27 craniometric characters in 12 extant OTUs from the bat genus Myotis (Chiroptera, Vespertilionidae). Squared Mahalanobis distance was used as a multivariate measure for amount of divergence, and squared Mahalanobis distance weighted by time was used as a measure for the rate of divergence. Estimates for the rates of divergence were found to be consistent with random walk hypothesis. Thus, the divergence in Myotis could be guided by random drift and mutations. The high dispersion in rate estimates suggests also a possible input of randomly fluctuating selection. The highest rates were recorded for divergence between M. myotis–M. blythii species group and the other OTUs. Rates of divergence between the subspecies of M. blythii occur to be lesser than rates of divergence between the earlier diverged species, their divergence could probably be slowed down by stabilizing selection. Size-adjusted data appeared to be lesser then the initial data, and it can be concluded that both size and shape were involved in divergence of Myotis species. The skeletal characters in bats are known to be extremely conservative during long-term evolution, however, the possibility for random walk at short time interval implies that bat evolution is constrained rather ecologically and biomechanically than genetically or developmentally.     

Evolution of placentally expressed cathepsins

Species and strain variants of a family of placentally expressed cathepsins (PECs) were cloned and sequenced in order to identify evolutionary conserved structural characteristics of this large family of cysteine proteases. Cathepsins M, P, Q, and R, are conserved in mice and rats but homologs of these genes are not found in human or rabbit placenta, showing that this family of proteases are probably restricted to rodents. Species-specific gene duplications have given rise to variants of cathepsin M in mice, and cathepsin Q in rats. Although the PECs have diverged at a greater rate than the other lysosomal cathepsins, residues around the specificity sub-sites of the individual enzymes are conserved. Strain-specific polymorphisms show that the evolutionary rate of divergence of cathepsins M and 3, the most recently duplicated pair of mouse genes, is even higher than the other PECs. In human placenta, critical functions of the PECs are probably performed by broader specificity proteases such as cathepsins B and L.     

The complex evolutionary dynamics of ancient and recent polyploidy in Leucaena (Leguminosae; Mimosoideae)

? Premise of the study: The evolutionary history of Leucaena has been impacted by polyploidy, hybridization, and divergent allopatric species diversification, suggesting that this is an ideal group to investigate the evolutionary tempo of polyploidy and the complexities of reticulation and divergence in plant diversification. ? Methods: Parsimony- and ML-based phylogenetic approaches were applied to 105 accessions sequenced for six sequence characterized amplified region-based nuclear encoded loci, nrDNA ITS, and four cpDNA regions. Hypotheses for the origin of tetraploid species were inferred using results derived from a novel species tree and established gene tree methods and from data on genome sizes and geographic distributions. ? Results: The combination of comprehensively sampled multilocus DNA sequence data sets and a novel methodology provide strong resolution and support for the origins of all five tetraploid species. A minimum of four allopolyploidization events are required to explain the origins of these species. The origin(s) of one tetraploid pair (L. involucrata/L. pallida) can be equally explained by two unique allopolyploidizations or a single event followed by divergent speciation. ? Conclusions: Alongside other recent findings, a comprehensive picture of the complex evolutionary dynamics of polyploidy in Leucaena is emerging that includes paleotetraploidization, diploidization of the last common ancestor to Leucaena, allopatric divergence among diploids, and recent allopolyploid origins for tetraploid species likely associated with human translocation of seed. These results provide insights into the role of divergence and reticulation in a well-characterized angiosperm lineage and into traits of diploid parents and derived tetraploids (particularly self-compatibility and year-round flowering) favoring the formation and establishment of novel tetraploids combinations.     

Understanding the functional difference between growth arrest-specific protein 6 and protein S: an evolutionary approach

Although protein S (PROS1) and growth arrest-specific protein 6 (GAS6) proteins are homologous with a high degree of structural similarity, they are functionally different. The objectives of this study were to identify the evolutionary origins from which these functional differences arose. Bioinformatics methods were used to estimate the evolutionary divergence time and to detect the amino acid residues under functional divergence between GAS6 and PROS1. The properties of these residues were analysed in the light of their three-dimensional structures, such as their stability effects, the identification of electrostatic patches and the identification potential protein–protein interaction. The divergence between GAS6 and PROS1 probably occurred during the whole-genome duplications in vertebrates. A total of 78 amino acid sites were identified to be under functional divergence. One of these sites, Asn463, is involved in N-glycosylation in GAS6, but is mutated in PROS1, preventing this post-translational modification. Sites experiencing functional divergence tend to express a greater diversity of stabilizing/destabilizing effects than sites that do not experience such functional divergence. Three electrostatic patches in the LG1/LG2 domains were found to differ between GAS6 and PROS1. Finally, a surface responsible for protein–protein interactions was identified. These results may help researchers to analyse disease-causing mutations in the light of evolutionary and structural constraints, and link genetic pathology to clinical phenotypes.     

Molecular distance and divergence time in carnivores and primates

Numerous studies have used indices of genetic distance between species to reconstruct evolutionary relationships and to estimate divergence time. However, the empirical relationship between molecular-based indices of genetic divergence and divergence time based on the fossil record is poorly known. To date, the results of empirical studies conflict and are difficult to compare because they differ widely in their choice of taxa, genetic techniques, or methods for calibrating rates of molecular evolution. We use a single methodology to analyze the relationship of molecular distance and divergence time in 86 taxa (72 carnivores and 14 primates). These taxa have divergence times of 0.01-55 Myr and provide a graded series of phylogenetic divergences such that the shape of the curve relating genetic distance and divergence time is often well defined. The techniques used to obtain genetic distance estimates include one- and two-dimensional protein electrophoresis, DNA hybridization, and microcomplement fixation. Our results suggest that estimates of molecular distance and divergence time are highly correlated. However, rates of molecular evolution are not constant; rather, in general they decline with increasing divergence time in a linear fashion. The rate of decline may differ according to technique and taxa. Moreover, in some cases the variability in evolutionary rates changes with increasing divergence time such that the accuracy of nodes in a phylogenetic tree varies predictably with time.     

The origin and evolution of Ebola and Marburg viruses

Molecular evolutionary analyses for Ebola and Marburg viruses were conducted with the aim of elucidating evolutionary features of these viruses. In particular, the rate of nonsynonymous substitutions for the glycoprotein gene of Ebola virus was estimated to be, on the average, 3.6 x 10(-5) per site per year. Marburg virus was also suggested to be evolving at a similar rate. Those rates were a hundred times slower than those of retroviruses and human influenza A virus, but were of the same order of magnitude as that of the hepatitis B virus. When these rates were applied to the degree of sequence divergence, the divergence time between Ebola and Marburg viruses was estimated to be more than several thousand years ago. Moreover, most of the nucleotide substitutions were transitions and synonymous for Marburg virus. This suggests that purifying selection has operated on Marburg virus during evolution.      

Hydraulic architecture and the evolution of shoot allometry in contrasting climates

We used pairs of congeneric shrub species from contrasting habitats to test for repeated evolutionary divergence in leaf-stem allometry and shoot hydraulic architecture in response to water availability. Allometric relationships and mean ratios between leaf size (individual and total area and mass per shoot) and stem cross-sectional area were compared between habitats using six species pairs representing three genera (Arctostaphylos, Baccharis, Ceanothus). We measured correlations among evolutionary changes in allometric, morphological, and physiological traits using phylogenetic independent contrasts. Allometric analysis revealed habitat differences: slopes were homogeneous among species (=1.46), but the more mesic-adapted species generally supported more leaf area at a common stem cross-sectional area. Reducing bivariate allometry to a ratio obscured this pattern because ratios varied with stem size, which was unrelated to habitat. Mean individual leaf size also was not correlated with either water availability or leaf-stem allometry. Stem hydraulic conductivity was generally lower in the xeric-adapted species of each pair, and its evolution mirrored changes in shoot allometry. This study provides evidence for repeated evolutionary divergence in shoot allometry and hydraulic architecture associated with water availability and demonstrates the importance of shoot allometry to water relations, independent of leaf size.