Degree of Functional Divergence in Duplicates Is Associated with Distinct Roles in Plant Evolution

Gene duplication is a major mechanism to create new genes. After gene duplication, some duplicated genes undergo functionalization, whereas others largely maintain redundant functions. Duplicated genes comprise various degrees of functional diversification in plants. However, the evolutionary fate of high and low diversified duplicates is unclear at genomic scale. To infer high and low diversified duplicates in Arabidopsis thaliana genome, we generated a prediction method for predicting whether a pair of duplicate genes was subjected to high or low diversification based on the phenotypes of knock-out mutants. Among 4,017 pairs of recently duplicated A. thaliana genes, 1,052 and 600 are high and low diversified duplicate pairs, respectively. The predictions were validated based on the phenotypes of generated knock-down transgenic plants. We determined that the high diversified duplicates resulting from tandem duplications tend to have lineage-specific functions, whereas the low diversified duplicates produced by whole-genome duplications are related to essential signaling pathways. To assess the evolutionary impact of high and low diversified duplicates in closely related species, we compared the retention rates and selection pressures on the orthologs of A. thaliana duplicates in two closely related species. Interestingly, high diversified duplicates resulting from tandem duplications tend to be retained in multiple lineages under positive selection. Low diversified duplicates by whole-genome duplications tend to be retained in multiple lineages under purifying selection. Taken together, the functional diversities determined by different duplication mechanisms had distinct effects on plant evolution.     

Fifteen Million Years of Evolution in the Oryza Genus Shows Extensive Gene Family Expansion

In analyzing gene families in the whole-genome sequences available for O. sativa （AA）, O. glaberrima （AA）, and O. brachyantha （FF）, we observed large size expansions in the AA genomes compared to FF genomes for the superfamilies F-box and NB-ARC, and five additional families： the Aspartic proteases, BTB/POZ proteins （BTB）, Glutaredoxins, Trypsin a-amylase inhibitor proteins, and Zf-Dof proteins. Their evolutionary dynamic was investigated to understand how and why such important size variations are observed between these closely related species. We show that expansions resulted from both amplification, largely by tandem duplications, and contraction by gene losses. For the F-box and NB-ARC gene families, the genes conserved in all species were under strong purifying selection while expanded orthologous genes were under more relaxed purifying selection. In F-box, NB-ARC, and BTB, the expanded groups were enriched in genes with little evidence of expression, in comparison with conserved groups. We also detected 87 loci under positive selection in the expanded groups. These results show that most of the duplicated copies in the expanded groups evolve neutrally after duplication because of functional redundancy but a fraction of these genes were preserved following neofunctionalization. Hence, the lineage-specific expansions observed between Oryza species were partly driven by directional selection.     

Eusociality Shapes Convergent Patterns of Molecular Evolution across Mitochondrial Genomes of Snapping Shrimps

Eusociality is a highly conspicuous and ecologically impactful behavioral syndrome that has evolved independently across multiple animal lineages. So far, comparative genomic analyses of advanced sociality have been mostly limited to insects. Here, we study the only clade of animals known to exhibit eusociality in the marine realm—lineages of socially diverse snapping shrimps in the genus Synalpheus. To investigate the molecular impact of sociality, we assembled the mitochondrial genomes of eight Synalpheus species that represent three independent origins of eusociality and analyzed patterns of molecular evolution in protein-coding genes. Synonymous substitution rates are lower and potential signals of relaxed purifying selection are higher in eusocial relative to noneusocial taxa. Our results suggest that mitochondrial genome evolution was shaped by eusociality-linked traits—extended generation times and reduced effective population sizes that are hallmarks of advanced animal societies. This is the first direct evidence of eusociality impacting genome evolution in marine taxa. Our results also strongly support the idea that eusociality can shape genome evolution through profound changes in life history and demography.     

重复基因的进化一回顾与进展

