Indel Evolution Research Articles

FLOWERING LOCUS T indel variants confer vernalization-independent and photoperiod-insensitive flowering of yellow lupin (Lupinus luteus L.).

Ongoing climate change has considerably reduced the seasonal window for crop vernalization, concurrently expanding cultivation area into northern latitudes with long-day photoperiod. To address these changes, cool season legume breeders need to understand molecular control of vernalization and photoperiod. A key floral transition gene integrating signals from these pathways is the Flowering locus T (FT). Here, a recently domesticated grain legume, yellow lupin (Lupinus luteus L.), was explored for potential involvement of FT homologues in abolition of vernalization and photoperiod requirements. Two FTa (LlutFTa1a and LlutFTa1b) and FTc (LlutFTc1 and LlutFTc2) homologues were identified and sequenced for two contrasting parents of a reference recombinant inbred line (RIL) population, an early-flowering cultivar Wodjil and a late-flowering wild-type P28213. Large deletions were detected in the 5' promoter regions of three FT homologues. Quantitative trait loci were identified for flowering time and vernalization response in the RIL population and in a diverse panel of wild and domesticated accessions. A 2227bp deletion found in the LlutFTc1 promoter was linked with early phenology and vernalization independence, whereas LlutFTa1a and LlutFTc2 indels with photoperiod responsiveness. Comparative mapping highlighted convergence of FTc1 indel evolution in two Old World lupin species, addressing both artificial selection during domestication and natural adaptation to short season environmental conditions. We concluded that rapid flowering in yellow lupin is associated with the de-repression of the LlutFTc1 homologue from the juvenile phase, putatively due to the elimination of all binding sites in the promoter region for the AGAMOUS-like 15 transcription factor.

Using the Evolutionary History of Proteins to Engineer Insertion-Deletion Mutants from Robust, Ancestral Templates Using Graphical Representation of Ancestral Sequence Predictions (GRASP).

Analyzing the natural evolution of proteins by ancestral sequence reconstruction (ASR) can provide valuable information about the changes in sequence and structure that drive the development of novel protein functions. However, ASR has also been used as a protein engineering tool, as it often generates thermostable proteins which can serve as robust and evolvable templates for enzyme engineering. Importantly, ASR has the potential to provide an insight into the history of insertions and deletions that have occurred in the evolution of a protein family. Indels are strongly associated with functional change during enzyme evolution and represent a largely unexplored source of genetic diversity for designing proteins with novel or improved properties. Current ASR methods differ in the way they handle indels; inclusion or exclusion of indels is often managed subjectively, based on assumptions the user makes about the likelihood of each recombination event, yet most currently available ASR tools provide limited, if any, opportunities for evaluating indel placement in a reconstructed sequence. Graphical Representation of Ancestral Sequence Predictions (GRASP) is an ASR tool that maps indel evolution throughout a reconstruction and enables the evaluation of indel variants. This chapter provides a general protocol for performing a reconstruction using GRASP and using the results to create indel variants. The method addresses protein template selection, sequence curation, alignment refinement, tree building, ancestor reconstruction, evaluation of indel variants and approaches to library development.

ProPIP: a tool for progressive multiple sequence alignment with Poisson Indel Process

BackgroundCurrent alignment tools typically lack an explicit model of indel evolution, leading to artificially short inferred alignments (i.e., over-alignment) due to inconsistencies between the indel history and the phylogeny relating the input sequences.ResultsWe present a new progressive multiple sequence alignment tool ProPIP. The process of insertions and deletions is described using an explicit evolutionary model—the Poisson Indel Process or PIP. The method is based on dynamic programming and is implemented in a frequentist framework. The source code can be compiled on Linux, macOS and Microsoft Windows platforms. The algorithm is implemented in C++ as standalone program. The source code is freely available on GitHub at https://github.com/acg-team/ProPIP and is distributed under the terms of the GNU GPL v3 license.ConclusionsThe use of an explicit indel evolution model allows to avoid over-alignment, to infer gaps in a phylogenetically consistent way and to make inferences about the rates of insertions and deletions. Instead of the arbitrary gap penalties, the parameters used by ProPIP are the insertion and deletion rates, which have biological interpretation and are contextualized in a probabilistic environment. As a result, indel rate settings may be optimised in order to infer phylogenetically meaningful gap patterns.

