Likelihood of using phytoremediation for mine-tailing management in Chile.

Chile has more than 750 mine tailings across the country, mainly distributed in the northern region, which also includes a biodiversity hotspot and the driest desert in the world. So far, tailing management has included chemical and physical stabilization of tailings, exclusively. This research examined the perceived likelihood of stakeholders, namely: Academia, Industry, environmental Non-Governmental Organizations, and Government officials, in the management of tailings and explored their attitudes toward the inclusion of plants and their associated microbes, as an additional stabilization technology, through an online questionnaire (n = 43). We aimed to answer What was the perceived likelihood of stakeholders for using phytoremediation for mine-tailing management in Chile? The consensus opinion across sectors is that there is general lack of proper tailings management in Chile. There is a critical willingness to incorporate phytoremediation, with a preference for using native plants. An absence of demonstrably successful phytoremediation projects, funding and access to suitable plants are perceived to be significant limitations. Local community involvement is considered key to successful restoration of tailings.

The novelty of this research work comes from finding consensus opinion across sectors about the lack of proper management of tailings in Chile and a critical willingness to incorporate phytoremediation as an additional stabilization technology. Furthermore, it provides insights from professional stakeholders on the need for more studies on native plants that can be used in phytoremediation, to diversify the representation of stakeholders, and to include input from local communities, directly affected by mines.

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 CO2 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 CO2↔HCO3- equilibrium. This reversible conversion has a clear role in photosynthesis, primarily through sustaining the CO2 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 CO2 diffusion inside mesophyll cells by facilitating CO2 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 CO2 -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 CO2 diffusivity through the mesophyll and supply of CO2 to photosynthetic reactions.

Growth response, uptake and mobilization of metals in native plant species on tailings at a Chilean copper mine.

Chile has many mine tailing deposits available for phytoremediation through the establishment of metal-tolerant plants. To guide such efforts, it is necessary to know whether roots exclude or take up metals, or if metals are mobilized to shoots. We evaluated a polyculture of ten native species 6 years after they were planted directly into tailings, amended with mycorrhiza before planting or planted with compost. All species were assessed for survival and hare damage. Growth, vigor, chlorophyll content index, and chlorophyll fluorescence were measured in seven species, and root and shoot concentrations of 11 metals and As were measured in five. All species had some level of stress, but there was no clear pattern related to foliar metal concentrations or amendments, the latter also having little effect on growth or survival. Copper, Fe, and Mo concentrations exceeded threshold toxicity levels for leaves and for animal feed in several species. Copper was most concentrated in the roots, and readily mobilized to leaves in three of the five species tested, but not in Prosopis chilensis (algarrobo) and Quillaja saponaria (quillay). Because of lower uptake overall, quillay is recommended for stabilizing tailings while also reducing transport to aerial tissues.

Proteomic analysis of metabolic mechanisms associated with fatty acid biosynthesis during Styrax tonkinensis kernel development.

BACKGROUND: Styrax tonkinensis is a white-flowered tree with considerable potential as a feedstock source for biodiesel production from the oily seed contained within its nutlike drupes. Transcriptome changes during oil accumulation have been previously reported, but not concurrent changes in the proteome. RESULTS: Using proteomic analysis of samples collected at 50, 70, 100 and 130 days after flowering (DAF), we identified 1472 differentially expressed proteins (DEPs). Based on their expression patterns, we grouped the DEPs into nine clusters and analyzed the pathway enrichment. Proteins related to starch and sucrose metabolism were most abundant at 50 DAF. Proteins involved in fatty acid (FA) biosynthesis were mainly grouped into a cluster that peaked at 70 DAF. Proteins related to protein processing in endoplasmic reticulum had two major patterns, trending either upwards or downwards, while proteins involved in amino acid biosynthesis showed more complex relationships. We identified 42 key enzymes involved in lipid accumulation during kernel development, including the acetyl-CoA carboxylase complex (ACC) and the pyruvate dehydrogenase complex (PDC). One oil body membrane protein, oleosin, continuously increased during kernel development. CONCLUSION: A regulatory network of oil accumulation processes was built based on protein and available transcriptome expression data, which were in good temporal agreement. This analysis placed ACC and PDC in the center of the network, suggesting that the glycolytic provision of substrate plays a central regulatory role in FA biosynthesis and oil accumulation. © 2021 Society of Chemical Industry.

