全基因组复制事件的绝对定年揭示莲座蕨属植物的迟滞演化

全基因组复制在维管植物的物种形成过程中普遍存在, 被认为是物种适应极端环境的重要机制之一。确定全基因组复制事件的发生时间对理解生物的适应性演化具有重要意义。然而, 在维管植物, 特别是蕨类植物中, 全基因组复制事件的发生时间及其演化意义仍知之甚少。本研究以蕨类植物重要基部类群——福建莲座蕨(Angiopteris fokiensis)为例, 基于不同采样点(广东、广西、上海)的3个转录组学数据, 利用同义替换率(Ks)和绝对定年的方法分析全基因组复制事件的发生时间和物种单位时间内的分子演化速率, 并对事件发生后保留下的基因进行基因功能注释和富集分析。结果表明, 福建莲座蕨在159‒165 Mya发生了一次全基因组复制事件, 该复制事件优先保留的基因主要与营养代谢、信号传导、适应调节和组织结构生长相关。另外, 福建莲座蕨的分子演化速率为1.66 × 10^‒9 (同义替换/位点/年), 是除裸子植物外, 陆生植物中已知演化速率最缓慢的类群。综合以上研究结果, 我们推测福建莲座蕨全基因组复制的发生可能与裸子植物繁盛、核心被子植物集中兴起或托阿尔阶灭绝事件有关。而复制后显著保留基因可能促进了莲座蕨属(Angiopteris)植物的遗传和形态创新, 从而帮助其快速适应环境的剧烈变化。进一步对该类群植物演化速率缓慢的原因进行讨论, 推测莲座蕨属缓慢的演化速率可能与其本身世代周期长、基因组较大及其生长环境稳定有关。本研究通过分析福建莲座蕨的全基因组复制历史和复制基因的保留模式, 推测全基因组复制事件对促进演化速率较慢的植物适应极端环境变化具有重要意义, 可为理解其他陆生植物的适应性演化提供更多启发。     

Sox genes in grass carp (Ctenopharyngodon idella) with their implications for genome duplication and evolution

The Sox gene family is found in a broad range of animal taxa and encodes important gene regulatory proteins involved in a variety of developmental processes. We have obtained clones representing the HMG boxes of twelve Sox genes from grass carp (Ctenopharyngodon idella), one of the four major domestic carps in China. The cloned Sox genes belong to group B1, B2 and C. Our analyses show that whereas the human genome contains a single copy of Sox4, Sox11 and Sox14, each of these genes has two co-orthologs in grass carp, and the duplication of Sox4 and Sox11 occurred before the divergence of grass carp and zebrafish, which support the "fish-specific whole-genome duplication" theory. An estimation for the origin of grass carp based on the molecular clock using Sox1, Sox3 and Sox11 genes as markers indicates that grass carp (subfamily Leuciscinae) and zebrafish (subfamily Danioninae) diverged approximately 60 million years ago. The potential uses of Sox genes as markers in revealing the evolutionary history of grass carp are discussed.     

Diversification of non-TIR class NB-LRR genes in relation to whole-genome duplication events in Arabidopsis

Arabidopsis thaliana is believed to have experienced at least two and possibly three whole-genome duplication events in its evolutionary history. In order to investigate the evolutionary relationships between these duplication events and diversification of disease resistance (R) genes, segmental-duplication events containing R genes belonging to the nucleotide binding-leucine rich repeat (NB-LRR) class were identified. Of 153 segmental-duplication events containing NB-LRR genes, only 22 contained NB-LRR genes in both members of the duplication pair, indicating a high frequency of NB-LRR gene loss after whole-genome duplication. The relative age of the duplication events was estimated based on the average synonymous substitution rate of the duplicated gene pairs in the segments. These data were combined with phylogenetic analyses. NB-LRR genes present in segment pairs derived from the most recent whole-genome duplication event, estimated to have occurred only 20 to 40 million years ago, occupy very distant branches of the NB-LRR phylogenetic tree. These data suggest that when NB-LRR clusters are duplicated as part of a whole-genome duplication, homoeologous NB-LRR genes are preferentially lost, either by eliminating one copy of the cluster or by eliminating individual genes such that only paralogous NB-LRR genes are maintained.     

Fugu genome analysis provides evidence for a whole-genome duplication early during the evolution of ray-finned fishes