A Probabilistic Model for Indel Evolution: Differentiating Insertions from Deletions.

Insertions and deletions (indels) are common molecular evolutionary events. However, probabilistic models for indel evolution are under-developed due to their computational complexity. Here, we introduce several improvements to indel modeling: 1) While previous models for indel evolution assumed that the rates and length distributions of insertions and deletions are equal, here we propose a richer model that explicitly distinguishes between the two; 2) we introduce numerous summary statistics that allow approximate Bayesian computation-based parameter estimation; 3) we develop a method to correct for biases introduced by alignment programs, when inferring indel parameters from empirical data sets; and 4) using a model-selection scheme, we test whether the richer model better fits biological data compared with the simpler model. Our analyses suggest that both our inference scheme and the model-selection procedure achieve high accuracy on simulated data. We further demonstrate that our proposed richer model better fits a large number of empirical data sets and that, for the majority of these data sets, the deletion rate is higher than the insertion rate.

A Model of Indel Evolution by Finite-State, Continuous-Time Machines.

We introduce a systematic method of approximating finite-time transition probabilities for continuous-time insertion-deletion models on sequences. The method uses automata theory to describe the action of an infinitesimal evolutionary generator on a probability distribution over alignments, where both the generator and the alignment distribution can be represented by pair hidden Markov models (HMMs). In general, combining HMMs in this way induces a multiplication of their state spaces; to control this, we introduce a coarse-graining operation to keep the state space at a constant size. This leads naturally to ordinary differential equations for the evolution of the transition probabilities of the approximating pair HMM. The TKF91 model emerges as an exact solution to these equations for the special case of single-residue indels. For the more general case of multiple-residue indels, the equations can be solved by numerical integration. Using simulated data, we show that the resulting distribution over alignments, when compared to previous approximations, is a better fit over a broader range of parameters. We also propose a related approach to develop differential equations for sufficient statistics to estimate the underlying instantaneous indel rates by expectation maximization. Our code and data are available at https://github.com/ihh/trajectory-likelihood.

Accelerating phylogeny-aware alignment with indel evolution using short time Fourier transform.

Recently we presented a frequentist dynamic programming (DP) approach for multiple sequence alignment based on the explicit model of indel evolution Poisson Indel Process (PIP). This phylogeny-aware approach produces evolutionary meaningful gap patterns and is robust to the ‘over-alignment’ bias. Despite linear time complexity for the computation of marginal likelihoods, the overall method’s complexity is cubic in sequence length. Inspired by the popular aligner MAFFT, we propose a new technique to accelerate the evolutionary indel based alignment. Amino acid sequences are converted to sequences representing their physicochemical properties, and homologous blocks are identified by multi-scale short-time Fourier transform. Three three-dimensional DP matrices are then created under PIP, with homologous blocks defining sparse structures where most cells are excluded from the calculations. The homologous blocks are connected through intermediate ‘linking blocks’. The homologous and linking blocks are aligned under PIP as independent DP sub-matrices and their tracebacks merged to yield the final alignment. The new algorithm can largely profit from parallel computing, yielding a theoretical speed-up estimated to be proportional to the cubic power of the number of sub-blocks in the DP matrices. We compare the new method to the original PIP approach and demonstrate it on real data.

Mutational Landscape of Pediatric Acute Myeloid Leukemia: A Report of the AML-BFM Study Group with a Targeted NGS Approach in 525 Patients Integrating De Novo, Relapsed and Secondary AML

