The phylogenetic reconstruction of the Neotropical cycad genus Ceratozamia (Zamiaceae) reveals disparate patterns of niche evolution.

The cycad genus Ceratozamia comprises 40 species from Mexico, Guatemala, Belize, and Honduras, where cycads occur throughout climatically varied montane habitats. Ceratozamia has the potential to reveal the history and processes of species diversification across diverse Neotropical habitats in this region. However, the species relationships within Ceratozamia and the ecological trends during its evolution remain unclear. Here, we aimed to clarify the phylogenetic relationships, the timing of clade and species divergences, and the niche evolution throughout the phylogenetic history of Ceratozamia. Genome-wide DNA sequences were obtained with MIG-seq, and multiple data-filtering steps were used to optimize the dataset used to construct an ultrametric species tree. Divergence times among branches and ancestral niches were estimated. The niche variation among species was evaluated, summarized into two principal components, and their ancestral states were reconstructed to test whether niche shifts among branches can be explained by random processes, under a Brownian Motion model. Ceratozamia comprises three main clades, and most species relationships within the clades were resolved. Ceratozamia has diversified since the Oligocene, with major branching events occurring during the Miocene. This timing is consistent with fossil evidence, the timing estimated for other Neotropical plant groups, and the major geological events that shaped the topographic and climatic variation in Mexico. Patterns of niche evolution in the genus do not accord with the Brownian Motion model. Rather, non-random evolution with shifts towards more seasonal environments at high latitudes, or shifts towards humid or dry environments at low latitudes explain the diversification of Ceratozamia. We present a comprehensive phylogenetic reconstruction for Ceratozamia and identify for the first time the environmental factors involved in clade and species diversification within the genus. This study alleviates the controversies regarding the species relationships in the genus and provides the first evidence that latitude-associated environmental factors may influence processes of niche evolution in cycads.

Transcriptome sequencing data provides a solid base to understand the phylogenetic relationships, biogeography and reticulated evolution of the genus Zamia L. (Cycadales: Zamiaceae).

BACKGROUND AND AIMS: Cycad is a key lineage to understand the early evolution of seed plants and their response to past environmental changes. However, tracing the evolutionary trajectory of cycad species is challenging when the robust relationships at inter- or infrageneric level are not well resolved. METHODS: Here, using 2,901 single-copy nuclear genes, we explored the species relationships and gene flow within the second largest genus of cycads, i.e., Zamia, based on phylotranscriptomic analyses of 90% extant Zamia species. Based on a well-resolved phylogenetic framework, we performed gene flow analyses, molecular dating, and biogeographical reconstruction to examine the spatiotemporal evolution of Zamia. We also performed ancestral state reconstruction (ASR) of a total of 62 traits of the genus to comprehensively investigate its morphological evolution. KEY RESULTS: Zamia is comprised of seven major clades corresponding to seven distinct distribution areas in the Americas, with at least three reticulation nodes revealed in this genus. Extant lineages of Zamia initially diversified around 18.4-32.6 (29.14) million years ago (MA) in the Mega-Mexico, and then expanded eastward into the Caribbean and southward into Central and South America. ASR revealed homoplasy in most of the morphological characters. CONCLUSIONS: This study revealed congruent phylogenetic relationships from comparative methods/datasets, with some conflicts being the result of incomplete lineage sorting and ancient/recent hybridization events. The strong association between the clades and the biogeographic areas suggested that ancient dispersal events shaped the modern distribution pattern, and regional climatic factors may have resulted in the following in-situ diversification. Climate cooling starting during the mid Miocene is associated with the global expansion of Zamia to the tropical South America that have dramatically driven lineage diversification in the New World flora, as well as the extinction of cycad species in the nowadays cooler regions of both hemispheres as indicated by the fossil records.

The evolution of red color vision is linked to coordinated rhodopsin tuning in lycaenid butterflies.

Color vision has evolved multiple times in both vertebrates and invertebrates and is largely determined by the number and variation in spectral sensitivities of distinct opsin subclasses. However, because of the difficulty of expressing long-wavelength (LW) invertebrate opsins in vitro, our understanding of the molecular basis of functional shifts in opsin spectral sensitivities has been biased toward research primarily in vertebrates. This has restricted our ability to address whether invertebrate Gq protein-coupled opsins function in a novel or convergent way compared to vertebrate Gt opsins. Here we develop a robust heterologous expression system to purify invertebrate rhodopsins, identify specific amino acid changes responsible for adaptive spectral tuning, and pinpoint how molecular variation in invertebrate opsins underlie wavelength sensitivity shifts that enhance visual perception. By combining functional and optophysiological approaches, we disentangle the relative contributions of lateral filtering pigments from red-shifted LW and blue short-wavelength opsins expressed in distinct photoreceptor cells of individual ommatidia. We use in situ hybridization to visualize six ommatidial classes in the compound eye of a lycaenid butterfly with a four-opsin visual system. We show experimentally that certain key tuning residues underlying green spectral shifts in blue opsin paralogs have evolved repeatedly among short-wavelength opsin lineages. Taken together, our results demonstrate the interplay between regulatory and adaptive evolution at multiple Gq opsin loci, as well as how coordinated spectral shifts in LW and blue opsins can act together to enhance insect spectral sensitivity at blue and red wavelengths for visual performance adaptation.