With about 24,000 extant species, teleosts are the largest group of vertebrates. They constitute more than 99% of the ray-finned fishes (Actinopterygii) that diverged from the lobe-finned fish lineage (Sarcopterygii) about 450 MYA. Although the role of genome duplication in the evolution of vertebrates is now established, its role in structuring the teleost genomes has been controversial. At least two hypotheses have been proposed: a whole-genome duplication in an ancient ray-finned fish and independent gene duplications in different lineages. These hypotheses are, however, based on small data sets and lack adequate statistical and phylogenetic support. In this study, we have made a systematic comparison of the draft genome sequences of Fugu and humans to identify paralogous chromosomal regions ("paralogons") in the Fugu that arose in the ray-finned fish lineage ("fish-specific"). We identified duplicate genes in the Fugu by phylogenetic analyses of the Fugu, human, and invertebrate sequences. Our analyses provide evidence for 425 fish-specific duplicate genes in the Fugu and show that at least 6.6% of the genome is represented by fish-specific paralogons. We estimated the ages of Fugu duplicate genes and paralogons using the molecular clock. Remarkably, the ages of duplicate genes and paralogons are clustered, with a peak around 350 MYA. These data strongly suggest a whole-genome duplication event early during the evolution of ray-finned fishes, probably before the origin of teleosts.     

Estimating the time to the whole-genome duplication and the duration of concerted evolution via gene conversion in yeast

A maximum-likelihood (ML) method is developed to estimate the duration of concerted evolution and the time to the whole-genome duplication (WGD) event in baker's yeast (Saccharomyces cerevisiae). The models with concerted evolution fit the data significantly better than the molecular clock model, indicating a crucial role of concerted evolution via gene conversion after gene duplication in yeast. Our ML estimate of the time to the WGD is nearly identical to the time to the speciation event between S. cerevisiae and Kluyveromyces waltii, suggesting that the WGD occurred in very early stages after speciation or the WGD might have been involved in the speciation event.     

Consistent patterns of rate asymmetry and gene loss indicate widespread neofunctionalization of yeast genes after whole-genome duplication

We investigated patterns of rate asymmetry in sequence evolution among the gene pairs (ohnologs) formed by whole-genome duplication (WGD) in yeast species. By comparing three species (Saccharomyces cerevisiae, Candida glabrata, and S. castellii) that underwent WGD to a nonduplicated outgroup (Kluyveromyces lactis), and by using a synteny framework to establish orthology and paralogy relationships at each duplicated locus, we show that 56% of ohnolog pairs show significantly asymmetric protein sequence evolution. For ohnolog pairs that remain duplicated in two species there is a strong tendency for the faster-evolving copy in one species to be orthologous to the faster copy in the other species, which indicates that the evolutionary rate differences were established before speciation and hence soon after the WGD. We also present evidence that in cases where one ohnolog has been lost from the genome of a post-WGD species, the lost copy was likely to have been the faster-evolving member of the pair prior to its loss. These results suggest that a significant fraction of the retained ohnologs in yeast species underwent neofunctionalization soon after duplication.     

Rapid evolution of goat and sheep globin genes following gene duplication

Statistical analyses of DNA sequences of globin genes (beta A, beta C, and gamma) from goat and sheep (including new sequence information for the second intron of sheep beta A and gamma, kindly provided by A. Davis and A. W. Nienhuis) indicate that the rates of nonsynonymous substitution in these genes have been greatly accelerated following the gene duplication separating gamma and the ancestor of beta A and beta C and the gene duplication separating beta A and beta C. In both cases the acceleration was apparently due to relaxation of purifying selection (functional constraints) rather than advantageous mutations because acceleration occurred only in less important parts of the beta globin chain. The rates of nonsynonymous substitution in these genes are estimated to be about 2.3 x 10(-9) per site per year, which is three times higher than that for the divergence between human beta and mouse beta major globin genes. Our analyses further suggest that the rate of synonymous substitution in functional genes and the rate of substitution in pseudogenes are approximately equal and are between 2.8 x 10(-9) and 5.0 x 10(-9) and that the rate of substitution in introns is about 3.0 x 10(-9). The divergence time between beta A and beta C and that between gamma and the beta A-beta C pair are about 12 and 30 million years, respectively. The proportion of transition mutations is estimated to be 64%, two times higher than expected under random mutation but considerably lower than the 96% estimated for animal mitochondrial DNA.      

