Potentially toxic elements in soil-plant-water-animal continuum in a mining area from Northwestern Mexico: animal exposure pathways and health risks for children

Mining increases environmental concentrations of potentially toxic elements (PTEs) accumulating in organisms and spreading in the human food chain-their presence in milk is of great human health concern. Pathways were identified by which these elements reach raw milk from farms within a mining area in Northwestern Mexico; health risks for dairy cattle and children were also evaluated. Water from river and cattle waterers, as well as, soils showed that PTE concentrations generally below the Mexican and international limits; cattle forage concentrations were above the World Health Organization limits. Al, Mg, Mo, Ni and Zn were recorded in raw milk samples from the mining area, showing that Cd, Co, Cr, Cu, Pb and V are transferred from soil to plants but not accumulated in raw milk. Zn concentrations in raw milk exceeded the permissible limit; milk from farms without mining operations (comparison site) showed the presence of Al, Cr and Cu. In cattle tail hair, PTE did not correlate with raw milk concentrations. Metal accumulation in milk was higher through water consumption than that accumulated through forage consumption. Daily intakes (DI) of Al, Mg and Zn in cows could represent a risk for their health. The observed biotransference was higher than in other parts of Mexico, and the calculated DI and hazard quotients indicate no adverse health effects for children. However, the hazard Index values indicate that exposure to multiple PTE represents a risk for children. Management measures should be performed to control the cumulative risks to protect young children's health.

Yeast diversity in open agave fermentations across Mexico

Yeasts are a diverse group of fungal microorganisms that are widely used to produce fermented foods and beverages. In Mexico, open fermentations are used to obtain spirits from agave plants. Despite the prevalence of this traditional practice throughout the country, yeasts have only been isolated and studied from a limited number of distilleries. To systematically describe the diversity of yeast species from open agave fermentations, here we generate the YMX-1.0 culture collection by isolating 4524 strains from 68 sites with diverse climatic, geographical, and biological contexts. We used MALDI-TOF mass spectrometry for taxonomic classification and validated a subset of the strains by ITS and D1/D2 sequencing, which also revealed two potential novel species of Saccharomycetales. Overall, the composition of yeast communities was weakly associated with local variables and types of climate, yet a core set of six species was consistently isolated from most producing regions. To explore the intraspecific variation of the yeasts from agave fermentations, we sequenced the genomes of four isolates of the nonconventional yeast Kazachstania humilis. The genomes of these four strains were substantially distinct from a European isolate of the same species, suggesting that they may belong to different populations. Our work contributes to the understanding and conservation of an open fermentation system of great cultural and economic importance, providing a valuable resource to study the biology and genetic diversity of microorganisms living at the interface of natural and human-associated environments.

Water resources affected by potentially toxic elements in an area under current and historical mining in northwestern Mexico

Rio Sonora watershed and its aquifer-located in northwest Mexico-have been influenced by mining operations for 140 years, possibly causing emissions of potentially toxic elements (PTE) and affecting health of exposed populations. On the basis of available data from governmental surveys (2014-2017) and recent sampling (2018), this study constructed reliable PTE total concentration database that allowed us to report temporal/spatial variations in surface and groundwater and their associated health risks to the population living in the central part of the Rio Sonora basin. The data clearly showed that a mining spill that took place in 2014 has had an adverse impact on total PTE concentrations in surface water. They also indicated the presence of different PTE point source locations that have continued to cause contamination of surface water at levels of health concern. Data also suggested slight impacts of the spill event on groundwater possibly related to soil neutralizing potential. Two metal groups were detected for surface waters (Pb-Cd-As-Ni-Cr and of Zn-Al-Cr) and groundwaters (Cr-As-Cu-Cd and Zn-Al), which suggest that they have different sources or are being released by different processes. The potential health impacts of PTE concentrations were associated with specific age groups, dates, and areas. The results indicate that in this complex semi-arid rural system, current and historical mining activities, as well as contrasting hydrological conditions, have impacted surface and groundwater quality with important ecological and human health risks.

Barriers for plant establishment in the abandoned tailings of Nacozari, Sonora, Mexico: the influence of compost addition on seedling performance and tailing properties

Past mining activities have left a legacy of abandoned mine tailing deposits whose metal contaminants poses serious risks to ecosystems and human health. While the development of a vegetated cover in mine tailings can help in mitigating these risks, the local factors limiting plant establishment in these sites are not well understood, restricting phytostabilization efforts. Here, we explore some of the barriers that limit seedling establishment of two species (Vachellia farnesiana and Prosopis velutina) in a mine tailing deposit located in Nacozari, Sonora, Mexico, and assess whether compost addition can help in overcoming these barriers in pot and field experiments. Our field observations found 20 times more carbon and at least 4 times more nitrogen concentration in areas under vegetated patches than in non-vegetated areas, while a previous study found no difference in metal concentrations and other physicochemical parameters. This suggests that organic matter and nutrients are a major limitation for plant establishment. In agreement with this, species failed to establish without compost addition in the field experiment. Compost addition also had a positive effect on biomass accumulation, pH and microbial activity, but increased the substrate soluble concentration of As, Cu, and Zn. Nonetheless, only Cu, K, and Mo in P. velutina accumulated in tissues at levels considered toxic for animal consumption. Our study documents that compost addition facilitated plant establishment for the phytostabilization of mine tailings and help to prevent the dispersion of most metal contaminants via animal consumption. We encourage the use of complementary strategies to minimize the risk of dispersion of metal contaminants.

