Testing the evolutionary potential of an alpine plant: Phenotypic plasticity in response to growth temperature far outweighs parental environmental effects and other genetic causes of variation.

Phenotypic plasticity and rapid evolution are fundamental processes by which organisms can maintain their function and fitness in the face of environmental changes. Here we quantified the plasticity and evolutionary potential of an alpine herb Wahlenbergia ceracea. Utilising its mixed-mating system, we generated outcrossed and self-pollinated families that were grown in either cool or warm environments, and that had parents that had also been grown in either cool or warm environments. We then analysed the contribution of environmental and genetic factors to variation in a range of phenotypic traits including phenology, leaf mass per area, photosynthetic function, thermal tolerance, and reproductive fitness. The strongest effect was that of current growth temperature, indicating strong phenotypic plasticity. All traits except thermal tolerance were plastic, whereby warm-grown plants flowered earlier, grew larger, produced more reproductive stems compared to cool-grown plants. Flowering onset and biomass were heritable and under selection, with early flowering and larger plants having higher relative fitness. There was little evidence for transgenerational plasticity, maternal effects, or genotype-by-environment interactions. Inbreeding delayed flowering and reduced reproductive fitness and biomass. Overall, we found that W. ceracea has the capacity to respond rapidly to climate warming via plasticity, and the potential for evolutionary change.

A novel and high-throughput approach to assess photosynthetic thermal tolerance of kelp using chlorophyll α fluorometry.

Foundation seaweed species are experiencing widespread declines and localized extinctions due to increased instability of sea surface temperature. Characterizing temperature thresholds are useful for predicting patterns of change and identifying species most vulnerable to extremes. Existing methods for characterizing seaweed thermal tolerance produce diverse metrics and are often time-consuming, making comparisons between species and techniques difficult, hindering insight into global patterns of change. Using three kelp species, we adapted a high-throughput method - previously used in terrestrial plant thermal biology - for use on kelps. This method employs temperature-dependent fluorescence (T-F0 ) curves under heating or cooling regimes to determine the critical temperature (Tcrit ) of photosystem II (PSII), i.e., the breakpoint between slow and fast rise fluorescence response to changing temperature, enabling rapid assays of photosynthetic thermal tolerance using a standardized metric. This method enables characterization of Tcrit for up to 48 samples per two-hour assay, demonstrating the capacity of T-F0 curves for high-throughput assays of thermal tolerance. Temperature-dependent fluorescence curves and their derived metric, Tcrit , may offer a timely and powerful new method for the field of phycology, enabling characterization and comparison of photosynthetic thermal tolerance of seaweeds across many populations, species, and biomes.

Inherent conflicts between reaction norm slope and plasticity indices when comparing plasticity: a conceptual framework and empirical test.

Phenotypic plasticity index (PI), the slope of reaction norm (K) and relative distances plasticity index (RDPI), the most commonly used estimators, have occasionally been found to generate different plasticity rankings between groups (species, populations, cultivars or genotypes). However, no effort has been made to determine how frequent this incongruence is, and the factors that influence the occurrence of the incongruence. To address these problems, we first proposed a conceptual framework and then tested the framework (its predictions) by reanalyzing 1248 sets of published data. Our framework reveals inherent conflicts between K and PI or RDPI when comparing plasticity between two groups, and the frequency of these conflicts increases with increasing inter-group initial trait difference and/or K values of the groups compared. More importantly, the estimators also affect the magnitude of the inter-group plasticity differences even when they do not change groups' plasticity rankings. The above-mentioned effects of plasticity estimators were confirmed by our empirical test using data from the literature, and the conflicts occur in 203 (16%) of the 1248 comparisons between K and indices, indicating that a considerable proportion of the comparative conclusions on plasticity in literature are estimator-dependent. The frequency of the conflicts is influenced by phylogenetic relatedness of the groups compared, being lower when comparing within relative to between species, but not by specific types of environments, traits and species. Our study indicates that care is needed to select estimator when comparing groups' plasticity, and that the conclusions in relevant literature should be treated with great caution.

Repeated extreme heatwaves result in higher leaf thermal tolerances and greater safety margins.