Transcriptome analysis of metabolic pathways associated with oil accumulation in developing seed kernels of Styrax tonkinensis, a woody biodiesel species.

BACKGROUND: Styrax tonkinensis (Pierre) Craib ex Hartwich has great potential as a woody biodiesel species having seed kernels with high oil content, excellent fatty acid composition and good fuel properties. However, no transcriptome information is available on the molecular regulatory mechanism of oil accumulation in developing S. tonkinensis kernels. RESULTS: The dynamic patterns of oil content and fatty acid composition at 11 time points from 50 to 150 days after flowering (DAF) were analyzed. The percent oil content showed an up-down-up pattern, with yield and degree of unsaturation peaking on or after 140 DAF. Four time points (50, 70, 100, and 130 DAF) were selected for Illumina transcriptome sequencing. Approximately 73 million high quality clean reads were generated, and then assembled into 168,207 unigenes with a mean length of 854 bp. There were 5916 genes that were differentially expressed between different time points. These differentially expressed genes were grouped into 9 clusters based on their expression patterns. Expression patterns of a subset of 12 unigenes were confirmed by qRT-PCR. Based on their functional annotation through the Basic Local Alignment Search Tool and publicly available protein databases, specific unigenes encoding key enzymes, transmembrane transporters, and transcription factors associated with oil accumulation were determined. Three main patterns of expression were evident. Most unigenes peaked at 70 DAF, coincident with a rapid increase in oil content during kernel development. Unigenes with high expression at 50 DAF were associated with plastid formation and earlier stages of oil synthesis, including pyruvate and acetyl-CoA formation. Unigenes associated with triacylglycerol biosynthesis and oil body development peaked at 100 or 130 DAF. CONCLUSIONS: Transcriptome changes during oil accumulation show a distinct temporal trend with few abrupt transitions. Expression profiles suggest that acetyl-CoA formation for oil biosynthesis is both directly from pyruvate and indirectly via acetaldehyde, and indicate that the main carbon source for fatty acid biosynthesis is triosephosphate originating from phosphohexose outside the plastid. Different sn-glycerol-3-phosphate acyltransferases are implicated in diacylglycerol biosynthesis at early versus late stages of oil accumulation. Triacylglycerol biosynthesis may be accomplished by both diacylglycerol and by phospholipid:diacylglycerol acyltransferases.

Isotopic composition and concentration of total nitrogen and nitrate in xylem sap under near steady-state hydroponics.

After root uptake, nitrate is effluxed back to the medium, assimilated locally, or translocated to shoots. Rooted black cottonwood (Populus trichocarpa) scions were supplied with a NO3- -based (0.5 mM) nutrient medium of known isotopic composition (δ15 N), and xylem sap was collected by pressure bombing. To establish a sampling protocol, sap was collected from lower and upper stem sections at 0.1-0.2 MPa above the balancing pressure, and after increasing the pressure by a further 0.5 MPa. Xylem sap from upper stem sections was partially diluted at higher pressure. Further analysis was restricted to sap obtained from intact shoots at low pressure. Total-, NO3- -N and, by difference, organic-N concentrations ranged from 6.1-11.0, 1.2-2.4, and 4.6-9.4 mM, while discrimination relative to the nutrient medium was -6.3 to 0.5‰, -23.3 to -11.5‰ and - 1.3 to 4.9‰, respectively. There was diurnal variation in δ15 N of total- and organic-N, but not NO3- . The difference in δ15 N between xylem NO3- and organic-N suggests that discrimination by nitrate reductase is near 25.1 ± 1.6‰. When this value was used in an isotope mass balance model, the predicted xylem sap NO3- -N to total-N ratio closely matched direct measurement.

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

Occurrence of stomata on both leaf surfaces (amphistomaty) promotes higher stomatal conductance and photosynthesis while simultaneously increasing exposure to potential disease agents in black cottonwood (Populus trichocarpa). A genome-wide association study (GWAS) with 2.2M single nucleotide polymorphisms generated through whole-genome sequencing found 280 loci associated with variation in adaxial stomatal traits, implicating genes regulating stomatal development and behavior. Strikingly, numerous loci regulating plant growth and response to biotic and abiotic stresses were also identified. The most significant locus was a poplar homologue of SPEECHLESS (PtSPCH1). Individuals possessing PtSPCH1 alleles associated with greater adaxial stomatal density originated primarily from environments with shorter growing seasons (e.g. northern latitudes, high elevations) or with less precipitation. PtSPCH1 was expressed in developing leaves but not developing stem xylem. In developing leaves, RNA sequencing showed patterns of coordinated expression between PtSPCH1 and other GWAS-identified genes. The breadth of our GWAS results suggests that the evolution of amphistomaty is part of a larger, complex response in plants. Suites of genes underpin this response, retrieved through genetic association to adaxial stomata, and show coordinated expression during development. We propose that the occurrence of amphistomaty in P. trichocarpa involves PtSPCH1 and reflects selection for supporting rapid growth over investment in immunity.

