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

Wenguo Yang; Dan Yao; Hainan Wu; Wei Zhao; Yuhua Chen; Chunfa Tong

doi:10.1371/journal.pone.0259278

Abstract

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 F₁ 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.

Citation: Yang W, Yao D, Wu H, Zhao W, Chen Y, Tong C (2021) Multivariate genome-wide association study of leaf shape in a Populus deltoides and P. simonii F₁ pedigree. PLoS ONE 16(10): e0259278. https://doi.org/10.1371/journal.pone.0259278

Editor: Karthikeyan Adhimoolam, Jeju National University, REPUBLIC OF KOREA

Received: May 3, 2021; Accepted: October 16, 2021; Published: October 28, 2021

Copyright: © 2021 Yang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The RADseq data is available in the SRA database at NCBI (http://www.ncbi.nlm.nih.gov/Traces/sra) with the accession numbers listed in S1 Table. Other relevant data are within the paper and its Supporting Information files.

Funding: Funding for this research was provided by the National Natural Science Foundation of China (No. 31870654 and 31270706) and the Priority Academic Program Development of the Jiangsu Higher Education Institutions (PAPD).The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Leaves are the most fundamental photosynthetic organs in plants; they are responsible for absorbing solar energy to generate power for plant growth and thus provide food for many species on earth [1, 2]. Leaf morphology exhibits tremendous diversity between or within species, such as the broad leaves of poplars and needle leaves of conifers. Leaf size and shape are evolutionarily adapted to environmental changes in response to water and light stress [3, 4], making it possible to reconstruct the paleoclimate [5, 6]. In model systems, several genes and networks have been identified to affect initial leaf development and pattern formation [2, 7, 8] as well as leaf length and width [9–11] using the mutagenesis screening method. Moreover, quantitative trait loci have also been detected for leaf morphological traits in species such as tomato [12], Arabidopsis [13], Brassica [14], maize [15], barley [16], and Populus [17, 18]. Despite advances made in these studies, the identified genes or loci may only cover a portion of the leaf morphological variation observed in nature because the variation is considered to be under polygenic control [11, 19].

The genus Populus (2n = 38) is an ecologically and economically important tree with a wide distribution in diverse environments of the Northern Hemisphere [20, 21]. The genus, comprising approximately 30 species, was grouped into six sections (i.e., Abaso, Aigeiros, Leucoides, Populus, Tacamahaca, and Turanga) according to morphological parameters [22]. Most species have several attractive biological characteristics, such as fast growth and easy asexual reproduction, so they are of particular interest to forest breeders for developing new cultivars to meet the needs of pulp, paper, lumber, and biofuels industries. Several studies have shown that leaf traits are highly related to growth and habitat and can be used as predictors of productivity and determinant factors in phylogenetic relationships [11, 23, 24]. Therefore, persistent efforts have been made to dissect the genetic mechanism of morphological traits in the genus. In the 1990s, Wu et al. [25] first conducted QTL mapping of leaf morphology in F₂ hybrids of P. trichocarpa × P. deltoides, with up to 3 QTLs identified for each trait, leaf area and the ratio of leaf length to width, at four crown positions. Recently, Mckown et al. [26] found 6 and 5 SNPs significantly associated with leaf length and width, respectively, in a GWAS on unrelated wild accessions of Populus trichocarpa. Drost et al. [11] detected 2 QTLs for lamina length, 2 for width, and 5 for their ratio in a pseudobackcross population of P. trichocarpa and P. deltoides. More recently, Chhetri et al. [17, 18] performed a GWAS of many traits with different genotypes from natural populations of P. trichocarpa; they did not detect any significant SNPs for single leaf traits, including leaf length, width, perimeter, area, and aspect ratio, but the detected up to 9 SNPs for leaf morphology multitraits. In contrast, Fu et al. [27] precisely described leaf shape with radii from the centroid to the contour at regular intervals and performed a marker-trait association analysis of principal components of the high-dimensional radius data, leading to several QTLs identified for leaf shape in a natural population of P. szechuanica var tibetica. They further modeled the leaf contour of a QTL genotype as a dynamic trajectory and identified a few significant QTLs for the variation of leaf shape in the same population [28]. These studies could be considered an initial stage for unravelling the genetic mechanism behind leaf size and shape in Populus. More powerful strategies, such as the utilization of novel statistical methods and generation of more accurate phenotype and genotype data, should be taken into account to ensure the accuracy and precision of such a tough task.

Herein, we report a genome-wide association study (GWAS) of leaf size and shape with a randomized complete block design (RCBD), which was established using clones from an F₁ hybrid population of P. deltoides and P. simonii [29] belonging to the sections Aigeiros and Tacamahaca, respectively. The leaves of the female P. deltoides are broad, while those of the male P. simonii parent are narrow. This sharp contrast led to diverse leaf shapes in their progeny (Fig 1). The leaf traits were digitally derived from scanned images of leaves, including the classical indices of leaf length, width, and area, as well as the high-dimensional data of regular ordered radii between the leaf centroid and edge points, as described in Fu et al. [27]. Single nucleotide polymorphism (SNP) genotypes of each clone were generated by mapping the paired-end (PE) reads from restriction site-associated DNA sequencing (RADseq) to the reference genome of P. trichocarpa [21]. With these SNP data, we applied a linear mixed model (LMM) to conduct GWAS for multiple or single leaf traits using the R package EMMREML (https://cran.rproject.org/web/packages/EMMREML). Consequently, we identified many more SNPs significantly associated with leaf shape than those detected in previous studies. Furthermore, candidate genes associated with these SNPs were investigated to show that most flanking regions of these significant SNPs harbored genes that were related to leaf growth and development and to the regulation of leaf morphology. The results enhanced our understanding of the genetic mechanism of leaf shape variation in Populus. We demonstrated that the combined use of the traditional experimental design and the multivariate linear mixed model (mvLMM) could greatly improve the power of GWAS for leaf shape. Additionally, the multivariate data of 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. This is contrary to the expectation that the high-dimensional regular radius data could well represent the leaf shape for GWAS, as indicated by Fu et al. [27].

Download:

Fig 1. Leaf shape variation among parents and progeny of P. deltoides × P. simonii.

https://doi.org/10.1371/journal.pone.0259278.g001

Materials and methods

Genetic experimental design and measurement of leaf traits

To obtain repeated phenotypic data for accurate QTL analysis, we established an RCBD in the spring of 2017 with clones from an F₁ hybrid population of the female P. deltoides and the male P. simonii, which was produced in the spring of 2009 to 2011 [29, 30]. The design consisted of 234 clones with 3 blocks, 6 cuttings for each clone per plot within a block, and a 50 × 60 cm spacing in Xiashu Forest Farm of Nanjing Forestry University, Jurong County, Jiangsu Province, China (32.1224°N, 119.2155°E). The sixth most apical mature leaf of each individual was sampled in mid-July 2017 by placing it into a paper bag and then scanned using a Hewlett-Packard Scanjet G2410 A4 flatbed scanner at a resolution of 300 dpi. The A4-sized images were saved as bmp files, and leaf size and shape traits were analyzed with the R package LeafShape (https://github.com/tongchf/LeafShape). These traits included leaf area (A), length (L), maximum width (W), and widths at one-third length (W1/3), half length (W1/2), and two-thirds length (W2/3) from the base, as well as 360 regular polar radii (RD360) between the leaf centroid and edge points, as shown in Fig 2. As our primary analytical step, we applied multivariate statistical methods to these leaf traits with SAS 9.3 software (SAS Institute, Cary, USA), including canonical correlation analysis and principal component analysis.

Download:

Fig 2.

Workflow of the R package LeafShape for extracting leaf shape traits in poplar hybrids: (A) a fresh leaf; (B) original (blue) and position-adjusted (red) edge points extracted from the scanned image of a leaf; (C) the red vertical line indicates the leaf length (L), the red horizontal line indicates the maximum leaf width (W), and the blue horizontal lines represent the leaf widths at one-third length (W1/3), half length (W1/2), and two-thirds length (W2/3); (D) 360 regular polar radii between the leaf centroid and edge points across −π to π.

https://doi.org/10.1371/journal.pone.0259278.g002

SNP genotyping

Since 2013, more than four hundred individuals in the poplar hybrid population have been sequenced with RADseq technology in several batches [29, 30]. The 163 clones from the RCBD experiment and their two parents were sequenced previously, and their RADseq paired-end (PE) reads were deposited in the NCBI SRA database with the accession numbers listed in S1 Table. The PE reads for each individual were filtered using the NGS QC toolkit with default parameters [31], and the resulting high-quality (HQ) read data were used for calling SNP genotypes. The SNP calling procedures were largely the same as those described in Mousavi et al. [30]. In brief, the PE reads of each clone or parent were first aligned to the reference sequence of P. trichocarpa with the software BWA v0.7.17 [32]. Second, the resulting SAM (sequence alignment/map) file was converted into a BAM (binary alignment/map) formatted file and then sorted and indexed with SAMtools v1.9 [33]. Third, the sorted BAM file was analyzed to generate a VCF file using BCFtools v1.9 software [33]. Finally, the VCF file was filtered to obtain SNP genotypes of each individual such that a heterozygous allele had a read coverage depth (DP) of at least 3 and the quality of each SNP genotype was greater than 30.

Because the 163 clones were from an F₁ hybrid population, the SNPs were expected to segregate in several different patterns, such as aa×ab and ab×ab, due to the characteristics of outbred forest species [34, 35]. We classified the SNPs into subsets according to the segregation patterns and kept those that did not seriously deviate from Mendelian segregation ratios (p > 0.01), including 1:1, 1:2:1, and 1:1:1:1. In addition, if an SNP had more than 10% missing genotypes across the clones, it was removed from the dataset. To obtain independent SNP markers and linkage disequilibrium (LD) blocks, we performed the LD-based SNP pruning procedure for all the SNPs using PLINK v1.07 software with a window size of 25 SNPs, a step size of 2 SNPs, and an r² threshold of 0.7 [36].

Statistical methods for GWAS

Since the poplar leaf is generally symmetrical, the polar radii on the right side largely contain the full information of all the radius data on both sides (Fig 2D). We used the multiple traits of different numbers of regular radii across −π/2 to π/2 to find SNPs associated with leaf shape, which was implemented with the mvLMM as follows: (1) where y_ijkl is the lth polar radius of the kth tree leaf of the jth clone in the ith block; μ_l is the overall mean of the lth polar radius; B_il is the effect of the ith block; M_jl is the genotype effect of the jth clone at any tested SNP; G_jl is the polygenic background effect of the jth clone; and e_ijkl is the residual effect. It is assumed that B_il and M_jl are fixed effects, while G_jl and e_ijkl are random effects with , and cov(G_jl, e_ijkl) = 0. In matrix form, model (1) can be written as (2) where Y is the n×t matrix whose (I, j)th element is the jth trait value of the ith individual, i.e., the ith row of Y is the multiple trait data for the ith individual; X is an n×p known design matrix of fixed effects, including the overall mean, block effects, and individual genotype effects at the tested SNP site; B is the p×t matrix of fixed effect coefficients; G is the c×t matrix whose (i, j)th element is the background random additive genetic effect of the ith clone and the jth trait with Vec(G)~N_c×t(0, V_G⊗A), where Vec denotes the matrix vectorization function [37], A is the additive relationship matrix for the c clones and V_G is the additive genetic variance matrix for the t traits; Z is the n×c coefficient matrix corresponding to the matrix G; E is the random residual matrix with Vec(E)~N_n×t(0, V_E⊗I_n). Based on the assumptions above, the covariance matrix of Vec(Y) can be derived as (3) Because the clones in our experiment belonged to a full-sib family and their parents were unrelated and not inbred, the kinship coefficient of any two clones is 0.25 in theory [38], leading to the relationship matrix A with ones on the diagonal and 0.5 elsewhere.

To test whether any single SNP was significantly associated with the leaf traits, an F statistic was used under the null hypothesis MVec(B) = 0 for a full-rank q×pt matrix M, as (4) with q numerator degrees of freedom and t(n−p) denominator degrees of freedom [39]. The p-values from the F statistics (4) for each SNP are prone to genomic inflation [40, 41]. It is necessary to calculate the genomic inflation factor (λ_GC) to evaluate the inflation level. When there was no genomic inflation, the p-value threshold for testing significant SNPs was determined based on Bonferroni correction at the 0.05 significance level.

The proportion of phenotype variance explained (PVE) by a single SNP was calculated as (5) where RSS₀ and RSS are the residual sums of squares under the null hypothesis model and the full model (2), respectively [42].

As a comparison with the regular radius data, we also used mvLMM (2) to perform GWAS for the multiple traits of L, W, W1/3, W1/2, and W2/3 (LWs) (Fig 2C). For a single leaf trait such as length, width, and area, the GWAS was conducted with univariate LMM, which can be derived by simplifying the multivariate model (2) and is expressed as (6) where y is a vector of trait values for n individuals; X is a design matrix of fixed effects; β is a vector of fixed effects; g is a vector of random genetic effects for each clone with ; Z is the coefficient matrix corresponding to the random vector g; e is the random residual vector with . Moreover, we calculated the narrow heritability of a single trait as without incorporating the fixed effects of SNPs in model (6).

To calculate the restricted maximum likelihood (REML) estimates of genetic parameters, we applied the function emmremlMultivariate for the multivariate model (2) and emmreml for the univariate model (6) in the R package EMMREML (https://cran.r-project.org/web/packages/EMMREML). After the genetic parameters were calculated, the p-value for testing each SNP was calculated according to the F statistic, as in Eq (4).

Investigation of candidate genes

In our previous study [43], the average LD block length was estimated to be ~650 bp in the same hybrid population, which is so short that it could not be properly used as a downstream or upstream range for investigating candidate genes for the significant SNPs. Alternatively, we took the strategy as described in Slaten et al [44]. In brief, we considered candidate genes that contain significant SNPs or are within a LD block harboring significant SNPs. If a significant SNP is within an intergenic region and does not form a LD block, both the closest downstream and upstream genes are considered as candidates. Because no annotation information for leaf shape is available in the gene annotation database of P. trichocarpa in Phytozome (v4.1; https://genome.jgi.doe.gov), the coding sequences (CDSs) of these genes were obtained for further annotating. We first performed BLAST searches with their CDSs against the nonredundant protein database [45, 46] and then mapped all BLAST hits to Gene Ontology (GO) terms based on ID mapping information from http://ftp.pir.georgetown.edu/databases/idmapping/idmapping.tb.gz. The descriptions of the blast hits and GO terms were saved in an Excel file in which we could search which genes were possibly related to leaf shape.

Results

Leaf trait data

We successfully obtained the leaf trait data for a total of 2,244 individual trees belonging to 163 clones in the RCBD (S2 Table). Some plots had missed samples due to the damage from pest, disease, poor rooting ability, or other unknown reasons. To validate these measurements, we measured 100 randomly chosen leaves with ImageJ [47] and LeafShape software separately. The average relative differences in the leaf length, width, and area values measured from the two software programs were 1.45 (±0.99)%, 4.76 (±0.76)%, and 5.05 (±0.92)%, respectively, indicating that the two measurements from both methods were largely consistent (S3 Table). The descriptive statistics for the traits L, W, W1/3, W1/2, W2/3, A, and the L/W ratio are presented in Table 1, including the mean, standard deviation, range, and coefficient of variation (CV). The CVs for the leaf length and different leaf widths were similar, ranging from 20.79% for L to 25.14% for W2/3, while the CV for leaf area reached a maximum value of 42.25% and the CV for the length/width ratio had a minimum value of 10.12%. The histograms showed that these leaf traits basically followed a normal distribution (S1 Fig). Furthermore, the heritabilities of leaf length and different leaf widths as well as leaf area were similar (40~50%), but the heritability of the length/width ratio was much higher at 64.74% (Table 1). In addition, correlation analysis showed that the leaf length, measures of leaf width, and leaf area were significantly positively correlated (p < 0.01) with each other, with most coefficients over 0.90; the minimum coefficient value was 0.8137 between L and W2/3 (S4 Table). However, the L/W ratio was significantly negatively correlated with each of the leaf length, different leaf widths, and leaf area traits, with coefficients between -0.6160 and -0.2389. Finally, analysis of variance showed that the effects of each leaf trait, L, W, W1/3, W1/2, W2/3, and A, were significantly different among blocks and clones (S5 Table).

Download:

Table 1. Variation in leaf length, different leaf widths, leaf area, and the ratio of length/width in the F₁ progeny of Populus deltoides × Populus simonii based on a randomized complete block design.

https://doi.org/10.1371/journal.pone.0259278.t001

The 360 polar radii (Fig 2D) of all leaves were obtained with the R package LeafShape as a full dataset denoted as RD360. We also extracted 5 reduced datasets denoted as RD06, RD09, RD11, RD16, and RD61 that contained 6, 9, 11, 16, and 61 regular polar radii of each leaf on the right side from −π/2 to π/2, respectively; these datasets were expected to represent the leaf shape characteristics despite the lower dimensionality of the data due to leaf symmetry. Canonical correlation analysis showed that each of the commonly measured leaf traits, such as length, width, and area, was highly correlated with the polar radii in the RD360, RD61, RD16, RD11, RD09, and RD06 datasets, with a correlation coefficient value of over 0.98 and a p-value less than 0.0001 (S6 Table). Additionally, the multiple traits of leaf length and different leaf widths (LWs: L, W, W1/3, W1/2, and W2/3) were extremely significantly correlated with the 6 radius datasets, with the first canonical correlation coefficient greater than 0.9996 (p<0.0001). Moreover, the principal component analysis of the polar radius traits revealed that the proportion of total variance for the first principal component was at least 95.59% for the 6 radius datasets. The first principal component for each radius dataset was highly correlated with leaf length, different leaf widths, and leaf area, with coefficients greater than 0.90 and p-values less than 0.0001 (S7 Table).

SNP genotype data

A total of 33,086 SNPs across the 163 clones were obtained by mapping their high-quality PE reads separately to the reference genome of P. trichocarpa (v4.1; https://genome.jgi.doe.gov). For a SNP genotype of each clone at each SNP site, the heterozygous allele was required to have a coverage depth of at least 3 reads, whereas the coverage depth for a homozygous allele was at least 5. Furthermore, the quality of each genotype needed a Phred score of at least 30, and the missing genotype rate at each SNP was set to less than 10%. All the SNPs were categorized into five segregation types, aa×ab, aa×ac, ab×aa, ab×ab, and ab×cc (Table 2). The majority of SNPs segregated at a ratio of 1:1 (p > 0.01) with aa×ab and ab×aa types.

Download:

Table 2. Summary of SNPs obtained across the 163 clones based on a randomized complete block design.

https://doi.org/10.1371/journal.pone.0259278.t002

The LD analysis of these SNPs was performed with the software PLINK [36], resulting in 10,735 independent SNP markers and LD blocks. Therefore, the p-value threshold for significant SNPs in our genome-wide analyses was set to 0.05/10735 = 4.66E-6 (-log₁₀(p-value) = 5.33) based on the Bonferroni correction at the 0.05 significance level.

Significant SNPs associated with leaf traits

mvLMM (2) was applied to perform the GWAS for the multiple traits of the regular polar radius datasets RD06, RD09, RD11, RD16, and RD61 separately, as well as the multiple traits of LWs. The quantile-quantile plots of the p-values on base 10 logarithm scale showed that there existed different levels of genomic inflation, with inflation factors greater than 1 for datasets RD06, RD09, and LWs; less than 1 for datasets RD16 and RD61; and almost equal to 1 for dataset RD11 (Fig 3). Because the result from dataset RD11 showed good genomic control, we used this result to determine the significant SNPs associated with leaf shape. Consequently, a total of 35 SNPs were found to be significantly associated with the multiple traits of leaf shape under the p-value threshold of 4.66E-6, each explaining 0.18–0.32% of the phenotypic variance (Table 3). Fig A shows the Manhattan plot of the negative base 10 logarithm of the p-value against the corresponding SNP position. These significant SNPs were distributed on 8 chromosomes, with a few SNPs on chromosomes 4, 6, 10, 15, 16, and 18 but more on chromosome 14. There were 11 significant SNPs detected on chromosome 1, of which the first 10 were within a 3-Mb region. More surprisingly, chromosome 14 harbored the most (16) significant SNPs, which could be divided into five regions according to the position where the negative logarithm of the p-value changed from decreasing to increasing (Fig 4B).

Download:

Fig 3. Quantile-quantile plots of observed p-values versus expected p-values on a base 10 logarithm scale with genomic inflation factors for GWAS of the different multiple traits representing the poplar leaf shape.

LWs indicate the multiple traits of leaf length and 4 different leaf widths. RD06, RD09, RD11, RD16, and RD61 indicate the multiple traits of 6, 9, 11, 16, and 61 regular polar radii between the leaf centroid and edge points across −π/2 to π/2, respectively.

https://doi.org/10.1371/journal.pone.0259278.g003

Download:

Fig 4. Manhattan plot of the association analysis for the 11 regular polar radii between the leaf centroid and edge points from -π/2 to π/2.

(A) The plot shows the 19 chromosomes of the reference genome of P. trichocarpa. The horizontal dashed line indicates the genome-wide significance threshold of 5.33, which is a base 10 logarithm of the p-value based on the Bonferroni correction at the 0.05 significance level. (B) The significant SNPs on chromosome 14 were divided into five regions roughly according to the position where the negative logarithm of the p-value changes from decreasing to increasing.

https://doi.org/10.1371/journal.pone.0259278.g004

Download:

Table 3. Summary of the significant SNPs associated with the leaf shape represented by the multiple trait dataset RD11.

https://doi.org/10.1371/journal.pone.0259278.t003

Moreover, LMM (6) was applied to detect associations with each single leaf trait, such as leaf length, width and area. The results showed that the genomic inflation factors for these traits ranged from 1.774 for W2/3 to 2.461 for the L/W ratio (Fig 5). After genomic control, we found that there were no significant SNPs associated with any single trait under the p-value threshold based on Bonferroni correction (Fig 5). However, without genomic control, various numbers of significant SNPs were found for these traits: only one significant SNP each was detected for W, W1/3, and A; 6, 8, and 33 SNPs were detected for W1/2, W2/3, and the L/W ratio, respectively; and no significant SNPs were detected for L (S8 and S9 Tables; S2 Fig). The SNP at position 10669990 on chromosome 10 was a common SNP significantly associated with the four different leaf widths, and the significant SNPs for the width at two-thirds length shared all but one significant SNP with the width at half length. In addition, the most significant SNP sites or regions for the ratio of leaf length to leaf width were consistent with those for the multiple traits of leaf length and four different leaf widths, except for the two significant SNPs on chromosomes 2 and 13.

Download:

Fig 5.

Manhattan and quantile-quantile (QQ) plots of the association analyses for each univariate trait, L (A), W (B), W31 (C), W21 (D), W32 (E), and A (F), and the ratio of L to W (G) across the 19 chromosomes of the reference genome of P. trichocarpa. The left panel presents the Manhattan plots under genomic control, while the right panel shows the corresponding QQ plots before (blue) and after (green) genomic control. The horizontal dashed line indicates the genome-wide significance threshold of 5.33, which is a base 10 logarithm of the p-value based on the Bonferroni correction at the 0.05 significance level.

https://doi.org/10.1371/journal.pone.0259278.g005

Candidate genes affecting leaf shape

The candidate genes of the significant SNPs for the multiple traits of the 11-dimensional regular polar radii data were annotated with the nonredundant protein database at the NCBI and GO databases (S10 Table). One significant SNP region on chromosome 1 and five on chromosome 14 were found to harbor a total of 40 candidate genes functionally related to leaf shape (Table 3). However, the rest 9 significant SNPs on chromosome 1, 4, 6, 10, 15, 16, and 18 had no candidate genes that have descriptions directly related to leaf shape, possibly due to the reason that each of them did not form a LD block with other SNPs and thus had at most two candidate genes. We found that there are 8 candidate genes in 5 significant SNP regions, which directly affect leaf growth and development, with descriptions such as “leaf development” and “regulation of leaf morphogenesis”. It was also noticed that there are 6 candidate genes on chromosomes 1 and 14 related to the hormone auxin, which plays important roles in initial leaf formation, lamina margin elaboration, and leaf vasculature patterning [48–51]. Moreover, 12 candidate genes were found to belong to MYB gene family, which was previously reported to be involved in leaf development in Arabidopsis [52] and maize [53]. Furthermore, 2 candidate genes on chromosomes 1 and 14 are related to TCP genes, which were found to be involved in leaf development and morphology in Arabidopsis [54, 55]. In addition, 14 candidate genes were related to light responses or photosynthesis in 5 significant SNP regions distributed on chromosomes 1 and14; these genes are involved in activities such as response to light intensity, light harvesting, and photosynthesis. Undoubtedly, these genes play important roles in leaf development and pattern formation.

Discussion

Leaf size and shape are the most important traits during the development and growth of Populus. Understanding the genetic mechanism of these traits is of great interest to many poplar breeders. In the present study, we successfully detected dozens of SNPs significantly associated with the multiple traits of the 11-dimensional regular leaf polar radii in a randomized complete block test with clones from the F₁ hybrids of P. deltoides and P. simonii. Multiple traits could be considered to represent the leaf shape because the regular polar radii on the right side largely reflect the two-dimensional pattern of the leaf. Compared with previous studies for identifying QTLs or SNPs associated with leaf shape in Populus (see Introduction), we were able to identify many more QTLs or significant SNP regions. One of the main reasons for the powerful ability to identify the associated SNPs may be attributed to the use of the RCBD in the current GWAS. This kind of test design provided replicates of clones not only at the block level but also at the plot level, allowing thousands of individuals to be used for the association analysis. From a statistical perspective, the repeated phenotype data for each genotype that originated from a single seed can control for the spatial effects in the field and reduce systematic errors, hence improving the accuracy and power of GWAS. In contrast, in previous GWAS or QTL mapping studies on poplar leaf traits, phenotype data were measured from single plants with different genotypes in natural populations or full-sib families, possibly limiting QTL detection power.

Another advantage of our association analysis strategy may be due to incorporating the multiple traits of leaf polar radii into the mvLMM for GWAS. Although mvLMMs have become increasingly important in GWAS because of their power gain over univariate analysis, the computation of genetic parameter estimates is nontrivial [56]. We successfully implemented the parameter and statistical calculations with the flexible R package EMMREML by adding or modifying some codes. Consequently, the mvLMM helped identify many more significant SNPs associated with leaf traits without genomic deflation. In contrast, after genomic control, the univariate LMM did not have the ability to detect any significant SNPs for any single trait, such as leaf length and width (Fig 5). Even if genomic deflation was permitted, we could see that fewer than 10 significant SNPs were detected for the single traits W, W1/3, W1/2, W2/3, and A, whereas no significant SNPs were detected for L (S8 Table). However, in such cases, the number of significant SNPs dramatically increased to 33 for the ratio of leaf length to width but was still less than the number detected based on the multiple traits of the 11-dimensional regular leaf polar radii dataset (S9 Table). The fact that more SNPs were detected for the ratio of leaf length to leaf width than for the other single traits may largely be due to much higher heritability of this trait (Table 1). This phenomenon can also be found in a previous study [11], where the authors identified 2 QTLs for leaf length and 2 for width but 5 for the ratio of the two traits.

Although our association analysis of the multiple traits based on the mvLMM was able to identify many more significant SNPs, it seems that the PVE of each SNP was much lower, ranging from 0.18 to 0.32% (Table 3). An intuitive explanation for this result is that leaf shape is possibly controlled by many genes with small effects, conforming to the infinitesimal model [57]. This explanation could further confirm that our strategy for GWAS in the current study is powerful for detecting such small-effect genes. This phenomenon may be the main reason why previous studies had a lower power for locating QTLs for single leaf traits, with only a few detected, although the PVEs of the QTLs were apparently larger than those estimated in this study [11, 25, 28]. However, the PVEs of SNPs or QTLs cannot be compared directly because they are calculated based on not only different population structures but also different statistical models. Even in the same study using the same statistical model, the PVE may or may not consistently increase or decrease with the corresponding statistical value for determining the significance of the hypothesis test. This is because the estimates of the environmental variance vary for different SNPs or QTLs, possibly leading to inconsistencies between the PVEs and statistical values. This phenomenon can be commonly found in the literature. For example, in Drost et al. [11], the first QTL for lamina length had a PVE value of 6.31% with a LOD value of 3.14, while the PVE of the second QTL was 8.10% with a lower LOD value of 2.68. In addition to these factors, the most important consideration is how to calculate the PVE based on a statistical model. For most fixed linear models with uncorrelated phenotype data, the R² statistic is generally used to measure the PVE in QTL mapping studies or GWAS. However, for mixed linear models, such measurements are not well established [58]. Here, we calculated the R² statistic as Eq (5) based on the weighted residual sum of squares [59].

It is worth emphasizing that the 11-dimensional multivariate data of the regular leaf polar radii can largely represent the poplar leaf shape and can be applied in association analyses with SNPs for such traits that are difficult to measure. Naturally, it was believed that the higher the dimensionality of the radius data between the leaf centroid and edge points is, the better the characteristics of leaf shape can be represented (Fig 2D). Fu et al. [27] first implemented such an idea by extracting 360 coordinates on leaf outlines from scanned images and performed a series of association analyses with the leaf shape [28, 60]. We also performed GWAS of leaf shape with different dimensions of the radius data (e.g., RD61, RD16), which were extracted by our own R package (https://github.com/tongchf/LeafShape) because Fu et al. did not provide public software for the task. However, our results showed that for the higher dimensional data (i.e., RD61 and RD16), genomic deflation existed with λ_GC≤0.820, while for the lower dimensional data (i.e., RD09 and RD06), genomic inflation existed with λ_GC≥1.120 (Fig 3). In contrast, the RD11 data presented a balanced result between genomic inflation and deflation, exhibiting the best performance regarding genomic control in the GWAS with different dimensional data of the regular leaf polar radii.

Compared with previous studies for poplar leaf shape, we found that there were a few overlapping regions (<5 Mb) containing significant SNPs or QTLs. S11 Table lists those significant SNPs or QTLs associated with leaf shape in the current study and in four recent studies [17, 18, 26, 61], excluding those previous QTL studies in which no physical QTL position information was available [11, 25, 27, 28]. The results in the previous studies for single leaf traits such leaf length and width were not considered because we thought that the leaf shape could not be described by a single leaf parameter. We found that there were 7 significant SNPs detected in our study very close (<5 Mb) to one or more SNPs identified in previous studies, of which 5 were consistent with Xia et al. [61], 4 with Chhetri et al. [17], 1 with McKown et al. [26], and 1 with Chhetri et al. [18]. In contrast, 5 overlapping regions were found between the four previous studies. It is interesting to find that 3 regions on chromosome 4, 6, and 8 were coincidentally detected for leaf shape in three studies. Although our GWAS findings have more consistent SNPs with the previous results, most SNPs identified in the current and previous studies did not share an overlapping region. This result may be due to many reasons, but one of the main reasons is that different methods were used to describe the complex trait of leaf shape in the GWAS or QTL studies. Drost et al. [11] described the leaf shape with the ratio of leaf length to width, while Chhetri et al. [17, 18] described it with the combination of leaf area (LA), leaf dry weight (LD), leaf length (LL) and leaf width (LW) or the combination of leaf aspect ratio (AR) and specific leaf area (SL). However, based on the method of Fu et al. [27, 28], we used high-dimensional regular polar radii data to describe the leaf shape.

Conclusion

The novel strategy for GWAS with direct integration of the traditional randomized complete block design and the multiple traits of regular leaf polar radii into the multivariate linear mixed model facilitated the identification of many more significant SNPs associated with leaf shape in Populus than previous studies have detected. Moreover, it was demonstrated that the multivariate linear mixed model was more powerful than the univariate linear mixed model in the association analyses for leaf traits such as leaf length, width, and area. Most flanking regions surrounding significant SNPs harbored potential candidate genes that were related to the growth and development of the poplar leaf. Our results enhance the understanding of the molecular mechanism underlying leaf morphological variation in Populus. In addition, the multivariate data from a moderate number of regular leaf polar radii could largely represent the leaf shape and exhibited better genomic control in the GWAS of poplar leaf shape.

Supporting information

S1 Fig.

Histograms with probability density curves (red) of normal distributions for each univariate trait of L (A), W (B), W31 (C), W21 (D), W32 (E), A (F), and the ratio of L to W (G) in the randomized complete block design derived from the F1 progeny of Populus deltoides × Populus simonii.

https://doi.org/10.1371/journal.pone.0259278.s001

(DOCX)

S2 Fig.

Manhattan plots of the association analyses without genomic control for each univariate trait of L (A), W (B), W31 (C), W21 (D), W32 (E), A (F), and the ratio of L to W (G) across the 19 chromosomes of the reference genome of P. trichocarpa. The horizontal dashed line indicates the genome-wide significant threshold of 5.33, a base 10 logarithm of p-value based on the Bonferroni correction at the 0.05 significant level.

https://doi.org/10.1371/journal.pone.0259278.s002

(DOCX)

S1 Table. The RADseq data information for the 2 parents and 163 progeny in the F1 hybrid population of Populus deltoides and Populus simonii.

https://doi.org/10.1371/journal.pone.0259278.s003

(DOCX)

S2 Table. The raw data of leaf length, different widths and area from a randomized complete block design in Populus.

https://doi.org/10.1371/journal.pone.0259278.s004

(XLSX)

S3 Table. Relative difference between leaf measurements with the two software ImageJ and LeafShape.

https://doi.org/10.1371/journal.pone.0259278.s005

(XLSX)

S4 Table. Correlation coefficients among the leaf traits of L, W, W31, W21, W32, A, and the ratio of L to W in the randomized complete block design derived from the F1 progeny of Populus deltoides × Populus simonii.

https://doi.org/10.1371/journal.pone.0259278.s006

(DOCX)

S5 Table. Analysis of variance for the leaf parameters of L, W, W1/3, W1/2, W2/3, and Area in the randomized complete block experiment derived from the F1 progeny of Populus deltoides × Populus simonii.

https://doi.org/10.1371/journal.pone.0259278.s007

(DOCX)

S6 Table. Canonical correlation coefficients among the leaf length, widths, area, the length/width ratio, and polar radii in the randomized complete block design derived from the F1 progeny of Populus deltoides × Populus simonii.

https://doi.org/10.1371/journal.pone.0259278.s008

(DOCX)

S7 Table. Correlation coefficients between the first principal component of different radius datasets and the leaf length, different widths, or area in the randomized complete block design derived from the F1 progeny of Populus deltoides × Populus simonii.

https://doi.org/10.1371/journal.pone.0259278.s009

(DOCX)

S8 Table. Summary of significant SNPs associated to each trait of the four different leaf widths and area without genomic control.

https://doi.org/10.1371/journal.pone.0259278.s010

(DOCX)

S9 Table. Summary of significant SNPs associated to the ratio of the leaf length to the maximum width without genomic control.

https://doi.org/10.1371/journal.pone.0259278.s011

(DOCX)

S10 Table. The annotation of candidate genes for the significant SNPs with the non-redundant protein database at NCBI and GO database.

https://doi.org/10.1371/journal.pone.0259278.s012

(XLSX)

S11 Table. Consistency between significant SNPs for poplar leaf shape identified in the current and previous studies.

The distance of two close SNPs between two different studies is presented in brackets.

https://doi.org/10.1371/journal.pone.0259278.s013

(XLSX)

Acknowledgments

We thank Professor Huogen Li in Nanjing Forestry University for his great help in establishing the randomized complete block design.

References

1. Field CB, Behrenfeld MJ, Randerson JT, Falkowski P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science. 1998;281(5374):237–40. pmid:9657713.
- View Article
- PubMed/NCBI
- Google Scholar
2. Fleming AJ. The control of leaf development. New Phytol. 2005;166(1):9–20. pmid:15760347.
- View Article
- PubMed/NCBI
- Google Scholar
3. Xu F, Guo W, Xu W, Wei Y, Wang R. Leaf morphology correlates with water and light availability: What consequences for simple and compound leaves? Progress in Natural Science. 2009;19(12):1789–98. WOS:000271902400017.
- View Article
- Google Scholar
4. Liao F, Peng J, Chen R. LeafletAnalyzer, an automated software for quantifying, comparing and classifying blade and serration features of compound leaves during development, and among induced mutants and natural variants in the legume Medicago truncatula. Front Plant Sci. 2017;8. WOS:000402386300001. pmid:28620405
- View Article
- PubMed/NCBI
- Google Scholar
5. Krieger JD, Guralnick RP, Smith DM. Generating empirically determined, continuous measures of leaf shape for paleoclimate reconstruction. PALAIOS. 2007;22(2):212–9. WOS:000248060700010.
- View Article
- Google Scholar
6. Peppe DJ, Royer DL, Cariglino B, Oliver SY, Newman S, Leight E, et al. Sensitivity of leaf size and shape to climate: global patterns and paleoclimatic applications. New Phytologist. 2011;190(3):724–39. pmid:21294735
- View Article
- PubMed/NCBI
- Google Scholar
7. Byrne ME. Networks in leaf development. Curr Opin Plant Biol. 2005;8(1):59–66. pmid:15653401.
- View Article
- PubMed/NCBI
- Google Scholar
8. Barkoulas M, Hay A, Kougioumoutzi E, Tsiantis M. A developmental framework for dissected leaf formation in the Arabidopsis relative Cardamine hirsuta. Nat Genet. 2008;40(9):1136–41. pmid:19165928.
- View Article
- PubMed/NCBI
- Google Scholar
9. Tsuge T, Tsukaya H, Uchimiya H. Two independent and polarized processes of cell elongation regulate leaf blade expansion in Arabidopsis thaliana (L.) Heynh. Development. 1996;122(5):1589–600. pmid:8625845.
- View Article
- PubMed/NCBI
- Google Scholar
10. Narita NN, Moore S, Horiguchi G, Kubo M, Demura T, Fukuda H, et al. Overexpression of a novel small peptide ROTUNDIFOLIA4 decreases cell proliferation and alters leaf shape in Arabidopsis thaliana. Plant J. 2004;38(4):699–713. pmid:15125775.
- View Article
- PubMed/NCBI
- Google Scholar
11. Drost DR, Puranik S, Novaes E, Novaes CR, Dervinis C, Gailing O, et al. Genetical genomics of Populus leaf shape variation. BMC Plant Biol. 2015;15:166. pmid:26122556; PubMed Central PMCID: PMC4486686.
- View Article
- PubMed/NCBI
- Google Scholar
12. Holtan HEE, Hake S. Quantitative trait locus analysis of leaf dissection in tomato using Lycopersicon pennellii segmental introgression lines. Genetics. 2003;165(3):1541–50. WOS:000187459100048. pmid:14668401
- View Article
- PubMed/NCBI
- Google Scholar
13. Juenger T, Perez-Perez JM, Bernal S, Micol JL. Quantitative trait loci mapping of floral and leaf morphology traits in Arabidopsis thaliana: evidence for modular genetic architecture. Evolution & Development. 2005;7(3):259–71. WOS:000228690100010. pmid:15876198
- View Article
- PubMed/NCBI
- Google Scholar
14. Li F, Kitashiba H, Inaba K, Nishio T. A Brassica rapa linkage map of EST-based SNP markers for identification of candidate genes controlling flowering time and leaf morphological traits. DNA Research. 2009;16(6):311–23. WOS:000272832700001. pmid:19884167
- View Article
- PubMed/NCBI
- Google Scholar
15. Fu Y, Xu G, Chen H, Wang X, Chen Q, Huang C, et al. QTL mapping for leaf morphology traits in a large maize-teosinte population. Molecular Breeding. 2019;39(7). WOS:000473170500002.
- View Article
- Google Scholar
16. Du B, Liu L, Wang Q, Sun G, Ren X, Li C, et al. Identification of QTL underlying the leaf length and area of different leaves in barley. Sci Rep. 2019;9(1):4431. pmid:30872632; PubMed Central PMCID: PMC6418291.
- View Article
- PubMed/NCBI
- Google Scholar
17. Chhetri HB, Macaya-Sanz D, Kainer D, Biswal AK, Evans LM, Chen JG, et al. Multitrait genome-wide association analysis of Populus trichocarpa identifies key polymorphisms controlling morphological and physiological traits. New Phytologist. 2019;223(1):293–309. WOS:308432139100027.
- View Article
- Google Scholar
18. Chhetri HB, Furches A, Macaya-Sanz D, Walker AR, Kainer D, Jones P, et al. Genome-Wide Association Study of Wood Anatomical and Morphological Traits inPopulus trichocarpa. Front Plant Sci. 2020;11. WOS:000574412000001. pmid:33013968
- View Article
- PubMed/NCBI
- Google Scholar
19. Bylesjo M, Segura V, Soolanayakanahally RY, Rae AM, Trygg J, Gustafsson P, et al. LAMINA: a tool for rapid quantification of leaf size and shape parameters. BMC Plant Biol. 2008;8:82. pmid:18647399; PubMed Central PMCID: PMC2500018.
- View Article
- PubMed/NCBI
- Google Scholar
20. Geraldes A, Pang J, Thiessen N, Cezard T, Moore R, Zhao Y, et al. SNP discovery in black cottonwood (Populus trichocarpa) by population transcriptome resequencing. Mol Ecol Resour. 2011;11(Suppl 1):81–92. pmid:21429165.
- View Article
- PubMed/NCBI
- Google Scholar
21. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006; 313:1596–604. pmid:16973872
- View Article
- PubMed/NCBI
- Google Scholar
22. Eckenwalder JE. Systematics and evolution of Populus. In: Stettler RF, Bradshaw HD, Heilman PE, Hinckley TM, editors. Biology of Populus and its implications for management and conservation. Ottawa: NRC Research Press, National Council of Canada; 1996. p. p. 7–32.
23. Stettler R, Bradshaw HD, Heilman PE, Hinckley TM. Biology of Populus and its implications for management and conservation. Ottawa: NRC Research Press; 1996.
24. Marron N, Ceulemans R. Genetic variation of leaf traits related to productivity in a Populus deltoides x Populus nigra family. Canadian Journal of Forest Research. 2006;36(2):390–400. WOS:000237199900012.
- View Article
- Google Scholar
25. Wu R, Bradshaw H, Stettler R. Molecular genetics of growth and development in Populus (Salicaceae). v. mapping quantitative trait loci affecting leaf variation. American Journal of Botany. 1997;84(2):143–53. MEDLINE:pmid:21712194.
- View Article
- PubMed/NCBI
- Google Scholar
26. McKown AD, Klapste J, Guy RD, Geraldes A, Porth I, Hannemann J, et al. Genome-wide association implicates numerous genes underlying ecological trait variation in natural populations of Populus trichocarpa. New Phytologist. 2014;203(2):535–53. WOS:000337639800019. pmid:24750093
- View Article
- PubMed/NCBI
- Google Scholar
27. Fu G, Bo W, Pang X, Wang Z, Chen L, Song Y, et al. Mapping shape quantitative trait loci using a radius-centroid-contour model. Heredity. 2013;110(6):511–9. WOS:000319112000002. pmid:23572125
- View Article
- PubMed/NCBI
- Google Scholar
28. Fu GF, Huang M, Bo WH, Hao H, Wu RL. Mapping morphological shape as a high-dimensional functional curve. Briefings in Bioinformatics. 2018;19(3):461–71. WOS:000432676200009. pmid:28062411
- View Article
- PubMed/NCBI
- Google Scholar
29. Tong CF, Li HG, Wang Y, Li XR, Ou JJ, Wang DY, et al. Construction of high-density linkage maps of Populus deltoides × P. simonii using restriction-site associated DNA sequencing. PLoS One. 2016;11(3):e0150692. pmid:26964097
- View Article
- PubMed/NCBI
- Google Scholar
30. Mousavi M, Tong C, Liu F, Tao S, Wu J, Li H, et al. De novo SNP discovery and genetic linkage mapping in poplar using restriction site associated DNA and whole-genome sequencing technologies. Bmc Genomics. 2016;17:656. pmid:27538483; PubMed Central PMCID: PMC4991039.
- View Article
- PubMed/NCBI
- Google Scholar
31. Patel RK, Jain M. NGS QC Toolkit: A toolkit for quality control of next generation sequencing data. PLoS One. 2012;7(2):e30619. pmid:22312429
- View Article
- PubMed/NCBI
- Google Scholar
32. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25:1754–60. pmid:19451168
- View Article
- PubMed/NCBI
- Google Scholar
33. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9. pmid:19505943
- View Article
- PubMed/NCBI
- Google Scholar
34. Maliepaard C, Jansen J, Van Ooijen JW. Linkage analysis in a full-sib family of an outbreeding plant species: overview and consequences for applications. Genet Res. 1997;70:237–50.
- View Article
- Google Scholar
35. Tong CF, Zhang B, Shi JS. A hidden Markov model approach to multilocus linkage analysis in a full-sib family. Tree Genet Genomes. 2010; 6(5): 651–62.
- View Article
- Google Scholar
36. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics. 2007;81(3):559–75. WOS:000249128200012. pmid:17701901
- View Article
- PubMed/NCBI
- Google Scholar
37. Searle SR, Casella G, McCulloch CE. Variance Components. New Jersey, USA: John Wiley & Sons, Inc.; 2006. 458 p.
38. Lynch M, Walsh B. Genetics and Analysis of Quantitative Traits. Sunderland, MA, USA: Sinauer Associates, Inc.; 1998. 143–5 p.
39. Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, Daly MJ, et al. Efficient control of population structure in model organism association mapping. Genetics. 2008;178(3):1709–23. pmid:18385116; PubMed Central PMCID: PMC2278096.
- View Article
- PubMed/NCBI
- Google Scholar
40. Devlin B, Roeder K. Genomic control for association studies. Biometrics. 1999;55(4):997–1004. MEDLINE:pmid:11315092.
- View Article
- PubMed/NCBI
- Google Scholar
41. van Iterson M, van Zwet EW, Heijmans BT, Consortium B. Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution. Genome Biology. 2017;18. WOS:000394827100001. pmid:28126036
- View Article
- PubMed/NCBI
- Google Scholar
42. Xu R. Measuring explained variation in linear mixed effects models. Stat Med. 2003;22:3527–41. pmid:14601017
- View Article
- PubMed/NCBI
- Google Scholar
43. Chen Y, Wu H, Yang W, Zhao W, Tong C. Multivariate linear mixed model enhanced the power of identifying genome-wide association to poplar tree heights in a randomized complete block design. G3. 2021;11(2). Epub 2021/02/20. pmid:33604666; PubMed Central PMCID: PMC8022933.
- View Article
- PubMed/NCBI
- Google Scholar
44. Slaten ML, Yobi A, Bagaza C, Chan YO, Shrestha V, Holden S, et al. mGWAS Uncovers Gln-Glucosinolate Seed-Specific Interaction and its Role in Metabolic Homeostasis. Plant Physiol. 2020;183(2):483–500. Epub 2020/04/23. pmid:32317360; PubMed Central PMCID: PMC7271782.
- View Article
- PubMed/NCBI
- Google Scholar
45. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. Journal of Molecular Biology. 1990;215(3):403–10. pmid:2231712.
- View Article
- PubMed/NCBI
- Google Scholar
46. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–5. pmid:17130148; PubMed Central PMCID: PMC1716718.
- View Article
- PubMed/NCBI
- Google Scholar
47. Schindelin J, Rueden CT, Hiner MC, Eliceiri KW. The ImageJ ecosystem: An open platform for biomedical image analysis. Molecular Reproduction And Development. 2015;82(7–8):518–29. WOS:000358486700004. pmid:26153368
- View Article
- PubMed/NCBI
- Google Scholar
48. Reinhardt D, Pesce ER, Stieger P, Mandel T, Baltensperger K, Bennett M, et al. Regulation of phyllotaxis by polar auxin transport. Nature. 2003;426(6964):255–60. WOS:000186660800035. pmid:14628043
- View Article
- PubMed/NCBI
- Google Scholar
49. Scanlon MJ. The polar auxin transport inhibitor N-1-naphthylphthalamic acid disrupts leaf initiation, KNOX protein regulation, and formation of leaf margins in maize. Plant Physiology. 2003;133(2):597–605. WOS:000185974800022. pmid:14500790
- View Article
- PubMed/NCBI
- Google Scholar
50. Hay A, Tsiantis M. The genetic basis for differences in leaf form between Arabidopsis thaliana and its wild relative Cardamine hirsuta. Nature Genetics. 2006;38(8):942–7. WOS:000239325700026. pmid:16823378
- View Article
- PubMed/NCBI
- Google Scholar
51. Scarpella E, Marcos D, Friml J, Berleth T. Control of leaf vascular patterning by polar auxin transport. Genes & Development. 2006;20(8):1015–27. WOS:000236951500011. pmid:16618807
- View Article
- PubMed/NCBI
- Google Scholar
52. Sun Y, Zhou Q, Zhang W, Fu Y, Huang H. ASYMMETRIC LEAVES1, an Arabidopsis gene that is involved in the control of cell differentiation in leaves. Planta. 2002;214:694–702. pmid:11882937
- View Article
- PubMed/NCBI
- Google Scholar
53. Theodoris G, Inada N, Freeling M. Conservation and molecular dissection of ROUGH SHEATH2 and ASYMMETRIC LEAVES1 function in leaf development. Proc Natl Acad Sci USA. 2003;100:6837–42. pmid:12750468
- View Article
- PubMed/NCBI
- Google Scholar
54. Kieffer M, Master V, Waites R, Davies B. TCP14 and TCP15 affect internode length and leaf shape in Arabidopsis. Plant J. 2011;68(1):147–58. pmid:21668538; PubMed Central PMCID: PMC3229714.
- View Article
- PubMed/NCBI
- Google Scholar
55. Aguilar-Martinez JA, Sinha N. Analysis of the role of Arabidopsis class I TCP genes AtTCP7, AtTCP8, AtTCP22, and AtTCP23 in leaf development. Front Plant Sci. 2013;4:406. WOS:241371718200001.
- View Article
- Google Scholar
56. Zhou X, Stephens M. Efficient multivariate linear mixed model algorithms for genome-wide association studies. Nat Methods. 2014;11(4):407–9. pmid:24531419; PubMed Central PMCID: PMC4211878.
- View Article
- PubMed/NCBI
- Google Scholar
57. Hu Z, Wang Z, Xu S. An infinitesimal model for quantitative trait genomic value prediction. PLoS One. 2012;7(7):e41336. pmid:22815992; PubMed Central PMCID: PMC3399838.
- View Article
- PubMed/NCBI
- Google Scholar
58. Sun G, Zhu C, Kramer MH, Yang SS, Song W, Piepho HP, et al. Variation explained in mixed-model association mapping. Heredity. 2010;105(4):333–40. pmid:20145669.
- View Article
- PubMed/NCBI
- Google Scholar
59. Buse A. Goodness of fit in generalized least-squares estimation. The American Statistician. 1973;27:106–8.
- View Article
- Google Scholar
60. Fu GF, Saunders G, Stevens J. Holm multiple correction for large-scale gene-shape association mapping. BMC Genetics. 2014;15 (Suppl 1):S5. WOS:000345651700007. pmid:25079623
- View Article
- PubMed/NCBI
- Google Scholar
61. Xia WX, Xiao ZA, Cao P, Zhang Y, Du KB, Wang N. Construction of a high-density genetic map and its application for leaf shape QTL mapping in poplar. Planta. 2018;248(5):1173–85. WOS:000447030900009. pmid:30088086
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Field CB, Behrenfeld MJ, Randerson JT, Falkowski P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science. 1998;281(5374):237–40. pmid:9657713.
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Fleming AJ. The control of leaf development. New Phytol. 2005;166(1):9–20. pmid:15760347.
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Xu F, Guo W, Xu W, Wei Y, Wang R. Leaf morphology correlates with water and light availability: What consequences for simple and compound leaves? Progress in Natural Science. 2009;19(12):1789–98. WOS:000271902400017.
View Article
Google Scholar

[10] View Article

[11] Google Scholar

[ref4] 4. Liao F, Peng J, Chen R. LeafletAnalyzer, an automated software for quantifying, comparing and classifying blade and serration features of compound leaves during development, and among induced mutants and natural variants in the legume Medicago truncatula. Front Plant Sci. 2017;8. WOS:000402386300001. pmid:28620405
View Article
PubMed/NCBI
Google Scholar

[13] View Article

[14] PubMed/NCBI

[15] Google Scholar

[ref5] 5. Krieger JD, Guralnick RP, Smith DM. Generating empirically determined, continuous measures of leaf shape for paleoclimate reconstruction. PALAIOS. 2007;22(2):212–9. WOS:000248060700010.
View Article
Google Scholar

[17] View Article

[18] Google Scholar

[ref6] 6. Peppe DJ, Royer DL, Cariglino B, Oliver SY, Newman S, Leight E, et al. Sensitivity of leaf size and shape to climate: global patterns and paleoclimatic applications. New Phytologist. 2011;190(3):724–39. pmid:21294735
View Article
PubMed/NCBI
Google Scholar

[20] View Article

[21] PubMed/NCBI

[22] Google Scholar

[ref7] 7. Byrne ME. Networks in leaf development. Curr Opin Plant Biol. 2005;8(1):59–66. pmid:15653401.
View Article
PubMed/NCBI
Google Scholar

[24] View Article

[25] PubMed/NCBI

[26] Google Scholar

[ref8] 8. Barkoulas M, Hay A, Kougioumoutzi E, Tsiantis M. A developmental framework for dissected leaf formation in the Arabidopsis relative Cardamine hirsuta. Nat Genet. 2008;40(9):1136–41. pmid:19165928.
View Article
PubMed/NCBI
Google Scholar

[28] View Article

[29] PubMed/NCBI

[30] Google Scholar

[ref9] 9. Tsuge T, Tsukaya H, Uchimiya H. Two independent and polarized processes of cell elongation regulate leaf blade expansion in Arabidopsis thaliana (L.) Heynh. Development. 1996;122(5):1589–600. pmid:8625845.
View Article
PubMed/NCBI
Google Scholar

[32] View Article

[33] PubMed/NCBI

[34] Google Scholar

[ref10] 10. Narita NN, Moore S, Horiguchi G, Kubo M, Demura T, Fukuda H, et al. Overexpression of a novel small peptide ROTUNDIFOLIA4 decreases cell proliferation and alters leaf shape in Arabidopsis thaliana. Plant J. 2004;38(4):699–713. pmid:15125775.
View Article
PubMed/NCBI
Google Scholar

[36] View Article

[37] PubMed/NCBI

[38] Google Scholar

[ref11] 11. Drost DR, Puranik S, Novaes E, Novaes CR, Dervinis C, Gailing O, et al. Genetical genomics of Populus leaf shape variation. BMC Plant Biol. 2015;15:166. pmid:26122556; PubMed Central PMCID: PMC4486686.
View Article
PubMed/NCBI
Google Scholar

[40] View Article

[41] PubMed/NCBI

[42] Google Scholar

[ref12] 12. Holtan HEE, Hake S. Quantitative trait locus analysis of leaf dissection in tomato using Lycopersicon pennellii segmental introgression lines. Genetics. 2003;165(3):1541–50. WOS:000187459100048. pmid:14668401
View Article
PubMed/NCBI
Google Scholar

[44] View Article

[45] PubMed/NCBI

[46] Google Scholar

[ref13] 13. Juenger T, Perez-Perez JM, Bernal S, Micol JL. Quantitative trait loci mapping of floral and leaf morphology traits in Arabidopsis thaliana: evidence for modular genetic architecture. Evolution & Development. 2005;7(3):259–71. WOS:000228690100010. pmid:15876198
View Article
PubMed/NCBI
Google Scholar

[48] View Article

[49] PubMed/NCBI

[50] Google Scholar

[ref14] 14. Li F, Kitashiba H, Inaba K, Nishio T. A Brassica rapa linkage map of EST-based SNP markers for identification of candidate genes controlling flowering time and leaf morphological traits. DNA Research. 2009;16(6):311–23. WOS:000272832700001. pmid:19884167
View Article
PubMed/NCBI
Google Scholar

[52] View Article

[53] PubMed/NCBI

[54] Google Scholar

[ref15] 15. Fu Y, Xu G, Chen H, Wang X, Chen Q, Huang C, et al. QTL mapping for leaf morphology traits in a large maize-teosinte population. Molecular Breeding. 2019;39(7). WOS:000473170500002.
View Article
Google Scholar

[56] View Article

[57] Google Scholar

[ref16] 16. Du B, Liu L, Wang Q, Sun G, Ren X, Li C, et al. Identification of QTL underlying the leaf length and area of different leaves in barley. Sci Rep. 2019;9(1):4431. pmid:30872632; PubMed Central PMCID: PMC6418291.
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref17] 17. Chhetri HB, Macaya-Sanz D, Kainer D, Biswal AK, Evans LM, Chen JG, et al. Multitrait genome-wide association analysis of Populus trichocarpa identifies key polymorphisms controlling morphological and physiological traits. New Phytologist. 2019;223(1):293–309. WOS:308432139100027.
View Article
Google Scholar

[63] View Article

[64] Google Scholar

[ref18] 18. Chhetri HB, Furches A, Macaya-Sanz D, Walker AR, Kainer D, Jones P, et al. Genome-Wide Association Study of Wood Anatomical and Morphological Traits inPopulus trichocarpa. Front Plant Sci. 2020;11. WOS:000574412000001. pmid:33013968
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref19] 19. Bylesjo M, Segura V, Soolanayakanahally RY, Rae AM, Trygg J, Gustafsson P, et al. LAMINA: a tool for rapid quantification of leaf size and shape parameters. BMC Plant Biol. 2008;8:82. pmid:18647399; PubMed Central PMCID: PMC2500018.
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref20] 20. Geraldes A, Pang J, Thiessen N, Cezard T, Moore R, Zhao Y, et al. SNP discovery in black cottonwood (Populus trichocarpa) by population transcriptome resequencing. Mol Ecol Resour. 2011;11(Suppl 1):81–92. pmid:21429165.
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref21] 21. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006; 313:1596–604. pmid:16973872
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref22] 22. Eckenwalder JE. Systematics and evolution of Populus. In: Stettler RF, Bradshaw HD, Heilman PE, Hinckley TM, editors. Biology of Populus and its implications for management and conservation. Ottawa: NRC Research Press, National Council of Canada; 1996. p. p. 7–32.

[ref23] 23. Stettler R, Bradshaw HD, Heilman PE, Hinckley TM. Biology of Populus and its implications for management and conservation. Ottawa: NRC Research Press; 1996.

[ref24] 24. Marron N, Ceulemans R. Genetic variation of leaf traits related to productivity in a Populus deltoides x Populus nigra family. Canadian Journal of Forest Research. 2006;36(2):390–400. WOS:000237199900012.
View Article
Google Scholar

[84] View Article

[85] Google Scholar

[ref25] 25. Wu R, Bradshaw H, Stettler R. Molecular genetics of growth and development in Populus (Salicaceae). v. mapping quantitative trait loci affecting leaf variation. American Journal of Botany. 1997;84(2):143–53. MEDLINE:pmid:21712194.
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref26] 26. McKown AD, Klapste J, Guy RD, Geraldes A, Porth I, Hannemann J, et al. Genome-wide association implicates numerous genes underlying ecological trait variation in natural populations of Populus trichocarpa. New Phytologist. 2014;203(2):535–53. WOS:000337639800019. pmid:24750093
View Article
PubMed/NCBI
Google Scholar

[91] View Article

[92] PubMed/NCBI

[93] Google Scholar

[ref27] 27. Fu G, Bo W, Pang X, Wang Z, Chen L, Song Y, et al. Mapping shape quantitative trait loci using a radius-centroid-contour model. Heredity. 2013;110(6):511–9. WOS:000319112000002. pmid:23572125
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref28] 28. Fu GF, Huang M, Bo WH, Hao H, Wu RL. Mapping morphological shape as a high-dimensional functional curve. Briefings in Bioinformatics. 2018;19(3):461–71. WOS:000432676200009. pmid:28062411
View Article
PubMed/NCBI
Google Scholar

