Phylogenomics of mulberries (Morus, Moraceae) inferred from plastomes and single copy nuclear genes.

Mulberries (genus Morus), belonging to the order Rosales, family Moraceae, are important woody plants due to their economic values in sericulture, as well as for nutritional benefits and medicinal values. However, the taxonomy and phylogeny of Morus, especially for the Asian species, remains challenging due to its wide geographical distribution, morphological plasticity, and interspecific hybridization. To better understand the evolutionary history of Morus, we combined plastomes and a large-scale nuclear gene analyses to investigate their phylogenetic relationships. We assembled the plastomes and screened 211 single-copy nuclear genes from 13 Morus species and related taxa. The plastomes of Morus species were relatively conserved in terms of genome size, gene content, synteny, IR boundary and codon usage. Using nuclear data, our results elucidated identical topologies based on coalescent and concatenation methods. The genus Morus was supported as monophyletic, with M. notabilis as the first diverging lineage and the two North American Morus species, M. microphylla and M. rubra, as sister to the other Asian species. In the Asian Morus species, interspecific relationships were completely resolved. However, cyto-nuclear discordances and gene tree-species tree conflicts were detected in the phylogenies of Morus, with multiple evidences supporting hybridization/introgression as the main cause of discordances between nuclear and plastid phylogenies, while gene tree-species tree conflicts were mainly caused by ILS.

Evolutionary gain and loss of a plant pattern-recognition receptor for HAMP recognition.

As a first step in innate immunity, pattern recognition receptors (PRRs) recognize the distinct pathogen and herbivore-associated molecular patterns and mediate activation of immune responses, but specific steps in the evolution of new PRR sensing functions are not well understood. We employed comparative genomic and functional analyses to define evolutionary events leading to the sensing of the herbivore-associated peptide inceptin (In11) by the PRR inceptin receptor (INR) in legume plant species. Existing and de novo genome assemblies revealed that the presence of a functional INR gene corresponded with ability to respond to In11 across ~53 million years (my) of evolution. In11 recognition is unique to the clade of Phaseoloid legumes, and only a single clade of INR homologs from Phaseoloids was functional in a heterologous model. The syntenic loci of several non-Phaseoloid outgroup species nonetheless contain non-functional INR-like homologs, suggesting that an ancestral gene insertion event and diversification preceded the evolution of a specific INR receptor function ~28 my ago. Chimeric and ancestrally reconstructed receptors indicated that 16 amino acid differences in the C1 leucine-rich repeat domain and C2 intervening motif mediate gain of In11 recognition. Thus, high PRR diversity was likely followed by a small number of mutations to expand innate immune recognition to a novel peptide elicitor. Analysis of INR evolution provides a model for functional diversification of other germline-encoded PRRs.

The health status of a plant depends on the immune system it inherits from its parents. Plants have many receptor proteins that can recognize distinct molecules from insects and microbes, and trigger an immune response. Inheriting the right set of receptors allows plants to detect certain threats and to cope with diseases and pests. Soybeans, chickpeas and other closely-related crop plants belong to a family of plants known as the legumes. Previous studies have found that, unlike other plants, some legumes are able to respond to oral secretions from caterpillars. These plants have a receptor known as INR that binds to a molecule called inceptin in the secretions. However, it remained unclear how or when INR evolved. To address this gap, Snoeck et al. tested immune responses to inceptin in the leaves of 22 species of legume. The experiments revealed that only members of a subgroup of legumes called the Phaseoloids were able to recognize the molecule. Analyzing the genomes of several legume species revealed that the gene encoding INR first emerged around 28 million years ago. Among the descendants of the legumes that first evolved this receptor, only the crop plant soybean and a few other species were unable to respond to inceptin. The genomic data indicated that these species had in fact lost the gene encoding INR over evolutionary time. Snoeck et al. then combined data from genes encoding modern-day receptors to reconstruct the sequence of building blocks that make up the 28-million-year-old version of INR. This ancestral receptor was able to respond to inceptin in the caterpillar secretion, whereas an older version of the protein, which had a slightly different set of building blocks, could not. This suggests that INR evolved the ability to respond to inceptin as a result of small mutations in the gene encoding a more ancient receptor. The work of Snoeck et al. reveals how the Phaseoloids evolved to respond to caterpillars, and how this ability has been lost in soybeans and other members of the subgroup. In the future, these findings may aid plant breeding or genetic engineering approaches for enhancing soybeans and other crops resistance to caterpillar pests.

