Next Article in Journal
Screening Evaluation of Antiproliferative, Antimicrobial and Antioxidant Activity of Lichen Extracts and Secondary Metabolites In Vitro
Previous Article in Journal
Analysis of Fast Fluorescence Kinetics of a Single Cyanobacterium Trapped in an Optical Microcavity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Functional Mapping of Genes Modulating Plant Shade Avoidance Using Leaf Traits

1
Center for Computational Biology, College of Biological Sciences and Technology, Beijing Forestry University, Beijing 100083, China
2
National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing 100083, China
3
Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, Beijing 100083, China
4
The Tree and Ornamental Plant Breeding and Biotechnology Laboratory of National Forestry and Grassland Administration, Beijing Forestry University, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Plants 2023, 12(3), 608; https://doi.org/10.3390/plants12030608
Submission received: 4 December 2022 / Revised: 19 January 2023 / Accepted: 23 January 2023 / Published: 30 January 2023
(This article belongs to the Section Plant Genetics, Genomics and Biotechnology)

Abstract

:
Shade avoidance syndrome (SAS) refers to a set of plant responses that increases light capture in dense stands. This process is crucial for plants in natural and agricultural environments as they compete for resources and avoid suboptimal conditions. Although the molecular, biochemical, and physiological mechanisms underlying the SAS response have been extensively studied, the genetic basis of developmental variation in leaves in regard to leaf area, petiole length, and leaf length (i.e., their allometric relationships) remains unresolved. In this study, with the recombinant inbred line (RIL) population, the developmental traits of leaves of Arabidopsis were investigated under two growth density conditions (high- and low-density plantings). The observed changes were then reconstructed digitally, and their allometric relationships were modelled. Taking the genome-wide association analysis, the SNP genotype and the dynamic phenotype of the leaf from both densities were combined to explore the allometry QTLs. Under different densities, leaf change phenotype was analyzed from two core ecological scenarios: (i) the allometric change of the leaf area with leaf length, and (ii) the change of the leaf length with petiole length. QTLs modulating these two scenarios were characterized as ‘leaf shape QTLs’ and ‘leaf position QTLs’. With functional mapping, results showed a total of 30 and 24 significant SNPs for shapeQTLs and positionQTLs, respectively. By annotation, immune pathway genes, photosensory receptor genes, and phytohormone genes were identified to be involved in the SAS response. Interestingly, genes modulating the immune pathway and salt tolerance, i.e., systemic acquired resistance (SAR) regulatory proteins (MININ-1-related) and salt tolerance homologs (STH), were reported to mediate the SAS response. By dissecting and comparing QTL effects from low- and high-density conditions, our results elucidate the genetic control of leaf formation in the context of the SAS response. The mechanism with leaf development × density interaction can further aid the development of density-tolerant crop varieties for agricultural practices.

1. Introduction

Over the course of evolution, plants have acquired an avoidance strategy to adjust their growth when exposed to suboptimal conditions, similar to an escape strategy in animals. Under highly shaded conditions, most plants tailor their photosynthetic strategy to decreased light levels while at the same time expressing shade-avoidance syndrome (SAS). Plants exhibit SAS by positioning leaves into well-lit microsites; accelerating the elongation of hypocotyls, internodes, and petioles; and elongating shoots away from the shade of neighboring plants to reach light. This plays out both in nature and in agricultural fields, for example, in high-density cropping systems. Thus, a detailed understanding of this process could help improve agricultural productivity.
The molecular and physiological underpinnings of the SAS response have been well studied and are covered in previous studies. Chlorophyll in leaf selectively absorb red light (R) and cannot absorb far-red light (FR), so when increasing the plant density, the unabsorbed FR will be reflected by leaves, causing the reduction of the R:FR ratio. To sense the density change, through the detection of FR reflection, the plant inactivates phytochrome (i.e., phyA to phyE). This subsequently modulates a network of transcription factors, hormones, and cell wall-modifying proteins to induce shade avoidance responses [1,2,3]. Upon exposure to R, the cytosolic phyB directly interacts with phytochrome-interacting factors (PIFs) [1]. In particular, PIF4, PIF5, and PIF7 are associated with SAS to induce the hypocotyl elongation [4]. In addition to PIFs, LONG HYPOCOTYL IN FAR-RED1 (HFR1), PHYTOCHROME RAPIDLY REGULATED1 (PAR1), and PAR2 are likely to be the negative regulators of hypocotyl elongation [5,6]. Additionally, in response to low R:FR, ATHB4 are transcriptionally induced to modulate hormone responses [7]. Indole acetic acid (IAA) level in the elongating hypocotyl can specifically be elevated in low R:FR, which appears to be conditional upon PIN3-mediated auxin transport [8]. Additionally, the low ratio of R:FR can induce the expression of GA20ox, which has a function of enhancing the level of GA, triggering the degradation of growth-inhibiting nuclear proteins of DELLA in cascade, and finally bring out an elongated plant phenotype [9]. All this knowledge is essential for promoting the further development of modern agriculture, e.g., increasing crop yields without agricultural expansion, whilst using the density-tolerant genotype under a higher planting density. Although key elements mediating SAS have been identified, we still argue that further work is needed to assess the genetic architecture of SAS, e.g., specific effects of genes in mediating the SAS-inferred phenotypic change, and the mechanistic basis of natural genetic variation in the SAS response.
Recent studies have reported quantitative trait loci (QTLs) that affect the SAS response. In soybean, a major QTL, named qSAR1, has been mapped to plant height [10]. In Arabidopsis, EODINDEX1 has been identified as a QTL that strongly impacts hypocotyl inhibition under white light [11]. ELF3 is a QTL involved in the shade avoidance response index based on leaf morphology traits [12]. Among phenotype changes of the SAS response, the leaf is one of the most important responsive organs, as leaves undergo photosynthesis and detect cue changes in light spectral composition. However, these studies were mainly based on static assessments of leaf traits, whereas leaf formation is a dynamic process that changes during the growth process. In addition, the allometric scaling relationships remain unknown, that is, changes in leaf length and/or shape are very flexible. During ontology, the leaf length varies with petiole length, thus causing the position change of the leaf, which facilitates the leaf to gain more light. These could potentially be modeled based on the allometric relationships between petiole length (P) and leaf length (L), and between leaf length (L) and leaf area (A).
To incorporate the biological scaling allometry law into the mathematical framework of QTL mapping, while retaining the dynamic phenotype of leaf development, we emphasize the superiority of functional mapping [13,14] in locating SAS genes from the aspect of allometry QTLs. Compared with other QTL mapping methods that use the static phenotype from only one timepoint, functional mapping using traits investigated at a series of timepoints can increase the model power for QTL detection, because the developmental process can be modelled well using mathematical equations with functional mapping. For allometry QTLs, through previous studies, genes mediating height-diameter allometry for the tree stem had been computationally identified, with the allometric growth equation equipped in the functional mapping framework [15]. Then, this information can be used with functional mapping to identify the underlying QTL controlling changes in these trait relationships and, thus, deepen our understanding of developmental regulation of SAS in leaves [13,14].
According to the advantage of functional mapping, we design two density conditions and incorporate the dynamic petiole and leaf traits from two density sets into a functional mapping framework. From a methodological viewpoint, this strongly support the detection of SAS QTL, but through a QTL × density interaction approach. Taking an Arabidopsis recombinant inbred line (RIL) as the mapping population, we propose the use of bivariate functional mapping to study the genetic architecture of SAS in this study. In particular, how genes modulate the dynamic covariation of petiole with leaf length, how leaf length covaries with leaf area, and how genes mediate the interaction of leaf development and density environment were all comprehensively analyzed in this study.