Characterization of AmphiF-spondin reveals the modular evolution of chordate F-spondin genes

The F-spondin genes are a family of extracellular matrix molecules united by two conserved domains, FS1 and FS2, at the amino terminus plus a variable number of thrombospondin repeats at the carboxy terminus. Currently, characterized members include a single gene in Drosophila and multiple genes in vertebrates. The vertebrate genes are expressed in the midline of the developing embryo, primarily in the floor plate of the neural tube. To investigate the evolution of chordate F-spondin genes, I have used the basal position in chordate phylogeny of the acraniate amphioxus. A single F-spondin-related gene, named AmphiF-spondin, was isolated from amphioxus. Based on molecular phylogenetics, AmphiF-spondin is closely related to a particular subgroup of vertebrate F-spondin genes that encode six thrombospondin repeats. However, unlike these genes, expression of AmphiF-spondin is not confined to the midline but is found through most of the central nervous system. Additionally, AmphiF-spondin has lost three thrombospondin repeats and gained two fibronectin type III repeats, one of which has strong identity to a fibronectin type III repeat from Deleted in Colorectal Cancer (DCC). Taken together, these results suggest a complex evolutionary history for chordate F-spondin genes that includes (1) domain loss, (2) domain gain by tandem duplication and divergence of existing domains, and (3) gain of heterologous domains by exon shuffling.      

Ordered structure of the transcription network inherited from the yeast whole-genome duplication

Background  

Gene duplication, a major evolutionary path to genomic innovation, can occur at the scale of an entire genome. One such "whole-genome duplication" (WGD) event among the Ascomycota fungi gave rise to genes with distinct biological properties compared to small-scale duplications.     

The evolution of pepsinogen C genes in vertebrates: duplication, loss and functional diversification

Background

Aspartic proteases comprise a large group of enzymes involved in peptide proteolysis. This collection includes prominent enzymes globally categorized as pepsins, which are derived from pepsinogen precursors. Pepsins are involved in gastric digestion, a hallmark of vertebrate physiology. An important member among the pepsinogens is pepsinogen C (Pgc). A particular aspect of Pgc is its apparent single copy status, which contrasts with the numerous gene copies found for example in pepsinogen A (Pga). Although gene sequences with similarity to Pgc have been described in some vertebrate groups, no exhaustive evolutionary framework has been considered so far.

Methodology/Principal Findings

By combining phylogenetics and genomic analysis, we find an unexpected Pgc diversity in the vertebrate sub-phylum. We were able to reconstruct gene duplication timings relative to the divergence of major vertebrate clades. Before tetrapod divergence, a single Pgc gene tandemly expanded to produce two gene lineages (Pgbc and Pgc2). These have been differentially retained in various classes. Accordingly, we find Pgc2 in sauropsids, amphibians and marsupials, but not in eutherian mammals. Pgbc was retained in amphibians, but duplicated in the ancestor of amniotes giving rise to Pgb and Pgc1. The latter was retained in mammals and probably in reptiles and marsupials but not in birds. Pgb was kept in all of the amniote clade with independent episodes of loss in some mammalian species. Lineage specific expansions of Pgc2 and Pgbc have also occurred in marsupials and amphibians respectively. We find that teleost and tetrapod Pgc genes reside in distinct genomic regions hinting at a possible translocation.

Conclusions

We conclude that the repertoire of Pgc genes is larger than previously reported, and that tandem duplications have modelled the history of Pgc genes. We hypothesize that gene expansion lead to functional divergence in tetrapods, coincident with the invasion of terrestrial habitats.     

Comparative evolution of photosynthetic genes in response to polyploid and nonpolyploid duplication

The likelihood of duplicate gene retention following polyploidy varies by functional properties (e.g. gene ontologies or protein family domains), but little is known about the effects of whole-genome duplication on gene networks related by a common physiological process. Here, we examined the effects of both polyploid and nonpolyploid duplications on genes encoding the major functional groups of photosynthesis (photosystem I, photosystem II, the light-harvesting complex, and the Calvin cycle) in the cultivated soybean (Glycine max), which has experienced two rounds of whole-genome duplication. Photosystem gene families exhibit retention patterns consistent with dosage sensitivity (preferential retention of polyploid duplicates and elimination of nonpolyploid duplicates), whereas Calvin cycle and light-harvesting complex gene families do not. We observed similar patterns in barrel medic (Medicago truncatula), which shared the older genome duplication with soybean but has evolved independently for approximately 50 million years, and in Arabidopsis (Arabidopsis thaliana), which experienced two nested polyploidy events independent from the legume duplications. In both soybean and Arabidopsis, Calvin cycle gene duplicates exhibit a greater capacity for functional differentiation than do duplicates within the photosystems, which likely explains the greater retention of ancient, nonpolyploid duplicates and larger average gene family size for the Calvin cycle relative to the photosystems.     