Ecological genomics of variation in bud-break phenology and mechanisms of response to climate warming in Populus trichocarpa.

Spring bud-break phenology is a critical adaptive feature common to temperate perennial woody plants. Understanding the molecular underpinnings of variation in bud-break is important for elucidating adaptive evolution and predicting outcomes relating to climate change. Field and controlled growth chamber tests were used to assess population-wide patterns in bud-break from wild-sourced black cottonwood (Populus trichocarpa) genotypes. We conducted a genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) derived from whole genome sequencing to test for loci underlying variation in bud-break. Bud-break had a quadratic relationship with latitude, where southern- and northern-most provenances generally broke bud earlier than those from central parts of the species' range. Reduced winter chilling increased population-wide variation in bud-break, whereas greater chilling decreased variation. GWAS uncovered 16 loci associated with bud-break. Phenotypic changes connected with allelic variation were replicated in an independent set of P. trichocarpa trees. Despite phenotypic similarities, genetic profiles between southern- and northern-most genotypes were dissimilar based on our GWAS-identified SNPs. We propose that the GWAS-identified loci underpin the geographical pattern in P. trichocarpa and that variation in bud-break reflects different selection for winter chilling and heat sum accumulation, both of which can be affected by climate warming.

Mechanosensitive Adhesion Explains Stepping Motility in Amoeboid Cells.

Cells employing amoeboid motility exhibit repetitive cycles of rapid expansion and contraction and apply coordinated traction forces to their environment. Although aspects of this process are well studied, it is unclear how the cell controls the coordination of cell length changes with adhesion to the surface. Here, we develop a simple model to mechanistically explain the emergence of periodic changes in length and spatiotemporal dynamics of traction forces measured in chemotaxing unicellular amoeba, Dictyostelium discoideum. In contrast to the biochemical mechanisms that have been implicated in the coordination of some cellular processes, we show that many features of amoeboid locomotion emerge from a simple mechanochemical model. The mechanism for interaction with the environment in Dictyostelium is unknown and thus, we explore different cell-environment interaction models to reveal that mechanosensitive adhesions are necessary to reproduce the spatiotemporal adhesion patterns. In this modeling framework, we find that the other motility modes, such as smooth gliding, arise naturally with variations in the physical properties of the surface. Thus, our work highlights the prominent role of biomechanics in determining the emergent features of amoeboid locomotion.

Substantial role for carbonic anhydrase in latitudinal variation in mesophyll conductance of Populus trichocarpa Torr. & Gray.

In Populus trichocarpa (black cottonwood), net photosynthesis (An ) varies with latitude and, in northern genotypes, is supported by higher stomatal conductance (gs ). We report here a parallel cline in mesophyll conductance (gm ) and link this variation to carbonic anhydrase (CA) activity. Using concurrent carbon isotope discrimination and chlorophyll fluorescence methods, we examined the effects of acetazolamide, an inhibitor of CA, on gm in six representative genotypes (three from either end of the north-south cline). Acetazolamide reduced CA activity, gm , gs , chloroplast CO2 concentration (Cc ) and An at normal CO2 (400 µmol mol-1 ), the latter being reversible at saturating CO2 . Absolute reductions in An , gm and CA activity were greater in northern genotypes than in southern genotypes (P < 0.025) but percent reductions were similar. In contrast, northern genotypes showed lower percent reduction in Cc compared to southern genotypes (P < 0.025). The northern genotypes had greater CA activity relative to both leaf area (two-fold) and mass (1.8-fold) (P < 0.016). The relationship between CA activity and gm was similar whether the variation was inherent or inhibitor induced. We suggest that greater CA activity contributes to higher gm in northern P. trichocarpa genotypes, but other diffusion pathway components may also be involved.