Niche conservatism promotes speciation in cycads: the case of Dioon merolae (Zamiaceae) in Mexico

Niche conservatism is the tendency of lineages to retain the same niche as their ancestors. It constrains biological groups and prevents ecological divergence. However, theory predicts that niche conservatism can hinder gene flow, strengthen drift and increase local adaptation: does it mean that it also can facilitate speciation? Why does this happen? We aim to answer these questions. We examined the variation of chloroplast DNA, genome-wide single nucleotide polymorphisms, morphological traits and environmental variables across the Dioon merolae cycad populations. We tested geographical structure, scenarios of demographic history, and niche conservatism between population groups. Lineage divergence is associated with the presence of a geographical barrier consisting of unsuitable habitats for cycads. There is a clear genetic and morphological distinction between the geographical groups, suggesting allopatric divergence. However, even in contrasting available environmental conditions, groups retain their ancestral niche, supporting niche conservatism. Niche conservatism is a process that can promote speciation. In D. merolae, lineage divergence occurred because unsuitable habitats represented a barrier against gene flow, incurring populations to experience isolated demographic histories and disparate environmental conditions. This study explains why cycads, despite their ancient lineage origin and biological stasis, have been able to diversify into modern ecosystems worldwide.

Mobility and accessibility of Zn, Pb, and As in abandoned mine tailings of northwestern Mexico

Generation, storage, and management of waste coming from industrial processes are a growing worldwide problem. One of the main contributors is the mining industry, in particular tailings generated by historical mining, which are barely maintained, especially in developing countries. Assessing the impact of a mining site to surrounding soils and ecosystems can be complex, especially when determining mobility and accessibility of the contaminants is required to perform ecological and human health risk assessment. As an effort to obtain information regarding mobility and accessibility of some potentially toxic elements (Zn, Pb, and As) from an historical mining site of northwestern Mexico, the abandoned mine tailings of San Felipe de Jesús in central Sonora and adjacent agricultural soils were investigated. Mobility and accessibility were assessed by means of sequential extraction procedures and using simulated physiological media. Additionally, an assessment of accidental oral intake was calculated considering the bioaccessible fractions. Results show that higher concentrations of contaminants were found in sulfide-rich tailings (Zn = 92,540; Pb = 21,288; As = 19,740 mg kg^-1) compared with oxide-rich tailings (Zn = 43,240; Pb = 14,763; As = 13,401 mg kg^-1). Concentrations in agricultural soils were on average Zn = 4755, Pb = 2840, and As = 103 mg kg^-1. Zinc was mainly recovered from labile fractions in oxide-rich tailings (~ 60%) and in a lower amount from sulfide-rich tailings (~ 30%). Pb and As were mainly associated with residual fractions (80-95%) in both types of tailings. The percentage of mobile fractions (sum of water-soluble, exchangeable, and bound to carbonate fractions) in agricultural soils was as follows: Zn ~ 60%, Pb ~ 15%, and As ~ 70%. Regarding the phytoaccessible fraction, the studied elements in mine tailings and agricultural soil samples exceeded the threshold limits, except for As in agricultural soils. According to data obtained, toxic effects were also calculated. As for daily oral intake for non-carcinogenic effects in adults and children, only Pb and As exceeded reference dose values, especially in children exposed to sulfide-rich tailings and agricultural soils. Regarding carcinogenic effects of Pb and As, most of the samples were above acceptable risk values.

Environmental assessment and historic erosion calculation of abandoned mine tailings from a semi-arid zone of northwestern Mexico: insights from geochemistry and unmanned aerial vehicles

Mining is known as one of the primary economic activities where exploitation of minerals and other materials have become essential for human development. However, this activity may represent a risk to the environment, starting from deforestation and ending with production of residues that might contain potentially toxic elements. Tailing deposits from historical mining are an example of waste that may represent an environmental concern when abandoned and exposed to environmental conditions. The town of Nacozari de Garcia, in northwestern Mexico, has three abandoned mine tailings (locally known as tailings I, II, and III) located around the urban area that represent important sources of dust and pollution. Images obtained using unmanned aerial vehicles (UAV) in conjunction with geochemical data are used to assess historic erosion calculation and pollution considering contamination and hazard indexes in tailings II and III. Digital elevation models of abandoned tailings were obtained using photogrammetry with UAV. A total of 37 surficial samples were collected from mine tailings to determine elemental concentrations (As, Cu, Pb, W, Zn) using portable X-ray fluorescence. Higher concentrations were found on samples from mine tailing II. Average concentrations followed the decreasing order of Cu > Zn > W > Pb > As for tailing II, whereas decreasing order of Cu > Zn > W > As > Pb was found for tailing III. Contamination Index (CI) values obtained from tailings II and III represent a low potential of pollution, whereas efflorescent crusts from these tailings represent a high potential of polluting soils and sediments by dust generation. Hazard Average Quotient (HAQ) values on both tailings suggest a very high potential of contamination if fluids infiltrate tailings and interact with surficial water and/or groundwater. Obtained surfaces of mine tailings II and III are 146,216 and 216,689 m², respectively, which represent around 11% of the urbanized area. A loss mass of 321,675 tons was determined for mine tailing II, whereas 634,062 tons for tailing III, accounting for 0.96 million tons of total eroded mass. Since abandonment, calculated erosion rates of 493 t ha^-1 year^-1 (tailing II) and 232 t ha^-1 year^-1 (tailing III) are in agreement with those determined in other mining areas. CI and HAQ indexes provide good estimates of pollution associated with abandoned mine tailings from Nacozari de García. Historic erosion determined in these tailings is an environmental concern since eroded material and polluted water have been incorporated into the Moctezuma River, which feeds several villages, whose major activities include agriculture and livestock raising.

