Turning the clock back on ancient genome duplication

A De Novo Genome Sequence Assembly of the Arabidopsis thaliana Accession Niederzenz-1 Displays Presence/Absence Variation and Strong Synteny

Arabidopsis thaliana is the most important model organism for fundamental plant biology. The genome diversity of different accessions of this species has been intensively studied, for example in the 1001 genome project which led to the identification of many small nucleotide polymorphisms (SNPs) and small insertions and deletions (InDels). In addition, presence/absence variation (PAV), copy number variation (CNV) and mobile genetic elements contribute to genomic differences between A. thaliana accessions. To address larger genome rearrangements between the A. thaliana reference accession Columbia-0 (Col-0) and another accession of about average distance to Col-0, we created a de novo next generation sequencing (NGS)-based assembly from the accession Niederzenz-1 (Nd-1). The result was evaluated with respect to assembly strategy and synteny to Col-0. We provide a high quality genome sequence of the A. thaliana accession (Nd-1, LXSY01000000). The assembly displays an N50 of 0.590 Mbp and covers 99% of the Col-0 reference sequence. Scaffolds from the de novo assembly were positioned on the basis of sequence similarity to the reference. Errors in this automatic scaffold anchoring were manually corrected based on analyzing reciprocal best BLAST hits (RBHs) of genes. Comparison of the final Nd-1 assembly to the reference revealed duplications and deletions (PAV). We identified 826 insertions and 746 deletions in Nd-1. Randomly selected candidates of PAV were experimentally validated. Our Nd-1 de novo assembly allowed reliable identification of larger genic and intergenic variants, which was difficult or error-prone by short read mapping approaches alone. While overall sequence similarity as well as synteny is very high, we detected short and larger (affecting more than 100 bp) differences between Col-0 and Nd-1 based on bi-directional comparisons. The de novo assembly provided here and additional assemblies that will certainly be published in the future will allow to describe the pan-genome of A. thaliana.

The E-Subgroup Pentatricopeptide Repeat Protein Family in Arabidopsis thaliana and Confirmation of the Responsiveness PPR96 to Abiotic Stresses

Pentatricopeptide repeat (PPR) proteins are extensive in all eukaryotes. Their functions remain as yet largely unknown. Mining potential stress responsive PPRs, and checking whether known PPR editing factors are affected in the stress treatments. It is beneficial to elucidate the regulation mechanism of PPRs involved in biotic and abiotic stress. Here, we explored the characteristics and origin of the 105 E subgroup PPRs in Arabidopsis thaliana. Phylogenetic analysis categorized the E subgroup PPRs into five discrete groups (Cluster I to V), and they may have a common origin in both A. thaliana and rice. An in silico expression analysis of the 105 E subgroup PPRs in A. thaliana was performed using available microarray data. Thirty-four PPRs were differentially expressed during A. thaliana seed imbibition, seed development stage(s), and flowers development processes. To explore potential stress responsive PPRs, differential expression of 92 PPRs was observed in A. thaliana seedlings subjected to different abiotic stresses. qPCR data of E subgroup PPRs under stress conditions revealed that the expression of 5 PPRs was responsive to abiotic stresses. In addition, PPR96 is involved in plant responses to salt, abscisic acid (ABA), and oxidative stress. The T-DNA insertion mutation inactivating PPR96 expression results in plant insensitivity to salt, ABA, and oxidative stress. The PPR96 protein is localized in the mitochondria, and altered transcription levels of several stress-responsive genes under abiotic stress treatments. Our results suggest that PPR96 may important function in a role connecting the regulation of oxidative respiration and environmental responses in A. thaliana.

Evolutionary consequences, constraints and potential of polyploidy in plants