Nuclear phylotranscriptomics and phylogenomics support numerous polyploidization events and hypotheses for the evolution of rhizobial nitrogen-fixing symbiosis in Fabaceae.

Fabaceae are the third largest angiosperm family, with 765 genera and ∼19 500 species. They are important both economically and ecologically, and global Fabaceae crops are intensively studied in part for their nitrogen-fixing ability. However, resolution of the intrasubfamilial Fabaceae phylogeny and divergence times has remained elusive, precluding a reconstruction of the evolutionary history of symbiotic nitrogen fixation in Fabaceae. Here, we report a highly resolved phylogeny using >1500 nuclear genes from newly sequenced transcriptomes and genomes of 391 species, along with other datasets, for a total of 463 legumes spanning all 6 subfamilies and 333 of 765 genera. The subfamilies are maximally supported as monophyletic. The clade comprising subfamilies Cercidoideae and Detarioideae is sister to the remaining legumes, and Duparquetioideae and Dialioideae are successive sisters to the clade of Papilionoideae and Caesalpinioideae. Molecular clock estimation revealed an early radiation of subfamilies near the K/Pg boundary, marked by mass extinction, and subsequent divergence of most tribe-level clades within ∼15 million years. Phylogenomic analyses of thousands of gene families support 28 proposed putative whole-genome duplication/whole-genome triplication events across Fabaceae, including those at the ancestors of Fabaceae and five of the subfamilies, and further analyses supported the Fabaceae ancestral polyploidy. The evolution of rhizobial nitrogen-fixing nodulation in Fabaceae was probed by ancestral character reconstruction and phylogenetic analyses of related gene families and the results support the hypotheses of one or two switch(es) to rhizobial nodulation followed by multiple losses. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Fabaceae.

Tandem gene duplication and recombination at the AT3 locus in the Solanaceae, a gene essential for capsaicinoid biosynthesis in Capsicum.

Capsaicinoids are compounds synthesized exclusively in the genus Capsicum and are responsible for the burning sensation experienced when consuming hot pepper fruits. To date, only one gene, AT3, a member of the BAHD family of acyltransferases, is currently known to have a measurable quantitative effect on capsaicinoid biosynthesis. Multiple AT3 paralogs exist in the Capsicum genome, but their evolutionary relationships have not been characterized well. Recessive alleles at this locus result in absence of capsaicinoids in pepper fruit. To explore the evolution of AT3 in Capsicum and the Solanaceae, we sequenced this gene from diverse Capsicum genotypes and species, along with a number of representative solanaceous taxa. Our results revealed that the coding region of AT3 is highly conserved throughout the family. Further, we uncovered a tandem duplication that predates the diversification of the Solanaceae taxa sampled in this study. This pair of tandem duplications were designated AT3-1 and AT3-2. Sequence alignments showed that the AT3-2 locus, a pseudogene, retains regions of amino acid conservation relative to AT3-1. Gene tree estimation demonstrated that AT3-1 and AT3-2 form well supported, distinct clades. In C. rhomboideum, a non-pungent basal Capsicum species, we describe a recombination event between AT3-1 and AT3-2 that modified the putative active site of AT3-1, also resulting in a frame-shift mutation in the second exon. Our data suggest that duplication of the original AT3 representative, in combination with divergence and pseudogene degeneration, may account for the patterns of sequence divergence and punctuated amino acid conservation observed in this study. Further, an early rearrangement in C. rhomboidium could account for the absence of pungency in this Capsicum species.

De novo transcriptome assembly of Pueraria montana var. lobata and Neustanthus phaseoloides for the development of eSSR and SNP markers: narrowing the US origin(s) of the invasive kudzu.

