The genetic mechanisms underlying the convergent evolution of pollination syndromes in the Neotropical radiation of Costus L.

Selection together with variation in floral traits can act to mold floral form, often driven by a plant's predominant or most effective pollinators. To investigate the evolution of traits associated with pollination, we developed a phylogenetic framework for evaluating tempo and mode of pollination shifts across the genus Costus L., known for its evolutionary toggle between traits related to bee and bird pollination. Using a target enrichment approach, we obtained 957 loci for 171 accessions to expand the phylogenetic sampling of Neotropical Costus. In addition, we performed whole genome resequencing for a subset of 20 closely related species with contrasting pollination syndromes. For each of these 20 genomes, a high-quality assembled transcriptome was used as reference for consensus calling of candidate loci hypothesized to be associated with pollination-related traits of interest. To test for the role these candidate genes may play in evolutionary shifts in pollinators, signatures of selection were estimated as dN/dS across the identified candidate loci. We obtained a well-resolved phylogeny for Neotropical Costus despite conflict among gene trees that provide evidence of incomplete lineage sorting and/or reticulation. The overall topology and the network of genome-wide single nucleotide polymorphisms (SNPs) indicate that multiple shifts in pollination strategy have occurred across Costus, while also suggesting the presence of previously undetected signatures of hybridization between distantly related taxa. Traits related to pollination syndromes are strongly correlated and have been gained and lost in concert several times throughout the evolution of the genus. The presence of bract appendages is correlated with two traits associated with defenses against herbivory. Although labellum shape is strongly correlated with overall pollination syndrome, we found no significant impact of labellum shape on diversification rates. Evidence suggests an interplay of pollination success with other selective pressures shaping the evolution of the Costus inflorescence. Although most of the loci used for phylogenetic inference appear to be under purifying selection, many candidate genes associated with functional traits show evidence of being under positive selection. Together these results indicate an interplay of phylogenetic history with adaptive evolution leading to the diversification of pollination-associated traits in Neotropical Costus.

Expanded COVID-19 phenotype definitions reveal distinct patterns of genetic association and protective effects.

Multiple COVID-19 genome-wide association studies (GWASs) have identified reproducible genetic associations indicating that there is a genetic component to susceptibility and severity risk. To complement these studies, we collected deep coronavirus disease 2019 (COVID-19) phenotype data from a survey of 736,723 AncestryDNA research participants. With these data, we defined eight phenotypes related to COVID-19 outcomes: four phenotypes that align with previously studied COVID-19 definitions and four 'expanded' phenotypes that focus on susceptibility given exposure, mild clinical manifestations and an aggregate score of symptom severity. We performed a replication analysis of 12 previously reported COVID-19 genetic associations with all eight phenotypes in a trans-ancestry meta-analysis of AncestryDNA research participants. In this analysis, we show distinct patterns of association at the 12 loci with the eight outcomes that we assessed. We also performed a genome-wide discovery analysis of all eight phenotypes, which did not yield new genome-wide significant loci but did suggest that three of the four 'expanded' COVID-19 phenotypes have enhanced power to capture protective genetic associations relative to the previously studied phenotypes. Thus, we conclude that continued large-scale ascertainment of deep COVID-19 phenotype data would likely represent a boon for COVID-19 therapeutic target identification.

A natural diversity screen in Arabidopsis thaliana reveals determinants for HopZ1a recognition in the ZAR1-ZED1 immune complex.

