A Population-Genetic Lens into the Process of Gene Loss Following Whole-Genome Duplication

Whole-genome duplications (WGDs) have occurred in many eukaryotic lineages. However, the underlying evolutionary forces and molecular mechanisms responsible for the long-term retention of gene duplicates created by WGDs are not well understood. We employ a population-genomic approach to understand the selective forces acting on paralogs and investigate ongoing duplicate-gene loss in multiple species of Paramecium that share an ancient WGD. We show that mutations that abolish protein function are more likely to be segregating in retained WGD paralogs than in single-copy genes, most likely because of ongoing nonfunctionalization post-WGD. This relaxation of purifying selection occurs in only one WGD paralog, accompanied by the gradual fixation of nonsynonymous mutations and reduction in levels of expression, and occurs over a long period of evolutionary time, “marking” one locus for future loss. Concordantly, the fitness effects of new nonsynonymous mutations and frameshift-causing indels are significantly more deleterious in the highly expressed copy compared with their paralogs with lower expression. Our results provide a novel mechanistic model of gene duplicate loss following WGDs, wherein selection acts on the sum of functional activity of both duplicate genes, allowing the two to wander in expression and functional space, until one duplicate locus eventually degenerates enough in functional efficiency or expression that its contribution to total activity is too insignificant to be retained by purifying selection. Retention of duplicates by such mechanisms predicts long times to duplicate-gene loss, which should not be falsely attributed to retention due to gain/change in function.     

Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

A comparative genomics study of nine Paramecium species reveals successful invasion of genes by transposable elements in their germline genomes, showing that the internal eliminated sequences (IESs) followed an evolutionary trajectory remarkably similar to that of spliceosomal introns.     

Screening of duplicated loci reveals hidden divergence patterns in a complex salmonid genome

A whole‐genome duplication (WGD) doubles the entire genomic content of a species and is thought to have catalysed adaptive radiation in some polyploid‐origin lineages. However, little is known about general consequences of a WGD because gene duplicates (i.e., paralogs) are commonly filtered in genomic studies; such filtering may remove substantial portions of the genome in data sets from polyploid‐origin species. We demonstrate a new method that enables genome‐wide scans for signatures of selection at both nonduplicated and duplicated loci by taking locus‐specific copy number into account. We apply this method to RAD sequence data from different ecotypes of a polyploid‐origin salmonid (Oncorhynchus nerka) and reveal signatures of divergent selection that would have been missed if duplicated loci were filtered. We also find conserved signatures of elevated divergence at pairs of homeologous chromosomes with residual tetrasomic inheritance, suggesting that joint evolution of some nondiverged gene duplicates may affect the adaptive potential of these genes. These findings illustrate that including duplicated loci in genomic analyses enables novel insights into the evolutionary consequences of WGDs and local segmental gene duplications.     

Help or hindrance? The evolutionary impact of whole‐genome duplication on immunogenetic diversity and parasite load

Whole‐genome duplication (WGD) events occur in all kingdoms and have been hypothesized to promote adaptability. WGDs identified in the early history of vertebrates, teleosts, and angiosperms have been linked to the large‐scale diversification of these lineages. However, the mechanics and full outcomes of WGD regarding potential evolutionary impacts remain a topic of debate. The Corydoradinae are a diverse subfamily of Neotropical catfishes with over 170 species described and a history of WGDs. They are divided into nine mtDNA lineages, with species coexisting in sympatric—and often mimetic—communities containing representatives of two or more of the nine lineages. Given their similar life histories, coexisting species of Corydoras might be exposed to similar parasite loads and because of their different histories of WGD and genome size they provide a powerful system for investigating the impacts of WGD on immune diversity and function in an animal system. Here, we compared parasite counts and the diversity of the immune‐related toll‐like receptors (TLR) in two coexisting species of Corydoras catfish (C. maculifer and C. araguaiaensis), one diploid and one putative tetraploid. In the putative tetraploid C. araguaiaensis, we found significantly lower numbers of parasites and significantly higher diversity (measured by both synonymous and nonsynonymous SNP counts) in two TLR genes than in the diploid C. maculifer. These results provide insight into how WGD may impact evolution, in this case by providing greater immunogenetic diversity.     

Comparison of the Esterases and Acid Phosphatases in Paramecium multimicronucleatum,Syngens 1–5, P. jenningsi,P. caudatum,and the P. aurelia Complex1