[99] View Article

[100] PubMed/NCBI

[101] Google Scholar

[ref29] 29. Tong CF, Li HG, Wang Y, Li XR, Ou JJ, Wang DY, et al. Construction of high-density linkage maps of Populus deltoides × P. simonii using restriction-site associated DNA sequencing. PLoS One. 2016;11(3):e0150692. pmid:26964097
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref30] 30. Mousavi M, Tong C, Liu F, Tao S, Wu J, Li H, et al. De novo SNP discovery and genetic linkage mapping in poplar using restriction site associated DNA and whole-genome sequencing technologies. Bmc Genomics. 2016;17:656. pmid:27538483; PubMed Central PMCID: PMC4991039.
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref31] 31. Patel RK, Jain M. NGS QC Toolkit: A toolkit for quality control of next generation sequencing data. PLoS One. 2012;7(2):e30619. pmid:22312429
View Article
PubMed/NCBI
Google Scholar

[111] View Article

[112] PubMed/NCBI

[113] Google Scholar

[ref32] 32. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25:1754–60. pmid:19451168
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref33] 33. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9. pmid:19505943
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref34] 34. Maliepaard C, Jansen J, Van Ooijen JW. Linkage analysis in a full-sib family of an outbreeding plant species: overview and consequences for applications. Genet Res. 1997;70:237–50.
View Article
Google Scholar