Pathogen pressure on hosts can lead to the evolution of genes regulating the innate immune response. By characterizing naturally occurring polymorphisms in immune receptors, we can better understand the molecular determinants of pathogen recognition. ZAR1 is an ancient Arabidopsis thaliana NLR (Nucleotide-binding [NB] Leucine-rich-repeat [LRR] Receptor) that recognizes multiple secreted effector proteins from the pathogenic bacteria Pseudomonas syringae and Xanthomonas campestris through its interaction with receptor-like cytoplasmic kinases (RLCKs). ZAR1 was first identified for its role in recognizing P. syringae effector HopZ1a, through its interaction with the RLCK ZED1. To identify additional determinants of HopZ1a recognition, we performed a computational screen for ecotypes from the 1001 Genomes project that were likely to lack HopZ1a recognition, and tested ~300 ecotypes. We identified ecotypes containing polymorphisms in ZAR1 and ZED1. Using our previously established Nicotiana benthamiana transient assay and Arabidopsis ecotypes, we tested for the effect of naturally occurring polymorphisms on ZAR1 interactions and the immune response. We identified key residues in the NB or LRR domain of ZAR1 that impact the interaction with ZED1. We demonstrate that natural diversity combined with functional assays can help define the molecular determinants and interactions necessary to regulate immune induction in response to pathogens.

Geographic restriction, genetic divergence, and morphological disparity in the Brazilian Atlantic Forests: Insights from Leposoma lizards (Gymnophthalmidae, Squamata).

Lineage differentiation, long-term persistence, and range limitation promote high levels of phylogenetic and phylogeographic endemisms and likely underlie the abundant morphologically cryptic diversity observed in the Brazilian Atlantic Forests (AF). We explore lineage differentiation and range restriction in the AF and ask if genetic divergence and morphological disparity are correlated by integrating coalescent-based species delimitation, molecular phylogenetic, and morphological analyses in the lizard genus Leposoma. We present the first species tree for Leposoma and of their tribe, the Ecpleopodini. The analyses are based on the largest dataset ever assembled for Leposoma in terms of number of species (all represented), genetic markers (12 loci), and geographic coverage (~2,500 km). The exercise allows us to robustly delimit species within the genus and phylogeographic lineages within all species. We find support for the monophyly of the genus and for the recognition of a yet undescribed species around the Baía de Todos-os-Santos, in the state of Bahia; this form is distinct from all other congeners, both genetically and morphologically. We find that L. baturitensis, from the northeastern state of Ceará, is basal to the genus - and sister to a clade of six species restricted to the AF across the eastern coast of Brazil. Relationships within this coastal clade are ((((L. annectans, Leposoma sp.), L. scincoides), L. puk) (L. nanodactylus, L. sinepollex)). Phylogenetic and phylogeographic analyses, together with precise distribution data, allowed us to update the ranges of species and phylogeographic lineages. We reveal pervasive geographic restriction of divergent lineages in Leposoma at and below species level and discuss how forest refuges and rivers might have contributed to it. We find that morphological disparity lags behind genetic divergence in the genus because although they are correlated, the first accumulates at a much slower rate than the latter. We hope to encourage new studies in the area of AF north of the Doce river; phylogeographic sampling in that region has been much less common relative to southern sites, yet it may hold the key to several important processes defining biodiversity patterns in eastern Brazil. This appears to specially apply to processes underlying geographic restriction of morphologically cryptic, yet genetic divergent lineages, as the case of Leposoma.

Unraveling the Spiraling Radiation: A Phylogenomic Analysis of Neotropical Costus L.

The family of pantropical spiral gingers (Costaceae Nakai; c. 125 spp.) can be used as a model to enhance our understanding of the mechanisms underlying Neotropical diversity. Costaceae has higher taxonomic diversity in South and Central America (c. 72 Neotropical species, c. 30 African, c. 23 Southeast Asian), particularly due to a radiation of Neotropical species of the genus Costus L. (c. 57 spp.). However, a well-supported phylogeny of the Neotropical spiral gingers including thorough sampling of proposed species encompassing their full morphologic and geographic variation is lacking, partly due to poor resolution recovered in previous analyses using a small sampling of loci. Here we use a phylogenomic approach to estimate the phylogeny of a sample of Neotropical Costus species using a targeted enrichment approach. Baits were designed to capture conserved elements' variable at the species level using available genomic sequences of Costus species and relatives. We obtained 832 loci (generating 791,954 aligned base pairs and 31,142 parsimony informative sites) for samples that encompassed the geographical and/or morphological diversity of some recognized species. Higher support values that improve the results of previous studies were obtained when including all the available loci, even those producing unresolved gene trees and having a low proportion of variable sites. Concatenation and coalescent-based species trees methods converge in almost the same topology suggesting a robust estimation of the relationships, even under the high levels of gene tree conflict presented here. The bait set design here presented made inferring a robust phylogeny to test taxonomic hypotheses possible and will improve our understanding of the origins of the charismatic diversity of the Neotropical spiral gingers.

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).