BACKGROUND: Kudzu, Pueraria montana var. lobata, is a woody vine native to Southeast Asia that has been introduced globally for cattle forage and erosion control. The vine is highly invasive in its introduced areas, including the southeastern US. Modern molecular marker resources are limited for the species, despite its importance. Transcriptomes for P. montana var. lobata and a second phaseoloid legume taxon previously ascribed to genus Pueraria, Neustanthus phaseoloides, were generated and mined for microsatellites and single nucleotide polymorphisms. RESULTS: Roche 454 sequencing of P. montana var. lobata and N. phaseoloides transcriptomes produced read numbers ranging from ~ 280,000 to ~ 420,000. Trinity assemblies produced an average of 17,491 contigs with mean lengths ranging from 639 bp to 994 bp. Transcriptome completeness, according to BUSCO, ranged between 64 and 77%. After vetting for primer design, there were 1646 expressed simple sequence repeats (eSSRs) identified in P. montana var. lobata and 1459 in N. phaseoloides. From these eSSRs, 17 identical primer pairs, representing inter-generic phaseoloid eSSRs, were created. Additionally, 13 primer pairs specific to P. montana var. lobata were also created. From these 30 primer pairs, a final set of seven primer pairs were used on 68 individuals of P. montana var. lobata for characterization across the US, China, and Japan. The populations exhibited from 20 to 43 alleles across the seven loci. We also conducted pairwise tests for high-confidence SNP discovery from the kudzu transcriptomes we sequenced and two previously sequenced P. montana var. lobata transcriptomes. Pairwise comparisons between P. montana var. lobata ranged from 358 to 24,475 SNPs, while comparisons between P. montana var. lobata and N. phaseoloides ranged from 5185 to 30,143 SNPs. CONCLUSIONS: The discovered molecular markers for kudzu provide a starting point for comparative genetic studies within phaseoloid legumes. This study both adds to the current genetic resources and presents the first available genomic resources for the invasive kudzu vine. Additionally, this study is the first to provide molecular evidence to support the hypothesis of Japan as a source of US kudzu and begins to narrow the origin of US kudzu to the central Japanese island of Honshu.

Targeting legume loci: A comparison of three methods for target enrichment bait design in Leguminosae phylogenomics.

PREMISE OF THE STUDY: The development of pipelines for locus discovery has spurred the use of target enrichment for plant phylogenomics. However, few studies have compared pipelines from locus discovery and bait design, through validation, to tree inference. We compared three methods within Leguminosae (Fabaceae) and present a workflow for future efforts. METHODS: Using 30 transcriptomes, we compared Hyb-Seq, MarkerMiner, and the Yang and Smith (Y&S) pipelines for locus discovery, validated 7501 baits targeting 507 loci across 25 genera via Illumina sequencing, and inferred gene and species trees via concatenation- and coalescent-based methods. RESULTS: Hyb-Seq discovered loci with the longest mean length. MarkerMiner discovered the most conserved loci with the least flagged as paralogous. Y&S offered the most parsimony-informative sites and putative orthologs. Target recovery averaged 93% across taxa. We optimized our targeted locus set based on a workflow designed to minimize paralog/ortholog conflation and thus present 423 loci for legume phylogenomics. CONCLUSIONS: Methods differed across criteria important for phylogenetic marker development. We recommend Hyb-Seq as a method that may be useful for most phylogenomic projects. Our targeted locus set is a resource for future, community-driven efforts to reconstruct the legume tree of life.

A Review on Current Status and Future Prospects of Winged Bean (Psophocarpus tetragonolobus) in Tropical Agriculture.

Winged bean, Psophocarpus tetragonolobus (L.) DC., is analogous to soybean in yield and nutritional quality, proving a valuable alternative to soybean in tropical regions of the world. The presence of anti-nutritional factors and high costs associated with indeterminate plant habit have been major concerns in this crop. But occurrence of good genetic variability in germplasm collections offers precious resources for winged bean breeding. However, lack of germplasm characterization is hindering such efforts. From a genomic standpoint, winged bean has been little studied despite rapid advancement in legume genomics in the last decade. Exploiting modern genomics/breeding approaches for genetic resource characterization and the breeding of early maturing, high yielding, determinate varieties which are disease resistant and free of anti-nutritional factors along with developing consumer friendly value-added products of local significance are great challenges and opportunities in the future that would boost cultivation of winged bean in the tropics. We review past efforts and future prospects towards winged bean improvement.

Transcriptomic resources for the medicinal legume Mucuna pruriens: de novo transcriptome assembly, annotation, identification and validation of EST-SSR markers.