Background: Recent studies in adults revealed the complexity of the genetic profiles of AML blasts that are used for fine-tuning of the risk-adapted treatment regimens. It is known, that recurrent translocations detected at diagnosis are largely preserved at relapse in the AML blasts. However, the molecular landscape and clonal evolution of small indels and missense alterations in the blasts assessed by NGS has not been well analyzed yet for pediatric AML patients, neither for primary and relapsed children and adolescents nor for patients with secondary AML. Here, we present the implications of our amplicon-based targeted sequencing efforts in a large cohort of pediatric patients with de novo AML or secondary AML at diagnosis and relapse. Patients and Methods: In total, samples from 525 AML patients (age <18 years) registered by the AML-BFM study group in Germany from 11/2000 until 04/2019 were analyzed here (de novo AML, n=449; secondary AML, n=13; first relapse n=63). Paired diagnostic and relapse samples were available from 45 patients. Sequencing was performed using the TruSight Myeloid panel (Illumina) on a MiseqDX device detecting somatic alterations in 54 genes associated with myeloid malignancies. On average, we achieved a 500X coverage in ≥95% of the target regions. We excluded intronic, synonymous and variants with an allele frequency below 5% and a read depth below 50 reads. Five-year probabilities of overall survival (pOS) and event-free survival (pEFS) were determined by Kaplan-Meier analysis in SPSS (IBM SPSS Statistics version 25). Results: In this large population of 525 pediatric patients, 632 mutations were detected in 40 out of 54 genes (74%). In total, 448 (80%) mutations in 9 genes (ASXL1, CEBPA, FLT3, KIT, KRAS, NPM1, NRAS, PTPN11 and WT1) were present in ≥3% of de novo AML patients. In ≥3% of relapsed patients, 48 (78%) mutations have been identified in 10 genes (CEBPA, FLT3, KIT, NPM1, NRAS, PHF6, PTPN11, RUNX1, TP53 and WT1). Although patient numbers were limited for secondary AML, data indicate a largely overlapping profile of affected genes. The median number of mutations per patient was low (n=1) in all 3 subgroups. However, 20% of the patients with de novo AML and 12% of relapsed patients had ≥ 3 mutations. Patients with more than 4 mutations have been only identified in de novo AML. The number of mutations per patient did not have a prognostic impact, analyzed in an uniformly treated subgroup of patients (pEFS: p=0.358; pOS: p=0.261). The comparison of paired samples (diagnosis/relapse) revealed the presence of three categories: Mutations only present at diagnosis for de novo AML samples, mutations that persisted through relapse and mutations that were newly gained at relapse. 39% of initially detected mutations (28 out of 72) persisted in the respective relapse samples with the highest stability in CEBPA, FLT3, KIT and WT1. In contrast, 61% (44 out of 72 initial mutations) were not detected and 45% (25 out of 54 relapse mutations) have been gained within the disease progress. The genetic profile changed from de novo AML to relapse in most patients (n = 34; pOS= 41 ± 15%). However, it has been found to be stable in 24% (n = 11; pOS 21±17%; p= 0.49) of the patients. Conclusions: This analysis elucidates the mutational landscape of both pediatric de novo and secondary AML as well as the genetic evolution of pediatric relapsed AML. Very high numbers of mutations per patient only occurred in de novo AML patients. In contrast to the assumed stability of recurrent translocations, the mutational landscape changed during disease progression in most of the patients. Whether these altered clones detectable at relapse were already present at the initial diagnosis, albeit at low frequency, or occurred de novo during treatment is very challenging to answer in pediatric AML due to the often very limited sample sizes. Either new single-cell sequencing approaches and/or analysis after expanding the cells in xenotransplantation models will provide more insights. Of note, modification of the sequencing panel would be beneficial for future studies due to the fact, that 26% of genes included were not affected by mutations. Disclosures Reinhardt: Jazz: Other: Participation in Advisory Boards, Research Funding; Novartis: Other: Participation in Advisory Boards; Roche: Research Funding; CSL Behring: Research Funding.

Progressive multiple sequence alignment with indel evolution

BackgroundSequence alignment is crucial in genomics studies. However, optimal multiple sequence alignment (MSA) is NP-hard. Thus, modern MSA methods employ progressive heuristics, breaking the problem into a series of pairwise alignments guided by a phylogeny. Changes between homologous characters are typically modelled by a Markov substitution model. In contrast, the dynamics of indels are not modelled explicitly, because the computation of the marginal likelihood under such models has exponential time complexity in the number of taxa. But the failure to model indel evolution may lead to artificially short alignments due to biased indel placement, inconsistent with phylogenetic relationship.ResultsRecently, the classical indel model TKF91 was modified to describe indel evolution on a phylogeny via a Poisson process, termed PIP. PIP allows to compute the joint marginal probability of an MSA and a tree in linear time. We present a new dynamic programming algorithm to align two MSAs –represented by the underlying homology paths– by full maximum likelihood under PIP in polynomial time, and apply it progressively along a guide tree. We have corroborated the correctness of our method by simulation, and compared it with competitive methods on an illustrative real dataset.ConclusionsOur MSA method is the first polynomial time progressive aligner with a rigorous mathematical formulation of indel evolution. The new method infers phylogenetically meaningful gap patterns alternative to the popular PRANK, while producing alignments of similar length. Moreover, the inferred gap patterns agree with what was predicted qualitatively by previous studies. The algorithm is implemented in a standalone C++ program: https://github.com/acg-team/ProPIP. Supplementary data are available at BMC Bioinformatics online.