2. Result

2.1. Shape Reconstruction

We used the shape reconstruction method described by Fu et al [16]. On principle, with the known size of a black/white square, the coordinate system can be build using the corner point of the checkerboard grid. Figure 1 shows the results for two types of Arabidopsis leaf, one with a very short petiole (Figure 1A) and the other with a tenuous petiole (Figure 1B). Using this shape reconstruction method, a series of coordinates along the leaf boundary were determined for both the leaves (Figure 1C,D), which match the photographic images. Figure 1C,D fully possessed the feature of the petiole base from Figure 1A,B. Additionally, the shape with the leaf tip and other parts are reserved, suggesting a strong power with a such method. After the shape reconstruction, the traits of leaf length, leaf area, and petiole length can further be quantified.

2.2. Longitudinal Phenotypes of Leaf Traits under Different Density Conditions

Next, using 1680 seedlings, we depicted the phenotypic changes in A, L, and P at the progeny level (Figure 2). Taking all eight timepoints for low-density and seven timepoints for high-density conditions, all three traits showed S-shaped curves during leaf development. In addition, curve-fitting at the population level and between the two conditions, indicated a high variance for all three traits. Overall, as compared with the early stage, the high-density condition led to higher growth rates and larger traits in later stages, supporting the general understanding that the elevation of plant density will cause the alteration of light in quality and quantity, to affect the plant growth or elongation synthetically.
Low-density conditions led to more variance among phenotypes at nearly all timepoints, suggesting an activated SAS response. As shown in Figure 2, by comparing the mean curve from the two-density condition, their crossover points were located at 4.05, 3.67, and 3.23 weeks for the leaf area (Figure 2A), leaf length (Figure 2B), and petiole length (Figure 2C), respectively. The earliest crossover for the petiole trait might imply the sensitivity of petiole growth regulation when facing the changing density condition, suggesting that petiole growth regulation is sensitive to density conditions, and this trait may be the first to trigger the process of SAS.

2.3. Mathematical Modeling of Leaf Allometry

To verify the fitness between the allometric growth equation and the measured phenotypes, the R-squared value for plants grown under each condition was evaluated using the least squares method. Under high-density conditions, more than 71% of progenies had R-squared values greater than 0.73 for the P-A relationship. For the L-A and P-L relationships, more than 93% and 92% of progenies had R-squared values greater than 0.96 and 0.92, respectively. Under low-density conditions, 75%, 92%, and 90% of progenies had goodness-of-fit values of 0.75, 0.94, and 0.95, respectively. These results indicate that the allometry between P and A was not strong. Therefore, we subsequently focused on analyzing allometry for the latter two pairs of traits, i.e., L-A and P-L.
Figure 3 shows the allometric patterns for L-A and P-L during development. L-A has a concave-shaped curve (Figure 3A,B), whereas that for P-L is convex (Figure 3D,E). Although the allometric pattern for each progeny varied, those associated with each density were remarkably similar. The data in Figure 3C indicate that, for a given lamina area, high-density conditions led to slightly longer leaves, indicating that leaf elongation occurs to ensure light acquisition in such environments. Nevertheless, petiole lengths were nearly equal between the two densities (Figure 3F). This suggests that, under high density, petiole elongation does not necessarily account for the increase in leaf length. This coincides with the nonsignificant goodness-of-fit of the allometry data for P-A. Considering the substantial role of the petiole in changing the leaf position and the considerable variation in P-L under both densities, the allometric relationships between P and L can provide important insights into the genetic dissection of the SAS mechanism. By estimating the power coefficients from Equation (1), how one trait scales with the other can be illustrated, and the scaling pattern between the two densities can be compared.

2.4. QTLs for SAS

By integrating phenotypes from the high- and low-density groups into the QTL mapping framework, a bivariate functional mapping model incorporating Equation (1) was implemented. We tested how QTL regulates the allometric scaling of L-A and P-L during leaf ontogeny. We argue that L-A and P-L allometry has different developmental significance for leaves. For L-A, at any given length, leaves with a smaller A are narrower, and this might lead to better photosynthetic ability than larger leaves, increasing competitiveness under crowded and shaded conditions.
To help the interpretation of QTLs, we performed a gene enrichment analysis using the combined genes of two set of analysis. The enrichment result based on Gene Ontology (GO) showed six classes of functions with p-values less than 0.05: red, far-red light phototransduction (GO:0009585, 2.56 × 10−7), gibberellin 20-oxidase activity (GO:0045544, 3.28 × 10−5), shade avoidance (GO:0009641, 0.0014), response to light stimulus (GO:0009416, 0.038), leaf shaping (GO:0010358, 0.021), regulation of defense response (GO:0031347, 0.015). In comparison with the background gene number, 8/32, 4/8, 4/26, 2/12, 2/28, 5/83 genes were involved in these functions, respectively. In particular, for leaf development, the detailed annotation of ‘leaf pavement cell development’, ‘leaf shaping’, and ‘leaf vascular tissue pattern formation’ were found. For class 6, the ‘defense response to bacterium’, ‘regulation of systemic acquired resistance’, ‘defense response to fungus’, ’response to UV-B’, and ‘defense response to bacterium, incompatible interaction’ were found. All these have pointed out the importance of leaf development, immune system, far-red light regulation to the SAS, and have also promoted the screening of key genes for the further dissection of SAS.

2.5. ShapeQTLs: How Leaf Area Scales with Leaf Length

