Interchromosomal linkage disequilibrium and linked fitness cost loci associated with selection for herbicide resistance

The adaptation of weeds to herbicide is both a significant problem in agriculture and a model of rapid adaptation. However, significant gaps remain in our knowledge of resistance controlled by many loci and the evolutionary factors that influence the maintenance of resistance. Here, using herbicide-resistant populations of the common morning glory (Ipomoea purpurea), we perform a multilevel analysis of the genome and transcriptome to uncover putative loci involved in nontarget-site herbicide resistance (NTSR) and to examine evolutionary forces underlying the maintenance of resistance in natural populations. We found loci involved in herbicide detoxification and stress sensing to be under selection and confirmed that detoxification is responsible for glyphosate (RoundUp) resistance using a functional assay. We identified interchromosomal linkage disequilibrium (ILD) among loci under selection reflecting either historical processes or additive effects leading to the resistance phenotype. We further identified potential fitness cost loci that were strongly linked to resistance alleles, indicating the role of genetic hitchhiking in maintaining the cost. Overall, our work suggests that NTSR glyphosate resistance in I. purpurea is conferred by multiple genes which are potentially maintained through generations via ILD, and that the fitness cost associated with resistance in this species is likely a by-product of genetic hitchhiking.

From comparative gene content and gene order to ancestral contigs, chromosomes and karyotypes

To reconstruct the ancestral genome of a set of phylogenetically related descendant species, we use the RACCROCHE pipeline for organizing a large number of generalized gene adjacencies into contigs and then into chromosomes. Separate reconstructions are carried out for each ancestral node of the phylogenetic tree for focal taxa. The ancestral reconstructions are monoploids; they each contain at most one member of each gene family constructed from descendants, ordered along the chromosomes. We design and implement a new computational technique for solving the problem of estimating the ancestral monoploid number of chromosomes x. This involves a "g-mer" analysis to resolve a bias due long contigs, and gap statistics to estimate x. We find that the monoploid number of all the rosid and asterid orders is [Formula: see text]. We show that this is not an artifact of our method by deriving [Formula: see text] for the metazoan ancestor.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

A derived ZW chromosome system in Amborella trichopoda, representing the sister lineage to all other extant flowering plants

The genetic basis and evolution of sex determination in dioecious plants is emerging as an active area of research with exciting advances in genome sequencing and analysis technologies. As the sole species within the sister lineage to all other extant flowering plants, Amborella trichopoda is an important model for understanding the evolution and development of flowers. Plants typically produce only male or female flowers, but sex determination mechanisms are unknown for the species. Sequence data derived from plants of natural origin and an F1 mapping population were used to identify sex-linked genes and the nonrecombining region. Amborella trichopoda has a ZW sex determination system. Analysis of genes in a 4 Mb nonrecombining sex-determination region reveals recent divergence of Z and W gametologs, and few Z- and W-specific genes. The sex chromosomes of A. trichopoda evolved less than 16.5 Myr ago, long after the divergence of the extant angiosperms.

Monocot plastid phylogenomics, timeline, net rates of species diversification, the power of multi-gene analyses, and a functional model for the origin of monocots

PREMISE OF THE STUDY

We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification.

METHODS

We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth.

KEY RESULTS

Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches.

CONCLUSIONS

We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

Brassicales phylogeny inferred from 72 plastid genes: A reanalysis of the phylogenetic localization of two paleopolyploid events and origin of novel chemical defenses

PREMISE OF THE STUDY

Previous phylogenetic studies employing molecular markers have yielded various insights into the evolutionary history across Brassicales, but many relationships between families remain poorly supported or unresolved. A recent phylotranscriptomic approach utilizing 1155 nuclear markers obtained robust estimates for relationships among 14 of 17 families. Here we report a complete family-level phylogeny estimated using the plastid genome.

METHODS

We conducted phylogenetic analyses on a concatenated data set comprising 44,926 bp from 72 plastid genes for species distributed across all 17 families. Our analysis includes three additional families, Tovariaceae, Salvadoraceae, and Setchellanthaceae, that were omitted in the previous phylotranscriptomic study.

KEY RESULTS

Our phylogenetic analyses obtained fully resolved and strongly supported estimates for all nodes across Brassicales. Importantly, these findings are congruent with the topology reported in the phylotranscriptomic study. This consistency suggests that future studies could utilize plastid genomes as markers for resolving relationships within some notoriously difficult clades across Brassicales. We used this new phylogenetic framework to verify the placement of the At-α event near the origin of Brassicaceae, with median date estimates of 31.8 to 42.8 million years ago and restrict the At-β event to one of two nodes with median date estimates between 85 to 92.2 million years ago. These events ultimately gave rise to novel chemical defenses and are associated with subsequent shifts in net diversification rates.

CONCLUSIONS