BACKGROUND: The medicinal legume Mucuna pruriens (L.) DC. has attracted attention worldwide as a source of the anti-Parkinson's drug L-Dopa. It is also a popular green manure cover crop that offers many agronomic benefits including high protein content, nitrogen fixation and soil nutrients. The plant currently lacks genomic resources and there is limited knowledge on gene expression, metabolic pathways, and genetics of secondary metabolite production. Here, we present transcriptomic resources for M. pruriens, including a de novo transcriptome assembly and annotation, as well as differential transcript expression analyses between root, leaf, and pod tissues. We also develop microsatellite markers and analyze genetic diversity and population structure within a set of Indian germplasm accessions. RESULTS: One-hundred ninety-one million two hundred thirty-three thousand two hundred forty-two bp cleaned reads were assembled into 67,561 transcripts with mean length of 626 bp and N50 of 987 bp. Assembled sequences were annotated using BLASTX against public databases with over 80% of transcripts annotated. We identified 7,493 simple sequence repeat (SSR) motifs, including 787 polymorphic repeats between the parents of a mapping population. 134 SSRs from expressed sequenced tags (ESTs) were screened against 23 M. pruriens accessions from India, with 52 EST-SSRs retained after quality control. Population structure analysis using a Bayesian framework implemented in fastSTRUCTURE showed nearly similar groupings as with distance-based (neighbor-joining) and principal component analyses, with most of the accessions clustering per geographical origins. Pair-wise comparison of transcript expression in leaves, roots and pods identified 4,387 differentially expressed transcripts with the highest number occurring between roots and leaves. Differentially expressed transcripts were enriched with transcription factors and transcripts annotated as belonging to secondary metabolite pathways. CONCLUSIONS: The M. pruriens transcriptomic resources generated in this study provide foundational resources for gene discovery and development of molecular markers. Polymorphic SSRs identified can be used for genetic diversity, marker-trait analyses, and development of functional markers for crop improvement. The results of differential expression studies can be used to investigate genes involved in L-Dopa synthesis and other key metabolic pathways in M. pruriens.

Parsing polyphyletic Pueraria: Delimiting distinct evolutionary lineages through phylogeny.

Several taxonomic and phylogenetic studies have hypothesized polyphyly within Pueraria DC., a genus comprising 19 species (24 with varieties) including the highly invasive Pueraria montana var. lobata (Kudzu) introduced to the U.S.A. about 150years ago. Previous efforts to investigate monophyly of the genus have been hampered by limited taxon sampling or a lack of comprehensive evolutionary context that would enable definitive taxonomic associations. This work presents a comprehensive phylogenetic investigation of Pueraria within the context of tribe Phaseoleae (Leguminosae). Polyphyly was found to be more extensive than previously thought, with five distinct lineages spread across the tribe and spanning over 25mya of divergence strongly supported by two chloroplast and one nuclear marker, AS2, presented here as a phylogenetic marker for the first time. Our phylogenies support taxonomic revisions to rectify polyphyly within Pueraria, including the resurrection of Neustanthus, moving one species to Teyleria, and the creation of two new genera, Haymondia and Toxicopueraria (taxonomic revisions published elsewhere).

Transcriptome sequencing and marker development in winged bean (Psophocarpus tetragonolobus; Leguminosae).

Winged bean, Psophocarpus tetragonolobus (L.) DC., is similar to soybean in yield and nutritional value but more viable in tropical conditions. Here, we strengthen genetic resources for this orphan crop by producing a de novo transcriptome assembly and annotation of two Sri Lankan accessions (denoted herein as CPP34 [PI 491423] and CPP37 [PI 639033]), developing simple sequence repeat (SSR) markers, and identifying single nucleotide polymorphisms (SNPs) between geographically separated genotypes. A combined assembly based on 804,757 reads from two accessions produced 16,115 contigs with an N50 of 889 bp, over 90% of which has significant sequence similarity to other legumes. Combining contigs with singletons produced 97,241 transcripts. We identified 12,956 SSRs, including 2,594 repeats for which primers were designed and 5,190 high-confidence SNPs between Sri Lankan and Nigerian genotypes. The transcriptomic data sets generated here provide new resources for gene discovery and marker development in this orphan crop, and will be vital for future plant breeding efforts. We also analyzed the soybean trypsin inhibitor (STI) gene family, important plant defense genes, in the context of related legumes and found evidence for radiation of the Kunitz trypsin inhibitor (KTI) gene family within winged bean.

