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.

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.

Author Correction: Sexual homomorphism in dioecious trees: extensive tests fail to detect sexual dimorphism in Populus.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

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.

Integrating morphological characters, molecular markers, and distribution patterns to assess the identity of Blepharis species from Jordan.

BACKGROUND: Blepharis constitutes an important part of the vegetation of the Jordanian arid and semi-arid regions, yet whether one or more species of this genus occurs in the Jordanian area is uncertain. We addressed this question by assessing morphological characters and testing Inter-Simple Sequence Repeat (ISSR) markers from three populations of Blepharis: two northern (lower slopes of Kufranjah valley and the Dead Sea region) and one southern (Wadi al Yutm). RESULTS: Shoots from randomly chosen Blepharis plants were harvested from each of the three populations for morphological and molecular analyses. In the northern populations, spikes were lax and bract width was significantly shorter than length of the longest lateral spine compared to the southern population. A multivariate linear discriminant analysis distinguished the northern populations from the southern one by internode length, bract width, longest lateral spine length, and bract width to spine length ratio. The ISSR analysis revealed that 44 markers across eight primers were polymorphic with major allele frequency of 83.6% and an average of 5.5 polymorphic markers per primer. The genetic resemblance among individuals ranged from 0.27 to 0.96. The three Blepharis populations were accordingly clustered into two distinct groups, similar to the analysis of morphological differences and corresponding with the "northern" and "southern" population designations. CONCLUSIONS: Our results strongly indicate the occurrence of two discrete Blepharis species in Jordan and reject the hypothesis that the genus is represented by only one species. We propose that the Blepharis species in Jordan are B. attenutata Napper (represented by the northern populations) and B. ciliaris (L.) B. L. Burtt (represented by the southern population). These findings are important for informing and revising floristic work within the region and an updated key has been included in our findings.

C4 photosynthesis and transition of Kranz anatomy in cotyledons and leaves of Tetraena simplex.

PREMISE OF THE STUDY: Tetraena simplex is an independently evolved C4 species in the Zygophylloideae (Zygophyllaceae) and a characteristic forb of saline flats in hot and sandy desert habitats. During early ontogeny, the species had a morphological shift from planar cotyledons (dorsiventral symmetry) to terete, succulent leaves (radial symmetry). We tested whether this shift had a corresponding change in internal Kranz anatomy and tissue patterning. METHODS: For a comprehensive characterization of C4 photosynthesis across early ontogeny in T. simplex, structural and ultrastructural anatomical properties and localization patterns, activities, and immunoblotting of key C4 photosynthetic enzymes were compared in mesophyll and bundle sheath tissues in cotyledons and leaves. KEY RESULTS: Cotyledons and leaves possessed different types of Kranz anatomy (atriplicoid type and a "Tetraena" variant of the kochioid type, respectively), reflecting the change in leaf morphology. In bundle sheath cells, key differences in ultrastructural features included increased organelle numbers and chloroplast thylakoid stacking. C4 enzymes had strict tissue-specific localization patterns within bundle sheath and mesophyll cells in both cotyledons and leaves. The decarboxylase NAD-ME maintained the highest activity, increasing from cotyledons to leaves. This classified T. simplex as fully C4 across ontogeny and a strictly NAD-ME biochemical subtype. CONCLUSIONS: Tetraena simplex cotyledons and leaves showed differences in Kranz type, with associated progression in ultrastructural features, and differing activities/expression levels of C4 enzymes. Furthermore, leaves characterized a new "Tetraena" variation of the kochioid Kranz anatomy.

Sexual homomorphism in dioecious trees: extensive tests fail to detect sexual dimorphism in Populus †.

