The worldwide allometric relationship in anatomical structures for plant roots

Yue Zhang,Jing-Jing Cao,Qing-Pei Yang,Ming-Zuo Wu,Yong Zhao,De-Liang Kong

College of Forestry,Henan Agricultural University,Zhengzhou 450002,China

Keywords: Allometric relationship Root dimeter Cortex thickness Stele radius

ABSTRACT The cortex(i.e.,absorptive tissue)and stele(transportive vascular tissue)are fundamental to the function of plant roots.Unraveling how these anatomical structures are assembled in absorptive roots is essential for our understanding of plant ecology,physiology,and plant responses to global environmental changes.In this review,we first compile a large data set on anatomical traits in absorptive roots,including cortex thickness and stele radius,across 698 observations and 512 species.Using this data set,we reveal a common root allometry in absorptive root structures,i.e.,cortex thickness increases much faster than stele radius with increasing root diameter(hereafter,root allometry).Root allometry is further validated within and across plant growth forms (woody,grass,and liana species),mycorrhiza types (arbuscular mycorrhiza,ectomycorrhiza,and orchid mycorrhizas),phylogenetic gradients (from ferns to Orchidaceae),and environmental change scenarios (e.g.,elevation of atmospheric CO2 concentration and nitrogen fertilization).These findings indicate that root allometry is common in plants.Importantly,root allometry varies greatly across species.We then summarize recent research on the mechanisms of root allometry and potential issues regarding these mechanisms.We further discuss ecological and evolutionary implications of root allometry.Finally,we propose several important research directions that should be pursued regarding root allometry.

1.Introduction

Plant roots play a crucial role in plant growth,vegetation dynamics,and ecosystem functioning (e.g.,net primary productivity,nutrient cycling) as well as plant responses to environmental changes(Carmona et al.,2021;Laughlin et al.,2021;Chandregowda et al.,2022;Encinas-Valero et al.,2022;Hong et al.,2022).Although extensive research has examined both the form and function of stems and leaves(Wright et al.,2004;Diaz et al.,2016;Joswig et al.,2022),our understanding of plant roots is greatly limited.The core function of plant roots is to acquire and transport water and nutrients,which is undertaken by only a few terminal root branch orders with primary root structures,i.e.,absorptive roots (Guo et al.,2008b).

Root absorption is facilitated by microbial symbiosis (e.g.,with mycorrhizal fungi) and a range of root traits,including root morphology(e.g.,root diameter),physiology(e.g.,root respiration),anatomy(e.g.,root cortex and stele),chemistry (e.g.,root nitrogen content),mechanics(e.g.,tensile strength)and microbial symbiosis(e.g.,with mycorrhizal fungi) (McCormack et al.,2017;Wambsganss et al.,2021;Wen et al.,2022;Yan et al.,2022;Betekhtina et al.,2023;Weemstra et al.,2023).Among these root functional traits,root diameter is considered the most important,given that it usually has the largest inter-specific variation (Kong et al.,2014;Valverde-Barrantes et al.,2017;Weigelt et al.,2021).Such great variation in root diameter controls a range of root and mycorrhizal traits(Eissenstat et al.,2015;Li et al.,2018;Bergmann et al.,2020;Wen et al.,2022).Importantly,absorptive root diameter is the most phylogenetically conserved trait,suggesting that the inter-specific variation of root diameter is closely related to plant evolution(Comas et al.,2012;Chen et al.,2013;Pineiro et al.,2020;Lugli et al.,2021).

Roots are mainly composed of two cylindrical components,i.e.,cortex and stele.This anatomy determines the absorptive diameter of roots.The cortex absorbs water and nutrients both directly and indirectly,through symbiotic mycorrhizal fungi (Brundrett,2002;Ma et al.,2018;Rich et al.,2021).Root steles are responsible for the transportation of water and nutrients to stems and leaves,and of photosynthates to meet energy demands during nutrient absorption.Recent research has discovered that root growth is allometric,i.e.,the thickness of the cortex increases linearly and at a higher rate than that of the stele (hereafter,“root allometry”) (Gu et al.,2014;Kong et al.,2014),root allometry has since been reported in over 200 plant species (Kong et al.,2019).