[123] View Article

[124] Google Scholar

[ref35] 35. Tong CF, Zhang B, Shi JS. A hidden Markov model approach to multilocus linkage analysis in a full-sib family. Tree Genet Genomes. 2010; 6(5): 651–62.
View Article
Google Scholar

[126] View Article

[127] Google Scholar

[ref36] 36. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics. 2007;81(3):559–75. WOS:000249128200012. pmid:17701901
View Article
PubMed/NCBI
Google Scholar

[129] View Article

[130] PubMed/NCBI

[131] Google Scholar

[ref37] 37. Searle SR, Casella G, McCulloch CE. Variance Components. New Jersey, USA: John Wiley & Sons, Inc.; 2006. 458 p.

[ref38] 38. Lynch M, Walsh B. Genetics and Analysis of Quantitative Traits. Sunderland, MA, USA: Sinauer Associates, Inc.; 1998. 143–5 p.

[ref39] 39. Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, Daly MJ, et al. Efficient control of population structure in model organism association mapping. Genetics. 2008;178(3):1709–23. pmid:18385116; PubMed Central PMCID: PMC2278096.
View Article
PubMed/NCBI
Google Scholar

[135] View Article

[136] PubMed/NCBI

[137] Google Scholar

[ref40] 40. Devlin B, Roeder K. Genomic control for association studies. Biometrics. 1999;55(4):997–1004. MEDLINE:pmid:11315092.
View Article
PubMed/NCBI
Google Scholar