Self-organized mechano-chemical dynamics in amoeboid locomotion of Physarum fragments.

The aim of this work is to quantify the spatio-temporal dynamics of flow-driven amoeboid locomotion in small (~100 µm) fragments of the true slime mold Physarum polycephalum. In this model organism, cellular contraction drives intracellular flows, and these flows transport the chemical signals that regulate contraction in the first place. As a consequence of these non-linear interactions, a diversity of migratory behaviors can be observed in migrating Physarum fragments. To study these dynamics, we measure the spatio-temporal distributions of the velocities of the endoplasm and ectoplasm of each migrating fragment, the traction stresses it generates on the substratum, and the concentration of free intracellular calcium. Using these unprecedented experimental data, we classify migrating Physarum fragments according to their dynamics, finding that they often exhibit spontaneously coordinated waves of flow, contractility and chemical signaling. We show that Physarum fragments exhibiting symmetric spatio-temporal patterns of endoplasmic flow migrate significantly slower than fragments with asymmetric patterns. In addition, our joint measurements of ectoplasm velocity and traction stress at the substratum suggest that forward motion of the ectoplasm is enabled by a succession of stick-slip transitions, which we conjecture are also organized in the form of waves. Combining our experiments with a simplified convection-diffusion model, we show that the convective transport of calcium ions may be key for establishing and maintaining the spatiotemporal patterns of calcium concentration that regulate the generation of contractile forces.

Neural mechanism of optimal limb coordination in crustacean swimming.

A fundamental challenge in neuroscience is to understand how biologically salient motor behaviors emerge from properties of the underlying neural circuits. Crayfish, krill, prawns, lobsters, and other long-tailed crustaceans swim by rhythmically moving limbs called swimmerets. Over the entire biological range of animal size and paddling frequency, movements of adjacent swimmerets maintain an approximate quarter-period phase difference with the more posterior limbs leading the cycle. We use a computational fluid dynamics model to show that this frequency-invariant stroke pattern is the most effective and mechanically efficient paddling rhythm across the full range of biologically relevant Reynolds numbers in crustacean swimming. We then show that the organization of the neural circuit underlying swimmeret coordination provides a robust mechanism for generating this stroke pattern. Specifically, the wave-like limb coordination emerges robustly from a combination of the half-center structure of the local central pattern generating circuits (CPGs) that drive the movements of each limb, the asymmetric network topology of the connections between local CPGs, and the phase response properties of the local CPGs, which we measure experimentally. Thus, the crustacean swimmeret system serves as a concrete example in which the architecture of a neural circuit leads to optimal behavior in a robust manner. Furthermore, we consider all possible connection topologies between local CPGs and show that the natural connectivity pattern generates the biomechanically optimal stroke pattern most robustly. Given the high metabolic cost of crustacean swimming, our results suggest that natural selection has pushed the swimmeret neural circuit toward a connection topology that produces optimal behavior.

Intracellular Pressure Dynamics in Blebbing Cells.

Blebs are pressure-driven protrusions that play an important role in cell migration, particularly in three-dimensional environments. A bleb is initiated when the cytoskeleton detaches from the cell membrane, resulting in the pressure-driven flow of cytosol toward the area of detachment and local expansion of the cell membrane. Recent experiments involving blebbing cells have led to conflicting hypotheses regarding the timescale of intracellular pressure propagation. The interpretation of one set of experiments supports a poroelastic model of the cytoplasm that leads to slow pressure equilibration when compared to the timescale of bleb expansion. A different study concludes that pressure equilibrates faster than the timescale of bleb expansion. To address this discrepancy, a dynamic computational model of the cell was developed that includes mechanics of and the interactions among the cytoplasm, the actin cortex, the cell membrane, and the cytoskeleton. The model results quantify the relationship among cytoplasmic rheology, pressure, and bleb expansion dynamics, and provide a more detailed picture of intracellular pressure dynamics. This study shows the elastic response of the cytoplasm relieves pressure and limits bleb size, and that both permeability and elasticity of the cytoplasm determine bleb expansion time. Our model with a poroelastic cytoplasm shows that pressure disturbances from bleb initiation propagate faster than the timescale of bleb expansion and that pressure equilibrates slower than the timescale of bleb expansion. The multiple timescales in intracellular pressure dynamics explain the apparent discrepancy in the interpretation of experimental results.