Shifting Quaternary migration patterns in the Bahamian archipelago: Evidence from the Zamia pumila complex at the northern limits of the Caribbean island biodiversity hotspot.

PREMISE OF THE STUDY: The Bahamas archipelago is formed by young, tectonically stable carbonate banks that harbor direct geological evidence of global ice-volume changes. We sought to detect signatures of major changes on gene flow patterns and reconstruct the phylogeographic history of the monophyletic Zamia pumila complex across the Bahamas. METHODS: Nuclear molecular markers with both high and low mutation rates were used to capture two different time scale signatures and test several gene flow and demographic hypotheses. KEY RESULTS: Single-copy nuclear genes unveiled apparent ancestral admixture on Andros, suggesting a significant role of this island as main hub of diversity of the archipelago. We detected demographic and spatial expansion of the Zamia pumila complex on both paleo-provinces around the Piacenzian (Pliocene)/Gelasian (Pleistocene). Populations evidenced signatures of different migration models that have occurred at two different times. Populations on Long Island (Z. lucayana) may either represent a secondary colonization of the Bahamas by Zamia or a rapid and early-divergence event of at least one population on the Bahamas. CONCLUSIONS: Despite changes in migration patterns with global climate, expected heterozygosity with both marker systems remains within the range reported for cycads, but with significant levels of increased inbreeding detected by the microsatellites. This finding is likely associated with reduced gene flow between and within paleo-provinces, accompanied by genetic drift, as rising seas enforced isolation. Our study highlights the importance of the maintenance of the predominant direction of genetic exchange and the role of overseas dispersion among the islands during climate oscillations.

Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony, likelihood and species tree inference methods.

BACKGROUND AND AIMS: Despite a recent new classification, a stable phylogeny for the cycads has been elusive, particularly regarding resolution of Bowenia, Stangeria and Dioon. In this study, five single-copy nuclear genes (SCNGs) are applied to the phylogeny of the order Cycadales. The specific aim is to evaluate several gene tree-species tree reconciliation approaches for developing an accurate phylogeny of the order, to contrast them with concatenated parsimony analysis and to resolve the erstwhile problematic phylogenetic position of these three genera. METHODS: DNA sequences of five SCNGs were obtained for 20 cycad species representing all ten genera of Cycadales. These were analysed with parsimony, maximum likelihood (ML) and three Bayesian methods of gene tree-species tree reconciliation, using Cycas as the outgroup. A calibrated date estimation was developed with Bayesian methods, and biogeographic analysis was also conducted. KEY RESULTS: Concatenated parsimony, ML and three species tree inference methods resolve exactly the same tree topology with high support at most nodes. Dioon and Bowenia are the first and second branches of Cycadales after Cycas, respectively, followed by an encephalartoid clade (Macrozamia-Lepidozamia-Encephalartos), which is sister to a zamioid clade, of which Ceratozamia is the first branch, and in which Stangeria is sister to Microcycas and Zamia. CONCLUSIONS: A single, well-supported phylogenetic hypothesis of the generic relationships of the Cycadales is presented. However, massive extinction events inferred from the fossil record that eliminated broader ancestral distributions within Zamiaceae compromise accurate optimization of ancestral biogeographical areas for that hypothesis. While major lineages of Cycadales are ancient, crown ages of all modern genera are no older than 12 million years, supporting a recent hypothesis of mostly Miocene radiations. This phylogeny can contribute to an accurate infrafamilial classification of Zamiaceae.

Cone humidity is a strong attractant in an obligate cycad pollination system.

Studies of pollination biology often focus on visual and olfactory aspects of attraction, with few studies addressing 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 furfuracea 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 produces 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.

Cycad phylogeny predicts host plant use of Eumaeus butterflies.

Eumaeus butterflies are obligate herbivores of Zamia, the most diverse neotropical genus of cycads. Eumaeus-Zamia interactions have been characterized mainly for species distributed in North and Central America. However, larval host plant use by the southern Eumaeus clade remains largely unknown, precluding a comprehensive study of co-evolution between the genera. Here, we combine fieldwork with museum and literature surveys to expand herbivory records for Eumaeus from 21 to 38 Zamia species. We inferred a time-calibrated phylogeny of Eumaeus to test for distinct macroevolutionary scenarios of larval host plant conservatism and co-evolution. We found a remarkable coincidence between Eumaeus and Zamia diversification, with the butterfly stem group diverging at the same time as the most recent radiation of Zamia in the Miocene. Cophylogenetic reconciliation analyses show a strong cophylogenetic signal between cycads and their butterfly herbivores. Bipartite model-based approaches indicate that this is because closely related Zamia species are used by the same Eumaeus species, suggesting larval host plant resource tracking by the butterfly herbivores. Our results highlight a case of tight evolution between Eumaeus butterflies and cycads, pointing to the generality of correlated evolution and phylogenetic tracking in plant-herbivore interactions across seed plants.