基因重复是普遍存在的生物学现象, 是基因组和遗传系统多样化的重要推动力量, 在生物进化过程中发挥着极其重要的作用。基因重复有何利弊, 基因发生重复后, 2个重复子拷贝的保留在基因功能方面是否存在偏好性, 子拷贝在表达和进化速率上如何分化, 以及重复基因为什么会被保留下来一直是进化生物学领域研究的热点问题之一。该文对以上重复基因研究的热点问题进行了介绍, 并对重复基因的进化机制和理论模型及其近年来的一些主要研究进展进行了综述。     

Chromosomal-level genome assembly of the orchid tree Bauhinia variegata (Leguminosae; Cercidoideae) supports the allotetraploid origin hypothesis of Bauhinia

Cercidoideae, one of the six subfamilies of Leguminosae, contains one genus Cercis with its chromosome number 2n = 14 and all other genera with 2n = 28. An allotetraploid origin hypothesis for the common ancestor of non-Cercis genera in this subfamily has been proposed; however, no chromosome-level genomes from Cercidoideae have been available to test this hypothesis. Here, we conducted a chromosome-level genome assembly of Bauhinia variegata to test this hypothesis. The assembled genome is 326.4 Mb with the scaffold N50 of 22.1 Mb and contains 37,996 protein-coding genes. The Ks distribution between gene pairs in the syntenic regions indicates two whole-genome duplications (WGDs): one is B. variegata-specific, and the other is shared among core eudicots. Although Ks between gene pairs generated by the recent WGD in Bauhinia is greater than that between Bauhinia and Cercis, the WGD was not detected in Cercis, which can be explained by an accelerated evolutionary rate in Bauhinia after divergence from Cercis. Ks distribution and phylogenetic analysis for gene pairs generated by the recent WGD in Bauhinia and their corresponding orthologs in Cercis support the allopolyploidy origin hypothesis of Bauhinia. The genome of B. variegata also provides a genomic resource for dissecting genetic basis of its ornamental traits.     

The evolutionary demography of duplicate genes

Although gene duplication has generally been viewed as a necessary source of material for the origin of evolutionary novelties, the rates of origin, loss, and preservation of gene duplicates are not well understood. Applying steady-state demographic techniques to the age distributions of duplicate genes censused in seven completely sequenced genomes, we estimate the average rate of duplication of a eukaryotic gene to be on the order of 0.01/gene/million years, which is of the same order of magnitude as the mutation rate per nucleotide site. However, the average half-life of duplicate genes is relatively small, on the order of 4.0 million years. Significant interspecific variation in these rates appears to be responsible for differences in species-specific genome sizes that arise as a consequence of a quasi-equilibrium birth-death process. Most duplicated genes experience a brief period of relaxed selection early in their history and a minority exhibit the signature of directional selection, but those that survive more than a few million years eventually experience strong purifying selection. Thus, although most theoretical work on the gene-duplication process has focused on issues related to adaptive evolution, the origin of a new function appears to be a very rare fate for a duplicate gene. A more significant role of the duplication process may be the generation of microchromosomal rearrangements through reciprocal silencing of alternative copies, which can lead to the passive origin of post-zygotic reproductive barriers in descendant lineages of incipient species.     

Insights into Three Whole-Genome Duplications Gleaned from the Paramecium caudatum Genome Sequence

Paramecium has long been a model eukaryote. The sequence of the Paramecium tetraurelia genome reveals a history of three successive whole-genome duplications (WGDs), and the sequences of P. biaurelia and P. sexaurelia suggest that these WGDs are shared by all members of the aurelia species complex. Here, we present the genome sequence of P. caudatum, a species closely related to the P. aurelia species group. P. caudatum shares only the most ancient of the three WGDs with the aurelia complex. We found that P. caudatum maintains twice as many paralogs from this early event as the P. aurelia species, suggesting that post-WGD gene retention is influenced by subsequent WGDs and supporting the importance of selection for dosage in gene retention. The availability of P. caudatum as an outgroup allows an expanded analysis of the aurelia intermediate and recent WGD events. Both the Guanine+Cytosine (GC) content and the expression level of preduplication genes are significant predictors of duplicate retention. We find widespread asymmetrical evolution among aurelia paralogs, which is likely caused by gradual pseudogenization rather than by neofunctionalization. Finally, cases of divergent resolution of intermediate WGD duplicates between aurelia species implicate this process acts as an ongoing reinforcement mechanism of reproductive isolation long after a WGD event.     