The frequency and severity of heatwave events are increasing, exposing species to conditions beyond their physiological limits. Species respond to heatwaves in different ways, however it remains unclear if plants have the adaptive capacity to successfully respond to hotter and more frequent heatwaves. We exposed eight tree populations from two climate regions grown under cool and warm temperatures to repeated heatwave events of moderate (40°C) and extreme (46°C) severity to assess adaptive capacity to heatwaves. Leaf damage and maximum quantum efficiency of photosystem II (Fv /Fm ) were significantly impacted by heatwave severity and growth temperatures, respectively; populations from a warm-origin avoided damage under moderate heatwaves compared to those from a cool-origin, indicating a degree of local adaptation. We found that plasticity to heatwave severity and repeated heatwaves contributed to enhanced thermal tolerance and lower leaf temperatures, leading to greater thermal safety margins (thermal tolerance minus leaf temperature) in a second heatwave. Notably, while we show that adaptation and physiological plasticity are important factors affecting plant adaptive capacity to thermal stress, plasticity of thermal tolerances and thermal safety margins provides the opportunity for trees to persist among fluctuating heatwave exposures.

The thermal tolerance of photosynthetic tissues: a global systematic review and agenda for future research.

Understanding plant thermal tolerance is fundamental to predicting impacts of extreme temperature events that are increasing in frequency and intensity across the globe. Extremes, not averages, drive species evolution, determine survival and increase crop performance. To better prioritize agricultural and natural systems research, it is crucial to evaluate how researchers are assessing the capacity of plants to tolerate extreme events. We conducted a systematic review to determine how plant thermal tolerance research is distributed across wild and domesticated plants, growth forms and biomes, and to identify crucial knowledge gaps. Our review shows that most thermal tolerance research examines cold tolerance of cultivated species; c. 5% of articles consider both heat and cold tolerance. Plants of extreme environments are understudied, and techniques widely applied in cultivated systems are largely unused in natural systems. Lastly, we find that lack of standardized methods and metrics compromises the potential for mechanistic insight. Our review provides an entry point for those new to the methods used in plant thermal tolerance research and bridges often disparate ecological and agricultural perspectives for the more experienced. We present a considered agenda of thermal tolerance research priorities to stimulate efficient, reliable and repeatable research across the spectrum of plant thermal tolerance.

Predicting species and community responses to global change using structured expert judgement: An Australian mountain ecosystems case study.

Conservation managers are under increasing pressure to make decisions about the allocation of finite resources to protect biodiversity under a changing climate. However, the impacts of climate and global change drivers on species are outpacing our capacity to collect the empirical data necessary to inform these decisions. This is particularly the case in the Australian Alps which have already undergone recent changes in climate and experienced more frequent large-scale bushfires. In lieu of empirical data, we use a structured expert elicitation method (the IDEA protocol) to estimate the change in abundance and distribution of nine vegetation groups and 89 Australian alpine and subalpine species by the year 2050. Experts predicted that most alpine vegetation communities would decline in extent by 2050; only woodlands and heathlands are predicted to increase in extent. Predicted species-level responses for alpine plants and animals were highly variable and uncertain. In general, alpine plants spanned the range of possible responses, with some expected to increase, decrease or not change in cover. By contrast, almost all animal species are predicted to decline or not change in abundance or elevation range; more species with water-centric life-cycles are expected to decline in abundance than other species. While long-term ecological data will always be the gold standard for informing the future of biodiversity, the method and outcomes outlined here provide a pragmatic and coherent basis upon which to start informing conservation policy and management in the face of rapid change and a paucity of data.

How to analyse plant phenotypic plasticity in response to a changing climate.

Contents Summary 1235 I. Introduction 1235 II. The many shapes of phenotypic plasticity 1236 III. Random regression mixed model framework 1237 IV. Conclusions 1240 Acknowledgements 1240 References 1240 SUMMARY: Plant biology is experiencing a renewed interest in the mechanistic underpinnings and evolution of phenotypic plasticity that calls for a re-evaluation of how we analyse phenotypic responses to a rapidly changing climate. We suggest that dissecting plant plasticity in response to increasing temperature needs an approach that can represent plasticity over multiple environments, and considers both population-level responses and the variation between genotypes in their response. Here, we outline how a random regression mixed model framework can be applied to plastic traits that show linear or nonlinear responses to temperature. Random regressions provide a powerful and efficient means of characterising plasticity and its variation. Although they have been used widely in other fields, they have only recently been implemented in plant evolutionary ecology. We outline their structure and provide an example tutorial of their implementation.

