Morphometrics Reveals Complex and Heritable Apple Leaf Shapes

Apple (Malus spp.) is a widely grown and valuable fruit crop. Leaf shape is important for flowering in apple and may also be an early indicator for other agriculturally valuable traits. We examined 9,000 leaves from 869 unique apple accessions using linear measurements and comprehensive morphometric techniques. We identified allometric variation as the result of differing length-to-width aspect ratios between accessions and species of apple. The allometric variation was due to variation in the width of the leaf blade, not the length. Aspect ratio was highly correlated with the first principal component (PC1) of morphometric variation quantified using elliptical Fourier descriptors (EFDs) and persistent homology (PH). While the primary source of variation was aspect ratio, subsequent PCs corresponded to complex shape variation not captured by linear measurements. After linking the morphometric information with over 122,000 genome-wide single nucleotide polymorphisms (SNPs), we found high SNP heritability values even at later PCs, indicating that comprehensive morphometrics can capture complex, heritable phenotypes. Thus, techniques such as EFDs and PH are capturing heritable biological variation that would be missed using linear measurements alone.

[1]  Zhiwu Zhang,et al.  Mixed linear model approach adapted for genome-wide association studies , 2010, Nature Genetics.

[2]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[3]  Pedro M. Valero-Mora,et al.  ggplot2: Elegant Graphics for Data Analysis , 2010 .

[4]  R. Varshney,et al.  Marker‐Assisted Backcrossing to Introgress Resistance to Fusarium Wilt Race 1 and Ascochyta Blight in C 214, an Elite Cultivar of Chickpea , 2014 .

[5]  Eugene W. Myers,et al.  Basic local alignment search tool. Journal of Molecular Biology , 1990 .

[6]  S. Myles,et al.  LinkImputeR: user-guided genotype calling and imputation for non-model organisms , 2017, BMC Genomics.

[7]  Sara Taskinen,et al.  smatr 3– an R package for estimation and inference about allometric lines , 2012 .

[8]  Peter J. Bradbury,et al.  Genome-wide association study of leaf architecture in the maize nested association mapping population , 2011, Nature Genetics.

[9]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[10]  E.E. Pissaloux,et al.  Image Processing , 1994, Proceedings. Second Euromicro Workshop on Parallel and Distributed Processing.

[11]  Robert J. Elshire,et al.  A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species , 2011, PloS one.

[12]  H. Kläring,et al.  ALLOMETRY TO ESTIMATE LEAF AREA OF TOMATO , 2001 .

[13]  T. Toldam-Andersen,et al.  Modification of leaf morphology and anatomy as a consequence of columnar architecture in domestic apple (Malus × domestica Borkh.) trees , 2013 .

[14]  Kevin P. Smith,et al.  PopVar: A Genome‐Wide Procedure for Predicting Genetic Variance and Correlated Response in Biparental Breeding Populations , 2015 .

[15]  J. Gurevitch,et al.  Sources of variation in leaf shape among two populations of Achillea lanulosa. , 1992, Genetics.

[16]  D. Garrick,et al.  Quantitative evaluation of apple (Malus × domestica Borkh.) fruit shape by principal component analysis of Fourier descriptors , 2000, Euphytica.

[17]  H. Iwata,et al.  SHAPE: a computer program package for quantitative evaluation of biological shapes based on elliptic Fourier descriptors. , 2002, The Journal of heredity.

[18]  Morphometric analysis of Passiflora leaves: the relationship between landmarks of the vasculature and elliptical Fourier descriptors of the blade , 2016 .

[19]  H. Tsukaya,et al.  Two independent and polarized processes of cell elongation regulate leaf blade expansion in Arabidopsis thaliana (L.) Heynh. , 1996, Development.

[20]  Allison J. Miller,et al.  Persistent homology demarcates a leaf morphospace , 2017, bioRxiv.

[21]  S. Myles,et al.  Genome to Phenome Mapping in Apple Using Historical Data , 2016, The plant genome.

[22]  Robert J. Elshire,et al.  TASSEL-GBS: A High Capacity Genotyping by Sequencing Analysis Pipeline , 2014, PloS one.

[23]  Allison J. Miller,et al.  Digital Morphometrics of Two North American Grapevines (Vitis: Vitaceae) Quantifies Leaf Variation between Species, within Species, and among Individuals , 2017, Front. Plant Sci..

[24]  Mao Li,et al.  Persistent homology and the branching topologies of plants. , 2017, American journal of botany.