Zamia (Cycadales: Zamiaceae) on Puerto Rico: asymmetric genetic differentiation and the hypothesis of multiple introductions.

PREMISE OF THE STUDY: This study of Zamia in Puerto Rico is the most intensive population genetics investigation of a cycad to date in terms of number of markers, and one of few microsatellite DNA studies of plants from the highly critical Caribbean biodiversity hotspot. Three distinctive Zamia taxa occur on the island: Z. erosa on the north coast, and Z. portoricensis and Z. pumila, both in the south. Their relationships are largely unknown. We tested three hypotheses about their genetic diversity, including the possibility of multiple introductions. METHODS: We used 31 microsatellite loci across 10 populations and analyzed the data with AMOVA, Bayesian clustering, and ABC coalescent modeling. KEY RESULTS: Puerto Rican zamias exhibit an amalgam of patterns of genetic differentiation that have been reported for cycads. Overall, the taxa are slightly inbred, with high infra-populational variation and little evidence of recent bottlenecks. Zamia erosa exhibits a more than threefold greater degree of population differentiation than the other two taxa. Admixture is evident only between Z. portoricensis and Z. pumila. Zamia portoricensis is inferred to be the youngest taxon on the island, on the basis of estimates of coalescence time and effective population size. A selective sweep may be underway in a small population of Z. erosa in a saline environment. CONCLUSIONS: Zamia erosa may represent an independent introduction into Puerto Rico; Z. portoricensis and Z. pumila fit a scenario of allopatric speciation. This will be explored further in the context of genetic analysis across the entire Caribbean region.

Cycad-Weevil Pollination Symbiosis Is Characterized by Rapidly Evolving and Highly Specific Plant-Insect Chemical Communication.

Coevolution between plants and insects is thought to be responsible for generating biodiversity. Extensive research has focused largely on antagonistic herbivorous relationships, but mutualistic pollination systems also likely contribute to diversification. Here we describe an example of chemically-mediated mutualistic species interactions affecting trait evolution and lineage diversification. We show that volatile compounds produced by closely related species of Zamia cycads are more strikingly different from each other than are other phenotypic characters, and that two distantly related pollinating weevil species have specialized responses only to volatiles from their specific host Zamia species. Plant transcriptomes show that approximately a fifth of genes related to volatile production are evolving under positive selection, but we find no differences in the relative proportion of genes under positive selection in different categories. The importance of phenotypic divergence coupled with chemical communication for the maintenance of this obligate mutualism highlights chemical signaling as a key mechanism of coevolution between cycads and their weevil pollinators.

Cyclone tolerance in new world arecaceae: biogeographic variation and abiotic natural selection.

BACKGROUND AND AIMS: Consistent abiotic factors can affect directional selection; cyclones are abiotic phenomena with near-discrete geographic limits. The current study investigates selective pressure of cyclones on plants at the species level, testing for possible natural selection. METHODS: New World Arecaceae (palms) are used as a model system, as plants with monopodial, unbranched arborescent form are most directly affected by the selective pressure of wind load. Living specimens of known provenance grown at a common site were affected by the same cyclone. Data on percentage mortality were compiled and analysed in biogeographic and phylogenetic contexts. KEY RESULTS: Palms of cyclone-prone provenance exhibited a much lower (one order of magnitude) range in cyclone tolerance, and significantly lower (P < 0.001) mean percentage mortality than collections from cyclone-free areas. Palms of cyclone-free provenance had much greater variation in tolerance, and significantly greater mean percentage mortality. A test for serial independence recovered no significant phylogenetic autocorrelation of percentage mortality. CONCLUSIONS: Variation in cyclone tolerance in New World Arecaceae correlates with biogeography, and is not confounded with phylogeny. These results suggest natural selection of cyclone tolerance in cyclone-prone areas.

Conserved genetic regions across angiosperms as tools to develop single-copy nuclear markers in gymnosperms: an example using cycads.

Several individuals of the Caribbean Zamia clade and other cycad genera were used to identify single-copy nuclear genes for phylogeographic and phylogenetic studies in Cycadales. Two strategies were employed to select target loci: (i) a tblastX search of Arabidopsis conserved ortholog sequence (COS) set and (ii) a tblastX search of Arabidopsis-Populus-Vitis-Oryza Shared Single-Copy genes (APVO SSC) against the EST Zamia databases in GenBank. From the first strategy, 30 loci were selected, and from the second, 16 loci. In both cases, the matching GenBank accessions of Zamia were used as a query for retrieving highly similar sequences from Cycas, Picea, Pinus species or Ginkgo biloba. After retrieving and aligning all the sequences in each locus, intron predictions were completed to assist in primer design. PCR was carried out in three rounds to detect paralogous loci. A total of 29 loci were successfully amplified as a single band of which 20 were likely single-copy loci. These loci showed different diversity and divergence levels. A preliminary screening allowed us to select 8 promising loci (40S, ATG2, BG, GroES, GTP, LiSH, PEX4 and TR) for the Zamia pumila complex and 4 loci (COS26, GroES, GTP and HTS) for all other cycad genera.