The evolution of sexual dimorphism and expansion of sex chromosomes are both driven through sexual conflict, arising from differing fitness optima between males and females. Here, we pair work in poplar (Populus) describing one of the smallest sex-determining regions known thus far in complex eukaryotes (~100 kbp) with comprehensive tests for sexual dimorphism using >1300 individuals from two Populus species and assessing 96 non-reproductive functional traits. Against expectation, we found sexual homomorphism (no non-reproductive trait differences between the sexes), suggesting that gender is functionally neutral with respect to non-reproductive features that affect plant survival and fitness. Combined with a small sex-determining region, we infer that sexual conflict may be effectively stymied or non-existent within these taxa. Both sexual homomorphism and the small sex-determining region occur against a background of strong environmental selection and local adaptation in Populus. This presents a powerful hypothesis for the evolution of dioecious species. Here, we suggest that environmental selection may be sufficient to suppress and stymy sexual conflict if it acts orthogonal to sexual selection, thereby placing limitations on the evolution of sexual dimorphism and genomic expansion of sex chromosomes.

Trait convergence and diversification arising from a complex evolutionary history in Hawaiian species of Scaevola.

Species variation in functional traits may reflect diversification relating to convergence and/or divergence depending on environmental pressures and phylogenetic history. We tested trait-environment relationships and their basis in finer-scale evolutionary processes among nine extant Hawaiian species of Scaevola L. (Goodeniaceae), a taxon with a complex history of three independent colonizations by different phylogenetic lineages, parallel ecological specialization, and homoploid hybridization events in Hawai'i. Using a wild population for each species, we evaluated traits related to plant function (morphology, leaf and wood anatomy, nutrient and carbon isotope composition). Hawaiian Scaevola species were distributed across coastal, dry forest and wet forest environments; multivariate environmental analysis using abiotic and biotic factors further showed that species from distantly related lineages inhabited similar environments. Many traits correlated with environment (based on the multivariate environmental analysis), considering both distantly related species and more closely related species. Scaevola species within shared habitats generally showed trait convergence across distantly related lineages, particularly among wet forest species. Furthermore, trait diversification through divergence was extensive among closely related Scaevola species that radiated into novel environments, especially in plant morphology and traits affecting water relations. Homoploid hybrid-origin species were "intermediate" compared to their ancestral parent species, and possessed trait combinations relevant for their current habitat. The diversity in functional traits reflected strong influences of both ecology and evolutionary history in native Hawaiian Scaevola species, and trait correspondence with environment was due to the combination of multiple processes within the taxon: trait pre-adaptation and filtering, evolutionary convergence, divergence, and hybridization.

Evolutionary Quantitative Genomics of Populus trichocarpa.

Forest trees generally show high levels of local adaptation and efforts focusing on understanding adaptation to climate will be crucial for species survival and management. Here, we address fundamental questions regarding the molecular basis of adaptation in undomesticated forest tree populations to past climatic environments by employing an integrative quantitative genetics and landscape genomics approach. Using this comprehensive approach, we studied the molecular basis of climate adaptation in 433 Populus trichocarpa (black cottonwood) genotypes originating across western North America. Variation in 74 field-assessed traits (growth, ecophysiology, phenology, leaf stomata, wood, and disease resistance) was investigated for signatures of selection (comparing QST-FST) using clustering of individuals by climate of origin (temperature and precipitation). 29,354 SNPs were investigated employing three different outlier detection methods and marker-inferred relatedness was estimated to obtain the narrow-sense estimate of population differentiation in wild populations. In addition, we compared our results with previously assessed selection of candidate SNPs using the 25 topographical units (drainages) across the P. trichocarpa sampling range as population groupings. Narrow-sense QST for 53% of distinct field traits was significantly divergent from expectations of neutrality (indicating adaptive trait variation); 2,855 SNPs showed signals of diversifying selection and of these, 118 SNPs (within 81 genes) were associated with adaptive traits (based on significant QST). Many SNPs were putatively pleiotropic for functionally uncorrelated adaptive traits, such as autumn phenology, height, and disease resistance. Evolutionary quantitative genomics in P. trichocarpa provides an enhanced understanding regarding the molecular basis of climate-driven selection in forest trees and we highlight that important loci underlying adaptive trait variation also show relationship to climate of origin. We consider our approach the most comprehensive, as it uncovers the molecular mechanisms of adaptation using multiple methods and tests. We also provide a detailed outline of the required analyses for studying adaptation to the environment in a population genomics context to better understand the species' potential adaptive capacity to future climatic scenarios.

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.