Polyploidy, the possession of more than 2 complete genomes, is a major force in plant evolution known to affect the genetic and genomic constitution and the phenotype of an organism, which will have consequences for its ecology and geography as well as for lineage diversification and speciation. In this review, we discuss phylogenetic patterns in the incidence of polyploidy including possible underlying causes, the role of polyploidy for diversification, the effects of polyploidy on geographical and ecological patterns, and putative underlying mechanisms as well as chromosome evolution and evolution of repetitive DNA following polyploidization. Spurred by technological advances, a lot has been learned about these aspects both in model and increasingly also in nonmodel species. Despite this enormous progress, long-standing questions about polyploidy still cannot be unambiguously answered, due to frequently idiosyncratic outcomes and insufficient integration of different organizational levels (from genes to ecology), but likely this will change in the near future. See also the sister article focusing on animals by Choleva and Janko in this themed issue.

Dramatic change in function and expression pattern of a gene duplicated by polyploidy created a paternal effect gene in the Brassicaceae

New gene formation by polyploidy has been an ongoing process during the evolution of various eukaryotes that has contributed greatly to the large number of genes in their genomes. After duplication, some genes that are retained can acquire new functions or expression patterns, or subdivide their functions or expression patterns between duplicates. Here, we show that SHORT SUSPENSOR (SSP) and Brassinosteroid Kinase 1 (BSK1) are paralogs duplicated by a polyploidy event that occurred in the Brassicaceae family about 23 Ma. SSP is involved in paternal control of zygote elongation in Arabidopsis thaliana by transcription in the sperm cells of pollen and then translation in the zygote, whereas BSK1 is involved in brassinosteroid signal transduction. Comparative analysis of expression in 63 different organs and developmental stages revealed that BSK1 and SSP have opposite expression patterns in pollen compared with all other parts of the plant. We determined that BSK1 retains the ancestral expression pattern and function. Thus, SSP has diverged in function after duplication from a component of the brassinosteroid signaling pathway to a paternal regulator of the timing of zygote elongation. The ancestral function of SSP was lost by deletions in the kinase domain. Our sequence rate analysis revealed that SSP but not BSK1 has experienced a greatly accelerated rate of amino acid sequence changes and relaxation of purifying selection. In addition, SSP has been duplicated to create a new gene (SSP-like1) with a completely different expression pattern, a shorter coding sequence that has lost a critical functional domain, and a greatly accelerated rate of amino acid sequence evolution along with evidence for positive selection, together indicative of neofunctionalization. This study illustrates two dramatic examples of neofunctionalization following gene duplication by complete changes in expression pattern and function. In addition, our findings indicate that paternal control of zygote elongation by SSP is an evolutionarily recent innovation in the Brassicaceae family.

Phylogenetic detection of numerous gene duplications shared by animals, fungi and plants

BACKGROUND

Gene duplication is considered a major driving force for evolution of genetic novelty, thereby facilitating functional divergence and organismal diversity, including the process of speciation. Animals, fungi and plants are major eukaryotic kingdoms and the divergences between them are some of the most significant evolutionary events. Although gene duplications in each lineage have been studied extensively in various contexts, the extent of gene duplication prior to the split of plants and animals/fungi is not clear.

RESULTS

Here, we have studied gene duplications in early eukaryotes by phylogenetic relative dating. We have reconstructed gene families (with one or more orthogroups) with members from both animals/fungi and plants by using two different clustering strategies. Extensive phylogenetic analyses of the gene families show that, among nearly 2,600 orthogroups identified, at least 300 of them still retain duplication that occurred before the divergence of the three kingdoms. We further found evidence that such duplications were also detected in some highly divergent protists, suggesting that these duplication events occurred in the ancestors of most major extant eukaryotic groups.

CONCLUSIONS

Our phylogenetic analyses show that numerous gene duplications happened at the early stage of eukaryotic evolution, probably before the separation of known major eukaryotic lineages. We discuss the implication of our results in the contexts of different models of eukaryotic phylogeny. One possible explanation for the large number of gene duplication events is one or more large-scale duplications, possibly whole genome or segmental duplication(s), which provides a genomic basis for the successful radiation of early eukaryotes.

