Differences in evolution rates among eudicotyledon species observed by analysis of protein divergence

The wild grape genome sequence provides insights into the transition from dioecy to hermaphroditism during grape domestication

Background

A key step in domestication of the grapevine was the transition from separate sexes (dioecy) in wild Vitis vinifera ssp. sylvestris (V. sylvestris) to hermaphroditism in cultivated Vitis vinifera ssp. sativa (V. vinifera). It is known that V. sylvestris has an XY system and V. vinifera a modified Y haplotype (Yh) and that the sex locus is small, but it has not previously been precisely characterized.

Results

We generate a high-quality de novo reference genome for V. sylvestris, onto which we map whole-genome re-sequencing data of a cross to locate the sex locus. Assembly of the full X, Y, and Yh haplotypes of V. sylvestris and V. vinifera sex locus and examining their gene content and expression profiles during flower development in wild and cultivated accessions show that truncation and deletion of tapetum and pollen development genes on the X haplotype likely causes male sterility, while the upregulation of a Y allele of a cytokinin regulator (APRT3) may cause female sterility. The downregulation of this cytokinin regulator in the Yh haplotype may be sufficient to trigger reversal to hermaphroditism. Molecular dating of X and Y haplotypes is consistent with the sex locus being as old as the Vitis genus, but the mechanism by which recombination was suppressed remains undetermined.

Conclusions

We describe the genomic and evolutionary characterization of the sex locus of cultivated and wild grapevine, providing a coherent model of sex determination in the latter and for transition from dioecy to hermaphroditism during domestication.

Different divergence events for three pairs of PEBPs in Gossypium as implied by evolutionary analysis

INTRODUCTION

The phosphatidylethanolamine-binding protein (PEBP) gene family plays a crucial role in seed germination, reproductive transformation, and other important developmental processes in plants, but its distribution in Gossypium genomes or species, evolutionary properties, and the fates of multiple duplicated genes remain unclear.

OBJECTIVES

The primary objectives of this study were to elucidate the distribution and characteristics of PEBP genes in Gossypium, as well as the evolutionary pattern of duplication and deletion, and functional differentiation of PEBPs in plants.

METHODS

Using the PEBP protein sequences in Arabidopsis thaliana as queries, blast alignment was carried out for the identification of PEBP genes in four sequenced cotton species. Using the primers designed according to the PEBP genome sequences, PEBP genes were cloned from 15 representative genomes of Gossypium genus, and the gene structure, CDS sequence, protein sequence and properties were predicted and phylogenetic analysis was performed. Taking PEBP proteins of grape as reference, grouping of orthologous gene, analysis of phylogeny and divergence of PEBPs in nine species were conducted to reconstruct the evolutionary pattern of PEBP genes in plants.

RESULTS

We identified and cloned 160 PEBPs from 15 cotton species, and the phylogenetic analysis showed that the genes could be classified into the following three subfamilies: MFT-like, FT-like and TFL1-like. There were eight single orthologous group (OG) members in each diploid and 16 double OG members in each tetraploid. An analysis of the expression and selective pressure indicated that expression divergence and strong purification selection within the same OG presented in the PEBP gene family.

CONCLUSION

An evolutionary pattern of duplication and deletion of the PEBP family in the evolutionary history of Gossypium was suggested, and three pairs of genes resulted from different whole-genome duplication events.

Diverse genome organization following 13 independent mesopolyploid events in Brassicaceae contrasts with convergent patterns of gene retention

Hybridization and polyploidy followed by genome-wide diploidization had a significant impact on the diversification of land plants. The ancient At-α whole-genome duplication (WGD) preceded the diversification of crucifers (Brassicaceae). Some genera and tribes also experienced younger, mesopolyploid WGDs concealed by subsequent genome diploidization. Here we tested if multiple base chromosome numbers originated due to genome diploidization after independent mesopolyploid WGDs and how diploidization affected post-polyploid gene retention. Sixteen species representing 10 Brassicaceae tribes were analyzed by comparative chromosome painting and/or whole-transcriptome analysis of gene age distributions and phylogenetic analyses of gene duplications. Overall, we found evidence for at least 13 independent mesopolyploidies followed by different degrees of diploidization across the Brassicaceae. New mesotetraploid events were uncovered for the tribes Anastaticeae, Iberideae and Schizopetaleae, and mesohexaploid WGDs for Cochlearieae and Physarieae. In contrast, we found convergent patterns of gene retention and loss among these independent WGDs. Our combined analyses of genomic data for Brassicaceae indicate that extant chromosome number variation in many plant groups, and especially monophyletic taxa with multiple base chromosome numbers, can result from clade-specific genome duplications followed by diploidization. Our observation of parallel gene retention and loss across multiple independent WGDs provides one of the first multi-species tests of the predictability of patterns of post-polyploid genome evolution.

Evolutionary Analyses of GRAS Transcription Factors in Angiosperms