Pleistocene diversification in an ancient lineage: a role for glacial cycles in the evolutionary history of Dioon Lindl. (Zamiaceae).

PREMISE OF THE STUDY: Recent estimates of crown ages for cycad genera (Late Miocene) challenge us to consider what processes have produced the extant diversity of this ancient group in such relatively little time. Pleistocene climate change has driven major shifts in species distributions in Mexico and may have led to speciation in the genus Dioon by forcing populations to migrate up in elevation, thereby becoming separated by topography. METHODS: We inferred orthologs from transcriptomes of five species and sequenced these in 42 individuals representing all Dioon species. From these data and published plastid sequences, we inferred dated species trees and lineage-specific diversification rates. KEY RESULTS: Analyses of 84 newly sequenced nuclear orthologs and published plastid data confirm four major clades within Dioon, all of Pleistocene age. Gene tree analysis, divergence dates, and an increase in diversification rate support very recent and rapid divergence of extant taxa. CONCLUSIONS: This study confirms the Pleistocene age of Dioon species and implicates Pleistocene climate change and established topography in lineage spitting. These results add to our understanding of the cycads as evolutionarily dynamic lineages, not relicts or evolutionary dead ends. We also find that well-supported secondary calibration points can be reliable in the absence of fossils. Our hypothesis of lineage splitting mediated by habitat shifts may be applicable to other taxa that are restricted to elevation specific ecotones.

Dense infraspecific sampling reveals rapid and independent trajectories of plastome degradation in a heterotrophic orchid complex.

Heterotrophic plants provide excellent opportunities to study the effects of altered selective regimes on genome evolution. Plastid genome (plastome) studies in heterotrophic plants are often based on one or a few highly divergent species or sequences as representatives of an entire lineage, thus missing important evolutionary-transitory events. Here, we present the first infraspecific analysis of plastome evolution in any heterotrophic plant. By combining genome skimming and targeted sequence capture, we address hypotheses on the degree and rate of plastome degradation in a complex of leafless orchids (Corallorhiza striata) across its geographic range. Plastomes provide strong support for relationships and evidence of reciprocal monophyly between C. involuta and the endangered C. bentleyi. Plastome degradation is extensive, occurring rapidly over a few million years, with evidence of differing rates of genomic change among the two principal clades of the complex. Genome skimming and targeted sequence capture differ widely in coverage depth overall, with depth in targeted sequence capture datasets varying immensely across the plastome as a function of GC content. These findings will help to fill a knowledge gap in models of heterotrophic plastid genome evolution, and have implications for future studies in heterotrophs.

Development of single-copy nuclear intron markers for species-level phylogenetics: Case study with Paullinieae (Sapindaceae).

PREMISE OF THE STUDY: We developed a bioinformatic pipeline that leverages a publicly available genome and published transcriptomes to design primers in conserved coding sequences flanking targeted introns of single-copy nuclear loci. Paullinieae (Sapindaceae) is used to demonstrate the pipeline. METHODS AND RESULTS: Transcriptome reads phylogenetically closer to the lineage of interest are aligned to the closest genome. Single-nucleotide polymorphisms are called, generating a "pseudoreference" closer to the lineage of interest. Several filters are applied to meet the criteria of single-copy nuclear loci with introns of a desired size. Primers are designed in conserved coding sequences flanking introns. Using this pipeline, we developed nine single-copy nuclear intron markers for Paullinieae. CONCLUSIONS: This pipeline is highly flexible and can be used for any group with available genomic and transcriptomic resources. This pipeline led to the development of nine variable markers for phylogenetic study without generating sequence data de novo.