Genes controlling legume nodule numbers affect phenotypic plasticity responses to nitrogen in the presence and absence of rhizobia.

We investigated the role of three autoregulation of nodulation (AON) genes in regulating of root and shoot phenotypes when responding to changing nitrogen availability in the model legume, Medicago truncatula. These genes, RDN1-1 (ROOT DETERMINED NODULATION1-1), SUNN (SUPER NUMERIC NODULES), and LSS (LIKE SUNN SUPERNODULAOR), act in a systemic signalling pathway that limits nodule numbers. This pathway is also influenced by nitrogen availability, but it is not well known if AON genes control root and shoot phenotypes other than nodule numbers in response to nitrogen. We conducted a controlled glasshouse experiment to compare root and shoot phenotypes of mutants and wild type plants treated with four nitrate concentrations. All AON mutants showed altered rhizobia-independent phenotypes, including biomass allocation, lateral root length, lateral root density, and root length ratio. In response to nitrogen, uninoculated AON mutants were less plastic than the wild type in controlling root mass ratio, root length ratio, and lateral root length. This suggests that AON genes control nodulation-independent root architecture phenotypes in response to nitrogen. The phenotypic differences between wild type and AON mutants were exacerbated by the presence of nodules, pointing to resource competition as an additional mechanism affecting root and shoot responses to nitrogen.

EGRINs (Environmental Gene Regulatory Influence Networks) in Rice That Function in the Response to Water Deficit, High Temperature, and Agricultural Environments.

Environmental gene regulatory influence networks (EGRINs) coordinate the timing and rate of gene expression in response to environmental signals. EGRINs encompass many layers of regulation, which culminate in changes in accumulated transcript levels. Here, we inferred EGRINs for the response of five tropical Asian rice (Oryza sativa) cultivars to high temperatures, water deficit, and agricultural field conditions by systematically integrating time-series transcriptome data, patterns of nucleosome-free chromatin, and the occurrence of known cis-regulatory elements. First, we identified 5447 putative target genes for 445 transcription factors (TFs) by connecting TFs with genes harboring known cis-regulatory motifs in nucleosome-free regions proximal to their transcriptional start sites. We then used network component analysis to estimate the regulatory activity for each TF based on the expression of its putative target genes. Finally, we inferred an EGRIN using the estimated transcription factor activity (TFA) as the regulator. The EGRINs include regulatory interactions between 4052 target genes regulated by 113 TFs. We resolved distinct regulatory roles for members of the heat shock factor family, including a putative regulatory connection between abiotic stress and the circadian clock. TFA estimation using network component analysis is an effective way of incorporating multiple genome-scale measurements into network inference.

Temperature variability drives within-species variation in germination strategy and establishment characteristics of an alpine herb.

Plant establishment and subsequent persistence are strongly influenced by germination strategy, especially in temporally and spatially heterogeneous environments. Germination strategy determines the plant's ability to synchronise germination timing and seedling emergence to a favourable growing season and thus variation in germination strategy within species may be key to persistence under more extreme and variable future climates. However, the determinants of variation in germination strategy are not well resolved. To understand the variation of germination strategy and the climate drivers, we assessed seed traits, germination patterns, and seedling establishment traits of Oreomyrrhis eriopoda from 29 populations across its range. Germination patterns were then analysed against climate data to determine the strongest climate correlates influencing the germination strategy. Oreomyrrhis eriopoda exhibits a striking range of germination strategies among populations: varying from immediate to staggered, postponed, and postponed-deep. Seeds from regions with lower temperature variability were more likely to exhibit an immediate germination strategy; however, those patterns depended on the timescale of climatic assessment. In addition, we show that these strategy differences extend to seedling establishment traits: autumn seedlings (from populations with an immediate or staggered germination strategy) exhibited a higher leaf production rate than spring seedlings (of staggered or postponed strategy). Our results demonstrate not only substantial within-species variation in germination strategy across the species distribution range, but also that this variation correlates with environmental drivers. Given that these differences also extend to establishment traits, they may reflect a critical mechanism for persistence in changing climate.