[139] View Article

[140] PubMed/NCBI

[141] Google Scholar

[ref41] 41. van Iterson M, van Zwet EW, Heijmans BT, Consortium B. Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution. Genome Biology. 2017;18. WOS:000394827100001. pmid:28126036
View Article
PubMed/NCBI
Google Scholar

[143] View Article

[144] PubMed/NCBI

[145] Google Scholar

[ref42] 42. Xu R. Measuring explained variation in linear mixed effects models. Stat Med. 2003;22:3527–41. pmid:14601017
View Article
PubMed/NCBI
Google Scholar

[147] View Article

[148] PubMed/NCBI

[149] Google Scholar

[ref43] 43. Chen Y, Wu H, Yang W, Zhao W, Tong C. Multivariate linear mixed model enhanced the power of identifying genome-wide association to poplar tree heights in a randomized complete block design. G3. 2021;11(2). Epub 2021/02/20. pmid:33604666; PubMed Central PMCID: PMC8022933.
View Article
PubMed/NCBI
Google Scholar

[151] View Article

[152] PubMed/NCBI

[153] Google Scholar

[ref44] 44. Slaten ML, Yobi A, Bagaza C, Chan YO, Shrestha V, Holden S, et al. mGWAS Uncovers Gln-Glucosinolate Seed-Specific Interaction and its Role in Metabolic Homeostasis. Plant Physiol. 2020;183(2):483–500. Epub 2020/04/23. pmid:32317360; PubMed Central PMCID: PMC7271782.
View Article
PubMed/NCBI
Google Scholar

