A role for SPEECHLESS in the integration of leaf stomatal patterning with the growth vs disease trade-off in poplar

Admixture mapping and selection scans identify genomic regions associated with stomatal patterning and disease resistance in hybrid poplars

Variation in fitness components can be linked in some cases to variation in key traits. Metric traits that lie at the intersection of development, defense, and ecological interactions may be expected to experience environmental selection, informing our understanding of evolutionary and ecological processes. Here, we use quantitative genetic and population genomic methods to investigate disease dynamics in hybrid and non-hybrid populations. We focus our investigation on morphological and ecophysiological traits which inform our understanding of physiology, growth, and defense against a pathogen. In particular, we investigate stomata, microscopic pores on the surface of a leaf that regulate gas exchange during photosynthesis and are sites of entry for various plant pathogens. Stomatal patterning traits were highly predictive of disease risk. Admixture mapping identified a polygenic basis of disease resistance. Candidate genes for stomatal and disease resistance map to the same genomic regions and experienced positive selection. Genes with functions to guard cell homeostasis, the plant immune system, components of constitutive defenses, and growth-related transcription factors were identified. Our results indicate positive selection acted on candidate genes for stomatal patterning and disease resistance, potentially acting in concert to structure their variation in naturally formed backcrossing hybrid populations.

Experimental validation of the mechanism of stomatal development diversification

Stomata are the structures responsible for gas exchange in plants. The established framework for stomatal development is based on the model plant Arabidopsis, but diverse patterns of stomatal development have been observed in other plant lineages and species. The molecular mechanisms behind these diversified patterns are still poorly understood. We recently proposed a model for the molecular mechanisms of the diversification of stomatal development based on the genus Callitriche (Plantaginaceae), according to which a temporal shift in the expression of key stomatal transcription factors SPEECHLESS and MUTE leads to changes in the behavior of meristemoids (stomatal precursor cells). In the present study, we genetically manipulated Arabidopsis to test this model. By altering the timing of MUTE expression, we successfully generated Arabidopsis plants with early differentiation or prolonged divisions of meristemoids, as predicted by the model. The epidermal morphology of the generated lines resembled that of species with prolonged or no meristemoid divisions. Thus, the evolutionary process can be reproduced by varying the SPEECHLESS to MUTE transition. We also observed unexpected phenotypes, which indicated the participation of additional factors in the evolution of the patterns observed in nature. This study provides novel experimental insights into the diversification of meristemoid behaviors.

Responses to Drought Stress in Poplar: What Do We Know and What Can We Learn?

Poplar (Populus spp.) is a high-value crop for wood and biomass production and a model organism for tree physiology and genomics. The early release, in 2006, of the complete genome sequence of P. trichocarpa was followed by a wealth of studies that significantly enriched our knowledge of complex pathways inherent to woody plants, such as lignin biosynthesis and secondary cell wall deposition. Recently, in the attempt to cope with the challenges posed by ongoing climate change, fundamental studies and breeding programs with poplar have gradually shifted their focus to address the responses to abiotic stresses, particularly drought. Taking advantage from a set of modern genomic and phenotyping tools, these studies are now shedding light on important processes, including embolism formation (the entry and expansion of air bubbles in the xylem) and repair, the impact of drought stress on biomass yield and quality, and the long-term effects of drought events. In this review, we summarize the status of the research on the molecular bases of the responses to drought in poplar. We highlight how this knowledge can be exploited to select more tolerant genotypes and how it can be translated to other tree species to improve our understanding of forest dynamics under rapidly changing environmental conditions.

Leaf physiology variations are modulated by natural variations that underlie stomatal morphology in Populus