Supporting the general understanding that the elevation of plant density will cause the alteration of light in quality and quantity, this will affect the plant growth or elongation synthetically. Of these, 280, 82, 86, 90, and 7 were distributed between chromosomes 1 to 5 (Figure 4A). All were annotated using the TAIR website (www.arabidopsis.org, accessed on 1 October 2022) according to their genomic position, resulting in 30 annotated functional proteins. In terms of protein function, there were fourteen photosensory receptor genes, four phytohormone genes, six immune pathway genes, and one leaf shape-related gene (Table 1), which implies that SAS is regulated through a multi-step mechanism. The photosensory receptor pathway had the greatest number of genes, suggesting its paramount role in regulating the SAS response. For these key genes, their phenotype variation explained (PVE) was calculated, ranging from 0.99% to 1.91%.
Next, we focused on an SNP located at AT4G01895.1, which encodes systemic acquired resistance (SAR) regulatory proteins (NIMIN-1-related) belonging to the immune regulation pathway. As shown in Figure 5A, each genotype at this shapeQTL showed a difference in the pattern of how leaf area and leaf length change with time. With the two-dimensional expansion of the leaf, A showed a period of linear growth from 2 to 4 weeks and increased substantially, compared to the increase in L. Indeed, the curves for A showed greater differences between density conditions than between genotypes. Over time, the allometric difference between the two densities narrowed (Figure 5B). Regarding genotypes, CC and GG showed significant differences between the two densities. In addition, all 10 shapeQTLs showed a weakening impact on the L-A relationship during the middle-late stages of growth after a slight increase in the early stage (Figure 6A). Once a leaf reaches 4 cm in length, the additive effects decrease. In comparing effects between the two densities, the inversive patterns of 10 additive effect curves were disclosed (Figure 6B).

2.6. PositionQTLs: How Leaf Length Scales with Petiole Length

There were 418 significant SNPs identified for P-L allometry, involving 280, 30, 17, 59, and 32 SNPs on chromosomes 1 to 5 (Figure 4B). Functional annotation identified 24 protein genes classified into the categories listed in Table 2. The pathways and gene categories were the same as those for shapeQTLs, clearly indicating their shared function in regulating SAS. Therefore, positionQTLs not only regulate leaf position by changing L or P, but also affect traits that are regulated by shapeQTLs. These key genes owned PVEs ranging from 1.01% to 1.88%.
Next, we focused on a QTL located at AT3G22170.1, which encodes far-red elongated hypocotyls 3 protein. Under both densities, P changed with time; although its growth rate was slower than that of L, both P and L changed more at the density level than at the genotype level (Figure 5C). High-density conditions led to a larger P than low-density conditions, which suggests that this QTL plays a role in developing a longer petiole to move leaves towards well-lit microsites. For L, the curves diverged between genotypes after 1.5 weeks under low density. Using the estimated parameters from Equation (1), we analyzed the allometry pattern at the genotype level. The positionQTLs fit the allometry scaling relationship for P-L well (Figure 5D). The allometric patterns among QTL genotypes under high density were indistinguishable. This suggests a G × E interaction property with this QTL. For genotypes AA and TT, the differences in allometry scaling between the two densities were significant. Calculations of the additive effects under low density indicated that QTLs 4, 5, and 7 showed a downward parabola-shaped additive effect. QTLs 3, 6, 9, and 10 showed a sharp decrease with an increase in P (Figure 6C). Under high density, QTLs 1, 2, 3, 6, 9, and 10 showed a decreasing role in controlling allometry in relation to L. QTLs 5, 7, 8, and 9 showed an increasing role with P (Figure 6D).
Figure 6. Additive effect with significant SNPs. (A,B) show the additive effect of shapeQTLs with 10 significant genes under low- and high-density conditions. Numbers 1–10 present genes AT1G17020.1, AT1G19850.1, AT2G03340.1, AT2G27110.1, AT3G59060.1, AT3G63300.1, AT4G16250.1, AT4G25420.1, AT5G28530.1 and AT5G40910.1, respectively. (C,D) show the additive effect of positionQTLs with 10 significant genes under low- and high-density conditions. Numbers 1-10 present AT1G29690.1, AT1G10240.1, AT1G51090.1, AT2G43010.1, AT3G22170.1, AT3G58850.1, AT4G18130.1, AT4G27430.1, AT5G38860.1 and AT5G51810.1, respectively.
Figure 6. Additive effect with significant SNPs. (A,B) show the additive effect of shapeQTLs with 10 significant genes under low- and high-density conditions. Numbers 1–10 present genes AT1G17020.1, AT1G19850.1, AT2G03340.1, AT2G27110.1, AT3G59060.1, AT3G63300.1, AT4G16250.1, AT4G25420.1, AT5G28530.1 and AT5G40910.1, respectively. (C,D) show the additive effect of positionQTLs with 10 significant genes under low- and high-density conditions. Numbers 1-10 present AT1G29690.1, AT1G10240.1, AT1G51090.1, AT2G43010.1, AT3G22170.1, AT3G58850.1, AT4G18130.1, AT4G27430.1, AT5G38860.1 and AT5G51810.1, respectively.
Plants 12 00608 g006

2.7. Pleiotropism among QTLs

Interestingly, 89 significant SNPs were identified for shapeQTLs that were also significant SNPs for positionQTLs, which suggests that these SNPs simultaneously govern allometric scaling of both L-A and P-L. There were 55, 6, 5, 21, and 2 shared SNPs located on chromosomes 1 to 5. As the two allometric relationships have different biological meanings, these QTLs appear to be pleiotropic. AT1G10240.1, AT1G69010.1, and AT5G38860.1 were directly associated with the SAS response. These SNPs encode FAR1-related sequence 11, BES1-interacting Myc-like protein 2, and BES1-interacting Myc-like protein, respectively, and were all involved in photosensory pathways. The QTL at AT1G10240.1 had a significant role in regulating L-A (Figure 7A,B) and P-L allometry under both densities (Figure 7C,D). Such a QTL may behave more sensitively to alter the phenotype of the leaf when density changes, so that plants can produce a flexible mechanism to better adapt to shade.

3. Discussion