[155] View Article

[156] PubMed/NCBI

[157] Google Scholar

[ref45] 45. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. Journal of Molecular Biology. 1990;215(3):403–10. pmid:2231712.
View Article
PubMed/NCBI
Google Scholar

[159] View Article

[160] PubMed/NCBI

[161] Google Scholar

[ref46] 46. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–5. pmid:17130148; PubMed Central PMCID: PMC1716718.
View Article
PubMed/NCBI
Google Scholar

[163] View Article

[164] PubMed/NCBI

[165] Google Scholar

[ref47] 47. Schindelin J, Rueden CT, Hiner MC, Eliceiri KW. The ImageJ ecosystem: An open platform for biomedical image analysis. Molecular Reproduction And Development. 2015;82(7–8):518–29. WOS:000358486700004. pmid:26153368
View Article
PubMed/NCBI
Google Scholar

[167] View Article

[168] PubMed/NCBI

[169] Google Scholar

[ref48] 48. Reinhardt D, Pesce ER, Stieger P, Mandel T, Baltensperger K, Bennett M, et al. Regulation of phyllotaxis by polar auxin transport. Nature. 2003;426(6964):255–60. WOS:000186660800035. pmid:14628043
View Article
PubMed/NCBI
Google Scholar

[171] View Article

[172] PubMed/NCBI