Variation in fluxes estimated from nitrogen isotope discrimination corresponds with independent measures of nitrogen flux in Populus balsamifera L.

Acquisition of mineral nitrogen by roots from the surrounding environment is often not completely efficient, in which a variable amount of leakage (efflux) relative to gross uptake (influx) occurs. The efflux/influx ratio (E/I) is, therefore, inversely related to the efficiency of nutrient uptake at the root level. Time-integrated estimates of E/I and other nitrogen-use traits may be obtainable from variation in stable isotope ratios or through compartmental analysis of tracer efflux (CATE) using radioactive or stable isotopes. To compare these two methods, Populus balsamifera L. genotypes were selected, a priori, for high or low nitrogen isotope discrimination. Vegetative cuttings were grown hydroponically, and E/I was calculated using an isotope mass balance model (IMB) and compared to E/I calculated using (15) N CATE. Both methods indicated that plants grown with ammonium had greater E/I than nitrate-grown plants. Genotypes with high or low E/I using CATE also had similarly high or low estimates of E/I using IMB, respectively. Genotype-specific means were linearly correlated (r = 0.77; P = 0.0065). Discrepancies in E/I between methods may reflect uncertainties in discrimination factors for the assimilatory enzymes, or temporal differences in uptake patterns. By utilizing genotypes with known variation in nitrogen isotope discrimination, a relationship between nitrogen isotope discrimination and bidirectional nitrogen fluxes at the root level was observed.

Extensive functional pleiotropy of REVOLUTA substantiated through forward genetics.

In plants, genes may sustain extensive pleiotropic functional properties by individually affecting multiple, distinct traits. We discuss results from three genome-wide association studies of approximately 400 natural poplar (Populus trichocarpa) accessions phenotyped for 60 ecological/biomass, wood quality, and rust fungus resistance traits. Single-nucleotide polymorphisms (SNPs) in the poplar ortholog of the class III homeodomain-leucine zipper transcription factor gene REVOLUTA (PtREV) were significantly associated with three specific traits. Based on SNP associations with fungal resistance, leaf drop, and cellulose content, the PtREV gene contains three potential regulatory sites within noncoding regions at the gene's 3' end, where alternative splicing and messenger RNA processing actively occur. The polymorphisms in this region associated with leaf abscission and cellulose content are suggested to represent more recent variants, whereas the SNP associated with leaf rust resistance may be more ancient, consistent with REV's primary role in auxin signaling and its functional evolution in supporting fundamental processes of vascular plant development.

Actin-myosin spatial patterns from a simplified isotropic viscoelastic model.

F-actin networks are involved in cell mechanical processes ranging from motility to endocytosis. The mesoscale architecture of assemblies of individual F-actin polymers that gives rise to micrometer-scale rheological properties is poorly understood, despite numerous in vivo and vitro studies. In vitro networks have been shown to organize into spatial patterns when spatially confined, including dense spherical shells inside spherical emulsion droplets. Here we develop a simplified model of an isotropic, compressible, viscoelastic material continually assembling and disassembling. We demonstrate that spherical shells emerge naturally when the strain relaxation rate (corresponding to internal network reorganization) is slower than the disassembly rate (corresponding to F-actin depolymerization). These patterns are consistent with recent experiments, including a collapse of shells to a central high-density focus of F-actin when either assembly or disassembly is reduced with drugs. Our results demonstrate how complex spatio-temporal patterns can emerge without spatially distributed force generation, polar alignment of F-actin polymers, or spatially nonuniform regulation of F-actin by upstream biochemical networks.

Geographical and environmental gradients shape phenotypic trait variation and genetic structure in Populus trichocarpa.

• Populus trichocarpa is widespread across western North America spanning extensive variation in photoperiod, growing season and climate. We investigated trait variation in P. trichocarpa using over 2000 trees from a common garden at Vancouver, Canada, representing replicate plantings of 461 genotypes originating from 136 provenance localities. • We measured 40 traits encompassing phenological events, biomass accumulation, growth rates, and leaf, isotope and gas exchange-based ecophysiology traits. With replicated plantings and 29,354 single nucleotide polymorphisms (SNPs) from 3518 genes, we estimated both broad-sense trait heritability (H(2)) and overall population genetic structure from principal component analysis. • Populus trichocarpa had high phenotypic variation and moderate/high H(2) for many traits. H(2) ranged from 0.3 to 0.9 in phenology, 0.3 to 0.8 in biomass and 0.1 to 0.8 in ecophysiology traits. Most traits correlated strongly with latitude, maximum daylength and temperature of tree origin, but not necessarily with elevation, precipitation or heat : moisture indices. Trait H(2) values reflected trait correlation strength with geoclimate variables. The population genetic structure had one significant principal component (PC1) which correlated with daylength and showed enrichment for genes relating to circadian rhythm and photoperiod. • Robust relationships between traits, population structure and geoclimate in P. trichocarpa reflect patterns which suggest that range-wide geographical and environment gradients have shaped its genotypic and phenotypic variability.