Tracking Host Trees: The Phylogeography of Endemic Psittacanthus sonorae (Loranthaceae) Mistletoe in the Sonoran Desert

The host dependence of mistletoes suggests that they track the distributions of their hosts. However, the factors that determine the geographic distribution of mistletoes are not well understood. In this study, the phylogeography of Psittacanthus sonorae was reconstructed by sequencing one nuclear (ITS) and two plastid (trnL-F and atpB-rbcL) regions of 148 plants from populations separated by the Sea of Cortez. Divergence time and gene flow were estimated to gain insight into the historical demography and geographic structuring of genetic variation. We also described and mapped the spatial distribution of suitable habitat occupied by P. sonorae and its most common host Bursera microphylla in the Sonoran Desert, along with their responses to Quaternary climate fluctuations using environmental data and ecological niche modeling (ENM). We detected environmental and genetic differentiation between the peninsular and continental P. sonorae populations. Population divergence occurred during the Pleistocene, around the time of the Last Glacial Maximum. No signals of population growth were detected, with net gene flow moving from the continent to the peninsula. ENM models indicate decoupled responses by the mistletoe and its main host to past climate changes. For the Last Interglacial to the present, most models produce only partial areas of overlap on both the peninsula and the continent. Our results support a scenario of Late-Pleistocene isolation and divergence with asymmetrical gene flow between peninsular and continental P. sonorae populations. Continental populations migrated to the peninsula and the spatial isolation probably produced genetic differentiation under different environmental conditions.

The phylogeography of the cycad genus Dioon (Zamiaceae) clarifies its Cenozoic expansion and diversification in the Mexican transition zone

Background and Aims

Biogeographic transition zones are promising areas to study processes of biogeographic evolution and its influence on biological groups. The Mexican transition zone originated due to the overlap of Nearctic and Neotropical biota, which promoted great biological diversification. However, since most previous studies in this area were focused on revealing the phylogeography of Nearctic plants, how historical biogeographic configuration influenced the expansion and diversification of the Neotropical flora remains almost unknown. Using the cycad genus Dioon (Zamiaceae), this study aimed to test whether the biogeographic provinciality of the Mexican transition zone reflects the history of diversification of Neotropical plants.

Methods

Two chloroplast DNA (cpDNA) regions were analysed from 101 specimens of 15 Dioon species to reveal the distribution of haplogroups. In addition, genome-wide single nucleotide polymorphisms (SNPs) from 84 specimens were used to test the concordance between phylogenetic clusters and the biogeographic provinces. An ultrametric tree was constructed from the sequences containing SNPs to reconstruct the biogeographic events of vicariance and dispersal of Dioon across the Neotropical biogeographic provinces.

Key Results

Four Dioon lineages with strong phylogeographic structures were recognized using both cpDNA and SNP data. The lineages correspond to two clades that originated from a common ancestor in Eastern Mexico. One clade expanded and diversified in South-east Mexico and Central America. Another clade diversified into three lineages that dispersed to North-east, South and North-west Mexico. Each lineage was biogeographically delimitated. Biogeographic provinces might have provided disparate ecological conditions that facilitated speciation in Dioon since the Miocene.

Conclusions

The current genetic structure and species diversity of Dioon depict the history of expansion and diversification of the northernmost Neotropical provinces. Past biogeographic connectivities were favoured by elevated topographies, since mountain systems served as corridors for the migration of Dioon and as refugia of tropical communities that diversified during the formation of modern Neotropical forests.

Aridification as a driver of biodiversity: a case study for the cycad genus Dioon (Zamiaceae)

Background and Aims

Aridification is considered a selective pressure that might have influenced plant diversification. It is suggested that plants adapted to aridity diversified during the Miocene, an epoch of global aridification (≈15 million years ago). However, evidence supporting diversification being a direct response to aridity is scarce, and multidisciplinary evidence, besides just phylogenetic estimations, is necessary to support the idea that aridification has driven diversification. The cycad genus Dioon (Zamiaceae), a tropical group including species occurring from humid forests to arid zones, was investigated as a promising study system to understand the associations among habitat shifts, diversification times, the evolution of leaf epidermal adaptations, and aridification of Mexico.

Methods

A phylogenetic tree was constructed from seven chloroplast DNA sequences and the ITS2 spacer to reveal the relationships among 14 Dioon species from habitats ranging from humid forests to deserts. Divergence times were estimated and the habitat shifts throughout Dioon phylogeny were detected. The epidermal anatomy among Dioon species was compared and correlation tests were performed to associate the epidermal variations with habitat parameters.

Key Results

Events of habitat shifts towards arid zones happened exclusively in one of the two main clades of Dioon. Such habitat shifts happened during the species diversification of Dioon, mainly during the Miocene. Comparative anatomy showed epidermal differences between species from arid and mesic habitats. The variation of epidermal structures was found to be correlated with habitat parameters. Also, most of the analysed epidermal traits showed significant phylogenetic signals.

Conclusions

The diversification of Dioon has been driven by the aridification of Mexico. The Miocene timing corresponds to the expansion of arid zones that embedded the ancestral Dioon populations. As response, species in arid zones evolved epidermal traits to counteract aridity stress. This case study provides a robust body of evidence supporting the idea that aridification is an important driver of biodiversity.

