Target capture and genome skimming for plant diversity studies.

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Diversity and divergence: evolution of secondary metabolism in the tropical tree genus Inga.

Plants are widely recognized as chemical factories, with each species producing dozens to hundreds of unique secondary metabolites. These compounds shape the interactions between plants and their natural enemies. We explore the evolutionary patterns and processes by which plants generate chemical diversity, from evolving novel compounds to unique chemical profiles. We characterized the chemical profile of one-third of the species of tropical rainforest trees in the genus Inga (c. 100, Fabaceae) using ultraperformance liquid chromatography-mass spectrometry-based metabolomics and applied phylogenetic comparative methods to understand the mode of chemical evolution. We show: each Inga species contain structurally unrelated compounds and high levels of phytochemical diversity; closely related species have divergent chemical profiles, with individual compounds, compound classes, and chemical profiles showing little-to-no phylogenetic signal; at the evolutionary time scale, a species' chemical profile shows a signature of divergent adaptation. At the ecological time scale, sympatric species were the most divergent, implying it is also advantageous to maintain a unique chemical profile from community members; finally, we integrate these patterns with a model for how chemical diversity evolves. Taken together, these results show that phytochemical diversity and divergence are fundamental to the ecology and evolution of plants.

Tracking of Host Defenses and Phylogeny During the Radiation of Neotropical Inga-Feeding Sawflies (Hymenoptera; Argidae).

Coevolutionary theory has long predicted that the arms race between plants and herbivores is a major driver of host selection and diversification. At a local scale, plant defenses contribute significantly to the structure of herbivore assemblages and the high alpha diversity of plants in tropical rain forests. However, the general importance of plant defenses in host associations and divergence at regional scales remains unclear. Here, we examine the role of plant defensive traits and phylogeny in the evolution of host range and species divergence in leaf-feeding sawflies of the family Argidae associated with Neotropical trees in the genus Inga throughout the Amazon, the Guiana Shield and Panama. Our analyses show that the phylogenies of both the sawfly herbivores and their Inga hosts are congruent, and that sawflies radiated at approximately the same time, or more recently than their Inga hosts. Analyses controlling for phylogenetic effects show that the evolution of host use in the sawflies associated with Inga is better correlated with Inga chemistry than with Inga phylogeny, suggesting a pattern of delayed host tracking closely tied to host chemistry. Finally, phylogenetic analyses show that sister species of Inga-sawflies are dispersed across the Neotropics, suggesting a role for allopatric divergence and vicariance in Inga diversification. These results are consistent with the idea that host defensive traits play a key role not only in structuring the herbivore assemblages at a single site, but also in the processes shaping host association and species divergence at a regional scale.

Chemocoding as an identification tool where morphological- and DNA-based methods fall short: Inga as a case study.

The need for species identification and taxonomic discovery has led to the development of innovative technologies for large-scale plant identification. DNA barcoding has been useful, but fails to distinguish among many species in species-rich plant genera, particularly in tropical regions. Here, we show that chemical fingerprinting, or 'chemocoding', has great potential for plant identification in challenging tropical biomes. Using untargeted metabolomics in combination with multivariate analysis, we constructed species-level fingerprints, which we define as chemocoding. We evaluated the utility of chemocoding with species that were defined morphologically and subject to next-generation DNA sequencing in the diverse and recently radiated neotropical genus Inga (Leguminosae), both at single study sites and across broad geographic scales. Our results show that chemocoding is a robust method for distinguishing morphologically similar species at a single site and for identifying widespread species across continental-scale ranges. Given that species are the fundamental unit of analysis for conservation and biodiversity research, the development of accurate identification methods is essential. We suggest that chemocoding will be a valuable additional source of data for a quick identification of plants, especially for groups where other methods fall short.

Transcriptome mining for phylogenetic markers in a recently radiated genus of tropical plants (Renealmia L.f., Zingiberaceae).

The reconstruction of relationships within species-rich groups that have recently evolved in biodiversity hotspots is hampered by a lack of phylogenetically informative markers. It is also made difficult by the lack of sampling necessary to reconstruct a species-level phylogeny. We use transcriptome mining to search for markers to reconstruct a phylogeny of the amphi-Atlantic genus Renealmia L. f. (Zingiberaceae). We recover seven introns from single copy genes and use them to reconstruct the phylogeny of the genus together with a commonly used phylogenetic marker, internal transcribed spacers of ribosomal DNA (ITS) that has previously been used to reconstruct the phylogeny of the genus. We targeted genes with low numbers of base pairs that improves sequencing success using highly degraded DNA from herbarium specimens. The use of herbarium specimens greatly increased the number of species in the study as these were readily available in historical collections. Data were obtained for 14 of the 17 African species and 54 of the 65 Neotropical species. The phylogeny was well-supported for a number of Renealmia subgroups although relationships among those clades remained poorly supported.

Ultrastructure and development of non-contiguous stomatal clusters and helicocytic patterning in Begonia.