Thick root cortex is physiologically disadvantageous for resource absorption(water and nutrients)because it increases the resistance resources encounter (Steudle and Carol,1998).In contrast,a thick stele,which bears thick vascular conduits(e.g.,vessels),usually has a high capacity to transport resources (Kong et al.,2014).Consequently,root allometry may increase the imbalance between resource absorption and transportation.Thick absorptive roots usually have low specific root length (root length per unit root biomass) and hence limited root surface area for resource absorption(Pregitzer et al.,2002).Root allometry,thus,does not enhance resource absorption by roots themselves;instead,the fast increase in nutrient absorption required to eventually balance nutrient absorption and transportation must be mediated by other agents,e.g.,mycorrhizal fungi (Kong et al.,2017).Thicker steles usually have thicker vascular conduits (e.g.,sieve tubes),which improve transportation of leaf photosynthates belowground(Jensen et al.,2016).This may,in turn,increase allometric root growth through the construction and maintenance of even thicker cortex (see Kong et al.,2021 for details).Thus,structural allometry in roots (cortex and stele) may reflect the functional coordination of the uptake,transportation and utilization of elements essential for plant life(water,nutrients and carbon).If so,root allometry paves a new way of understanding how root form and function are related to plant evolution and adaptation to environmental changes (McCormack et al.,2017;Kong et al.,2019;Bergmann et al.,2020;Zhou et al.,2022).

Although several studies have examined allometric growth in roots,various plant groups have been neglected,including Orchidaceae and vine plants.Furthermore,little is known about whether and how root allometry varies among different mycorrhizal types,e.g.,arbuscular mycorrhiza (AM),ectomycorrhiza(EM),and Orchid mycorrhiza (OM) (Brundrett,2002;Martin et al.,2017).AM fungi and EM fungi are known to differ structurally and functionally in ways that may be affected by root allometry.For instance,AM fungi penetrate root cortical cells and form arbuscules to facilitate nutrient transfer,whereas EM fungi create a sheathing mantle on the exterior of the root epidermis.Thus,increased root cortical space would provide a greater habitat for AM fungi,whereas EM fungi habitat is proportional to the root epidermal surface area(Martin et al.,2008;Terrer et al.,2016,2017,2018).It is important to learn whether the structural and functional differences between AM and EM affect their host root structures,i.e.,the allometric assembly of root cortex and stele.

Absorptive root diameter shows a strong phylogenetic signal(Kong et al.,2014;Valverde-barrantes et al.,2017;Wang et al.,2018),indicating that root anatomical traits are more influenced by ancestry than by environmental conditions.This is important for our understanding of how roots and even entire plants respond to climate change.However,we lack empirical evidence that indicates whether and/or how environmental changes (e.g.,nitrogen deposition,elevation of atmospheric CO2concentration) influence root allometry.It is necessary to test the generality of root allometry using a large species pool that includes neglected plant groups and that considers environmental changes.

Here,we aim to (1) test the generality of root allometry in different plant growth forms,mycorrhizal types,and in response to different environmental change scenarios;(2) summarize current knowledge on how root allometry is formed and on the ecological and evolutionary implications of the root allometry;and (3) propose promising directions that should be pursued in future studies on root allometry.

2.Methods