Genome-wide association implicates numerous genes underlying ecological trait variation in natural populations of Populus trichocarpa.

In order to uncover the genetic basis of phenotypic trait variation, we used 448 unrelated wild accessions of black cottonwood (Populus trichocarpa) from much of its range in western North America. Extensive data from large-scale trait phenotyping (with spatial and temporal replications within a common garden) and genotyping (with a 34 K Populus single nucleotide polymorphism (SNP) array) of all accessions were used for gene discovery in a genome-wide association study (GWAS). We performed GWAS with 40 biomass, ecophysiology and phenology traits and 29,355 filtered SNPs representing 3518 genes. The association analyses were carried out using a Unified Mixed Model accounting for population structure effects among accessions. We uncovered 410 significant SNPs using a Bonferroni-corrected threshold (P<1.7×10(-6)). Markers were found across 19 chromosomes, explained 1-13% of trait variation, and implicated 275 unique genes in trait associations. Phenology had the largest number of associated genes (240 genes), followed by biomass (53 genes) and ecophysiology traits (25 genes). The GWAS results propose numerous loci for further investigation. Many traits had significant associations with multiple genes, underscoring their genetic complexity. Genes were also identified with multiple trait associations within and/or across trait categories. In some cases, traits were genetically correlated while in others they were not.

Association genetics, geography and ecophysiology link stomatal patterning in Populus trichocarpa with carbon gain and disease resistance trade-offs.

Stomata are essential for diffusive entry of gases to support photosynthesis, but may also expose internal leaf tissues to pathogens. To uncover trade-offs in range-wide adaptation relating to stomata, we investigated the underlying genetics of stomatal traits and linked variability in these traits with geoclimate, ecophysiology, condensed foliar tannins and pathogen susceptibility in black cottonwood (Populus trichocarpa). Upper (adaxial) and lower (abaxial) leaf stomatal traits were measured from 454 accessions collected throughout much of the species range. We calculated broad-sense heritability (H(2) ) of stomatal traits and, using SNP data from a 34K Populus SNP array, performed a genome-wide association studies (GWAS) to uncover genes underlying stomatal trait variation. H(2) values for stomatal traits were moderate (average H(2) = 0.33). GWAS identified genes associated primarily with adaxial stomata, including polarity genes (PHABULOSA), stomatal development genes (BRASSINOSTEROID-INSENSITIVE 2) and disease/wound-response genes (GLUTAMATE-CYSTEINE LIGASE). Stomatal traits correlated with latitude, gas exchange, condensed tannins and leaf rust (Melampsora) infection. Latitudinal trends of greater adaxial stomata numbers and guard cell pore size corresponded with higher stomatal conductance (gs ) and photosynthesis (Amax ), faster shoot elongation, lower foliar tannins and greater Melampsora susceptibility. This suggests an evolutionary trade-off related to differing selection pressures across the species range. In northern environments, more adaxial stomata and larger pore sizes reflect selection for rapid carbon gain and growth. By contrast, southern genotypes have fewer adaxial stomata, smaller pore sizes and higher levels of condensed tannins, possibly linked to greater pressure from natural leaf pathogens, which are less significant in northern ecosystems.

Mechanisms of elastic enhancement and hindrance for finite-length undulatory swimmers in viscoelastic fluids.

A computational model of finite-length undulatory swimmers is used to examine the physical origin of the effect of elasticity on swimming speed. We explore two distinct target swimming strokes: one derived from the motion of Caenorhabditis elegans, with large head undulations, and a contrasting stroke with large tail undulations. We show that both favorable stroke asymmetry and swimmer elasticity contribute to a speed-up, but a substantial boost results only when these two effects work together. We reproduce conflicting results from the literature simply by changing relevant physical parameters.