Trees tolerate an extreme heatwave via sustained transpirational cooling and increased leaf thermal tolerance.

Heatwaves are likely to increase in frequency and intensity with climate change, which may impair tree function and forest C uptake. However, we have little information regarding the impact of extreme heatwaves on the physiological performance of large trees in the field. Here, we grew Eucalyptus parramattensis trees for 1 year with experimental warming (+3°C) in a field setting, until they were greater than 6 m tall. We withheld irrigation for 1 month to dry the surface soils and then implemented an extreme heatwave treatment of 4 consecutive days with air temperatures exceeding 43°C, while monitoring whole-canopy exchange of CO2 and H2 O, leaf temperatures, leaf thermal tolerance, and leaf and branch hydraulic status. The heatwave reduced midday canopy photosynthesis to near zero but transpiration persisted, maintaining canopy cooling. A standard photosynthetic model was unable to capture the observed decoupling between photosynthesis and transpiration at high temperatures, suggesting that climate models may underestimate a moderating feedback of vegetation on heatwave intensity. The heatwave also triggered a rapid increase in leaf thermal tolerance, such that leaf temperatures observed during the heatwave were maintained within the thermal limits of leaf function. All responses were equivalent for trees with a prior history of ambient and warmed (+3°C) temperatures, indicating that climate warming conferred no added tolerance of heatwaves expected in the future. This coordinated physiological response utilizing latent cooling and adjustment of thermal thresholds has implications for tree tolerance of future climate extremes as well as model predictions of future heatwave intensity at landscape and global scales.

Population and phylogenomic decomposition via genotyping-by-sequencing in Australian Pelargonium.

Species delimitation has seen a paradigm shift as increasing accessibility of genomic-scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.

The presence of nodules on legume root systems can alter phenotypic plasticity in response to internal nitrogen independent of nitrogen fixation.

All higher plants show developmental plasticity in response to the availability of nitrogen (N) in the soil. In legumes, N starvation causes the formation of root nodules, where symbiotic rhizobacteria fix atmospheric N2 for the host in exchange for fixed carbon (C) from the shoot. Here, we tested whether plastic responses to internal [N] of legumes are altered by their symbionts. Glasshouse experiments compared root phenotypes of three legumes, Medicago truncatula, Medicago sativa and Trifolium subterraneum, inoculated with their compatible symbiont partners and grown under four nitrate levels. In addition, six strains of rhizobia, differing in their ability to fix N2 in M. truncatula, were compared to test if plastic responses to internal [N] were dependent on the rhizobia or N2 -fixing capability of the nodules. We found that the presence of rhizobia affected phenotypic plasticity of the legumes to internal [N], particularly in root length and root mass ratio (RMR), in a plant species-dependent way. While root length responses of M. truncatula to internal [N] were dependent on the ability of rhizobial symbionts to fix N2 , RMR response to internal [N] was dependent only on initiation of nodules, irrespective of N2 -fixing ability of the rhizobia strains.

Assessing the components of adaptive capacity to improve conservation and management efforts under global change.

Natural-resource managers and other conservation practitioners are under unprecedented pressure to categorize and quantify the vulnerability of natural systems based on assessment of the exposure, sensitivity, and adaptive capacity of species to climate change. Despite the urgent need for these assessments, neither the theoretical basis of adaptive capacity nor the practical issues underlying its quantification has been articulated in a manner that is directly applicable to natural-resource management. Both are critical for researchers, managers, and other conservation practitioners to develop reliable strategies for assessing adaptive capacity. Drawing from principles of classical and contemporary research and examples from terrestrial, marine, plant, and animal systems, we examined broadly the theory behind the concept of adaptive capacity. We then considered how interdisciplinary, trait- and triage-based approaches encompassing the oft-overlooked interactions among components of adaptive capacity can be used to identify species and populations likely to have higher (or lower) adaptive capacity. We identified the challenges and value of such endeavors and argue for a concerted interdisciplinary research approach that combines ecology, ecological genetics, and eco-physiology to reflect the interacting components of adaptive capacity. We aimed to provide a basis for constructive discussion between natural-resource managers and researchers, discussions urgently needed to identify research directions that will deliver answers to real-world questions facing resource managers, other conservation practitioners, and policy makers. Directing research to both seek general patterns and identify ways to facilitate adaptive capacity of key species and populations within species, will enable conservation ecologists and resource managers to maximize returns on research and management investment and arrive at novel and dynamic management and policy decisions.