[173] Google Scholar

[ref49] 49. Scanlon MJ. The polar auxin transport inhibitor N-1-naphthylphthalamic acid disrupts leaf initiation, KNOX protein regulation, and formation of leaf margins in maize. Plant Physiology. 2003;133(2):597–605. WOS:000185974800022. pmid:14500790
View Article
PubMed/NCBI
Google Scholar

[175] View Article

[176] PubMed/NCBI

[177] Google Scholar

[ref50] 50. Hay A, Tsiantis M. The genetic basis for differences in leaf form between Arabidopsis thaliana and its wild relative Cardamine hirsuta. Nature Genetics. 2006;38(8):942–7. WOS:000239325700026. pmid:16823378
View Article
PubMed/NCBI
Google Scholar

[179] View Article

[180] PubMed/NCBI

[181] Google Scholar

[ref51] 51. Scarpella E, Marcos D, Friml J, Berleth T. Control of leaf vascular patterning by polar auxin transport. Genes & Development. 2006;20(8):1015–27. WOS:000236951500011. pmid:16618807
View Article
PubMed/NCBI
Google Scholar

[183] View Article

[184] PubMed/NCBI

[185] Google Scholar

[ref52] 52. Sun Y, Zhou Q, Zhang W, Fu Y, Huang H. ASYMMETRIC LEAVES1, an Arabidopsis gene that is involved in the control of cell differentiation in leaves. Planta. 2002;214:694–702. pmid:11882937
View Article
PubMed/NCBI
Google Scholar