A multilocus phylogenetic analysis reveals the monophyly of a recircumscribed papilionoid legume tribe Diocleae with well-supported generic relationships.

Deciphering the phylogenetic relationships within the species-rich Millettioid clade has persisted as one of the major challenges in the systematics and evolutionary history of papilionoid legumes (Leguminosae, Papilionoideae). Historically, the predominantly neotropical lianas of subtribe Diocleinae in the Millettioid legumes have been taxonomically tangled together with the largely heterogeneous tribe Phaseoleae. This work presents a comprehensive molecular phylogenetic analysis based on nuclear and chloroplast markers and includes all genera ever referred to Diocleae except for the monospecific Philippine Luzonia, resolving several key generic relationships within the Millettioid legumes. The first of two separate analyses includes 310 matK accessions and strongly supports the reestablishment of tribe Diocleae as a branch of the Millettioid clade. This work sheds greater light on the higher-level phylogeny of Diocleae and allows the recognition of three major lineages: the Canavalia, Dioclea, and Galactia clades. The second set of phylogenetic analyses utilized nuclear (ITS/5.8S and ETS) and plastid (matK and trnT-Y) DNA sequences to reveal (i) the monophyly of Canavalia and Cleobulia; (ii) the monophyly of Bionia with the exclusion of Bionia bella; (iii) the paraphyly of Dioclea with respect to Cleobulia, Cymbosema, and Macropsychanthus; (iv) the paraphyly of Cratylia with respect to the broadly polyphyletic Camptosema; and (v) the polyphyly of Galactia with species scattered widely across the tree.

Species delimitation and recognition in the Pediomelummegalanthum complex (Fabaceae) via multivariate morphometrics.

Pediomelum is a genus endemic to North America comprising about 26 species, including the megalanthum complex, which consists of Pediomelummegalanthum and its varieties retrorsum and megalanthum, Pediomelummephiticum, and the recently described Pediomelumverdiense and Pediomelumpauperitense. Historically, species of the megalanthum complex have been variably recognized at the species or variety levels, dependent upon the relative importance of morphological characters as diagnostic of species. Ten quantitative morphological characters regarded as diagnostic at the species level were analyzed using multivariate morphometrics across these taxa in order to examine the discriminatory power of these characters to delineate species and to aid in species delimitation. The analyses support the recognition of Pediomelummegalanthum, Pediomelummephiticum, and Pediomelumverdiense at the species level, Pediomelumretrorsum as a variety under Pediomelummegalanthum, and suggest the sinking of Pediomelumpauperitense into Pediomelumverdiense. The findings of the present study help quantify the power of certain characters at delimiting taxa and provide a basis for taxonomic revision of the Pediomelummegalanthum complex.

The wild side of a major crop: soybean's perennial cousins from Down Under.

The accumulation of over 30 years of basic research on the biology, genetic variation, and evolution of the wild perennial relatives of soybean (Glycine max) provides a foundation to improve cultivated soybean. The cultivated soybean and its wild progenitor, G. soja, have a center of origin in eastern Asia and are the only two species in the annual subgenus Soja. Systematic and evolutionary studies of the ca. 30 perennial species of subgenus Glycine, native to Australia, have benefited from the availability of the G. max genomic sequence. The perennial species harbor many traits of interest to soybean breeders, among them resistance to major soybean pathogens such as cyst nematode and leaf rust. New species in the Australian subgenus continue to be described, due to the collection of new material and to insights gleaned through systematic studies of accessions in germplasm collections. Ongoing studies in perennial species focus on genomic regions that contain genes for key traits relevant to soybean breeding. These comparisons also include the homoeologous regions that are the result of polyploidy in the common ancestor of all Glycine species. Subgenus Glycine includes a complex of recently formed allopolyploids that are the focus of studies aimed at elucidating genomic, transcriptomic, physiological, taxonomic, morphological, developmental, and ecological processes related to polyploid evolution. Here we review what has been learned over the past 30 years and outline ongoing work on photosynthesis, nitrogen fixation, and floral biology, much of it drawing on new technologies and resources.

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.

