Phylogenetic signal from rearrangements in 18 Anopheles species by joint scaffolding extant and ancestral genomes

Background

Genomes rearrangements carry valuable information for phylogenetic inference or the elucidation of molecular mechanisms of adaptation. However, the detection of genome rearrangements is often hampered by current deficiencies in data and methods: Genomes obtained from short sequence reads have generally very fragmented assemblies, and comparing multiple gene orders generally leads to computationally intractable algorithmic questions.

Results

We present a computational method, ADseq, which, by combining ancestral gene order reconstruction, comparative scaffolding and de novo scaffolding methods, overcomes these two caveats. ADseq provides simultaneously improved assemblies and ancestral genomes, with statistical supports on all local features. Compared to previous comparative methods, it runs in polynomial time, it samples solutions in a probabilistic space, and it can handle a significantly larger gene complement from the considered extant genomes, with complex histories including gene duplications and losses. We use ADseq to provide improved assemblies and a genome history made of duplications, losses, gene translocations, rearrangements, of 18 complete Anopheles genomes, including several important malaria vectors. We also provide additional support for a differentiated mode of evolution of the sex chromosome and of the autosomes in these mosquito genomes.

Conclusions

We demonstrate the method’s ability to improve extant assemblies accurately through a procedure simulating realistic assembly fragmentation. We study a debated issue regarding the phylogeny of the Gambiae complex group of Anopheles genomes in the light of the evolution of chromosomal rearrangements, suggesting that the phylogenetic signal they carry can differ from the phylogenetic signal carried by gene sequences, more prone to introgression.

     

Evolution of genome size in the angiosperms

Genome size varies extensively across the flowering plants, which has stimulated speculation regarding the ancestral genome size of these plants and trends in genome evolution. We investigated the evolution of C-values across the angiosperms using a molecular phylogenetic framework and C-values not previously available for crucial basal angiosperms, including Amborella, Illiciaceae, and Austrobaileya. Reconstructions of genome size across the angiosperms and extant gymnosperms indicate that the ancestral genome size for angiosperms is very small (1C ≤ 1.4 pg), in agreement with an earlier analysis of Leitch et al. (1998). Furthermore, a very small genome size (1C ≤ 1.4 pg) is ancestral not only for the angiosperms in general, but also for most major clades of flowering plants, including the monocots and the eudicots. The ancestral genome of core eudicots may also have been very small given that very low 1C-values appear to be ancestral for major clades of core eudicots, such as Caryophyllales, Saxifragales, and asterids. Very large genomes occur in clades that occupy derived positions within the monocots and Santalales.     

Shrinking genomes? Evidence from genome size variation in Crepis (Compositae)

Large-scale surveys of genome size evolution in angiosperms show that the ancestral genome was most likely small, with a tendency towards an increase in DNA content during evolution. Due to polyploidisation and self-replicating DNA elements, angiosperm genomes were considered to have a 'one-way ticket to obesity' (Bennetzen & Kellogg 1997). New findings on how organisms can lose DNA challenged the hypotheses of unidirectional evolution of genome size. The present study is based on the classical work of Babcock (1947a) on karyotype evolution within Crepis and analyses karyotypic diversification within the genus in a phylogenetic context. Genome size of 21 Crepis species was estimated using flow cytometry. Additional data of 17 further species were taken from the literature. Within 30 diploid Crepis species there is a striking trend towards genome contraction. The direction of genome size evolution was analysed by reconstructing ancestral character states on a molecular phylogeny based on ITS sequence data. DNA content is correlated to distributional aspects as well as life form. Genome size is significantly higher in perennials than in annuals. Within sampled species, very small genomes are only present in Mediterranean or European species, whereas their Central and East Asian relatives have larger 1C values.     

Chromosome evolution in perissodactyla

Comparative painting has provided a wealth of useful information and helped to reconstruct the pathways of karyotype evolution within major eutherian phylogenetic clades. New data have come from gene localizations, BAC mapping and high throughout sequencing projects that enrich and provide new details of genome evolution. Extensive research on perissodactyl genomes has revealed not only increased rates of chromosomal rearrangements, but also an exceptionally high number of centromere repositioning events in equids. Here were combined new physical mapping, comparative painting and genome sequencing data to refine the putative ancestral karyotype maps and to revise the previously proposed scenario of perissodactyl karyotype evolution.     

