Population structure and hybridization under contemporary and future climates in a heteroploid foundational shrub species (Artemisia tridentata).

Current and past climatic changes can shift plant climatic niches, which may cause spatial overlap or separation between related taxa. The former often leads to hybridization and introgression, which may generate novel variation and influence the adaptive capacity of plants. An additional mechanism facilitating adaptations to novel environments and an important evolutionary driver in plants is polyploidy as the result of whole genome duplication. Artemisia tridentata (big sagebrush) is a landscape-dominating foundational shrub in the western United States which occupies distinct ecological niches, exhibiting diploid and tetraploid cytotypes. Tetraploids have a large impact on the species' landscape dominance as they occupy a preponderance of the arid spectrum of A. tridentata range. Three distinct subspecies are recognized, which co-occur in ecotones - the transition zone between two or more distinct ecological niches - allowing for hybridization and introgression. Here we assess the genomic distinctiveness and extent of hybridization among subspecies at different ploidies under both contemporary and predicted future climates. We sampled five transects throughout the western United States where a subspecies overlap was predicted using subspecies-specific climate niche models. Along each transect, we sampled multiple plots representing the parental and the potential hybrid habitats. We performed reduced representation sequencing and processed the data using a ploidy-informed genotyping approach. Population genomic analyses revealed distinct diploid subspecies and at least two distinct tetraploid gene pools, indicating independent origins of the tetraploid populations. We detected low levels of hybridization (2.5%) between the diploid subspecies, while we found evidence for increased admixture between ploidy levels (18%), indicating hybridization has an important role in the formation of tetraploids. Our analyses highlight the importance of subspecies co-occurrence within these ecotones to maintain gene exchange and potential formation of tetraploid populations. Genomic confirmations of subspecies in the ecotones support the subspecies overlap predicted by the contemporary climate niche models. However, future mid-century projections of subspecies niches predict a substantial loss in range and subspecies overlap. Thus, reductions in hybridization potential could affect new recruitment of genetically variable tetraploids that are vital to this species' ecological role. Our results underscore the importance of ecotone conservation and restoration.

Growth, drought response, and climate-associated genomic structure in whitebark pine in the Sierra Nevada of California.

Whitebark pine (Pinus albicaulis Engelm.) has experienced rapid population declines and is listed as threatened under the Endangered Species Act in the United States. Whitebark pine in the Sierra Nevada of California represents the southernmost end of the species' distribution and, like other portions of its range, faces threats from an introduced pathogen, native bark beetles, and a rapidly warming climate. Beyond these chronic stressors, there is also concern about how this species will respond to acute stressors, such as drought. We present patterns of stem growth from 766 large (average diameter at breast height >25 cm), disease-free whitebark pine across the Sierra Nevada before and during a recent period of drought. We contextualize growth patterns using population genomic diversity and structure from a subset of 327 trees. Sampled whitebark pine generally had positive to neutral stem growth trends from 1970 to 2011, which was positively correlated with minimum temperature and precipitation. Indices of stem growth during drought years (2012 to 2015) relative to a predrought interval were mostly positive to neutral at our sampled sites. Individual tree growth response phenotypes appeared to be linked to genotypic variation in climate-associated loci, suggesting that some genotypes can take better advantage of local climatic conditions than others. We speculate that reduced snowpack during the 2012 to 2015 drought years may have lengthened the growing season while retaining sufficient moisture to maintain growth at most study sites. Growth responses may differ under future warming, however, particularly if drought severity increases and modifies interactions with pests and pathogens.

Subspecies differentiation and range-wide genetic structure are driven by climate in the California gnatcatcher, a flagship species for coastal sage scrub conservation.