Change of Fate and Staminodial Laminarity as Potential Agents of Floral Diversification in the Zingiberales.

The evolution of floral morphology in the monocot order Zingiberales shows a trend in which androecial whorl organs are progressively modified into variously conspicuous "petaloid" structures with differing degrees of fertility. Petaloidy of androecial members results from extensive laminarization of an otherwise radially symmetric structure. The genetic basis of the laminarization of androecial members has been addressed through recent candidate gene studies focused on understanding the spatiotemporal expression patterns of genes known to be necessary to floral organ formation. Here, we explore the correlation between gene duplication events and floral and inflorescence morphological diversification across the Zingiberales by inferring ancestral character states and gene copy number using the most widely accepted phylogenetic hypotheses. Our results suggest that the duplication and differential loss of GLOBOSA (GLO) copies is correlated with a change in the degree of the laminarization of androecial members. We also find an association with increased diversification in most families. We hypothesize that retention of paralogs in flower development genes could have led to a developmental shift affecting androecial organs with potential adaptive consequences, thus favoring diversification in some lineages but not others.

The phylogeny of Heliconia (Heliconiaceae) and the evolution of floral presentation.

Heliconia (Heliconiaceae, order Zingiberales) is among the showiest plants of the Neotropical rainforest and represent a spectacular co-evolutionary radiation with hummingbirds. Despite the attractiveness and ecological importance of many Heliconia, the genus has been the subject of limited molecular phylogenetic studies. We sample seven markers from the plastid and nuclear genomes for 202 samples of Heliconia. This represents ca. 75% of accepted species and includes coverage of all taxonomic subgenera and sections. We date this phylogeny using fossils associated with other families in the Zingiberales; in particular we review and evaluate the Eocene fossil Ensete oregonense. We use this dated phylogenetic framework to evaluate the evolution of two components of flower orientation that are hypothesized to be important for modulating pollinator discrimination and pollen placement: resupination and erect versus pendant inflorescence habit. Our phylogenetic results suggest that the monophyletic Melanesian subgenus Heliconiopsis and a small clade of Ecuadorian species are together the sister group to the rest of Heliconia. Extant diversity of Heliconia originated in the Late Eocene (39Ma) with rapid diversification through the Early Miocene, making it the oldest known clade of hummingbird-pollinated plants. Most described subgenera and sections are not monophyletic, though closely related groups of species, often defined by shared geography, mirror earlier morphological cladistic analyses. Evaluation of changes in resupination and inflorescence habit suggests that these characters are more homoplasious than expected, and this largely explains the non-monophyly of previously circumscribed subgenera, which were based on these characters. We also find strong evidence for the correlated evolution of resupination and inflorescence habit. The correlated model suggests that the most recent common ancestor of all extant Heliconia had resupinate flowers and erect inflorescences. Finally, we note our nearly complete species sampling and dated phylogeny allow for an assessment of taxonomic history in terms of phylogenetic diversity. We find approximately half of the currently recognized species, corresponding to half of the phylogenetic diversity, have been described since 1975, highlighting the continued importance of basic taxonomic research and conservation initiatives to preserve both described and undiscovered species of Heliconia.

Revisiting the Zingiberales: using multiplexed exon capture to resolve ancient and recent phylogenetic splits in a charismatic plant lineage.

The Zingiberales are an iconic order of monocotyledonous plants comprising eight families with distinctive and diverse floral morphologies and representing an important ecological element of tropical and subtropical forests. While the eight families are demonstrated to be monophyletic, phylogenetic relationships among these families remain unresolved. Neither combined morphological and molecular studies nor recent attempts to resolve family relationships using sequence data from whole plastomes has resulted in a well-supported, family-level phylogenetic hypothesis of relationships. Here we approach this challenge by leveraging the complete genome of one member of the order, Musa acuminata, together with transcriptome information from each of the other seven families to design a set of nuclear loci that can be enriched from highly divergent taxa with a single array-based capture of indexed genomic DNA. A total of 494 exons from 418 nuclear genes were captured for 53 ingroup taxa. The entire plastid genome was also captured for the same 53 taxa. Of the total genes captured, 308 nuclear and 68 plastid genes were used for phylogenetic estimation. The concatenated plastid and nuclear dataset supports the position of Musaceae as sister to the remaining seven families. Moreover, the combined dataset recovers known intra- and inter-family phylogenetic relationships with generally high bootstrap support. This is a flexible and cost effective method that gives the broader plant biology community a tool for generating phylogenomic scale sequence data in non-model systems at varying evolutionary depths.