Testing a 'genes-to-ecosystems' approach to understanding aquatic-terrestrial linkages.

A 'genes-to-ecosystems' approach has been proposed as a novel avenue for integrating the consequences of intraspecific genetic variation with the underlying genetic architecture of a species to shed light on the relationships among hierarchies of ecological organization (genes → individuals → communities → ecosystems). However, attempts to identify genes with major effect on the structure of communities and/or ecosystem processes have been limited and a comprehensive test of this approach has yet to emerge. Here, we present an interdisciplinary field study that integrated a common garden containing different genotypes of a dominant, riparian tree, Populus trichocarpa, and aquatic mesocosms to determine how intraspecific variation in leaf litter alters both terrestrial and aquatic communities and ecosystem functioning. Moreover, we incorporate data from extensive trait screening and genome-wide association studies estimating the heritability and genes associated with litter characteristics. We found that tree genotypes varied considerably in the quality and production of leaf litter, which contributed to variation in phytoplankton abundances, as well as nutrient dynamics and light availability in aquatic mesocosms. These 'after-life' effects of litter from different genotypes were comparable to the responses of terrestrial communities associated with the living foliage. We found that multiple litter traits corresponding with aquatic community and ecosystem responses differed in their heritability. Moreover, the underlying genetic architecture of these traits was complex, and many genes contributed only a small proportion to phenotypic variation. Our results provide further evidence that genetic variation is a key component of aquatic-terrestrial linkages, but challenge the ability to predict community or ecosystem responses based on the actions of one or a few genes.

Landscape genomics of Populus trichocarpa: the role of hybridization, limited gene flow, and natural selection in shaping patterns of population structure.

Populus trichocarpa is an ecologically important tree across western North America. We used a large population sample of 498 accessions over a wide geographical area genotyped with a 34K Populus SNP array to quantify geographical patterns of genetic variation in this species (landscape genomics). We present evidence that three processes contribute to the observed patterns: (1) introgression from the sister species P. balsamifera, (2) isolation by distance (IBD), and (3) natural selection. Introgression was detected only at the margins of the species' distribution. IBD was significant across the sampled area as a whole, but no evidence of restricted gene flow was detected in a core of drainages from southern British Columbia (BC). We identified a large number of FST outliers. Gene Ontology analyses revealed that FST outliers are overrepresented in genes involved in circadian rhythm and response to red/far-red light when the entire dataset is considered, whereas in southern BC heat response genes are overrepresented. We also identified strong correlations between geoclimate variables and allele frequencies at FST outlier loci that provide clues regarding the selective pressures acting at these loci.

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.

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.

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.

Genome-wide association mapping for wood characteristics in Populus identifies an array of candidate single nucleotide polymorphisms.

Establishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits. Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29,233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 0.05) associated with 16 wood chemistry/ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation. Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Significant involvement of PEP-CK in carbon assimilation of C4 eudicots.