Pioneering polyploids: the impact of whole-genome duplication on biome shifting in New Zealand Coprosma (Rubiaceae) and Veronica (Plantaginaceae)

The role of whole-genome duplication (WGD) in facilitating shifts into novel biomes remains unknown. Focusing on two diverse woody plant groups in New Zealand, Coprosma (Rubiaceae) and Veronica (Plantaginaceae), we investigate how biome occupancy varies with ploidy level, and test the hypothesis that WGD increases the rate of biome shifting. Ploidy levels and biome occupancy (forest, open and alpine) were determined for indigenous species in both clades. The distribution of low-ploidy (Coprosma: 2x, Veronica: 6x) versus high-ploidy (Coprosma: 4–10x, Veronica: 12–18x) species across biomes was tested statistically. Estimation of the phylogenetic history of biome occupancy and WGD was performed using time-calibrated phylogenies and the R package BioGeoBEARS. Trait-dependent dispersal models were implemented to determine support for an increased rate of biome shifting among high-ploidy lineages. We find support for a greater than random portion of high-ploidy species occupying multiple biomes. We also find strong support for high-ploidy lineages showing a three- to eightfold increase in the rate of biome shifts. These results suggest that WGD promotes ecological expansion into new biomes.     

Natural selection on gene order in the genome reorganization process after whole-genome duplication of yeast

A genome must locate its coding genes on the chromosomes in a meaningful manner with the help of natural selection, but the mechanism of gene order evolution is poorly understood. To explore the role of selection in shaping the current order of coding genes and their cis-regulatory elements, a comparative genomic approach was applied to the baker's yeast Saccharomyces cerevisiae and its close relatives. S. cerevisiae have experienced a whole-genome duplication followed by an extensive reorganization process of gene order, during which a number of new adjacent gene pairs appeared. We found that the proportion of new adjacent gene pairs in divergent orientation is significantly reduced, suggesting that such new divergent gene pairs may be disfavored most likely because their coregulation may be deleterious. It is also found that such new divergent gene pairs have particularly long intergenic regions. These observations suggest that selection specifically worked against deletions in intergenic regions of new divergent gene pairs, perhaps because they should be physically kept away so that they are not coregulated. It is indicated that gene regulation would be one of the major factors to determine the order of coding genes.     

On the origins of Mendelian disease genes in man: the impact of gene duplication

Over 3,000 human diseases are known to be linked to heritable genetic variation, mapping to over 1,700 unique genes. Dating of the evolutionary age of these disease-associated genes has suggested that they have a tendency to be ancient, specifically coming into existence with early metazoa. The approach taken by past studies, however, assumes that the age of a disease is the same as the age of its common ancestor, ignoring the fundamental contribution of duplication events in the evolution of new genes and function. Here, we date both the common ancestor and the duplication history of known human disease-associated genes. We find that the majority of disease genes (80%) are genes that have been duplicated in their evolutionary history. Periods for which there are more disease-associated genes, for example, at the origins of bony vertebrates, are explained by the emergence of more genes at that time, and the majority of these are duplicates inferred to have arisen by whole-genome duplication. These relationships are similar for different disease types and the disease-associated gene's cellular function. This indicates that the emergence of duplication-associated diseases has been ongoing and approximately constant (relative to the retention of duplicate genes) throughout the evolution of life. This continued until approximately 390 Ma from which time relatively fewer novel genes came into existence on the human lineage, let alone disease genes. For single-copy genes associated with disease, we find that the numbers of disease genes decreases with recency. For the majority of duplicates, the disease-associated mutation is associated with just one of the duplicate copies. A universal explanation for heritable disease is, thus, that it is merely a by-product of the evolutionary process; the evolution of new genes (de novo or by duplication) results in the potential for new diseases to emerge.