Adaptive horizontal transfer of a bacterial gene to an invasive insect pest of coffee.

Horizontal gene transfer (HGT) involves the nonsexual transmission of genetic material across species boundaries. Although often detected in prokaryotes, examples of HGT involving animals are relatively rare, and any evolutionary advantage conferred to the recipient is typically obscure. We identified a gene (HhMAN1) from the coffee berry borer beetle, Hypothenemus hampei, a devastating pest of coffee, which shows clear evidence of HGT from bacteria. HhMAN1 encodes a mannanase, representing a class of glycosyl hydrolases that has not previously been reported in insects. Recombinant HhMAN1 protein hydrolyzes coffee berry galactomannan, the major storage polysaccharide in this species and the presumed food of H. hampei. HhMAN1 was found to be widespread in a broad biogeographic survey of H. hampei accessions, indicating that the HGT event occurred before radiation of the insect from West Africa to Asia and South America. However, the gene was not detected in the closely related species H. obscurus (the tropical nut borer or "false berry borer"), which does not colonize coffee beans. Thus, HGT of HhMAN1 from bacteria represents a likely adaptation to a specific ecological niche and may have been promoted by intensive agricultural practices.

Applications of next-generation sequencing in plant biology.

The last several years have seen revolutionary advances in DNA sequencing technologies with the advent of next-generation sequencing (NGS) techniques. NGS methods now allow millions of bases to be sequenced in one round, at a fraction of the cost relative to traditional Sanger sequencing. As costs and capabilities of these technologies continue to improve, we are only beginning to see the possibilities of NGS platforms, which are developing in parallel with online availability of a wide range of biological data sets and scientific publications and allowing us to address a variety of questions not possible before. As techniques and data sets continue to improve and grow, we are rapidly moving to the point where every organism, not just select "model organisms", is open to the power of NGS. This volume presents a brief synopsis of NGS technologies and the development of exemplary applications of such methods in the fields of molecular marker development, hybridization and introgression, transcriptome investigations, phylogenetic and ecological studies, polyploid genetics, and applications for large genebank collections.

A comparison of global, gene-specific, and relaxed clock methods in a comparative genomics framework: dating the polyploid history of soybean (Glycine max).

It is widely recognized that many genes and lineages do not adhere to a molecular clock, yet molecular clocks are commonly used to date divergences in comparative genomic studies. We test the application of a molecular clock across genes and lineages in a phylogenetic framework utilizing 12 genes linked in a 1-Mb region on chromosome 13 of soybean (Glycine max); homoeologous copies of these genes formed by polyploidy in Glycine; and orthologous copies in G. tomentella, Phaseolus vulgaris, and Medicago truncatula. We compare divergence dates estimated by two methods each in three frameworks: a global molecular clock with a single rate across genes and lineages using full and approximate likelihood methods based on synonymous substitutions, a gene-specific clock assuming rate constancy over lineages but allowing a different rate for each gene, and a relaxed molecular clock where rates may vary across genes and lineages estimated under penalized likelihood and Bayesian inference. We use the cumulative variance across genes as a means of quantifying precision. Our results suggest that divergence dating methods produce results that are correlated, but that older nodes are more variable and more difficult to estimate with precision and accuracy. We also find that models incorporating less rate heterogeneity estimate older dates of divergence than more complex models, as node age increases. A mixed model nested analysis of variance testing the effects of framework, method, and gene found that framework had a significant effect on the divergence date estimates but that most variation among dates is due to variation among genes, suggesting a need to further characterize and understand the evolutionary phenomena underlying rate variation within genomes, among genes, and across lineages.

Dating the origins of polyploidy events.

