Genome Duplication and Gene Loss Affect the Evolution of Heat Shock Transcription Factor Genes in Legumes

Whole-genome duplication events (polyploidy events) and gene loss events have played important roles in the evolution of legumes. Here we show that the vast majority of Hsf gene duplications resulted from whole genome duplication events rather than tandem duplication, and significant differences in gene retention exist between species. By searching for intraspecies gene colinearity (microsynteny) and dating the age distributions of duplicated genes, we found that genome duplications accounted for 42 of 46 Hsf-containing segments in Glycine max, while paired segments were rarely identified in Lotus japonicas, Medicago truncatula and Cajanus cajan. However, by comparing interspecies microsynteny, we determined that the great majority of Hsf-containing segments in Lotus japonicas, Medicago truncatula and Cajanus cajan show extensive conservation with the duplicated regions of Glycine max. These segments formed 17 groups of orthologous segments. These results suggest that these regions shared ancient genome duplication with Hsf genes in Glycine max, but more than half of the copies of these genes were lost. On the other hand, the Glycine max Hsf gene family retained approximately 75% and 84% of duplicated genes produced from the ancient genome duplication and recent Glycine-specific genome duplication, respectively. Continuous purifying selection has played a key role in the maintenance of Hsf genes in Glycine max. Expression analysis of the Hsf genes in Lotus japonicus revealed their putative involvement in multiple tissue-/developmental stages and responses to various abiotic stimuli. This study traces the evolution of Hsf genes in legume species and demonstrates that the rates of gene gain and loss are far from equilibrium in different species.     

Genome Duplication in Soybean (Glycine Subgenus Soja)

Restriction fragment length polymorphism mapping data from nine populations (Glycine max X G. soja and G. max X G. max) of the Glycine subgenus soja genome led to the identification of many duplicated segments of the genome. Linkage groups contained up to 33 markers that were duplicated on other linkage groups. The size of homoeologous regions ranged from 1.5 to 106.4 cM, with an average size of 45.3 cM. We observed segments in the soybean genome that were present in as many as six copies with an average of 2.55 duplications per segment. The presence of nested duplications suggests that at least one of the original genomes may have undergone an additional round of tetraploidization. Tetraploidization, along with large internal duplications, accounts for the highly duplicated nature of the genome of the subgenus. Quantitative trait loci for seed protein and oil showed correspondence across homoeologous regions, suggesting that the genes or gene families contributing to seed composition have retained similar functions throughout the evolution of the chromosomes.     

Polyploidy and genome evolution in plants

Genome doubling (polyploidy) has been and continues to be a pervasive force in plant evolution. Modern plant genomes harbor evidence of multiple rounds of past polyploidization events, often followed by massive silencing and elimination of duplicated genes. Recent studies have refined our inferences of the number and timing of polyploidy events and the impact of these events on genome structure. Many polyploids experience extensive and rapid genomic alterations, some arising with the onset of polyploidy. Survivorship of duplicated genes are differential across gene classes, with some duplicate genes more prone to retention than others. Recent theory is now supported by evidence showing that genes that are retained in duplicate typically diversify in function or undergo subfunctionalization. Polyploidy has extensive effects on gene expression, with gene silencing accompanying polyploid formation and continuing over evolutionary time.     

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. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s11103-005-4414-1     

A Comprehensive Analysis of the Cupin Gene Family in Soybean (Glycine max)

Cupin superfamily of proteins, including germin and germin-like proteins (GLPs) from higher plants, is known to play crucial roles in plant development and defense. To date, no systematic analysis has been conducted in soybean (Glycine max) incorporating genome organization, gene structure, expression compendium. In this study, 69 putative Cupin genes were identified from the whole-genome of soybean, which were non-randomly distributed on 17 of the 20 chromosomes. These Gmcupin proteins were phylogenetically clustered into ten distinct subgroups among which the gene structures were highly conserved. Eighteen pairs (52.2%) of duplicate paralogous genes were preferentially retained in duplicated regions of the soybean genome. The distributions of GmCupin genes implied that long segmental duplications contributed significantly to the expansion of the GmCupin gene family. According to the RNA-seq data analysis, most of the Gmcupins were differentially expressed in tissue-specific expression pattern and the expression of some duplicate genes were partially redundant while others showed functional diversity, suggesting the Gmcupins have been retained by substantial subfunctionalization during soybean evolutionary processes. Selective analysis based on single nucleotide polymorphisms (SNPs) in cultivated and wild soybeans revealed sixteen Gmcupins had selected site(s), with all SNPs in Gmcupin10.3 and Gmcupin07.2 genes were selected sites, which implied these genes may have undergone strong selection effects during soybean domestication. Taken together, our results contribute to the functional characterization of Gmcupin genes in soybean.     