The temporal distribution of gene duplication events in a set of highly conserved human gene families

Using a data set of protein translations associated with map positions in the human genome, we identified 1520 mapped highly conserved gene families. By comparing sharing of families between genomic windows, we identified 92 potentially duplicated blocks in the human genome containing 422 duplicated members of these families. Using branching order in the phylogenetic trees, we timed gene duplication events in these families relative to the primate-rodent divergence, the amniote-amphibian divergence, and the deuterostome-protostome divergence. The results showed similar patterns of gene duplication times within duplicated blocks and outside duplicated blocks. Both within and outside duplicated blocks, numerous duplications were timed prior to the deuterostome-protostome divergence, whereas others occurred after the amniote-amphibian divergence. Thus, neither gene duplication in general nor duplication of genomic blocks could be attributed entirely to polyploidization early in vertebrate history. The strongest signal in the data was a tendency for intrachromosomal duplications to be more recent than interchromosomal duplications, consistent with a model whereby tandem duplication-whether of single genes or of genomic blocks-may be followed by eventual separation of duplicates due to chromosomal rearrangements. The rate of separation of tandemly duplicated gene pairs onto separated chromosomes in the human lineage was estimated at 1.7 x 10(-9) per gene-pair per year.     

Linked by Ancestral Bonds: Multiple Whole-Genome Duplications and Reticulate Evolution in a Brassicaceae Tribe

Pervasive hybridization and whole-genome duplications (WGDs) influenced genome evolution in several eukaryotic lineages. Although frequent and recurrent hybridizations may result in reticulate phylogenies, the evolutionary events underlying these reticulations, including detailed structure of the ancestral diploid and polyploid genomes, were only rarely reconstructed. Here, we elucidate the complex genomic history of a monophyletic clade from the mustard family (Brassicaceae), showing contentious relationships to the early-diverging clades of this model plant family. Genome evolution in the crucifer tribe Biscutelleae (∼60 species, 5 genera) was dominated by pervasive hybridizations and subsequent genome duplications. Diversification of an ancestral diploid genome into several divergent but crossable genomes was followed by hybridizations between these genomes. Whereas a single genus (Megadenia) remained diploid, the four remaining genera originated by allopolyploidy (Biscutella, Lunaria, Ricotia) or autopolyploidy (Heldreichia). The contentious relationships among the Biscutelleae genera, and between the tribe and other early diverged crucifer lineages, are best explained by close genomic relatedness among the recurrently hybridizing ancestral genomes. By using complementary cytogenomics and phylogenomics approaches, we demonstrate that the origin of a monophyletic plant clade can be more complex than a parsimonious assumption of a single WGD spurring postpolyploid cladogenesis. Instead, recurrent hybridization among the same and/or closely related parental genomes may phylogenetically interlink diploid and polyploid genomes despite the incidence of multiple independent WGDs. Our results provide new insights into evolution of early-diverging Brassicaceae lineages and elucidate challenges in resolving the contentious relationships within and between land plant lineages with pervasive hybridization and WGDs.     

The American Paddlefish Genome Provides Novel Insights into Chromosomal Evolution and Bone Mineralization in Early Vertebrates