Genome trees from conservation profiles

The concept of the genome tree depends on the potential evolutionary significance in the clustering of species according to similarities in the gene content of their genomes. In this respect, genome trees have often been identified with species trees. With the rapid expansion of genome sequence data it becomes of increasing importance to develop accurate methods for grasping global trends for the phylogenetic signals that mutually link the various genomes. We therefore derive here the methodological concept of genome trees based on protein conservation profiles in multiple species. The basic idea in this derivation is that the multi-component "presence-absence" protein conservation profiles permit tracking of common evolutionary histories of genes across multiple genomes. We show that a significant reduction in informational redundancy is achieved by considering only the subset of distinct conservation profiles. Beyond these basic ideas, we point out various pitfalls and limitations associated with the data handling, paving the way for further improvements. As an illustration for the methods, we analyze a genome tree based on the above principles, along with a series of other trees derived from the same data and based on pair-wise comparisons (ancestral duplication-conservation and shared orthologs). In all trees we observe a sharp discrimination between the three primary domains of life: Bacteria, Archaea, and Eukarya. The new genome tree, based on conservation profiles, displays a significant correspondence with classically recognized taxonomical groupings, along with a series of departures from such conventional clusterings.     

Reconstructing an ancestral mammalian immune supercomplex from a marsupial major histocompatibility complex

The first sequenced marsupial genome promises to reveal unparalleled insights into mammalian evolution. We have used theMonodelphis domestica (gray short-tailed opossum) sequence to construct the first map of a marsupial major histocompatibility complex (MHC). The MHC is the most gene-dense region of the mammalian genome and is critical to immunity and reproductive success. The marsupial MHC bridges the phylogenetic gap between the complex MHC of eutherian mammals and the minimal essential MHC of birds. Here we show that the opossum MHC is gene dense and complex, as in humans, but shares more organizational features with non-mammals. The Class I genes have amplified within the Class II region, resulting in a unique Class I/II region. We present a model of the organization of the MHC in ancestral mammals and its elaboration during mammalian evolution. The opossum genome, together with other extant genomes, reveals the existence of an ancestral “immune supercomplex” that contained genes of both types of natural killer receptors together with antigen processing genes and MHC genes.     

Polyploidy and angiosperm diversification

Polyploidy has long been recognized as a major force in angiosperm evolution. Recent genomic investigations not only indicate that polyploidy is ubiquitous among angiosperms, but also suggest several ancient genome-doubling events. These include ancient whole genome duplication (WGD) events in basal angiosperm lineages, as well as a proposed paleohexaploid event that may have occurred close to the eudicot divergence. However, there is currently no evidence for WGD in Amborella, the putative sister species to other extant angiosperms. The question is no longer "What proportion of angiosperms are polyploid?", but "How many episodes of polyploidy characterize any given lineage?" New algorithms provide promise that ancestral genomes can be reconstructed for deep divergences (e.g., it may be possible to reconstruct the ancestral eudicot or even the ancestral angiosperm genome). Comparisons of diversification rates suggest that genome doubling may have led to a dramatic increase in species richness in several angiosperm lineages, including Poaceae, Solanaceae, Fabaceae, and Brassicaceae. However, additional genomic studies are needed to pinpoint the exact phylogenetic placement of the ancient polyploidy events within these lineages and to determine when novel genes resulting from polyploidy have enabled adaptive radiations.     

Genome size evolution of the extant lycophytes and ferns

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.     

Use of genomic history to improve phylogeny and understanding of births and deaths in a gene family

