Extensive gene tree discordance and hemiplasy shaped the genomes of North American columnar cacti.

Few clades of plants have proven as difficult to classify as cacti. One explanation may be an unusually high level of convergent and parallel evolution (homoplasy). To evaluate support for this phylogenetic hypothesis at the molecular level, we sequenced the genomes of four cacti in the especially problematic tribe Pachycereeae, which contains most of the large columnar cacti of Mexico and adjacent areas, including the iconic saguaro cactus (Carnegiea gigantea) of the Sonoran Desert. We assembled a high-coverage draft genome for saguaro and lower coverage genomes for three other genera of tribe Pachycereeae (Pachycereus, Lophocereus, and Stenocereus) and a more distant outgroup cactus, Pereskia We used these to construct 4,436 orthologous gene alignments. Species tree inference consistently returned the same phylogeny, but gene tree discordance was high: 37% of gene trees having at least 90% bootstrap support conflicted with the species tree. Evidently, discordance is a product of long generation times and moderately large effective population sizes, leading to extensive incomplete lineage sorting (ILS). In the best supported gene trees, 58% of apparent homoplasy at amino sites in the species tree is due to gene tree-species tree discordance rather than parallel substitutions in the gene trees themselves, a phenomenon termed "hemiplasy." The high rate of genomic hemiplasy may contribute to apparent parallelisms in phenotypic traits, which could confound understanding of species relationships and character evolution in cacti.

Nonlinear modeling of saguaro growth rates reveals the importance of temperature for size-dependent growth.

PREMISE: The saguaro cactus is an iconic species of the Sonoran Desert. Its individual growth rates have been investigated for over 100 years. Its growth dynamics have been studied using phenomenological models intended to estimate growth, but not to understand the underlying biological processes. Most studies have suggested summer rainfall as the sole factor determining saguaro growth, overlooking the influence of other factors related to the process of growth. METHODS: We analyzed the annual growth rates for 13 saguaro populations in the Sonoran Desert using nonlinear models. These are better suited to analyze growth since they consider the fact that maximum growth rates diminish just before the onset of reproduction. We related model parameters to the local climate. RESULTS: The most parsimonious model was the Ricker function that described growth considering cactus decline with age. Variance in temperature, rather than precipitation, was more closely related to growth. Higher variance in temperature at the beginning of the warm season was detrimental to saguaro growth. CONCLUSIONS: Simple nonlinear equations modeled growth rate using biologically interpretable parameters related to climate factors. Because the temperature is projected to increase in both mean and variance by climate change, the population dynamics of this iconic cactus are likely to be affected.

Exceptional reduction of the plastid genome of saguaro cactus (Carnegiea gigantea): Loss of the ndh gene suite and inverted repeat.

UNLABELLED: • PREMISE OF THE STUDY: Land-plant plastid genomes have only rarely undergone significant changes in gene content and order. Thus, discovery of additional examples adds power to tests for causes of such genome-scale structural changes.• METHODS: Using next-generation sequence data, we assembled the plastid genome of saguaro cactus and probed the nuclear genome for transferred plastid genes and functionally related nuclear genes. We combined these results with available data across Cactaceae and seed plants more broadly to infer the history of gene loss and to assess the strength of phylogenetic association between gene loss and loss of the inverted repeat (IR).• KEY RESULTS: The saguaro plastid genome is the smallest known for an obligately photosynthetic angiosperm (∼113 kb), having lost the IR and plastid ndh genes. This loss supports a statistically strong association across seed plants between the loss of ndh genes and the loss of the IR. Many nonplastid copies of plastid ndh genes were found in the nuclear genome, but none had intact reading frames; nor did three related nuclear-encoded subunits. However, nuclear pgr5, which functions in a partially redundant pathway, was intact.• CONCLUSIONS: The existence of an alternative pathway redundant with the function of the plastid NADH dehydrogenase-like complex (NDH) complex may permit loss of the plastid ndh gene suite in photoautotrophs like saguaro. Loss of these genes may be a recurring mechanism for overall plastid genome size reduction, especially in combination with loss of the IR.

Are 150 km of open sea enough? Gene flow and population differentiation in a bat-pollinated columnar cactus.

Genetic differentiations and phylogeographical patterns are controlled by the interplay between spatial isolation and gene flow. To assess the extent of gene flow across an oceanic barrier, we explored the effect of the separation of the peninsula of Baja California on the evolution of mainland and peninsular populations of the long-lived columnar cactus Stenocereus thurberi. We analyzed twelve populations throughout the OPC distribution range to assess genetic diversity and structure using chloroplast DNA sequences. Genetic diversity was higher (Hd = 0.81), and genetic structure was lower (GST = 0.143) in mainland populations vs peninsular populations (Hd = 0.71, GST = 0.358 respectively). Genetic diversity was negatively associated with elevation but positively with rainfall. Two mainland and one peninsular ancestral haplotypes were reconstructed. Peninsular populations were as isolated among them as with mainland populations. Peninsular haplotypes formed a group with one mainland coastal population, and populations across the gulf shared common haplotypes giving support to regular gene flow across the Gulf. Gene flow is likely mediated by bats, the main pollinators and seed dispersers. Niche modeling suggests that during the Last Glacial Maximum (c. 130 ka), OPC populations shrank to southern locations. Currently, Stenocereus thurberi populations are expanding, and the species is under population divergence despite ongoing gene flow. Ancestral populations are located on the mainland and although vicariant peninsular populations cannot be ruled out, they are likely the result of gene flow across the seemingly formidable barrier of the Gulf of California. Still, unique haplotypes occur in the peninsula and the mainland, and peninsular populations are more structured than those on the mainland.