Using the Neandertal genome to study the evolution of small insertions and deletions in modern humans

BackgroundSmall insertions and deletions occur in humans at a lower rate compared to nucleotide changes, but evolve under more constraint than nucleotide changes. While the evolution of insertions and deletions have been investigated using ape outgroups, the now available genome of a Neandertal can shed light on the evolution of indels in more recent times.ResultsWe used the Neandertal genome together with several primate outgroup genomes to differentiate between human insertion/deletion changes that likely occurred before the split from Neandertals and those that likely arose later. Changes that pre-date the split from Neandertals show a smaller proportion of deletions than those that occurred later. The presence of a Neandertal-shared allele in Europeans or Asians but the absence in Africans was used to detect putatively introgressed indels in Europeans and Asians. A larger proportion of these variants reside in intergenic regions compared to other modern human variants, and some variants are linked to SNPs that have been associated with traits in modern humans.ConclusionsOur results are in agreement with earlier results that suggested that deletions evolve under more constraint than insertions. When considering Neandertal introgressed variants, we find some evidence that negative selection affected these variants more than other variants segregating in modern humans. Among introgressed variants we also identify indels that may influence the phenotype of their carriers. In particular an introgressed deletion associated with a decrease in the time to menarche may constitute an example of a former Neandertal-specific trait contributing to modern human phenotypic diversity.

The Evolution of Small Insertions and Deletions in the Coding Genes of Drosophila melanogaster

Studies of protein evolution have focused on amino acid substitutions with much less systematic analysis on insertion and deletions (indels) in protein coding genes. We hence surveyed 7,500 genes between Drosophila melanogaster and D. simulans, using D. yakuba as an outgroup for this purpose. The evolutionary rate of coding indels is indeed low, at only 3% of that of nonsynonymous substitutions. As coding indels follow a geometric distribution in size and tend to fall in low-complexity regions of proteins, it is unclear whether selection or mutation underlies this low rate. To resolve the issue, we collected genomic sequences from an isogenic African line of D. melanogaster (ZS30) at a high coverage of 70× and analyzed indel polymorphism between ZS30 and the reference genome. In comparing polymorphism and divergence, we found that the divergence to polymorphism ratio (i.e., fixation index) for smaller indels (size ≤ 10 bp) is very similar to that for synonymous changes, suggesting that most of the within-species polymorphism and between-species divergence for indels are selectively neutral. Interestingly, deletions of larger sizes (size ≥ 11 bp and ≤ 30 bp) have a much higher fixation index than synonymous mutations and 44.4% of fixed middle-sized deletions are estimated to be adaptive. To our surprise, this pattern is not found for insertions. Protein indel evolution appear to be in a dynamic flux of neutrally driven expansion (insertions) together with adaptive-driven contraction (deletions), and these observations provide important insights for understanding the fitness of new mutations as well as the evolutionary driving forces for genomic evolution in Drosophila species.

Protein Insertions and Deletions Enabled by Neutral Roaming in Sequence Space