We anticipate that these findings will aid future comparative evolutionary studies across Brassicales, including selecting candidates for whole-genome sequencing projects.

Shifts in gene expression profiles are associated with weak and strong Crassulacean acid metabolism

PREMISE OF THE STUDY

The relative ease of high throughput sequencing is facilitating comprehensive phylogenomic and gene expression studies, even for nonmodel groups. To date, however, these two approaches have not been merged; while phylogenomic methods might use transcriptome sequences to resolve relationships, assessment of gene expression patterns in a phylogenetic context is less common. Here we analyzed both carbon assimilation and gene expression patterns of closely related species within the Agavoideae (Asparagaceae) to elucidate changes in gene expression across weak and strong phenotypes for Crassulacean acid metabolism (CAM).

METHODS

Gene expression patterns were compared across four genera: Agave (CAM), which is paraphyletic with Polianthes (weak CAM) and Manfreda (CAM), and Beschorneria (weak CAM). RNA-sequencing was paired with measures of gas exchange and titratable acidity. Climate niche space was compared across the four lineages to examine abiotic factors and their correlation to CAM.

KEY RESULTS

Expression of homologous genes showed both shared and variable patterns in weak and strong CAM species. Network analysis highlights that despite shared expression patterns, highly connected genes differ between weak and strong CAM, implicating shifts in regulatory gene function as key for the evolution of CAM. Variation in carbohydrate metabolism between weak and strong CAM supports the importance of sugar turnovers for CAM physiology.

CONCLUSIONS

Integration of phylogenetics and RNA-sequencing provides a powerful tool to study the evolution of CAM photosynthesis across closely related but photosynthetically variable species. Our findings regarding shared or shifted gene expression and regulation of CAM via carbohydrate metabolism have important implications for efforts to engineer the CAM pathway into C₃ food and biofuel crops.

Insights into the Evolution of Hydroxyproline-Rich Glycoproteins from 1000 Plant Transcriptomes

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.

Concomitant loss of NDH complex-related genes within chloroplast and nuclear genomes in some orchids

The chloroplast NAD(P)H dehydrogenase-like (NDH) complex consists of about 30 subunits from both the nuclear and chloroplast genomes and is ubiquitous across most land plants. In some orchids, such as Phalaenopsis equestris, Dendrobium officinale and Dendrobium catenatum, most of the 11 chloroplast genome-encoded ndh genes (cp-ndh) have been lost. Here we investigated whether functional cp-ndh genes have been completely lost in these orchids or whether they have been transferred and retained in the nuclear genome. Further, we assessed whether both cp-ndh genes and nucleus-encoded NDH-related genes can be lost, resulting in the absence of the NDH complex. Comparative analyses of the genome of Apostasia odorata, an orchid species with a complete complement of cp-ndh genes which represents the sister lineage to all other orchids, and three published orchid genome sequences for P. equestris, D. officinale and D. catenatum, which are all missing cp-ndh genes, indicated that copies of cp-ndh genes are not present in any of these four nuclear genomes. This observation suggests that the NDH complex is not necessary for some plants. Comparative genomic/transcriptomic analyses of currently available plastid genome sequences and nuclear transcriptome data showed that 47 out of 660 photoautotrophic plants and all the heterotrophic plants are missing plastid-encoded cp-ndh genes and exhibit no evidence for maintenance of a functional NDH complex. Our data indicate that the NDH complex can be lost in photoautotrophic plant species. Further, the loss of the NDH complex may increase the probability of transition from a photoautotrophic to a heterotrophic life history.

Orchid phylogenomics and multiple drivers of their extraordinary diversification

Orchids are the most diverse family of angiosperms, with over 25 000 species,more than mammals, birds and reptiles combined. Tests of hypotheses to account for such diversity have been stymied by the lack of a fully resolved broad-scale phylogeny. Here,we provide such a phylogeny, based on 75 chloroplast genes for 39 species representing all orchid subfamilies and 16 of 17 tribes, time-calibrated against 17 angiosperm fossils. Asupermatrix analysis places an additional 144 species based on three plastid genes. Orchids appear to have arisen roughly 112 million years ago (Mya); the subfamilies Orchidoideae and Epidendroideae diverged from each other at the end of the Cretaceous; and the eight tribes and three previously unplaced subtribes of the upper epidendroids diverged rapidly from each other between 37.9 and 30.8 Mya. Orchids appear to have undergone one significant acceleration of net species diversification in the orchidoids, and two accelerations and one deceleration in the upper epidendroids. Consistent with theory, such accelerations were correlated with the evolution of pollinia, the epiphytic habit, CAM photosynthesis, tropical distribution (especially in extensive cordilleras),and pollination via Lepidoptera or euglossine bees. Deceit pollination appears to have elevated the number of orchid species by one-half but not via acceleration of the rate of net diversification. The highest rate of net species diversification within the orchids (0.382 sp sp(-1) My(-1)) is 6.8 times that at the Asparagales crown.