Changes in twelve homoeologous genomic regions in soybean following three rounds of polyploidy

With the advent of high-throughput sequencing, the availability of genomic sequence for comparative genomics is increasing exponentially. Numerous completed plant genome sequences enable characterization of patterns of the retention and evolution of genes within gene families due to multiple polyploidy events, gene loss and fractionation, and differential evolutionary pressures over time and across different gene families. In this report, we trace the changes that have occurred in 12 surviving homoeologous genomic regions from three rounds of polyploidy that contributed to the current Glycine max genome: a genome triplication before the origin of the rosids (~130 to 240 million years ago), a genome duplication early in the legumes (~58 million years ago), and a duplication in the Glycine lineage (~13 million years ago). Patterns of gene retention following the genome triplication event generally support predictions of the Gene Balance Hypothesis. Finally, we find that genes in networks with a high level of connectivity are more strongly conserved than those with low connectivity and that the enrichment of these highly connected genes in the 12 highly conserved homoeologous segments may in part explain their retention over more than 100 million years and repeated polyploidy events.     

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.     

Sequencing wheat chromosome arm 7BS delimits the 7BS/4AL translocation and reveals homoeologous gene conservation

Complex Triticeae genomes pose a challenge to genome sequencing efforts due to their size and repetitive nature. Genome sequencing can reveal details of conservation and rearrangements between related genomes. We have applied Illumina second generation sequencing technology to sequence and assemble the low copy and unique regions of Triticum aestivum chromosome arm 7BS, followed by the construction of a syntenic build based on gene order in Brachypodium. We have delimited the position of a previously reported translocation between 7BS and 4AL with a resolution of one or a few genes and report approximately 13% genes from 7BS having been translocated to 4AL. An additional 13 genes are found on 7BS which appear to have originated from 4AL. The gene content of the 7DS and 7BS syntenic builds indicate a total of ~77,000 genes in wheat. Within wheat syntenic regions, 7BS and 7DS share 740 genes and a common gene conservation rate of ~39% of the genes from the corresponding regions in Brachypodium, as well as a common rate of colinearity with Brachypodium of ~60%. Comparison of wheat homoeologues revealed ~84% of genes previously identified in 7DS have a homoeologue on 7BS or 4AL. The conservation rates we have identified among wheat homoeologues and with Brachypodium provide a benchmark of homoeologous gene conservation levels for future comparative genomic analysis. The syntenic build of 7BS is publicly available at http://www.wheatgenome.info.     

Gene duplication within the Green Lineage: the case of TEL genes

Recent years have witnessed a breathtaking increase in the availability of genome sequence data, providing evidence of the highly duplicate nature of eukaryotic genomes. Plants are exceptional among eukaryotic organisms in that duplicate loci compose a large fraction of their genomes, partly because of the frequent occurrence of polyploidy (or whole-genome duplication) events. Tandem gene duplication and transposition have also contributed to the large number of duplicated genes in plant genomes. Evolutionary analyses allowed the dynamics of duplicate gene evolution to be studied and several models were proposed. It seems that, over time, many duplicated genes were lost and some of those that were retained gained new functions and/or expression patterns (neofunctionalization) or subdivided their functions and/or expression patterns between them (subfunctionalization). Recent studies have provided examples of genes that originated by duplication with successive diversification within plants. In this review, we focused on the TEL (TERMINAL EAR1-like) genes to illustrate such mechanisms. Emerged from the mei2 gene family, these TEL genes are likely to be land plant-specific. Phylogenetic analyses revealed one or two TEL copies per diploid genome. TEL gene degeneration and loss in several Angiosperm species such as in poplar and maize seem to have occurred. In Arabidopsis thaliana, whose genome experienced at least three polyploidy events followed by massive gene loss and genomic reorganization, two TEL genes were retained and two new shorter TEL-like (MCT) genes emerged. Molecular and expression analyses suggest for these genes sub- and neofunctionalization events, but confirmation will come from their functional characterization.     

Pericentromeric effects shape the patterns of divergence, retention, and expression of duplicated genes in the paleopolyploid soybean