Unique haplotypes of cacao trees as revealed by trnH-psbA chloroplast DNA

Cacao trees have been cultivated in Mesoamerica for at least 4,000 years. In this study, we analyzed sequence variation in the chloroplast DNA trnH-psbA intergenic spacer from 28 cacao trees from different farms in the Soconusco region in southern Mexico. Genetic relationships were established by two analysis approaches based on geographic origin (five populations) and genetic origin (based on a previous study). We identified six polymorphic sites, including five insertion/deletion (indels) types and one transversion. The overall nucleotide diversity was low for both approaches (geographic = 0.0032 and genetic = 0.0038). Conversely, we obtained moderate to high haplotype diversity (0.66 and 0.80) with 10 and 12 haplotypes, respectively. The common haplotype (H1) for both networks included cacao trees from all geographic locations (geographic approach) and four genetic groups (genetic approach). This common haplotype (ancient) derived a set of intermediate haplotypes and singletons interconnected by one or two mutational steps, which suggested directional selection and event purification from the expansion of narrow populations. Cacao trees from Soconusco region were grouped into one cluster without any evidence of subclustering based on AMOVA (F_ST = 0) and SAMOVA (F_ST = 0.04393) results. One population (Mazatán) showed a high haplotype frequency; thus, this population could be considered an important reservoir of genetic material. The indels located in the trnH-psbA intergenic spacer of cacao trees could be useful as markers for the development of DNA barcoding.

Genetic, morphological, geographical and ecological approaches reveal phylogenetic relationships in complex groups, an example of recently diverged pinyon pine species (Subsection Cembroides)

Elucidating phylogenetic relationships and species boundaries within complex taxonomic groups is challenging for intrinsic and extrinsic (i.e., technical) reasons. Mexican pinyon pines are a complex group whose phylogenetic relationships and species boundaries have been widely studied but poorly resolved, partly due to intrinsic ecological and evolutionary features such as low morphological and genetic differentiation caused by recent divergence, hybridization and introgression. Extrinsic factors such as limited sampling and difficulty in selecting informative molecular markers have also impeded progress. Some of the Mexican pinyon pines are of conservation concern but others may remain unprotected because the species boundaries have not been established. In this study we combined approaches to resolve the phylogenetic relationships in this complex group and to establish species boundaries in four recently diverged taxa: P. discolor, P. johannis, P. culminicola and P. cembroides. We performed phylogenetic analyses using the chloroplast markers matK and psbA-trnH as well as complete and partial chloroplast genomes of species of Subsection Cembroides. Additionally, we performed a phylogeographic analysis combining genetic data (18 chloroplast markers), morphological data and geographical data to define species boundaries in four recently diverged taxa. Ecological divergence was supported by differences in climate among localities for distinct genetic lineages. Whereas the phylogenetic analysis inferred with matK and psbA-trnH was unable to resolve the relationships in this complex group, we obtained a resolved phylogeny with the use of the chloroplast genomes. The resolved phylogeny was concordant with a haplotype network obtained using chloroplast markers. In species with potential for recent divergence, hybridization or introgression, nonhierarchical network-based approaches are probably more appropriate to protect against misclassification due to incomplete lineage sorting. The boundaries among genetic lineages were delimited by the inclusion of morphological, geographical and ecological data in the haplotype network. These multiple lines of evidence help to assign species boundaries in this complex group. P. johannis, P. discolor, P. culminicola and P. cembroides are different species based on their genetic, morphological and ecological niche differences. We suggest a reevaluation of the conservation status of these species considering the information generated in this study.

Sexual stability in the nearly dioecious Pinus johannis (Pinaceae)

PREMISE OF THE STUDY

Even though dioecy is a dominant sexual system among gymnosperms, little is known about its evolutionary history. Pinus johannis may represent a model system because unisexual and monoecious individuals compose its populations. The presence of unisexual individuals in other Pinus species is a consequence of sexual lability. Here we determined whether P. johannis represents the first example of a dioecious or nearly dioecious reproductive system in conifers by evaluating its sexual stability. •

METHODS

To assess the stability of sexual expression, we quantified the proportion of male vs. female reproductive structures produced by trees over multiple years and tested for the presence of sexual dimorphism. Sexual lability hypotheses were also examined by looking at the relationship between environmental factors and sexual expression and by comparing the reproductive behavior of P. johannis with its closest labile relative, P. edulis. •

KEY RESULTS

Pinus johannis is nearly dioecious: ~99% of individuals are unisexual or express a low proportion of the opposite gender with few changes in sexual expression through time. We found sexual dimorphism consistent with sexual stability. Sexual expression did not vary with tree size/age, abiotic environment, or herbivore removal, providing evidence against sexual lability. Individuals of P. johannis tended to produce only male or female strobili, whereas those of P. edulis were mainly monoecious with a gradient in the female to male strobili ratio. •

CONCLUSIONS

This study represents the first report of a nearly stable dioecious Pinus species. The variety of sexual morphs coexisting in the same population makes P. johannis a model for studying the evolution of dioecy in gymnosperms.

Consequences of buffelgrass pasture development for primary productivity, perennial plant richness, and vegetation structure in the drylands of Sonora, Mexico