The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change.

Species are the unit of analysis in many global change and conservation biology studies; however, species are not uniform entities but are composed of different, sometimes locally adapted, populations differing in plasticity. We examined how intraspecific variation in thermal niches and phenotypic plasticity will affect species distributions in a warming climate. We first developed a conceptual model linking plasticity and niche breadth, providing five alternative intraspecific scenarios that are consistent with existing literature. Secondly, we used ecological niche-modeling techniques to quantify the impact of each intraspecific scenario on the distribution of a virtual species across a geographically realistic setting. Finally, we performed an analogous modeling exercise using real data on the climatic niches of different tree provenances. We show that when population differentiation is accounted for and dispersal is restricted, forecasts of species range shifts under climate change are even more pessimistic than those using the conventional assumption of homogeneously high plasticity across a species' range. Suitable population-level data are not available for most species so identifying general patterns of population differentiation could fill this gap. However, the literature review revealed contrasting patterns among species, urging greater levels of integration among empirical, modeling and theoretical research on intraspecific phenotypic variation.

Environmental and Biogeographic Drivers behind Alpine Plant Thermal Tolerance and Genetic Variation.

In alpine ecosystems, elevation broadly functions as a steep thermal gradient, with plant communities exposed to regular fluctuations in hot and cold temperatures. These conditions lead to selective filtering, potentially contributing to species-level variation in thermal tolerance and population-level genetic divergence. Few studies have explored the breadth of alpine plant thermal tolerances across a thermal gradient or the underlying genetic variation thereof. We measured photosystem heat (Tcrit-hot) and cold (Tcrit-cold) thresholds of ten Australian alpine species across elevation gradients and characterised their neutral genetic variation. To reveal the biogeographical drivers of present-day genetic signatures, we also reconstructed temporal changes in habitat suitability across potential distributional ranges. We found intraspecific variation in thermal thresholds, but this was not associated with elevation, nor underpinned by genetic differentiation on a local scale. Instead, regional population differentiation and considerable homozygosity within populations may, in part, be driven by distributional contractions, long-term persistence, and migrations following habitat suitability. Our habitat suitability models suggest that cool-climate-distributed alpine plants may be threatened by a warming climate. Yet, the observed wide thermal tolerances did not reflect this vulnerability. Conservation efforts should seek to understand variations in species-level thermal tolerance across alpine microclimates.

The host bias of three epiphytic Aeridinae orchid species is reflected, but not explained, by mycorrhizal fungal associations.

PREMISE OF THE STUDY: The three co-occurring epiphytic orchid species, Sarcochilus hillii, Plectorrhiza tridentata, and Sarcochilus parviflorus vary in host specificity; all are found predominantly on the tree Backhousia myrtifolia but some also associate with a broad range of species. Despite this specialization, no fitness advantage has been detected for adult orchid plants growing on the preferred host. Therefore, we predicted that the host specialization of these orchid species is a consequence of a bias toward particular orchid mycorrhizal fungi, which are in turn biased toward particular woody plant species. METHODS: To test this hypothesis, we sampled representatives of each orchid species on B. myrtifolia and other host species across sites. Rhizoctonia-like fungi were isolated from orchid roots and identified using molecular markers. KEY RESULTS: Three groups of fungi were identified, and the orchid species varied in their specificity for these. All fungal groups were found on the host B. myrtifolia; yet at all sites, only one orchid species, S. hillii, associated with all three groups. CONCLUSIONS: Our results demonstrate that these orchid species did vary in their mycorrhizal specificity; however, the distribution of their mycorrhizal associates did not directly explain their host associations. Rather, we propose that the mycorrhizal relationship of these orchid species is complex and have suggested future avenues of research.

Phylogenetic influences on leaf trait integration in Pelargonium (Geraniaceae): convergence, divergence, and historical adaptation to a rapidly changing climate.