is a widespread speciation mechanism, particularly in plants. Estimating the time of origin of polyploid species is important for understanding issues such as gene loss and changes in regulation and expression among homoeologous copies that coexist in a single genome owing to polyploidy. Polyploid species can originate in various ways; the effects of mode of origin, genetic system, and sampling on estimates of the age of polyploid origin using distances between alleles of polyploids and their diploid progenitors, or between homoeologous loci in a polyploid genome, are explored. Even in the simplest cases, simulations confirm that different loci are expected to give very different estimates of the date of origin. The time of polyploid origin is at least as old as the time estimated from comparison of an allele sampled from the polyploid with the most closely related allele in the diploid progenitor. The polyploidy literature often does not make clear the longstanding observation that the divergence of homoeologous copies in an allopolyploid tracks the divergence of diploid species, not the origin of the polyploid. Estimating the date of origin of a polyploid is difficult, and in some circumstances impossible. Skepticism about dates of polyploid origins is clearly warranted.

Divergence and diversification in North American Psoraleeae (Fabaceae) due to climate change.

BACKGROUND: Past studies in the legume family (Fabaceae) have uncovered several evolutionary trends including differential mutation and diversification rates across varying taxonomic levels. The legume tribe Psoraleeae is shown herein to exemplify these trends at the generic and species levels. This group includes a sizable diversification within North America dated at approximately 6.3 million years ago with skewed species distribution to the most recently derived genus, Pediomelum, suggesting a diversification rate shift. We estimate divergence dates of North American (NAm) Psoraleeae using Bayesian MCMC sampling in BEAST based on eight DNA regions (ITS, waxy, matK, trnD-trnT, trnL-trnF, trnK, trnS-trnG, and rpoB-trnC). We also test the hypothesis of a diversification rate shift within NAm Psoraleeae using topological and temporal methods. We investigate the impact of climate change on diversification in this group by (1) testing the hypothesis that a shift from mesic to xeric habitats acted as a key innovation and (2) investigating diversification rate shifts along geologic time, discussing the impact of Quaternary climate oscillations on diversification. RESULTS: NAm Psoraleeae represents a recent, rapid radiation with several genera originating during the Pleistocene, 1 to 2 million years ago. A shift in diversification rate is supported by both methods with a 2.67-fold increase suggested around 2 million years ago followed by a 8.73-fold decrease 440,000 years ago. The hypothesis that a climate regime shift from mesic to xeric habitats drove increased diversification in affected taxa was not supported. Timing of the diversification rate increase supports the hypothesis that glaciation-induced climate changes during the Quaternary influenced diversification of the group. Nonrandom spatial diversification also exists, with greater species richness in the American Southwest. CONCLUSION: This study outlines NAm Psoraleeae as a model example of a recent, rapid radiation. Diversification rate shifts in NAm Psoraleeae are not due to current climate regimes as represented by habitat, but instead to past global climate change resulting from Quaternary glaciations. NAm Psoraleeae diversification is a good example of how earthly dynamics including global climate change and topography work together to shape biodiversity.

Replication of nonautonomous retroelements in soybean appears to be both recent and common.

Retrotransposons and their remnants often constitute more than 50% of higher plant genomes. Although extensively studied in monocot crops such as maize (Zea mays) and rice (Oryza sativa), the impact of retrotransposons on dicot crop genomes is not well documented. Here, we present an analysis of retrotransposons in soybean (Glycine max). Analysis of approximately 3.7 megabases (Mb) of genomic sequence, including 0.87 Mb of pericentromeric sequence, uncovered 45 intact long terminal repeat (LTR)-retrotransposons. The ratio of intact elements to solo LTRs was 8:1, one of the highest reported to date in plants, suggesting that removal of retrotransposons by homologous recombination between LTRs is occurring more slowly in soybean than in previously characterized plant species. Analysis of paired LTR sequences uncovered a low frequency of deletions relative to base substitutions, indicating that removal of retrotransposon sequences by illegitimate recombination is also operating more slowly. Significantly, we identified three subfamilies of nonautonomous elements that have replicated in the recent past, suggesting that retrotransposition can be catalyzed in trans by autonomous elements elsewhere in the genome. Analysis of 1.6 Mb of sequence from Glycine tomentella, a wild perennial relative of soybean, uncovered 23 intact retroelements, two of which had accumulated no mutations in their LTRs, indicating very recent insertion. A similar pattern was found in 0.94 Mb of sequence from Phaseolus vulgaris (common bean). Thus, autonomous and nonautonomous retrotransposons appear to be both abundant and active in Glycine and Phaseolus. The impact of nonautonomous retrotransposon replication on genome size appears to be much greater than previously appreciated.