Mapping morphological shape as a high-dimensional functional curve

Multivariate genome-wide association study of leaf shape in a Populus deltoides and P. simonii F1 pedigree

Leaf morphology exhibits tremendous diversity between and within species, and is likely related to adaptation to environmental factors. Most poplar species are of great economic and ecological values and their leaf morphology can be a good predictor for wood productivity and environment adaptation. It is important to understand the genetic mechanism behind variation in leaf shape. Although some initial efforts have been made to identify quantitative trait loci (QTLs) for poplar leaf traits, more effort needs to be expended to unravel the polygenic architecture of the complex traits of leaf shape. Here, we performed a genome-wide association analysis (GWAS) of poplar leaf shape traits in a randomized complete block design with clones from F1 hybrids of Populus deltoides and Populus simonii. A total of 35 SNPs were identified as significantly associated with the multiple traits of a moderate number of regular polar radii between the leaf centroid and its edge points, which could represent the leaf shape, based on a multivariate linear mixed model. In contrast, the univariate linear mixed model was applied as single leaf traits for GWAS, leading to genomic inflation; thus, no significant SNPs were detected for leaf length, measures of leaf width, leaf area, or the ratio of leaf length to leaf width under genomic control. Investigation of the candidate genes showed that most flanking regions of the significant leaf shape-associated SNPs harbored genes that were related to leaf growth and development and to the regulation of leaf morphology. The combined use of the traditional experimental design and the multivariate linear mixed model could greatly improve the power in GWAS because the multiple trait data from a large number of individuals with replicates of clones were incorporated into the statistical model. The results of this study will enhance the understanding of the genetic mechanism of leaf shape variation in Populus. In addition, a moderate number of regular leaf polar radii can largely represent the leaf shape and can be used for GWAS of such a complicated trait in Populus, instead of the higher-dimensional regular radius data that were previously considered to well represent leaf shape.

Heterophylly Quantitative Trait Loci Respond to Salt Stress in the Desert Tree Populus euphratica

Heterophylly, or leaf morphological changes along plant shoot axes, is an important indicator of plant eco-adaptation to heterogeneous microenvironments. Despite extensive studies on the genetic control of leaf shape, the genetic architecture of heterophylly remains elusive. To identify genes related to heterophylly and their associations with plant saline tolerance, we conducted a leaf shape mapping experiment using leaves from a natural population of Populus euphratica. We included 106 genotypes grown under salt stress and salt-free (control) conditions using clonal seedling replicates. We developed a shape tracking method to monitor and analyze the leaf shape using principal component (PC) analysis. PC1 explained 42.18% of the shape variation, indicating that shape variation is mainly determined by the leaf length. Using leaf length along shoot axes as a dynamic trait, we implemented a functional mapping-assisted genome-wide association study (GWAS) for heterophylly. We identified 171 and 134 significant quantitative trait loci (QTLs) in control and stressed plants, respectively, which were annotated as candidate genes for stress resistance, auxin, shape, and disease resistance. Functions of the stress resistance genes ABSCISIC ACIS-INSENSITIVE 5-like (ABI5), WRKY72, and MAPK3 were found to be related to many tolerance responses. The detection of AUXIN RESPONSE FACTOR17-LIKE (ARF17) suggests a balance between auxin-regulated leaf growth and stress resistance within the genome, which led to the development of heterophylly via evolution. Differentially expressed genes between control and stressed plants included several factors with similar functions affecting stress-mediated heterophylly, such as the stress-related genes ABC transporter C family member 2 (ABCC2) and ABC transporter F family member (ABCF), and the stomata-regulating and reactive oxygen species (ROS) signaling gene RESPIRATORY BURST OXIDASE HOMOLOG (RBOH). A comparison of the genetic architecture of control and salt-stressed plants revealed a potential link between heterophylly and saline tolerance in P. euphratica, which will provide new avenues for research on saline resistance-related genetic mechanisms.