PREMISE OF THE STUDY: Trait integration may improve prediction of species and lineage responses to future climate change more than individual traits alone, particularly when analyses incorporate effects of phylogenetic relationships. The South African genus Pelargonium contains divergent major clades that have radiated along the same seasonal aridity gradient, presenting the opportunity to ask whether patterns of evolution in mean leaf trait values are achieved through the same set of coordinated changes among traits in each clade. METHODS: Seven leaf traits were measured on field-collected leaves from one-third of the species (98) of the genus. Trait relationships were examined using phylogenetic regression within major clades. Disparity analysis determined whether the course of trait evolution paralleled historical climate change events. KEY RESULTS: Divergence in mean trait values between sister clades A1 and A2 was consistent with expectations for leaves differing in longevity, despite strong similarity between clades in trait interactions. No traits in either clade exhibited significant relationships with multivariate climate axes, with one exception. Species in clades C and A2 included in this study occupied similar environments. These clades had similar values of individual trait means, except for δ(13)C, but they exhibited distinctive patterns of trait integration. CONCLUSIONS: Differing present-day patterns of trait integration are consistent with interpretations of adaptive responses to the prevailing climate at the time of each clade's origin. These differing patterns of integration are likely to exert strong effects on clade-level responses to future climate change in the winter rainfall region of South Africa.

The impact of beneficial plant-associated microbes on plant phenotypic plasticity.

Plants show phenotypic plasticity in response to changing or extreme abiotic environments; but over millions of years they also have co-evolved to respond to the presence of soil microbes. Studies on phenotypic plasticity in plants have focused mainly on the effects of the changing environments on plants' growth and survival. Evidence is now accumulating that the presence of microbes can alter plant phenotypic plasticity in a broad range of traits in response to a changing environment. In this review, we discuss the effects of microbes on plant phenotypic plasticity in response to changing environmental conditions, and how this may affect plant fitness. By using a range of specific plant-microbe interactions as examples, we demonstrate that one way that microbes can alleviate the effect of environmental stress on plants and thus increase plant fitness is to remove the stress, e.g., nutrient limitation, directly. Furthermore, microbes indirectly affect plant phenotypic plasticity and fitness through modulation of plant development and defense responses. In doing so, microbes affect fitness by both increasing or decreasing the degree of phenotypic plasticity, depending on the phenotype and the environmental stress studied, with no clear difference between the effect of prokaryotic and eukaryotic microbes in general. Additionally, plants have the ability to modulate microbial behaviors, suggesting that they manipulate bacteria, enhancing interactions that help them cope with stressful environments. Future challenges remain in the identification of the many microbial signals that modulate phenotypic plasticity, the characterization of plant genes, e.g. receptors, that mediate the microbial effects on plasticity, and the elucidation of the molecular mechanisms that link phenotypic plasticity with fitness. The characterization of plant and microbial mutants defective in signal synthesis or perception, together with carefully designed glasshouse or field experiments that test various environmental stresses will be necessary to understand the link between molecular mechanisms controlling plastic phenotypes with the resulting effects on plant fitness.

An innovative approach to using an intensive field course to build scientific and professional skills.

This paper reports on the design and evaluation of Field Studies in Functional Ecology (FSFE), a two-week intensive residential field course that enables students to master core content in functional ecology alongside skills that facilitate their transition from "student" to "scientist." We provide an overview of the course structure, showing how the constituent elements have been designed and refined over successive iterations of the course. We detail how FSFE students: (1) Work closely with discipline specialists to develop a small group project that tests an hypothesis to answer a genuine scientific question in the field; (2) Learn critical skills of data management and communication; and (3) Analyze, interpret, and present their results in the format of a scientific symposium. This process is repeated in an iterative "cognitive apprenticeship" model, supported by a series of workshops that name and explicitly instruct the students in "hard" and "soft" skills (e.g., statistics and teamwork, respectively) critically relevant for research and other careers. FSFE students develop a coherent and nuanced understanding of how to approach and execute ecological studies. The sophisticated knowledge and ecological research skills that they develop during the course is demonstrated through high-quality presentations and peer-reviewed publications in an open-access, student-led journal. We outline our course structure and evaluate its efficacy to show how this novel combination of field course elements allows students to gain maximum value from their educational journey, and to develop cognitive, affective, and reflective tools to help apply their skills as scientists.