The problem of the eukaryotic genome size

The current state of knowledge concerning the unsolved problem of the huge interspecific eukaryotic genome size variations not correlating with the species phenotypic complexity (C-value enigma also known as C-value paradox) is reviewed. Characteristic features of eukaryotic genome structure and molecular mechanisms that are the basis of genome size changes are examined in connection with the C-value enigma. It is emphasized that endogenous mutagens, including reactive oxygen species, create a constant nuclear environment where any genome evolves. An original quantitative model and general conception are proposed to explain the C-value enigma. In accordance with the theory, the noncoding sequences of the eukaryotic genome provide genes with global and differential protection against chemical mutagens and (in addition to the anti-mutagenesis and DNA repair systems) form a new, third system that protects eukaryotic genetic information. The joint action of these systems controls the spontaneous mutation rate in coding sequences of the eukaryotic genome. It is hypothesized that the genome size is inversely proportional to functional efficiency of the anti-mutagenesis and/or DNA repair systems in a particular biological species. In this connection, a model of eukaryotic genome evolution is proposed.

Persistence of duplicated PAC1 receptors in the teleost, Sparus auratus

Background:

Duplicated genes are common in vertebrate genomes. Their persistence is assumed to be either a consequence of gain of novel function (neofunctionalisation) or partitioning of the function of the ancestral molecule (sub-functionalisation). Surprisingly few studies have evaluated the extent of such modifications despite the numerous duplicated receptor and ligand genes identified in vertebrate genomes to date. In order to study the importance of function in the maintenance of duplicated genes, sea bream (Sparus auratus) PAC₁ receptors, sequence homologues of the mammalian receptor specific for PACAP (Pituitary Adenylate Cyclase-Activating Polypeptide), were studied. These receptors belong to family 2 GPCRs and most of their members are duplicated in teleosts although the reason why both persist in the genome is unknown.

Results:

Duplicate sea bream PACAP receptor genes (sbPAC₁A and sbPAC₁B), members of family 2 GPCRs, were isolated and share 77% amino acid sequence identity. RT-PCR with specific primers for each gene revealed that they have a differential tissue distribution which overlaps with the distribution of the single mammalian receptor. Furthermore, in common with mammals, the teleost genes undergo alternative splicing and a PAC₁Ahop1 isoform has been characterised. Duplicated orthologous receptors have also been identified in other teleost genomes and their distribution profile suggests that function may be species specific. Functional analysis of the paralogue sbPAC₁s in Cos7 cells revealed that they are strongly stimulated in the presence of mammalian PACAP₂₇ and PACAP₃₈ and far less with VIP (Vasoactive Intestinal Peptide). The sbPAC₁ receptors are equally stimulated (LOG_EC50 values for maximal cAMP production) in the presence of PACAP₂₇ (-8.74 ± 0.29 M and -9.15 ± 0.21 M, respectively for sbPAC₁A and sbPAC₁B, P > 0.05) and PACAP₃₈ (-8.54 ± 0.18 M and -8.92 ± 0.24 M, respectively for sbPAC₁A and sbPAC₁B, P > 0.05). Human VIP was found to stimulate sbPAC₁A (-7.23 ± 0.20 M) more strongly than sbPAC₁B (-6.57 ± 0.14 M, P < 0.05) and human secretin (SCT), which has not so far been identified in fish genomes, caused negligible stimulation of both receptors.

Conclusion:

The existence of functionally divergent duplicate sbPAC₁ receptors is in line with previously proposed theories about the origin and maintenance of duplicated genes. Sea bream PAC₁ duplicate receptors resemble the typical mammalian PAC₁, and PACAP peptides were found to be more effective than VIP in stimulating cAMP production, although sbPAC₁A was more responsive for VIP than sbPAC₁B. These results together with the highly divergent pattern of tissue distribution suggest that a process involving neofunctionalisation occurred after receptor duplication within the fish lineage and probably accounts for their persistence in the genome. The characterisation of further duplicated receptors and their ligands should provide insights into the evolution and function of novel protein-protein interactions associated with the vertebrate radiation.