Polyploidy events have played an important role in the evolution of angiosperm genomes. Here, we demonstrate how genomic histories can increase phylogenetic resolution in a gene family, specifically the expansin superfamily of cell wall proteins. There are 36 expansins in Arabidopsis and 58 in rice. Traditional sequence-based phylogenetic trees yield poor resolution below the family level. To improve upon these analyses, we searched for gene colinearity (microsynteny) between Arabidopsis and rice genomic segments containing expansin genes. Multiple rounds of genome duplication and extensive gene loss have obscured synteny. However, by simultaneously aligning groups of up to 10 potentially orthologous segments from the two species, we traced the history of 49 out of 63 expansin-containing segments back to the ancestor of monocots and eudicots. Our results indicate that this ancestor had 15-17 expansin genes, each ancestral to an extant clade. Some clades have strikingly different growth patterns in the rice and Arabidopsis lineages, with more than half of all rice expansins arising from two ancestral genes. Segmental duplications, most of them part of polyploidy events, account for 12 out of 21 new expansin genes in Arabidopsis and 16 out of 44 in rice. Tandem duplications explain most of the rest. We were also able to estimate a minimum of 28 gene deaths in the Arabidopsis lineage and nine in rice. This analysis greatly clarifies expansin evolution since the last common ancestor of monocots and eudicots and the method should be broadly applicable to many other gene families.     

The complete mitochondrial DNA sequence of the green alga Pseudendoclonium akinetum (Ulvophyceae) highlights distinctive evolutionary trends in the chlorophyta and suggests a sister-group relationship between the Ulvophyceae and Chlorophyceae

The mitochondrial genome has undergone radical changes in both the Chlorophyta and Streptophyta, yet little is known about the dynamics of mtDNA evolution in either of these lineages. In the Chlorophyta, which comprises four of the five recognized classes of green algae (Prasinophyceae, Trebouxiophyceae, Ulvophyceae, and Chlorophyceae), the mitochondrial genome varies from 16 to 55 kb. This genome has retained a compact gene organization and a relatively complex gene repertoire ("ancestral" pattern) in the basal lineages represented by the Trebouxiophyceae and Prasinophyceae, whereas it has been reduced in size and gene complement and tends to evolve much more rapidly at the sequence level ("reduced-derived" pattern of evolution) in the Chlorophyceae and the lineage leading to the enigmatic chlorophyte Pedinomonas. To gain information about the evolutionary trends of mtDNA in the Ulvophyceae and also to gain insights into the phylogenetic relationships between ulvophytes and other chlorophytes, we have determined the mtDNA sequence of Pseudendoclonium akinetum. At 95,880 bp, Pseudendoclonium mtDNA is the largest green-algal mitochondrial genome sequenced to date and has the lowest gene density. These derived features are reminiscent of the "expanded" pattern exhibited by embryophyte mtDNAs, indicating that convergent evolution towards genome expansion has occurred independently in the Chlorophyta and Streptophyta. With 57 conserved genes, the gene repertoire of Pseudendoclonium mtDNA is slightly smaller than those of the prasinophyte Nephroselmis olivacea and the trebouxiophyte Prototheca wickerhamii. This ulvophyte mtDNA contains seven group I introns, four of which have homologs in green-algal mtDNAs displaying an "ancestral" or a "reduced-derived" pattern of evolution. Like its counterpart in the chlorophycean green alga Scenedesmus obliquus, it features numerous small, dispersed repeats in intergenic regions and introns. Its overall rate of sequence evolution appears to be accelerated to an intermediary level as compared with the rates observed in "ancestral" and "reduced-derived" mtDNAs. In agreement with the finding that Pseudendoclonium mtDNA exhibits features typical of both the "ancestral" and "reduced-derived" patterns of evolution, phylogenetic analyses of seven mtDNA-encoded proteins revealed a sister-group relationship between this ulvophyte and chlorophytes displaying "reduced-derived" mtDNAs.     

Spontaneous symmetry breaking in genome evolution

The quest for evolutionary mechanisms providing separation between the coding (exons) and noncoding (introns) parts of genomic DNA remains an important focus of genetics. This work combines an analysis of the most recent achievements of genomics and fundamental concepts of random processes to provide a novel point of view on genome evolution. Exon sizes in sequenced genomes show a lognormal distribution typical of a random Kolmogoroff fractioning process. This implies that the process of intron incretion may be independent of exon size, and therefore could be dependent on intron-exon boundaries. All genomes examined have two distinctive classes of exons, each with different evolutionary histories. In the framework proposed in this article, these two classes of exons can be derived from a hypothetical ancestral genome by (spontaneous) symmetry breaking. We note that one of these exon classes comprises mostly alternatively spliced exons.     