ABSTRACT. Forty-eight stocks in Paramecium jenningsi, syngens 1–5 of P. multimicronucleatum, P. caudatum, P. primaurelia, P. biaurelia, and P. tetraurelia were grown axenically and tested for their esterases and acid phosphatases using starch gel electrophoresis. The five esterases and the acid phosphatases previously characterized in species of the P. aurelia complex were also found in P. jenningsi, and three to four of the esterases and the acid phosphatases were found in the P. multimicronucleatum species complex and in P. caudatum. Additional subtypes were observed for each of the enzyme phenotypes in these new (though here unnamed) species of Paramecium. Two of the new acid phosphatase subtypes, which depart radically in mobility and in pattern, were found in syngen 3 of P. multimicronucleatum and in P. caudatum. Except for syngens 1 and 5 in P. multimicronucleatum, the degree of similarity between syngens 1, 5 and 2, 3, and 4 appears to be very low—perhaps even lower than that seen for species in the aurelia complex. More realistically, the syngens of P. multimicronucleatum should be considered as separate species although they are not here given separate taxonomic names. Limited sharing of subtypes occurred between species in different species complexes. This observation suggests that the molecular distances between species complexes may be even greater than between species within a complex.     

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.     

Preferential retention of the slowly evolving gene in pairs of duplicates in angiosperm genomes

Gene duplication provides raw material for functional innovation, but gene duplicability varies considerably. Previous studies have found widespread asymmetrical sequence evolution between paralogs. However, it remains unknown whether the rate of evolution among paralogs affects their propensity of being retained after another round of whole-genome duplication (WGD). In this study, we investigated gene groups that have experienced two successive WGDs to determine which of two older duplicates with different evolutionary rates was more likely to retain both younger duplicates. To uncouple the measurement of evolutionary rates from any assignment of duplicate or singleton status, we measured the evolutionary rates of singleton genes in out-lineages but classified these singleton genes according to whether they are retained or not in a crown group of species. We found that genes that retained younger duplicates in the crown group of genomes were more constrained prior to the younger duplication event than those that failed to leave duplicates. In addition, we also found that the retained clades have more genes in out-lineages. Subsequent analyses showed that genes in the retained clades were expressed more broadly and highly than genes in the singleton clades. We concluded that the set of repeatedly retained genes after two WGDs is biased toward slowly evolving genes in angiosperms, suggesting that the potential of genes for both functional conservation and divergence likely affects their propensity of being retained after WGD in angiosperms.     

Intraspecific differentiation of Paramecium novaurelia strains (Ciliophora,Protozoa) inferred from phylogenetic analysis of ribosomal and mitochondrial DNA variation

Paramecium novaurelia Beale and Schneller, 1954, was first found in Scotland and is known to occur mainly in Europe, where it is the most common species of the P. aurelia complex. In recent years, two non-European localities have been described: Turkey and the United States of America. This article presents the analysis of intraspecific variability among 25 strains of P. novaurelia with the application of ribosomal and mitochondrial loci (ITS1-5.8S-ITS2, 5′ large subunit rDNA (5′LSU rDNA) and cytochrome c oxidase subunit 1 (COI) mtDNA). The mean distance observed for all of the studied P. novaurelia sequence pairs was p = 0.008/0.016/0.092 (ITS1-5.8S-ITS2/5′LSU rDNA/COI). Phylogenetic trees (NJ/MP/BI) based on a comparison of all of the analysed sequences show that the studied strains of P. novaurelia form a distinct clade, separate from the P. caudatum outgroup, and are divided into two clusters (A and B) and two branches (C and D). The occurrence of substantial genetic differentiation within P. novaurelia, confirmed by the analysed DNA fragments, indicates a rapid evolution of particular species within the Paramecium genus.     

Ancestral gene duplications in mosses characterized by integrated phylogenomic analyses

Mosses (Bryophyta) are a key group occupying an important phylogenetic position in land plant (embryophyte) evolution. The class Bryopsida represents the most diversified lineage, containing more than 95% of modern mosses, whereas other classes are species‐poor. Two branches with large numbers of gene duplications were elucidated by phylogenomic analyses, one in the ancestry of all mosses and another before the separation of the Bryopsida, Polytrichopsida, and Tetraphidopsida. The analysis of the phylogenetic progression of duplicated paralogs retained on genomic syntenic regions in the Physcomitrella patens genome confirmed that the whole‐genome duplication events WGD1 and WGD2 were re‐recognized as the ψ event and the Funarioideae duplication event, respectively. The ψ polyploidy event was tightly associated with the early diversification of Bryopsida, in the ancestor of Bryidae, Dicranidae, Timmiidae, and Funariidae. Together, four branches with large numbers of gene duplications were unveiled in the evolutionary past of P. patens. Gene retention patterns following the four large‐scale duplications in different moss lineages were analyzed and discussed. Recurrent significant retention of stress‐related genes may have contributed to their adaption to distinct ecological environments and the evolutionary success of this early‐diverging land plant lineage.     

An Experimental Study on the Food Chain Among Bacteria,Paramecium and Daphnia