Stomata are essential for photosynthesis and abiotic stress tolerance. Here, we used multiomics approaches to dissect the genetic architecture and adaptive mechanisms that underlie stomatal morphology in Populus tomentosa juvenile natural population (303 accessions). We detected 46 candidate genes and 15 epistatic gene-pairs, associated with 5 stomatal morphologies and 18 leaf development and photosynthesis traits, through genome-wide association studies. Expression quantitative trait locus mapping revealed that stomata-associated gene loci were significantly associated with the expression of leaf-related genes; selective sweep analysis uncovered significant differentiation in the allele frequencies of genes that underlie stomatal variations. An allelic regulatory network operating under drought stress and adequate precipitation conditions, with three key regulators (DUF538, TRA2 and AbFH2) and eight interacting genes, was identified that might regulate leaf physiology via modulation of stomatal shape and density. Validation of candidate gene variations in drought-tolerant and F₁ hybrid populations of P. tomentosa showed that the DUF538, TRA2 and AbFH2 loci cause functional stabilisation of spatiotemporal regulatory, whose favourable alleles can be faithfully transmitted to offspring. This study provides insights concerning leaf physiology and stress tolerance via the regulation of stomatal determination in perennial plants.

Genetic architecture of the metabolic pathway of salicylic acid biosynthesis in Populus

Salicylic acid (SA) is a vital hormone for adaptive responses to biotic and abiotic stresses, which facilitates growth-immunity trade-offs in plants. However, the genetic regulatory networks underlying the metabolic pathway of SA biosynthesis in perennial species remain unclear. Here, we integrated genome-wide association study (GWAS) with metabolite and expression profiling methodologies to dissect the genetic architecture of SA biosynthesis in Populus. First, we quantified nine intermediate metabolites of SA biosynthesis in 300 unrelated Populus tomentosa Carr. individuals. Then, we used a systematic genetic strategy to identify candidate genes for constructing the genetic regulatory network of SA biosynthesis. We focused on WRKY70, an efficient transcription factor, as the key causal gene in the regulatory network, and combined the novel genes coordinating the accumulation of SA. Finally, we identified eight GWAS signals and eight expression quantitative trait loci situated in a selective sweep, and showed the presence of large allele frequency differences among the three geographic populations, revealing that candidate genes subject to selection were involved in SA biosynthesis. This study provides an integrated strategy for dissecting the genetic architecture of the metabolic pathway of SA biosynthesis in Populus, thereby enhancing our understanding of genetic regulation of SA biosynthesis in trees, and accelerating marker-assisted breeding efforts toward high-resistance elite varieties of Populus.

Callitriche as a potential model system for evolutionary studies on the dorsiventral distribution of stomata

Controlling the distribution of stomata is crucial for the adaptation of plants to new, or changing environments. While many plant species produce stomata predominantly on the abaxial leaf surface (hypostomy), some produce stomata on both surfaces (amphistomy), and the remaining few produce them only on the adaxial surface (hyperstomy). Various selective pressures have driven the evolution of these three modes of stomatal distribution. Despite recent advances in our understanding of stomatal development and dorsiventral leaf polarity, the genetic basis for the evolution of different stomatal distributions is still unclear. Here, we propose the genus Callitriche as a new model system to investigate patterns in the evolution of stomatal distribution. Callitriche comprises species with diverse lifestyles, including terrestrial, amphibious, and obligately aquatic plants. We found that species in this genus cover all three modes of dorsiventral stomatal distribution, making it a desirable model for comparative and evolutionary analyses on distribution modes. We further characterized the genetic basis of the different distribution modes, focusing on the stomatal key transcription factor SPEECHLESS. Future research using the promising model system Callitriche would open a new direction for evolutionary developmental biology studies on stomata.

Growth-defense trade-offs masked in unadmixed populations are revealed by hybridization

