The rhizosphere microbiome and host plant glucosinolates exhibit feedback cycles in Brassica rapa

The rhizosphere microbiome influences many aspects of plant fitness, including production of secondary compounds and defence against insect herbivores. Plants also modulate the composition of the microbial community in the rhizosphere via secretion of root exudates. We tested both the effect of the rhizosphere microbiome on plant traits, and host plant effects on rhizosphere microbes using recombinant inbred lines (RILs) of Brassica rapa that differ in production of glucosinolates (GLS), secondary metabolites that contribute to defence against insect herbivores. First, we investigated the effect of genetic variation in GLS production on the composition of the rhizosphere microbiome. Using a Bayesian Dirichlet-multinomial regression model (DMBVS), we identified both negative and positive associations between bacteria from six genera and the concentration of five GLS compounds produced in plant roots. Additionally, we tested the effects of microbial inoculation (an intact vs. disrupted soil microbiome) on GLS production and insect damage in these RILs. We found a significant microbial treatment × genotype interaction, in which total GLS was higher in the intact relative to the disrupted microbiome treatment in some RILs. However, despite differences in GLS production between microbial treatments, we observed no difference in insect damage between treatments. Together, these results provide evidence for a full feedback cycle of plant-microbe interactions mediated by GLS; that is, GLS compounds produced by the host plant "feed-down" to influence rhizosphere microbial community and rhizosphere microbes "feed-up" to influence GLS production.

Time outweighs the effect of host developmental stage on microbial community composition

Thousands of microbial taxa in the soil form symbioses with host plants, and due to their contribution to plant performance, these microbes are often considered an extension of the host genome. Given microbial effects on host performance, it is important to understand factors that govern microbial community assembly. Host developmental stage could affect rhizosphere microbial diversity while, alternatively, microbial assemblages could change simply as a consequence of time and the opportunity for microbial succession. Previous studies suggest that rhizosphere microbial assemblages shift across plant developmental stages, but time since germination is confounded with developmental stage. We asked how elapsed time and potential microbial succession relative to host development affected microbial diversity in the rhizosphere using monogenic flowering-time mutants of Arabidopsis thaliana. Under our experimental design, different developmental stages were present among host genotypes after the same amount of time following germination, e.g. at 76 days following germination some host genotypes were flowering while others were fruiting or senescing. We found that elapsed time was a strong predictor of microbial diversity whereas there were few differences among developmental stages. Our results support the idea that time and, likely, microbial succession more strongly affect microbial community assembly than host developmental stage.

Soil Microsite Outweighs Cultivar Genotype Contribution to Brassica Rhizobacterial Community Structure

Microorganisms residing on root surfaces play a central role in plant development and performance and may promote growth in agricultural settings. Studies have started to uncover the environmental parameters and host interactions governing their assembly. However, soil microbial communities are extremely diverse and heterogeneous, showing strong variations over short spatial scales. Here, we quantify the relative effect of meter-scale variation in soil bacterial community composition among adjacent field microsites, to better understand how microbial communities vary by host plant genotype as well as soil microsite heterogeneity. We used bacterial 16S rDNA amplicon sequencing to compare rhizosphere communities from four Brassica rapa cultivars grown in three contiguous field plots (blocks) and evaluated the relative contribution of resident soil communities and host genotypes in determining rhizosphere community structure. We characterize concomitant meter-scale variation in bacterial community structure among soils and rhizospheres and show that this block-scale variability surpasses the influence of host genotype in shaping rhizosphere communities. We identified biomarker amplicon sequence variants (ASVs) associated with bulk soil and rhizosphere habitats, each block, and three of four cultivars. Numbers and percent abundances of block-specific biomarkers in rhizosphere communities far surpassed those from bulk soils. These results highlight the importance of fine-scale variation in the pool of colonizing microorganisms during rhizosphere assembly and demonstrate that microsite variation may constitute a confounding effect while testing biotic and abiotic factors governing rhizosphere community structure.

A Nested Association Mapping Panel in Arabidopsis thaliana for Mapping and Characterizing Genetic Architecture