BACKGROUND AND AIMS: C4 eudicot species are classified into biochemical sub-types of C4 photosynthesis based on the principal decarboxylating enzyme. Two sub-types are known, NADP-malic enzyme (ME) and NAD-ME; however, evidence for the occurrence or involvement of the third sub-type (phosphoenolpyruvate carboxykinase; PEP-CK) is emerging. In this study, the presence and activity of PEP-CK in C4 eudicot species of Trianthema and Zaleya (Sesuvioideae, Aizoaceae) is clarified through analysis of key anatomical features and C4 photosynthetic enzymes. METHODS: Three C4 species (T. portulacastrum, T. sheilae and Z. pentandra) were examined with light and transmission electron microscopy for leaf structural properties. Activities and immunolocalizations of C4 enzymes were measured for biochemical characteristics. KEY RESULTS: Leaves of each species possess atriplicoid-type Kranz anatomy, but differ in ultrastructural features. Bundle sheath organelles are centripetal in T. portulacastrum and Z. pentandra, and centrifugal in T. sheilae. Bundle sheath chloroplasts in T. portulacastrum are almost agranal, whereas mesophyll counterparts have grana. Both T. sheilae and Z. pentandra are similar, where bundle sheath chloroplasts contain well-developed grana while mesophyll chloroplasts are grana deficient. Cell wall thickness is significantly greater in T. sheilae than in the other species. Biochemically, T. portulacastrum is NADP-ME, while T. sheilae and Z. pentandra are NAD-ME. Both T. portulacastrum and Z. pentandra exhibit considerable PEP-CK activity, and immunolocalization studies show dense and specific compartmentation of PEP-CK in these species, consistent with high PEP-CK enzyme activity. CONCLUSIONS: Involvement of PEP-CK in C4 NADP-ME T. portulacastrum and NAD-ME Z. petandra occurs irrespective of biochemical sub-type, or the position of bundle sheath chloroplasts. Ultrastructural traits, including numbers of bundle sheath peroxisomes and mesophyll chloroplasts, and degree of grana development in bundle sheath chloroplasts, coincide more directly with PEP-CK recruitment. Discovery of high PEP-CK activity in C4 Sesuvioideae species offers a unique opportunity for evaluating PEP-CK expression and suggests the possibility that PEP-CK recruitment may exist elsewhere in C4 eudicots.

Seasonality and phenology alter functional leaf traits.

In plant ecophysiology, functional leaf traits are generally not assessed in relation to phenological phase of the canopy. Leaf traits measured in deciduous perennial species are known to vary between spring and summer seasons, but there is a knowledge gap relating to the late-summer phase marked by growth cessation and bud set occurring well before fall leaf senescence. The effects of phenology on canopy physiology were tested using a common garden of over 2,000 black cottonwood (Populus trichocarpa) individuals originating from a wide geographical range (44-60ºN). Annual phenological events and 12 leaf-based functional trait measurements were collected spanning the entire summer season prior to, and following, bud set. Patterns of seasonal trait change emerged by synchronizing trees using their date of bud set. In particular, photosynthetic, mass, and N-based traits increased substantially following bud set. Most traits were significantly different between pre-bud set and post-bud set phase trees, with many traits showing at least 25% alteration in mean value. Post-bud set, both the significance and direction of trait-trait relationships could be modified, with many relating directly to changes in leaf mass. In Populus, these dynamics in leaf traits throughout the summer season reflected a shift in whole plant physiology, but occurred long before the onset of leaf senescence. The marked shifts in measured trait values following bud set underscores the necessity to include phenology in trait-based ecological studies or large-scale phenotyping efforts, both at the local level and larger geographical scale.

Developmentally based scaling of leaf venation architecture explains global ecological patterns.

Leaf size and venation show remarkable diversity across dicotyledons, and are key determinants of plant adaptation in ecosystems past and present. Here we present global scaling relationships of venation traits with leaf size. Across a new database for 485 globally distributed species, larger leaves had major veins of larger diameter, but lower length per leaf area, whereas minor vein traits were independent of leaf size. These scaling relationships allow estimation of intact leaf size from fragments, to improve hindcasting of past climate and biodiversity from fossil remains. The vein scaling relationships can be explained by a uniquely synthetic model for leaf anatomy and development derived from published data for numerous species. Vein scaling relationships can explain the global biogeographical trend for smaller leaves in drier areas, the greater construction cost of larger leaves and the ability of angiosperms to develop larger and more densely vascularised lamina to outcompete earlier-evolved plant lineages.