Phylotranscriptomic analysis of the origin and early diversification of land plants

Reconstructing the origin and evolution of land plants and their algal relatives is a fundamental problem in plant phylogenetics, and is essential for understanding how critical adaptations arose, including the embryo, vascular tissue, seeds, and flowers. Despite advances in molecular systematics, some hypotheses of relationships remain weakly resolved. Inferring deep phylogenies with bouts of rapid diversification can be problematic; however, genome-scale data should significantly increase the number of informative characters for analyses. Recent phylogenomic reconstructions focused on the major divergences of plants have resulted in promising but inconsistent results. One limitation is sparse taxon sampling, likely resulting from the difficulty and cost of data generation. To address this limitation, transcriptome data for 92 streptophyte taxa were generated and analyzed along with 11 published plant genome sequences. Phylogenetic reconstructions were conducted using up to 852 nuclear genes and 1,701,170 aligned sites. Sixty-nine analyses were performed to test the robustness of phylogenetic inferences to permutations of the data matrix or to phylogenetic method, including supermatrix, supertree, and coalescent-based approaches, maximum-likelihood and Bayesian methods, partitioned and unpartitioned analyses, and amino acid versus DNA alignments. Among other results, we find robust support for a sister-group relationship between land plants and one group of streptophyte green algae, the Zygnematophyceae. Strong and robust support for a clade comprising liverworts and mosses is inconsistent with a widely accepted view of early land plant evolution, and suggests that phylogenetic hypotheses used to understand the evolution of fundamental plant traits should be reevaluated.

Young, intact and nested retrotransposons are abundant in the onion and asparagus genomes

BACKGROUND AND AIMS

Although monocotyledonous plants comprise one of the two major groups of angiosperms and include >65 000 species, comprehensive genome analysis has been focused mainly on the Poaceae (grass) family. Due to this bias, most of the conclusions that have been drawn for monocot genome evolution are based on grasses. It is not known whether these conclusions apply to many other monocots.

METHODS

To extend our understanding of genome evolution in the monocots, Asparagales genomic sequence data were acquired and the structural properties of asparagus and onion genomes were analysed. Specifically, several available onion and asparagus bacterial artificial chromosomes (BACs) with contig sizes >35 kb were annotated and analysed, with a particular focus on the characterization of long terminal repeat (LTR) retrotransposons.

KEY RESULTS

The results reveal that LTR retrotransposons are the major components of the onion and garden asparagus genomes. These elements are mostly intact (i.e. with two LTRs), have mainly inserted within the past 6 million years and are piled up into nested structures. Analysis of shotgun genomic sequence data and the observation of two copies for some transposable elements (TEs) in annotated BACs indicates that some families have become particularly abundant, as high as 4-5 % (asparagus) or 3-4 % (onion) of the genome for the most abundant families, as also seen in large grass genomes such as wheat and maize.

CONCLUSIONS

Although previous annotations of contiguous genomic sequences have suggested that LTR retrotransposons were highly fragmented in these two Asparagales genomes, the results presented here show that this was largely due to the methodology used. In contrast, this current work indicates an ensemble of genomic features similar to those observed in the Poaceae.

Evaluating methods for isolating total RNA and predicting the success of sequencing phylogenetically diverse plant transcriptomes

Next-generation sequencing plays a central role in the characterization and quantification of transcriptomes. Although numerous metrics are purported to quantify the quality of RNA, there have been no large-scale empirical evaluations of the major determinants of sequencing success. We used a combination of existing and newly developed methods to isolate total RNA from 1115 samples from 695 plant species in 324 families, which represents >900 million years of phylogenetic diversity from green algae through flowering plants, including many plants of economic importance. We then sequenced 629 of these samples on Illumina GAIIx and HiSeq platforms and performed a large comparative analysis to identify predictors of RNA quality and the diversity of putative genes (scaffolds) expressed within samples. Tissue types (e.g., leaf vs. flower) varied in RNA quality, sequencing depth and the number of scaffolds. Tissue age also influenced RNA quality but not the number of scaffolds ≥1000 bp. Overall, 36% of the variation in the number of scaffolds was explained by metrics of RNA integrity (RIN score), RNA purity (OD 260/230), sequencing platform (GAIIx vs HiSeq) and the amount of total RNA used for sequencing. However, our results show that the most commonly used measures of RNA quality (e.g., RIN) are weak predictors of the number of scaffolds because Illumina sequencing is robust to variation in RNA quality. These results provide novel insight into the methods that are most important in isolating high quality RNA for sequencing and assembling plant transcriptomes. The methods and recommendations provided here could increase the efficiency and decrease the cost of RNA sequencing for individual labs and genome centers.

Quality and quantity of data recovered from massively parallel sequencing: Examples in Asparagales and Poaceae