Reconstruction of the vertebrate ancestral genome reveals dynamic genome reorganization in early vertebrates

Although several vertebrate genomes have been sequenced, little is known about the genome evolution of early vertebrates and how large-scale genomic changes such as the two rounds of whole-genome duplications (2R WGD) affected evolutionary complexity and novelty in vertebrates. Reconstructing the ancestral vertebrate genome is highly nontrivial because of the difficulty in identifying traces originating from the 2R WGD. To resolve this problem, we developed a novel method capable of pinning down remains of the 2R WGD in the human and medaka fish genomes using invertebrate tunicate and sea urchin genes to define ohnologs, i.e., paralogs produced by the 2R WGD. We validated the reconstruction using the chicken genome, which was not considered in the reconstruction step, and observed that many ancestral proto-chromosomes were retained in the chicken genome and had one-to-one correspondence to chicken microchromosomes, thereby confirming the reconstructed ancestral genomes. Our reconstruction revealed a contrast between the slow karyotype evolution after the second WGD and the rapid, lineage-specific genome reorganizations that occurred in the ancestral lineages of major taxonomic groups such as teleost fishes, amphibians, reptiles, and marsupials.

Evolution of secretin family GPCR members in the metazoa

Background

Comparative approaches using protostome and deuterostome data have greatly contributed to understanding gene function and organismal complexity. The family 2 G-protein coupled receptors (GPCRs) are one of the largest and best studied hormone and neuropeptide receptor families. They are suggested to have arisen from a single ancestral gene via duplication events. Despite the recent identification of receptor members in protostome and early deuterostome genomes, relatively little is known about their function or origin during metazoan divergence. In this study a comprehensive description of family 2 GPCR evolution is given based on in silico and expression analyses of the invertebrate receptor genes.

Results

Family 2 GPCR members were identified in the invertebrate genomes of the nematodes C. elegans and C. briggsae, the arthropods D. melanogaster and A. gambiae (mosquito) and in the tunicate C. intestinalis. This suggests that they are of ancient origin and have evolved through gene/genome duplication events. Sequence comparisons and phylogenetic analyses have demonstrated that the immediate gene environment, with regard to gene content, is conserved between the protostome and deuterostome receptor genomic regions. Also that the protostome genes are more like the deuterostome Corticotrophin Releasing Factor (CRF) and Calcitonin/Calcitonin Gene-Related Peptide (CAL/CGRP) receptors members than the other family 2 GPCR members. The evolution of family 2 GPCRs in deuterostomes is characterised by acquisition of new family members, with SCT (Secretin) receptors only present in tetrapods. Gene structure is characterised by an increase in intron number with organismal complexity with the exception of the vertebrate CAL/CGRP receptors.

Conclusion

The family 2 GPCR members provide a good example of gene duplication events occurring in tandem with increasing organismal complexity during metazoan evolution. The putative ancestral receptors are proposed to be more like the deuterostome CAL/CGRP and CRF receptors and this may be associated with their fundamental role in calcium regulation and the stress response, both of which are essential for survival.

Coding limits on the number of transcription factors

Background

Transcription factor proteins bind specific DNA sequences to control the expression of genes. They contain DNA binding domains which belong to several super-families, each with a specific mechanism of DNA binding. The total number of transcription factors encoded in a genome increases with the number of genes in the genome. Here, we examined the number of transcription factors from each super-family in diverse organisms.

Results