Linkage and association mapping populations are crucial public resources that facilitate the characterization of trait genetic architecture in natural and agricultural systems. We define a large nested association mapping panel (NAM) from 14 publicly available recombinant inbred line populations (RILs) of Arabidopsis thaliana, which share a common recurrent parent (Col-0). Using a genotype-by-sequencing approach (GBS), we identified single nucleotide polymorphisms (SNPs; range 563-1525 per population) and subsequently built updated linkage maps in each of the 14 RIL sets. Simulations in individual RIL populations indicate that our GBS markers have improved power to detect small effect QTL and enhanced resolution of QTL support intervals in comparison to original linkage maps. Using these robust linkage maps, we imputed a common set of publicly available parental SNPs into each RIL linkage map, generating overlapping markers across all populations. Though ultimately depending on allele frequencies at causal loci, simulations of the NAM panel suggest that surveying between 4 to 7 of the 14 RIL populations provides high resolution of the genetic architecture of complex traits, relative to a single mapping population.

QTL Underlying Circadian Clock Parameters Under Seasonally Variable Field Settings in Arabidopsis thaliana

The circadian clock facilitates coordination of the internal rhythms of an organism to daily environmental conditions, such as the light-dark cycle of one day. Circadian period length (the duration of one endogenous cycle) and phase (the timing of peak activity) exhibit quantitative variation in natural populations. Here, we measured circadian period and phase in June, July and September in three Arabidopsis thaliana recombinant inbred line populations. Circadian period and phase were estimated from bioluminescence of a genetic construct between a native circadian clock gene (COLD CIRCADIAN RHYTHM RNA BINDING 2) and the reporter gene (LUCIFERASE) after lines were entrained under field settings. Using a Bayesian mapping approach, we estimated the median number and effect size of genomic regions (Quantitative Trait Loci, QTL) underlying circadian parameters and the degree to which these regions overlap across months of the growing season. We also tested for QTL associations between the circadian clock and plant morphology. The genetic architecture of circadian phase was largely independent across months, as evidenced by the fact that QTL determining phase values in one month of the growing season were different from those determining phase in a second month. QTL for circadian parameters were shared with both cauline and rosette branching in at least one mapping population. The results provide insights into the QTL architecture of the clock under field settings, and suggest that the circadian clock is highly responsive to changing environments and that selection can act on clock phase in a nuanced manner.

Circadian rhythms are associated with shoot architecture in natural settings

Circadian rhythms are key regulators of diverse biological processes under controlled settings. Yet, the phenotypic and fitness consequences of quantitative variation in circadian rhythms remain largely unexplored in the field. As with other pathways, phenotypic characterization of circadian outputs in the field may reveal novel clock functions. Across consecutive growing seasons, we test for associations between clock variation and flowering phenology, plant size, shoot architecture, and fruit set in clock mutants and segregating progenies of Arabidopsis thaliana expressing quantitative variation in circadian rhythms. Using structural equation modeling, we find that genotypic variation in circadian rhythms within a growing season is associated directly with branching, which in turn affects fruit production. Consistent with direct associations between the clock and branching in segregating progenies, cauline branch number is lower and rosette branch number higher in a short-period mutant relative to wild-type and long-period genotypes, independent of flowering time. Differences in branching arise from variation in meristem fate as well as leaf production rate before flowering and attendant increases in meristem number. Our results suggest that clock variation directly affects shoot architecture in the field, suggesting a novel clock function and means by which the clock affects performance.

Bayesian estimation and use of high-throughput remote sensing indices for quantitative genetic analyses of leaf growth

We develop Bayesian function-valued trait models that mathematically isolate genetic mechanisms underlying leaf growth trajectories by factoring out genotype-specific differences in photosynthesis. Remote sensing data can be used instead of leaf-level physiological measurements. Characterizing the genetic basis of traits that vary during ontogeny and affect plant performance is a major goal in evolutionary biology and agronomy. Describing genetic programs that specifically regulate morphological traits can be complicated by genotypic differences in physiological traits. We describe the growth trajectories of leaves using novel Bayesian function-valued trait (FVT) modeling approaches in Brassica rapa recombinant inbred lines raised in heterogeneous field settings. While frequentist approaches estimate parameter values by treating each experimental replicate discretely, Bayesian models can utilize information in the global dataset, potentially leading to more robust trait estimation. We illustrate this principle by estimating growth asymptotes in the face of missing data and comparing heritabilities of growth trajectory parameters estimated by Bayesian and frequentist approaches. Using pseudo-Bayes factors, we compare the performance of an initial Bayesian logistic growth model and a model that incorporates carbon assimilation (A _max) as a cofactor, thus statistically accounting for genotypic differences in carbon resources. We further evaluate two remotely sensed spectroradiometric indices, photochemical reflectance (pri2) and MERIS Terrestrial Chlorophyll Index (mtci) as covariates in lieu of A _max, because these two indices were genetically correlated with A _max across years and treatments yet allow much higher throughput compared to direct leaf-level gas-exchange measurements. For leaf lengths in uncrowded settings, including A _max improves model fit over the initial model. The mtci and pri2 indices also outperform direct A _max measurements. Of particular importance for evolutionary biologists and plant breeders, hierarchical Bayesian models estimating FVT parameters improve heritabilities compared to frequentist approaches.