In large parts of northern Mexico native plant communities are being converted to non-native buffelgrass (Pennisetum ciliare) pastures, and this conversion could fundamentally alter primary productivity and species richness. In Sonora, Mexico land conversion is occurring at a regional scale along a rainfall-driven gradient of primary productivity, across which native plant communities transition from desert scrub to thorn scrub. We used a paired sampling design to compare a satellite-derived index of primary productivity, richness of perennial plant species, and canopy-height profiles of native plant communities with buffelgrass pastures. We sampled species richness across a gradient of primary productivity in desert scrub and thorn scrub vegetation to examine the influence of site productivity on the outcomes of land conversion. We also examined the influence of pasture age on species richness of perennial plants. Index values of primary productivity were lower in buffelgrass pastures than in native vegetation, which suggests a reduction in primary productivity. Land conversion reduced species richness by approximately 50% at local and regional scales, reduced tree and shrub cover by 78%, and reduced canopy height. Land conversion disproportionately reduced shrub species richness, which reflects the common practice among Sonoran ranchers of conserving certain tree and cactus species. Site productivity did not affect the outcomes of land conversion. The age of a buffelgrass pasture was unrelated to species richness within the pasture, which suggests that passive recovery of species richness to preconversion levels is unlikely. Our findings demonstrate that land conversion can result in large losses of plant species richness at local and regional scales and in substantial changes to primary productivity and vegetation structure, which casts doubt on the feasibility of restoring native plant communities without active intervention on the part of land managers.

Variation in sex ratio, morph-specific reproductive ecology and an experimental test of frequency-dependence in the gynodioecious Kallstroemia grandiflora (Zygophyllaceae)

An enduring puzzle in gynodioecious species is the great variation in female frequency seen among populations. We quantified sex ratio in 44 populations of gynodioecious Kallstroemia grandiflora. Then, we measured pollinator visitation, pollen deposition, autonomous selfing rate and pollen limitation of females. Finally, using experimental populations, we tested whether female fitness responds to the frequency of female plants. We found broad variability in sex ratio among populations (0-44% female). Hermaphrodite flowers received more pollinator visits and pollen grains than females, and bagged hermaphrodite flowers produced fruits. However, we found no evidence of pollen limitation in females. In experimental populations, female plants showed no evidence of frequency-dependent pollinator visitation, fruit set, seed set or total seed mass. These results do not support frequency-dependent variation in fitness as a major mechanism affecting female frequencies in K. grandiflora. Within the context of this study, pollinators are abundant and pollinator movement appears to operate at a large enough scale to overcome the potential reproductive disadvantages of producing solely female flowers.

The genetic structure of the gynodioecious Kallstroemia grandiflora (Zygophyllaceae): the role of male sterility and colonization history

In gynodioecious populations, the frequency of females is expected to have a strong influence on the contemporary genetic structure of populations. Historical patterns of range contraction and expansion are also known to influence the genetic diversity of plant populations. We explore the influence of male sterility and colonization history on the genetic diversity in populations of Kallstroemia grandiflora along the Pacific of México. Both the overall population Fis and Fis values of hermaphrodites showed a negative association with female frequency. Genetic diversity declined with latitude. Our results provide evidence that females have a significant effect on the genetic structure as predicted by theoretical models and provide support for the hypothesis that historical processes have modified the genetic structure of K. grandiflora.

Local genetic differentiation among populations of the mass-flowering tropical shrub Erythroxylum havanense (Erythroxylaceae)

This study is aimed to understand the role of life-history traits in determining the genetic structure of populations. We used Erythroxylum havanense, a distylous shrub with synchronous and massive patterns of flowering and fruiting. We suggest that the high concentration of ephemeral resources produced by mass flowering satiates both pollinators and frugivores, restricting gene flow and leading to genetic differentiation among populations. Using random amplified polymorphic DNAs as genetic markers, we estimated genetic diversity and structure statistics to quantify the amount and distribution of genetic variation within and among five populations from the Pacific coast of Mexico. High levels of genetic variation within populations and significant differentiation among populations located very near to each other were found. Furthermore, spatial autocorrelation analysis indicated the presence of significant genetic structure at short spatial distances. We suggest that by influencing the foraging behavior of pollinators and frugivores, mass flowering may produce the observed patterns of genetic structure, while small differences in flowering or fruiting phenology could further reinforce the isolation of nearby populations.

Clonal diversity and distribution in Stenocereus eruca (Cactaceae), a narrow endemic cactus of the Sonoran Desert

Stenocereus eruca (Cactaceae), a prostrate cactus endemic to the Sonoran Desert, is thought to be highly clonal. We examined its clonal diversity and distribution: (1) at the population level, in four distinct populations along its distribution range; and (2) at a micro scale level, within a single population. Our objective was to evaluate the importance of sexual versus clonal recruitment through the use of RAPD markers. Contrary to previous field observations, clonal diversity was relatively high across the distribution range. This finding suggests that sexual recruitment is an important regeneration mechanism. The proportions of distinguishable genotypes (G/N = 0.83) and genotypic diversity (D = 0.987) were greater than in other clonal cacti, suggesting that clonal propagation is not the major regeneration mechanism. Autocorrelation analyses revealed a spatial genetic structure that may be the result of restricted gene flow (via pollen or seeds) and clonal propagation. A molecular variance analysis (AMOVA) indicated that most of the variation (66.3%) was found within and not across populations. Future studies on pollen and seed dispersal are needed to understand the role of the clonal habit in the mating system of S. eruca.

Geographic differentiation in the pollination system of the columnar cactus Pachycereus pecten -aboriginum