Organisms are constantly challenged by pathogens and pests, which can drive the evolution of growth-defense strategies. Plant stomata are essential for gas exchange during photosynthesis and conceptually lie at the intersection of the physiological demands of growth and exposure to foliar fungal pathogens. Generations of natural selection for locally adapted growth-defense strategies can eliminate variation between traits, potentially masking trade-offs and selection conflicts that may have existed in the past. Hybrid populations offer a unique opportunity to reset the clock on selection and to study potentially maladaptive trait variation before selection removes it. We study the interactions of growth, stomatal, ecopysiological, and disease resistance traits in poplars (Populus) after infection by the leaf rust Melampsora medusae. Phenotypes were measured in a common garden and genotyped at 227K SNPs. We isolate the effects of hybridization on trait variance, discover correlations between stomatal, ecophysiology, and disease resistance, examine trade-offs and selection conflicts, and explore the evolution of growth-defense strategies potentially mediated by selection for stomatal traits on the upper leaf surface. These results suggest an important role for stomata in determining growth-defense strategies in organisms susceptible to foliar pathogens, and reinforces the contribution of hybridization studies toward our understanding of trait evolution.

A Stomatal Model of Anatomical Tradeoffs Between Gas Exchange and Pathogen Colonization

Stomatal pores control leaf gas exchange and are one route for infection of internal plant tissues by many foliar pathogens, setting up the potential for tradeoffs between photosynthesis and pathogen colonization. Anatomical shifts to lower stomatal density and/or size may also limit pathogen colonization, but such developmental changes could permanently reduce the gas exchange capacity for the life of the leaf. I developed and analyzed a spatially explicit model of pathogen colonization on the leaf as a function of stomatal size and density, anatomical traits which partially determine maximum rates of gas exchange. The model predicts greater stomatal size or density increases the probability of colonization, but the effect is most pronounced when the fraction of leaf surface covered by stomata is low. I also derived scaling relationships between stomatal size and density that preserves a given probability of colonization. These scaling relationships set up a potential anatomical conflict between limiting pathogen colonization and minimizing the fraction of leaf surface covered by stomata. Although a connection between gas exchange and pathogen defense has been suggested empirically, this is the first mathematical model connecting gas exchange and pathogen defense via stomatal anatomy. A limitation of the model is that it does not include variation in innate immunity and stomatal closure in response to pathogens. Nevertheless, the model makes predictions that can be tested with experiments and may explain variation in stomatal size and density among plants. The model is generalizable to many types of pathogens, but lacks significant biological realism that may be needed for precise predictions.

Stomatal development in the grasses: lessons from models and crops (and crop models)

When plants emerged from their aquatic origins to colonise land, they needed to avoid desiccation while still enabling gas and water exchange with the environment. The solution was the development of a waxy cuticle interrupted by epidermal pores, known as stomata. Despite the importance of stomata in plant physiology and their contribution to global water and carbon cycles, our knowledge of the genetic basis of stomatal development is limited mostly to the model dicot, Arabidopsis thaliana. This limitation is particularly troublesome when evaluating grasses, whose members represent our most agriculturally significant crops. Grass stomatal development follows a trajectory strikingly different from Arabidopsis and their uniquely shaped four-celled stomatal complexes are especially responsive to environmental inputs. Thus, understanding the development and regulation of these efficient complexes is of particular interest for the purposes of crop engineering. This review focuses on genetic regulation of grass stomatal development and prospects for the future, highlighting discoveries enabled by parallel comparative investigations in cereal crops and related genetic model species such as Brachypodium distachyon.

Discovery of salicyl benzoate UDP-glycosyltransferase, a central enzyme in poplar salicinoid phenolic glycoside biosynthesis

The salicinoids are anti-herbivore phenolic glycosides unique to the Salicaceae (Populus and Salix). They consist of a salicyl alcohol glucoside core, which is usually further acylated with benzoic, cinnamic or phenolic acids. While salicinoid structures are well known, their biosynthesis remains enigmatic. Recently, two enzymes from poplar, salicyl alcohol benzoyl transferase and benzyl alcohol benzoyl transferase, were shown to catalyze the production of salicyl benzoate, a predicted potential intermediate in salicinoid biosynthesis. Here, we used transcriptomics and co-expression analysis with these two genes to identify two UDP-glucose-dependent glycosyltransferases (UGT71L1 and UGT78M1) as candidate enzymes in this pathway. Both recombinant enzymes accepted only salicyl benzoate, salicylaldehyde and 2-hydroxycinnamic acid as glucose acceptors. Knocking out the UGT71L1 gene by CRISPR/Cas9 in poplar hairy root cultures led to the complete loss of salicortin, tremulacin and tremuloidin, and a partial reduction of salicin content. This demonstrated that UGT71L1 is required for synthesis of the major salicinoids, and suggested that an additional route can lead to salicin. CRISPR/Cas9 knockouts for UGT78M1 were not successful, and its in vivo role thus remains to be determined. Although it has a similar substrate preference and predicted structure as UGT71L1, it appears not to contribute to the synthesis of salicortin, tremulacin and tremuloidin, at least in roots. The demonstration of UGT71L1 as an enzyme of salicinoid biosynthesis will open up new avenues for the elucidation of this pathway.

