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.     

Phylogenetic perspectives on the origins of nodulation

Recent refinements to the phylogeny of rosid angiosperms support the conclusion that nodulation has evolved several times in the so-called N(2)-fixing clade (NFC), and provide dates for these origins. The hypothesized predisposition that enabled the evolution of nodulation occurred approximately 100 million years ago (MYA), was retained in the various lineages that radiated rapidly shortly thereafter, and was functional in its non-nodulation role for at least an additional 30 million years in each nodulating lineage. Legumes radiated rapidly shortly after their origin approximately 60 MYA, and nodulation most likely evolved several times during this radiation. The major lineages of papilionoid legumes diverged close to the time of origin of nodulation, accounting for the diversity of nodule biology in the group. Nodulation symbioses exemplify the concept of "deep homology," sharing various homologous components across nonhomologous origins of nodulation, largely due to recruitment from existing functions, notably the older arbuscular mycorrhizal symbiosis. Although polyploidy may have played a role in the origin of papilionoid legume nodules, it did not do so in other legumes, nor did the prerosid whole-genome triplication lead directly to the predisposition of nodulation.     

The dalbergioid legumes (Fabaceae): delimitation of a pantropical monophyletic clade

A monophyletic pantropical group of papilionoid legumes, here referred to as the "dalbergioid" legumes, is circumscribed to include all genera previously referred to the tribes Aeschynomeneae and Adesmieae, the subtribe Bryinae of the Desmodieae, and tribe Dalbergieae except Andira, Hymenolobium, Vatairea, and Vataireopsis. This previously undetected group was discovered with phylogenetic analysis of DNA sequences from the chloroplast trnK (including matK) and trnL introns, and the nuclear ribosomal 5.8S and flanking internal transcribed spacers 1 and 2. All dalbergioids belong to one of three well-supported subclades, the Adesmia, Dalbergia, and Pterocarpus clades. The dalbergioid clade and its three main subclades are cryptic in the sense that they are genetically distinct but poorly, if at all, distinguished by nonmolecular data. Traditionally important taxonomic characters, such as arborescent habit, free stamens, and lomented pods, do not provide support for the major clades identified by the molecular analysis. Short shoots, glandular-based trichomes, bilabiate calyces, and aeschynomenoid root nodules, in contrast, are better indicators of relationship at this hierarchical level. The discovery of the dalbergioid clade prompted a re-analysis of root nodule structure and the subsequent finding that the aeschynomenoid root nodule is synapomorphic for the dalbergioids.     

A phylogenetic strategy based on a legume-specific whole genome duplication yields symbiotic cytokinin type-A response regulators

Legumes host their Rhizobium spp. symbiont in novel root organs called nodules. Nodules originate from differentiated root cortical cells that dedifferentiate and subsequently form nodule primordia, a process controlled by cytokinin. A whole-genome duplication has occurred at the root of the legume Papilionoideae subfamily. We hypothesize that gene pairs originating from this duplication event and are conserved in distinct Papilionoideae lineages have evolved symbiotic functions. A phylogenetic strategy was applied to search for such gene pairs to identify novel regulators of nodulation, using the cytokinin phosphorelay pathway as a test case. In this way, two paralogous type-A cytokinin response regulators were identified that are involved in root nodule symbiosis. Response Regulator9 (MtRR9) and MtRR11 in medicago (Medicago truncatula) and an ortholog in lotus (Lotus japonicus) are rapidly induced upon Rhizobium spp. Nod factor signaling. Constitutive expression of MtRR9 results in arrested primordia that have emerged from cortical, endodermal, and pericycle cells. In legumes, lateral root primordia are not exclusively formed from pericycle cells but also require the involvement of the root cortical cell layer. Therefore, the MtRR9-induced foci of cell divisions show a strong resemblance to lateral root primordia, suggesting an ancestral function of MtRR9 in this process. Together, these findings provide a proof of principle for the applied phylogenetic strategy to identify genes with a symbiotic function in legumes.     

Crown group Oxyphotobacteria postdate the rise of oxygen