Circadian rhythms vary over the growing season and correlate with fitness components

Circadian clocks have evolved independently in all three domains of life, suggesting that internal mechanisms of time-keeping are adaptive in contemporary populations. However, the performance consequences of either discrete or quantitative clock variation have rarely been tested in field settings. Clock sensitivity of diverse segregating lines to the environment remains uncharacterized as do the statistical genetic parameters that determine evolutionary potential. In field studies with Arabidopsis thaliana, we found that major perturbations to circadian cycle length (referred to as clock period) via mutation reduce both survival and fecundity. Subtler adjustments via genomic introgression of naturally occurring alleles indicated that clock periods slightly >24 hr were adaptive, consistent with prior models describing how well the timing of biological processes is adjusted within a diurnal cycle (referred to as phase). In segregating recombinant inbred lines (RILs), circadian phase varied up to 2 hr across months of the growing season, and both period and phase expressed significant genetic variances. Performance metrics including developmental rate, size and fruit set were described by principal components (PC) analyses and circadian parameters correlated with the first PC, such that period lengths slightly >24 hr were associated with improved performance in multiple RIL sets. These experiments translate functional analyses of clock behaviour performed in controlled settings to natural ones, demonstrating that quantitative variation in circadian phase is highly responsive to seasonally variable abiotic factors. The results expand upon prior studies in controlled settings, showing that discrete and quantitative variation in clock phenotypes correlates with performance in nature.

Allocation to male vs female floral function varies by currency and responds differentially to density and moisture stress

Allocation of finite resources to separate reproductive functions is predicted to vary across environments and affect fitness. Biomass is the most commonly measured allocation currency; however, in comparison with nutrients it may be less limited and express different environmental and evolutionary responses. Here, we measured carbon, nitrogen, phosphorus, and biomass allocation among floral whorls in recombinant inbred lines of Brassica rapa in multiple environments to characterize the genetic architecture of floral allocation, including its sensitivity to environmental heterogeneity and to choice of currency. Mass, carbon, and nitrogen allocation to female whorls (pistils and sepals) decreased under high density, whereas nitrogen allocation to male organs (stamens) decreased under drought. Phosphorus allocation decreased by half in pistils under drought, while stamen phosphorus was unaffected by environment. While the contents of each currency were positively correlated among whorls, selection to improve fitness through female (or male) function typically favored increased allocation to pistils (or stamens) but decreased allocation to other whorls. Finally, genomic regions underlying correlations among allocation metrics were mapped, and loci related to nitrogen uptake and floral organ development were located within mapped quantitative trait loci. Our candidate gene identification suggests that nutrient uptake may be a limiting step in maintaining male allocation. Taken together, allocation to male vs female function is sensitive to distinct environmental stresses, and the choice of currency affects the interpretation of floral allocation responses to the environment. Further, genetic correlations may counter the evolution of allocation patterns that optimize fitness through female or male function.

FREQUENCY-DEPENDENT FITNESS IN SILENE VULGARIS, A GYNODIOECIOUS PLANT

In gynodioecious plants the selective processes that determine the relative number of female and hermaphroditic individuals are often frequency dependent. Frequency-dependent fitness can occur in the two sexes through a variety of mechanisms, especially given pollen limitation and inbreeding depression when hermaphrodites are rare. Frequency dependence in several components of the fitness of female and hermaphroditic Silene vulgaris was tested in experiments in which the relative numbers of the two sexes was varied among 12 artificial populations. In females, the proportion of flowers that set fruit covaried positively among populations with the frequency of hermaphrodites in two separate experiments, whereas the number of flowers/plant covaried negatively in one case. In hermaphrodites, the number of seeds/fruit covaried positively with the frequency of hermaphrodites, whereas the fitness of hermaphrodites estimated through pollen transfer covaried negatively. The results are discussed as they relate to the selective maintenance of gynodioecy in S. vulgaris and in light of a recent model of the effect of population structure on selection in gynodioecious systems.

Genetic architecture, biochemical underpinnings and ecological impact of floral UV patterning