[25]  M. Province,et al.  Avoiding the high Bonferroni penalty in genome‐wide association studies , 2009, Genetic epidemiology.

[26]  C. Topp,et al.  Persistent homology: a tool to universally measure plant morphologies across organs and scales , 2017, bioRxiv.

[27]  P. Visscher,et al.  GCTA: a tool for genome-wide complex trait analysis. , 2011, American journal of human genetics.

[28]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[29]  Michael F. Covington,et al.  A Modern Ampelography: A Genetic Basis for Leaf Shape and Venation Patterning in Grape1[C][W][OPEN] , 2013, Plant Physiology.

[30]  M. Lenhard,et al.  The INDETERMINATE DOMAIN Protein BROAD LEAF1 Limits Barley Leaf Width by Restricting Lateral Proliferation , 2016, Current Biology.

[31]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[32]  Michael F. Covington,et al.  A Quantitative Genetic Basis for Leaf Morphology in a Set of Precisely Defined Tomato Introgression Lines[C][W][OPEN] , 2013, Plant Cell.

[33]  Eden R Martin,et al.  A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms , 2008, Genetic epidemiology.

[34]  Characterizing Water Use Efficiency and Water Deficit Responses in Apple (Malus ×domestica Borkh. and Malus sieversii Ledeb.) M. Roem. , 2011 .

[35]  Christopher R. Cooney,et al.  Mega-evolutionary dynamics of the adaptive radiation of birds , 2017, Nature.

[36]  G. P. Lohmann Eigenshape analysis of microfossils: A general morphometric procedure for describing changes in shape , 1983 .

[37]  Don C. Jones,et al.  Modifications to a LATE MERISTEM IDENTITY1 gene are responsible for the major leaf shapes of Upland cotton (Gossypium hirsutum L.) , 2016, Proceedings of the National Academy of Sciences.

[38]  Ping Zheng,et al.  The Genome Database for Rosaceae (GDR): year 10 update , 2013, Nucleic Acids Res..

[39]  Roger E Bumgarner,et al.  The genome of the domesticated apple (Malus × domestica Borkh.) , 2010, Nature Genetics.

[40]  John W. Palmer,et al.  Effects of crop load on fruiting and gas-exchange characteristics of 'Braeburn'/M.26 apple trees at full canopy. , 2000 .

[41]  Achim Tresch,et al.  Heterochrony underpins natural variation in Cardamine hirsuta leaf form , 2015, Proceedings of the National Academy of Sciences.

[42]  C. Gessler,et al.  Mapping quantitative physiological traits in apple (Malus × domestica Borkh.) , 2003, Plant Molecular Biology.

[43]  R. Hellens,et al.  Asthma and the environment. , 1997, BMC Genomics.

[44]  R. Hellens,et al.  Multiple Repeats of a Promoter Segment Causes Transcription Factor Autoregulation in Red Apples[W] , 2009, The Plant Cell Online.

[45]  D. Bates,et al.  Fitting Linear Mixed-Effects Models Using lme4 , 2014, 1406.5823.

[46]  Julin N. Maloof,et al.  The Developmental Trajectory of Leaflet Morphology in Wild Tomato Species[C][W][OA] , 2012, Plant Physiology.

[47]  K. Olsen,et al.  Fruit Quality Traits Have Played Critical Roles in Domestication of the Apple , 2014 .

[48]  Takashi Akagi,et al.  Quantitative characterization of fruit shape and its differentiation pattern in diverse persimmon ( Diospyros kaki ) cultivars , 2018 .

[49]  D. Ferree,et al.  Flowering, pollination and fruit set and development. , 2003 .

[50]  Edward S. Buckler,et al.  TASSEL: software for association mapping of complex traits in diverse samples , 2007, Bioinform..

[51]  Russell V. Lenth,et al.  Least-Squares Means: The R Package lsmeans , 2016 .

[52]  Alan N. Lakso,et al.  The Relationship Between Leaf Area and Light Interception by Spur and Extension Shoot Leaves and Apple Orchard Productivity , 2000 .

[53]  Daniel Koenig,et al.  Natural Variation in Leaf Morphology Results from Mutation of a Novel KNOX Gene , 2008, Current Biology.

[54]  Cédric Gaucherel,et al.  Momocs: Outline Analysis Using R , 2014 .

[55]  Michael D. Abràmoff,et al.  Image processing with ImageJ , 2004 .

[56]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[57]  Xianzhong Feng,et al.  Evolution through genetically controlled allometry space. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[58]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..