The rise of oxygen ca. 2.3 billion years ago (Ga) is the most distinct environmental transition in Earth history. This event was enabled by the evolution of oxygenic photosynthesis in the ancestors of Cyanobacteria. However, long‐standing questions concern the evolutionary timing of this metabolism, with conflicting answers spanning more than one billion years. Recently, knowledge of the Cyanobacteria phylum has expanded with the discovery of non‐photosynthetic members, including a closely related sister group termed Melainabacteria, with the known oxygenic phototrophs restricted to a clade recently designated Oxyphotobacteria. By integrating genomic data from the Melainabacteria, cross‐calibrated Bayesian relaxed molecular clock analyses show that crown group Oxyphotobacteria evolved ca. 2.0 billion years ago (Ga), well after the rise of atmospheric dioxygen. We further estimate the divergence between Oxyphotobacteria and Melainabacteria ca. 2.5–2.6 Ga, which—if oxygenic photosynthesis is an evolutionary synapomorphy of the Oxyphotobacteria—marks an upper limit for the origin of oxygenic photosynthesis. Together, these results are consistent with the hypothesis that oxygenic photosynthesis evolved relatively close in time to the rise of oxygen.     

Conservation of Arabidopsis flowering genes in model legumes

The model plants Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) have provided a wealth of information about genes and genetic pathways controlling the flowering process, but little is known about the corresponding pathways in legumes. The garden pea (Pisum sativum) has been used for several decades as a model system for physiological genetics of flowering, but the lack of molecular information about pea flowering genes has prevented direct comparison with other systems. To address this problem, we have searched expressed sequence tag and genome sequence databases to identify flowering-gene-related sequences from Medicago truncatula, soybean (Glycine max), and Lotus japonicus, and isolated corresponding sequences from pea by degenerate-primer polymerase chain reaction and library screening. We found that the majority of Arabidopsis flowering genes are represented in pea and in legume sequence databases, although several gene families, including the MADS-box, CONSTANS, and FLOWERING LOCUS T/TERMINAL FLOWER1 families, appear to have undergone differential expansion, and several important Arabidopsis genes, including FRIGIDA and members of the FLOWERING LOCUS C clade, are conspicuously absent. In several cases, pea and Medicago orthologs are shown to map to conserved map positions, emphasizing the closely syntenic relationship between these two species. These results demonstrate the potential benefit of parallel model systems for an understanding of flowering phenology in crop and model legume species.     

Computational identification and characterization of novel genes from legumes

The Fabaceae, the third largest family of plants and the source of many crops, has been the target of many genomic studies. Currently, only the grasses surpass the legumes for the number of publicly available expressed sequence tags (ESTs). The quantity of sequences from diverse plants enables the use of computational approaches to identify novel genes in specific taxa. We used BLAST algorithms to compare unigene sets from Medicago truncatula, Lotus japonicus, and soybean (Glycine max and Glycine soja) to nonlegume unigene sets, to GenBank's nonredundant and EST databases, and to the genomic sequences of rice (Oryza sativa) and Arabidopsis. As a working definition, putatively legume-specific genes had no sequence homology, below a specified threshold, to publicly available sequences of nonlegumes. Using this approach, 2,525 legume-specific EST contigs were identified, of which less than three percent had clear homology to previously characterized legume genes. As a first step toward predicting function, related sequences were clustered to build motifs that could be searched against protein databases. Three families of interest were more deeply characterized: F-box related proteins, Pro-rich proteins, and Cys cluster proteins (CCPs). Of particular interest were the >300 CCPs, primarily from nodules or seeds, with predicted similarity to defensins. Motif searching also identified several previously unknown CCP-like open reading frames in Arabidopsis. Evolutionary analyses of the genomic sequences of several CCPs in M. truncatula suggest that this family has evolved by local duplications and divergent selection.     

Placing paleopolyploidy in relation to taxon divergence: a phylogenetic analysis in legumes using 39 gene families

Young polyploid events are easily diagnosed by various methods, but older polyploid events become increasingly difficult to identify as chromosomal rearrangements, tandem gene or partial chromosome duplications, changes in substitution rates among duplicated genes, pseudogenization or locus loss, and interlocus interactions complicate the means of inferring past genetic events. Genomic data have provided valuable information about the polyploid history of numerous species, but on their own fail to show whether related species, each with a polyploid past, share a particular polyploid event. A phylogenetic approach provides a powerful method to determine this but many processes may mislead investigators. These processes can affect individual gene trees, but most likely will not affect all genes, and almost certainly will not affect all genes in the same way. Thus, a multigene approach, which combines the large-scale aspect of genomics with the resolution of phylogenetics, has the power to overcome these difficulties and allow us to infer genomic events further into the past than would otherwise be possible. Previous work using synonymous distances among gene pairs within species has shown evidence for large-scale duplications in the legumes Glycine max and Medicago truncatula. We present a case study using 39 gene families, each with three or four members in G. max and the putative orthologues in M. truncatula, rooted using Arabidopsis thaliana. We tested whether the gene duplications in these legumes occurred separately in each lineage after their divergence (Hypothesis 1), or whether they share a round of gene duplications (Hypothesis 2). Many more gene family topologies supported Hypothesis 2 over Hypothesis 1 (11 and 2, respectively), even after synonymous distance analysis revealed that some topologies were providing misleading results. Only ca. 33% of genes examined support either hypothesis, which strongly suggests that single gene family approaches may be insufficient when studying ancient events with nuclear DNA. Our results suggest that G. max and M. truncatula, along with approximately 7000 other legume species from the same clade, share an ancient round of gene duplications, either due to polyploidy or to some other process.     