Sequence capture using PCR-generated probes: a cost-effective method of targeted high-throughput sequencing for nonmodel organisms.

Recent advances in high-throughput sequencing library preparation and subgenomic enrichment methods have opened new avenues for population genetics and phylogenetics of nonmodel organisms. To multiplex large numbers of indexed samples while sequencing predominantly orthologous, targeted regions of the genome, we propose modifications to an existing, in-solution capture that utilizes PCR products as target probes to enrich library pools for the genomic subset of interest. The sequence capture using PCR-generated probes (SCPP) protocol requires no specialized equipment, is highly flexible and significantly reduces experimental costs for projects where a modest scale of genetic data is optimal (25-100 genomic loci). Our alterations enable application of this method across a wider phylogenetic range of taxa and result in higher capture efficiencies and coverage at each locus. Efficient and consistent capture over multiple SCPP experiments and at various phylogenetic distances is demonstrated, extending the utility of this method to both phylogeographic and phylogenomic studies.

Phylogenetic estimation of the core Bromelioids with an emphasis on the genus Aechmea (Bromeliaceae).

We developed a phylogeny of the core Bromelioideae including Aechmea and related genera, with the specific goals of investigating the monophyly of Aechmea and its allied genera, redefining monophyletic lineages for taxonomic revision, and investigating the biogeographic history of the group. Chloroplast, nuclear ribosomal, and low copy nuclear DNA sequences from 150 species within the Bromelioideae were used to develop the phylogeny. Phylogenies constructed with the combined four gene dataset provided sufficient resolution for investigating evolutionary relationships among species. Many genera are nested within Aechmea, or are rendered para- or polyphyletic by inclusion of Aechmea species. Several genera and subgenera of Aechmea with species in disjunct geographic locations are found to be polyphyletic, divided into separate clades that reflect geographic distribution rather than morphological similarity. This suggests that certain morphological characteristics thought to be indicative of common ancestry have instead evolved multiple times in parallel (i.e. ecological conservatism), possibly indicative of local adaptations to an epiphytic habit across the range of the Bromelioideae. These apparently homoplastic morphological characters used to assign species to genera or subgenera may be useful taxonomically when geography is also taken into account.

DNA barcoding in the cycadales: testing the potential of proposed barcoding markers for species identification of cycads.

Barcodes are short segments of DNA that can be used to uniquely identify an unknown specimen to species, particularly when diagnostic morphological features are absent. These sequences could offer a new forensic tool in plant and animal conservation-especially for endangered species such as members of the Cycadales. Ideally, barcodes could be used to positively identify illegally obtained material even in cases where diagnostic features have been purposefully removed or to release confiscated organisms into the proper breeding population. In order to be useful, a DNA barcode sequence must not only easily PCR amplify with universal or near-universal reaction conditions and primers, but also contain enough variation to generate unique identifiers at either the species or population levels. Chloroplast regions suggested by the Plant Working Group of the Consortium for the Barcode of Life (CBoL), and two alternatives, the chloroplast psbA-trnH intergenic spacer and the nuclear ribosomal internal transcribed spacer (nrITS), were tested for their utility in generating unique identifiers for members of the Cycadales. Ease of amplification and sequence generation with universal primers and reaction conditions was determined for each of the seven proposed markers. While none of the proposed markers provided unique identifiers for all species tested, nrITS showed the most promise in terms of variability, although sequencing difficulties remain a drawback. We suggest a workflow for DNA barcoding, including database generation and management, which will ultimately be necessary if we are to succeed in establishing a universal DNA barcode for plants.