SAS triggers phenotypic changes in plants, as they attempt to avoid shade and seek enhanced light interception. For plants, shade varies both in the vertical and horizontal directions and also temporally, and shade avoidance can involve any one or more of these axes. All these hold important implications for agriculture practice. For example, in the context of decreasingly available arable land and the increasing population, rational close cultivation is one of the effective strategies to increase crop yield per unit area. To date, although studies have identified dominant components of SAS regulation, the dissection on genetic control of SAS plasticity had largely not factored the co-varying relationships among distinct phenotypes of leaf.
In this study, we focused on the allometric relationships among leaf traits. Specific traits were investigated longitudinally and nondestructively by modeling the changes observed in laboratory experiments. We also modeled the allometric relationships among traits and determined how QTLs may control allometry. The association analysis for each SNP have been used to determine the genetic architecture of complex traits [17,18]. Functional mapping, a QTL identification framework that utilizes mathematical curves to model biological growth, was conceptualized two decades ago [13]. Although only one mapping population used in this study, we hope to construct a second mapping population to validate the QTL results to improve the reliability. Compared to GWAS using static traits, the use of dynamic phenotypes can provide deeper biological insights by correlating genotypes with phenotypes. Thus, since this method was introduced, it has been extensively used, for example, to dissect gravitropism and phototropism based on dynamic phenotypes of seedlings [19], to study leaf allometry between leaf area and leaf mass in the common bean [20], and to model phenotypic plasticity using growth trajectories of rice [21].
In this study, L-A and P-L relationships were used as a predictor of the capacity of Arabidopsis to respond to shade. Through functional mapping, the genetic architecture underlying allometry was dissected. Focusing on traits of leaf development and incorporating phenotypes from plants grown under two different densities, we were able to map shapeQTLs and positionQTLs. We also identified pleiotropic genes. Among genes modulating the change in A as a function of L, we identified 545 functional genes in the shapeQTLs group, 30 of which encode a molecular function: these included photosensory receptor genes, phytohormone genes, and genes regulating leaf development and leaf shape. Most photosensory receptors were phytochromes, which respond to low R:FR ratios caused by dense vegetation [22]. AT3G63300.1, encoding FORKED 1 (FKD1), plays a role in patterning leaf vascular and can coordinate leaf size with vein density [23]. Interestingly, AT1G75540.1, which encodes salt tolerance homolog2 (STH2), was also associated with shade avoidance. It interacts with COP1 to control de-etiolation and positively regulates photosynthesis [24]. Interestingly, high-density growth conditions lead to more severe plant diseases [25] which may be due to the change in light environment (e.g., a decrease in UV-B radiation typically decreases plant resistance to biotic stressors).
Genes modulating changes in L as a function of P (i.e., positionQTLs) fell into the same pathway categories as those for shapeQTLs but the gene sets were different. This implies that changes in leaf position and leaf shape are linked to the SAS response. The overlap of gene categories between L-A and P-L may reflect their combined role in driving changes in leaf position and shape. On the other hand, this overlapping further illustrated the need for QTL mapping of SAS using distinct allometric trait pairs; indeed, doing so identified another phenotype regulating SAS under high-density conditions, namely, the A trait. This indicates that there are QTL–density interactions. Remarkably, a previous study reported the role of AT3G22170 in mediating the elongation of hypocotyls in response to FR light in concert with FAR1 [26,27]. However, in our study, this QTL was related with P-L, revealing another role during later development.
The allometry underlying complex traits is pervasive. For example, the metabolic theory of ecology explaining metabolic rate scales with the three-fourths power of body mass [28]. Another example is the quarter-power scaling relationship between separate organs and overall body size [29]. As a practical case, in forms of height-diameter allometry [15] identified genes modulating the tree stem growth. Using Reeve and Huxley’s allometric equation [30], the increase in the dry weight of soybean seed was examined using a functional mapping method [31]. Therefore, the allometry principle deserves a widespread application to the dissection of more complex traits. In this study, we used functional mapping to map allometry-related QTLs, with an emphasis on SAS. We utilized an RIL population, which can give sufficient replicate samples, as each progeny within the RIL population can propagate many replicates through self-reproduction. This can reduce spurious associations when mapping QTLs. In addition, we reduce other factors by growing plants in chambers with controllable conditions. All these steps improved the precision of QTL detection.
For rosette species such as Arabidopsis, vegetative internode elongation is highly suppressed and we, therefore, used leaf traits to analyze the genetic architecture of SAS. However, for other non-rosette plants, accelerated elongation of internodes and the resulting increase in stem length have been reported as the most striking SAS component [32]. This might also be accompanied by enhanced apical dominance and stimulated vertical growth of the main stem. In such species, it may be better to adjust the model to consider the growth relationship between internode length and stem height. The further modeling of such potential allometry could reveal molecular mechanisms underlying dwarf or semidwarf phenotypes, which could aid the development of modern crop species that can grow at high densities.

4. Materials and Methods

4.1. Plant Materials