A quantitative study is made of the feeding of Paramecium caudatum on bacteria and of the feeding of Daphnia pulex on Paramecium and bacteria. The growth yield of P. caudatum on a diet of an unidentified yellow-colored bacterium was about 55% on a dry weight basis. Paramecium would not grow if the concentration of bacteria is less than 1 . 10₇ cells ml. Several other bacterial species were tested. Paramecium could not grow successfully on all of them, but could on Rhodopseudomonas spheroides. Adults of Daphnia pulex could grow and reproduce on a diet of P. caudatum alone, but young grew successfully only when bacteria were present.     

Alternative Transposition Generates New Chimeric Genes and Segmental Duplications at the Maize p1 Locus

The maize Ac/Ds transposon family was the first transposable element system identified and characterized by Barbara McClintock. Ac/Ds transposons belong to the hAT family of class II DNA transposons. We and others have shown that Ac/Ds elements can undergo a process of alternative transposition in which the Ac/Ds transposase acts on the termini of two separate, nearby transposons. Because these termini are present in different elements, alternative transposition can generate a variety of genome alterations such as inversions, duplications, deletions, and translocations. Moreover, Ac/Ds elements transpose preferentially into genic regions, suggesting that structural changes arising from alternative transposition may potentially generate chimeric genes at the rearrangement breakpoints. Here we identified and characterized 11 independent cases of gene fusion induced by Ac alternative transposition. In each case, a functional chimeric gene was created by fusion of two linked, paralogous genes; moreover, each event was associated with duplication of the ∼70-kb segment located between the two paralogs. An extant gene in the maize B73 genome that contains an internal duplication apparently generated by an alternative transposition event was also identified. Our study demonstrates that alternative transposition-induced duplications may be a source for spontaneous creation of diverse genome structures and novel genes in maize.     

The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment

Entodinium caudatum is an anaerobic binucleated ciliate representing the most dominant protozoal species in the rumen. However, its biological features are largely unknown due to the inability to establish an axenic culture. In this study, we primally sequenced its macronucleus (MAC) genome to aid the understanding of its metabolism, physiology, ecology. We isolated the MAC of E. caudatum strain MZG-1 and sequenced the MAC genome using Illumina MiSeq, MinION, and PacBio RSII systems. De novo assembly of the MiSeq sequence reads followed with subsequent scaffolding with MinION and PacBio reads resulted in a draft MAC genome about 117 Mbp. A large number of carbohydrate-active enzymes were likely acquired through horizontal gene transfer. About 8.74% of the E. caudatum predicted proteome was predicted as proteases. The MAC genome of E. caudatum will help better understand its important roles in rumen carbohydrate metabolism, and interaction with other members of the rumen microbiome.     

Did Gause Have a Yeast Infection?

We planned to develop predator–prey models using Paramecium and yeast, but they have not been empirically examined since work by Gause in the 1930s. Therefore, we evaluated if Paramecium aurelia ingests and grows on eight yeasts. Recognising that it ingested yeasts but could not grow, we assessed if it might grow on other yeasts, by empirically parameterising a predator–prey model that relies on ingestion, not growth. Simulations were compared to P. aurelia‐yeast time‐series data, from Gause. We hypothesised that if the model simulated predator–prey dynamics that mimicked the original data, then possibly P. aurelia could grow on yeast; simulations did not mimic the original data. Reviewing works by Gause exposed two issues: experiments were undoubtedly contaminated with bacteria, allowing growth on bacteria, not yeast; and the population cycle data cannot be considered a self‐sustaining time series, as they were manipulated by adding yeast and ciliates. We conclude that past and future work should not rely on this system, for either empirical or theoretical evaluations. Finally, although we show that P. aurelia, P. caudatum, Euplotes patella, and Blepharisma sp. cannot grow on yeast, Tetrahymena pyriformis and Colpidium striatum can; these may provide models to explore predator–prey dynamics.     

Comparative analysis of Gossypium and Vitis genomes indicates genome duplication specific to the Gossypium lineage

Genetic mapping studies have suggested that diploid cotton (Gossypium) might be an ancient polyploid. However, further evidence is lacking due to the complexity of the genome and the lack of sequence resources. Here, we used the grape (Vitis vinifera) genome as an out-group in two different approaches to further explore evidence regarding ancient genome duplication (WGD) event(s) in the diploid Gossypium lineage and its (their) effects: a genome-level alignment analysis and a local-level sequence component analysis. Both studies suggest that at least one round of genome duplication occurred in the Gossypium lineage. Also, gene densities in corresponding regions from Gossypium raimondii, V. vinifera, Arabidopsis thaliana and Carica papaya genomes are similar, despite the huge difference in their genome sizes and the different number of WGDs each genome has experienced. These observations fit the model that differences in plant genome sizes are largely explained by transposon insertions into heterochromatic regions.