GRAS transcription factors (TFs) play critical roles in plant growth and development such as gibberellin and mycorrhizal signaling. Proteins belonging to this gene family contain a typical GRAS domain in the C-terminal sequence, whereas the N-terminal region is highly variable. Although, GRAS genes have been characterized in a number of plant species, their classification is still not completely resolved. Based on a panel of eight representative species of angiosperms, we identified 29 orthologous groups or orthogroups (OGs) for the GRAS gene family, suggesting that at least 29 ancestor genes were present in the angiosperm lineage before the "Amborella" evolutionary split. Interestingly, some taxonomic groups were missing members of one or more OGs. The gene number expansion usually observed in transcription factors was not observed in GRAS while the genome triplication ancestral to the eudicots (γ hexaploidization event) was detectable in a limited number of GRAS orthogroups. We also found conserved OG-specific motifs in the variable N-terminal region. Finally, we could regroup OGs in 17 subfamilies for which names were homogenized based on a literature review and described 5 new subfamilies (DLT, RAD1, RAM1, SCLA, and SCLB). This study establishes a consistent framework for the classification of GRAS members in angiosperm species, and thereby a tool to correctly establish the orthologous relationships of GRAS genes in most of the food crops in order to facilitate any subsequent functional analyses in the GRAS gene family. The multi-fasta file containing all the sequences used in our study could be used as database to perform diagnostic BLASTp to classify GRAS genes from other non-model species.

Regulatory divergence between parental alleles determines gene expression patterns in hybrids

Both hybridization and allopolyploidization generate novel phenotypes by conciliating divergent genomes and regulatory networks in the same cellular context. To understand the rewiring of gene expression in hybrids, the total expression of 21,025 genes and the allele-specific expression of over 11,000 genes were quantified in interspecific hybrids and their parental species, Coffea canephora and Coffea eugenioides using RNA-seq technology. Between parental species, cis- and trans-regulatory divergences affected around 32% and 35% of analyzed genes, respectively, with nearly 17% of them showing both. The relative importance of trans-regulatory divergences between both species could be related to their low genetic divergence and perennial habit. In hybrids, among divergently expressed genes between parental species and hybrids, 77% was expressed like one parent (expression level dominance), including 65% like C. eugenioides. Gene expression was shown to result from the expression of both alleles affected by intertwined parental trans-regulatory factors. A strong impact of C. eugenioides trans-regulatory factors on the upregulation of C. canephora alleles was revealed. The gene expression patterns appeared determined by complex combinations of cis- and trans-regulatory divergences. In particular, the observed biased expression level dominance seemed to be derived from the asymmetric effects of trans-regulatory parental factors on regulation of alleles. More generally, this study illustrates the effects of divergent trans-regulatory parental factors on the gene expression pattern in hybrids. The characteristics of the transcriptional response to hybridization appear to be determined by the compatibility of gene regulatory networks and therefore depend on genetic divergences between the parental species and their evolutionary history.

The coffee genome provides insight into the convergent evolution of caffeine biosynthesis

Coffee is a valuable beverage crop due to its characteristic flavor, aroma, and the stimulating effects of caffeine. We generated a high-quality draft genome of the species Coffea canephora, which displays a conserved chromosomal gene order among asterid angiosperms. Although it shows no sign of the whole-genome triplication identified in Solanaceae species such as tomato, the genome includes several species-specific gene family expansions, among them N-methyltransferases (NMTs) involved in caffeine production, defense-related genes, and alkaloid and flavonoid enzymes involved in secondary compound synthesis. Comparative analyses of caffeine NMTs demonstrate that these genes expanded through sequential tandem duplications independently of genes from cacao and tea, suggesting that caffeine in eudicots is of polyphyletic origin.

Shifts in the evolutionary rate and intensity of purifying selection between two Brassica genomes revealed by analyses of orthologous transposons and relics of a whole genome triplication

Recent sequencing of the Brassica rapa and Brassica oleracea genomes revealed extremely contrasting genomic features such as the abundance and distribution of transposable elements between the two genomes. However, whether and how these structural differentiations may have influenced the evolutionary rates of the two genomes since their split from a common ancestor are unknown. Here, we investigated and compared the rates of nucleotide substitution between two long terminal repeats (LTRs) of individual orthologous LTR-retrotransposons, the rates of synonymous and non-synonymous substitution among triplicated genes retained in both genomes from a shared whole genome triplication event, and the rates of genetic recombination estimated/deduced by the comparison of physical and genetic distances along chromosomes and ratios of solo LTRs to intact elements. Overall, LTR sequences and genic sequences showed more rapid nucleotide substitution in B. rapa than in B. oleracea. Synonymous substitution of triplicated genes retained from a shared whole genome triplication was detected at higher rates in B. rapa than in B. oleracea. Interestingly, non-synonymous substitution was observed at lower rates in the former than in the latter, indicating shifted densities of purifying selection between the two genomes. In addition to evolutionary asymmetry, orthologous genes differentially regulated and/or disrupted by transposable elements between the two genomes were also characterized. Our analyses suggest that local genomic and epigenomic features, such as recombination rates and chromatin dynamics reshaped by independent proliferation of transposable elements and elimination between the two genomes, are perhaps partially the causes and partially the outcomes of the observed inter-specific asymmetric evolution.