Effects of plant size and weather on the flowering phenology of the organ pipe cactus (Stenocereus thurberi).

BACKGROUND AND AIMS: Flowering phenology is a critical life-history trait that influences reproductive success. It has been shown that genetic, climatic and other factors such as plant size affect the timing of flowering and its duration. The spatial and temporal variation in the reproductive phenology of the columnar cactus Stenocereus thurberi and its association with plant size and environmental cues was studied. METHODS: Flowering was monitored during 3 years in three populations of S. thurberi along a latitudinal gradient. Plant size was related to phenological parameters. The actual and past weather were used for each site and year to investigate the environmental correlates of flowering. KEY RESULTS: There was significant variation in the timing of flowering within and among populations. Flowering lasted 4 months in the southern population and only 2 months in the northern population. A single flowering peak was evident in each population, but ocurred at different times. Large plants produced more flowers, and bloomed earlier and for a longer period than small plants. Population synchrony increased as the mean duration of flowering per individual decreased. The onset of flowering is primarily related to the variance in winter minimum temperatures and the duration to the autumn-winter mean maximum temperature, whereas spring mean maximum temperature is best correlated with synchrony. CONCLUSIONS: Plant size affects individual plant fecundity as well as flowering time. Thus the population structure strongly affects flowering phenology. Indications of clinal variation in the timing of flowering and reproductive effort suggest selection pressures related to the arrival of migrating pollinators, climate and resource economy in a desert environment. These pressures are likely to be relaxed in populations where individual plants can attain large sizes.

Population Genetic Structure of a Widespread Bat-Pollinated Columnar Cactus.

Bats are the main pollinators and seed dispersers of Stenocereus thurberi, a xenogamous columnar cactus of northwestern Mexico and a good model to illustrate spatial dynamics of gene flow in long-lived species. Previous studies in this cactus showed differences among populations in the type and abundance of pollinators, and in the timing of flowering and fruiting. In this study we analyzed genetic variability and population differentiation among populations. We used three primers of ISSR to analyze within and among populations genetic variation from eight widely separated populations of S. thurberi in Sonora, Mexico. Sixty-six out of 99 of the ISSR bands (P = 66.7%) were polymorphic. Total heterozygosity for all populations sampled revealed high genetic diversity (Hsp = 0.207, HBT = 0.224). The AMOVA showed that most of the genetic variation was within populations (80.5%). At the species level, estimates of population differentiation, θ = 0.175 and θB = 0.194, indicated moderate gene flow among populations. The absence of a significant correlation between genetic and geographic distances indicated little isolation by geographic distance. The large genetic variation and diversity found in S. thurberi is consistent with its open reproductive system and the high mobility of bats, a major pollinator. However, small changes in number or kind of pollinators and seed dispersal agents, in the directionality of migratory routes, and/or in the timing of flowering and fruiting among populations, can critically affect gene flow dynamics.

Geographic variation in reproductive success of Stenocereus thurberi (Cactaceae): Effects of pollination timing and pollinator guild.

PREMISE OF THE STUDY: It has been proposed that species of columnar cacti from dry tropical areas depend on bats for their reproduction, whereas species from dry subtropical areas are also pollinated by other species. To test this hypothesis, we examined the effects of pollinator guild and of variation in time and space on the reproductive success of a widespread species. • METHODS: Changes in fruit set, seed set, and pollinator activity through time were recorded in three widely separated populations of Stenocereus thurberi. Breeding system and sources of pollination limitation were determined by controlled pollinator exclusions in each population. • KEY RESULTS: Significant differences were found in the timing of activity and in the effectiveness of pollinators among sites. In the northern and central populations, reproductive success depends on bats, whereas in the southern population a combination of pollinators was more effective. No difference between open and hand cross-pollination treatments was found in the northern and central populations, which suggests no pollen limitation. However, significant differences were detected in the southern population, which indicates temporal differences in pollinator abundance or arrival time. • CONCLUSIONS: Local variation in pollinator assemblages and reproductive success could greatly affect the evolution of pollination systems. The pattern of generalist pollination in the southernmost populations and specialized pollination in the central and northern populations contradicts the hypothesis of latitudinal variation. In the absence of nocturnal pollinators, the accumulated nectar can sustain visits by diurnal pollinators, a bet-hedging strategy that increases the chances of fruit set in some populations.