Backbone modifications via insertions and deletions (InDels) may exert dramatic effects, for better (mediating new functions) and for worse (causing loss of structure and/or function). However, contrary to point mutations (substitutions), our knowledge of the evolution and structural-functional effects of InDels is limited and so is our capability to engineer them. We sought to assess how deleterious InDels are relative to point mutations and understand the mechanisms that mediate their acceptance. Analysis of the evolution of InDels in orthologous protein phylogenies indicated that their rate of purging is 9- to 100-fold higher than for point mutations. In yeast, for example, the substitutions-to-InDels ratio is approximately 14-fold higher in protein coding than in noncoding regions. The incorporation of InDels relative to substitutions is not only slow but also nonlinear. On average, ≥50 substitutions accumulate before the appearance of the first InDel. We also found enriched substitutions in sequential and spatial proximity to InDels, suggesting that certain substitutions are correlated with InDels. As indicated by the lag in InDels accumulation, some of these correlated substitutions may have occurred first, as apparently neutral mutations, and later enabled the accumulation of InDels that would be otherwise purged. Thus, compensatory substitutions may follow InDels in an "adaptive walk" as traditionally assumed, but might also accumulate first, by "neutral roaming." The dynamics of InDels accumulation also depends on their genomic frequencies-InDels in flies are 4-fold more frequent than in yeast and tend to be compensated rather than enabled.

Evolution of protein indels in plants, animals and fungi

BackgroundInsertions/deletions (indels) in protein sequences are useful as drug targets, protein structure predictors, species diagnostics and evolutionary markers. However there is limited understanding of indel evolutionary patterns. We sought to characterize indel patterns focusing first on the major groups of multicellular eukaryotes.ResultsComparisons of complete proteomes from a taxonically broad set of primarily Metazoa, Fungi and Viridiplantae yielded 299 substantial (>250aa) universal, single-copy (in-paralog only) proteins, from which 901 simple (present/absent) and 3,806 complex (multistate) indels were extracted. Simple indels are mostly small (1-7aa) with a most frequent size class of 1aa. However, even these simple looking indels show a surprisingly high level of hidden homoplasy (multiple independent origins). Among the apparently homoplasy-free simple indels, we identify 69 potential clade-defining indels (CDIs) that may warrant closer examination. CDIs show a very uneven taxonomic distribution among Viridiplante (13 CDIs), Fungi (40 CDIs), and Metazoa (0 CDIs). An examination of singleton indels shows an excess of insertions over deletions in nearly all examined taxa. This excess averages 2.31 overall, with a maximum observed value of 7.5 fold.ConclusionsWe find considerable potential for identifying taxon-marker indels using an automated pipeline. However, it appears that simple indels in universal proteins are too rare and homoplasy-rich to be used for pure indel-based phylogeny. The excess of insertions over deletions seen in nearly every genome and major group examined maybe useful in defining more realistic gap penalties for sequence alignment. This bias also suggests that insertions in highly conserved proteins experience less purifying selection than do deletions.

Phylogenetic utility and evolution of indels: A study in neognathous birds

Indels are increasingly used in phylogenetics and play a major role in genome size evolution, and yet both the phylogenetic information content of indels and their evolutionary significance remain to be better assessed. Using three presumably independently evolving nuclear gene fragments (28S rDNA, β-fibrinogen, ornithine decarboxylase) from 29 families of neognathous birds, we have obtained a topology that is in general agreement with the current molecular consensus tree, supports the monophyly of Metaves, and provides evidence for the unresolved relationships within the Charadriiformes. Based on the retrieved topology, we assess the relative impact of indels and nucleotide substitutions and demonstrate that the superposition of the two kinds of data yields a topology that could not be obtained from either data set alone. Although only two out of three gene fragments reveal the deletion bias, the combined nucleotide insertion-to-deletion ratio is 0.22, indicating a rapid decrease of intron length. The average indel fixation rate in the neognaths is 2.5 times faster than that in therian (placental) mammals of similar geologic age. As in mammals, there is a considerable variation of indel fixation rate that is 1.5 times higher in Galloanseres compared to Neoaves, and 2.4 times higher in the Rallidae compared to the average for Neoaves (8.2 times higher compared to the related Gruidae). Our results add to the evidence that indel fixation rates correlate with lineage-specific evolutionary rates.

Comprehensive Phylogenetic Reconstruction of Amoebozoa Based on Concatenated Analyses of SSU-rDNA and Actin Genes

Evolutionary relationships within Amoebozoa have been the subject of controversy for two reasons: 1) paucity of morphological characters in traditional surveys and 2) haphazard taxonomic sampling in modern molecular reconstructions. These along with other factors have prevented the erection of a definitive system that resolves confidently both higher and lower-level relationships. Additionally, the recent recognition that many protosteloid amoebae are in fact scattered throughout the Amoebozoa suggests that phylogenetic reconstructions have been excluding an extensive and integral group of organisms. Here we provide a comprehensive phylogenetic reconstruction based on 139 taxa using molecular information from both SSU-rDNA and actin genes. We provide molecular data for 13 of those taxa, 12 of which had not been previously characterized. We explored the dataset extensively by generating 18 alternative reconstructions that assess the effect of missing data, long-branched taxa, unstable taxa, fast evolving sites and inclusion of environmental sequences. We compared reconstructions with each other as well as against previously published phylogenies. Our analyses show that many of the morphologically established lower-level relationships (defined here as relationships roughly equivalent to Order level or below) are congruent with molecular data. However, the data are insufficient to corroborate or reject the large majority of proposed higher-level relationships (above the Order-level), with the exception of Tubulinea, Archamoebae and Myxogastrea, which are consistently recovered. Moreover, contrary to previous expectations, the inclusion of available environmental sequences does not significantly improve the Amoebozoa reconstruction. This is probably because key amoebozoan taxa are not easily amplified by environmental sequencing methodology due to high rates of molecular evolution and regular occurrence of large indels and introns. Finally, in an effort to facilitate future sampling of key amoebozoan taxa, we provide a novel methodology for genome amplification and cDNA extraction from single or a few cells, a method that is culture-independent and allows both photodocumentation and extraction of multiple genes from natural samples.

Tracing the Most Parsimonious Indel History

Sequence alignment (the grouping of homologous bases into one column) is fundamental to almost any task in comparative genomics. This translates to positing gaps in the genomic sequences to account for events of insertions and deletions (indels). The interrelationship between sequence alignment and phylogenetic reconstruction has drawn substantial attention recently with works showing the significance of differences in alignments. One of the plausible approaches in this direction is to grade the suitability of a tree to an associated alignment and vice verse. We here present a combinatorial (as opposed to statistical) approach based on the indel history. We show--both by simulations and by using real biological data from the Encyclopedia of DNA Elements (ENCODE)--that this criterion is sound. The novelty of our approach is the distinguishing between insertions and deletions, and augmenting the analysis with a dimension of "depth," extending it from the sequence space to the phylogenetic space. Using this approach, we perform a comprehensive study of indel characteristic behavior among mammals in both coding and non-coding regions. Our results show significant differences in indel patterns between coding and non-coding regions. We also show other characteristic patterns of indel evolution in the depth of the underlying phylogeny.

Parameter Estimation in Multiple-Hidden I.I.D. Models from Biological Multiple Alignment

In this work we deal with parameter estimation in a latent variable model, namely the multiple-hidden i.i.d. model, which is derived from multiple alignment algorithms. We first provide a rigorous formalism for the homology structure of k sequences related by a star-shaped phylogenetic tree in the context of multiple alignment based on indel evolution models. We discuss possible definitions of likelihoods and compare them to the criterion used in multiple alignment algorithms. Existence of two different Information divergence rates is established and a divergence property is shown under additional assumptions. This would yield consistency for the parameter in parametrization schemes for which the divergence property holds. We finally extend the definition of the multiple-hidden i.i.d. model and the results obtained to the case in which the sequences are related by an arbitrary phylogenetic tree. Simulations illustrate different cases which are not covered by our results.

Phylogeny and diversification of diving beetles (Coleoptera: Dytiscidae).

Dytiscidae is the most diverse family of beetles in which both adults and larvae are aquatic, with examples of extreme morphological and ecological adaptations. Despite continuous attention from systematists and ecologists, existing phylogenetic hypotheses remain unsatisfactory because of limited taxon sampling or low node support. Here we provide a phylogenetic tree inferred from four gene fragments (cox1, rrnL, H3 and SSU, ≈ 4000 aligned base pairs), including 222 species in 116 of 174 known genera and 25 of 26 tribes. We aligned ribosomal genes prior to tree building with parsimony and Bayesian methods using three approaches: progressive pair-wise alignment with refinement, progressive alignment modeling the evolution of indels, and deletion of hypervariable sites. Results were generally congruent across alignment and tree inference methods. Basal relationships were not well defined, although we identified 28 well supported lineages corresponding to recognized tribes or groups of genera, among which the most prominent novel results were the polyphyly of Dytiscinae; the grouping of Pachydrini with Bidessini, Peschetius with Methlini and Coptotomus within Copelatinae; the monophyly of all Australian Hydroporini (Necterosoma group), and their relationship with the Graptodytes and Deronectes groups plus Hygrotini. We found support for a clade formed by Hydroporinae plus Laccophilini, and their sister relationship with Cybistrini and Copelatinae. The tree provided a framework for the analysis of species diversification in Dytiscidae. We found a positive correlation between the number of species in a lineage and the age of the crown group as estimated through a molecular clock approach, but the correlation with the stem age was non-significant. Imbalances between sister clades were significant for several nodes, but the residuals of the regression of species numbers with the crown age of the group identified only Bidessini and the Coptotomus + Agaporomorphus clade as lineages with, respectively, above and below expected levels of species diversity. © The Willi Hennig Society 2008.