On the need for integrative phylogenomics,and some steps toward its creation

Recently improved understanding of evolutionary processes suggests that tree-based phylogenetic analyses of evolutionary change cannot adequately explain the divergent evolutionary histories of a great many genes and gene complexes. In particular, genetic diversity in the genomes of prokaryotes, phages, and plasmids cannot be fit into classic tree-like models of evolution. These findings entail the need for fundamental reform of our understanding of molecular evolution and the need to devise alternative apparatus for integrated analysis of these genomes. We advocate the development of integrative phylogenomics for analyzing these genomes and their histories, with tools suited to analyzing the importance of lateral gene transfer (LGT) and of DNA evolution in extra-cellular mobile genetic elements (e.g., viruses, plasmids). These phenomena greatly increase the complexity of relationships among interacting genetic partners, as they exchange functional genetic units. We examine the ontology of functional genetic units, interacting genetic partners, and emergent genetic associations, argue that these three categories of entities are required for a successful integrated phylogenomics. We conclude with arguments to suggest that the proposed new perspective and associated tools are suitable, and perhaps required, as a replacement for the bifurcating trees that have dominated evolutionary thinking for the last 150 years.     

Evolution of photosynthetic prokaryotes: a maximum-likelihood mapping approach

Reconstructing the early evolution of photosynthesis has been guided in part by the geological record, but the complexity and great antiquity of these early events require molecular genetic techniques as the primary tools of inference. Recent genome sequencing efforts have made whole genome data available from representatives of each of the five phyla of bacteria with photosynthetic members, allowing extensive phylogenetic comparisons of these organisms. Here, we have undertaken whole genome comparisons using maximum likelihood to compare 527 unique sets of orthologous genes from all five photosynthetic phyla. Substantiating recent whole genome analyses of other prokaryotes, our results indicate that horizontal gene transfer (HGT) has played a significant part in the evolution of these organisms, resulting in genomes with mosaic evolutionary histories. A small plurality phylogenetic signal was observed, which may be a core of remnant genes not subject to HGT, or may result from a propensity for gene exchange between two or more of the photosynthetic organisms compared.     

Untangling complex histories of genome mergings in high polyploids

Polyploidy, the duplication of entire genomes, plays a major role in plant evolution. In allopolyploids, genome duplication is associated with hybridization between two or more divergent genomes. Successive hybridization and polyploidization events can build up species complexes of allopolyploids with complicated network-like histories, and the evolutionary history of many plant groups cannot be adequately represented by phylogenetic trees because of such reticulate events. The history of complex genome mergings within a high-polyploid species complex in the genus Cerastium (Caryophyllaceae) is here untangled by the use of a network algorithm and noncoding sequences of a low-copy number gene. The resulting network illustrates how hybridization and polyploidization have acted as key evolutionary processes in creating a plant group where high-level allopolyploids clearly outnumber extant parental genomes.     

Evolution of plant mitochondrial genomes via substoichiometric intermediates

Comparison of the modern fertile maize mitochondrial genome (N) with an ancestral maize mitochondrial genome (RU) reveals a 12 kb duplication (containing the atpA gene) in the modern genome that is absent from the ancestor. Cloning, mapping, and sequencing of the relevant portions of the ancestral genome shows that this duplication probably arose via a three-stage recombination process involving substoichiometric intermediates. Comparison with analogous observations on yeast mitochondrial genomes suggests that this three-stage model of genome reorganization can be generally applied to plant mitochondrial genomes to explain both deletions and the creation of novel repeats, common features of plant mitochondrial genome evolution.     

Deciphering the evolutionary interplay between subgenomes following polyploidy: A paleogenomics approach in grasses

How did plant species emerge from their most recent common ancestors (MRCAs) 250 million years ago? Modern plant genomes help to address such key questions in unveiling precise species genealogies. The field of paleogenomics is undergoing a paradigm shift for investigating species evolution from the study of ancestral genomes from extinct species to deciphering the evolutionary forces (in terms of duplication, fusion, fission, deletion, and translocation) that drove present‐day plant diversity (in terms of chromosome/gene number and genome size). In this review, inferred ancestral karyotype genomes are shown to be powerful tools to (1) unravel the past history of extant species by recovering the variations of ancestral genomic compartments and (2) accelerate translational research by facilitating the transfer of genomic information from model systems to species of agronomic interest.     