Sturgeons and paddlefishes (Acipenseriformes) occupy the basal position of ray-finned fishes, although they have cartilaginous skeletons as in Chondrichthyes. This evolutionary status and their morphological specializations make them a research focus, but their complex genomes (polyploidy and the presence of microchromosomes) bring obstacles and challenges to molecular studies. Here, we generated the first high-quality genome assembly of the American paddlefish (Polyodon spathula) at a chromosome level. Comparative genomic analyses revealed a recent species-specific whole-genome duplication event, and extensive chromosomal changes, including head-to-head fusions of pairs of intact, large ancestral chromosomes within the paddlefish. We also provide an overview of the paddlefish SCPP (secretory calcium-binding phosphoprotein) repertoire that is responsible for tissue mineralization, demonstrating that the earliest flourishing of SCPP members occurred at least before the split between Acipenseriformes and teleosts. In summary, this genome assembly provides a genetic resource for understanding chromosomal evolution in polyploid nonteleost fishes and bone mineralization in early vertebrates.     

Evolution of RAG-1 in polyploid clawed frogs

Possible genetic fates of a gene duplicate are silencing, redundancy, subfunctionalization, or novel function. These different fates can be realized at the DNA, RNA, or protein level, and their genetic determinants are poorly understood. We explored molecular evolution of duplicated RAG-1 genes in African clawed frogs (Xenopus and Silurana) (1) to examine the fate of paralogs of this gene at the DNA level in terms of recombination, positive selection, and gene degeneration and in the absence of extensive recombination among alleles at different paralogs, (2) to test phylogenetic hypotheses about the origins of polyploid species. We found that recombination between different RAG-1 paralogs is infrequent, that degeneration of some paralogs has occurred via stop codons and frameshift mutations, and that this degeneration occurred in paralogs inherited from only one diploid progenitor species. Simulations and phylogenetic analyses of RAG-1 and mitochondrial DNA support one origin of extant tetraploids in Xenopus and at least one origin in Silurana, five allopolyploid origins of extant octoploids, and two allopolyploid origins of extant dodecaploids. In allopolyploid species, which inherit a complete genome from two different ancestors, genes inherited from the same ancestor have a longer period of coevolution than genes inherited from different ancestors. Because of this, gene ancestry could potentially influence gene fate: interacting paralogs derived from the same lower ploidy ancestor might have similar genetic destinies.     

Chromosome-level genome of Tibetan naked carp (Gymnocypris przewalskii) provides insights into Tibetan highland adaptation

Gymnocypris przewalskii, a cyprinid fish endemic to the Qinghai-Tibetan Plateau, has evolved unique morphological, physiological and genetic characteristics to adapt to the highland environment. Herein, we assembled a high-quality G. przewalskii tetraploid genome with a size of 2.03 Gb and scaffold N50 of 44.93 Mb, which was anchored onto 46 chromosomes. The comparative analysis suggested that gene families related to highland adaptation were significantly expanded in G. przewalskii. According to the G. przewalskii genome, we evaluated the phylogenetic relationship of 13 schizothoracine fishes, and inferred that the demographic history of G. przewalskii was strongly associated with geographic and eco-environmental alterations. We noticed that G. przewalskii experienced whole-genome duplication, and genes preserved post duplication were functionally associated with adaptation to high salinity and alkalinity. In conclusion, a chromosome-scale G. przewalskii genome provides an important genomic resource for teleost fish, and will particularly promote our understanding of the molecular evolution and speciation of fish in the highland environment.     

直翅目昆虫线粒体基因组研究进展

