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.     

Differential accumulation of retroelements and diversification of NB-LRR disease resistance genes in duplicated regions following polyploidy in the ancestor of soybean

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.     

Functional specialization of duplicated AP3‐like genes in Medicago truncatula

The B–class of MADS box genes has been studied in a wide range of plant species, but has remained largely uncharacterized in legumes. Here we investigate the evolutionary fate of the duplicated AP3‐like genes of a legume species. To obtain insight into the extent to which B‐class MADS box gene functions are conserved or have diversified in legumes, we isolated and characterized the two members of the AP3 lineage in Medicago truncatula: MtNMH7 and MtTM6 (euAP3 and paleoAP3 genes, respectively). A non‐overlapping and complementary expression pattern of both genes was observed in petals and stamens. MtTM6 was expressed predominantly in the outer cell layers of both floral organs, and MtNMH7 in the inner cell layers of petals and stamens. Functional analyses by reverse genetics approaches (RNAi and Tnt1 mutagenesis) showed that the contribution of MtNMH7 to petal identity is more important than that of MtTM6, whereas MtTM6 plays a more important role in stamen identity than its paralog MtNMH7. Our results suggest that the M. truncatula AP3‐like genes have undergone a functional specialization process associated with complete partitioning of gene expression patterns of the ancestral gene lineage. We provide information regarding the similarities and differences in petal and stamen development among core eudicots.     

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     

Divergence in expression between duplicated genes in Arabidopsis

New genes may arise through tandem duplication, dispersed small-scale duplication, and polyploidy, and patterns of divergence between duplicated genes may vary among these classes. We have examined patterns of gene expression and coding sequence divergence between duplicated genes in Arabidopsis thaliana. Due to the simultaneous origin of polyploidy-derived gene pairs, we can compare covariation in the rates of expression divergence and sequence divergence within this group. Among tandem and dispersed duplicates, much of the divergence in expression profile appears to occur at or shortly after duplication. Contrary to findings from other eukaryotic systems, there is little relationship between expression divergence and synonymous substitutions, whereas there is a strong positive relationship between expression divergence and nonsynonymous substitutions. Because this pattern is pronounced among the polyploidy-derived pairs, we infer that the strength of purifying selection acting on protein sequence and expression pattern is correlated. The polyploidy-derived pairs are somewhat atypical in that they have broader expression patterns and are expressed at higher levels, suggesting differences among polyploidy- and nonpolyploidy-derived duplicates in the types of genes that revert to single copy. Finally, within many of the duplicated pairs, 1 gene is expressed at a higher level across all assayed conditions, which suggests that the subfunctionalization model for duplicate gene preservation provides, at best, only a partial explanation for the patterns of expression divergence between duplicated genes.     

Structure and expression analysis of rice paleo duplications

Having a well-known history of genome duplication, rice is a good model for studying structural and functional evolution of paleo duplications. Improved sequence alignment criteria were used to characterize 10 major chromosome-to-chromosome duplication relationships associated with 1440 paralogous pairs, covering 47.8% of the rice genome, with 12.6% of genes that are conserved within sister blocks. Using a micro-array experiment, a genome-wide expression map has been produced, in which 2382 genes show significant differences of expression in root, leaf and grain. By integrating both structural (1440 paralogous pairs) and functional information (2382 differentially expressed genes), we identified 115 paralogous gene pairs for which at least one copy is differentially expressed in one of the three tissues. A vast majority of the 115 paralogous gene pairs have been neofunctionalized or subfunctionalized as 88%, 89% and 96% of duplicates, respectively, expressed in grain, leaf and root show distinct expression patterns. On the basis of a Gene Ontology analysis, we have identified and characterized the gene families that have been structurally and functionally preferentially retained in the duplication showing that the vast majority (>85%) of duplicated have been either lost or have been subfunctionalized or neofunctionalized during 50–70 million years of evolution.     

Gene duplication and divergence produce divergent MHC genotypes without disassortative mating

Genes of the major histocompatibility complex (MHC) exhibit heterozygote advantage in immune defence, which in turn can select for MHC‐disassortative mate choice. However, many species lack this expected pattern of MHC‐disassortative mating. A possible explanation lies in evolutionary processes following gene duplication: if two duplicated MHC genes become functionally diverged from each other, offspring will inherit diverse multilocus genotypes even under random mating. We used locus‐specific primers for high‐throughput sequencing of two expressed MHC Class II B genes in Leach's storm‐petrels, Oceanodroma leucorhoa, and found that exon 2 alleles fall into two gene‐specific monophyletic clades. We tested for disassortative vs. random mating at these two functionally diverged Class II B genes, using multiple metrics and different subsets of exon 2 sequence data. With good statistical power, we consistently found random assortment of mates at MHC. Despite random mating, birds had MHC genotypes with functionally diverged alleles, averaging 13 amino acid differences in pairwise comparisons of exon 2 alleles within individuals. To test whether this high MHC diversity in individuals is driven by evolutionary divergence of the two duplicated genes, we built a phylogenetic permutation model. The model showed that genotypic diversity was strongly impacted by sequence divergence between the most common allele of each gene, with a smaller additional impact of monophyly of the two genes. Divergence of allele sequences between genes may have reduced the benefits of actively seeking MHC‐dissimilar mates, in which case the evolutionary history of duplicated genes is shaping the adaptive landscape of sexual selection.     

Expression partitioning between genes duplicated by polyploidy under abiotic stress and during organ development

Allopolyploidy has been a prominent mode of speciation and a recurrent process during plant evolution and has contributed greatly to the large number of duplicated genes in plant genomes [1-4]. Polyploidy often leads to changes in genome organization and gene expression [5-9]. The expression of genes that are duplicated by polyploidy (termed homeologs) can be partitioned between the duplicates so that one copy is expressed and functions only in some organs and the other copy is expressed only in other organs, indicative of subfunctionalization [10]. To determine how homeologous-gene expression patterns change during organ development and in response to abiotic stress conditions, we have examined expression of the alcohol dehydrogenase gene AdhA in allopolyploid cotton (Gossypium hirsutum). Expression ratios of the two homeologs vary considerably during the development of organs from seedlings and fruits. Abiotic stress treatments, including cold, dark, and water submersion, altered homeologous-gene expression. Most notably, only one copy is expressed in hypocotyls during a water-submersion treatment, and only the other copy is expressed during cold stress. These results imply that subfunctionalization of genes duplicated by polyploidy has occurred in response to abiotic stress conditions. Partitioning of duplicate gene expression in response to environmental stress may lead to duplicate gene retention during subsequent evolution.