Floral attraction traits can significantly affect pollinator visitation patterns, but adaptive evolution of these traits may be constrained by correlations with other traits. In some cases, molecular pathways contributing to floral attraction are well characterized, offering the opportunity to explore loci potentially underlying variation among individuals. Here, we quantify the range of variation in floral UV patterning (i.e. UV 'bulls-eye nectar guides) among crop and wild accessions of Brassica rapa. We then use experimental crosses to examine the genetic architecture, candidate loci and biochemical underpinnings of this patterning as well as phenotypic manipulations to test the ecological impact. We find qualitative variation in UV patterning between wild (commonly lacking UV patterns) and crop (commonly exhibiting UV patterns) accessions. Similar to the majority of crops, recombinant inbred lines (RILs) derived from an oilseed crop × WI fast-plant® cross exhibit UV patterns, the size of which varies extensively among genotypes. In RILs, we further observe strong statistical-genetic and QTL correlations within petal morphological traits and within measurements of petal UV patterning; however, correlations between morphology and UV patterning are weak or nonsignificant, suggesting that UV patterning is regulated and may evolve independently of overall petal size. HPLC analyses reveal a high concentration of sinapoyl glucose in UV-absorbing petal regions, which, in concert with physical locations of UV-trait QTLs, suggest a regulatory and structural gene as candidates underlying observed quantitative variation. Finally, insects prefer flowers with UV bulls-eye patterns over those that lack patterns, validating the importance of UV patterning in pollen-limited populations of B. rapa.

Modeling development and quantitative trait mapping reveal independent genetic modules for leaf size and shape

Improved predictions of fitness and yield may be obtained by characterizing the genetic controls and environmental dependencies of organismal ontogeny. Elucidating the shape of growth curves may reveal novel genetic controls that single-time-point (STP) analyses do not because, in theory, infinite numbers of growth curves can result in the same final measurement. We measured leaf lengths and widths in Brassica rapa recombinant inbred lines (RILs) throughout ontogeny. We modeled leaf growth and allometry as function valued traits (FVT), and examined genetic correlations between these traits and aspects of phenology, physiology, circadian rhythms and fitness. We used RNA-seq to construct a SNP linkage map and mapped trait quantitative trait loci (QTL). We found genetic trade-offs between leaf size and growth rate FVT and uncovered differences in genotypic and QTL correlations involving FVT vs STPs. We identified leaf shape (allometry) as a genetic module independent of length and width and identified selection on FVT parameters of development. Leaf shape is associated with venation features that affect desiccation resistance. The genetic independence of leaf shape from other leaf traits may therefore enable crop optimization in leaf shape without negative effects on traits such as size, growth rate, duration or gas exchange.

QTL architecture of reproductive fitness characters in Brassica rapa

BACKGROUND

Reproductive output is critical to both agronomists seeking to increase seed yield and to evolutionary biologists interested in understanding natural selection. We examine the genetic architecture of diverse reproductive fitness traits in recombinant inbred lines (RILs) developed from a crop (seed oil) × wild-like (rapid cycling) genotype of Brassica rapa in field and greenhouse environments.

RESULTS

Several fitness traits showed strong correlations and QTL-colocalization across environments (days to bolting, fruit length and seed color). Total fruit number was uncorrelated across environments and most QTL affecting this trait were correspondingly environment-specific. Most fitness components were positively correlated, consistent with life-history theory that genotypic variation in resource acquisition masks tradeoffs. Finally, we detected evidence of transgenerational pleiotropy, that is, maternal days to bolting was negatively correlated with days to offspring germination. A QTL for this transgenerational correlation was mapped to a genomic region harboring one copy of FLOWERING LOCUS C, a genetic locus known to affect both days to flowering as well as germination phenotypes.

CONCLUSIONS

This study characterizes the genetic structure of important fitness/yield traits within and between generations in B. rapa. Several identified QTL are suitable candidates for fine-mapping for the improvement of yield in crop Brassicas. Specifically, brFLC1, warrants further investigation as a potential regulator of phenology between generations.

Quantitative trait loci × environment interactions for plant morphology vary over ontogeny in Brassica rapa

Growth in plants occurs via the addition of repeating modules, suggesting that the genetic architecture of similar subunits may vary between earlier- and later-developing modules. These complex environment × ontogeny interactions are not well elucidated, as studies examining quantitative trait loci (QTLs) expression over ontogeny have not included multiple environments. Here, we characterized the genetic architecture of vegetative traits and onset of reproduction over ontogeny in recombinant inbred lines of Brassica rapa in the field and glasshouse. The magnitude of genetic variation in plasticity of seedling internodes was greater than in those produced later in ontogeny. We correspondingly detected that QTLs for seedling internode length were environment-specific, whereas later in ontogeny the majority of QTLs affected internode lengths in all treatments. The relationship between internode traits and onset of reproduction varied with environment and ontogenetic stage. This relationship was observed only in the glasshouse environment and was largely attributable to one environment-specific QTL. Our results provide the first evidence of a QTL × environment × ontogeny interaction, and provide QTL resolution for differences between early- and later-stage plasticity for stem elongation. These results also suggest potential constraints on morphological evolution in early vs later modules as a result of associations with reproductive timing.

Natural variation in GA1 associates with floral morphology in Arabidopsis thaliana

• The genetic architecture of floral traits is evolutionarily important due to the fitness consequences of quantitative variation in floral morphology. Yet, little is known about the genes underlying these traits in natural populations. Using Arabidopsis thaliana, we examine molecular variation at GIBBERELLIC ACID REQUIRING 1 (GA1) and test for associations with floral morphology. • We examined full-length sequence in 32 accessions and describe two haplotypes (comprising four nonsynonymous polymorphisms) in GA1 that segregate at intermediate frequencies. In 133 A. thaliana accessions, we test for genotype-phenotype associations and corroborate these findings in segregating progenies. • The two common GA1 haplotypes were associated with the length of petals, stamens, and to a lesser extent style-stigma length. Associations were confirmed in a segregating progeny developed from 19 accessions. We find analogous results in recombinant inbred lines of the Bayreuth × Shahdara cross, which differ only at one of 4 SNPs, suggesting that this SNP may contribute to the observed association. • Assuming GA1 causally affects floral organ size, it is interesting that adjacent petal and stamen whorls are most strongly affected. This pattern suggests that GA1 could contribute to the greater strength of petal-stamen correlations relative to other floral-length correlations observed in some Brassicaceous species.

Genetic architecture of the circadian clock and flowering time in Brassica rapa

The circadian clock serves to coordinate physiology and behavior with the diurnal cycles derived from the daily rotation of the earth. In plants, circadian rhythms contribute to growth and yield and, hence, to both agricultural productivity and evolutionary fitness. Arabidopsis thaliana has served as a tractable model species in which to dissect clock mechanism and function, but it now becomes important to define the extent to which the Arabidopsis model can be extrapolated to other species, including crops. Accordingly, we have extended our studies to the close Arabidopsis relative and crop species, Brassica rapa. We have investigated natural variation in circadian function and flowering time among multiple B. rapa collections. There is wide variation in clock function, based on a robust rhythm in cotyledon movement, within a collection of B. rapa accessions, wild populations and recombinant inbred lines (RILs) derived from a cross between parents from two distinct subspecies, a rapid cycling Chinese cabbage (ssp. pekinensis) and a Yellow Sarson oilseed (ssp. trilocularis). We further analyzed the RILs to identify the quantitative trait loci (QTL) responsible for this natural variation in clock period and temperature compensation, as well as for flowering time under different temperature and day length settings. Most clock and flowering-time QTL mapped to overlapping chromosomal loci. We have exploited micro-synteny between the Arabidopsis and B. rapa genomes to identify candidate genes for these QTL.

Genes underlying quantitative variation in ecologically important traits: PIF4 (phytochrome interacting factor 4) is associated with variation in internode length, flowering time, and fruit set in Arabidopsis thaliana

Association studies utilize the action of recombination over numerous generations to identify loci that underlie quantitative traits. We use a candidate-gene association approach, segregation analyses and analyses of local linkage disequilibrium (LD) to evaluate the potentially causal effects of molecular variation at PIF4 (PHYTOCROME INTERACTING FACTOR 4) on ecologically important traits in Arabidopsis thaliana. A preliminary analysis of sequence diversity in 14 natural genotypes revealed one intermediate-frequency replacement polymorphism at PIF4. A sample of 161 natural accessions was genotyped at PIF4 and screened for average length of early internodes, inflorescence length, days to flowering and flowering interval (days between bolting and flowering) under high- and low-density environments to test for genotype-phenotype associations. PIF4 was associated with early internode lengths, while the PIF4x treatment interaction was associated with flowering interval in the panel of 161 accessions. Further, in a set of recombinant inbred lines that segregate for the PIF4 polymorphism, nucleotide substitutions at PIF4 co-segregated with early internode lengths, days to flowering and fruit set, suggesting that cryptic population structure in the association-mapping panel and attendant LD with a physically distant locus do not account for the observed association. Finally, in a panel of pseudochromosomes from 20 re-sequenced genotypes, LD appeared to decay rapidly in the immediate vicinity of PIF4, suggesting that flanking loci contribute little to the observed association. In sum, the results suggest that PIF4 causally affects early internode lengths on the primary inflorescence, potentially via effects on reproductive timing and that these traits in turn affect fitness.

Genetic variation in tolerance of competition and neighbour suppression in Arabidopsis thaliana

Intraspecific competitive interactions can profoundly influence phenotypic evolution. However, prior studies have rarely evaluated the evolutionary potential of the two components of competitive ability, tolerance of competition and suppression of neighbours. Here, we grow a set of 20 Arabidopsis thaliana recombinant inbred lines in three competitive treatments (noncompetitive, intra-genotypic competition and inter-genotypic competition) to examine if there is genetic variation for the components of competitive ability and whether neighbour relatedness has an effect on fitness. We find evidence for genetic variation in tolerance of competition and neighbour suppression and that these two competitive strategies are correlated, such that genotypes that tolerate competition will also strongly suppress neighbours. We further observe that the effect of neighbour relatedness on fitness of target individuals depends on neighbour identity, i.e. whether target individuals perform better when competing against self vs. nonself individuals depends on the genotypic identity of the nonself neighbour. The results are particularly relevant to evolutionary responses under multi-level selection.

Indirect effects of FRIGIDA: floral trait (co)variances are altered by seasonally variable abiotic factors associated with flowering time

Reproductive timing is a critical life-history event that could influence the (co)variation of traits developing later in ontogeny by regulating exposure to seasonally variable factors. In a field experiment with Arabidopsis thaliana, we explore whether allelic variation at a flowering-time gene of major effect (FRIGIDA) affects (co)variation of floral traits by regulating exposure to photoperiod, temperature, and moisture levels. We detect a positive latitudinal cline in floral organ size among plants with putatively functional FRI alleles. Statistically controlling for bolting day removes the cline, suggesting that seasonal abiotic variation affects floral morphology. Both photoperiod and precipitation at bolting correlate positively with the length of petals, stamens, and pistils. Additionally, floral (co)variances differ significantly across FRI backgrounds, such that the sign of some floral-trait correlations reverses. Subsequent experimental manipulations of photoperiod and water availability demonstrate direct effects of these abiotic factors on floral traits. In sum, these results highlight how the timing of life-history events can affect the expression of traits developing later in ontogeny, and provide some of the first empirical evidence for the effects of major genes on evolutionary potential.

Polymorphic genes of major effect: consequences for variation, selection and evolution in Arabidopsis thaliana

The importance of genes of major effect for evolutionary trajectories within and among natural populations has long been the subject of intense debate. For example, if allelic variation at a major-effect locus fundamentally alters the structure of quantitative trait variation, then fixation of a single locus can have rapid and profound effects on the rate or direction of subsequent evolutionary change. Using an Arabidopsis thaliana RIL mapping population, we compare G-matrix structure between lines possessing different alleles at ERECTA, a locus known to affect ecologically relevant variation in plant architecture. We find that the allele present at ERECTA significantly alters G-matrix structure-in particular the genetic correlations between branch number and flowering time traits-and may also modulate the strength of natural selection on these traits. Despite these differences, however, when we extend our analysis to determine how evolution might differ depending on the ERECTA allele, we find that predicted responses to selection are similar. To compare responses to selection between allele classes, we developed a resampling strategy that incorporates uncertainty in estimates of selection that can also be used for statistical comparisons of G matrices.

Plasticity and environment-specific covariances: an investigation of floral-vegetative and within flower correlations

Floral traits are commonly thought to be more canalized than vegetative ones. In addition, floral and vegetative traits are hypothesized to be genetically decoupled, enabling vegetative structures to respond plastically to environmental heterogeneity, and to evolve in response to selection without disrupting the reproductive function of flowers. To test these hypotheses, we evaluate the genetic architecture of floral and vegetative traits in natural populations of Arabidopsis thaliana raised under variable light-quality environments. Plants were grown either under high or low ratios of red to far-red (R:FR) light, an aspect of light quality that varies with neighbor proximity and regulates competitive shade-avoidance responses. Across environments, we detected significant genetic variation for the average expression of all measured floral traits (petal length and width, stamen length, pistil length, stigma-anther separation, and exsertion of both the stamen and pistil beyond the corolla). Light quality significantly influenced the absolute size of several floral traits as well as the allometry (i.e., relative scaling) of all floral traits, and genotypes differed in the plasticity of floral traits to the light treatments. Exposure to low relative to high R:FR resulted in significantly greater elongation in the vegetative trait, petiole length, and genotypes again differed in the plasticity of this trait to R:FR. Consistent with prior studies, most floral traits were less plastic than the vegetative trait; herkogamy (i.e., stigma-anther separation) was the exception and expressed more variable trait values across environments than petiole length, apparently as a consequence of the independent responses of stamens and pistils. Flowers also showed strong phenotypic integration; genotypic correlations were significantly positive among floral traits within each light treatment. Although floral-vegetative correlations were not significant in the high R:FR light treatment, significant correlations were detected between petal traits, pistil length, and petiole length under low R:FR, in contrast to the widely held hypothesis that floral and vegetative traits are genetically independent. Finally, we detected selection for reduced herkogamy in the low R:FR light treatment. The observed correlation between functional trait groups suggest that vegetative plasticity may affect the expression of floral traits in some environments, and that environment-specific constraints may exist on the evolution of floral and vegetative traits.

Sequence diversity and haplotype associations with phenotypic responses to crowding: GIGANTEA affects fruit set in Arabidopsis thaliana

Identifying the molecular genetic basis of intraspecific variation in quantitative traits promises to provide novel insight into their evolutionary history as well as genetic mechanisms of adaptation. In an attempt to identify genes responsible for natural variation in competitive responses in Arabidopsis thaliana, we examined DNA sequence diversity at seven loci previously identified as members of the phytochrome B signalling network. For one gene, GIGANTEA (GI), we detected significant haplotype structure. To test for GI haplogroup-phenotype associations, we genotyped 161 A. thaliana accessions at GI and censused the same accessions for total fruit set and the expression of three phenotypic traits (days to flowering, petiole length, and inflorescence height) in a greenhouse experiment where plants were grown in crowded and uncrowded environments. We detected a significant association between GI and total fruit set that resulted in a 14% difference in average fruit set among GI haplogroups. Given that fruit set is an important component of fitness in this species and given the magnitude of the effect, the question arises as to how variation at this locus is maintained. Our observation of frequent and significant epistasis between GI and background single nucleotide polymorphisms (SNP), where the fitness ranking of the GI allele either reverses or does not differ depending on the allele at the interacting SNP, suggests that epistatic selection may actively maintain or at least slow the loss of variation at GI. This result is particularly noteworthy in the light of the ongoing debate regarding the genetic underpinnings of phenotypic evolution and recent observations that epistasis for phenotypic traits and components of fitness is common in A. thaliana.

Constraints on the evolution of adaptive plasticity: costs of plasticity to density are expressed in segregating progenies

Phenotypic plasticity, the ability of a genotype to express different phenotypes across environments, is an adaptive strategy expected to evolve in heterogeneous environments. One widely held hypothesis is that the evolutionary benefits of plasticity are reduced by its costs, but when compared with the number of traits tested, the evidence for costs is limited. Selection gradients were calculated for traits and trait plasticities to test for costs of plasticity to density in a field study using recombinant inbred lines (RILs) of Brassica rapa. Significant costs of putatively adaptive plasticity were found in three out of six measured traits. For one trait, petiole length, a cost of plasticity was detected in both environments tested; such global costs are expected to more strongly constrain the evolution of plasticity than local costs expressed in a single environment. These results, in combination with evidence from studies in segregating progenies of Arabidopsis thaliana, suggest that the potential for genetic costs of plasticity exists in natural populations. Detection of costs in previous studies may have been limited because historical selection has purged genotypes with costly plasticity, and experimental conditions often lack environmental stresses.

Resolving the genetic basis of invasiveness and predicting invasions

Considerable effort has been invested in determining traits underlying invasiveness. Yet, identifying a set of traits that commonly confers invasiveness in a range of species has proven elusive, and almost nothing is known about genetic loci affecting invasive success. Incorporating genetic model organisms into ecologically relevant studies is one promising avenue to begin dissecting the genetic underpinnings of invasiveness. Molecular biologists are rapidly characterizing genes mediating developmental responses to diverse environmental cues, i.e., genes for plasticity, as well as to environmental factors likely to impose strong selection on invading species, e.g., resistance to herbivores and competitors, coordination of life-history events with seasonal changes, and physiological tolerance of heat, drought, or cold. Here, we give an overview of molecular genetic tools increasingly used to characterize the genetic basis of adaptation and that may be used to begin identifying genetic mechanisms of invasiveness. Given the divergent traits that affect invasiveness, "invasiveness genes" common to many clades are unlikely, but the combination of developmental genetic advances with further evolutionary studies and modeling may provide a framework for identifying genes that account for invasiveness in related species.

Drought tolerance in the alpine dandelion, Taraxacum ceratophorum (Asteraceae), its exotic congener T. officinale, and interspecific hybrids under natural and experimental conditions

We compared water relations and adaptations to drought stress in native and invasive exotic dandelions, Taraxacum ceratophorum and T. officinale. Photosynthesis (A), transpiration (E), and water use efficiency (WUE; carbon gained/water lost) were measured for the two species under extreme drought in the alpine tundra of Colorado, USA. We also subjected both species and F(1) hybrids to a dry-down experiment to determine how relative physiological performance varied with water availability. Photosynthesis and transpiration in the field were low and did not differ between Taraxacum congeners; however, native T. ceratophorum had higher WUE than T. officinale. After 6 days of greenhouse drought, photosynthesis and transpiration were reduced in T. officinale compared to T. ceratophorum. Taraxacum ceratophorum maintained high WUE under control and drought treatments. Conversely, WUE in T. officinale was highly plastic between watered (low WUE) and dry-down (high WUE) treatments. Hybrids did not exhibit heterosis; instead, they were similar to T. officinale in A and E and intermediate to the parental species in WUE. Overall, results suggest that native dandelions are more drought tolerant than invasive congeners or their hybrids, but have less plasticity in WUE. Arid habitats and occasional drought in mesic sites may provide native dandelions with refugia from negative interactions with invasives.

A comparison of phenotypic plasticity in the native dandelion Taraxacum ceratophorum and its invasive congener T. officinale

We compared plastic responses to variation in the light environment for sympatric populations of native and exotic dandelion species, Taraxacum ceratophorum and Taraxacum officinale. Plasticity in leaf size, inflorescence height, reproductive phenology and dispersal-related traits were measured under experimentally altered light quality (red : far-red light ratio, R : FR) and light intensity (photosynthetically active radiation, PAR). To test whether differences in means and reaction norms of dispersal-related traits between species affected colonization potential, we created seed-dispersal models based on seed-fall rate and release height. Differences in plasticity between species were not systematic, but varied in direction and magnitude among traits. Taraxacum officinale produced larger leaves that exhibited greater plasticity in size under variable light intensity than T. ceratophorum. Plasticity in scape length at flowering occurred in relation to R : FR ratio in both species, but tended to be greater in T. ceratophorum. Seed-bearing scapes of T. officinale were taller and more canalized in height across light regimes than scapes of T. ceratophorum. Seeds of T. officinale were smaller than seeds of T. ceratophorum. Models predict greater dispersal in T. officinale within open and vegetated habitats. In contrast to the idea that plasticity promotes invasiveness, results suggest that the lack of plasticity in dispersal-related traits enhances the colonization potential of T. officinale.

The potential for genetic assimilation of a native dandelion species, Taraxacum ceratophorum (Asteraceae), by the exotic congener T. officinale

Exotic plant species can threaten closely related native congeners through asymmetric hybridization and subsequent backcrossing, the process known as genetic assimilation. I explore the initial stages of this process in Taraxacum ceratophorum (Asteraceae), the native alpine dandelion, and the invasive apomict T. officinale. In central Colorado, seven T. ceratophorum populations all occur in sympatry with T. officinale. In one large population on Pennsylvania Mountain, surveys further revealed that flowering phenologies and visiting insect taxa overlap almost completely for both Taraxacum species. Together these results indicated that heterospecific pollen transfer is likely. Crossing experiments showed that T. ceratophorum is an obligate outcrosser, and interspecific hand pollinations resulted in 37.3% seed set. However, molecular analysis of the F1 offspring indicated that only 33.2% of germinating seeds were hybrids; the remainder were selfed offspring produced from a breakdown in self-incompatibility (the mentor effect). Although the mentor effect helps reduce the production of hybrids, the asymmetrical direction of hybridization creates the potential for genetic assimilation of T. ceratophorum by T. officinale.