We find that the number of transcription factors from most super-families appears to be bounded. For example, the number of winged helix factors does not generally exceed 300, even in very large genomes. The magnitude of the maximal number of transcription factors from each super-family seems to correlate with the number of DNA bases effectively recognized by the binding mechanism of that super-family. Coding theory predicts that such upper bounds on the number of transcription factors should exist, in order to minimize cross-binding errors between transcription factors. This theory further predicts that factors with similar binding sequences should tend to have similar biological effect, so that errors based on mis-recognition are minimal. We present evidence that transcription factors with similar binding sequences tend to regulate genes with similar biological functions, supporting this prediction.

Conclusion

The present study suggests limits on the transcription factor repertoire of cells, and suggests coding constraints that might apply more generally to the mapping between binding sites and biological function.

Ancient duplicated conserved noncoding elements in vertebrates: a genomic and functional analysis

Fish-mammal genomic comparisons have proved powerful in identifying conserved noncoding elements likely to be cis-regulatory in nature, and the majority of those tested in vivo have been shown to act as tissue-specific enhancers associated with genes involved in transcriptional regulation of development. Although most of these elements share little sequence identity to each other, a small number are remarkably similar and appear to be the product of duplication events. Here, we searched for duplicated conserved noncoding elements in the human genome, using comparisons with Fugu to select putative cis-regulatory sequences. We identified 124 families of duplicated elements, each containing between two and five members, that are highly conserved within and between vertebrate genomes. In 74% of cases, we were able to assign a specific set of paralogous genes with annotation relating to transcriptional regulation and/or development to each family, thus removing much of the ambiguity in identifying associated genes. We find that duplicate elements have the potential to up-regulate reporter gene expression in a tissue-specific manner and that expression domains often overlap, but are not necessarily identical, between family members. Over two thirds of the families are conserved in duplicate in fish and appear to predate the large-scale duplication events thought to have occurred at the origin of vertebrates. We propose a model whereby gene duplication and the evolution of cis-regulatory elements can be considered in the context of increased morphological diversity and the emergence of the modern vertebrate body plan.

Diagnosing duplications – can it be done?

New genes arise through duplication and modification of DNA sequences on a range of scales: single gene duplication, duplication of large chromosomal fragments and whole-genome duplication. Each duplication mechanism has specific characteristics that influence the fate of the resulting duplicates, such as the size of the duplicated fragment, the potential for dosage imbalance, the preservation or disruption of regulatory control and genomic context. The ability to diagnose or identify the mechanism that produced a pair of paralogs has the potential to increase our ability to reconstruct evolutionary history, to understand the processes that govern genome evolution and to make functional predictions based on paralogy. The recent availability of large amounts of whole-genome sequence, often from several closely related species, has stimulated a wealth of new computational methods to diagnose gene duplications.

Yeast evolution and comparative genomics

Recent sequencing efforts and experiments have advanced our understanding of genome evolution in yeasts, particularly the Saccharomyces yeasts. The ancestral genome of the Saccharomyces sensu stricto complex has been subject to both whole-genome duplication, followed by massive sequence loss and divergence, and segmental duplication. In addition the subtelomeric regions are subject to further duplications and rearrangements via ectopic exchanges. Translocations and other gross chromosomal rearrangements that break down syntenic relationships occur; however, they do not appear to be a driving force of speciation. Analysis of single genomes has been fruitful for hypothesis generation such as the whole-genome duplication, but comparative genomics between close and more distant species has proven to be a powerful tool in testing these hypotheses as well as elucidating evolutionary processes acting on the genome. Future work on population genomics and experimental evolution will keep yeast at the forefront of studies in genome evolution.

Molecular evidence for asymmetric evolution of sister duplicated blocks after cereal polyploidy