L1 (LINE-1) retrotransposable elements provide a "fossil" record of the phylogenetic history of murid rodents

The single most difficult problem in phylogenetic analysis is deciding whether a shared taxonomic character is due to common ancestry or one that appeared independently due to convergence, parallelism, or reversion to an ancestral state. Mammalian L1 retrotransposons undergo periodic amplifications in which multiple copies of the elements are interspersed in the genome. Because these elements apparently are transmitted only by inheritance and are retained in the genome, a shared L1 amplification event can only be an inherited ancestral character. We propose that L1 amplification events can be an excellent tool for analyzing mammalian evolution and demonstrate here how we addressed several refractory problems in rodent systematics using L1 DNA as a taxonomic character.      

Evolution of parental ITS regions of nuclear rDNA in allopolyploid Aegilops (Poaceae) species

The genus Aegilops comprises approximately 25 diploid, tetraploid and hexaploid species, in which the genome types of all allopolyploids involve either U or D genome, or both of them. The internal transcribed spacer (ITS) region of 18S-26S nuclear ribosomal DNA (rDNA) from 11 allopolyploid species and 7 related diploid species in the genus were directly sequenced by pooled PCR products. Phylogenetic analyses for tracing evolutionary patterns of parental rDNA in allopolyploid species were performed using the neighbor-joining method. The D genome involved tree included three clades (CC-DDCC, DDMM-DDMMSS-DDMMUU, and MM-MhMh-DDNN), but did not include Ae. squarrosa (DD). It indicated that the rDNA of ancestral D genome had been somewhat differentiated in allopolyploids. The U genome involved tree showed that the allopolyploids and their common ancestor, Ae. umbellulata, formed a clade, suggesting that rDNA in UUMM and UUSS genomes has been homogenizing toward that of ancestral U genome. The phylogenetic pattern of U genome based on ITS sequences also supported the "pivotal-differential" hypothesis.     

Mosaic genomes of the six major primate lentivirus lineages revealed by phylogenetic analyses

To clarify the origin and evolution of the primate lentiviruses (PLVs), which include human immunodeficiency virus types 1 and 2 as well as their simian relatives, simian immunodeficiency viruses (SIVs), isolated from several host species, we investigated the phylogenetic relationships among the six supposedly nonrecombinant PLV lineages for which the full genome sequences are available. Employing bootscanning as an exploratory tool, we located several regions in the PLV genome that seem to have uncertain or conflicting phylogenetic histories. Phylogeny reconstruction based on distance and maximum-likelihood algorithms followed by a number of statistical tests confirms the existence of at least five putative recombinant fragments in the PLV genome with different clustering patterns. Split decomposition analysis also shows that phylogenetic relationships among PLVs may be better represented by network-based graphs, such as the ones produced by SplitsTree. Our findings not only imply that the six so-called pure PLV lineages have in fact mosaic genomes but also make more unlikely the hypothesis of cospeciation of SIVs and their simian hosts.     

Effects of Life Histories on Genome Size Variation in Squamata

Genome size changes significantly among taxonomic levels, and this variation is often related to the patterns shaped by the phylogeny, life histories and ecological factors. However, there are mixed evidences on the main factors affecting molecular evolution in animals.In this study, we used phylogenetic comparative analysis to investigate the evolutionary rate of genome size and the relationships between genome size and life histories(i.e.,hatchling mass, clutch size, clutches per year, age at sexual maturity, lifespan and body mass) among 199 squamata species. Our results showed that the evolutionary rate of genome size in Lacertilia was significantly faster than Serpentes. Moreover, we also found that larger species showed larger hatchling mass, more clutches per year and clutch size and longer lifespan. However, genome size was negatively associated with clutch size and clutches per year, but not associated with body mass we looked at.The findings suggest that larger species do not possess the evolution of large genomes in squamata.