本文总结了本实验室对40余种直翅目昆虫的线粒体基因组序列的研究方法和主要结果.直翅目线粒体基因组研究中最重要的发现包括：（1）在直翅目昆虫线粒体基因组中发现了3种基因排列次序.蝗亚目除蜢总科外都具有DK排列.蜢总科的变色乌蜢为KD 排列,与蝗亚目其他总科不同,而与螽亚目昆虫的排序方式相同.已测出的螽亚目大多数昆虫的KD 排列顺序与典型节肢动物的完全相同,但在黄脸油葫芦Teleogryllus emma发生了tRNA^Glu,tRNA^Ser和tRNA^Asn的倒置;（2）在疑钩额螽Ruspolia dubia中发现了一种到目前为止具有最短控制区（70 bp）的线粒体基因组;（3）采用多种方法分析了昆虫A+T富集区存在的调控序列和二级结构特征,获得了昆虫A+T富集区保守序列的一致结构.采用Z曲线分析蝗虫的A+T富集区,表明也存在与原核生物复制起点类似的信号;（4）构建了30种蝗虫12S rRNA和16S rRNA的二级结构.在昆虫线粒体基因组非编码链中发现了一些类tRNA结构和tRNA异构体;（5）构建了基于线粒体基因组数据的直翅目昆虫主要亚科以上分类单元之间的系统发育关系.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

Evolution of follistatin in teleosts revealed through phylogenetic, genomic and expression analyses

Follistatin (Fst) inhibits transforming growth factor-β (TGF-B) proteins and is a known regulator of amniote myogenesis. Here, we used phylogenetic, genomic and experimental approaches to study its evolution in teleosts. Phylogenetic analyses suggested that one fst gene (fst1) is common to euteleosts, but a second gene (fst2) is conserved specifically within the Ostariophysi. Zebrafish fst1/2 respectively appear on chromosomes 5 and 10 in two genomic regions, each with conserved synteny to a single region in tetrapods. Interestingly, other teleosts have two corresponding chromosomal regions with a similar repertoire of paralogues. Phylogenetic reconstruction clustered these gene duplicates into two sister clades branching from tetrapod sequences. We suggest that an ancestral fst-containing chromosome was duplicated during the teleost whole genome duplication, but that fst2 was lost in lineages external to the Ostariophysi. We show that Fst1 of teleosts/mammals has evolved under strong purifying selection, but the N-terminal of Fst2 may have evolved under positive selection. Furthermore, the tissue-specific expression of zebrafish fst2 was restricted to fewer tissues compared to its paralogue and the single fst1 orthologue of Atlantic salmon (Salmo salar). Zebrafish fst1/2 may have subfunctionalized relative to non-duplicated vertebrate lineages, as several regions in the fst promoter of tetrapods were conserved with one paralogue, but not both. Finally, we examined the embryonic expression of fst1 in a teleost outside the Ostariophysi (Atlantic salmon). During segmentation, fst1 was expressed in the anterior somite compartment but was excluded from muscle progenitors that strongly expressed myogenic regulatory factors (MRFs). Later, fst1 was expressed in myogenic progenitors of the pectoral fin buds and also within the pax7 ⁺ cell layer external to the myotome. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

One Hundred and One Years of Chemotaxis. Pfeffer,Pheromones, and Fertilization

Chemotaxis of flagellated unicellular gametes was discovered by Wilhelm Pfeffer. In most known systems the lure is a metabolic waste product ritualized for sexual attraction. A special case are the female gametes of brown algae. They secrete polyunsaturated C-8 or C-11 hydrocarbons mostly with straight chain or unsubstituted tri-, penta- or heptacyclic structures. Their biogenetic origin stems from polyunsaturated fatty acids. The reactions leading to the different structures comprise lipoxygenase-type radical formations and dehydrogenating decarboxylations. — Most gametes produce a variety of such hydrocarbons; however, only one of them is the active principle. Biological activity is highly stereospecific; but frequently the bouquets contain lures for other species. Although cross-fertilization is excluded, they may be used to defend the habitat against overgrowth by other seaweeds. Besides chemoattractive action, some secretions also synchronize hatching of the spermatozoids. One of the compounds responsible is characterized by an epoxy substitution. — For most of the structures, stereospecific and versatile syntheses are available, which were used to investigate the binding site. This forms a template for the signal molecule with specifically situated polarizing groups that interact with the π-electrons of the double bonds. The message is transformed into changes of flagellar beat patterns which ensure meeting of the partners followed by fertilization.