Emerging roles for carbonic anhydrase in mesophyll conductance and photosynthesis

Carbonic anhydrase (CA) is an abundant protein in most photosynthesizing organisms and higher plants. This review paper considers the physiological importance of the more abundant CA isoforms in photosynthesis, through their effects on CO₂ diffusion and other processes in photosynthetic organisms. In plants, CA has multiple isoforms in three different families (α, β and γ) and is mainly known to catalyze the CO₂↔HCO3- equilibrium. This reversible conversion has a clear role in photosynthesis, primarily through sustaining the CO₂ concentration at the site of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco). Despite showing the same major reaction mechanism, the three main CA families are evolutionarily distinct. For different CA isoforms, cellular localization and total gene expression as a function of developmental stage are predicted to determine the role of each family in relation to the net assimilation rate. Reaction-diffusion modeling and observational evidence support a role for CA activity in reducing resistance to CO₂ diffusion inside mesophyll cells by facilitating CO₂ transfer in both gas and liquid phases. In addition, physical and/or biochemical interactions between CAs and other membrane-bound compartments, for example aquaporins, are suggested to trigger a CO₂ -sensing response by stomatal movement. In response to environmental stresses, changes in the expression level of CAs and/or stimulated deactivation of CAs may correspond with lower photosynthetic capacity. We suggest that further studies should focus on the dynamics of the relationship between the activity of CAs (with different subcellular localization, abundance and gene expression) and limitations due to CO₂ diffusivity through the mesophyll and supply of CO₂ to photosynthetic reactions.

Genome-Wide Association Study of Wood Anatomical and Morphological Traits in Populus trichocarpa

To understand the genetic mechanisms underlying wood anatomical and morphological traits in Populus trichocarpa, we used 869 unrelated genotypes from a common garden in Clatskanie, Oregon that were previously collected from across the distribution range in western North America. Using GEMMA mixed model analysis, we tested for the association of 25 phenotypic traits and nine multitrait combinations with 6.741 million SNPs covering the entire genome. Broad-sense trait heritabilities ranged from 0.117 to 0.477. Most traits were significantly correlated with geoclimatic variables suggesting a role of climate and geography in shaping the variation of this species. Fifty-seven SNPs from single trait GWAS and 11 SNPs from multitrait GWAS passed an FDR threshold of 0.05, leading to the identification of eight and seven nearby candidate genes, respectively. The percentage of phenotypic variance explained (PVE) by the significant SNPs for both single and multitrait GWAS ranged from 0.01% to 6.18%. To further evaluate the potential roles of candidate genes, we used a multi-omic network containing five additional data sets, including leaf and wood metabolite GWAS layers and coexpression and comethylation networks. We also performed a functional enrichment analysis on coexpression nearest neighbors for each gene model identified by the wood anatomical and morphological trait GWAS analyses. Genes affecting cell wall composition and transport related genes were enriched in wood anatomy and stomatal density trait networks. Signaling and metabolism related genes were also common in networks for stomatal density. For leaf morphology traits (leaf dry and wet weight) the networks were significantly enriched for GO terms related to photosynthetic processes as well as cellular homeostasis. The identified genes provide further insights into the genetic control of these traits, which are important determinants of the suitability and sustainability of improved genotypes for lignocellulosic biofuel production.