Arabidopsis ecotype Ler (Landsbergerecta) and Sha (Shakdara) were used as parental material. F1 seeds were obtained by manual crossing and subsequent self-fertilization for 10 generations to develop a RIL population. In all, 84 progenies were used, with each progeny, 20 replicates were taken for both the high- and low-density conditions. The genome DNA for the whole RIL population was extracted with a TIANGEN DP305 DNA extraction kit, using 100 mg of young leaf. DNA quality was inspected by agarose electrophoresis and imaged with an ultraviolet imager. DNA samples without degradation were sequenced subsequently. The RIL population was genotyped at Total Genomics Solution (TGS) Institute, China. Using the Illumina Hiseq 2000 platform, the 125 bp paired-end reads were produced at a depth of 9.91×. After filtering the low-quality sequence, 1.10 GB data per sample were generated on average. The percentage of Q30 and Q20 bases were more than 90% and 95%, respectively, with a normal GC distribution. Taking the reference genome from the TAIR website (https://www.arabidopsis.org/download/index-auto.jsp?dir=%2Fdownload_files%2FSequences%2FAssemblies, accessed on 1 October 2022), SNP detection was performed with GATK, a widely used mutation detection software. The SNP screening criteria were customized as follows: the sequencing depth ≥ 4, the SNP base quality ≥ 20, the variation detection quality ≥ 50. SNPs with the allele frequency ≥ 5%, or with the missing rate ≥ 50% were filtered out. With these criteria, a total of 417,495 high-quality SNPs were determined for QTL mapping [33].

4.2. Experimental Design

From the 84 progenies, 40 seeds were sowed, from which 20 seedlings were used for high- and low-density planting (with 3 × 3 cm and 5 × 5 cm row spacing for the high and low conditions, respectively). For each density condition, a randomized block design of 20 replicates with each progeny, 5 blocks each with 4 replicates, was used. Plants were grown in an artificial climate chamber at Beijing Forestry University, China. Light and temperature conditions in the chamber were 16 h light and 8 h dark per day, with a light intensity of 150 μmol·s–1·m–2 and a constant room temperature of 22 °C with 80% relative humidity. A mixture of turfy soil, perlite, and vermiculite at a 1:1:1 ratio was used for cultivation.

4.3. Phenotype Collection

One week after transplanting, the third true leaf of each plant was chosen for phenotypic investigation. Selected leaves were photographed weekly using a digital camera (NIKON COOLPIX A10). To assess phenotypes, an image calibration card (a printed checkerboard, each square having a size of 10 × 10 mm) was used as illustrated in Figure 1. At the time of taking the photograph, the leaf was flattened against the blank space of the card, making the checkerboard and leaf integrity show on photo (Figure 1A,B). The information was imported into the CameraCalibrator module of MATLAB for further analysis. Then, following a previously described method [16], all leaf shapes were reconstructed. Furthermore, the three traits under study (P, L, and A) were assessed via MATLAB.

4.4. Statistical Models

Mean phenotypes were calculated for each progeny in the RIL population, and each of the three traits were calculated. Equation (1) was used to model allometric growth:
μ j = α j x i β j
where xi is the longitudinal phenotype of one trait for progeny i, μj is the mean vector of another trait for SNP genotype j, and αj and βj are parameters for the allometric law.
To incorporate the dynamic phenotype from the two densities into one model, taking the scaling relationship of L-A as an example, let xi = (xih; xil) = (xih(t1),(xih(t2),…, (xih(tTh); xil(t1),(xil(t2),…, (xil(tTl)) denote the joint phenotype of L for progeny i, where h represents the high-density environment and l represents the low-density environment. Similarly, yi = (yih; yil) = (yih(t1),(yih(t2),…, (yih(tTh); yil(t1),(yil(t2),…, (yil(tTl)) expresses the joint phenotype of A. For the two-density condition, we allowed the starting measurement time to vary depending on the growth after seedling transplantation.
Next, assuming two genotypes for each SNP, n progenies can be grouped into two sets. The likelihood for the dynamic trait from the two-density condition is formulated as follows:
L y = i = 1 n f j | i y i h ;   y i l   | μ j , Σ i  
According to this equation, at a specific SNP, the given progeny i carrying QTL genotype j, the function f j | i y i h ;   y i l   | μ j , Σ i   is a multivariate normal distribution, where the expected mean vector μj is expressed as
μ j = ( μ j h ;   μ j l ) = ( μ j h t 1 , μ j h t 2 , , μ j h t T h ;   μ j l ( t 1 ) ,   μ j l ( t 2 ) ,   μ j l ( t T l )
For (co-)variance matrix i, because xi in Equation (1) denotes one trait within the allometric trait pair, xi has progeny-specific values. To synthesize the two density conditions, i can be expressed as
Σ i = Σ i h h Σ i h l Σ i l h Σ i l l
where i(hh) and i(ll) represent a (Th × Th) and (Tl × Tl) residual covariance matrix of growth over the y trait, respectively, and i(hl) = (i(lh))’ is a (Th × Tl) covariance matrix between the two density conditions. Based on a previous derivation of multi-trait functional mapping [34], combined with a one-order structured ante-dependence model (SAD (1)) [35] we further rewrite (4) as
Σ i = Λ Σ ε Λ
where
Λ = I 0 0 0 V T 2 T 1 I 0 0 V T 3 T 1 V T 3 T 2 I 0 V T T T 1 V T T T 2 V T T T T 1 I
and V = φ 1 0 0 φ 2 ;   φ 1 and φ 2   refer to the parameters of the two-dimensional SAD (1), Σ ε = diag( σ 2 ).
A maximum likelihood estimation was used to estimate curve parameters ( α j , β j ) = ( α j h , α j l , β j h , β j l ) for the mean vector, and (co-)variance parameters ( φ 1 , φ 2   , σ 2 ). Under likelihood (2), whether a QTL for SAS regulation exists can be tested by formulating the null and alternative hypotheses as follows:
H 0 :   α j = α ,   β j = β H 1 :   At   least   one   of   the   above   expression   does   not   hold
For this hypothesis, we used the likelihood ratio test to determine QTLs, the threshold for which can be determined through 1000 replicates of a permutation test, using a significance level of 0.01.
L R = 2 l o g L 0 l o g L 1
We can also test how QTLs may change allometric scaling relationships under different densities to impact the SAS response. This can give deeper insights into the genetic control of leaf traits to optimize photosynthetic properties. The equation can be expressed as follows:
H 0 :   α j h = α j l ,   β j h = β j l H 1 :   At   least   one   of   the   above   expressions   does   not   hold
We performed the same calculations for P and L to determine their allometric relationship. Then, considering the three correlated leaf traits (P, A, and L), the two allometric scaling relationships, and the two (high/low) density conditions, the QTLs responsible for the SAS response were comprehensively assessed.

4.5. Functional Annotation

Using the models above, all SNPs with a high LOD value larger than the threshold were taken as QTLs. All SNPs were aligned against the Arabidopsis genome; thus, the annotation result with all SNPs can be determined, including their genomic position, and their reported function if SNPs were located at functional proteins.

5. Conclusions

Using the developmental trait of leaf, this study performed the QTL mapping of plant SAS with the Arabidopsis RIL population. With the aid of bivariate functional mapping, leaf traits investigated under a high- and low-density planting environment were mathematically modeled from the viewpoint of allometry. For shapeQTL and positionQTL results using the L-A and P-L allometry, a total of 30 and 24 functional SNPs were determined to mediate the plant SAS, respectively. Through functional annotation, immune pathway genes, photosensory receptor genes, and phytohormone genes were identified to be involved in the SAS response, among which, systemic acquired resistance (SAR) regulatory proteins (NIMIN-1-related) and salt tolerance homolog2 (STH2) were identified. By dissecting and comparing QTL effects, the results elucidate the genetic control of leaf formation in the context of the SAS response.

Author Contributions

Conceptualization, H.Z., Y.C. and M.Y.; methodology, H.Z. and Y.C.; validation, H.Z., Y.C. and M.Y.; formal analysis, H.Z. and Z.W.; software, H.Z. and Z.W.; Supervision, M.Y. and R.W.; writing- original draft, H.Z.; writing - review& editing, H.Z., Z.W., M.Y. and R.W. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (32071796).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, M.; Chory, J. Phytochrome signaling mechanisms and the control of plant development. Trends Cell Biol. 2011, 21, 664–671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Casal, J.J. Photoreceptor signaling networks in plant responses to shade. Annu. Rev. Plant Biol. 2013, 64, 403–427. [Google Scholar] [CrossRef] [PubMed]
  3. Fraser, D.P.; Hayes, S.; Franklin, K.A. Photoreceptor crosstalk in shade avoidance. Curr. Opin. Plant Biol. 2016, 33, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zhang, Y.; Mayba, O.; Pfeiffer, A.; Shi, H.; Tepperman, J.M.; Speed, T.P.; Quail, P.H. A quartet of PIF bHLH factors provides a transcriptionally centered signaling hub that regulates seedling morphogenesis through differential expression-patterning of shared target genes in Arabidopsis. PLoS Genet. 2013, 9, e1003244. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Roig-Villanova, I.; Bou-Torrent, J.; Galstyan, A.; Carretero-Paulet, L.; Portoles, S.; Rodriguez-Concepcion, M.; Martinez-Garcia, J.F. Interaction of shade avoidance and auxin responses: A role for two novel atypical bHLH proteins. Embo J. 2007, 26, 4756–4767. [Google Scholar] [CrossRef] [Green Version]
  6. Galstyan, A.; Cifuentes-Esquivel, N.; Bou-Torrent, J.; Martinez-Garcia, J.F. The shade avoidance syndrome in Arabidopsis: A fundamental role for atypical basic helix-loop-helix proteins as transcriptional cofactors. Plant J. 2011, 66, 258–267. [Google Scholar] [CrossRef]
  7. Sorin, C.; Salla-Martret, M.; Bou-Torrent, J.; Roig-Villanova, I.; Martinez-Garcia, J.F. ATHB4, a regulator of shade avoidance, modulates hormone response in Arabidopsis seedlings. Plant J. 2009, 59, 266–277. [Google Scholar] [CrossRef]
  8. Keuskamp, D.H.; Pollmann, S.; Voesenek, L.A.; Peeters, A.J.; Pierik, R. Auxin transport through PIN-FORMED 3 (PIN3) controls shade avoidance and fitness during competition. Proc. Natl. Acad. Sci. USA 2010, 107, 22740–22744. [Google Scholar] [CrossRef] [Green Version]
  9. Schwechheimer, C. Understanding gibberellic acid signaling--are we there yet? Curr. Opin. Plant Biol. 2008, 11, 9–15. [Google Scholar] [CrossRef]
  10. Zhao, J.; Shi, X.; Chen, L.; Chen, Q.; Tian, X.; Ai, L.; Zhao, H.; Yang, C.; Yan, L.; Zhang, M. Genetic and transcriptome analyses reveal the candidate genes and pathways involved in the inactive shade-avoidance response enabling high-density planting of soybean. Front. Plant Sci. 2022, 13, 973643. [Google Scholar] [CrossRef]
  11. Jiang, Y.; Yang, C.; Huang, S.; Xie, F.; Xu, Y.; Liu, C.; Li, L. The ELF3-PIF7 interaction mediates the circadian gating of the shade response in Arabidopsis. Iscience 2019, 22, 288–298. [Google Scholar] [CrossRef]
  12. Jimenez-Gomez, J.M.; Wallace, A.D.; Maloof, J.N. Network analysis identifies ELF3 as a QTL for the shade avoidance response in Arabidopsis. PLoS Genet. 2010, 6, e1001100. [Google Scholar] [CrossRef] [Green Version]
  13. Ma, C.X.; Casella, G.; Wu, R. Functional mapping of quantitative trait loci underlying the character process: A theoretical framework. Genetics 2002, 161, 1751–1762. [Google Scholar] [CrossRef]
  14. Wu, R.; Lin, M. Functional mapping—how to map and study the genetic architecture of dynamic complex traits. Nat. Rev. Genet. 2006, 7, 229–237. [Google Scholar] [CrossRef]
  15. Jiang, L.; Ye, M.; Zhu, S.; Zhai, Y.; Xu, M.; Huang, M.; Wu, R. Computational identification of genes modulating stem height-diameter allometry. Plant Biotechnol. J. 2016, 14, 2254–2264. [Google Scholar] [CrossRef] [Green Version]
  16. Fu, Y.; Li, F.; Mu, S.; Jiang, L.; Ye, M.; Wu, R. Heterophylly quantitative trait loci respond to salt stress in the desert tree Populus euphratica. Front. Plant Sci. 2021, 12, 692494. [Google Scholar] [CrossRef]
  17. Chang, M.; He, L.; Cai, L. An overview of genome-wide association studies. Methods Mol. Biol. 2018, 1754, 97–108. [Google Scholar] [CrossRef]
  18. Gondro, C.; Lee, S.H.; Lee, H.K.; Porto-Neto, L.R. Quality control for genome-wide association studies. Methods Mol. Biol. 2013, 1019, 129–147. [Google Scholar] [CrossRef]
  19. Jiang, P.; Ma, N.; Sun, F.; Ye, M.; Wu, R. Functional mapping of gravitropism and phototropism for a desert tree, Populus euphratica. Front. Biosci. 2021, 26, 988–1000. [Google Scholar] [CrossRef]
  20. Zhang, M.; Zhang, S.; Ye, M.; Jiang, L.; Vallejos, C.E.; Wu, R. The genetic control of leaf allometry in the common bean, Phaseolus vulgaris. Bmc Genet. 2020, 21, 29. [Google Scholar] [CrossRef]
  21. Wang, Z.; Pang, X.; Wu, W.; Wang, J.; Wang, Z.; Wu, R. Modeling phenotypic plasticity in growth trajectories: A statistical framework. Evolution 2014, 68, 81–91. [Google Scholar] [CrossRef] [PubMed]
  22. de Wit, M.; Keuskamp, D.H.; Bongers, F.J.; Hornitschek, P.; Gommers, C.; Reinen, E.; Martinez-Ceron, C.; Fankhauser, C.; Pierik, R. Integration of phytochrome and cryptochrome signals determines plant growth during competition for light. Curr. Biol. 2016, 26, 3320–3326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Prabhakaran, M.N.; Clarke, K.J.; Yu, H.; Wilton, E.E.; Van Dyk, J.; Hou, H.; Schultz, E.A. Members of the Arabidopsis FORKED1-LIKE gene family act to localize PIN1 in developing veins. J. Exp. Bot. 2018, 69, 4773–4790. [Google Scholar] [CrossRef] [PubMed]
  24. Datta, S.; Hettiarachchi, C.; Johansson, H.; Holm, M. SALT TOLERANCE HOMOLOG2, a B-box protein in Arabidopsis that activates transcription and positively regulates light-mediated development. Plant Cell 2007, 19, 3242–3255. [Google Scholar] [CrossRef] [Green Version]
  25. Ballare, C.L.; Mazza, C.A.; Austin, A.T.; Pierik, R. Canopy light and plant health. Plant Physiol. 2012, 160, 145–155. [Google Scholar] [CrossRef] [Green Version]
  26. Liu, L.; Li, B.; Liu, X. FAR-RED ELONGATED HYPOCOTYL3 promotes floral meristem determinacy in Arabidopsis. Plant Signal. Behav. 2016, 11, e1238545. [Google Scholar] [CrossRef] [Green Version]
  27. Lin, R.; Ding, L.; Casola, C.; Ripoll, D.R.; Feschotte, C.; Wang, H. Transposase-derived transcription factors regulate light signaling in Arabidopsis. Science 2007, 318, 1302–1305. [Google Scholar] [CrossRef] [Green Version]
  28. Gillooly, J.F.; Brown, J.H.; West, G.B.; Savage, V.M.; Charnov, E.L. Effects of size and temperature on metabolic rate. Science 2001, 293, 2248–2251. [Google Scholar] [CrossRef] [Green Version]
  29. Brown, J.H.; West, G.B.; Enquist, B.J. Scaling in Biology: Patterns and Processes, Causes and Consequences; Oxford University Press, Inc.: New York, NY, USA, 2000; pp. 1–24. [Google Scholar]
  30. Reeve, E.; Huxley Js, L.G.C.W.; Pb, M. Some problems in the Study of Allometric Growth, Essays on “Growth and Form” Presented to D’Arcy Wentworth Thompson; Oxford Clarendon Press: Oxford, UK, 1945; pp. 121–156. [Google Scholar]
  31. Huang, Z.; Tong, C.; Bo, W.; Pang, X.; Wang, Z.; Xu, J.; Gai, J.; Wu, R. An allometric model for mapping seed development in plants. Brief. Bioinform. 2014, 15, 562–570. [Google Scholar] [CrossRef] [Green Version]
  32. Pierik, R.; Testerink, C. The art of being flexible: How to escape from shade, salt, and drought. Plant Physiol. 2014, 166, 5–22. [Google Scholar] [CrossRef]
  33. Jiang, L.; Shi, C.; Ye, M.; Xi, F.; Cao, Y.; Wang, L.; Zhang, M.; Sang, M.; Wu, R. A computational-experimental framework for mapping plant coexistence. Methods Ecol. Evol. 2018, 9, 1335–1352. [Google Scholar] [CrossRef]
  34. Mu, S.; Zhu, X.; Ye, M.; Wu, R. Application of multi-trait functional mapping of QTL in Populus euphratica Oliv.seedlings. Plant Sci. J. 2022, 40, 365–377. [Google Scholar]
  35. Zhao, W.; Hou, W.; Littell, R.C.; Wu, R. Structured antedependence models for functional mapping of multiple longitudinal traits. Stat. Appl. Genet. Mol. Biol. 2005, 4, e33. [Google Scholar] [CrossRef]
Figure 1. Effect with the shape reconstruction procedure. (A,B) are the original images taken by a digital camera, (C,D) are the reconstructed shapes manifested in coordinate systems.
Figure 1. Effect with the shape reconstruction procedure. (A,B) are the original images taken by a digital camera, (C,D) are the reconstructed shapes manifested in coordinate systems.
Plants 12 00608 g001
Figure 2. Dynamic growth pattern of leaf traits under high-density and low-density conditions. (AC), respectively, present the ontogenic growth curve of the leaf area, leaf length, and petiole length. Yellow curves denote trait patterns from a low-density environment while blue curves denote that from a high-density environment. The two thick lines represent the fitted mean growth curve, while thin lines represents the observed growth phenotype with each progeny in the population. The vertical line indicates the timepoint when the two mean curves for the two densities crossed.
Figure 2. Dynamic growth pattern of leaf traits under high-density and low-density conditions. (AC), respectively, present the ontogenic growth curve of the leaf area, leaf length, and petiole length. Yellow curves denote trait patterns from a low-density environment while blue curves denote that from a high-density environment. The two thick lines represent the fitted mean growth curve, while thin lines represents the observed growth phenotype with each progeny in the population. The vertical line indicates the timepoint when the two mean curves for the two densities crossed.
Plants 12 00608 g002
Figure 3. Fitness of allometric growth with different traits of leaf. (A,B) are the allometric pattern between the leaf length and leaf area, respectively, for high- and low-density conditions. (D,E) are the pattern between the petiole length and leaf length for high- and low- density conditions. (C,F) are the comparison of the mean allometric curve fitted for leaf length–leaf area and petiole length–leaf length, respectively.
Figure 3. Fitness of allometric growth with different traits of leaf. (A,B) are the allometric pattern between the leaf length and leaf area, respectively, for high- and low-density conditions. (D,E) are the pattern between the petiole length and leaf length for high- and low- density conditions. (C,F) are the comparison of the mean allometric curve fitted for leaf length–leaf area and petiole length–leaf length, respectively.
Plants 12 00608 g003
Figure 4. Manhattan plot of association results with shade avoidance syndrome. (A,B) are association results for the allometry association results for leaf length–leaf area and petiole length–leaf length, respectively. The dashed horizontal line represents the LOD threshold at the significance level of 0.05.
Figure 4. Manhattan plot of association results with shade avoidance syndrome. (A,B) are association results for the allometry association results for leaf length–leaf area and petiole length–leaf length, respectively. The dashed horizontal line represents the LOD threshold at the significance level of 0.05.
Plants 12 00608 g004
Figure 5. Expression patterns of allometry QTLs. With QTL located at AT4G01895.1, (A) shows the ontogenic growth curve of leaf length and leaf area, and (B) shows their allometry growth curve. With QTL located at AT3G22170.1, (C) shows the ontogenic growth curve of petiole length and leaf length, and (D) shows their allometry growth curve. The solid line represents low density, the dotted line represents high density, and red and black color express different genotypes with the QTL.
Figure 5. Expression patterns of allometry QTLs. With QTL located at AT4G01895.1, (A) shows the ontogenic growth curve of leaf length and leaf area, and (B) shows their allometry growth curve. With QTL located at AT3G22170.1, (C) shows the ontogenic growth curve of petiole length and leaf length, and (D) shows their allometry growth curve. The solid line represents low density, the dotted line represents high density, and red and black color express different genotypes with the QTL.
Plants 12 00608 g005
Figure 7. Expression patterns of pleiotropic QTL. With QTL located at AT1G10240.1, (A) shows the ontogenic growth curve of leaf length and leaf area, and (B) shows their allometry growth curve. (C) shows the ontogenic growth curve of petiole length and leaf length, and (D) shows their allometry growth curve. The solid line represents low density, the dotted line represents high density, and red and black color express different genotypes with the QTL.
Figure 7. Expression patterns of pleiotropic QTL. With QTL located at AT1G10240.1, (A) shows the ontogenic growth curve of leaf length and leaf area, and (B) shows their allometry growth curve. (C) shows the ontogenic growth curve of petiole length and leaf length, and (D) shows their allometry growth curve. The solid line represents low density, the dotted line represents high density, and red and black color express different genotypes with the QTL.
Plants 12 00608 g007
Table 1. Functional annotation of QTLs associated with allometry between leaf length and leaf area.
Table 1. Functional annotation of QTLs associated with allometry between leaf length and leaf area.
TypeChromosomeLocationGene IDGene DescriptionPVE
Immune pathway15821947AT1G17020.1Senescence-related gene 1 (SRG1)1.312%
111605859AT1G32210.1Death resistance protein1.415%
21021553AT2G03340.1WRKY DNA-binding protein (WRKY33)1.502%
27308099AT2G16870.1TIR-NBS-LRR-like resistance protein1.004%
4818434AT4G01895.1Systemic acquired resistance (SAR) regulatory proteins NIMIN-1-related1.628%
516392619AT5G40910.1TIR-NBS-LRR-like resistance protein1.012%
Photosensory receptor genes13359050AT1G10240.1FAR1-related sequence 11 (FRS11)1.911%
119567276AT1G52520.1FAR1-related sequence 6 (FRS6)1.328%
211581554AT2G27110.1FAR1-related sequence 3 (FRS3)1.401%
510525858AT5G28530.1FAR1-related sequence 10 (FRS10)1.408%
13100411AT1G09570.1Phytochrome A1.909%
28144368AT2G18790.1Phytochrome B1.650%
49203572AT4G16250.1Phytochrome D1.341%
28708022AT2G20180.1Phytochrome interacting factor 3-like 5 (PIL5)1.406%
321829916AT3G59060.1Phytochrome interacting factor 3-like 6 (PIL6)1.326%
125946396AT1G69010.1BES1-interacting Myc-like protein 2 (BIM2)1.073%
515558746AT5G38860.1BES1-interacting Myc-like protein 3 (BIM1)0.991%
217837651AT2G42870.1PHY rapidly regulated 1 (PAR1)1.212%
2344412AT2G01800.1COP1-interacting protein-related (CIP)1.303%
514727698AT5G37190.1COP1-interacting protein 4 (CIP4)1.531%
Phytohormones12218953AT1G07240.1UDP-glucosyl transferase 71C51.307%
125051840AT1G67080.1Abscisic acid (aba)-deficient 41.721%
412993626AT4G25420.12-Oxoglutarate (2OG) and Fe(II)-dependent oxygenase superfamily protein1.405%
413221364AT4G26080.1Protein phosphatase 2C family protein (PP2C)1.702%
Leaf development-related11245253AT1G04550.1AUX/IAA transcriptional regulator family protein (AXR3)1.356%
15883293AT1G17220.1Translation initiation factor 2, small GTP-binding protein (FUG1)1.377%
16152210AT1G17870.1Ethylene-dependent gravitropism-deficient and yellow-green-like 3 (EGY3)1.398%
16882081AT1G19850.1Transcriptional factor B3 family protein/auxin-responsive factor AUX/IAA-related (MP)1.395%
Leaf morphology correlation323386503AT3G63300.1FORKED 1 (FKD1)
1.281%
Others128362736AT1G75540.1Salt tolerance homolog2 (BBX21)1.341%
Table 2. Functional annotation of QTLs associated with allometry between the petiole length and leaf length.
Table 2. Functional annotation of QTLs associated with allometry between the petiole length and leaf length.
TypeChromosomeLocationGene IDGene DescriptionPVE
Immune pathway110379838AT1G29690.1MAC/Perforin domain-containing protein (NSL1)1.310%
118933188AT1G51090.1Heavy metal transport/detoxification superfamily protein (NAKR1)1.202%
124067759AT1G64790.1ILITYHIA (ILA)1.344%
Photosensory receptor genes13359050AT1G10240.1FAR1-related sequence 11 (FRS11)1.011%
410832385AT4G19990.1FAR1-related sequence 1 (FRS1)1.142%
417905395AT4G38170.1FAR1-related sequence 9 (FRS9)1.149%
125946396AT1G69010.1BES1-interacting Myc-like protein 2 (BIM2)1.261%
515558746AT5G38860.1BES1-interacting Myc-like protein 3 (BIM3)1.447%
217881666AT2G43010.1Phytochrome interacting factor 4 (PIF4)1.282%
37826998AT3G22170.1Far-red elongated hypocotyls 3 (FHY3)1.889%
321759608AT3G58850.1PHY rapidly regulated 2 (PAR2)1.156%
413723181AT4G27430.1COP1-interacting protein 7 (CIP7)1.600%
410049797AT4G18130.1Phytochrome E1.717%
514015733AT5G35840.1Phytochrome C1.456%
Phytohormones112448321AT1G34170.1Auxin response factor 131.335%
121976267AT1G59750.1Auxin response factor 11.112%
127657184AT1G73590.1Auxin efflux carrier family protein (PIN1)1.287%
128970875AT1G77110.1Auxin efflux carrier family protein (PIN1)1.271%
521055662AT5G51810.1Gibberellin 20 oxidase 2 (GA20OX2)1.122%
Leaf development related1453280AT1G02280.1Translocon at the outer envelope membrane of chloroplasts 33 (TOC33)1.574%
11245256AT1G04550.1AUX/IAA transcriptional regulator family protein (AXR3)1.725%
11610455AT1G05470.1DNAse I-like superfamily protein1.625%
11682406AT1G05630.1Endonuclease/exonuclease/phosphatase family protein (Eepd1)1.351%
Leaf morphology correlation11105105AT1G04180.1YUCCA 91.381%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, H.; Cao, Y.; Wang, Z.; Ye, M.; Wu, R. Functional Mapping of Genes Modulating Plant Shade Avoidance Using Leaf Traits. Plants 2023, 12, 608. https://doi.org/10.3390/plants12030608

AMA Style

Zhang H, Cao Y, Wang Z, Ye M, Wu R. Functional Mapping of Genes Modulating Plant Shade Avoidance Using Leaf Traits. Plants. 2023; 12(3):608. https://doi.org/10.3390/plants12030608

Chicago/Turabian Style

Zhang, Han, Yige Cao, Zijian Wang, Meixia Ye, and Rongling Wu. 2023. "Functional Mapping of Genes Modulating Plant Shade Avoidance Using Leaf Traits" Plants 12, no. 3: 608. https://doi.org/10.3390/plants12030608

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop