The origin and maintenance of nuclear endosperms: viewing development through a phylogenetic lens

The endosperm develops in fertilized ovules of angiosperms following fertilization of the central cell and nuclei in the female gametophyte. Endosperms differ in whether, and which, nuclear divisions are followed by cellular divisions; the variants are classified as cellular, nuclear or helobial. Functional correlates of this variation are little understood. Phylogenetic methods provide a powerful means of exploring taxonomic variation and phylogenetic patterns, to frame questions regarding biological processes. Data on endosperms across angiosperms were analysed in a phylogenetic context in order to determine homologies and detect biases in the direction of evolutionary transitions. Analyses confirm that neither all nuclear nor all helobial endosperms are homologous, raise the possibility that cellular development is a reversal in some derived angiosperms (e.g. asterids) and show that a statistically significant bias towards evolution of nuclear endosperms (and against reversals) prevails in angiosperms as a whole. This bias suggests strong selective advantages to having nuclear endosperm, developmental constraints to reversals or both. Homologies suggest that the microtubular cycle and cellularization pattern characteristic of reproductive cells across land plants may have been independently co-opted during multiple origins of nuclear endosperms, but information on cellular endosperms is essential to investigate further.     

Detection of dispersed short tandem repeats using reversible jump Markov chain Monte Carlo

Tandem repeats occur frequently in biological sequences. They are important for studying genome evolution and human disease. A number of methods have been designed to detect a single tandem repeat in a sliding window. In this article, we focus on the case that an unknown number of tandem repeat segments of the same pattern are dispersively distributed in a sequence. We construct a probabilistic generative model for the tandem repeats, where the sequence pattern is represented by a motif matrix. A Bayesian approach is adopted to compute this model. Markov chain Monte Carlo (MCMC) algorithms are used to explore the posterior distribution as an effort to infer both the motif matrix of tandem repeats and the location of repeat segments. Reversible jump Markov chain Monte Carlo (RJMCMC) algorithms are used to address the transdimensional model selection problem raised by the variable number of repeat segments. Experiments on both synthetic data and real data show that this new approach is powerful in detecting dispersed short tandem repeats. As far as we know, it is the first work to adopt RJMCMC algorithms in the detection of tandem repeats.     

Inferring phylogenetic networks by the maximum parsimony criterion: a case study

Horizontal gene transfer (HGT) may result in genes whose evolutionary histories disagree with each other, as well as with the species tree. In this case, reconciling the species and gene trees results in a network of relationships, known as the "phylogenetic network" of the set of species. A phylogenetic network that incorporates HGT consists of an underlying species tree that captures vertical inheritance and a set of edges which model the "horizontal" transfer of genetic material. In a series of papers, Nakhleh and colleagues have recently formulated a maximum parsimony (MP) criterion for phylogenetic networks, provided an array of computationally efficient algorithms and heuristics for computing it, and demonstrated its plausibility on simulated data. In this article, we study the performance and robustness of this criterion on biological data. Our findings indicate that MP is very promising when its application is extended to the domain of phylogenetic network reconstruction and HGT detection. In all cases we investigated, the MP criterion detected the correct number of HGT events required to map the evolutionary history of a gene data set onto the species phylogeny. Furthermore, our results indicate that the criterion is robust with respect to both incomplete taxon sampling and the use of different site substitution matrices. Finally, our results show that the MP criterion is very promising in detecting HGT in chimeric genes, whose evolutionary histories are a mix of vertical and horizontal evolution. Besides the performance analysis of MP, our findings offer new insights into the evolution of 4 biological data sets and new possible explanations of HGT scenarios in their evolutionary history.     

Sorting by restricted-length-weighted reversals

Classical sorting by reversals uses the unit-cost model, that is, each reversal consumes an equal cost. This model limits the biological meaning of sorting by reversal. Bender and his colleagues extended it by assigning a cost function f（1） = l^a for all a≥ 0, where l is the length of the reversed subsequence. In this paper, we extend their results by considering a model in which long reversals are prohibited. Using the same cost function above for permitted reversals, we present tight or nearly tight bounds for the worst-case cost of sorting by reversals. Then we develop algorithms to approximate the optimal cost to sort a given 0/1 sequence as well as a given permutation. Our proposed problems are more biologically meaningful and more algorithmically general and challenging than the problem considered by Bender et al. Furthermore, our bounds are tight and nearly tight, whereas our algorithms provide good approximation ratios compared to the optimal cost to sort 0/1 sequences or permutations by reversals.     

Molecular evolutionary processes and conflicting gene trees: The hominoid case

Molecular evolutionary processes modify DNA over time, creating both newly derived substitutions shared by related descendant lineages (phylogenetic signal) and “false” similarities which confound phylogenetic reconstruction (homoplasy). However, some types of DNA regions, for example those containing tandem duplicate repeats, are preferentially subject to homoplasy-inducing processes such as sporadically occurring concerted evolution and DNA insertion/deletion. This added level of homoplasic “noise” can make DNA regions with repeats less reliable in phylogenetic reconstruction than those without repeats. Most molecular datasets which distinguish among African hominoids support a human-chimpanzee clade; the most notable exception is from the involucrin gene. However, phylogenetic resolution supporting a chimpanzee-gorilla clade is based entirely on involucrin DNA repeat regions. This is problematic because (1) involucrin repeats are difficult to align, and published alignments are contradictory; (2) involucrin repeats are subject to DNA insertion/deletion; (3) gorillas are polymorphic in that some do not have repeats reported to be synapomorphies linking chimpanzees and gorillas. Gene tree/species tree conflicts can occur due to the sorting of ancestrally polymorphic alleles during speciation. Because hominoid females transfer between groups, mitochondrial and nuclear gene flow occur to the same extent, and the probability of conflict between mitochondrial and nuclear gene trees is theoretically low. When hominoid intraspecific mitochondrial variability is taken into account [based on cytochrome oxidase subunit II (COII) gene sequences], humans and chimpanzees are most closely related, showing the same relative degree of separation from gorillas as when single individuals representing species are analyzed. Conflicting molecular phylogenies can be explained in terms of molecular evolutionary processes and sorting of ancient polymorphisms. This perspective can enhance our understanding of hominoid molecular phylogenies. © 1994 Wiley-Liss, Inc.     

Dispersal between shallow and abyssal seas and evolutionary loss and regain of compound eyes in cylindroleberidid ostracods: conflicting conclusions from different comparative methods

Complex organs such as eyes are commonly lost during evolution, but the timescale on which lost phenotypes could be reactivated is a matter of long-standing debate, with important implications for the molecular mechanisms of trait loss. Two phylogenetic approaches have been used to test whether regain of traits has occurred. One way is by comparison of nested, continuous-time Markov models of trait evolution, approaches that we term tree-based tests. A second way to demonstrate statistical support for trait regain is through use of node-based tests that employ explicit estimation of ancestral node states. Here, we estimate new molecular and morphological phylogenies and use them to examine the possibility of eye regain and dispersal between abyssal and shallow seas during the history of cylindroleberidid ostracods, a family of about 200 species, comprising both eyeless and sighted species. First, we confirmed that eye presence/absence is correlated with habitat depth. Parameter estimates from a phylogenetic model indicate that speciation is more rapid in deep-sea eyeless clades compared with shallow-water sighted clades. In addition, we found that tree-based statistical tests usually indicated reversals, including both transitions from deep to shallow seas and regain of eyes. In contrast, node-based statistical tests usually failed to show significant support for reversals. These results also hold for simulated phylogenies, indicating that they are not unique to the current data set. We recommend that both tree-based and node-based tests should be examined before making conclusions about character reversal and that ideally, alternative character histories should be tested using additional data, besides just the phylogenetic distribution of presence/absence of the characters.     

Solving the preserving reversal median problem

Genomic rearrangement operations can be very useful to infer the phylogenetic relationship of gene orders representing species. We study the problem of finding potential ancestral gene orders for the gene orders of given taxa, such that the corresponding rearrangement scenario has a minimal number of reversals, and where each of the reversals has to preserve the common intervals of the given input gene orders. Common intervals identify sets of genes that occur consecutively in all input gene orders. The problem of finding such an ancestral gene order is called the preserving reversal median problem (pRMP). A tree-based data structure for the representation of the common intervals of all input gene orders is used in our exact algorithm TCIP for solving the pRMP. It is known that the minimum number of reversals to transform one gene order into another can be computed in polynomial time, whereas the corresponding problem with the restriction that common intervals should not be destroyed is already NP-hard. It is shown theoretically that TCIP can solve a large class of pRMP instances in polynomial time. Empirically we show the good performance of TCIP on biological and artificial data.     

Multiple ITS Haplotypes in the Genome of the Lichenized Basidiomycete Cora inversa (Hygrophoraceae): Fact or Artifact?

The internal transcribed spacer region (ITS) of the nuclear rDNA cistron represents the barcoding locus for Fungi. Intragenomic variation of this multicopy gene can interfere with accurate phylogenetic reconstruction of biological entities. We investigated the amount and nature of this variation for the lichenized fungus Cora inversa in the Hygrophoraceae (Basidiomycota: Agaricales), analyzing base call and length variation in ITS1 454 pyrosequencing data of three samples of the target mycobiont, for a total of 16,665 reads obtained from three separate repeats of the same samples under different conditions. Using multiple fixed alignment methods (PaPaRa) and maximum likelihood phylogenetic analysis (RAxML), we assessed phylogenetic relationships of the obtained reads, together with Sanger ITS sequences from the same samples. Phylogenetic analysis showed that all ITS1 reads belonged to a single species, C. inversa. Pyrosequencing data showed 266 insertion sites in addition to the 325 sites expected from Sanger sequences, for a total of 15,654 insertions (0.94 insertions per read). An additional 3,279 substitutions relative to the Sanger sequences were detected in the dataset, out of 5,461,125 bases to be called. Up to 99.3 % of the observed indels in the dataset could be interpreted as 454 pyrosequencing errors, approximately 65 % corresponding to incorrectly recovered homopolymer segments, and 35 % to carry-forward-incomplete-extension errors. Comparison of automated clustering and alignment-based phylogenetic analysis demonstrated that clustering of these reads produced a 35-fold overestimation of biological diversity in the dataset at the 95 % similarity threshold level, whereas phylogenetic analysis using a maximum likelihood approach accurately recovered a single biological entity. We conclude that variation detected in 454 pyrosequencing data must be interpreted with great care and that a combination of a sufficiently large number of reads per taxon, a set of Sanger references for the same taxon, and at least two runs under different emulsion PCR and sequencing conditions, are necessary to reliably separate biological variation from 454 sequencing errors. Our study shows that clustering methods are highly sensitive to artifactual sequence variation and inadequate to properly recover biological diversity in a dataset, if sequencing errors are substantial and not removed prior to clustering analysis.     

Evolution of Alu repeats. Imitation model

At present, nucleotide sequences of 100 different Alu repeats are known, i.e. 0.01% of the total number of Alu repeats in the genome. It is clear that one can not refer the evolutionary characteristics of Alu repeats obtained from the analysis of the available limited sample to all Alu repeats comprised in the genome, without additional statistical estimations. For supplementary investigation of such average evolutionary characteristics as the extent of intraspecific divergence, the rate of Alu repeats transposition (insertion, excision), we used the method of imitation simulation of the process of Alu repeats transposition in the genome. As a result of simulation, phylogenetic relations were obtained among all Alu repeats. It was shown that the evolutionary characteristics evaluated for different samples of repeats were similar. It was proved that the extent of divergence of Alu repeats in the model is twice as small as that evaluated, according to the real data (0.15, instead of 0.3). Possible reasons for such discrepancy have been discussed.     

THE LIMITS OF AMINO ACID SEQUENCE DATA IN ANGIOSPERM PHYLOGENETIC RECONSTRUCTION

Amino acid sequence data are available for ribulose biphosphate carboxylase, plastocyanin, cytochrome c, and ferredoxin for a number of angiosperm families. Cladistic analysis of the data, including evaluation of all equally or almost equally parsimonious cladograms, shows that much homoplasy (parallelisms and reversals) is present and that few or no well supported monophyletic groups of families can be demonstrated. In one analysis of nine angiosperm families and 40 variable amino acid positions from three proteins, the most parsimonious cladograms were 151 steps long and contained 63 parallelisms and reversals (consistency index = 0.583). In another analysis of six families and 53 variable amino acid positions from four proteins, the most parsimonious cladogram was 161 steps long and contained 50 parallelisms and reversals (consistency index = 0.689). Single changes in both data matrices could yield most parsimonious cladograms with quite different topologies and without common monophyletic groups. Presently, amino acid sequence data are not comprehensive enough for phylogenetic reconstruction among angiosperms. More informative positions are needed, either from sequencing longer parts of the proteins or from sequencing more proteins from the same taxa.     

Evolution and phylogenetic information content of mitochondrial genomic structural features illustrated with acrodont lizards

DNA sequences from 195 squamate reptiles indicate that mitochondrial gene order is the most reliable phylogenetic character establishing monophyly of acrodont lizards and of the snake families Boidae, Colubridae, and Viperidae. Gene order shows no evidence of evolutionary parallelisms or reversals in these taxa. Derived secondary structures of mitochondrial tRNAs also prove to be useful phylogenetic characters showing no reversals. Parallelisms for secondary structures of tRNAs are restricted to deep lineages that are separated by at least 200 million years of independent evolution. Presence of a stem-and-loop structure between the genes encoding tRNA(Asn) and tRNA(Cys), where the replication origin for light-strand synthesis is typically located in vertebrate mitochondrial genomes, is found to undergo at least three and possibly as many as seven evolutionary shifts, most likely parallel losses. This character is therefore a less desirable phylogenetic marker than the other structural changes examined. Sequencing regions that contain multiple genes, including tRNA genes, may be preferable to the common practice of obtaining single-gene fragments for phylogenetic inference because it permits observation of major structural changes in the mitochondrial genome. Such characters may occasionally provide phylogenetic information on relatively short internal branches for which base substitutional changes are expected to be relatively uninformative.     

The Tandem Repeats Enabling Reversible Switching between the Two Phases of β-Lactamase Substrate Spectrum

Expansion or shrinkage of existing tandem repeats (TRs) associated with various biological processes has been actively studied in both prokaryotic and eukaryotic genomes, while their origin and biological implications remain mostly unknown. Here we describe various duplications (de novo TRs) that occurred in the coding region of a β-lactamase gene, where a conserved structure called the omega loop is encoded. These duplications that occurred under selection using ceftazidime conferred substrate spectrum extension to include the antibiotic. Under selective pressure with one of the original substrates (amoxicillin), a high level of reversion occurred in the mutant β-lactamase genes completing a cycle back to the original substrate spectrum. The de novo TRs coupled with reversion makes a genetic toggling mechanism enabling reversible switching between the two phases of the substrate spectrum of β-lactamases. This toggle exemplifies the effective adaptation of de novo TRs for enhanced bacterial survival. We found pairs of direct repeats that mediated the DNA duplication (TR formation). In addition, we found different duos of sequences that mediated the DNA duplication. These novel elements—that we named SCSs (same-strand complementary sequences)—were also found associated with β-lactamase TR mutations from clinical isolates. Both direct repeats and SCSs had a high correlation with TRs in diverse bacterial genomes throughout the major phylogenetic lineages, suggesting that they comprise a fundamental mechanism shaping the bacterial evolution.     

Deep homology: A view from systematics

Over the past decade, it has been discovered that disparate aspects of morphology – often of distantly related groups of organisms – are regulated by the same genetic regulatory mechanisms. Those discoveries provide a new perspective on morphological evolutionary change. A conceptual framework for exploring these research findings is termed ‘deep homology’. A comparative framework for morphological relations of homology is provided that distinguishes analogy, homoplasy, plesiomorphy and synapomorphy. Four examples – three from plants and one from animals – demonstrate that homologous developmental mechanisms can regulate a range of morphological relations including analogy, homoplasy and examples of uncertain homology. Deep homology is part of a much wider range of phenomena in which biological (genes, regulatory mechanisms, morphological traits) and phylogenetic levels of homology can both be disassociated. Therefore, to understand homology, precise, comparative, independent statements of both biological and phylogenetic levels of homology are necessary.     

Potentials and limitations of histone repeat sequences for phylogenetic reconstruction of Sophophora.

Simplified DNA sequence acquisition has provided many new data sets that are useful for phylogenetic reconstruction, including single- and multiple-copy nuclear and organellar genes. Although transcribed regions receive much attention, nontranscribed regions have recently been added to the repertoire of sequences suitable for phylogenetic studies, especially for closely related taxa. We evaluated the efficacy of a small portion of the histone repeat for phylogenetic reconstruction among Drosophila species. Histone repeats in invertebrates offer distinct advantages similar to those of widely used ribosomal repeats. First, the units are tandemly repeated and undergo concerted evolution. Second, histone repeats include both highly conserved coding and variable intergenic regions. This composition facilitates application of "universal" primers spanning potentially informative sites. We examined a small region of the histone repeat, including the intergenic spacer segments of coding regions from the divergently transcribed H2A and H2B histone genes. The spacer (about 230 bp) exists as a mosaic with highly conserved functional motifs interspersed with rapidly diverging regions; the former aid in alignment of the spacer. There are no ambiguities in alignment of coding regions. Coding and noncoding regions were analyzed together and separately for phylogenetic information. Parsimony, distance, and maximum-likelihood methods successfully retrieve the corroborated phylogeny for the taxa examined. This study demonstrates the resolving power of a small histone region which may now be added to the growing collection of phylogenetically useful DNA sequences.     

Distribution of phylogenetic diversity under random extinction

Phylogenetic diversity is a measure for describing how much of an evolutionary tree is spanned by a subset of species. If one applies this to the unknown subset of current species that will still be present at some future time, then this ‘future phylogenetic diversity’ provides a measure of the impact of various extinction scenarios in biodiversity conservation. In this paper, we study the distribution of future phylogenetic diversity under a simple model of extinction (a generalized ‘field of bullets’ model). We show that the distribution of future phylogenetic diversity converges to a normal distribution as the number of species grows, under mild conditions, which are necessary. We also describe an algorithm to compute the distribution efficiently, provided the edge lengths are integral, and briefly outline the significance of our findings for biodiversity conservation.     

Mining functional microsatellites in legume unigenes

Highly polymorphic and transferable microsatellites (SSRs) are important for comparative genomics, genome analysis and phylogenetic studies. Development of novel species-specific microsatellite markers remains a costly and labor-intensive project. Therefore, interest has been shifted from genomic to genic markers owing to their high inter-species transferability as they are developed from conserved coding regions of the genome. This study concentrates on comparative analysis of genic microsatellites in nine important legume (Arachis hypogaea, Cajanus cajan, Cicer arietinum, Glycine max, Lotus japonicus, Medicago truncatula, Phaseolus vulgaris, Pisum sativum and Vigna unguiculata) and two model plant species (Oryza sativa and Arabidopsis thaliana). Screening of a total of 228090 putative unique sequences spanning 219610522 bp using a microsatellite search tool, MISA, identified 12.18% of the unigenes containing 36248 microsatellite motifs excluding mononucleotide repeats. Frequency of legume unigene-derived SSRs was one SSR in every 6.0 kb of analyzed sequences. The trinucleotide repeats were predominant in all the unigenes with the exception of C. cajan, which showed prevalence of dinucleotide repeats over trinucleotide repeats. Dinucleotide repeats along with trinucleotides counted for more than 90% of the total microsatellites. Among dinucleotide and trinucleotide repeats, AG and AAG motifs, respectively, were the most frequent. Microsatellite positive chickpea unigenes were assigned Gene Ontology (GO) terms to identify the possible role of unigenes in various molecular and biological functions. These unigene based microsatellite markers will prove valuable for recording allelic variance across germplasm collections, gene tagging and searching for putative candidate genes.     

A Phylogenetic Perspective on Sequence Evolution in Microsatellite Loci

We examined the evolution of the repeat regions of three noncoding microsatellite loci in 58 species of the Polistinae, a subfamily of wasps that diverged over 140 million years ago. A phylogenetic approach allows two new kinds of approaches to studying microsatellite evolution: character mapping and comparative analysis. The basic repeat structure of the loci was highly conserved, but was often punctuated with imperfections that appear to be phylogenetically informative. Repeat numbers evolved more rapidly than other changes in the repeat region. Changes in number of repeats among species seem consistent with the stepwise mutation model, which is based on slippage during replication as the main source of mutations. Changes in repeat numbers can occur even when there are very few tandem repeats but longer repeats, especially perfect repeats led to greater rates of evolutionary change. Species phylogenetically closer to the one from which we identified the loci had longer stretches of uninterrupted repeats and more different motifs, but not longer total repeat regions. The number of perfect repeats increased more often than it decreased. However, there was no evidence that some species have consistently greater numbers of repeats across loci than other species have, once ascertainment bias is eliminated. We also found no evidence for a population size effect posited by one form of the directionality hypothesis. Overall, phylogenetic variation in repeat regions can be explained by adding neutral evolution to what is already known about the mutation process. The life cycle of microsatellites appears to reflect a balance between growth by slippage and degradation by an essentially irreversible accumulation of imperfections. Received: 13 April 1999 / Accepted: 8 September 1999     

Complete chloroplast genome sequences of Dipterygium glaucum and Cleome chrysantha and other Cleomaceae Species,comparative analysis and phylogenetic relationships

This current study presents, for the first time, the complete chloroplast genome of two Cleomaceae species: Dipterygium glaucum and Cleome chrysantha in order to evaluate the evolutionary relationship. The cp genome is 158,576 bp in length with 35.74% GC content in D. glaucum and 158,111 bp with 35.96% GC in C. chrysantha. Inverted repeats IR 26,209 bp, 26,251 bp each, LSC of 87,738 bp, 87,184 bp and SSC of 18,420 bp, 18,425 bp respectively. There are 136 genes in the genome, which includes 80 protein coding genes, 31 tRNA genes and four rRNA genes were observed in both chloroplast genomes. 117 genes are unique while the remaining 19 genes are duplicated in IR regions. The analysis of repeats shows that the cp genome includes all types of repeats with more frequent occurrences of palindromic; Also, this analysis indicates that the total number of simple sequence repeats (SSR) were 323 in D. glaucum, and 313 in C. chrysantha, of which the majority of the SSRs in these plastid genomes were mononucleotide repeats A/T which are located in the intergenic spacer. Moreover, the comparative analysis of the four cp sequences revealed four hotspot genes (atpF, rpoC2, rps19, and ycf1), these variable regions could be used as molecular makers for the species authentication as well as resources for inferring phylogenetic relationships of the species. All the relationships in the phylogenetic tree are with high support, this indicate that the complete chloroplast genome is a useful data for inferring phylogenetic relationship within the Cleomaceae and other families. The simple sequence repeats identified will be useful for identification, genetic diversity, and other evolutionary studies of the species. This study reported the first cp genome of the genus Dipterygium and Cleome. The finding of this study will be beneficial for biological disciplines such as evolutionary and genetic diversity studies of the species within the core Cleomaceae.