[187] View Article

[188] PubMed/NCBI

[189] Google Scholar

[ref53] 53. Theodoris G, Inada N, Freeling M. Conservation and molecular dissection of ROUGH SHEATH2 and ASYMMETRIC LEAVES1 function in leaf development. Proc Natl Acad Sci USA. 2003;100:6837–42. pmid:12750468
View Article
PubMed/NCBI
Google Scholar

[191] View Article

[192] PubMed/NCBI

[193] Google Scholar

[ref54] 54. Kieffer M, Master V, Waites R, Davies B. TCP14 and TCP15 affect internode length and leaf shape in Arabidopsis. Plant J. 2011;68(1):147–58. pmid:21668538; PubMed Central PMCID: PMC3229714.
View Article
PubMed/NCBI
Google Scholar

[195] View Article

[196] PubMed/NCBI

[197] Google Scholar

[ref55] 55. Aguilar-Martinez JA, Sinha N. Analysis of the role of Arabidopsis class I TCP genes AtTCP7, AtTCP8, AtTCP22, and AtTCP23 in leaf development. Front Plant Sci. 2013;4:406. WOS:241371718200001.
View Article
Google Scholar

[199] View Article

[200] Google Scholar

[ref56] 56. Zhou X, Stephens M. Efficient multivariate linear mixed model algorithms for genome-wide association studies. Nat Methods. 2014;11(4):407–9. pmid:24531419; PubMed Central PMCID: PMC4211878.
View Article
PubMed/NCBI
Google Scholar