Polyploidy (genome duplication) is thought to have contributed to the evolution of the eukaryotic genome, but complex genome structures and massive gene loss during evolution has complicated detection of these ancestral duplication events. The major factors determining the fate of duplicated genes are currently unclear, as are the processes by which duplicated genes evolve after polyploidy. Fine-scale analysis between homologous regions may allow us to better understand post-polyploidy evolution. Here, using gene-by-gene and gene-by-genome strategies, we identified the S5 region and four homologous regions within the japonica genome. Additional phylogenomic analyses of the comparable duplicated blocks indicate that four successive duplication events gave rise to these five regions, allowing us to propose a model for this local chromosomal evolution. According to this model, gene loss may play a major role in post-duplication genetic evolution at the segmental level. Moreover, we found molecular evidence that one of the sister duplicated blocks experienced more gene loss and a more rapid evolution subsequent to two recent duplication events. Given that these two recent duplication events were likely involved in polyploidy, this asymmetric evolution (gene loss and gene divergence) may be one possible mechanism accounting for the diploidization at the segmental level.

Genome evolution and biodiversity in teleost fish

Teleost fish, which roughly make up half of the extant vertebrate species, exhibit an amazing level of biodiversity affecting their morphology, ecology and behaviour as well as many other aspects of their biology. This huge variability makes fish extremely attractive for the study of many biological questions, particularly of those related to evolution. New insights gained from different teleost species and sequencing projects have recently revealed several peculiar features of fish genomes that might have played a role in fish evolution and speciation. There is now substantial evidence that a round of tetraploidization/rediploidization has taken place during the early evolution of the ray-finned fish lineage, and that hundreds of duplicate pairs generated by this event have been maintained over hundreds of millions of years of evolution. Differential loss or subfunction partitioning of such gene duplicates might have been involved in the generation of fish variability. In contrast to mammalian genomes, teleost genomes also contain multiple families of active transposable elements, which might have played a role in speciation by affecting hybrid sterility and viability. Finally, the amazing diversity of sex determination systems and the plasticity of sex chromosomes observed in teleost might have been involved in both pre- and postmating reproductive isolation. Comparison of data generated by current and future genome projects as well as complementary studies in other species will allow one to approach the molecular and evolutionary mechanisms underlying genome diversity in fish, and will certainly significantly contribute to our understanding of gene evolution and function in humans and other vertebrates.

Duplication and DNA segmental loss in the rice genome: implications for diploidization

* Large-scale duplication events have been recently uncovered in the rice genome, but different interpretations were proposed regarding the extent of the duplications. * Through analysing the 370 Mb genome sequences assembled into 12 chromosomes of Oryza sativa subspecies indica, we detected 10 duplicated blocks on all 12 chromosomes that contained 47% of the total predicted genes. Based on the phylogenetic analysis, we inferred that this was a result of a genome duplication that occurred c. 70 million years ago, supporting the polyploidy origin of the rice genome. In addition, a segmental duplication was also identified involving chromosomes 11 and 12, which occurred c. 5 million years ago. * Following the duplications, there have been large-scale chromosomal rearrangements and deletions. About 30-65% of duplicated genes were lost shortly after the duplications, leading to a rapid diploidization. * Together with other lines of evidence, we propose that polyploidization is still an ongoing process in grasses of polyploidy origins.

Computational approaches to unveiling ancient genome duplications

Recent analyses of complete genome sequences have revealed that many genomes have been duplicated in their evolutionary past. Such events have been associated with important biological transitions, major leaps in evolution and adaptive radiations of species. Here, we consider recently developed computational methods to detect such ancient large-scale gene duplication events. Several new approaches have been used to show that large-scale gene duplications are more common than previously thought.

Recent developments in computational approaches for uncovering genomic homology

Identifying genomic homology within and between genomes is essential when studying genome evolution. In the past years, different computational techniques have been developed to detect homology even when the actual similarity between homologous segments is low. Depending on the strategy used, these methods search for pairs of chromosomal segments between which either both gene content and order are conserved or gene content only. However, due to fact that, after their divergence, homologous segments can lose a different set of genes, these methods still often fail to detect genomic homology. Recently, more advanced approaches have been developed that can combine gene order and content information of multiple genomic segments.