A novel family of lectins evolutionarily related to class V chitinases: an example of neofunctionalization in legumes

A lectin has been identified in black locust (Robinia pseudoacacia) bark that shares approximately 50% sequence identity with plant class V chitinases but is essentially devoid of chitinase activity. Specificity studies indicated that the black locust chitinase-related agglutinin (RobpsCRA) preferentially binds to high-mannose N-glycans comprising the proximal pentasaccharide core structure. Closely related orthologs of RobpsCRA could be identified in the legumes Glycine max, Medicago truncatula, and Lotus japonicus but in no other plant species, suggesting that this novel lectin family most probably evolved in an ancient legume species or possibly an earlier ancestor. This identification of RobpsCRA not only illustrates neofunctionalization in plants, but also provides firm evidence that plants are capable of developing a sugar-binding domain from an existing structural scaffold with a different activity and accordingly sheds new light on the molecular evolution of plant lectins.     

The chicken or the egg? Plastome evolution and an independent loss of the inverted repeat in papilionoid legumes

The plastid genome (plastome), while surprisingly constant in gene order and content across most photosynthetic angiosperms, exhibits variability in several unrelated lineages. During the diversification history of the legume family Fabaceae, plastomes have undergone many rearrangements, including inversions, expansion, contraction and loss of the typical inverted repeat (IR), gene loss and repeat accumulation in both shared and independent events. While legume plastomes have been the subject of study for some time, most work has focused on agricultural species in the IR-lacking clade (IRLC) and the plant model Medicago truncatula. The subfamily Papilionoideae, which contains virtually all of the agricultural legume species, also comprises most of the plastome variation detected thus far in the family. In this study three non-papilioniods were included among 34 newly sequenced legume plastomes, along with 33 publicly available sequences, to assess plastome structural evolution in the subfamily. In an effort to examine plastome variation across the subfamily, approximately 20% of the sampling represents the IRLC with the remainder selected to represent the early-branching papilionoid clades. A number of IR-related and repeat-mediated changes were identified and examined in a phylogenetic context. Recombination between direct repeats associated with ycf2 resulted in intraindividual plastome heteroplasmy. Although loss of the IR has not been reported in legumes outside of the IRLC, one genistoid taxon was found to completely lack the typical plastome IR. The role of the IR and non-IR repeats in the progression of plastome change is discussed.     

Phylogeny and character evolution in Medicago (Leguminosae): Evidence from analyses of plastid trnK/matK and nuclear GA3ox1 sequences

? Premise of the study: The genus Medicago, with about 87 species, includes the model legume species M. truncatula, and a number of important forage species such as M. sativa (alfalfa), M. scutellata (snail medic), and M. lupulina (black medic). Relationships within the genus are not yet sufficiently resolved, contributing to difficulty in understanding the evolution of a number of distinguishing characteristics such as aneuploidy and polyploidy, life history, structure of cotyledons, and number of seeds per fruit. ? Methods: Phylogenetic relationships of 70-73 species of Medicago and its sister genus Trigonella (including Melilotus) were reconstructed from nucleotide sequences of the plastid trnK/matK region and the nuclear-encoded GA3ox1 gene (gibberellin 3-β-hydroxylase) using maximum parsimony and Bayesian inference methods. ? Key results: Our results support certain currently recognized taxonomic groups, e.g., sect. Medicago (with M. sativa) and sect. Buceras. However, other strongly supported clades-the "reduced subsection Leptospireae clade" that includes M. lupulina, the "polymorpha clade" that includes M. murex and M. polymorpha and the "subsection Pachyspireae clade" that includes M. truncatula-each of which includes species presently in different subsections of sect. Spirocarpos, contradict the current classification. ? Conclusions: These results support the hypothesis that some characters considered important in existing taxonomies, for example, single-seeded fruits that have arisen more than once in both Medicago and Trigonella, are indeed homoplastic. Others, such as the 2n = 14 chromosome number, have also arisen independently within the genus. In addition, we demonstrate support for the utility of GA3ox1 sequences for phylogenetic analysis among and within closely related genera of legumes.     