The evolutionary forces that govern the divergence and retention of duplicated genes in polyploids are poorly understood. In this study, we first investigated the rates of nonsynonymous substitution (Ka) and the rates of synonymous substitution (Ks) for a nearly complete set of genes in the paleopolyploid soybean (Glycine max) by comparing the orthologs between soybean and its progenitor species Glycine soja and then compared the patterns of gene divergence and expression between pericentromeric regions and chromosomal arms in different gene categories. Our results reveal strong associations between duplication status and Ka and gene expression levels and overall low Ks and low levels of gene expression in pericentromeric regions. It is theorized that deleterious mutations can easily accumulate in recombination-suppressed regions, because of Hill-Robertson effects. Intriguingly, the genes in pericentromeric regions-the cold spots for meiotic recombination in soybean-showed significantly lower Ka and higher levels of expression than their homoeologs in chromosomal arms. This asymmetric evolution of two members of individual whole genome duplication (WGD)-derived gene pairs, echoing the biased accumulation of singletons in pericentromeric regions, suggests that distinct genomic features between the two distinct chromatin types are important determinants shaping the patterns of divergence and retention of WGD-derived genes.     

Allopolyploid origin of Cardamine asarifolia (Brassicaceae): incongruence between plastid and nuclear ribosomal DNA sequences solved by a single-copy nuclear gene

Interspecific hybridization and polyploidization have played central roles in plant diversification. However, technical difficulties in the analyses of low-copy genes have limited the study of the origins of hybrid and polyploid plants. Here, we present a phylogenetic analysis of the hexaploid Cardamine asarifolia, distributed in the southern European Alps and northern Apennines. Our study included all relevant taxa of the genus found in Europe. A marked discrepancy was revealed between the trnL-trnF region of cpDNA and internal transcribed spacer (nrDNA ITS) sequences. To solve the incongruence, we sequenced a single-copy nuclear CHS gene (chalcone synthase) using a novel method to design homoeologue-specific PCR primers to bypass artefacts caused by artificial recombination of homoeologues during PCR and/or cloning. Three homoeologues were isolated from C. asarifolia, providing evidence for its allopolyploid origin. One homoeologue, showing the same phylogenetic position as the ITS sequences, most likely originated from an extinct parent. Furthermore, we documented recurrent polytopic hybridizations between C. asarifolia and diploid C. amara. The allohexaploidization and the following hybridization with a diploid species exemplify the ongoing dynamic processes of speciation in the genus Cardamine.     

Functional conservation and divergence of GmCHLI genes in polyploid soybean

Polyploidy is prevalent in nature. As the fate of duplicated genes becomes more complicated when the encoded proteins function as oligomers, functional investigations into duplicated oligomer‐encoding genes in polyploid genomes will facilitate our understanding of how traits are expressed. In this study, we identified GmCHLI1, a gene encoding the I subunit of magnesium (Mg)‐chelatase, which functions in hexamers as responsible for the semi‐dominant etiolation phenotype in soybean. Four GmCHLI copies derived from two polyploidy events were identified in the soybean genome. Further investigation with regard to expression patterns indicated that these four copies have diverged into two pairs; mutation in the other copy of the pair that includes GmCHLI1 also resulted in a chlorophyll‐deficient phenotype. Protein interaction assays showed that these four GmCHLIs can interact with each other, but stronger interactions were found with mutated subunits. The results indicate that, in polyploidy, deficiency in each copy of duplicated oligomer‐encoding genes could result in a mutant phenotype due to hetero‐oligomer formation, which is different from the model of allelic dosage or functional redundancy. In addition, we interestingly found an increase in isoflavonoids in the heterozygous etiolated plants, which might be useful for improving soybean seed quality.     

Genome organization in dicots. II. Arabidopsis as a ’bridging species’ to resolve genome evolution events among legumes

Analysis of molecular linkage groups within the soybean (Glycine max L. Merr.) genome reveals many homologous regions, reflecting the ancient polyploidy of soybean. The fragmented arrangement of the duplicated regions suggests that extensive rearrangements, as well as additional duplications, have occurred since the initial polyploidization event. In this study we used comparisons between homoeologous regions in soybean, and the homologous regions in the related diploids Phaseolus vulgaris and Vigna radiata, to elucidate the evolutionary history of the three legume genomes. Our results show that there is not only conservation of large regions of the genomes but that these conserved linkage blocks are also represented twice in the soybean genome. To gain a better understanding of the process of genome evolution in dicots, molecular comparisons have been extended to another well-studied species, Arabidopsis thaliana. Interestingly, the conserved regions we identified in the legume species are also relatively conserved in Arabidopsis. Our results suggest that there is conservation of blocks of DNA between species as distantly related as legumes and brassicas, representing 90 million years of divergence. We also present evidence for an additional, presumably earlier, genome duplication in soybean. These duplicated regions were only recognized by using Arabidopsis as a ’bridging’ species in the genome comparisons. Received: 10 October 2000 / Accepted: 13 January 2001