The pollination biology of the cactus Pachycereus pecten-aboriginum was studied in a tropical location in western Mexico (ca. 18° N latitude) to compare with data from a northern population (ca. 28° N latitude). Throughout this range, the nectar-feeding bat Leptonycteris curasoae is resident within the tropics but migratory in its northern range. The hypothesis was tested that if a predictable bat presence has been an important force in the evolution of pollination systems in columnar cacti, P. pecten-aboriginum will have a specialized pollination system within the tropics and a generalized pollination system in northern populations. In both areas, pollination experiments showed that P. pecten-aboriginum has a self-incompatible, hermaphroditic breeding system. In the tropical area, flowers open at night and close early in the morning. Nectar is secreted only during the night, and flowers are exclusively pollinated by three species of nectar-feeding bats, indicating a specialized pollination system. In contrast, anthesis and nectar secretion in northern populations occur during the night and day, allowing visitation and effective pollination by both nocturnal and diurnal pollinators. This study provides evidence of divergence mediated by pollinator predictability (resident vs. migrant), through shifts from short to long anthesis and nectar production periods from southern to northern populations.

The pollination biology of two paniculate agaves (Agavaceae) from northwestern Mexico: contrasting roles of bats as pollinators

Paniculate agaves from tropical deserts depend on nectar-feeding bats for their reproduction, while species from extratropical areas near the northern limit of Leptonycteris curasoae are pollinated by diurnal and nocturnal insects and birds. Agave angustifolia is a paniculate agave with a broad distribution in Mexico whose range coincides with the distribution of L. curasoae, while A. subsimplex has a narrow distribution in Sonora within the range of nectar-feeding bats. We studied the pollination biology of A. angustifolia and A. subsimplex in northwestern Mexico to evaluate the relative importance of bat pollination in a region where L. curasoae is seasonal. Flower visitors included a wide variety of animals, with bats having greater visitation rates in A. angustifolia. A pollinator-exclusion experiment revealed that bats were responsible for most of the control fruit set in A. angustifolia, whereas for A. subsimplex, diurnal and nocturnal visitors were equally effective. Overall, our data indicate that in central Sonora, A. angustifolia depends on nectar-feeding bats for its sexual reproductive success, while A. subsimplex relies on both diurnal and nocturnal pollinators. Given the contrasting distribution ranges of paniculate agaves, our results seem to support the trend from specialization to generalization along the continuum of tropical to extratropical deserts.

The genetic structure of a columnar cactus with a disjunct distribution: Stenocereus gummosus in the Sonoran desert

Stenocereus gummosus is a columnar cactus endemic to the Sonoran desert that exhibits a disjunct distribution: it is widely distributed in Baja California and restricted to a small coastal area in mainland Sonora. In this paper, we examine the genetic structure and the mating system of this species in order to explore the origin of the disjunction and describe aspects of the pollination biology. Flowers are nocturnal, pollinated mainly by sphingids and self-incompatible. Polymorphism for allozymes (11 loci) was relatively high (P=75%) but moderate levels of heterozygosity were detected (H(o)=0.103 and H(e)=0.261). Sonoran populations exhibited higher levels of genetic variation than peninsular populations. H(e) declined with latitude when just peninsular and two island populations are included. Substantial levels of inbreeding within populations (f=0.60), moderate differentiation among populations (theta=0.10), and no evidence of isolation by distance were detected. The neighbor-joining phenogram showed Sonoran and island populations nested within peninsular populations. Mainland populations showed greater genetic similarity to island populations, supporting a dispersal hypothesis for the origin of the disjunction. Future studies using DNA markers are suggested in order to better understand the forces that have shaped the genetic structure of this species.

Genetic differentiation in the Agave deserti (Agavaceae) complex of the Sonoran desert