Depression of legume growth by liming

Summary Early growth or nodulation of certain species was depressed when a Hawaiian Oxisol was limed at rates above 6 tons/ha (pH 6). In 8 legumes, the depression later gave way to positive response. This was evident in plant weights of Desmodium intortum and Glycine wightii var. Cooper, and in pod weights of Phaseolus vulgaris. A transient depression was observed visually in Desmodium canum, Dolichos axillaris, Glycine wightii var. Tinaroo, and Trifolium subterraneum. In Stylosanthes guyanensis and S. fruticosa, the depression persisted throughout the experiment (6 months). Growth was not depressed in Arachis hypogea, Coronilla varia, Glycine max, Leucaena leucocephala, Medicago sativa, Trifolium repens, or Vigna sinensis. Journal Series No. 1953 of the Hawaii Agr. Exp. Station Journal Series No. 1953 of the Hawaii Agr. Exp. Station     

Mining EST databases to resolve evolutionary events in major crop species.

Using plant EST collections, we obtained 1392 potential gene duplicates across 8 plant species: Zea mays, Oryza sativa, Sorghum bicolor, Hordeum vulgare, Solanum tuberosum, Lycopersicon esculentum, Medicago truncatula, and Glycine max. We estimated the synonymous and nonsynonymous distances between each gene pair and identified two to three mixtures of normal distributions corresponding to one to three rounds of genome duplication in each species. Within the Poaceae, we found a conserved duplication event among all four species that occurred approximately 50-60 million years ago (Mya); an event that probably occurred before the major radiation of the grasses. In the Solanaceae, we found evidence for a conserved duplication event approximately 50-52 Mya. A duplication in soybean occurred approximately 44 Mya and a duplication in Medicago about 58 Mya. Comparing synonymous and nonsynonymous distances allowed us to determine that most duplicate gene pairs are under purifying, negative selection. We calculated Pearson's correlation coefficients to provide us with a measure of how gene expression patterns have changed between duplicate pairs, and compared this across evolutionary distances. This analysis showed that some duplicates seemed to retain expression patterns between pairs, whereas others showed uncorrelated expression.     

Legume agglutinins that bind to Rhizobium meliloti.

A protein found in seeds and roots of alfalfa (Medicago sativa) was implicated in the specificity of the infection process, based on its binding to the symbiont Rhizobium meliloti. We found an agglutinin with similar properties in seeds and roots of sweet clover (Melilotis alba). The sweet clover differed from alfalfa in nodulation by a mutant strain of R. meliloti, but the agglutinins were indistinguishable by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, Rhizobium agglutination, and cross-reactivity to antibodies. Similar agglutinins binding R. meliloti were found in seeds of legumes from different cross-inoculation groups, including soybean (Glycine max), cowpea (Vigna unguiculata), pea (Pisum sativum L), and mung bean (Vigna mungo). The agglutinins from these legumes were recognized by antibodies raised against the agglutinins of alfalfa and sweet clover. Seeds of corn (Zea mays) and tomato (Lycopersicon esculentum) contained a protein similar to the legume agglutinin, but it did not react with the antibodies. We conclude that the alfalfa agglutinin is representative of a common legume protein and that there is no evidence for its role in specificity or nodule initiation.     

Root nodulation studies inAeschynomene aspera

Summary Fifteen isolates of nodule bacteria were isolated from root and stem nodules ofAeschynomene aspera and they were characterized as Rhizobium by well known laboratory tests. All these isolates together with other efficient strains of known rhizobia belonging to different cross-inoculation groups were evaluated for their nodulation abilities onAeschynomene aspera, Cajanus cajan (pigeon pea),Cicer arietinum (chickpea),Pisum sativum (pea),Trifolium repens (clover),Medicago sativa (lucerne),Lens culinaris (lentil),Glycine max (soybean),Vigna sinensis (cowpea),Vigna radiata (mung bean),Vigna mungo (urd bean) andArachis hypogea (peanut). The results demonstrated that Rhizobium fromAeschynomene could form nodules only on its homologous host (Aeschynomene) but not on other legumes tested. Secondly, none of the rhizobia of other cross-inoculation groups could nodulateA. aspera.