We screened data on cortex thickness and stele radius for absorptive roots in Web of Science,Google Scholar,FRED 3.0(http://roots.ornl.gov/) and CNKI (China's national knowledge infrastructure,https://www.cnki.net/) using the following keywords: “cortex”,“stele”,“anatomic structure”,“allometric relationship” or “root diameter”.Our search yielded 498,754 papers and reports.We then refined these results according to additional criteria: (1) Studies included in our analysis must be empirical rather than a review or perspective;(2)Data on root diameter and stele radius must be accessible.We also included some unpublished data (Supplementary data 2) on root anatomical traits that were measured at the same sites and following the same sampling procedures as in our previous study (Kong et al.,2014).

In total,our data set included 32 empirical studies (Supplementary data 1) with 698 observations of 512 species at 41 sites globally(Fig.1).Specifically,there were 272 woody species and 240 non-woody species(78 grass,92 herb,37 fern,28 Orchidaceae and 5 non-woody liana species).In addition,13 liana species were included in the data set.For the same species sampled in different studies,we used the mean value of the root traits across studies as the trait value of this species.For some studies with only data on stele radius and root diameter,we used the difference between root diameter and stele radius,namely the stele tissue and the tissue outside the stele (ToS),including the epidermis,exodermis and cortex.The thickness of the cortex,i.e.,tToS,has been used as a proxy approximate to the cortex thickness (Kong et al.,2019).For studies with root trait data displayed in figures,or only photos of root anatomical structures,we digitalized root trait data using the software “SigmaScan Pro software (v.5.0,SPSS Inc.,Chicago,USA)”and the software “IMAGE J (NIH Image,Bethesda,MD,USA)”,respectively.All 204 plant species in Kong’s study were included in our study.

Fig.1.A global map of 32 studies reporting data on the cortex thickness and stele radius in absorptive roots.Each study is represented by a red circle,and the size of the circle is proportional to the number of species in the study.

The existence of root allometry was tested by the following steps.First,we calculated the slope of regression of the stele radius with root diameter(Slope(RD-SR))and the slope of the regression of the cortex thickness with root diameter(Slope(RD-CT)).Then,a slope comparison was conducted between Slope (RD-SR) and Slope (RD-CT)using the method of standardized major axis(Falster et al.,2006).A significant difference in the slope comparison indicates the existence of the root allometry.Finally,we tested whether root allometry exists across the total observations and total species,among different plant growth forms (woody vs.non-woody),mycorrhizal types (AM,EM and OM),and environmental change scenarios (e.g.,nitrogen deposition,elevation of atmospheric CO2concentration).

3.Generality of root allometry in absorptive roots

Root allometry was observed in the root anatomical structures of 512 species(Table 1;Fig.2a and b).Root allometry was observed in 22 of 26 studies in which root anatomical structures were examined in at least three species(Table 2;Fig.S1;Supplementary data 1).

Table 1 The relationships between root dimeter (RD) and cortex thickness (CT) and stele radius (SR) for different plant forms,mycorrhizal types,and different environmental treatments.Slope (RD-CT),the slope for the regression of cortex thickness with root diameter;Slope (RD-SR),the slope for the regression of stele radius with root diameter;Slope(RD-CT): Slope (RD-SR, ratio of the Slope (RD-CT) to Slope (RD-SR);AM,arbuscular mycorrhiza;EM,ectomycorrhiza;OM,Orchid mycorrhiza.

Table 2 The relationships between root diameter(RD)and cortex thickness(CT)and stele radius(SR)in each of the 26 studies with at least 3 plant species.Slope(RD-CT),the slope for the regression of cortex thickness with root diameter;Slope (RD-SR),the slope for the regression of stele radius with root diameter;Slope (RD-CT):Slope (RD-SR),ratio of Slope (RD-CT)to Slope (RD-SR).

Fig.2.The pattern of root allometry across 698 observations (a) and 512 species (b).CT: cortex thickness;SR: stele radius.

Root allometry occurs in different plant growth forms.For example,in woody species,Slope(RD-CT): Slope(RD-SR)was 7.6,and in non-woody species the Slope(RD-CT)was 3.3(Table 1;Fig.S3).The presence of root allometry in woody and non-woody species is very similar to that reported by Kong et al.(2019),which used less than half the plant species as in this current study.Interestingly,we also observed that the Slope(RD-CT): Slope (RD-SR) ratio was higher in woody vines than in non-woody vines (4.2 vs.2.5,Table S1;Fig.S5a-c).The difference in root allometry between woody and non-woody species arise mainly from higher Slope(RD-CT)of the woody species than that of the non-woody species (Table 1).This could be because woody plants rely more on mycorrhizal fungi than do non-woody species(Ma et al.,2018);hence,thicker cortex(i.e.,higher Slope(RD-CT)) is needed to accommodate more mycorrhizal fungi.

Root allometry was observed in ferns,grasses,herbs,and even the Orchidaceae bearing the thickest absorptive roots (up to 4.6 mm) (Zhu et al.,2016) among terrestrial plants (Table 1;Fig.S4d).Interestingly,we observed a slope ratio of 5.1 for Slope(RD-CT):Slope(RD-SR)in Orchidaceae species,higher than the slope ratio in ferns (1.8),grasses (2.1) and herbs (3.5) (Table 1).The difference in root allometry (indicated by the above slope ratios)in these non-woody species is mainly owed to the lower Slope(RD-SR)of the Orchidaceae species than that of other non-woody species (Table 1).This could be because Orchidaceae species,as the most evolutionarily advanced angiosperms,could have vascular conduits with higher transport efficiency than other nonwoody species;this may ensure that Orchidaceae species have a higher increase of matter transport efficiency by a smaller increase in conduit and stele diameter,i.e.,low Slope(RD-SR).Another possibility may be that Orchidaceae species usually grow in shady habitats,which have low transpiration demand,resulting in a low increase in conduit and stele diameter with increasing root diameter,i.e.,low Slope(RD-SR).

Among dominant mycorrhizal types,root allometry is observed in AM(Fig.S6a),in both woody and non-woody plants(Fig.S7a and 7b),and in the dual mycorrhizas of AM and EM plants(i.e.,roots are colonized by both AM and EM) but not in EM plants (Table 1;Fig.S5b and S5c).We observed root allometry in broadleaf EM trees but not the coniferous EM trees (Table 2;Fig.S8).Thus,the inclusion of the coniferous EM trees greatly influenced our statistical analysis of root allometry across the EM plants.However,only six coniferous EM plant species were included in our data set,and we cannot rule out the existence of the root allometry in other coniferous EM trees.If root allometry is absent in all coniferous EM trees,it might be interesting to explore the underlying mechanisms,given the distinct anatomical structures of coniferous and broadleaf tree roots (Guo et al.,2008b;Chen et al.,2016).For example,coniferous EM trees usually have vascular conduits (tracheid and sieve cells) with lower matter (water and photosynthates) transport efficiency relative to those of broadleaf EM plants,which have conduits with higher transport efficiency(vessels and sieve tubes)(Guo et al.,2008b).However,coniferous species usually have a special type of pit on the side walls of conduits,i.e.,bordered pits,which confer the conduits with higher transport efficiency than that expected only by the conduit diameter.It is possible that bordered pits rather than conduit and stele radius vary in tandem with root diameter.This would decouple stele radius (and cortex thickness) from root diameter,and hence result in no root allometry.

Root anatomical structures have been measured sparsely under different environmental change scenarios (e.g.,soil nitrogen or phosphorus fertilization,elevation of atmospheric CO2concentration,and seasonality in precipitation) (Table 1;Fig.S9),which provides an opportunity to test the consistency of root allometry.Although root allometry in some scenarios does not hold statistically(Fig.S9),this is apparently due to the inclusion of a few species with “exceptionally”sized root cortex and stele or the inclusion of some species with “exceptional” responses to environmental changes(e.g.,Fig.S9a2,a3,b2,b3,d2,d3).Here,we found that root allometry does occur after exclusion of these exceptional species(Fig.S9).These results suggest that root allometry is relatively insensitive to environmental changes.Unexpectedly,we note that cortex thickness increases slower and stele radius increases faster with increasing root diameter in the rainy season than in the dry season,causing an equal increase rate of the cortex thickness and stele radius in the rainy season (Fig.S8d2,d3).It would be interesting to test whether seasonality impacts root allometry in other plant species.

Together,the widespread existence of root allometry between root cortex and stele across different climatic zones (tropical,subtropical and temperate),ecosystem types (forests,grasslands,deserts and mangroves),mycorrhizal types (AM,some EM,AM&EM,OM),phylogenetic gradients(from ferns to Orchidaceae),and some environmental change scenarios supports the general existence of the root allometry in terrestrial plants.

4.Mechanisms of the root allometric relationship

Currently,two mechanisms have been proposed to explain why root allometry is formed.One is the nutrient absorptiontransportation balance theory (Kong et al.,2017) and the other is the carbon supply-consumption balance theory(Kong et al.,2021;Colombi et al.,2022) (Fig.3).Both theories run according to the principle of functional balance of matter (nutrients,photosynthates)transport within the root stele as well as the physical law of fluid transport in the conduits,namely the“Hagen-Poiseuille law”(Jensen et al.,2016).Here,we briefly explain the two mechanisms.

There are two paralleled vascular systems in root steles,i.e.,vessels responsible for transporting water and nutrients upward to stems and leaves,and sieves responsible for transporting photosynthates downward to meet the carbon demands of absorptive roots.According to the Hagen-Poiseuille law,both volumetric flow rates in the conduits (i.e.,water and nutrient transportation via vessels and photosynthate transportation via sieves) scale to the fourth power of the root radius (∝r4,r is radius of the absorptive roots);in contrast,even the maximum nutrient absorption (via mycorrhizal fungi in the cortex) and carbon consumption (in the cortex)scales to less than the second power of root radius(∝r2).In this case,only a much faster increase in cortex thickness than in stele radius(i.e.,the allometric relationship)can lead to afunctional balance between nutrient absorption(via mycorrhizal fungi in the cortex) andnutrient transportation(via vessels in the stele) andabalance between the carbon supply(via sieves in the stele)and carbon consumption(via the cortex and mycorrhizal fungi).

Although the above two mechanisms have been widely accepted,we should keep in mind that there are important limitations to these theories.First,thenutrient absorptiontransportation balancetheory is based on a universal mycorrhizal association in terrestrial plant roots,although many plant species have no mycorrhizal associations (Van der Heijden et al.,2015;Brundrett and Tedersoo,2018;Correia et al.,2018).Second,the two theories seem to be independent of one another,although they are both based on two interconnected vascular conduits,i.e.,vessels and sieves in roots.It is also interesting to learn about how the two theories are linked with leaves.This is because leaves function as an important nutrient sink and the source of carbon,both of which are functionally linked to absorptive roots.Thirdly,empirical evidence is still needed to test the predictions of the functional balance that underlies the formation of root allometry.In testing the predictions,we should also consider the temporal dynamics of carbon supply from leaves as well as the carbon cost in root exudates,mycorrhizal symbiosis and root turnover.

5.Implications of root allometry

5.1.The “root economics spectrum"

Similar to the well-known leaf economics spectrum (Wright et al.,2004),absorptive roots have long been assumed to vary following the root economics spectrum,a leading trait dimension in roots conveying a trade-off between nutrient acquisition and conservation (Freschet et al.,2010;Reich,2014;Kramer-Walter et al.,2016;Bergmann et al.,2020).If the nutrient acquisitionconservation tradeoff (i.e.,root economics spectrum) exists in absorptive roots,there should be a positive correlation between root diameter and root tissue density,given that thicker absorptive roots have longer lifespan(i.e.,longer conservation of the nutrients in the roots) (Guo et al.,2008a;Gu et al.,2017;Liese et al.,2019),and roots with higher root tissue density usually have lower activity in nutrient acquisition (Zadworny et al.,2017;Stock et al.,2021).Given the higher tissue density of the stele than the cortex(Valverde-barrantes et al.,2021),root allometry theoretically leads to a negative relationship between root tissue density and root diameter (Han and Zhu,2021;Kong et al.,2016;Weemstra et al.,2016).Such a prediction has been confirmed using global root trait data,which runs against the prediction from the root economics spectrum(Kong et al.,2019).This in turn suggests that root allometry may constitute indirect evidence against the existence of the traditionally recognized root economics spectrum (Kong et al.,2019).

5.2.Plant evolution

The evolution of angiosperms is closely related to the decline of atmospheric CO2concentrations since the mid-Cretaceous(Beerling and Berner,2005;Gerhart and Ward,2010).For example,the reduction of atmospheric CO2concentrations often lowers leaf photosynthesis,consequently causing a “carbon starvation” in plants.To survive in this carbon limitation condition,plants tend to increase stomatal conductance to compensatively improve leaf CO2fixation (Zhou et al.,2022;Holtta et al.,2017).However,large stomatal conductance enhances transpiration water loss,causing a physiological drought stress to plants (Khan et al.,2007;Wang et al.,2018).The surge of leaf vein density in angiosperms since the mid-Cretaceous,and hence more efficient water supply to the mesophyll cell for photosynthesis,has been considered evidence for plant adaptation to physiological drought(Baraloto et al.,2010;Feild et al.,2011;Baird et al.,2021;Yan et al.,2022).

Coordinated with the evolutionary change in leaves,thinning of the absorptive roots is regarded as an adaptation to physiological drought stress caused by the decline of atmospheric CO2concentrations (Fig.4a) (Comas et al.,2012;Chen et al.,2013;Ma et al.,2018).Root allometry may also be an adaptation to drought stress,as it effectively reduces the resistance of water and nutrients entering the root tissues and reduces carbon consumption by the root cortex (Kong et al.,2017,2021) (Fig.4b).Therefore,even though root allometry is present in non-angiosperms,e.g.,ferns,whose origin is much earlier than the Cretaceous,this adaptation has likely been instrumental in helping angiosperms survive in carbon-and water-limited environments for millions of years.From this point of view,root allometry provides insights into how roots,whole plants and even ecosystems respond and adapt to on-going geological and environmental changes.Further,the Orchidaceae plants,as the most advanced angiosperms,have the thickest absorptive roots.Therefore,the formation of root allometry in these plants may not be linked with plant adaptation to the decline of atmospheric CO2.It would be interesting to discover what root allometry implies for the plant ecology and evolution of these plants.

Fig.4.The co-variation of atmospheric CO2 concentration(a)and anatomical structures of absorptive roots in angiosperms(b)since the Cretaceous.The pattern for the change of atmospheric CO2 concentration is redrawn from Beerling and Franks (2010).The resultant stress by the change in atmospheric CO2 concentration (Comas et al.,2012) and the adaptive responses of the absorptive roots are indicated as inlets(Kong et al.,2017,2021)(see the text for the details).RD:absorptive root diameter;CT:cortex thickness,SR:stele radius.

6.Future directions

6.1.Variations of the root allometry

To date,studies have concentrated on the general pattern of root allometry,i.e.,higher Slope(RD-CT)than Slope(RD-SR).However,our analysis has shown a great variation in the Slope(RD-CT):Slope(RD-SR)ratio from the minimum of 1.2 to the maximum of 15.7(Table 2).We are still unclear why there is so large a difference in root allometry among studies,across species or even within species.This is a fascinating area that could stimulate new insights into our understanding of the root form and function.

6.2.More plant species need to be examined

Root allometry has been observed in over 500 plant species,much smaller than the total number of the vascular plant species(about 390,000) on Earth (Cantwell-Jones et al.,2022).Even in global root trait data sets,such as FRED 3.0 (Iversen et al.,2021),measurements of root anatomical traits are far less common than those of other root traits such as root diameter and root tissue density.Therefore,it is necessary to examine root anatomical structures in more plant species,especially the families with a large number of species like Orchidaceae (>20,000 species).

6.3.Environmental changes and root allometry

Few studies have examined whether environmental changes alter root allometry in absorptive roots.Moreover,these studies have only examined a few environmental changes,failing to examine interactions among these environmental factors.In addition to studies that examine how environmental changes relate to root allometry,we need to pay more attention to plants growing in naturally stressed conditions such as alpine forests,deserts,and coast environments in high salinity (Cao et al.,2023).Comparison of the root allometry among these conditions could be instructive for our understanding and prediction of vegetation dynamics under global climate change.

6.4.Plant above-ground linkage and root allometry

Plant growth and evolution depends on the functional coordination between plant above-and below-ground organs(Aritsara et al.,2022;Weigelt et al.,2021;Zhou et al.,2022).Root allometry may be accompanied by a faster supply of leaf photosynthates to the roots.However,we know little about how plants with thick absorptive roots shape their leaf structures to meet the functional balance of water,nutrients and carbon between roots and leaves.Another interesting question is whether and how root allometry is coordinated with plant reproductive organs like flowers,fruits and seeds given that the reproductive organs can compete with roots for photosynthates and with leaves for water and nutrients.Therefore,linking root allometry with plant above-ground organs is an intriguing area to deepen our understanding of the co-evolution among plant organs and the co-evolution between plants and animals for pollination or seed dispersal.

Author contributions

Y.Z.and D.K.conceived the ideas of this review paper,D.K.and Y.Z.performed the data analysis,Y.Z.wrote the first draft of this manuscript and Y.Z.,D.K.,J.C.,Q.Y.,M.W.,and Y.Z.all contributed to the editing and revision of the final version of the manuscript.

Declaration of competing interest

We declare that we have no conflict interest.

Acknowledgements

We are grateful for Dr.Meng Zhou,Oscar J.Valverde-Barrantes and WenHua Xiang for their assistance in data collection.Thanks are also given to Dr.J.Aaron Hogan (Department of Biology,University of Florida) and the other anonymous reviewer for their insightful comments and inputs for improvement of this review paper.This study was funded by the National Natural Science Foundation of China (32171746,42077450,31870522 and 31670550),Funding for Characteristic and Backbone Forestry Discipline Group of Henan Province,and the Scientific Research Foundation of Henan Agricultural University(30500854),Research Funds for overseas returnee in Henan Province,China.

Appendix A.Supplementary data

Supplementary data to this article can be found online at https://doi.org/10.1016/j.pld.2023.05.002.