PREMISE OF THE STUDY

Genome survey sequences (GSS) from massively parallel sequencing have potential to provide large, cost-effective data sets for phylogenetic inference, replace single gene or spacer regions as DNA barcodes, and provide a plethora of data for other comparative molecular evolution studies. Here we report on the application of this method to estimating the molecular phylogeny of core Asparagales, investigating plastid gene losses, assembling complete plastid genomes, and determining the type and quality of assembled genomic data attainable from Illumina 80-120-bp reads.

METHODS

We sequenced total genomic DNA from samples in two lineages of monocotyledonous plants, Poaceae and Asparagales, on the Illumina platform in a multiplex arrangement. We compared reference-based assemblies to de novo contigs, evaluated consistency of assemblies resulting from use of various references sequences, and assessed our methods to obtain sequence assemblies in nonmodel taxa.

KEY RESULTS

Our method returned reliable, robust organellar and nrDNA sequences in a variety of plant lineages. High quality assemblies are not dependent on genome size, amount of plastid present in the total genomic DNA template, or relatedness of available reference sequences for assembly. Phylogenetic results revealed familial and subfamilial relationships within Asparagales with high bootstrap support, although placement of the monotypic genus Aphyllanthes was placed with moderate confidence.

CONCLUSIONS

The well-supported molecular phylogeny provides evidence for delineation of subfamilies within core Asparagales. With advances in technology and bioinformatics tools, the use of massively parallel sequencing will continue to become easier and more affordable for phylogenomic and molecular evolutionary biology investigations.

The iPlant Collaborative: Cyberinfrastructure for Plant Biology

The iPlant Collaborative (iPlant) is a United States National Science Foundation (NSF) funded project that aims to create an innovative, comprehensive, and foundational cyberinfrastructure in support of plant biology research (PSCIC, 2006). iPlant is developing cyberinfrastructure that uniquely enables scientists throughout the diverse fields that comprise plant biology to address Grand Challenges in new ways, to stimulate and facilitate cross-disciplinary research, to promote biology and computer science research interactions, and to train the next generation of scientists on the use of cyberinfrastructure in research and education. Meeting humanity's projected demands for agricultural and forest products and the expectation that natural ecosystems be managed sustainably will require synergies from the application of information technologies. The iPlant cyberinfrastructure design is based on an unprecedented period of research community input, and leverages developments in high-performance computing, data storage, and cyberinfrastructure for the physical sciences. iPlant is an open-source project with application programming interfaces that allow the community to extend the infrastructure to meet its needs. iPlant is sponsoring community-driven workshops addressing specific scientific questions via analysis tool integration and hypothesis testing. These workshops teach researchers how to add bioinformatics tools and/or datasets into the iPlant cyberinfrastructure enabling plant scientists to perform complex analyses on large datasets without the need to master the command-line or high-performance computational services.

Gene loss under neighborhood selection following whole genome duplication and the reconstruction of the ancestral Populus genome

We develop criteria to detect neighborhood selection effects on gene loss following whole genome duplication, and apply them to the recently sequenced poplar (Populus trichocarpa) genome. We improve on guided genome halving algorithms so that several thousand gene sets, each containing two paralogs in the descendant T of the doubling event and their single ortholog from an undoubled reference genome R, can be analyzed to reconstruct the ancestor A of T at the time of doubling. At the same time, large numbers of defective gene sets, either missing one paralog from T or missing their ortholog in R, may be incorporated into the analysis in a consistent way. We apply this genomic rearrangement distance-based approach to the poplar and grapevine (Vitis vinifera) genomes, as T and R respectively. We conclude that, after chromosome doubling, the "choice" of which paralogous gene pairs will lose copies is random, but that the retention of strings of single-copy genes on one chromosome versus the other is decidedly non-random.

Pattern and timing of diversification in Yucca (Agavaceae): specialized pollination does not escalate rates of diversification

The yucca-yucca moth interaction is one of the most well-known and remarkable obligate pollination mutualisms, and is an important study system for understanding coevolution. Previous research suggests that specialist pollinators can promote rapid diversification in plants, and theoretical work has predicted that obligate pollination mutualism promotes cospeciation between plants and their pollinators, resulting in contemporaneous, parallel diversification. However, a lack of information about the age of Yucca has impeded efforts to test these hypotheses. We used analyses of 4322 AFLP markers and cpDNA sequence data representing six non-protein-coding regions (trnT-trnL, trnL, trnL intron, trnL-trnF, rps16 and clpP intron 2) from all 34 species to recover a consensus organismal phylogeny, and used penalized likelihood to estimate divergence times and speciation rates in Yucca. The results indicate that the pollination mutualism did not accelerate diversification, as Yucca diversity (34 species) is not significantly greater than that of its non-moth-pollinated sister group, Agave sensu latissimus (240 species). The new phylogenetic estimates also corroborate the suggestion that the plant-moth pollination mutualism has at least two origins within the Agavaceae. Finally, age estimates show significant discord between the age of Yucca (ca 6-10Myr) and the current best estimates for the age of their pollinators (32-40Myr).

Analysis of 81 genes from 64 plastid genomes resolves relationships in angiosperms and identifies genome-scale evolutionary patterns

Angiosperms are the largest and most successful clade of land plants with >250,000 species distributed in nearly every terrestrial habitat. Many phylogenetic studies have been based on DNA sequences of one to several genes, but, despite decades of intensive efforts, relationships among early diverging lineages and several of the major clades remain either incompletely resolved or weakly supported. We performed phylogenetic analyses of 81 plastid genes in 64 sequenced genomes, including 13 new genomes, to estimate relationships among the major angiosperm clades, and the resulting trees are used to examine the evolution of gene and intron content. Phylogenetic trees from multiple methods, including model-based approaches, provide strong support for the position of Amborella as the earliest diverging lineage of flowering plants, followed by Nymphaeales and Austrobaileyales. The plastid genome trees also provide strong support for a sister relationship between eudicots and monocots, and this group is sister to a clade that includes Chloranthales and magnoliids. Resolution of relationships among the major clades of angiosperms provides the necessary framework for addressing numerous evolutionary questions regarding the rapid diversification of angiosperms. Gene and intron content are highly conserved among the early diverging angiosperms and basal eudicots, but 62 independent gene and intron losses are limited to the more derived monocot and eudicot clades. Moreover, a lineage-specific correlation was detected between rates of nucleotide substitutions, indels, and genomic rearrangements.

The phylogeny of yuccas

The genus Yucca is widely recognized for its pollination mutualism with yucca moths. Analysis of diversification in this interaction has been hampered by the lack of a robust phylogeny for the genus. Here we attempt the first extensive nuclear DNA based assessment of the phylogenetic relationships of Yucca. We used AFLP markers to recover the phylogeny of 87 samples representing 38 Yucca taxa. An analysis based on 4322 markers strongly supported a topology consistent with morphological classification at the section level (capsular-fruited Chaenocarpa, fleshy-fruited Sarcocarpa, and spongy-fruited Clistocarpa). Within Sarcocarpa, all but two of the traditional species were monophyletic. Within Chaenocarpa, the morphologically distinct series Rupicolae was strongly supported. In the remaining Chaenocarpa, a western group (Colorado Plateau southward) and an eastern group (Great Plains, central Texas east to Florida) were recovered. Within these groups, where taxonomic circumscriptions are narrow and historically contested, there was at most limited monophyly of traditional taxa, suggesting rapid recent diversification, introgression, or non-monophyletically circumscribed taxa.

Comparative genomic analyses in Asparagus

Garden asparagus (Asparagus officinalis L.) belongs to the monocot family Asparagaceae in the order Asparagales. Onion (Allium cepa L.) and Asparagus officinalis are 2 of the most economically important plants of the core Asparagales, a well supported monophyletic group within the Asparagales. Coding regions in onion have lower GC contents than the grasses. We compared the GC content of 3374 unique expressed sequence tags (ESTs) from A. officinalis with Lycoris longituba and onion (both members of the core Asparagales), Acorus americanus (sister to all other monocots), the grasses, and Arabidopsis. Although ESTs in A. officinalis and Acorus had a higher average GC content than Arabidopsis, Lycoris, and onion, all were clearly lower than the grasses. The Asparagaceae have the smallest nuclear genomes among all plants in the core Asparagales, which typically have huge genomes. Within the Asparagaceae, European Asparagus species have approximately twice the nuclear DNA of that of southern African Asparagus species. We cloned and sequenced 20 genomic amplicons from European A. officinalis and the southern African species Asparagus plumosus and observed no clear evidence for a recent genome doubling in A. officinalis relative to A. plumosus. These results indicate that members of the genus Asparagus with smaller genomes may be useful genomic models for plants in the core Asparagales.

Complete plastid genome sequences of Drimys, Liriodendron, and Piper: implications for the phylogenetic relationships of magnoliids

Background

The magnoliids with four orders, 19 families, and 8,500 species represent one of the largest clades of early diverging angiosperms. Although several recent angiosperm phylogenetic analyses supported the monophyly of magnoliids and suggested relationships among the orders, the limited number of genes examined resulted in only weak support, and these issues remain controversial. Furthermore, considerable incongruence resulted in phylogenetic reconstructions supporting three different sets of relationships among magnoliids and the two large angiosperm clades, monocots and eudicots. We sequenced the plastid genomes of three magnoliids, Drimys (Canellales), Liriodendron (Magnoliales), and Piper (Piperales), and used these data in combination with 32 other angiosperm plastid genomes to assess phylogenetic relationships among magnoliids and to examine patterns of variation of GC content.

Results

The Drimys, Liriodendron, and Piper plastid genomes are very similar in size at 160,604, 159,886 bp, and 160,624 bp, respectively. Gene content and order are nearly identical to many other unrearranged angiosperm plastid genomes, including Calycanthus, the other published magnoliid genome. Overall GC content ranges from 34–39%, and coding regions have a substantially higher GC content than non-coding regions. Among protein-coding genes, GC content varies by codon position with 1st codon > 2nd codon > 3rd codon, and it varies by functional group with photosynthetic genes having the highest percentage and NADH genes the lowest. Phylogenetic analyses using parsimony and likelihood methods and sequences of 61 protein-coding genes provided strong support for the monophyly of magnoliids and two strongly supported groups were identified, the Canellales/Piperales and the Laurales/Magnoliales. Strong support is reported for monocots and eudicots as sister clades with magnoliids diverging before the monocot-eudicot split. The trees also provided moderate or strong support for the position of Amborella as sister to a clade including all other angiosperms.

Conclusion

Evolutionary comparisons of three new magnoliid plastid genome sequences, combined with other published angiosperm genomes, confirm that GC content is unevenly distributed across the genome by location, codon position, and functional group. Furthermore, phylogenetic analyses provide the strongest support so far for the hypothesis that the magnoliids are sister to a large clade that includes both monocots and eudicots.

The genome of black cottonwood, Populus trichocarpa (Torr. & Gray)

We report the draft genome of the black cottonwood tree, Populus trichocarpa. Integration of shotgun sequence assembly with genetic mapping enabled chromosome-scale reconstruction of the genome. More than 45,000 putative protein-coding genes were identified. Analysis of the assembled genome revealed a whole-genome duplication event; about 8000 pairs of duplicated genes from that event survived in the Populus genome. A second, older duplication event is indistinguishably coincident with the divergence of the Populus and Arabidopsis lineages. Nucleotide substitution, tandem gene duplication, and gross chromosomal rearrangement appear to proceed substantially more slowly in Populus than in Arabidopsis. Populus has more protein-coding genes than Arabidopsis, ranging on average from 1.4 to 1.6 putative Populus homologs for each Arabidopsis gene. However, the relative frequency of protein domains in the two genomes is similar. Overrepresented exceptions in Populus include genes associated with lignocellulosic wall biosynthesis, meristem development, disease resistance, and metabolite transport.

Patterns of Speciation in the Yucca Moths: Parallel Species Radiations within the Tegeticula yuccasella Species Complex

The interaction between yuccas and yucca moths has been central to understanding the origin and loss of obligate mutualism and mutualism reversal. Previous systematic research using mtDNA sequence data and characters associated with genitalic morphology revealed that a widespread pollinator species in the genus Tegeticula was in fact a complex of pollinator species that differed in host use and the placement of eggs into yucca flowers. Within this mutualistic clade two nonpollinating "cheater" species evolved. Cheaters feed on yucca seeds but lack the tentacular mouthparts necessary for yucca pollination. Previous work suggested that the species complex formed via a rapid radiation within the last several million years. In this study, we use an expanded mtDNA sequence data set and AFLP markers to examine the phylogenetic relationships among this rapidly diverging clade of moths and compare these relationships to patterns in genitalic morphology. Topologies obtained from analyses of the mtDNA and AFLP data differed significantly. Both data sets, however, corroborated the hypothesis of a rapid species radiation and suggested that there were likely two independent species radiations. Morphological analyses based on oviposition habit produced species groupings more similar to the AFLP topology than the mtDNA topology and suggested the two radiations coincided with differences in oviposition habit. The evolution of cheating was reaffirmed to have evolved twice and the closest pollinating relative for one cheater species was identified by both mtDNA and AFLP markers. For the other cheater species, however, the closest pollinating relative remains ambiguous, and mtDNA, AFLP, and morphological data suggest this cheater species may be diverged based on host use. Much of the divergence in the species complex can be explained by geographic isolation associated with the evolution of two oviposition habits.

To B or Not to B a flower: the role of DEFICIENS and GLOBOSA orthologs in the evolution of the angiosperms

DEFICIENS (DEF) and GLOBOSA (GLO) function in petal and stamen organ identity in Antirrhinum and are orthologs of APETALA3 and PISTILLATA in Arabidopsis. These genes are known as B-function genes for their role in the ABC genetic model of floral organ identity. Phylogenetic analyses show that DEF and GLO are closely related paralogs, having originated from a gene duplication event after the separation of the lineages leading to the extant gymnosperms and the extant angiosperms. Several additional gene duplications followed, providing multiple potential opportunities for functional divergence. In most angiosperms studied to date, genes in the DEF/GLO MADS-box subfamily are expressed in the petals and stamens during flower development. However, in some angiosperms, the expression of DEF and GLO orthologs are occasionally observed in the first and fourth whorls of flowers or in nonfloral organs, where their function is unknown. In this article we review what is known about function, phylogeny, and expression in the DEF/GLO subfamily to examine their evolution in the angiosperms. Our analyses demonstrate that although the primary role of the DEF/GLO subfamily appears to be in specifying the stamens and inner perianth, several examples of potential sub- and neofunctionalization are observed.

Methods for obtaining and analyzing whole chloroplast genome sequences

During the past decade, there has been a rapid increase in our understanding of plastid genome organization and evolution due to the availability of many new completely sequenced genomes. There are 45 complete genomes published and ongoing projects are likely to increase this sampling to nearly 200 genomes during the next 5 years. Several groups of researchers including ours have been developing new techniques for gathering and analyzing entire plastid genome sequences and details of these developments are summarized in this chapter. The most important developments that enhance our ability to generate whole chloroplast genome sequences involve the generation of pure fractions of chloroplast genomes by whole genome amplification using rolling circle amplification, cloning genomes into Fosmid or bacterial artificial chromosome (BAC) vectors, and the development of an organellar annotation program (Dual Organellar GenoMe Annotator [DOGMA]). In addition to providing details of these methods, we provide an overview of methods for analyzing complete plastid genome sequences for repeats and gene content, as well as approaches for using gene order and sequence data for phylogeny reconstruction. This explosive increase in the number of sequenced plastid genomes and improved computational tools will provide many insights into the evolution of these genomes and much new data for assessing relationships at deep nodes in plants and other photosynthetic organisms.

EST clustering error evaluation and correction

MOTIVATION

The gene expression intensity information conveyed by (EST) Expressed Sequence Tag data can be used to infer important cDNA library properties, such as gene number and expression patterns. However, EST clustering errors, which often lead to greatly inflated estimates of obtained unique genes, have become a major obstacle in the analyses. The EST clustering error structure, the relationship between clustering error and clustering criteria, and possible error correction methods need to be systematically investigated.

RESULTS

We identify and quantify two types of EST clustering error, namely, Type I and II in EST clustering using CAP3 assembling program. A Type I error occurs when ESTs from the same gene do not form a cluster whereas a Type II error occurs when ESTs from distinct genes are falsely clustered together. While the Type II error rate is <1.5% for both 5' and 3' EST clustering, the Type I error in the 5' EST case is approximately 10 times higher than the 3' EST case (30% versus 3%). An over-stringent identity rule, e.g., P >/= 95%, may even inflate the Type I error in both cases. We demonstrate that approximately 80% of the Type I error is due to insufficient overlap among sibling ESTs (ISO error) in 5' EST clustering. A novel statistical approach is proposed to correct ISO error to provide more accurate estimates of the true gene cluster profile.

Patterns of Genetic Structure among Populations of an Oligophagous Pollinating Yucca Moth (Tegeticula yuccasella)

Plant-insect associations have served as models for investigations of coevolution and the influence of biotic interactions on diversification. The pollination association between yuccas and yucca moths is a classic example of an obligate mutualism often suggested to have been affected by coevolution. Recent work has shown high host specificity in pollinating yucca moths, and here we use Tegeticula yuccasella, the species with the widest diet breadth, to ask how host specificity and isolation by distance contribute to specialization. Isolation by distance at a regional scale was observed in nucleotide variation within the mitochondrial gene cytochrome oxidase I (COI) (r =.294; P =.003). Host-related genetic structure (F(ct) = 0.08) was found to be slightly lower than the level of structure observed between eastern and western moth populations (F(ct) = 0.096). However, 56% of the COI haplotypes sampled from moths on Yucca filamentosa mapped to a host-specific clade in the haplotype network. Taken together, these results suggest that differentiation among T. yuccasella populations on alternative hosts is slight, but gene flow is influenced by both host association and geographic distance.

Pollen dispersal in Yucca filamentosa (Agavaceae): the paradox of self-pollination behavior by Tegeticula yuccasella (Prodoxidae)

We investigated pollen dispersal in an obligate pollination mutualism between Yucca filamentosa and Tegeticula yuccasella. Yucca moths are the only documented pollinator of yuccas, and moth larvae feed solely on developing yucca seeds. The quality of pollination by a female moth affects larval survival because flowers receiving small amounts of pollen or self-pollen have a high abscission probability, and larvae die in abscised flowers. We tested the prediction that yucca moths primarily perform outcross pollinations by using fluorescent dye to track pollen dispersal in five populations of Y. filamentosa. Dye transfers within plants were common in all populations (mean ± 1 SE, 55 ± 3.0%), indicating that moths frequently deposit self-pollen. Distance of dye transfers ranged from 0 to 50 m, and the mean number of flowering plants between the pollen donor and recipient was 5 (median = 0), suggesting that most pollen was transferred among near neighbors. A multilocus genetic estimate of outcrossing based on seedlings matured from open-pollinated fruits at one site was 94 ± 6% (mean ± 1 SD). We discuss why moths frequently deposit self-pollen to the detriment of their offspring and compare the yucca-yucca moth interaction with other obligate pollinator mutualisms in which neither pollinator nor plant benefit from self-pollination.

Forty million years of mutualism: evidence for eocene origin of the yucca-yucca moth association

The obligate mutualism between yuccas and yucca moths is a major model system for the study of coevolving species interactions. Exploration of the processes that have generated current diversity and associations within this mutualism requires robust phylogenies and timelines for both moths and yuccas. Here we establish a molecular clock for the moths based on mtDNA and use it to estimate the time of major life history events within the yucca moths. Colonization of yuccas had occurred by 41.5 +/- 9.8 million years ago (Mya), with rapid life history diversification and the emergence of pollinators within 0-6 My after yucca colonization. A subsequent burst of diversification 3.2 +/- 1.8 Mya coincided with evolution of arid habitats in western North America. Derived nonpollinating cheater yucca moths evolved 1.26 +/- 0.96 Mya. The estimated age of the moths far predates the host fossil record, but is consistent with suggested host age based on paleobotanical, climatological, biogeographical, and geological data, and a tentative estimation from an rbcL-based molecular clock for yuccas. The moth data are used to establish three alternative scenarios of how the moths and plants have coevolved. They yield specific predictions that can be tested once a robust plant phylogeny becomes available.

Pollination biology in hybridizing Baptisia (Fabaceae) populations

In their classic study, Alston and Turner (American Journal of Botany, vol. 50, 159-173, 1963) documented extensive hybridization among four morphologically distinct Baptisia species native to East Texas. While Alston and Turner found putative F1 hybrids in great numbers, they found no evidence of backcrossing. In this study prezygotic and postzygotic reproductive barriers between two of these species, B. leucophaea and B. sphaerocarpa, were investigated and found to be quite weak. Flowering times overlap and bumble bees were observed visiting both species and intermediate hybrids. While pollinator constancy in flights between B. leucophaea and B. sphaerocarpa was moderately strong, significant levels of constancy were not observed in flights involving hybrids and either parental species. Thus, backcrossing was not impeded by pollinator behavior. Further, hybrid pollen was highly stainable (93.5%) and able to effectively set seed in crossing experiments with both parental species. Pollinator behavior was compared in experimental populations with and without hybrid ramets and found to differ between these two treatments. Hybrids were found to facilitate pollinator movement between species. In total, these results suggest that reproductive isolation is not responsible for the rarity of backcrossing in naturally hybridizing B. leucophaea and B. sphaerocarpa populations.

Nuclear DNA-based markers for plant evolutionary biology

While DNA-based markers can provide a wealth of information for the study of plant evolutionary biology, progress is limited by the lack of primers available for PCR. To overcome this limitation, we outline a protocol for developing oligonucleotide primers targeting regions of low copy-number nuclear genes. This protocol is intended to lead to universally useful primer sets. To test our approach, we designed eight primer sets and tested their abilities to amplify targets from representatives of each dicot and one monocot subclass. Five of the eight primer sets amplified targets from at least five of the seven taxa and thus exhibited broad taxonomic usefulness; the remaining primers were rather specific, however, and amplified targets from at most three taxa. In only one primer-taxon combination was a complex multiple-banded amplification produced. Overall, the protocol outlined proved quite useful at identifying broadly applicable primers targeted to low copy-number nuclear genes. Wider application of this approach should be effective at greatly increasing the amount of genetic information available for a diversity of plant nuclear genomes.

Non-mutualistic yucca moths and their evolutionary consequences

Interspecific mutualisms are regarded as having evolved from antagonistic or commensalistic interactions, with most mutualisms remaining facultative but some having coevolved into obligate reciprocal dependency. Underlying mutualism is an intrinsic conflict between the parties, in that each is under selection for increased exploitation of the other. Theoretical models suggest that this conflict is a source of evolutionary instability, and that evolution of 'cheating' by one party may lead to reciprocal extinction. Here we present phylogenic evidence for reversal of an obligate mutualism: within the yucca moth complex, distinct cheater species derived from obligate pollinators inflict a heavy cost on their yucca hosts by laying their eggs but not pollinating the yucca. Phylogenetic data show the cheaters to have existed for a long time. Coexisting pollinators and cheaters are not sister taxa, supporting predictions that evolution of cheating within a single pollinator is evolutionarily unstable. Several lines of evidence support a hypothesis that host shifts preceded the reversal of obligate mutualism. Host or partner shifts is a mechanism that can provide a route of evolutionary escape among obligate mutualists in general.