[202] View Article

[203] PubMed/NCBI

[204] Google Scholar

[ref57] 57. Hu Z, Wang Z, Xu S. An infinitesimal model for quantitative trait genomic value prediction. PLoS One. 2012;7(7):e41336. pmid:22815992; PubMed Central PMCID: PMC3399838.
View Article
PubMed/NCBI
Google Scholar

[206] View Article

[207] PubMed/NCBI

[208] Google Scholar

[ref58] 58. Sun G, Zhu C, Kramer MH, Yang SS, Song W, Piepho HP, et al. Variation explained in mixed-model association mapping. Heredity. 2010;105(4):333–40. pmid:20145669.
View Article
PubMed/NCBI
Google Scholar

[210] View Article

[211] PubMed/NCBI

[212] Google Scholar

[ref59] 59. Buse A. Goodness of fit in generalized least-squares estimation. The American Statistician. 1973;27:106–8.
View Article
Google Scholar

[214] View Article

[215] Google Scholar

[ref60] 60. Fu GF, Saunders G, Stevens J. Holm multiple correction for large-scale gene-shape association mapping. BMC Genetics. 2014;15 (Suppl 1):S5. WOS:000345651700007. pmid:25079623
View Article
PubMed/NCBI
Google Scholar

[217] View Article

[218] PubMed/NCBI

[219] Google Scholar

[ref61] 61. Xia WX, Xiao ZA, Cao P, Zhang Y, Du KB, Wang N. Construction of a high-density genetic map and its application for leaf shape QTL mapping in poplar. Planta. 2018;248(5):1173–85. WOS:000447030900009. pmid:30088086
View Article
PubMed/NCBI
Google Scholar

[221] View Article

[222] PubMed/NCBI

[223] Google Scholar

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

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

Figures

Abstract

Introduction

Materials and methods

Genetic experimental design and measurement of leaf traits

SNP genotyping

Statistical methods for GWAS

Investigation of candidate genes

Results

Leaf trait data

SNP genotype data

Significant SNPs associated with leaf traits

Candidate genes affecting leaf shape

Discussion

Conclusion

Supporting information

S1 Fig.

S2 Fig.

S1 Table. The RADseq data information for the 2 parents and 163 progeny in the F1 hybrid population of Populus deltoides and Populus simonii.

S2 Table. The raw data of leaf length, different widths and area from a randomized complete block design in Populus.

S3 Table. Relative difference between leaf measurements with the two software ImageJ and LeafShape.

S4 Table. Correlation coefficients among the leaf traits of L, W, W31, W21, W32, A, and the ratio of L to W in the randomized complete block design derived from the F1 progeny of Populus deltoides × Populus simonii.

S5 Table. Analysis of variance for the leaf parameters of L, W, W1/3, W1/2, W2/3, and Area in the randomized complete block experiment derived from the F1 progeny of Populus deltoides × Populus simonii.

S6 Table. Canonical correlation coefficients among the leaf length, widths, area, the length/width ratio, and polar radii in the randomized complete block design derived from the F1 progeny of Populus deltoides × Populus simonii.

S7 Table. Correlation coefficients between the first principal component of different radius datasets and the leaf length, different widths, or area in the randomized complete block design derived from the F1 progeny of Populus deltoides × Populus simonii.

S8 Table. Summary of significant SNPs associated to each trait of the four different leaf widths and area without genomic control.

S9 Table. Summary of significant SNPs associated to the ratio of the leaf length to the maximum width without genomic control.

S10 Table. The annotation of candidate genes for the significant SNPs with the non-redundant protein database at NCBI and GO database.

S11 Table. Consistency between significant SNPs for poplar leaf shape identified in the current and previous studies.

Acknowledgments

References