The Agave deserti complex, comprising A. deserti, A. cerulata and A. subsimplex, represents a group of species and subspecies with a near allopatric distribution and clear differences in morphology. Genetic differentiation and taxonomic status with respect to spatial distribution of 14 populations of the complex were analyzed in an effort to understand the evolution and speciation process within the genus. Allelic frequencies, levels of genetic variation, expected heterozygosity (H(S)), proportion of polymorphic loci (P), and genetic differentiation (theta and Nei's genetic distance) were estimated using 41 putative RAPD loci. All three species show high levels of genetic variation (H(S)=0.12-0.29, P=63.4-95.1), and low genetic differentiation between populations and species (theta populations=0.14+/-0.02 (SE); G(st)=0.11+/-0.02). Accordingly, gene flow among populations was estimated as high by three different methods (N(m)=2.91-6.14). Nei's genetic distances between the three species were low compared to the values obtained from other Agavaceae, and there was no clear correlation with taxonomic divisions. In a UPGMA analysis, A. subsimplex and A. cerulata formed exclusive monospecific clusters, whereas the A. deserti populations appear in more than one cluster together with other species. The results were consistent with a pattern of genetic isolation by distance.

The biomechanics of Pachycereus pringlei root systems

We report on the root system of the large columnar cactus species Pachycereus pringlei to explore the hypothesis that increasing plant size decreases the ability to resist wind-throw but increases the capacity to absorb and store nutrients in roots (i.e., plant size limits the performance of these functions and may shift the performance of one function in favor of another as size increases). Based on 18 plants differing in size, the root system is characterized by a broad and deep bayonet-like root central to a shallow and extensive lateral system of root elements bearing sinker roots near the stem base. All root types have a living secondary cortex and contain wood with a large volume fraction of ray tissues that increases toward the stem base. Wood stiffness and tensile strength are correlated negatively with the ray tissue volume fraction and thus decrease toward the stem base in lateral and bayonet roots. Calculations show that the ability of the bayonet and proximal lateral root elements to resist wind-throw decreases with increasing plant size, whereas the nutrient absorption/storage capacity of the total root system increases with plant size (i.e., a size-dependent shift between these two root functions occurs).

Flower orientation in Pachycereus pringlei

The orientation of flowers in columnar cacti has been a subject of great interest to plant biologists. The interpretation of this pattern has invoked warmer temperatures as the underlying factor. In this paper, we describe flower orientation in two populations of Pachycereus pringlei (S. Watson) Britton & Rose and propose a hypothesis of the underlying mechanism. Stems from the two populations showed a significant mean direction of their flowers, with most flowers (70-77%) present between 90 and 270°. Photosynthetic photon flux density interception and stem temperature reached maximum values on south-facing ribs, showing concordance with flower orientation. We suggest that PFD interception, through its influence on CO2 uptake and stem temperature, is the major factor underlying the observed orientation of flowers in P. pringlei.Key words: columnar cacti, flower orientation, Pachycereus pringlei, PFD interception.

Wood biomechanics and anatomy of PACHYCEREUS PRINGLEI

We report the longitudinal, biomechanical, and anatomical trends observed for tissue samples drawn from the parallel aligned, prismatic woody vascular bundles running the length of a Pachycereus pringlei plant measuring 5.22 m in height. The main vertical stem of this plant was cut into five segments (labeled A through E in the acropetal direction) measuring ∼1.02 m in length. Four of the 14 vascular bundles in each segment were surgically removed to obtain 20 vascular bundle segments that were tested in bending to determine their stiffness measured in the radial E(R) and tangential E(T) direction. We also determined the lignin content of representative samples of wood.A nonlinear trend in stiffness was observed: E(R) and E(T) were highest in segments B or C (1.67 GN/m and 1.09 GN/m, respectively), lower in segment A (E(R) = 1.18 GN/m and E(T) = 0.35 GN/m), and lowest in segment E (E(R) = 0.03 GN/m and E(T) = 0.20 GN/m). Similar longitudinal trends were seen for axial tissue volume fraction and fiber wall thickness, which achieved their highest values in segment B (69.8% and 6.59 μm, respectively). Wood stiffness also correlated significantly with cell wall lignin content: with respect to segment B (which had the highest lignin content, and was thus used as the standard reference for percent lignin content), lignin content, was 15, 60, 85, and 43% in segments E, D, C, and A, respectively. Fiber cell length, which increased toward the base of the stem and toward the vascular cambium in the most proximal vascular bundle segment, did not correlate with E(R) or E(T).Basic engineering principles were used to calculate stem stresses resulting from self-loading and any wind-induced bending moment (produced by drag forces). Calculations indicated that the less stiff wood produced in segment A eliminates a rapid and potentially dangerous increase in stresses that would otherwise occur in segments B or C. The less stiff wood in segment A also reduces the probability of shear failure at the cellular interface between the wood and surrounding tissues in this portion of the stem.We conclude that P. pringlei wood stiffness is dependent on the volume fraction and lignification of axial tissues, less so on fiber wall thickness, and that wood development in this species is adaptively responsive to self-loading and differentially applied external mechanical forces.

Biomechanics of the columnar cactus Pachycereus pringlei

We report the longitudinal variations in stiffness and bulk density of tissue samples drawn from along the length of two Pachycereus pringlei plants measuring 3.69 and 5.9 m in height to determine how different tissues contribute to the mechanical stability of these massive vertical organs. Each of the two stems was cut into segments of uniform length and subsequently dissected to obtain and mechanically test portions of xylem strands, stem ribs, and a limited number of pith and cortex samples. In each case, morphometric measurements were taken to determine the geometric contribution each tissue likely made to the ability of whole stems to resist bending forces. The stiffness of each xylem strand increased basipetally toward the base of each plant where stiffness sharply decreased, reaching a magnitude comparable to that of strands 1 m beneath the stem apex. The xylem was anisotropic in behavior, i.e., its stiffness measured in the radial and in the tangential directions differed significantly. Despite the abrupt decrease in xylem strand stiffness at the stem base, the contribution made by this tissue to resist bending forces increased exponentially from the tip to the base of each plant due to the accumulation of wood. A basipetal increase in the stiffness of the pith (and, to limited extent, that of the cortex) was also observed. In contrast, the stiffness of stem rib tissues varied little as a function of stem length. These tissues were stiffer than the xylem in the corresponding portions of the stem along the upper two-fifths of the length of either plant. Tissue stiffness and bulk density were not significantly correlated within or across tissue types. However, a weak inverse relationship was observed for these properties in the case of the xylem and stem rib tissues. We present a simple formula that predicts when stem ribs rather than the xylem strands serve as the principal stiffening agents in stems. This formula successfully predicted the observed aspect ratio of the stem ribs (the average quotient of the radial and tangential dimensions of rib transections), and thus provided circumstantial evidence that the ribs are important for mechanical stability for the distal and younger regions of the stems examined.

Stem biomechanics of three columnar cacti from the Sonoran Desert

The allometric relationship of stem length L with respect to mean stem diameter D was determined for 80 shoots of each of three columnar cactus species (Stenocereus thurberi, Lophocereus schottii, and S. gummosus) to determine whether this relationship accords with that predicted by each of three contending models purporting to describe the mechanical architecture of vertical shoots (i.e., geometric, stress, and elastic similitude, which predict L proportional to D(alpha), with alpha = 1/1, 1/2, and 2/3, respectively). In addition, anatomical, physical, and biomechanical stem properties were measured to determine how the stems of these three species maintain their elastic stability as they increase in size. Reduced major axis regression of L with respect to D showed that alpha = 2.82 ± 0.14 for S. thurberi, 2.32 ± 0.19 for L. schottii, and 4.21 ± 0.31 for S. gummosus. Thus, the scaling exponents for the allometry of L differed significantly from that predicted by each of the three biomechanical models. In contrast, these exponents were similar to that for the allometry previously reported for saguaro. Analyses of biomechanical data derived from bending tests performed on 30 stems selected from each of the three species indicated that the bulk stem tissue stiffness was roughly proportional to L2, while stem flexural rigidity (i.e., the ability to resist a bending force) scaled roughly as L3. Stem length was significantly and positively correlated with the volume fraction of wood, while regression analysis of the pooled data from the three species (i.e., 90 stems) indicated that bulk tissue stiffness scaled roughly as the 5/3-power of the volume fraction of wood in stems. These data were interpreted to indicate that wood served as the major stiffening agent in stems and that this tissue accumulates at a sufficient rate to afford unusually high scaling exponents tot stem length with respect to stem diameter (i.e., disproportionately large increments of stem length with respect to increments in stem diameter). Nevertheless, the safety factor against the elastic failure of stems (computed on the basis of the critical buckling height divided by actual stem length) decreased with increasing stem size tot each species, even though each species maintained an average safety factor equal to two. We speculate that the apparent upper limit to plant height calculated for each species may serve as a biomechanical mechanism for vegetative propagation and the establishment of dense plant colonies by means of extreme stem flexure and ultimate breakage, especially for S. gummosus.

INBREEDING EFFECTS IN A GYNODIOECIOUS POPULATION OF THE COLONIZING SPECIES TRIFOLIUM HIRTUM ALL

The demographic consequences of one and two generations of selfing were examined in a gynodioecious population of rose clover (Trifolium hirtum). Seeds collected from a roadside population were used to create an experimental population that consisted of seeds of open-pollinated individuals (S₀ ), and seeds obtained after one (S₁ ) and two (S₂ ) generations of selfing. Seeds from the three groups were equally likely to germinate. However, inbreeding effects were observed in the vigor and survivorship of the seedlings. The proportion of the transplanted seedlings that reached the reproductive stage was 0.42, 0.31, and 0.14 for the S₀ , S₁ , and S₂ groups, indicating a reduction of 26% and 66% by one and two generations of selfing, respectively. Inbreeding effects on reproductive components of fitness were nonsignificant. Groups were compared by multiplicative fitness functions based on germination, survival to reproduction, and fecundity. The relative fitness of the S₁ and S₂ groups were 0.64 and 0.23 relative to the S₀ group, representing a reduction in lifetime fitness of 36% and 77%, respectively. The differences in relative fitness among the groups were caused mainly by differences in survivorship. Thus, inbreeding depression is apparently an important factor in the maintenance of the gynodioecious breeding system in rose clover.

Female frequencies and fitness components between sex phenotypes among gynodioecious populations of the colonizing species Trifolium hirtum All. in California

Male sterility has been recently discovered in Californian populations of rose clover (Trifolium hirtum). This study describes the frequency distribution of male sterility in Turkish and Californian populations, and compares fitness components of hermaphrodites and females. As male-steriles were found in Turkey, it is likely that they were introduced to California during the 1940's with the original material derived from Turkey. The spread of male-sterile genotypes in California has given rise to an asymmetrical frequency distribution of male sterility with positive skewness. The frequency of females has not exceeded fifty percent and it appears to be temporally stable in most of the Californian populations. The hypothesis that female frequencies and fitness differences between phenotypes are correlated was tested by comparing sex phenotypes in seven populations with contrasting levels of male sterility. The analysis of those populations showed no evidence for such a correlation as no significant differences were found between sex phenotypes for fecundity and seed germination. The hypothesis that females are maintained due to fitness differences in the progeny of hermaphrodites and females was experimentally tested in the population with maximum frequency of male-steriles. The results showed no significant differences in the demographic performance of the progenies of hermaphrodites and females. The present results are discussed in terms of the possible mechanism of maintenance of gynodioecy in rose clover.

Breeding systems of hermaphroditic and gynodioecious populations of the colonizing species Trifolium hirtum All. in California

A multilocus procedure was used to estimate outcrossing rates in ten roadside populations of Trifolium hirtum in California. Three groups of populations were studied: cultivars, hermaphroditic, and gynodioecious (sexually dimorphic) populations. The multilocus outcrossing rate (tm) varied from 0.05 to 0.43 among populations. Population level tm estimates were significantly correlated with the observed heterozygosity in gynodioecious populations but not in hermaphroditic populations. The outcrossing rate of hermaphrodites and females was estimated in three gynodioecious populations; the estimates of tm varied from 0.09 to 0.23 for hermaphrodites and from 0.73 to 0.80 for females. The distribution of outcrossing rates in gynodioecious populations is bimodal. Our results indicate that for the levels of selfing observed among hermaphrodites, inbreeding depression is likely to be a major factor in the maintenance of females in gynodioecious populations.