Search for: All records

There is a general lack of consensus on the best practices for filtering of single‐nucleotide polymorphisms (SNPs) and whether it is better to use SNPs or include flanking regions (full “locus”) in phylogenomic analyses and subsequent comparative methods. Using genotyping‐by‐sequencing data from 22Glycinespecies, we assessed the effects of SNP vs. locus usage and SNP retention stringency. We compared branch length, node support, and divergence time estimation across 16 datasets with varying amounts of missing data and total size. Our results revealed five aspects of phylogenomic data usage that may be generally applicable: (1) tree topology is largely congruent across analyses; (2) filtering strictly for SNP retention (e.g., 90–100%) reduces support and can alter some inferred relationships; (3) absolute branch lengths vary by two orders of magnitude between SNP and locus datasets; (4) data type and branch length variation have little effect on divergence time estimation; and (5) phylograms alter the estimation of ancestral states and rates of morphological evolution. Using SNP or locus datasets does not alter phylogenetic inference significantly, unless researchers want or need to use absolute branch lengths. We recommend against using excessive filtering thresholds for SNP retention to reduce the risk of producing inconsistent topologies and generating low support. 

Zamia integrifolia L.f. (Cycadales), a threatened cycad native to Florida, depends on 2 native beetle species for pollination: Rhopalotria slossoni (Chevrolat; Coleoptera: Belidae) and Pharaxanotha floridana (Casey; Coleoptera: Erotylidae). Both insects are brood-site pollina- tion mutualists, known to live and feed within the pollen (male) cone. However, for pollination to occur, beetles must also visit ovulate (fe- male) cones, which have been assumed to offer no benefits to them as food or nurseries. We tested the potential for beetle pollinator use of ovulate cones by performing no-choice behavior and feeding trials for adults of both beetle species on both ovulate cones and pollen cones of Z. integrifolia. Rhopalotria slossoni beetles showed greater survival on ovulate cone tissues despite showing no significant difference in to- tal tissue mass consumed between cone sexes. Conversely, P. floridana consumed more tissue mass from ovulate cone scales yet showed no difference in survivorship on ovulate vs. pollen cone scales. Although neither beetle species is found in large numbers on ovulate cones in the field, our laboratory study suggests that both species could po- tentially benefit from feeding on ovulate cone tissues, questioning the standing hypothesis that Z. integrifolia pollination occurs by deceit. 

Abstract Identifying along which lineages shifts in diversification rates occur is a central goal of comparative phylogenetics; these shifts may coincide with key evolutionary events such as the development of novel morphological characters, the acquisition of adaptive traits, polyploidization or other structural genomic changes, or dispersal to a new habitat and subsequent increase in environmental niche space. However, while multiple methods now exist to estimate diversification rates and identify shifts using phylogenetic topologies, the appropriate use and accuracy of these methods are hotly debated. Here we test whether five Bayesian methods—Bayesian Analysis of Macroevolutionary Mixtures (BAMM), two implementations of the Lineage-Specific Birth–Death–Shift model (LSBDS and PESTO), the approximate Multi-Type Birth–Death model (MTBD; implemented in BEAST2), and the Cladogenetic Diversification Rate Shift model (ClaDS2)—produce comparable results. We apply each of these methods to a set of 65 empirical time-calibrated phylogenies and compare inferences of speciation rate, extinction rate, and net diversification rate. We find that the five methods often infer different speciation, extinction, and net-diversification rates. Consequently, these different estimates may lead to different interpretations of the macroevolutionary dynamics. The different estimates can be attributed to fundamental differences among the compared models. Therefore, the inference of shifts in diversification rates is strongly method dependent. We advise biologists to apply multiple methods to test the robustness of the conclusions or to carefully select the method based on the validity of the underlying model assumptions to their particular empirical system. 

Abstract Geological events such as mountain uplift affect how, when, and where species diversify, but measuring those effects is a longstanding challenge. Andean orogeny impacted the evolution of regional biota by creating barriers to gene flow, opening new habitats, and changing local climate. B⁢o⁢m⁢a⁢r⁢e⁢a (Alstroemeriaceae) are tropical plants with (often) small, isolated ranges; in total, B⁢o⁢m⁢a⁢r⁢e⁢a species occur from central Mexico to central Chile. This genus appears to have evolved rapidly and quite recently, and rapid radiations are often challenging to resolve with traditional phylogenetic inference. In this study, we apply phylogenomics—with hundreds of loci, gene-tree-based data curation, and a multispecies-coalescent approach—to infer the phylogeny of B⁢o⁢m⁢a⁢r⁢e⁢a. We use this phylogeny to untangle the potential drivers of diversification and biogeographic history. In particular, we test if Andean orogeny contributed to the diversification of B⁢o⁢m⁢a⁢r⁢e⁢a. We find that B⁢o⁢m⁢a⁢r⁢e⁢a originated in the central Andes during the mid-Miocene, then spread north, following the trajectory of mountain uplift. Furthermore, Andean lineages diversified faster than non-Andean relatives. B⁢o⁢m⁢a⁢r⁢e⁢a thus demonstrates that—at least in some cases—geological change rather than environmental stability has driven high species diversity in a tropical biodiversity hotspot. These results also demonstrate the utility (and danger) of genome-scale data for making macroevolutionary inferences. 