Inference of higher-order conifer relationships from a multi-locus plastid data setThis paper is one of a selection of papers published in the Special Issue on Systematics Research.

We reconstructed the broad backbone of conifer phylogeny from a survey of 15–17 plastid loci and associated noncoding regions from exemplar conifer species. Parsimony and likelihood analyses recover the same higher-order relationships, and we find strong support for most of the deep splits in conifer phylogeny, including those within our two most heavily sampled families, Araucariaceae and Cupressaceae. Our findings are broadly congruent with other recent studies, and are inferred with comparable or improved bootstrap support. The deepest phylogenetic split in conifers is inferred to be between Pinaceae and all other conifers (Cupressophyta). Our current gene and taxon sampling does not support a relationship between Pinaceae and Gnetales, observed in some published studies. Within the Cupressophyta clade, we infer well-supported relationships among Cephalotaxaceae, Cupressaceae, Sciadopityaceae, and Taxaceae. Our data support recent moves to recognize Cephalotaxus under Taxaceae, and we find strong support for a sister-group relationship between the two predominantly southern hemisphere conifer families, Araucariaceae and Podocarpaceae. A local hotspot of indel evolution shared by the latter two conifer families is identified in the coding portion of one of the plastid ribosomal protein genes. The removal of the most rapidly evolving plastid characters, as defined using a likelihood-based classification of substitution rates for the taxa considered here, is shown to have little to no effect on our inferences of higher-order conifer relationships.

Fifty mosses on five trees: comparing phylogenetic information in three types of non-coding mitochondrial DNA and two chloroplast loci

Given the frequent genomic recombinations in plant mitochondrial DNA, intergenic regions of this organelle genome had so far not been considered as loci of potential phylogenetic information. Based on the recent evidence that an evolutionary ancient mitochondrial nad5-nad4 gene continuum is conserved in bryophytes we have compiled a dataset for a phylogenetically wide sampling of 50 mosses covering this intergenic region. The length of the intergenic region was generally in the range of 585 (Diphyscium) to 646 bp (Tomentypnum) with rare exceptions, for example all Polytrichales taxa showing a 200 bp deletion as an apparent synapomorphy of this order. Phylogenetic information in the novel marker sequence was compared with that of a nad5 gene region containing a group I intron and a nad2 gene region containing a group II intron as well as with two widely sampled chloroplast data sets, rbcL and rps4. Indel evolution in the three types of non-coding mitochondrial sequences is obviously more taxon-dependent than locus-dependent, indicating lineage-specific insertion/deletion rates. For example, larger sequence deletions are a general feature in Schistostega and Tetraphis. Although confidence for particular nodes in the phylogeny was found to vary among data sets, gene trees were essentially without conflict with respect to well-supported ones and add up in information towards a reasonably well-resolved moss phylogeny. However, while a consensus on the latter is clearly emerging, sufficient confidence is still lacking for the first dichotomies among the arthrodontous mosses leading into subclasses Bryidae, Dicranidae and Funariidae and the relative placement of nematodontous mosses (Polytrichales and Tetraphidales) on the backbone of early moss phylogeny.

Whole-Gene Positive Selection, Elevated Synonymous Substitution Rates, Duplication, and Indel Evolution of the Chloroplast clpP1 Gene

BackgroundSynonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.Methodology/Principle FindingsWe found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.Conclusions/SignificanceWe found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.