Understanding genetic structure and diversity within species can uncover associations with environmental and geographic attributes that highlight adaptive potential and inform conservation and management. The California gnatcatcher, Polioptila californica, is a small songbird found in desert and coastal scrub habitats from the southern end of Baja California Sur to Ventura County, California. Lack of congruence among morphological subspecies hypotheses and lack of measurable genetic structure found in a few genetic markers led to questions about the validity of subspecies within P. californica and the listing status of the coastal California gnatcatcher, P. c. californica. As a U.S. federally threatened subspecies, P. c. californica is recognized as a flagship for coastal sage scrub conservation throughout southern California. We used restriction site-associated DNA sequencing to develop a genomic dataset for the California gnatcatcher. We sampled throughout the species' range, examined genetic structure, gene-environment associations, and demographic history, and tested for concordance between genetic structure and morphological subspecies groups. Our data support two distinct genetic groups with evidence of restricted movement and gene flow near the U.S.- Mexico international border. We found that climate-associated outlier loci were more strongly differentiated than climate neutral loci, suggesting that local climate adaptation may have helped to drive differentiation after Holocene range expansions. Patterns of habitat loss and fragmentation are also concordant with genetic substructure throughout the southern California portion of the range. Finally, our genetic data supported the morphologically defined P. c. californica as a distinct group, but there was little evidence of genetic differentiation among other previously hypothesized subspecies in Baja California. Our data suggest that retaining and restoring connectivity, and protecting populations, particularly at the northern range edge, could help preserve existing adaptive potential to allow for future range expansion and long-term persistence of the California gnatcatcher.

Quantitative trait loci for cell wall composition traits measured using near-infrared spectroscopy in the model C4 perennial grass Panicum hallii.

BACKGROUND: Biofuels derived from lignocellulosic plant material are an important component of current renewable energy strategies. Improvement efforts in biofuel feedstock crops have been primarily focused on increasing biomass yield with less consideration for tissue quality or composition. Four primary components found in the plant cell wall contribute to the overall quality of plant tissue and conversion characteristics, cellulose and hemicellulose polysaccharides are the primary targets for fuel conversion, while lignin and ash provide structure and defense. We explore the genetic architecture of tissue characteristics using a quantitative trait loci (QTL) mapping approach in Panicum hallii, a model lignocellulosic grass system. Diversity in the mapping population was generated by crossing xeric and mesic varietals, comparative to northern upland and southern lowland ecotypes in switchgrass. We use near-infrared spectroscopy with a primary analytical method to create a P. hallii specific calibration model to quickly quantify cell wall components. RESULTS: Ash, lignin, glucan, and xylan comprise 68% of total dry biomass in P. hallii: comparable to other feedstocks. We identified 14 QTL and one epistatic interaction across these four cell wall traits and found almost half of the QTL to localize to a single linkage group. CONCLUSIONS: Panicum hallii serves as the genomic model for its close relative and emerging biofuel crop, switchgrass (P. virgatum). We used high throughput phenotyping to map genomic regions that impact natural variation in leaf tissue composition. Understanding the genetic architecture of tissue traits in a tractable model grass system will lead to a better understanding of cell wall structure as well as provide genomic resources for bioenergy crop breeding programs.

The Genetic Basis of Upland/Lowland Ecotype Divergence in Switchgrass (Panicum virgatum).

The evolution of locally adapted ecotypes is a common phenomenon that generates diversity within plant species. However, we know surprisingly little about the genetic mechanisms underlying the locally adapted traits involved in ecotype formation. The genetic architecture underlying locally adapted traits dictates how an organism will respond to environmental selection pressures, and has major implications for evolutionary ecology, conservation, and crop breeding. To understand the genetic architecture underlying the divergence of switchgrass (Panicum virgatum) ecotypes, we constructed a genetic mapping population through a four-way outbred cross between two northern upland and two southern lowland accessions. Trait segregation in this mapping population was largely consistent with multiple independent loci controlling the suite of traits that characterizes ecotype divergence. We assembled a joint linkage map using ddRADseq, and mapped quantitative trait loci (QTL) for traits that are divergent between ecotypes, including flowering time, plant size, physiological processes, and disease resistance. Overall, we found that most QTL had small to intermediate effects. While we identified colocalizing QTL for multiple traits, we did not find any large-effect QTL that clearly controlled multiple traits through pleiotropy or tight physical linkage. These results indicate that ecologically important traits in switchgrass have a complex genetic basis, and that similar loci may underlie divergence across the geographic range of the ecotypes.

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata).

Plant-pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator-mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (FST ) for nuclear and chloroplast microsatellite loci and used a half-sib design to calculate quantitative trait divergence (QST ) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean FST to be 0.032. QST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean FST . We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non-neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.