We used nuclear genomic data and statistical models to evaluate the ecological and evolutionary processes shaping spatial variation in species richness inCalochortus(Liliaceae, 74 spp.).Calochortusoccupies diverse habitats in the western United States and Mexico and has a center of diversity in the California Floristic Province, marked by multiple orogenies, winter rainfall, and highly divergent climates and substrates (including serpentine). We used sequences of 294 low-copy nuclear loci to produce a time-calibrated phylogeny, estimate historical biogeography, and test hypotheses regarding drivers of present-day spatial patterns in species number. Speciation and species coexistence require reproductive isolation and ecological divergence, so we examined the roles of chromosome number, environmental heterogeneity, and migration in shaping local species richness. Six major clades—inhabiting different geographic/climatic areas, and often marked by different base chromosome numbers (n = 6 to 10)—began diverging from each other ~10.3 Mya. As predicted, local species number increased significantly with local heterogeneity in chromosome number, elevation, soil characteristics, and serpentine presence. Species richness is greatest in the Transverse/Peninsular Ranges where clades with different chromosome numbers overlap, topographic complexity provides diverse conditions over short distances, and several physiographic provinces meet allowing immigration by several clades. Recently diverged sister-species pairs generally have peri-patric distributions, and maximum geographic overlap between species increases over the first million years since divergence, suggesting that chromosomal evolution, genetic divergence leading to gametic isolation or hybrid inviability/sterility, and/or ecological divergence over small spatial scales may permit species co-occurrence. 

Studies of pollination biology often focus on visual and olfactory aspects of attraction, with few studies ad- dressing behavioral responses and morphological adaptation to primary metabolic attributes. As part of an in-depth study of obligate nursery pollination of cycads, we find that Rhopalotria furfuracea weevils show a strong physiological response and behavioral orientation to the cone humidity of the host plant Zamia furfur- acea in an equally sensitive manner to their responses to Z. furfuracea-produced cone volatiles. Our results demonstrate that weevils can perceive fine-scale differences in relative humidity (RH) and that individuals exhibit a strong behavioral preference for higher RH in binary choice assays. Host plant Z. furfuracea pro- duces a localized cloud of higher than ambient humidity around both pollen and ovulate cones, and R. furfuracea weevils preferentially land at the zone of maximum humidity on ovulate cones, i.e., the cracks between rows of megasporophylls that provide access to the ovules. Moreover, R. furfuracea weevils exhibit striking antennal morphological traits associated with RH perception, suggesting the importance of humidity sensing in the evolution of this insect lineage. Results from this study suggest that humidity functions in a signal-like fashion in this highly specialized pollination system and help to characterize a key pollination- mediating trait in an ancient plant lineage. 

Abstract Grass leaves develop from a ring of primordial initial cells within the periphery of the shoot apical meristem, a pool of organogenic stem cells that generates all of the organs of the plant shoot. At maturity, the grass leaf is a flattened, strap-like organ comprising a proximal supportive sheath surrounding the stem and a distal photosynthetic blade. The sheath and blade are partitioned by a hinge-like auricle and the ligule, a fringe of epidermally derived tissue that grows from the adaxial (top) leaf surface. Together, the ligule and auricle comprise morphological novelties that are specific to grass leaves. Understanding how the planar outgrowth of grass leaves and their adjoining ligules is genetically controlled can yield insight into their evolutionary origins. Here we use single-cell RNA-sequencing analyses to identify a ‘rim’ cell type present at the margins of maize leaf primordia. Cells in the leaf rim have a distinctive identity and share transcriptional signatures with proliferating ligule cells, suggesting that a shared developmental genetic programme patterns both leaves and ligules. Moreover, we show that rim function is regulated by genetically redundant Wuschel-like homeobox3 (WOX3) transcription factors. Higher-order mutations in maizeWox3genes greatly reduce leaf width and disrupt ligule outgrowth and patterning. Together, these findings illustrate the generalizable use of a rim domain during planar growth of maize leaves and ligules, and suggest a parsimonious model for the homology of the grass ligule as a distal extension of the leaf sheath margin.