<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="1.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title-group>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2026.1747538</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Research</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Quantifying the effects of cultivation mode and sprouting stage on tea bud morphology and chemical quality</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Zhu</surname><given-names>Wei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/3280112/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Funding acquisition" vocab-term-identifier="https://credit.niso.org/contributor-roles/funding-acquisition/">Funding acquisition</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="visualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/visualization/">Visualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Formal analysis" vocab-term-identifier="https://credit.niso.org/contributor-roles/formal-analysis/">Formal analysis</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Yang</surname><given-names>Jiayi</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="visualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/visualization/">Visualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Xiao</surname><given-names>Suqing</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Hu</surname><given-names>Wenjing</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="resources" vocab-term-identifier="https://credit.niso.org/contributor-roles/resources/">Resources</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Tian</surname><given-names>Jiamin</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Zeng</surname><given-names>Yuchen</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/3282637/overview"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Huang</surname><given-names>Zhiming</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Yuan</surname><given-names>Ziyi</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &amp; editing</role>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label><institution>Jingdezhen Key Laboratory of Tea Plant Germplasm Resources and Ecological Cultivation, Jingdezhen University</institution>, <city>Jingdezhen</city>,&#xa0;<country country="cn">China</country></aff>
<aff id="aff2"><label>2</label><institution>Ministry of Education Key Laboratory of Silviculture and Conservation, Beijing Forestry University</institution>, <city>Beijing</city>,&#xa0;<country country="cn">China</country></aff>
<aff id="aff3"><label>3</label><institution>Fuliang Changnan Tea Co., Ltd.</institution>, <city>Jingdezhen</city>,&#xa0;<country country="cn">China</country></aff>
<author-notes>
<corresp id="c001"><label>*</label>Correspondence: Wei Zhu, <email xlink:href="mailto:zhuwei@jdzu.edu.cn">zhuwei@jdzu.edu.cn</email>; <email xlink:href="mailto:18811375082@163.com">18811375082@163.com</email></corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-01-30">
<day>30</day>
<month>01</month>
<year>2026</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2026</year>
</pub-date>
<volume>17</volume>
<elocation-id>1747538</elocation-id>
<history>
<date date-type="received">
<day>16</day>
<month>11</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>01</month>
<year>2026</year>
</date>
<date date-type="rev-recd">
<day>10</day>
<month>01</month>
<year>2026</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2026 Zhu, Yang, Xiao, Hu, Tian, Zeng, Huang and Yuan.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>Zhu, Yang, Xiao, Hu, Tian, Zeng, Huang and Yuan</copyright-holder>
<license>
<ali:license_ref start_date="2026-01-30">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</license-p>
</license>
</permissions>
<abstract>
<p>While the influences of environment and genetics on plant traits are widely recognized, quantitative analysis of their relative contributions remains limited in crop systems, constraining precise environmental manipulation of agronomic traits. This study investigated the effects of cultivation mode (representing external environment) and sprouting stage (reflecting internal developmental program) on tea bud morphology and chemical quality, with particular focus on their intrinsic relationships. Through dynamic monitoring at three sprouting period across five contrasting cultivation modes in Fuliang County, Jiangxi Province, China, key morphological and chemical traits were obtained. The effects of cultivation mode and sprouting stage were assessed using linear mixed-effects models and partial &#x3c9;&#xb2;. Results revealed a consistent developmental trend over the observed sprouting period: across all cultivation modes, bud length and number increased continuously, while width decreased, and size first declined then increased. Critically, cultivation mode had a stronger influence on bud width, number, and size (effect sizes: 0.739&#x2014;0.768) than on length (0.132), which was predominantly governed by the sprouting stage (effect size: 0.812). Cultivation mode exerted a stronger effect on morphological composite traits (0.753) than on chemical composite traits (0.530). The correlations between bud length and width and between width and number shifted from positive to negative during sprouting (<italic>P</italic> &lt; 0.05), suggesting increasing resource competition among these traits. Most importantly, no significant correlations were observed between morphological traits and chemical quality at the final sprouting stage. Based on data from this single yet harvest-critical time point, our results directly challenge the traditional belief that bud size predicts tea quality. These findings provide new insights into gene-environment interactions, phenotypic plasticity, and ecological trade-offs in agroecosystems, establishing a scientific basis for precision management strategies that simultaneously optimize tea production and quality.</p>
</abstract>
<kwd-group>
<kwd>cultivation modes</kwd>
<kwd>sprouting dynamics</kwd>
<kwd>tea bud morphology</kwd>
<kwd>chemical quality</kwd>
<kwd>correlation</kwd>
</kwd-group>
<funding-group>
<funding-statement>The author(s) declared that financial support was received for this work and/or its publication. This research was funded by the Science and Technology Research Project of Jiangxi Provincial Department of Education [GJJ2402302] and the Science and Technology Development Project of Jingdezhen City [20234NY005].</funding-statement>
</funding-group>
<counts>
<fig-count count="8"/>
<table-count count="1"/>
<equation-count count="2"/>
<ref-count count="55"/>
<page-count count="16"/>
<word-count count="9353"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Crop and Product Physiology</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Plant traits represent the ultimate phenotype resulting from the interaction between genotype and environmental factors. The mechanisms underlying their formation constitute a central topic in genetics, ecology, and crop breeding (<xref ref-type="bibr" rid="B9">Chapin et&#xa0;al., 1987</xref>; <xref ref-type="bibr" rid="B31">Lynch and Walsh, 1998</xref>). Disentangling the relative contributions of intrinsic developmental programs and external environments to trait formation is not only crucial for understanding species adaptation and ecological differentiation (<xref ref-type="bibr" rid="B45">Tsukaya, 2005</xref>; <xref ref-type="bibr" rid="B7">Carbone et&#xa0;al., 2009</xref>; <xref ref-type="bibr" rid="B29">Li et&#xa0;al., 2020b</xref>), but also provides a theoretical foundation for the precise regulation of agricultural traits. However, in many crop systems, knowledge regarding the quantitative decomposition of their effects and the inherent trade-off relationships remains limited (<xref ref-type="bibr" rid="B46">Wachira et&#xa0;al., 2002</xref>; <xref ref-type="bibr" rid="B14">Dutta et&#xa0;al., 2010</xref>; Owuor et&#xa0;al., 2011; <xref ref-type="bibr" rid="B21">Hua et&#xa0;al., 2022</xref>). This significantly constrains our ability to optimize agronomic traits through targeted environmental management.</p>
<p>Tea (<italic>Camellia sinensis</italic>), as one of the world&#x2019;s most important economic crops, is cultivated in over 60 countries and regions and consumed by more than 2 billion people globally, demonstrating a profound global influence (<xref ref-type="bibr" rid="B4">International Institute for Sustainable Development, 2024</xref>). In premium tea production, harvesting standards such as &#x201c;a single bud&#x201d;, &#x201c;one bud with one leaf&#x201d;, and &#x201c;one bud with two leaves&#x201d; are widely adopted (<xref ref-type="bibr" rid="B1">Aaqil et&#xa0;al., 2023</xref>). Tea bud size is a critical agronomic and biological trait, not only directly affecting fresh leaf yield but also closely related to the processing suitability and quality of the final tea product (<xref ref-type="bibr" rid="B28">Li et&#xa0;al., 2020a</xref>; <xref ref-type="bibr" rid="B3">An et&#xa0;al., 2021</xref>). Therefore, the morphogenesis of tea buds and the formation of their chemical quality are core processes determining economic value. However, existing research has predominantly focused on the macroscopic growth rhythm of tea shoots (<xref ref-type="bibr" rid="B42">Smith et&#xa0;al., 1990</xref>; Botwright, 1997; <xref ref-type="bibr" rid="B52">Zhang et&#xa0;al., 2024a</xref>) or genetic analyses of bud static traits at specific developmental stages (<xref ref-type="bibr" rid="B50">Zhang et&#xa0;al., 2025b</xref>), while systematic investigation into the dynamic developmental patterns of tea bud morphology, their response to cultivation environments, and the correlation mechanisms between morphology and intrinsic chemical quality has long remained insufficient.</p>
<p>In recent years, researchers have begun to employ emerging technologies to explore the genetic basis of tea bud morphogenesis. For instance, <xref ref-type="bibr" rid="B50">Zhang et&#xa0;al. (2025b)</xref> utilized digital phenotyping technology to conduct association analysis on a large natural population, revealing extensive variation in tea bud size and providing preliminary insights into its genetic architecture. Such studies offer valuable perspectives for understanding the genetic potential of trait variation. However, this paradigm possesses inherent limitations: relying on static data from a single time point, it fundamentally fails to capture the dynamic trajectory and plastic responses of tea buds during continuous development. More critically, in agricultural ecosystems, phenotypes result from the continuous interaction between genetic programs and environmental signals; systematically deconstructing the effects of cultivation mode (as a manageable environmental factor) and sprouting stage (representing the internal developmental program) on trait formation, and quantitatively comparing their relative importance, is essential for elucidating genotype-by-environment (G&#xd7;E) interaction mechanisms (<xref ref-type="bibr" rid="B32">Malosetti et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B15">El-Soda et&#xa0;al., 2014</xref>). Nevertheless, in tea research, the quantitative dissection of the influence of these two factors remains an underexplored area, significantly impeding a deeper understanding of trait formation mechanisms and constraining the optimization of precision cultivation strategies.</p>
<p>The chemical quality of tea is determined by a diverse array of endogenous chemical components and is known to be intricately regulated by multiple factors, including cultivar, season, leaf position, climate, and agronomic practices (<xref ref-type="bibr" rid="B48">Wu et&#xa0;al., 2025</xref>). Numerous studies have sought to elucidate the dynamic patterns of these chemical constituents. For example, <xref ref-type="bibr" rid="B30">Li et&#xa0;al. (2016)</xref> observed that as leaves mature, carbon and nitrogen metabolic pathways undergo stage-specific shifts, leading to a significant decline in the content of nitrogen-containing compounds such as amino acids. <xref ref-type="bibr" rid="B49">Xu et&#xa0;al. (2021)</xref> further demonstrated that gallic acid, caffeine, and certain esterified catechins also decrease with increasing bud and leaf maturity, reflecting a clear developmental stage specificity. At the genetic level, <xref ref-type="bibr" rid="B53">Zhang et&#xa0;al. (2025a)</xref> found that the free amino acid content is not only influenced by season and germplasm but also exhibits genetic correlations with sprouting traits, suggesting that quality and morphological traits may share a partial genetic regulatory basis. Although these studies have enhanced our understanding of quality formation from various perspectives, none have directly and systematically examined the intrinsic relationship between the external morphology and internal chemical quality of tea buds. Consequently, the long-held traditional belief in both tea science and industry that &#x201c;larger buds signify superior quality&#x201d; still lacks rigorous scientific validation. The persistence of this critical knowledge gap significantly hinders the development of technologies for efficient quality prediction and precision breeding based on appearance traits.</p>
<p>To address the knowledge gaps identified above, this study established five tea garden cultivation modes with significant environmental gradients in Fuliang County, Jiangxi Province, China. Based on the biological principles of tea plant growth and development and the theoretical framework of secondary metabolic regulation, we conducted continuous dynamic monitoring and multi-dimensional quantitative analysis. This research aims to achieve the following three core objectives: (1) investigate the growth dynamics of tea buds under different cultivation modes, quantifying the relative effect size of cultivation mode and sprouting stage on bud morphology. (2) analyze the variation characteristics of chemical quality indicators in tea buds across different cultivation modes, quantify the effect size of cultivation methods on chemical quality, and determine whether cultivation mode has a greater influence on morphological traits or chemical quality. (3) uncover the dynamic correlations among morphological indicators of tea buds, as well as the intrinsic relationships between final bud morphology, growth rate, and chemical quality. By deconstructing the effects of cultivation mode and sprouting stage on tea bud phenotypes, this study will provide quantitative evidence for elucidating genotype-by-environment interaction mechanisms in tea plant trait formation and challenge the conventional belief that &#x201c;larger buds indicate better quality,&#x201d; thereby establishing a solid scientific foundation for developing precision agronomic strategies that simultaneously optimize tea yield and quality.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Study site and experimental tea plantation</title>
<p>The study site was located in a tea plantation in Fuliang County, Jingdezhen City, Jiangxi Province (E 116&#xb0;57&#x2032;--117&#xb0;42&#x2032;, N 28&#xb0;44&#x2032;--29&#xb0;56&#x2032;). The study area experiences a subtropical monsoon climate with distinct seasons and concurrent rain and heat. Precipitation is concentrated in spring and summer, while autumn and winter are mild. The annual average temperature is 17.36&#xb0;C, with historical extreme maximum temperatures exceeding 40&#xb0;C and extreme minimum temperatures below -5&#xb0;C. The average annual precipitation is 1763.5 mm, unevenly distributed in time and space, with the rainy season from April to June accounting for about 46% of the annual total. The terrain is predominantly low mountains and hills, with layered tea terraces, significant diurnal temperature variation, and high air humidity. The soil is primarily red and yellow earth, rich in organic matter and minerals. With a forest coverage rate exceeding 80% and minimal industrial pollution, it is a typical high-quality tea production area, providing an ideal environment for premium tea growth. The experiment was conducted from March to April 2025, encompassing five different cultivation modes, whose basic profiles are shown in <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>List of tea plantations under different cultivation modes.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Cultivation mode</th>
<th valign="middle" align="center">Abbreviation</th>
<th valign="middle" align="center">Altitude (m)</th>
<th valign="middle" align="center">Growth Environment</th>
<th valign="middle" align="center">Tea Cultivar</th>
<th valign="middle" align="center">Intercropped plants and their planting density</th>
<th valign="middle" align="center">Age (years)</th>
<th valign="middle" align="center">Total N (g&#xb7;kg<sup>&#x2212;1</sup>)</th>
<th valign="middle" align="center">pH</th>
<th valign="middle" align="center">Air temperature (&#xb0;C)</th>
<th valign="middle" align="center">Air humidity (%)</th>
<th valign="middle" align="center">Total radiation (W&#xb7;m<sup>&#x2212;2</sup>)</th>
<th valign="middle" align="center">Soil moisture (%)</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">Hilly Tea Monoculture</td>
<td valign="middle" align="center">HT</td>
<td valign="middle" align="center">92.11</td>
<td valign="middle" align="center">Hilltop</td>
<td valign="middle" align="center"><italic>Camellia sinensis</italic> (L.) O. Kuntze cv. Zhuyeqi</td>
<td valign="middle" align="center">&#x2014;</td>
<td valign="middle" align="center">6</td>
<td valign="middle" align="center">1.65 &#xb1; 0.08</td>
<td valign="middle" align="center">4.74 &#xb1; 0.04</td>
<td valign="middle" align="center">17.42 &#xb1; 3.72</td>
<td valign="middle" align="center">74.58 &#xb1; 8.20</td>
<td valign="middle" align="center">150.68 &#xb1; 57.61</td>
<td valign="middle" align="center">31.46 &#xb1; 2.29</td>
</tr>
<tr>
<td valign="middle" align="center">Hilly Tea-Forest Intercropping</td>
<td valign="middle" align="center">HTF</td>
<td valign="middle" align="center">87.77</td>
<td valign="middle" align="center">Hilltop</td>
<td valign="middle" align="center"><italic>C. sinensis</italic> (L.) O. Kuntze cv. Zhuyeqi</td>
<td valign="middle" align="center">Chinese fir, planted in shelterbelts</td>
<td valign="middle" align="center">6</td>
<td valign="middle" align="center">1.75 &#xb1; 0.12</td>
<td valign="middle" align="center">4.69 &#xb1; 0.21</td>
<td valign="middle" align="center">16.69 &#xb1; 3.71</td>
<td valign="middle" align="center">77.96 &#xb1; 7.66</td>
<td valign="middle" align="center">107.27 &#xb1; 42.23</td>
<td valign="middle" align="center">42.10 &#xb1; 1.67</td>
</tr>
<tr>
<td valign="middle" align="center">Plain Tea Monoculture</td>
<td valign="middle" align="center">PT</td>
<td valign="middle" align="center">64.53</td>
<td valign="middle" align="center">Riverbank</td>
<td valign="middle" align="center"><italic>C. sinensis</italic> (L.) O. Kuntze cv. Zhuyeqi</td>
<td valign="middle" align="center">&#x2014;</td>
<td valign="middle" align="center">6</td>
<td valign="middle" align="center">1.79 &#xb1; 0.06</td>
<td valign="middle" align="center">5.03 &#xb1; 0.05</td>
<td valign="middle" align="center">17.18 &#xb1; 3.70</td>
<td valign="middle" align="center">78.32 &#xb1; 8.45</td>
<td valign="middle" align="center">188.56 &#xb1; 75.19</td>
<td valign="middle" align="center">34.33 &#xb1; 3.02</td>
</tr>
<tr>
<td valign="middle" align="center">Plain Tea-Forest Intercropping</td>
<td valign="middle" align="center">PTF</td>
<td valign="middle" align="center">65.63</td>
<td valign="middle" align="center">Riverbank</td>
<td valign="middle" align="center"><italic>C. sinensis</italic> (L.) O. Kuntze cv. Zhuyeqi</td>
<td valign="middle" align="center">Cherry trees (canopy width approx. 1.8 m), spaced at 5 m &#xd7; 20 m</td>
<td valign="middle" align="center">6</td>
<td valign="middle" align="center">1.70 &#xb1; 0.03</td>
<td valign="middle" align="center">4.89 &#xb1; 0.05</td>
<td valign="middle" align="center">&#x2014;</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
</tr>
<tr>
<td valign="middle" align="center">Plain Tea-Soybean Intercropping</td>
<td valign="middle" align="center">PTS</td>
<td valign="middle" align="center">70.03</td>
<td valign="middle" align="center">Riverbank</td>
<td valign="middle" align="center"><italic>C. sinensis</italic> (L.) O. Kuntze cv. Zhuyeqi</td>
<td valign="middle" align="center">Soybean, interplanted within tea rows</td>
<td valign="middle" align="center">6</td>
<td valign="middle" align="center">1.83 &#xb1; 0.18</td>
<td valign="middle" align="center">5.63 &#xb1; 0.16</td>
<td valign="middle" align="center">&#x2014;</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>The meteorological data in the table represent the mean &#xb1; standard deviation during the experimental period from March to April 2025.</p></fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Microclimate monitoring in tea plantations</title>
<p>Microclimate data were collected using miniature environmental monitoring systems (ME40, Beijing Shiyu Tong Technology Co., Ltd.). One monitoring device was installed in each of the HT, HTF, and PT tea gardens to continuously track parameters such as air temperature, air humidity, total solar radiation, and topsoil moisture content at the tea canopy level. Data were recorded at 10-minute intervals. After quality control processing, daily average values were calculated. The mean values and standard deviations of each meteorological factor during the entire experimental period (March to April 2025) were compiled for subsequent analysis.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Measurement of tea bud growth dynamics</title>
<sec id="s2_3_1">
<label>2.3.1</label>
<title>Monitoring design</title>
<p><italic>Camellia sinensis</italic> (L.) O. Kuntze cv. Zhuyeqi is an intermediate-sprouting cultivar. In its region of origin, the &#x201c;one bud with three leaves&#x201d; stage typically occurs in early April (This cultivar was developed through individual plant selection from the Anhua population by the Tea Research Institute of Hunan Province; see <uri xlink:href="http://hntri.hunaas.cn/?p=30&amp;a=view&amp;r=373">http://hntri.hunaas.cn/?p=30&amp;a=view&amp;r=373</uri> for details). The period from bud dormancy release to the formation of &#x201c;one bud with three leaves&#x201d; spans approximately one month, representing the crucial sprouting phase for tea buds. To capture the dynamics during this critical period, we conducted three rounds of bud morphological monitoring, with the final round in early April coinciding with chemical quality assessment. Notably, early April (around the Qingming Festival) also marks the important agricultural milestone when local tea farmers begin harvesting the first batch of premium teas, which primarily consist of single buds. Specifically, three dynamic observations were carried out during the crucial tea bud sprouting period in 2025 on March 6 (the early stage), March 19 (the middle stage), and April 7 (the late stage, proximate to Qingming), across all experimental plots under the five cultivation modes (HT, HTF, PT, PTF, and PTS).</p>
</sec>
<sec id="s2_3_2">
<label>2.3.2</label>
<title>Quadrat setting</title>
<p>for each cultivation mode, five independent observation points (i.e., five biological replicates, n = 5 per mode) were randomly selected. A fixed standard quadrat of 30 cm &#xd7; 30 cm was established at each point, defining the experimental unit for all subsequent measurements and analyses.</p>
</sec>
<sec id="s2_3_3">
<label>2.3.3</label>
<title>Bud morphology and number</title>
<p>During each observation, the longitudinal (length, mm) and transverse (width, mm) diameters of all vigorous tea buds (with unexpanded apical leaves) within the quadrat were measured using a digital caliper (accuracy: 0.01 mm). Bud size (mm<sup>3</sup>) was expressed as the bud size index, calculated as the square of length multiplied by width. For each quadrat, mean values for length, width, and size were calculated from all measured buds. The total bud count was recorded to obtain the bud number per unit area (buds m<sup>&#x2212;2</sup>). Growth dynamics curves were constructed based on these three time-series observations.</p>
</sec>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Determination of tea chemical quality</title>
<sec id="s2_4_1">
<label>2.4.1</label>
<title>Sampling and pre-treatment</title>
<p>Sampling was conducted concurrently with the final tea bud growth measurement (April 7). Within the above mentioned fixed quadrats (in section 2.3), all disease-free and fresh tea bud samples were harvested (i.e., five biological replicates per cultivation mode, with at least 80 buds collected per replicate). The samples were immediately placed in an icebox and transported to the laboratory. They were fixed using a microwave oven (medium-high power for 2&#x2013;2.5 minutes), dried in an oven at 80&#xb0;C, ground, and sieved through a 40-mesh sieve. The processed powder was sealed and stored protected from light in a -40&#xb0;C ultra-low temperature freezer for subsequent analysis.</p>
</sec>
<sec id="s2_4_2">
<label>2.4.2</label>
<title>Determination of free amino acids and tea polyphenols</title>
<p>The contents of free amino acids and tea polyphenols were determined according to the Chinese National Standards GB/T 8314&#x2013;2013 and GB/T 8313-2018, respectively, using a UV spectrophotometer (UV-8000, Shanghai Metash Instruments CO., LTD., Shanghai, China). Samples were extracted with hot water (for amino acid determination) or 70% methanol (for tea polyphenol determination), followed by filtration or centrifugation and volume adjustment. An aliquot of the extract was then reacted with ninhydrin solution or Folin-Ciocalteu reagent, and the absorbance was measured at wavelengths of 570 nm or 765 nm, respectively. Quantification was performed using standard curves prepared with theanine or gallic acid as the standard, respectively.</p>
</sec>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Data analysis</title>
<sec id="s2_5_1">
<label>2.5.1</label>
<title>Part I: effect analysis of cultivation mode and sprouting stage on tea bud morphology and growth rate</title>
<p>To quantify the effects of cultivation mode and sprouting stage on tea bud morphology, we employed linear mixed-effects models for analysis. For the four morphological indicators of tea buds (length, width, number, and size), the model fixed effects included cultivation mode (5 levels: HT, HTF, PT, PTF, PTS), sprouting stage (3 levels: early, middle, late, set as an ordered factor), and their interaction. The random effect was quadrat, to control for the non-independence among repeated measurements from the same quadrat (n = 5). To handle temporal autocorrelation common in repeated measures data, we systematically compared six residual covariance structures: Independent (INDEP), Compound Symmetry (CS), First-Order Autoregressive (AR1), Heterogeneous Compound Symmetry (CSH), Toeplitz (TOEP), and Unstructured (UN). Specific formula expressions can be found in the supplementary material section &#x201c;1. Model Formulations&#x201d;.</p>
<p>The optimal covariance structure for each morphological indicator was selected based on the Akaike Information Criterion (AIC) (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S1</bold></xref>). Model assumption tests (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S2</bold></xref>) showed that, except for tea bud width residuals not fully meeting the normality assumption (Shapiro-Wilk <italic>p</italic> = 0.023), the residuals for the other indicators satisfied the assumptions of normality and homogeneity of variance. Given that linear mixed models are robust to slight deviations from normality, subsequent inferences were all based on the selected best models. We reported the Type III ANOVA results of the best models (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S3</bold></xref>), fixed effect coefficient estimates, standard errors and significance (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S4</bold></xref>), and residual diagnostic plots (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figures S1</bold></xref>&#x2013;<xref ref-type="supplementary-material" rid="SM1"><bold>S4</bold></xref>). Based on the model results, we used the emmeans package for <italic>post-hoc</italic> multiple comparisons, with <italic>p</italic>-values adjusted using the Tukey method. To quantify the relative importance of each factor, we calculated the relative contribution (%) of cultivation mode, sprouting stage, and their interaction to the variation of each morphological indicator, based on the <italic>F</italic>-values from the model Type III ANOVA, and calculated partial &#x3c9;&#xb2; as a standardized effect size to indicate the proportion of variance attributable to each factor after accounting for other effects in the model.</p>
<p>To analyze the dynamic growth of the morphological indicators, we calculated the growth rate for each indicator between two adjacent sprouting stages (<inline-formula>
<mml:math display="inline" id="im1"><mml:mrow><mml:mi>G</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>g</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>), i.e., the daily absolute change:</p>
<disp-formula>
<mml:math display="block" id="M1"><mml:mrow><mml:mi>G</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>s</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>g</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mi>Y</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>&#x2212;</mml:mo><mml:msub><mml:mi>Y</mml:mi><mml:mi>t</mml:mi></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>&#x394;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:mfrac></mml:mrow></mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im2"><mml:mi>Y</mml:mi></mml:math></inline-formula> represents the measured value of the morphological indicator, and <inline-formula>
<mml:math display="inline" id="im3"><mml:mrow><mml:mi>&#x394;</mml:mi><mml:mi>t</mml:mi></mml:mrow></mml:math></inline-formula> is the number of days between adjacent stages (early to middle: 13 days, middle to late: 19 days). Subsequently, we used paired <italic>t</italic>-tests to compare the differences in growth rates for each indicator between the &#x201c;early&#x2013;middle&#x201d; and &#x201c;middle&#x2013;late&#x201d; stages (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S5</bold></xref>, <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S5</bold></xref>).</p>
<p>Calculation of overall growth rate: To obtain a single indicator characterizing the growth intensity over the entire sprouting process (early to late) for subsequent correlation analysis in Part II, we calculated the time-weighted average overall growth rate (<inline-formula>
<mml:math display="inline" id="im4"><mml:mrow><mml:mi>O</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>G</mml:mi><mml:mi>R</mml:mi></mml:mrow></mml:math></inline-formula>). This method considered the influence of different stage durations, avoiding bias from simple arithmetic averaging. The calculation formula is as follows:</p>
<disp-formula>
<mml:math display="block" id="M2"><mml:mrow><mml:mi>O</mml:mi><mml:mi>v</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mi>l</mml:mi><mml:mi>G</mml:mi><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>G</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>l</mml:mi><mml:mi>y</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo>&#xd7;</mml:mo><mml:mi>&#x394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>l</mml:mi><mml:mi>y</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>+</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>G</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>&#xd7;</mml:mo><mml:mi>&#x394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>&#x394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>l</mml:mi><mml:mi>y</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>&#x394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:mfrac></mml:mrow></mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im5"><mml:mrow><mml:mi>G</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>l</mml:mi><mml:mi>y</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula>
<mml:math display="inline" id="im6"><mml:mrow><mml:mi>G</mml:mi><mml:msub><mml:mi>R</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> are the growth rates for the early-middle and middle-late stages, respectively, <inline-formula>
<mml:math display="inline" id="im7"><mml:mrow><mml:mi>&#x394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>e</mml:mi><mml:mi>a</mml:mi><mml:mi>r</mml:mi><mml:mi>l</mml:mi><mml:mi>y</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> = 13 days, <inline-formula>
<mml:math display="inline" id="im8"><mml:mrow><mml:mi>&#x394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>d</mml:mi><mml:mo>&#x2212;</mml:mo><mml:mi>l</mml:mi><mml:mi>a</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> = 19 days, <inline-formula>
<mml:math display="inline" id="im9"><mml:mrow><mml:mi>&#x394;</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi><mml:mi>t</mml:mi><mml:mi>a</mml:mi><mml:mi>l</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> = 32 days. The calculation was based on the raw data from each quadrat, ultimately obtaining the overall growth rates for length, width, number, and size for each quadrat replicate over the entire sprouting period.</p>
</sec>
<sec id="s2_5_2">
<label>2.5.2</label>
<title>Part II: effect analysis of cultivation mode on chemical quality, final morphology, and growth rate</title>
<p>To assess the comprehensive impact of cultivation mode on the chemical quality of mature tea buds, their final morphology, and the overall growth rate during the sprouting period, we established a series of linear mixed models. The model fixed effect was cultivation mode, and the random effect was quadrat (n = 5), used to analyze the following three categories of indicators: 1. Chemical quality indicators: Free amino acid content, tea polyphenol content, and the polyphenol-to-amino acid ratio determined at the late sprouting stage. 2. Final morphological indicators: Tea bud length, width, number, and size measured at the late sprouting stage. 3. Overall growth rate indicators: The overall growth rates for each morphological indicator over the entire sprouting period, calculated based on the method from Part I.</p>
<p>All models underwent diagnostic checks for residual normality (Shapiro-Wilk test) and homogeneity of variance (Levene&#x2019;s test) (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S6</bold></xref>), and residual diagnostic plots were provided (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figures S6</bold></xref>-<xref ref-type="supplementary-material" rid="SM1"><bold>S16</bold></xref>). Based on models meeting the basic assumptions, we performed Type III ANOVA to test the main effect of cultivation mode (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S7</bold></xref>) and reported the model parameter estimates table (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S8</bold></xref>). Simultaneously, we calculated the partial &#x3c9;&#xb2; for cultivation mode&#x2019;s effect on each indicator as a measure of its effect size.</p>
<p>To comprehensively evaluate the overall influence intensity of cultivation mode on each set of indicators, we performed Principal Component Analysis (PCA). PCA was conducted separately on the chemical indicators, morphological indicators, and growth rate indicator, and the scores of the first principal component (PC1) were extracted. These scores were named the &#x201c;Chemical Composite Index,&#x201d; &#x201c;Morphological Composite Index,&#x201d; and &#x201c;Growth Composite Index,&#x201d; respectively, summarizing most of the variation within each category (variance explained &gt; 60%, <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Tables S9</bold></xref>, <xref ref-type="supplementary-material" rid="SM1"><bold>S10</bold></xref>). Subsequently, we again used linear mixed models (fixed effect: cultivation mode, random effect: quadrat) to analyze the effect of cultivation mode on these composite indices (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Tables S7</bold></xref>, <xref ref-type="supplementary-material" rid="SM1"><bold>S8</bold></xref>, <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figures S17</bold></xref>, <xref ref-type="supplementary-material" rid="SM1"><bold>S19</bold></xref>) and calculated the partial &#x3c9;&#xb2;, thereby comparing the differences in the influence of cultivation mode at a comprehensive, dimension-reduced level.</p>
</sec>
<sec id="s2_5_3">
<label>2.5.3</label>
<title>Part III: inter-indicator correlations</title>
<p>To explore the intrinsic relationships among various indicators and their dynamic changes during the sprouting process, we conducted systematic non-parametric correlation analyses. First, we analyzed the pairwise correlations among the three morphological indicators (length, width, number) within each sprouting stage. Second, we performed a comprehensive correlation analysis on all indicators involved in Part II (3 late-stage chemical quality indicators, 4 late-stage morphological indicators, 4 overall growth rate indicators).</p>
<p>As some indicator data distributions did not conform to normality based on the Shapiro-Wilk test (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Tables S11</bold></xref>, <xref ref-type="supplementary-material" rid="SM1"><bold>S12</bold></xref>), and scatter plots suggested potential non-linear trends for some variable pairs, we chose Spearman&#x2019;s rank correlation coefficient (<inline-formula>
<mml:math display="inline" id="im10"><mml:mrow><mml:msub><mml:mi>&#x3c1;</mml:mi><mml:mi>s</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) for assessment. This method does not rely on data distribution assumptions and is robust to outliers.</p>
<p>To test whether the correlation coefficient for the same pair of morphological indicators (e.g., length <italic>vs</italic>. width) differed statistically between different sprouting stages (e.g., early <italic>vs</italic>. middle), we employed a non-parametric Bootstrap permutation test. The specific steps were: 1) Combine the data from the two stages; 2) Perform random resampling with replacement to create Bootstrap samples (1000 repetitions); 3) For each repetition, calculate the difference in <inline-formula>
<mml:math display="inline" id="im11"><mml:mrow><mml:msub><mml:mi>&#x3c1;</mml:mi><mml:mi>s</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> between the two stages based on the Bootstrap samples; 4) Calculate the 95% confidence interval based on the distribution of the differences. If this interval did not contain 0, the correlation coefficients for the two stages were considered significantly different (<italic>P</italic> &lt; 0.05; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S13</bold></xref>). This method avoids the strict data distribution requirements of parametric tests.</p>
<p>All statistical analyses and graphics were generated using R software (v. 4.5.2, <xref ref-type="bibr" rid="B36">R Core Team, 2025</xref>). The following packages were used: nlme, lme4, lmerTest, emmeans, tidyverse, ggplot2, MuMIn, car, performance, effectsize, FactoMineR, factoextra, psych, and cocor.</p>
</sec>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Tea bud growth followed similar dynamic trends across cultivation modes</title>
<p>The dynamic trends of various morphological indicators of tea buds were generally consistent across different cultivation modes (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>). Tea bud length and number continuously increased throughout the sprouting process (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f1"><bold>Figures&#xa0;1A, C</bold></xref>; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S3</bold></xref>). The relative growth for bud length was 20% and 106% during the early-to-middle and middle-to-late stages, respectively. The corresponding relative growth for bud number were 74% and 44%. In contrast, bud width gradually decreased (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1B</bold></xref>; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S3</bold></xref>), with a significant reduction of -21% during the early-to-middle stage (<italic>P</italic> &lt; 0.05) and no significant change thereafter (<italic>P</italic> &gt; 0.05). Bud size overall showed an initial decrease (relative change: -24%) followed by an increase (relative change: 92%) (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1D</bold></xref>; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S3</bold></xref>).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Dynamics of tea bud morphological traits across three sprouting stages under different cultivation modes (n = 5). Colored solid lines represent smoothed trends (loess fitting) for each cultivation mode; shaded areas indicate 95% confidence intervals. The black dashed line represents the overall mean trend across all cultivation modes. Different lowercase letters above the shaded areas indicate significant differences among stages based on Tukey&#x2019;s HSD post&#x2212;hoc tests (<italic>P</italic> &lt; 0.05). Abbreviations for cultivation modes (HT, HTF, PT, PTF, PTS) are defined in <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g001.tif">
<alt-text content-type="machine-generated">Line graphs illustrating bud characteristics at different sprouting stages. Graph A shows bud length increasing from early to late stages. Graph B depicts bud width slightly decreasing. Graph C displays bud number increasing over time. Graph D illustrates bud size slightly decreasing before increasing. Different colored lines represent various cultivation modes: HT, HTF, PT, PTF, and PTS. Shaded areas indicate variability.</alt-text>
</graphic></fig>
<p>Furthermore, the growth rates for tea bud length, width, and size during the middle-to-late stage were significantly higher than those during the early-to-middle stage (<italic>P</italic> &lt; 0.05, <xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2</bold></xref>), whereas the growth rate for bud number did not differ significantly between the two stages (<italic>P</italic> &gt; 0.05, <xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2</bold></xref>). This pattern was consistent across most cultivation modes (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S5</bold></xref>). Across all cultivation modes, the growth rates for length and number were consistently positive, whereas the rate for width remained negative. The growth rate for bud size shifted from negative to positive. These patterns of growth rates corroborate the consistent overall developmental trends described above (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S5</bold></xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Comparison of daily growth rates of tea bud morphological traits between two successive sprouting periods across all cultivation modes (n = 25). Growth rates were calculated as the daily absolute change of each trait from early&#x2212;to&#x2212;mid (13 d) and mid&#x2212;to&#x2212;late (19 d) stages. Boxes show the interquartile range (IQR), horizontal lines denote medians, and points represent individual observations. Asterisks above the boxes indicate significance levels from paired&#x2212;sample <italic>t</italic>-tests comparing the two periods (*<italic>, P</italic> &lt; 0.05; ns, not significant).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g002.tif">
<alt-text content-type="machine-generated">Four box plots compare bud growth rates between two sprouting stages: early to mid (blue) and mid to late (orange). The graphs show bud length, width, size, and number. Significant differences are marked with red asterisks except for bud number, marked as not significant. Data points are scattered around the boxes.</alt-text>
</graphic></fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Differential effects of cultivation modes on tea bud morphology and chemical quality</title>
<p>Cultivation mode had no significant effect on tea bud length at any of the three observed sprouting stages (<italic>P</italic> &gt; 0.05; <xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3A</bold></xref>). In contrast, it exhibited a significant regulatory effect on bud width and size (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B, D</bold></xref>). Specifically, Hilly Tea Monoculture (HT) and Plain Tea-Soybean Intercropping (PTS) maintained larger bud width and size throughout the sprouting period, with their average values across stages being significantly higher by 24%&#x2013;109% compared to other modes (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B, D</bold></xref>). The response of bud number to cultivation mode varied depending on the sprouting stage (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3C</bold></xref>). No cultivation mode showed a clear advantage at the early stage. However, by the middle and late stages, Plain Tea Monoculture (PT) and Plain Tea-Forest Intercropping (PTF) supported significantly higher densities than other modes, with this advantage further increasing to 72% at the late stage (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3C</bold></xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Effects of different cultivation modes on tea bud length, width, size, and number at different sprouting stages (n = 5). Different lowercase letters indicate significant differences among cultivation modes within the same sprouting stage according to Tukey&#x2019;s HSD post&#x2212;hoc test (<italic>P</italic> &lt; 0.05) following linear mixed&#x2212;effects models (fixed effect: cultivation &#xd7; stage; random effect: quadrat). Absence of letters indicates no significant differences among modes at <italic>P</italic> &gt; 0.05. Error bars represent standard errors of the mean. Cultivation mode abbreviations (HT, HTF, PT, PTF, PTS) are detailed in <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g003.tif">
<alt-text content-type="machine-generated">Bar charts labeled A to D compare bud characteristics across sprouting stages and cultivation modes. Chart A shows bud length increasing in the late stage. Chart B displays variable bud width. Chart C indicates an increase in bud number during the late stage. Chart D depicts varying bud size, with early and late stages having larger sizes. Cultivation modes include HT, HTF, PT, PTF, and PTS, each represented by different colors. Error bars indicate variability, and letters denote statistical significance.</alt-text>
</graphic></fig>
<p>Regarding chemical quality, cultivation mode significantly affected all three chemical indicators (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f4"><bold>Figures&#xa0;4A&#x2013;C</bold></xref>). Specifically, PT and PTS were associated with a &#x201c;refreshing&#x201d; quality profile, characterized by the highest free amino acid content (4.0%; <xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>) and the lowest tea polyphenol content and polyphenol-to-amino acid ratio (16.4% and 4.26; <xref ref-type="fig" rid="f4"><bold>Figures&#xa0;4B, C</bold></xref>). Conversely, HT and Hilly Tea-Forest Intercropping (HTF) were linked to a &#x201c;mellow and robust&#x201d; profile, exhibiting the lowest free amino acid content (2.5%; <xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>) and the highest tea polyphenol content and polyphenol-to-amino acid ratio (18.0% and 7.51; <xref ref-type="fig" rid="f4"><bold>Figures&#xa0;4B, C</bold></xref>). The chemical indicators for PTF were intermediate between these two distinct categories.</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Effects of different cultivation modes on tea leaf chemical composition (n = 5). Different lowercase letters indicate significant differences among cultivation modes according to Tukey&#x2019;s HSD post&#x2212;hoc test (<italic>P</italic> &lt; 0.05) following linear mixed&#x2212;effects models (fixed effect: cultivation mode; random effect: quadrat). Error bars represent standard errors of the mean. Cultivation mode abbreviations (HT, HTF, PT, PTF, PTS) are defined in <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g004.tif">
<alt-text content-type="machine-generated">Bar charts labeled A, B, and C compare different cultivation modes (HT, HTF, PT, PTF, PTS) with respect to free amino acids, tea polyphenols, and phenolic/amino ratio. Chart A shows PTS with the highest free amino acids. Chart B shows HT has the highest tea polyphenols. Chart C shows HT and HTF with the highest phenolic/amino ratio. Error bars indicate variability.</alt-text>
</graphic></fig>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Sprouting stage predominantly drove bud length, while cultivation mode more strongly influenced width and number</title>
<p>Except for its non-significant effect on bud length, cultivation mode, sprouting stage, and their interaction all had significant effects on the various tea bud morphological indicators (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5</bold></xref>; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S3</bold></xref>). Specifically, the sprouting stage accounted for a substantially higher proportion of variance in bud length variation than cultivation mode did (82.9% <italic>vs</italic>. 3.5%; <xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>). Conversely, bud size was more sensitive to cultivation mode than to sprouting stage (66.6% <italic>vs</italic>. 19.5%; <xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>). The relative contributions of these two factors to bud width variation were similar (40.7% <italic>vs</italic>. 43.1%; <xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>). Notably, in the model for bud number, the interaction effect accounted for a greater proportion of variance (59.2%) than either factor alone (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Contribution and effect size of cultivation mode, sprouting stage, and their interaction on tea bud morphological traits (n = 5). <bold>(A)</bold> Stacked bars show the relative contribution (%) of each factor to the total variation, calculated based on model III ANOVA from linear mixed&#x2212;effects models (fixed effects: cultivation &#xd7; stage; random effect: quadrat). <bold>(B)</bold> Grouped bars represent partial &#x3c9;&#xb2; effect sizes for each factor; asterisks denote that the factor had a significant effect on the corresponding indicator (*, <italic>P</italic> &lt; 0.05; ns, not significant).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g005.tif">
<alt-text content-type="machine-generated">Bar charts labeled A and B. Chart A shows the relative contribution of cultivation mode, sprouting stage, and interaction to morphological traits: bud length, bud width, bud number, and bud size. Chart B displays effect sizes in partial omega squared for the same traits, with significance indicated by asterisks.</alt-text>
</graphic></fig>
<p>Effect size analysis revealed that cultivation mode had the smallest effect on bud length (partial &#x3c9;&#xb2; = 0.132), while its effects on the other three morphological indicators were similar (0.739&#x2013;0.768) (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>). Sprouting stage had the largest effect on bud length (0.812), followed by its effect on bud width (0.676), and its effects on bud number and size were similar (0.341&#x2013;0.349) (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>). The interaction effect was greatest on bud number (0.735), which was 2&#x2013;4 times larger than its effects on the other three indicators (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>).</p>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Cultivation mode had a stronger effect on morphological than on chemical quality traits</title>
<p>Cultivation mode had significant effects on the chemical and morphological indicators measured at the late sprouting stage, as well as on the overall growth rates of tea buds during the entire sprouting period (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6</bold></xref>; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Table S7</bold></xref>). Among these, cultivation mode exerted the strongest influence on late-stage bud number and its overall growth rate (partial &#x3c9;&#xb2;: 0.831&#x2013;0.864), followed by late-stage bud width and its growth rate (0.601&#x2013;0.620). Its effects on late-stage free amino acid content and the polyphenol-to-amino acid ratio ranked third (0.508&#x2013;0.534) (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6A</bold></xref>). The effect of cultivation mode on the overall growth rate of length (0.500) was greater than its effects on late-stage bud size and its overall growth rate (0.403&#x2013;0.490) (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6A</bold></xref>). Late-stage bud length and tea polyphenol content were the least responsive to cultivation mode (0.360&#x2013;0.397) (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6A</bold></xref>). Overall, cultivation mode had a greater effect on the composite morphological index (0.753) than on the composite growth rate index (0.699) or the composite chemical quality index (0.530) (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6B</bold></xref>).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Effect sizes of cultivation modes on individual and composite indicators of tea chemistry, morphology, and growth (n = 5). <bold>(A)</bold> Partial &#x3c9;&#xb2; (partial omega&#x2212;squared) effect sizes for individual indicators. <bold>(B)</bold> Partial &#x3c9;&#xb2; effect sizes for the first principal component (PC1) of each composite indicator group. Both panels present effect sizes estimated from separate linear mixed&#x2212;effects models (fixed effect: cultivation mode; random effect: quadrat). Asterisks denote that the cultivation mode had a significant effect on the corresponding indicator (*, <italic>P</italic> &lt; 0.05; ns, not significant).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g006.tif">
<alt-text content-type="machine-generated">Bar charts labeled A and B show effect sizes of various indicators. Chart A lists indicators like bud number and free amino acids, with effect sizes ranging from 0.360 to 0.864. Chart B displays composite indicators&#x2014;Chemical PC1 at 0.530, Growth PC1 at 0.699, and Morphological PC1 at 0.753. Color-coded groups are chemical, growth, and morphological indicators.</alt-text>
</graphic></fig>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>Correlations among morphological traits shifted dynamically across sprouting stages</title>
<p>Bud length and width were significantly positively correlated at the early sprouting stage (r = 0.58, <italic>P</italic> &lt; 0.05) but showed no significant correlation at the middle or late stages (<italic>P</italic> &gt; 0.05). The correlation coefficient between length and width shifted from positive to negative from the early to the late stage (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7A</bold></xref>). Bud length and number were not significantly correlated at any of the three stages. Their correlation coefficient first decreased and then increased during sprouting, with the coefficients at the early and late stages both being significantly different from that at the middle stage (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7B</bold></xref>). Bud width and number were significantly negatively correlated at the late stage (r = -0.81, <italic>P</italic> &lt; 0.05), but not at the early or middle stages (<italic>P</italic> &gt; 0.05). The correlation coefficient between width and number also declined from positive to negative from the early to the late stage (<italic>P</italic> &lt; 0.05; <xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7C</bold></xref>).</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Changes in correlation between tea bud morphological traits across three sprouting stages (n = 5 per stage). Each panel shows the Spearman correlation coefficient (<italic>r</italic><sub>s</sub>) for a specific trait pair. Points represent <italic>r</italic><sub>s</sub> values with 95% confidence intervals; asterisks indicate correlation significance (*, <italic>P</italic> &lt; 0.05; ns, not significant). Different lowercase letters above points denote statistically significant differences between stages, as determined by non&#x2212;parametric bootstrap permutation tests for the same trait pair. Horizontal dashed line denotes zero correlation.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g007.tif">
<alt-text content-type="machine-generated">Line graphs showing Spearman correlation coefficients at different sprouting stages. Panel A depicts Length-Width: positive in early stage, decreasing to negative. Panel B shows Length-Number: positive early and late, negative in middle. Panel C illustrates Width-Number: positive early, strongly negative late. Significance markers (*) are included.</alt-text>
</graphic></fig>
</sec>
<sec id="s3_6">
<label>3.6</label>
<title>Chemical quality was not significantly correlated with final bud morphology or growth rates</title>
<p>Systematic correlations were observed among the chemical indicators (<xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8</bold></xref>): free amino acids were significantly negatively correlated with tea polyphenols (r = -0.69;<italic>P</italic> &lt; 0.05); the negative correlation between the polyphenol-to-amino acid ratio and free amino acids (r = -0.97;<italic>P</italic> &lt; 0.05) was stronger than its positive correlation with tea polyphenols (r = 0.81;<italic>P</italic> &lt; 0.05). Notably, no significant correlations were detected between any chemical indicators and the tea bud morphological indicators measured at the final stage or growth rates throughout the sprouting period (<xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8</bold></xref>). Further analysis revealed a positive correlation between the growth rates of bud number and length (r = 0.41; <italic>P</italic> &lt; 0.05), while the growth rate of number was negatively correlated with both the late-stage bud width (r = -0.73; <italic>P</italic> &lt; 0.05) and the growth rate of bud size (r = -0.70; <italic>P</italic> &lt; 0.05). The growth rate of bud size was positively correlated with the growth rates of both length (r = 0.59; <italic>P</italic> &lt; 0.05) and width (r = 0.83; <italic>P</italic> &lt; 0.05) (<xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8</bold></xref>).</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Spearman correlation analysis among chemical quality, final morphology, and growth rates (n = 5). Cells display correlation coefficients (<italic>r</italic><sub>s</sub>) with significance stars (*, <italic>P</italic> &lt; 0.05; **, <italic>P</italic> &lt; 0.01; <bold>***,</bold><italic>P</italic> &lt; 0.001; ns, not significant). Red and blue shades represent positive and negative correlations, respectively.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-17-1747538-g008.tif">
<alt-text content-type="machine-generated">Correlation matrix heatmap showing relationships between various growth and biochemical measures, such as size growth, bud dimensions, phenol/amino acid ratio, tea polyphenols, and free amino acids. Colors range from blue for negative correlations to red for positive ones, with correlation values overlaid on the cells.</alt-text>
</graphic></fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<sec id="s4_1">
<label>4.1</label>
<title>Inherent program of bud morphological development and its phenotypic plasticity</title>
<p>As the primary harvesting target of premium teas, the dynamic morphological development of tea buds directly affects the yield of fresh leaves and is a crucial objective of tea garden agronomic management. Within the observational framework of three time points in this study, we found that despite the diversity in cultivation modes, the morphological development of tea buds exhibits an inherent conservatism: their length and number continuously increased throughout the sprouting period, width gradually decreased, while bud size showed an initial decrease followed by an increase. This overall trend remained consistent across different cultivation modes (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>; <xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S5</bold></xref>). This consistent pattern indicates that tea bud development follows an inherent, genetically programmed blueprint (<xref ref-type="bibr" rid="B11">De Costa et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B51">Zhang et&#xa0;al., 2016</xref>, <xref ref-type="bibr" rid="B50">2025b</xref>). However, in-depth analysis revealed significant differences in the responses of different morphological traits to external cultivation environments, uncovering the complex interaction between genetic control and agronomic plasticity.</p>
<p>Specifically, the development of tea bud length exhibited strong genetic constraint. This trait had an extremely high effect size for the sprouting stage (0.812) and a very low effect size for cultivation mode (0.132) (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>), indicating that its orderly elongation is a core event in the sprouting process, primarily driven by internal developmental signals, with limited scope for external cultivation regulation. In contrast, the development of tea bud width, size, and number showed higher phenotypic plasticity. Their response to cultivation mode was stronger than that of bud length (effect sizes 0.739&#x2014;0.768, <xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>). Previous research has mostly focused on the genetic analysis of static morphology; for example, <xref ref-type="bibr" rid="B50">Zhang et&#xa0;al. (2025b)</xref> revealed key genes regulating tea bud size through digital phenotyping, GWAS, and transcriptomic analysis. Based on dynamic observations at limited time points, this study preliminarily outlines the plasticity patterns of tea bud morphological development, suggesting that width and number may be key traits responsive to cultivation regulation. This finding provides initial evidence and potential targets for considering tea bud width and number as key plastic traits for agronomic intervention, although the underlying patterns and mechanisms require further validation through higher-frequency continuous monitoring in the future.</p>
<p>During the early sprouting stage, tea bud length and width exhibited a significant positive correlation. Although the relationships between the other two pairs of morphological indicators did not reach statistical significance, they also showed positive associations (<xref ref-type="fig" rid="f7"><bold>Figures&#xa0;7A&#x2013;C</bold></xref>). This suggests that under conditions of sufficient resources in the initial sprouting phase, there is synergistic resource allocation among different morphological dimensions. As sprouting progressed, tea bud length and width, as well as width and number, showed negative correlations at the late sprouting stage, and these correlation coefficients were significantly different from those observed at the early stage (<xref ref-type="fig" rid="f7"><bold>Figures&#xa0;7A, C</bold></xref>). This dynamic pattern clearly reveals a shift in tea bud developmental strategy: from synergistic growth in the early stage to a gradual shift towards competitive allocation of limited resources in the late stage between individual lateral expansion (width) and population density (number), as well as longitudinal extension (length). It is noteworthy that length and number showed no significant correlation throughout all periods (<xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7B</bold></xref>), while a significant positive correlation existed between the growth rates of length and number (<xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8</bold></xref>), suggesting a synergistic promotion between vertical elongation and population expansion. This finding aligns well with the classic plant size&#x2013;density relationship theory in ecology and the dynamic optimization theory of plant growth (<xref ref-type="bibr" rid="B6">Callaway and Walker, 1997</xref>; <xref ref-type="bibr" rid="B22">Iwasa, 2000</xref>), and systematically verifies for the first time the existence of such a trade-off mechanism during tea bud sprouting. This intrinsic trade-off relationship provides a key mechanistic perspective for understanding the dynamic plasticity of morphological construction in tea plants.</p>
<p>Different cultivation modes possibly achieve differentiated stage-specific regulation of these plastic traits by altering the source-sink relationship and micro-environment of the tea plants. For instance, the Hilly Tea Monoculture (HT) and Plain Tea-Soybean Intercropping (PTS) modes maintained larger bud width and size throughout the sprouting period (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B, D</bold></xref>). The width advantage in the HT mode might benefit from its relatively higher temperature and moderate radiation (<xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>), conditions considered conducive to promoting cell lateral expansion (<xref ref-type="bibr" rid="B5">Brown and Rickless, 1949</xref>; <xref ref-type="bibr" rid="B8">Casal and Balasubramanian, 2019</xref>). The advantage of the PTS mode likely stems from the unique biological nitrogen fixation of leguminous crops, providing a sufficient nitrogen nutrition base for the lateral growth of tea buds (<xref ref-type="bibr" rid="B13">Duan et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B34">Pokharel et&#xa0;al., 2023</xref>). On the other hand, the bud number advantage shown by Plain Tea Monoculture (PT) and Plain Tea-Forest Intercropping (PTF) in the mid-to-late stages (<xref ref-type="fig" rid="f3"><bold>Figure&#xa0;3C</bold></xref>) might be related to the promotion of axillary bud differentiation by higher solar radiation in the PT mode (<xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref> shows its highest solar radiation) (<xref ref-type="bibr" rid="B27">Leduc et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B40">Samuoliene et&#xa0;al., 2017</xref>).</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Chemical quality divergence and potential micro-environmental drivers</title>
<p>Previous research generally holds that high-altitude environments, due to higher air humidity and diffuse light, usually favor increased free amino acid content, reduced tea polyphenol levels, and a lower polyphenol-to-amino acid ratio in tea, thereby forming a &#x201c;refreshing&#x201d; quality style (<xref ref-type="bibr" rid="B19">Han et&#xa0;al., 2017</xref>; Li et&#xa0;al., 2022; <xref ref-type="bibr" rid="B37">Ren et&#xa0;al., 2025</xref>). However, this study observed an opposite trend: plain tea gardens (PT/PTS) exhibited typical &#x201c;refreshing&#x201d; quality characteristics, while hilly tea gardens (HT/HTF) leaned towards a &#x201c;mellow and robust&#x201d; style. In this study, the core driver of chemical quality differentiation lies in the significant difference in free amino acid content, directly reflected in the significant change of the polyphenol-to-amino acid ratio (<xref ref-type="fig" rid="f4"><bold>Figures&#xa0;4A, C</bold></xref>). In contrast, tea polyphenol content remained relatively stable across treatments (<xref ref-type="fig" rid="f4"><bold>Figures&#xa0;4B</bold></xref>, <xref ref-type="fig" rid="f6"><bold>6A</bold></xref>), indicating that under the conditions of this study, the difference in tea chemical quality was primarily dominated by nitrogen metabolism products&#x2014;especially the content of free amino acids&#x2014;rather than the carbon metabolism product, tea polyphenols.</p>
<p>This free amino acid-dominated quality differentiation pattern profoundly reflects the differential regulation of tea plant nitrogen metabolism pathways by the micro-environments created by different cultivation modes. In the plain modes (PT/PTS), milder light, higher air humidity, and good soil moisture conditions (<xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>) collectively created a micro-environment favorable for nitrogen assimilation and amino acid accumulation (<xref ref-type="bibr" rid="B23">Jeyaramraja et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B26">Kc et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B35">Qifang et&#xa0;al., 2025</xref>). Among them, the PTS mode, further enhanced by the biological nitrogen fixation of leguminous crops, significantly increased rhizosphere nitrogen availability, strengthening its &#x201c;refreshing&#x201d; quality characteristics (<xref ref-type="bibr" rid="B39">Ruan et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B13">Duan et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B47">Wang et&#xa0;al., 2022</xref>). In contrast, amino acid accumulation was generally limited in the hilly modes (HT/HTF), but the limiting mechanisms differed: the lower soil water content (31.46%) in the HT mode might induce water stress, inhibiting the transport of amino acids from roots to new shoots (<xref ref-type="bibr" rid="B18">Good and Zaplachinski, 1994</xref>; <xref ref-type="bibr" rid="B23">Jeyaramraja et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B35">Qifang et&#xa0;al., 2025</xref>). Although the HTF mode had higher soil moisture and lower radiation, its lowest average temperature (16.69&#xb0;C) was below the optimal growth temperature range for tea plants (20&#x2013;25&#xb0;C), significantly inhibiting the efficiency of theanine synthesis in roots and its transport to the apical parts (<xref ref-type="bibr" rid="B48">Wu et&#xa0;al., 2025</xref>), reflecting that water and temperature factors differentially dominated the suppression of nitrogen metabolism in different hilly modes.</p>
<p>This study further quantified, through effect sizes, that the regulatory effect of cultivation mode on the final morphology of tea buds monitored at the late sprouting stage (0.753) was greater than its effect on the simultaneously measured chemical quality indicators (0.530). Even the response of morphological growth rates (0.699) to cultivation mode was stronger than that of chemical quality (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6B</bold></xref>). This result clearly demonstrates that, on an agronomic timescale, the external morphological traits of tea buds exhibit higher cultivation controllability than their internal chemical quality traits.</p>
<p>To cope with unstable environments, plants can adapt by adjusting their morphological structure or by rapidly altering their metabolite composition (<xref ref-type="bibr" rid="B12">Dong and Lin, 2021</xref>; <xref ref-type="bibr" rid="B54">Zhang et&#xa0;al., 2022</xref>). However, these two adaptive strategies differ fundamentally in the underlying mechanisms of their plasticity. The production and turnover of secondary metabolites can be very rapid; their biochemical response to and cessation following environmental signals are far more sensitive and reversible than morphological adjustments (<xref ref-type="bibr" rid="B33">Metlen et&#xa0;al., 2009</xref>). In contrast, the morphogenesis of leaves or buds involves complex developmental programs such as cell division, expansion, and tissue differentiation. Changes in morphology are often irreversible or accumulate slowly (<xref ref-type="bibr" rid="B16">Fritz et&#xa0;al., 2018</xref>). Therefore, although metabolic traits are environmentally sensitive, their strong dynamic fluctuations make it difficult to achieve stable, directional, and cumulative changes through single or short-term agronomic practices. In contrast, once morphological traits undergo directional changes in response to environmental influences, they are more easily fixed and observed, thus exhibiting greater potential for agronomic manipulation.</p>
<p>This study is the first to quantitatively verify this principle within the tea plant system, revealing that the shaping effect of cultivation mode on morphology (particularly width and number) is significant and stable. Although environmental factors also drive changes in key tea metabolites (e.g., tea polyphenols, amino acids) (<xref ref-type="bibr" rid="B48">Wu et&#xa0;al., 2025</xref>), their responses are likely more susceptible to interference from multiple factors such as transient climatic fluctuations and developmental stage transitions. This may result in their effect size in response to cultivation mode being lower than that of morphological traits at a given point in time. This understanding not only deepens the comprehension of the hierarchical regulation of &#x201c;morphology-metabolism&#x201d; in tea garden ecosystems but also provides a key theoretical basis for prioritizing the directional optimization of tea bud appearance traits through precise microenvironment management (e.g., targeted regulation of light and temperature resources via intercropping or shading). Future research could further integrate continuous time-series dynamic monitoring to analyze the differential response trajectories of morphological and metabolic plasticity over longer cycles, thereby refining the process-based theoretical system for precision cultivation.</p>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Stage-specific decoupling between morphology and chemistry</title>
<p>Through dynamic monitoring, this study found no significant correlations between all chemical quality indicators (free amino acids, tea polyphenols, and the polyphenol-to-amino acid ratio) and the morphological indicators of tea buds (length, width, size, number, and their growth rates) at the late sprouting stage (<xref ref-type="fig" rid="f8"><bold>Figure&#xa0;8</bold></xref>). This clear &#x201c;morphology-chemistry&#x201d; decoupling phenomenon challenges the traditional empirical belief that intrinsic tea quality can be inferred from appearance traits, such as bud size. The underlying mechanism may stem from a fundamental hierarchical separation between the biological processes of morphogenesis and metabolite accumulation, in terms of their regulatory logic, temporal scales, and resource allocation.</p>
<p>From the perspective of molecular regulation, morphological development and chemical quality formation may be controlled by distinct genetic programs and environmental response networks. The morphogenesis of tea buds, especially key traits like size, is primarily governed by a relatively fixed developmental program. Research indicates that key genes such as CsKNOX6 play a central role in regulating tea bud cell division and expansion, and variation in their expression can significantly affect bud size, reflecting the strong genetic constraint on morphological development (<xref ref-type="bibr" rid="B50">Zhang et&#xa0;al., 2025b</xref>). More importantly, the morphological plasticity of leaves (buds) itself is a highly integrated, modular process controlled by complex gene regulatory networks (GRNs). Environmental signals (e.g., light, temperature) must influence key components within these networks, such as transcription factors, chromatin states, and hormone signaling pathways, to trigger amplified yet relatively slow morphological adjustments (<xref ref-type="bibr" rid="B17">Gonzalez et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B16">Fritz et&#xa0;al., 2018</xref>). This regulatory mode determines the &#x201c;inertia&#x201d; and cumulative effects of morphological traits in response to the environment.</p>
<p>In stark contrast, the biosynthesis and accumulation of key components determining tea flavor chemical quality&#x2014;amino acids (e.g., theanine) and tea polyphenols (primarily catechins)&#x2014;are governed by a highly dynamic and plastic metabolic network. These secondary metabolites are core components of the plant&#x2019;s chemical defense system against environmental stresses (both biotic and abiotic), and their synthesis is characterized by rapid induction and reversibility (<xref ref-type="bibr" rid="B33">Metlen et&#xa0;al., 2009</xref>; <xref ref-type="bibr" rid="B41">Scheres and van der Putten, 2017</xref>). Substantial evidence shows that a range of environmental factors, including seasonal shifts, water stress, temperature fluctuations, and light changes, can rapidly and significantly (with variation amplitudes up to 50%) alter the content of secondary metabolites in tea leaves. However, the direction and extent of these effects vary depending on the specific environment, cultivar, and management practices, and there are even contradictory reports, highlighting the extreme complexity and context-dependency of metabolic network responses (<xref ref-type="bibr" rid="B2">Ahmed et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B48">Wu et&#xa0;al., 2025</xref>). In short, morphological development addresses &#x201c;how to build a stable structure,&#x201d; while chemical quality formation addresses &#x201c;how to perform dynamic biochemical defense and adaptation.&#x201d; There is an essential difference in their regulatory objectives and response dynamics.</p>
<p>This difference in regulatory hierarchy may lead to a temporal misalignment between the two processes on the developmental timescale. Morphological development (e.g., cell division, tissue differentiation) is a continuous, resource-intensive process, and its changes, once initiated, are not easily reversed. In contrast, the synthesis and turnover of secondary metabolites may fluctuate at a faster pace in response to transient environmental signals. The observations in this study were concentrated in the specific time window of the late sprouting stage, potentially capturing a &#x201c;phase difference&#x201d; between the key stage of morphological development and the dynamic equilibrium period of chemical substances. This temporal asynchrony might be an important reason for observing phenotypic decoupling at a single time point. Other studies also hint at the complexity of morphology-metabolism relationships. For instance, in tea leaves, photosynthetic and respiratory parameters were significantly correlated with tea polyphenol content but not with total amino acid content, indicating that the coupling strength between different metabolic pathways and physiological functions inherently differs (<xref ref-type="bibr" rid="B30">Li et&#xa0;al., 2016</xref>).</p>
<p>Furthermore, internal resource allocation trade-offs within the plant may provide a further explanation for the observed stage-specific decoupling phenomenon. The lack of correlation between tea bud morphogenesis (growth) and chemical quality formation (defense) can be interpreted through the classic &#x201c; theory of plant defense&#x201d; in plant ecology. This theory posits that under limited resource conditions, plants cannot simultaneously maximize investment in both primary growth (e.g., tissue expansion) and secondary metabolic defense (e.g., synthesis of phenolic compounds) (<xref ref-type="bibr" rid="B10">Coley et&#xa0;al., 1985</xref>; <xref ref-type="bibr" rid="B20">Herms and Mattson, 1992</xref>; <xref ref-type="bibr" rid="B43">Stamp, 2003</xref>). As mentioned above, our study indicates that trade-offs already exist among different morphological components during plant growth by the late sprouting stage (e.g., between bud &#x201c;width&#x201d; and &#x201c;number&#x201d;) (<xref ref-type="fig" rid="f7"><bold>Figure&#xa0;7</bold></xref>). Therefore, during the rapid sprouting phase, the developmental program likely prioritizes the allocation of resources to core morphogenetic processes, while the allocation for synthesizing key secondary metabolites (e.g., theanine, catechins) is relatively constrained. This differential resource allocation at critical developmental stages can lead to independent variation trajectories between morphological and chemical traits&#x2014;that is, the observed stage-specific decoupling phenomenon. In summary, this suggests that morphological construction and chemical quality formation may operate as partially independent modules. Consequently, agronomic measures targeting one module may not induce proportional or linear changes in the other.</p>
<p>The stage-specific decoupling between morphology and chemistry implies that yield- and quality-related traits in tea may be regulated independently. However, this study&#x2019;s sampling of chemical quality only at a single, late developmental stage limits our ability to characterize dynamic relationships or transient coordination between these trait categories across the entire sprouting period. Therefore, future research should integrate time&#x2212;series multi&#x2212;omics (transcriptomics, metabolomics) with continuous, precise phenotyping to dynamically track the trajectories of morphological and chemical traits and to identify the key regulatory nodes underlying their independent or coupled expression.</p>
</sec>
<sec id="s4_4">
<label>4.4</label>
<title>Implications for precision cultivation practices</title>
<p>The systematic understanding described above lays an important foundation for developing quality-oriented tea garden management strategies. The results of this study can provide references for integrating multi-level strategies from macro-mode selection to field micro-management. The following preliminary ideas are elaborated to offer insights for the targeted regulation of tea quality.</p>
<p>First, in the strategic selection of cultivation modes, it should be closely aligned with the quality requirements of the target tea category. The results of this study indicate that Plain Tea-Soybean Intercropping (PTS) is the preferred mode for producing &#x201c;refreshing&#x201d; premium green tea, as it combines the highest free amino acid content (4.3%) with significantly larger tea bud width (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B</bold></xref>, <xref ref-type="fig" rid="f4"><bold>4A</bold></xref>), providing ideal raw material for shaping high-end green tea with &#x201c;good taste and beautiful appearance&#x201d;. Conversely, Hilly Tea Monoculture (HT), due to its higher polyphenol-to-amino acid ratio and maintained good bud width (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3B</bold></xref>, <xref ref-type="fig" rid="f4"><bold>4A</bold></xref>), is more suitable for developing &#x201c;mellow and robust&#x201d; black tea products. This quality potential-based mode matching is a key decision for maximizing tea garden ecosystem services and economic benefits at the source.</p>
<p>Secondly, during the dynamic process of tea bud sprouting, stage-specific precise morphological regulation should be implemented. Based on the inherent programmability and plasticity windows of tea bud development, management measures need to vary with time. In the early sprouting stage, when tea bud width shows a significant decrease (-21%) (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1B</bold></xref>), the management core should focus on &#x201c;preserving width&#x201d;, potentially through strategies like foliar application of calcium fertilizer to enhance cell wall stability and delay lateral contraction. Entering the mid-to-late sprouting stage, when the increase in tea bud length significantly accelerates and the growth rate of tea bud number decreases (<xref ref-type="fig" rid="f1"><bold>Figures&#xa0;1A, C</bold></xref>), the focus should shift to &#x201c;promoting length and stabilizing density&#x201d;, potentially by ensuring timely nitrogen supply to promote longitudinal elongation.</p>
<p>Finally, in response to the stage-specific &#x201c;morphology-chemistry&#x201d; decoupling phenomenon, a new paradigm of synergistic management must be developed. We propose establishing a &#x201c;dual-track monitoring&#x201d; system to achieve parallel optimization of yield and quality. Track one: use image analysis technology to track key morphological indicators (e.g., width/number ratio) in real-time (<xref ref-type="bibr" rid="B55">Zhang et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B24">Jiang et&#xa0;al., 2025a</xref>; <xref ref-type="bibr" rid="B44">Sun et&#xa0;al., 2025</xref>). When the absolute value of their correlation (|r|) exceeds the competition threshold of 0.5, guide corresponding &#x201c;width-preserving&#x201d; or &#x201c;density-increasing&#x201d; agronomic operations promptly. Track two, which is key to breaking the decoupling bottleneck, is to actively explore the use of multispectral or hyperspectral remote sensing technology to directly and non-destructively monitor the chemical components of young shoots (e.g., amino acid accumulation dynamics) (<xref ref-type="bibr" rid="B38">Ren et&#xa0;al., 2024</xref>; <xref ref-type="bibr" rid="B25">Jiang et&#xa0;al., 2025b</xref>). This would bypass the traditional misconception of indirectly inferring quality from morphology, ultimately enabling direct and precise regulation of tea chemical quality.</p>
</sec>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusions</title>
<p>This study demonstrates that tea bud development follows a conserved program yet exhibits distinct trait-specific plasticity. While bud length was strongly governed by the sprouting stage, width, and number were more responsive to cultivation practices. Morphological traits were significantly more malleable to cultivation mode than chemical quality. Furthermore, the correlations among morphological traits shifted dynamically during sprouting: bud length and width transitioned from positive to negative correlation, as did width and number, indicating a progressive shift in resource allocation priorities during development.</p>
<p>Notably, at the final sprouting stage measured in this study, no significant correlations were detected between morphological traits and chemical quality indicators. This observed absence of correlation is based on analysis from a single, late developmental time point, which may not capture potential transient or stage-specific relationships that could occur earlier in bud development. Nevertheless, the clear lack of association at this key harvest-relevant stage challenges the widespread traditional assumption that larger buds inherently indicate better tea quality. Further research employing continuous temporal monitoring of both morphological and chemical dynamics will be important to fully clarify their relationship across the entire developmental progression.</p>
<p>Collectively, this work advances our understanding of gene-environment interactions in perennial crops and provides an actionable framework for precision management to simultaneously optimize tea yield and quality.</p>
</sec>
</body>
<back>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p></sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>WZ: Conceptualization, Data curation, Methodology, Writing &#x2013; review &amp; editing, Funding acquisition, Visualization, Formal Analysis. JY: Visualization, Investigation, Methodology, Writing &#x2013; original draft. SX: Writing &#x2013; original draft, Investigation. WH: Writing &#x2013; review &amp; editing, Investigation, Resources. JT: Investigation, Writing &#x2013; review &amp; editing. YZ: Writing &#x2013; review &amp; editing, Investigation. ZH: Writing &#x2013; review &amp; editing, Investigation. ZY: Investigation, Writing &#x2013; review &amp; editing.</p></sec>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>Author WH was employed by the company Fuliang Changnan Tea Co., Ltd.</p>
<p>The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p></sec>
<sec id="s10" sec-type="ai-statement">
<title>Generative AI statement</title>
<p>The author(s) declared that generative AI was not used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p></sec>
<sec id="s11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p></sec>
<sec id="s12" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2026.1747538/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2026.1747538/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="DataSheet1.docx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/></sec>
<ref-list>
<title>References</title>
<ref id="B1">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Aaqil</surname> <given-names>M.</given-names></name>
<name><surname>Peng</surname> <given-names>C.</given-names></name>
<name><surname>Kamal</surname> <given-names>A.</given-names></name>
<name><surname>Nawaz</surname> <given-names>T.</given-names></name>
<name><surname>Zhang</surname> <given-names>F.</given-names></name>
<name><surname>Gong</surname> <given-names>J.</given-names></name>
</person-group>. (<year>2023</year>). 
<article-title>Tea harvesting and processing techniques and its effect on phytochemical profile and final quality of black tea: A review</article-title>. <source>Foods</source> <volume>12</volume>, <fpage>4467</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/foods12244467</pub-id>, PMID: <pub-id pub-id-type="pmid">38137271</pub-id>
</mixed-citation>
</ref>
<ref id="B2">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Ahmed</surname> <given-names>S.</given-names></name>
<name><surname>Griffin</surname> <given-names>T. S.</given-names></name>
<name><surname>Kraner</surname> <given-names>D.</given-names></name>
<name><surname>Schaffner</surname> <given-names>M. K.</given-names></name>
<name><surname>Sharma</surname> <given-names>D.</given-names></name>
<name><surname>Hazel</surname> <given-names>M.</given-names></name>
<etal/>
</person-group>. (<year>2019</year>). 
<article-title>Environmental factors variably impact tea secondary metabolites in the context of climate change</article-title>. <source>Front. Plant Sci.</source> <volume>10</volume>, <elocation-id>939</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2019.00939</pub-id>, PMID: <pub-id pub-id-type="pmid">31475018</pub-id>
</mixed-citation>
</ref>
<ref id="B3">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>An</surname> <given-names>Y.</given-names></name>
<name><surname>Chen</surname> <given-names>L.</given-names></name>
<name><surname>Tao</surname> <given-names>L.</given-names></name>
<name><surname>Liu</surname> <given-names>S.</given-names></name>
<name><surname>Wei</surname> <given-names>C.</given-names></name>
</person-group> (<year>2021</year>). 
<article-title>QTL mapping for leaf area of tea plants (Camellia sinensis) based on a high-quality genetic map constructed by whole genome resequencing</article-title>. <source>Front. Plant Sci.</source> <volume>12</volume>, <elocation-id>705285</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2021.705285</pub-id>, PMID: <pub-id pub-id-type="pmid">34394160</pub-id>
</mixed-citation>
</ref>
<ref id="B4">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name><surname>Berm&#xfa;dez</surname> <given-names>S.</given-names></name>
<name><surname>Voora</surname> <given-names>V.</given-names></name>
<name><surname>Larrea</surname> <given-names>C.</given-names></name>
<name><surname>Luna</surname> <given-names>E.</given-names></name>
</person-group> (<year>2024</year>). <source>Global Market Report: Tea prices and sustainability</source>. 
<publisher-name>International Institute for Sustainable Development</publisher-name>. Available online at: <uri xlink:href="https://www.iisd.org/system/files/2024-01/2024-global-market-report-tea.pdf">https://www.iisd.org/system/files/2024-01/2024-global-market-report-tea.pdf</uri>
</mixed-citation>
</ref>
<ref id="B5">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Brown</surname> <given-names>R.</given-names></name>
<name><surname>Rickless</surname> <given-names>P.</given-names></name>
</person-group> (<year>1949</year>). 
<article-title>A new method for the study of cell division and cell extension with some preliminary observations on the effect of temperature and of nutrients</article-title>. <source>Proc. R. Soc. London. Ser. B-Biological Sci.</source> <volume>136</volume>, <fpage>110</fpage>&#x2013;<lpage>125</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1098/rspb.1949.0008</pub-id>, PMID: <pub-id pub-id-type="pmid">18130566</pub-id>
</mixed-citation>
</ref>
<ref id="B6">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Callaway</surname> <given-names>R. M.</given-names></name>
<name><surname>Walker</surname> <given-names>L. R.</given-names></name>
</person-group> (<year>1997</year>). 
<article-title>Competition and facilitation: a synthetic approach to interactions in plant communities</article-title>. <source>Ecology</source> <volume>78</volume>, <fpage>1958</fpage>&#x2013;<lpage>1965</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/0012-9658(1997)078[1958:CAFASA]2.0.CO;2</pub-id>
</mixed-citation>
</ref>
<ref id="B7">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Carbone</surname> <given-names>F.</given-names></name>
<name><surname>Preuss</surname> <given-names>A.</given-names></name>
<name><surname>De Vos</surname> <given-names>R. C.</given-names></name>
<name><surname>D&#x2019;amico</surname> <given-names>E.</given-names></name>
<name><surname>Perrotta</surname> <given-names>G.</given-names></name>
<name><surname>Bovy</surname> <given-names>A. G.</given-names></name>
<etal/>
</person-group>. (<year>2009</year>). 
<article-title>Developmental, genetic and environmental factors affect the expression of flavonoid genes, enzymes and metabolites in strawberry fruits</article-title>. <source>Plant Cell Environ.</source> <volume>32</volume>, <fpage>1117</fpage>&#x2013;<lpage>1131</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-3040.2009.01994.x</pub-id>, PMID: <pub-id pub-id-type="pmid">19422609</pub-id>
</mixed-citation>
</ref>
<ref id="B8">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Casal</surname> <given-names>J. J.</given-names></name>
<name><surname>Balasubramanian</surname> <given-names>S.</given-names></name>
</person-group> (<year>2019</year>). 
<article-title>Thermomorphogenesis</article-title>. <source>Annu. Rev. Plant Biol.</source> <volume>70</volume>, <fpage>321</fpage>&#x2013;<lpage>346</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev-arplant-050718-095919</pub-id>, PMID: <pub-id pub-id-type="pmid">30786235</pub-id>
</mixed-citation>
</ref>
<ref id="B9">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Chapin</surname> <given-names>F. S.</given-names></name>
<name><surname>Bloom</surname> <given-names>A. J.</given-names></name>
<name><surname>Field</surname> <given-names>C. B.</given-names></name>
<name><surname>Waring</surname> <given-names>R. H.</given-names></name>
</person-group> (<year>1987</year>). 
<article-title>Plant responses to multiple environmental factors</article-title>. <source>Bioscience</source> <volume>37</volume>, <fpage>49</fpage>&#x2013;<lpage>57</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/1310177</pub-id>
</mixed-citation>
</ref>
<ref id="B10">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Coley</surname> <given-names>P. D.</given-names></name>
<name><surname>Bryant</surname> <given-names>J. P.</given-names></name>
<name><surname>Chapin</surname> <given-names>F. S.</given-names> <suffix>III</suffix></name>
</person-group> (<year>1985</year>). 
<article-title>Resource availability and plant antiherbivore defense</article-title>. <source>Science</source> <volume>230</volume>, <fpage>895</fpage>&#x2013;<lpage>899</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.230.4728.895</pub-id>, PMID: <pub-id pub-id-type="pmid">17739203</pub-id>
</mixed-citation>
</ref>
<ref id="B11">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>De Costa</surname> <given-names>W.</given-names></name>
<name><surname>Mohotti</surname> <given-names>A. J.</given-names></name>
<name><surname>Wijeratne</surname> <given-names>M. A.</given-names></name>
</person-group> (<year>2007</year>). 
<article-title>Ecophysiology of tea</article-title>. <source>Braz. J. Plant Physiol.</source> <volume>19</volume>, <fpage>299</fpage>&#x2013;<lpage>332</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1590/S1677-04202007000400005</pub-id>
</mixed-citation>
</ref>
<ref id="B12">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Dong</surname> <given-names>N. Q.</given-names></name>
<name><surname>Lin</surname> <given-names>H. X.</given-names></name>
</person-group> (<year>2021</year>). 
<article-title>Contribution of phenylpropanoid metabolism to plant development and plant&#x2013;environment interactions</article-title>. <source>J. Integr. Plant Biol.</source> <volume>63</volume>, <fpage>180</fpage>&#x2013;<lpage>209</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/jipb.13054</pub-id>, PMID: <pub-id pub-id-type="pmid">33325112</pub-id>
</mixed-citation>
</ref>
<ref id="B13">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Duan</surname> <given-names>Y.</given-names></name>
<name><surname>Shang</surname> <given-names>X.</given-names></name>
<name><surname>Liu</surname> <given-names>G.</given-names></name>
<name><surname>Zou</surname> <given-names>Z.</given-names></name>
<name><surname>Zhu</surname> <given-names>X.</given-names></name>
<name><surname>Ma</surname> <given-names>Y.</given-names></name>
<etal/>
</person-group>. (<year>2021</year>). 
<article-title>The effects of tea plants-soybean intercropping on the secondary metabolites of tea plants by metabolomics analysis</article-title>. <source>BMC Plant Biol.</source> <volume>21</volume>, <fpage>1</fpage>&#x2013;<lpage>14</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12870-021-03258-1</pub-id>, PMID: <pub-id pub-id-type="pmid">34686144</pub-id>
</mixed-citation>
</ref>
<ref id="B14">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Dutta</surname> <given-names>R.</given-names></name>
<name><surname>Stein</surname> <given-names>A.</given-names></name>
<name><surname>Smaling</surname> <given-names>E.</given-names></name>
<name><surname>Bhagat</surname> <given-names>R.</given-names></name>
<name><surname>Hazarika</surname> <given-names>M.</given-names></name>
</person-group> (<year>2010</year>). 
<article-title>Effects of plant age and environmental and management factors on tea yield in Northeast India</article-title>. <source>Agron. J.</source> <volume>102</volume>, <fpage>1290</fpage>&#x2013;<lpage>1301</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2134/agronj2010.0091</pub-id>
</mixed-citation>
</ref>
<ref id="B15">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>El-Soda</surname> <given-names>M.</given-names></name>
<name><surname>Malosetti</surname> <given-names>M.</given-names></name>
<name><surname>Zwaan</surname> <given-names>B. J.</given-names></name>
<name><surname>Koornneef</surname> <given-names>M.</given-names></name>
<name><surname>Aarts</surname> <given-names>M. G.</given-names></name>
</person-group> (<year>2014</year>). 
<article-title>Genotype&#xd7; environment interaction QTL mapping in plants: lessons from Arabidopsis</article-title>. <source>Trends Plant Sci.</source> <volume>19</volume>, <fpage>390</fpage>&#x2013;<lpage>398</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tplants.2014.01.001</pub-id>, PMID: <pub-id pub-id-type="pmid">24491827</pub-id>
</mixed-citation>
</ref>
<ref id="B16">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Fritz</surname> <given-names>M. A.</given-names></name>
<name><surname>Rosa</surname> <given-names>S.</given-names></name>
<name><surname>Sicard</surname> <given-names>A.</given-names></name>
</person-group> (<year>2018</year>). 
<article-title>Mechanisms underlying the environmentally induced plasticity of leaf morphology</article-title>. <source>Front. Genet.</source> <volume>9</volume>, <elocation-id>478</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fgene.2018.00478</pub-id>, PMID: <pub-id pub-id-type="pmid">30405690</pub-id>
</mixed-citation>
</ref>
<ref id="B17">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Gonzalez</surname> <given-names>N.</given-names></name>
<name><surname>Vanhaeren</surname> <given-names>H.</given-names></name>
<name><surname>Inz&#xe9;</surname> <given-names>D.</given-names></name>
</person-group> (<year>2012</year>). 
<article-title>Leaf size control: complex coordination of cell division and expansion</article-title>. <source>Trends Plant Sci.</source> <volume>17</volume>, <fpage>332</fpage>&#x2013;<lpage>340</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tplants.2012.02.003</pub-id>, PMID: <pub-id pub-id-type="pmid">22401845</pub-id>
</mixed-citation>
</ref>
<ref id="B18">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Good</surname> <given-names>A. G.</given-names></name>
<name><surname>Zaplachinski</surname> <given-names>S. T.</given-names></name>
</person-group> (<year>1994</year>). 
<article-title>The effects of drought stress on free amino acid accumulation and protein synthesis in Brassica napus</article-title>. <source>Physiologia plantarum</source> <volume>90</volume>, <fpage>9</fpage>&#x2013;<lpage>14</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1399-3054.1994.tb02185.x</pub-id>
</mixed-citation>
</ref>
<ref id="B19">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Han</surname> <given-names>W.-Y.</given-names></name>
<name><surname>Huang</surname> <given-names>J.-G.</given-names></name>
<name><surname>Li</surname> <given-names>X.</given-names></name>
<name><surname>Li</surname> <given-names>Z.-X.</given-names></name>
<name><surname>Ahammed</surname> <given-names>G. J.</given-names></name>
<name><surname>Yan</surname> <given-names>P.</given-names></name>
<etal/>
</person-group>. (<year>2017</year>). 
<article-title>Altitudinal effects on the quality of green tea in east China: a climate change perspective</article-title>. <source>Eur. Food Res. Technol.</source> <volume>243</volume>, <fpage>323</fpage>&#x2013;<lpage>330</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00217-016-2746-5</pub-id>
</mixed-citation>
</ref>
<ref id="B20">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Herms</surname> <given-names>D. A.</given-names></name>
<name><surname>Mattson</surname> <given-names>W. J.</given-names></name>
</person-group> (<year>1992</year>). 
<article-title>The dilemma of plants: to grow or defend</article-title>. <source>Q. Rev. Biol.</source> <volume>67</volume>, <fpage>283</fpage>&#x2013;<lpage>335</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1086/417659</pub-id>
</mixed-citation>
</ref>
<ref id="B21">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Hua</surname> <given-names>S.</given-names></name>
<name><surname>Dal-Bianco</surname> <given-names>M.</given-names></name>
<name><surname>Chen</surname> <given-names>Z.-H.</given-names></name>
</person-group> (<year>2022</year>). 
<article-title>Crop yield and quality response to the interaction between environment and genetic factors</article-title>. <source>Front. Media SA</source>, <elocation-id>823279</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fgene.2022.823279</pub-id>, PMID: <pub-id pub-id-type="pmid">35368712</pub-id>
</mixed-citation>
</ref>
<ref id="B22">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Iwasa</surname> <given-names>Y.</given-names></name>
</person-group> (<year>2000</year>). 
<article-title>Dynamic optimization of plant growth</article-title>. <source>Evolutionary Ecol. Res.</source> <volume>2</volume>, <fpage>437</fpage>&#x2013;<lpage>455</lpage>.
</mixed-citation>
</ref>
<ref id="B23">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Jeyaramraja</surname> <given-names>P.</given-names></name>
<name><surname>Pius</surname> <given-names>P.</given-names></name>
<name><surname>Raj Kumar</surname> <given-names>R.</given-names></name>
<name><surname>Jayakumar</surname> <given-names>D.</given-names></name>
</person-group> (<year>2003</year>). 
<article-title>Soil moisture stress-induced alterations in bioconstituents determining tea quality</article-title>. <source>J. Sci. Food Agric.</source> <volume>83</volume>, <fpage>1187</fpage>&#x2013;<lpage>1191</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/jsfa.1440</pub-id>
</mixed-citation>
</ref>
<ref id="B24">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Jiang</surname> <given-names>J.</given-names></name>
<name><surname>Ji</surname> <given-names>H.</given-names></name>
<name><surname>Zhou</surname> <given-names>G.</given-names></name>
<name><surname>Pan</surname> <given-names>R.</given-names></name>
<name><surname>Zhao</surname> <given-names>L.</given-names></name>
<name><surname>Duan</surname> <given-names>Z.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>a). 
<article-title>Non-destructive monitoring of tea plant growth through UAV spectral imagery and meteorological data using machine learning and parameter optimization algorithms</article-title>. <source>Comput. Electron. Agric.</source> <volume>229</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.compag.2024.109795</pub-id>
</mixed-citation>
</ref>
<ref id="B25">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Jiang</surname> <given-names>J.</given-names></name>
<name><surname>Zhou</surname> <given-names>G.</given-names></name>
<name><surname>Ji</surname> <given-names>H.</given-names></name>
<name><surname>Pan</surname> <given-names>R.</given-names></name>
<name><surname>Yang</surname> <given-names>L.</given-names></name>
<name><surname>Zhao</surname> <given-names>S.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>b). 
<article-title>Advanced hyper-spectral feature engineering for real-time prediction of tea quality in support of high-grade cultivation</article-title>. <source>Comput. Electron. Agric.</source> <volume>236</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.compag.2025.110485</pub-id>
</mixed-citation>
</ref>
<ref id="B26">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Kc</surname> <given-names>S.</given-names></name>
<name><surname>Long</surname> <given-names>L.</given-names></name>
<name><surname>Liu</surname> <given-names>M.</given-names></name>
<name><surname>Zhang</surname> <given-names>Q.</given-names></name>
<name><surname>Ruan</surname> <given-names>J.</given-names></name>
</person-group> (<year>2021</year>). 
<article-title>Light intensity modulates the effect of phosphate limitation on carbohydrates, amino acids, and catechins in tea plants (Camellia sinensis L.)</article-title>. <source>Front. Plant Sci.</source> <volume>12</volume>, <elocation-id>743781</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2021.743781</pub-id>, PMID: <pub-id pub-id-type="pmid">34691121</pub-id>
</mixed-citation>
</ref>
<ref id="B27">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Leduc</surname> <given-names>N.</given-names></name>
<name><surname>Roman</surname> <given-names>H.</given-names></name>
<name><surname>Barbier</surname> <given-names>F.</given-names></name>
<name><surname>P&#xe9;ron</surname> <given-names>T.</given-names></name>
<name><surname>Huch&#xe9;-Th&#xe9;lier</surname> <given-names>L.</given-names></name>
<name><surname>Lothier</surname> <given-names>J.</given-names></name>
<etal/>
</person-group>. (<year>2014</year>). 
<article-title>Light signaling in bud outgrowth and branching in plants</article-title>. <source>Plants</source> <volume>3</volume>, <fpage>223</fpage>&#x2013;<lpage>250</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/plants3020223</pub-id>, PMID: <pub-id pub-id-type="pmid">27135502</pub-id>
</mixed-citation>
</ref>
<ref id="B28">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Li</surname> <given-names>J.</given-names></name>
<name><surname>Wang</surname> <given-names>J.</given-names></name>
<name><surname>Yao</surname> <given-names>Y.</given-names></name>
<name><surname>Hua</surname> <given-names>J.</given-names></name>
<name><surname>Zhou</surname> <given-names>Q.</given-names></name>
<name><surname>Jiang</surname> <given-names>Y.</given-names></name>
<etal/>
</person-group>. (<year>2020</year>a). 
<article-title>Phytochemical comparison of different tea (Camellia sinensis) cultivars and its association with sensory quality of finished tea</article-title>. <source>Lwt</source> <volume>117</volume>, <fpage>108595</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.lwt.2019.108595</pub-id>
</mixed-citation>
</ref>
<ref id="B29">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Li</surname> <given-names>Y.</given-names></name>
<name><surname>Kong</surname> <given-names>D.</given-names></name>
<name><surname>Fu</surname> <given-names>Y.</given-names></name>
<name><surname>Sussman</surname> <given-names>M. R.</given-names></name>
<name><surname>Wu</surname> <given-names>H.</given-names></name>
</person-group> (<year>2020</year>b). 
<article-title>The effect of developmental and environmental factors on secondary metabolites in medicinal plants</article-title>. <source>Plant Physiol. Biochem.</source> <volume>148</volume>, <fpage>80</fpage>&#x2013;<lpage>89</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.plaphy.2020.01.006</pub-id>, PMID: <pub-id pub-id-type="pmid">31951944</pub-id>
</mixed-citation>
</ref>
<ref id="B30">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Li</surname> <given-names>Z.-X.</given-names></name>
<name><surname>Yang</surname> <given-names>W.-J.</given-names></name>
<name><surname>Ahammed</surname> <given-names>G. J.</given-names></name>
<name><surname>Shen</surname> <given-names>C.</given-names></name>
<name><surname>Yan</surname> <given-names>P.</given-names></name>
<name><surname>Li</surname> <given-names>X.</given-names></name>
<etal/>
</person-group>. (<year>2016</year>). 
<article-title>Developmental changes in carbon and nitrogen metabolism affect tea quality in different leaf position</article-title>. <source>Plant Physiol. Biochem.</source> <volume>106</volume>, <fpage>327</fpage>&#x2013;<lpage>335</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.plaphy.2016.06.027</pub-id>, PMID: <pub-id pub-id-type="pmid">27380366</pub-id>
</mixed-citation>
</ref>
<ref id="B31">
<mixed-citation publication-type="book">
<person-group person-group-type="author">
<name><surname>Lynch</surname> <given-names>M.</given-names></name>
<name><surname>Walsh</surname> <given-names>B.</given-names></name>
</person-group> (<year>1998</year>). <source>Genetics and analysis of quantitative traits</source> (
<publisher-name>Sinauer Sunderland</publisher-name>). <publisher-loc>MA: Sinauer</publisher-loc>. <volume>1</volume>:<fpage>535</fpage>&#x2013;<lpage>557</lpage>.
</mixed-citation>
</ref>
<ref id="B32">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Malosetti</surname> <given-names>M.</given-names></name>
<name><surname>Ribaut</surname> <given-names>J.-M.</given-names></name>
<name><surname>van Eeuwijk</surname> <given-names>F. A.</given-names></name>
</person-group> (<year>2013</year>). 
<article-title>The statistical analysis of multi-environment data: modeling genotype-by-environment interaction and its genetic basis</article-title>. <source>Front. Physiol.</source> <volume>4</volume>, <elocation-id>44</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fphys.2013.00044</pub-id>, PMID: <pub-id pub-id-type="pmid">23487515</pub-id>
</mixed-citation>
</ref>
<ref id="B33">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Metlen</surname> <given-names>K. L.</given-names></name>
<name><surname>Aschehoug</surname> <given-names>E. T.</given-names></name>
<name><surname>Callaway</surname> <given-names>R. M.</given-names></name>
</person-group> (<year>2009</year>). 
<article-title>Plant behavioural ecology: dynamic plasticity in secondary metabolites</article-title>. <source>Plant Cell Environ.</source> <volume>32</volume>, <fpage>641</fpage>&#x2013;<lpage>653</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-3040.2008.01910.x</pub-id>, PMID: <pub-id pub-id-type="pmid">19021888</pub-id>
</mixed-citation>
</ref>
<ref id="B34">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Pokharel</surname> <given-names>S. S.</given-names></name>
<name><surname>Yu</surname> <given-names>H.</given-names></name>
<name><surname>Fang</surname> <given-names>W.</given-names></name>
<name><surname>Parajulee</surname> <given-names>M. N.</given-names></name>
<name><surname>Chen</surname> <given-names>F.</given-names></name>
</person-group> (<year>2023</year>). 
<article-title>Intercropping cover crops for a vital ecosystem service: A review of the biocontrol of insect pests in tea agroecosystems</article-title>. <source>Plants</source> <volume>12</volume>, <fpage>2361</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/plants12122361</pub-id>, PMID: <pub-id pub-id-type="pmid">37375986</pub-id>
</mixed-citation>
</ref>
<ref id="B35">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Qifang</surname> <given-names>S.</given-names></name>
<name><surname>Luhan</surname> <given-names>Z.</given-names></name>
<name><surname>Wenjie</surname> <given-names>Y.</given-names></name>
<name><surname>Fang</surname> <given-names>D.</given-names></name>
<name><surname>Yan</surname> <given-names>L.</given-names></name>
<name><surname>Xiaoqing</surname> <given-names>F.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>). 
<article-title>Effects of meteorological factors on content of free amino acids and caffeine in green tea</article-title>. <source>J. Appl. Meteorological Sci.</source> <volume>36</volume>, <fpage>307</fpage>&#x2013;<lpage>315</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.11898/1001-7313.20250305</pub-id>
</mixed-citation>
</ref>
<ref id="B36">
<mixed-citation publication-type="book">
<person-group person-group-type="author"><collab>R Core Team</collab>
</person-group>. (<year>2021</year>). <source>R: A language and environment for statistical computing</source>. 
<publisher-name>R Foundation for Statistical Computing</publisher-name>. Available online at: <uri xlink:href="https://www.r-project.org/">https://www.r-project.org/</uri>
</mixed-citation>
</ref>
<ref id="B37">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Ren</surname> <given-names>X.</given-names></name>
<name><surname>Lin</surname> <given-names>M.</given-names></name>
<name><surname>Liu</surname> <given-names>J.</given-names></name>
<name><surname>Khan</surname> <given-names>W.</given-names></name>
<name><surname>Zhao</surname> <given-names>H.</given-names></name>
<name><surname>Sun</surname> <given-names>B.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>). 
<article-title>Effects of altitude on tea composition: dual regulation by soil physicochemical properties and microbial communities</article-title>. <source>Plants</source> <volume>14</volume>, <fpage>1642</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/plants14111642</pub-id>, PMID: <pub-id pub-id-type="pmid">40508317</pub-id>
</mixed-citation>
</ref>
<ref id="B38">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Ren</surname> <given-names>G.</given-names></name>
<name><surname>Wu</surname> <given-names>R.</given-names></name>
<name><surname>Yin</surname> <given-names>L.</given-names></name>
<name><surname>Zhang</surname> <given-names>Z.</given-names></name>
<name><surname>Ning</surname> <given-names>J.</given-names></name>
</person-group> (<year>2024</year>). 
<article-title>Description of tea quality using deep learning and multi-sensor feature fusion</article-title>. <source>J. Food Composition Anal.</source> <volume>126</volume>, <fpage>105924</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jfca.2023.105924</pub-id>
</mixed-citation>
</ref>
<ref id="B39">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Ruan</surname> <given-names>J.</given-names></name>
<name><surname>Gerend&#xe1;s</surname> <given-names>J.</given-names></name>
<name><surname>H&#xe4;rdter</surname> <given-names>R.</given-names></name>
<name><surname>Sattelmacher</surname> <given-names>B.</given-names></name>
</person-group> (<year>2007</year>). 
<article-title>Effect of nitrogen form and root-zone pH on growth and nitrogen uptake of tea (Camellia sinensis) plants</article-title>. <source>Ann. Bot.</source> <volume>99</volume>, <fpage>301</fpage>&#x2013;<lpage>310</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/aob/mcl258</pub-id>, PMID: <pub-id pub-id-type="pmid">17204540</pub-id>
</mixed-citation>
</ref>
<ref id="B40">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Samuoliene</surname> <given-names>G.</given-names></name>
<name><surname>Virsile</surname> <given-names>A.</given-names></name>
<name><surname>Brazaityte</surname> <given-names>A.</given-names></name>
<name><surname>Jankauskiene</surname> <given-names>J.</given-names></name>
<name><surname>Sakalauskiene</surname> <given-names>S.</given-names></name>
<name><surname>Vastakaite</surname> <given-names>V.</given-names></name>
<etal/>
</person-group>. (<year>2017</year>). 
<article-title>Blue light dosage affects carotenoids and tocopherols in microgreens</article-title>. <source>Food Chem.</source> <volume>228</volume>, <fpage>50</fpage>&#x2013;<lpage>56</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.foodchem.2017.01.144</pub-id>, PMID: <pub-id pub-id-type="pmid">28317756</pub-id>
</mixed-citation>
</ref>
<ref id="B41">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Scheres</surname> <given-names>B.</given-names></name>
<name><surname>van der Putten</surname> <given-names>W. H.</given-names></name>
</person-group> (<year>2017</year>). 
<article-title>The plant perceptron connects environment to development</article-title>. <source>Nature</source> <volume>543</volume>, <fpage>337</fpage>&#x2013;<lpage>345</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nature22010</pub-id>, PMID: <pub-id pub-id-type="pmid">28300110</pub-id>
</mixed-citation>
</ref>
<ref id="B42">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Smith</surname> <given-names>R.</given-names></name>
<name><surname>Harvey</surname> <given-names>F.</given-names></name>
<name><surname>Cannell</surname> <given-names>M.</given-names></name>
</person-group> (<year>1990</year>). 
<article-title>Pattern of tea shoot growth</article-title>. <source>Exp. Agric.</source> <volume>26</volume>, <fpage>197</fpage>&#x2013;<lpage>208</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1017/S0014479700018196</pub-id>
</mixed-citation>
</ref>
<ref id="B43">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Stamp</surname> <given-names>N.</given-names></name>
</person-group> (<year>2003</year>). 
<article-title>Out of the quagmire of plant defense hypotheses</article-title>. <source>Q. Rev. Biol.</source> <volume>78</volume>, <fpage>23</fpage>&#x2013;<lpage>55</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1086/367580</pub-id>, PMID: <pub-id pub-id-type="pmid">12661508</pub-id>
</mixed-citation>
</ref>
<ref id="B44">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Sun</surname> <given-names>L.</given-names></name>
<name><surname>Shen</surname> <given-names>J.</given-names></name>
<name><surname>Mao</surname> <given-names>Y.</given-names></name>
<name><surname>Li</surname> <given-names>X.</given-names></name>
<name><surname>Fan</surname> <given-names>K.</given-names></name>
<name><surname>Qian</surname> <given-names>W.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>). 
<article-title>Discrimination of tea varieties and bud sprouting phenology using UAV-based RGB and multispectral images</article-title>. <source>Int. J. Remote Sens.</source> <volume>46</volume>, <fpage>6214</fpage>&#x2013;<lpage>6234</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/01431161.2025.2534992</pub-id>
</mixed-citation>
</ref>
<ref id="B45">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Tsukaya</surname> <given-names>H.</given-names></name>
</person-group> (<year>2005</year>). 
<article-title>Leaf shape: genetic controls and environmental factors</article-title>. <source>Int. J. Dev. Biol.</source> <volume>49</volume>, <fpage>547</fpage>&#x2013;<lpage>555</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1387/ijdb.041921ht</pub-id>, PMID: <pub-id pub-id-type="pmid">16096964</pub-id>
</mixed-citation>
</ref>
<ref id="B46">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Wachira</surname> <given-names>F.</given-names></name>
<name><surname>Ng&#x2019;etich</surname> <given-names>W.</given-names></name>
<name><surname>Omolo</surname> <given-names>J.</given-names></name>
<name><surname>Mamati</surname> <given-names>G.</given-names></name>
</person-group> (<year>2002</year>). 
<article-title>Genotype&#xd7; environment interactions for tea yields</article-title>. <source>Euphytica</source> <volume>127</volume>, <fpage>289</fpage>&#x2013;<lpage>297</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1023/A:1020273616349</pub-id>
</mixed-citation>
</ref>
<ref id="B47">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Wang</surname> <given-names>T.</given-names></name>
<name><surname>Duan</surname> <given-names>Y.</given-names></name>
<name><surname>Liu</surname> <given-names>G.</given-names></name>
<name><surname>Shang</surname> <given-names>X.</given-names></name>
<name><surname>Liu</surname> <given-names>L.</given-names></name>
<name><surname>Zhang</surname> <given-names>K.</given-names></name>
<etal/>
</person-group>. (<year>2022</year>). 
<article-title>Tea plantation intercropping green manure enhances soil functional microbial abundance and multifunctionality resistance to drying-rewetting cycles</article-title>. <source>Sci. Total Environ.</source> <volume>810</volume>, <fpage>151282</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.scitotenv.2021.151282</pub-id>, PMID: <pub-id pub-id-type="pmid">34757096</pub-id>
</mixed-citation>
</ref>
<ref id="B48">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Wu</surname> <given-names>T.</given-names></name>
<name><surname>He</surname> <given-names>J.</given-names></name>
<name><surname>Deng</surname> <given-names>X.</given-names></name>
<name><surname>Wang</surname> <given-names>X.</given-names></name>
<name><surname>Yuan</surname> <given-names>W.</given-names></name>
<name><surname>Wang</surname> <given-names>Q.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>). 
<article-title>Influencing factors and regulatory mechanisms of fresh tea leaf quality: A review</article-title>. <source>Foods</source> <volume>14</volume>, <fpage>3268</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/foods14183268</pub-id>, PMID: <pub-id pub-id-type="pmid">41008240</pub-id>
</mixed-citation>
</ref>
<ref id="B49">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Xu</surname> <given-names>C.</given-names></name>
<name><surname>Liang</surname> <given-names>L.</given-names></name>
<name><surname>Li</surname> <given-names>Y.</given-names></name>
<name><surname>Yang</surname> <given-names>T.</given-names></name>
<name><surname>Fan</surname> <given-names>Y.</given-names></name>
<name><surname>Mao</surname> <given-names>X.</given-names></name>
<etal/>
</person-group>. (<year>2021</year>). 
<article-title>Studies of quality development and major chemical composition of green tea processed from tea with different shoot maturity</article-title>. <source>Lwt</source> <volume>142</volume>, <fpage>111055</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.lwt.2021.111055</pub-id>
</mixed-citation>
</ref>
<ref id="B50">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Zhang</surname> <given-names>S.</given-names></name>
<name><surname>Chen</surname> <given-names>S.</given-names></name>
<name><surname>Fu</surname> <given-names>Z.</given-names></name>
<name><surname>Li</surname> <given-names>F.</given-names></name>
<name><surname>Chen</surname> <given-names>Q.</given-names></name>
<name><surname>Ma</surname> <given-names>J.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>b). 
<article-title>Integration of digital phenotyping, GWAS, and transcriptomic analysis revealed a key gene for bud size in tea plant (Camellia sinensis)</article-title>. <source>Horticulture Res.</source> <volume>12</volume>, <fpage>uhaf051</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/hr/uhaf051</pub-id>, PMID: <pub-id pub-id-type="pmid">40271457</pub-id>
</mixed-citation>
</ref>
<ref id="B51">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Zhang</surname> <given-names>H.</given-names></name>
<name><surname>Li</surname> <given-names>H.</given-names></name>
<name><surname>Lai</surname> <given-names>B.</given-names></name>
<name><surname>Xia</surname> <given-names>H.</given-names></name>
<name><surname>Wang</surname> <given-names>H.</given-names></name>
<name><surname>Huang</surname> <given-names>X.</given-names></name>
</person-group> (<year>2016</year>). 
<article-title>Morphological characterization and gene expression profiling during bud development in a tropical perennial, Litchi chinensis Sonn</article-title>. <source>Front. Plant Sci.</source> <volume>7</volume>, <elocation-id>1517</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2016.01517</pub-id>, PMID: <pub-id pub-id-type="pmid">27833615</pub-id>
</mixed-citation>
</ref>
<ref id="B52">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Zhang</surname> <given-names>J.</given-names></name>
<name><surname>Sun</surname> <given-names>K.</given-names></name>
<name><surname>Wang</surname> <given-names>Y.</given-names></name>
<name><surname>Qian</surname> <given-names>W.</given-names></name>
<name><surname>Sun</surname> <given-names>L.</given-names></name>
<name><surname>Shen</surname> <given-names>J.</given-names></name>
<etal/>
</person-group>. (<year>2024</year>a). 
<article-title>Integrated metabolomic and transcriptomic analyses reveal the molecular mechanism of amino acid transport between source and sink during tea shoot development</article-title>. <source>Plant Cell Rep.</source> <volume>43</volume>, <fpage>28</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00299-023-03110-w</pub-id>, PMID: <pub-id pub-id-type="pmid">38177567</pub-id>
</mixed-citation>
</ref>
<ref id="B53">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Zhang</surname> <given-names>D.</given-names></name>
<name><surname>Wei</surname> <given-names>X.</given-names></name>
<name><surname>Zhang</surname> <given-names>J.</given-names></name>
<name><surname>Cui</surname> <given-names>D.</given-names></name>
<name><surname>Zhang</surname> <given-names>P.</given-names></name>
<name><surname>Chen</surname> <given-names>S.</given-names></name>
<etal/>
</person-group>. (<year>2025</year>a). 
<article-title>Variation analysis and quantitative trait loci mapping of 16 free amino acid traits in the tea plant (Camellia sinensis)</article-title>. <source>BMC Plant Biol.</source> <volume>25</volume>, <fpage>194</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12870-024-06038-9</pub-id>, PMID: <pub-id pub-id-type="pmid">39953411</pub-id>
</mixed-citation>
</ref>
<ref id="B54">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Zhang</surname> <given-names>J.</given-names></name>
<name><surname>Yang</surname> <given-names>X.</given-names></name>
<name><surname>Zhang</surname> <given-names>X.</given-names></name>
<name><surname>Zhang</surname> <given-names>L.</given-names></name>
<name><surname>Zhang</surname> <given-names>Z.</given-names></name>
<name><surname>Zhang</surname> <given-names>Y.</given-names></name>
<etal/>
</person-group>. (<year>2022</year>). 
<article-title>Linking environmental signals to plant metabolism: The combination of field trials and environment simulators</article-title>. <source>Mol. Plant</source> <volume>15</volume>, <fpage>213</fpage>&#x2013;<lpage>215</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2021.12.017</pub-id>, PMID: <pub-id pub-id-type="pmid">34973455</pub-id>
</mixed-citation>
</ref>
<ref id="B55">
<mixed-citation publication-type="journal">
<person-group person-group-type="author">
<name><surname>Zhang</surname> <given-names>L.</given-names></name>
<name><surname>Zhang</surname> <given-names>H.</given-names></name>
<name><surname>Chen</surname> <given-names>Y.</given-names></name>
<name><surname>Dai</surname> <given-names>S.</given-names></name>
<name><surname>Li</surname> <given-names>X.</given-names></name>
<name><surname>Kenji</surname> <given-names>I.</given-names></name>
<etal/>
</person-group>. (<year>2019</year>). 
<article-title>Real-time monitoring of optimum timing for harvesting fresh tea leaves based on machine vision</article-title>. <source>Int. J. Agric. Biol. Eng.</source> <volume>12</volume>, <fpage>6</fpage>&#x2013;<lpage>9</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.25165/j.ijabe.20191201.3418</pub-id>
</mixed-citation>
</ref>
</ref-list>
<fn-group>
<fn id="n1" fn-type="custom" custom-type="edited-by">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2026725">Yuxiang Zeng</ext-link>, Chinese Academy of Agricultural Sciences, China</p></fn>
<fn id="n2" fn-type="custom" custom-type="reviewed-by">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1496047">Huan Si</ext-link>, Chinese Academy of Agricultural Sciences, China</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1669735">Ning Hao</ext-link>, The University of Tokyo, Japan</p></fn>
</fn-group>
</back>
</article>