Background and Aims: Helicocytic stomata are characterized by an inward spiral of mesogenous cells surrounding a central stomatal pore. They represent a relatively rare feature that occurs in some drought-tolerant angiosperm species. In some Begonia species with thick leaves, the stomata are not only helicocytic but also clustered into groups that are spaced apart by at least one cell. This paper presents a detailed ontogenetic study of this characteristic non-contiguous stomatal patterning in a developmental and phylogenetic context. Methods: Light microscopy and both scanning and transmission electron microscopy were used to examine stomatal development in several species of Begonia. Published reports of stomatal development in Begonia and other angiosperms were reviewed to provide a comprehensive discussion of the evolution of stomatal patterning. Key Results: Helicocytic stomata develop from meristemoids that undergo a series of oriented asymmetric divisions to produce a spiral of mesogene stomatal lineage ground cells (SLGCs) surrounding a stoma. A clear developmental similarity between anisocytic and helicocytic stomata is positively correlated with the number of iterations of amplifying divisions that result in SLGCs. Stomatal clusters develop from asymmetric divisions in neighbouring SLGCs. Within each cluster, non-contiguous spacing of meristemoids is maintained by asymmetric divisions oriented away from each developing meristemoid. Conclusions: Formation of non-contiguous stomatal clusters in Begonia relies on two primary developmental factors in the epidermis: an inwardly spiralling series of amplifying divisions that result in helicocytic stomata, and the development of a variable number of meristemoids from neighbouring SLGCs within each cluster. Optimization of these features on an angiosperm phylogeny indicates that the occurrence of amplifying divisions could be pre-adaptive for these factors. Both factors have been thoroughly studied in terms of developmental genetics in Arabidopsis, suggesting gene orthologues that could be implicated in Begonia stomatal patterning.

Comparative Analysis of Begonia Plastid Genomes and Their Utility for Species-Level Phylogenetics.

Recent, rapid radiations make species-level phylogenetics difficult to resolve. We used a multiplexed, high-throughput sequencing approach to identify informative genomic regions to resolve phylogenetic relationships at low taxonomic levels in Begonia from a survey of sixteen species. A long-range PCR method was used to generate draft plastid genomes to provide a strong phylogenetic backbone, identify fast evolving regions and provide informative molecular markers for species-level phylogenetic studies in Begonia.

Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae).

Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.

Maintenance of species boundaries in a Neotropical radiation of Begonia.

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by-product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later-generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.

The evolution of sex ratio differences and inflorescence architectures in Begonia (Begoniaceae).

PREMISE OF THE STUDY: A major benefit conferred by monoecy is the ability to alter floral sex ratio in response to selection. In monoecious species that produce flowers of a given sex at set positions on the inflorescence, floral sex ratio may be related to inflorescence architecture. We studied the loci underlying differences in inflorescence architecture between two monoecious Begonia species and related this to floral sex ratios. METHODS: We performed trait comparisons and quantitative trait locus (QTL) mapping in a segregating backcross population between Central American Begonia plebeja and B. conchifolia. We focused on traits related to inflorescence architecture, sex ratios, and other reproductive traits. KEY RESULTS: The inflorescence branching pattern of B. conchifolia was more asymmetric than B. plebeja, which in turn affects the floral sex ratio. Colocalizing QTLs of moderate effect influenced both the number of male flowers and the fate decisions of axillary meristems, demonstrating the close link between inflorescence architecture and sex ratio. Additional QTLs were found for stamen number (30% variance explained, VE) and pollen sterility (12.3% VE). CONCLUSIONS: One way in which Begonia species develop different floral sex ratios is through modifications of their inflorescence architecture. The potential pleiotropic action of QTL on inflorescence branching and floral sex ratios may have major implications for trait evolution and responses to selection. The presence of a single QTL of large effect on stamen number may allow rapid divergence for this key floral trait. We propose candidate loci for stamen number and inflorescence branching for future characterization.

Development and characterization of microsatellite markers for Central American Begonia sect. Gireoudia (Begoniaceae).

PREMISE OF THE STUDY: Transcriptome sequence data were used to design microsatellite primers for two widespread Central American Begonia species, B. heracleifolia and B. nelumbiifolia, to investigate population structure and hybridization. • METHODS AND RESULTS: The transcriptome from vegetative meristem tissue from the related B. plebeja was mined for microsatellite loci, and 31 primer pairs amplified in the target species. Fifteen primer pairs were combined in two multiplex PCR reactions, which amplified an average of four alleles per locus. • CONCLUSIONS: The markers developed will be a valuable genetic resource for medium-throughput genotyping of Central American species of Begonia sect. Gireoudia. A subset of these markers have perfect sequence matches to Asian B. venusta, and are promising for studies in other Begonia sections.

Signaling sides adaxial-abaxial patterning in leaves.

Most leaves are dorsiventrally flattened and develop clearly defined upper and lower surfaces. Light capturing is the specialization of the adaxial or upper surface and the abaxial or lower surface is specialized for gas exchange (Fig. 5.1). This division into adaxial and abaxial domains is also key for the outgrowth of the leaf blade or lamina, which occurs along the boundary between the upper and lower sides. How this polarity is set up is not clear but genetic analysis in a range of species suggests that several highly conserved interlocking pathways are involved. Positional information from the meristem is reinforced by signaling through the epidermal layer as the meristem grows away from the leaf primordium. Opposing ta-siRNA and miRNA gradients help refine distinct adaxial and abaxial sides, and mutual inhibition between the genes expressed on each side stabilizes the boundary. In this review we consider how recent work in a range of species is clarifying our understanding of these processes.

The many roles of small RNAs in leaf development.

Leaf development involves many complex genetic interactions, signals between adjacent cells or between more distant tissues and consequent changes in cell fate. This review describes three stages in leaf development where regulation by small RNAs have been used to modulate gene expression patterns.

Mixing and matching pathways in leaf polarity.

Leaves extend a large, porous surface to the environment to catch light and exchange gasses. The extension of the lamina is produced by the interaction of an upper (adaxial) and a lower (abaxial) domain in the developing leaf primordium. Recent studies have revealed that conserved genetic pathways, involving small regulatory RNAs and several distinct transcription factor families, have key roles in adaxial-abaxial patterning, suggesting candidate signals that convey positional information within the shoot to the newly initiated leaf. The interactions of the polarity pathways are distinguished by mutual antagonism and by redundancies. Analysis of these pathways in different model organisms reveals a surprising diversity in the genetic control of such a fundamental developmental process.

The role of ARGONAUTE1 (AGO1) in meristem formation and identity.

The ARGONAUTE gene family is involved in the regulation of gene expression via the RNAi Silencing Complex (RISC). microRNA (miRNA) are 20-22bp RNAs that direct RISC to target genes. Several miRNA have been characterized in plants. Their roles include control of flowering time, floral organ identity, cell division patterns, and leaf polarity. ARGONAUTE1 (AGO1) is required for stem cell function and organ polarity, as is the closely related protein PINHEAD/ZWILLE (PNH/ZLL). Through phenotypic and double mutant analysis, we show that AGO1 regulates stem cell function via SHOOT MERISTEMLESS (STM). CUPSHAPED COTYLEDONS1 and 2 (CUC1 and CUC2) positively regulate STM and are targets of miRNA. The effect of AGO1 on leaf polarity is dependent, in part, on its role in meristem function revealed by interactions with ASYMMETRIC LEAVES1(AS1). AGO1 is required for full expression of LEAFY (LFY), APETALA1 (AP1) and AGAMOUS (AG). Flowering time is unaffected but floral meristem identity is partially restored in a curlyleaf (clf) background and this is not due to clf's affects on AG expression. CLF is over expressed in ago1, showing that the RNAi pathway regulates polycomb-type epigenetic modifiers.

The developmental role of microRNA in plants.

MicroRNAs (miRNAs) are single-stranded RNA molecules of around 22 nucleotides (nt) in length that are associated with the RNA-induced silencing complex (RISC). They play an important role in plant development, either by targeting mRNA for cleavage or by inhibiting translation. Over the past year, the list of known miRNAs, confirmed targets and developmental effects has expanded, as has the realization that they are conserved during evolution and that small RNAs can play a direct role in cell-cell signaling.

Spatially restricted microRNA directs leaf polarity through ARGONAUTE1.

Gene regulation by RNA interference requires the functions of the PAZ domain protein Argonaute. In plants, mutations in ARGONAUTE1 (AGO1) are associated with distinctive developmental defects that suggest a role for microRNA (miRNA) in organ polarity. Potential targets of miRNA regulation are the homeodomain/leucine zipper genes PHABULOSA (PHB) and PHAVOLUTA (PHV). These genes are expressed in a polar fashion in leaf primordia and are required for adaxial cell fate. Here we show that a 21-nucleotide miRNA that directs cleavage of PHB/PHV messenger RNA accumulates first in the embryonic meristem, and then in the abaxial domain of the developing leaf. miRNA distribution is disrupted by mutations in AGO1, indicating that AGO1 affects the regulation of miRNA. In addition, interactions between homeodomain/leucine zipper genes and an allelic series of ago1 indicate that miRNA acts as a signal to specify leaf polarity.

Plant stem cells: divergent pathways and common themes in shoots and roots.

Stem cells in plant shoot and root meristems are maintained throughout the life of the plant and produce somatic daughter cells that make up the body of the plant. Plant stem cells can also be derived from somatic cells in vivo and in vitro. Recent findings are refining our knowledge of signaling pathways that define stem cell fate and specify either shoot or root stem cell function. New evidence also highlights a role for epigenetic mechanisms in controlling stem cell fate.

Macro effects of microRNAs in plants.

MicroRNAs (miRNAs) are 20- to 22-nucleotide fragments that regulate expression of mRNAs that have complementary sequences. They are numerous and widespread among eukaryotes, being conserved throughout evolution. The few miRNAs that have been fully characterized were found in Caenorhabditis elegans and are required for development. Recently, a study of miRNAs isolated from Arabidopsis showed that here also developmental genes are putative regulatory targets. A role for miRNAs have in plant development is supported by the developmental phenotypes of mutations in the genes required for miRNA processing.