Are yarn quality prediction tools useful in the breeding of high yielding and better fibre quality cotton (Gossypium hirsutum L.)?

The approach of directly testing yarn quality to define fibre quality breeding objectives and progress the selection is attractive but difficult when considering the need for time and labour. The question remains whether yarn prediction tools from textile research can serve as an alternative. In this study, using a dataset from three seasons of field testing recombinant inbred line population, Cottonspec, a software developed by the Commonwealth Scientific and Industrial Research Organisation (CSIRO) for predicting ring spun yarn quality from fibre properties measured by High Volume Instrument (HVI), was used to select improved fibre quality and lint yield in the population. The population was derived from an advanced generation inter-crossing of four CSIRO conventional commercial varieties. The Cottonspec program was able to provide an integrated index of the fibre qualities affecting yarn properties. That was compared with selection based on HVI-measured fibre properties, and two composite fibre quality variables, namely, fibre quality index (FQI), and premium and discount (PD) points. The latter represents the net points of fibre length, strength, and micronaire based on the Premiums and Discounts Schedule used in the market while modified by the inclusion of elongation. The population had large variations for lint yield, fibre properties, predicted yarn properties, and composite fibre quality values. Lint yield with all fibre quality traits was not correlated. When the selection was conducted first to keep those with improved fibre quality, and followed for high yields, a large proportion in the resultant populations was the same between selections based on Cottonspec predicted yarn quality and HVI-measured fibre properties. They both exceeded the selection based on FQI and PD points. The population contained elite segregants with improved yield and fibre properties, and Cottonspec predicted yarn quality is useful to effectively capture these elites. There is a need to further develop yarn quality prediction tools through collaborative efforts with textile mills, to draw better connectedness between fibre and yarn quality. This connection will support the entire cotton value chain research and evolution.


Introduction
Cotton is primarily produced for fibres supplied to textile mills where they are processed into yarn and various textile products.Over the value chain, fibre quality is of mutual interest to both cotton growers and processors (May and Taylor 1998;Bradow and Davidonis 2000;Kelly et al. 2015).Growers wish to produce a fibre quality that receives either no discount or a small premium.Poor quality fibre affects the value translation of field productivity into crop gross margin and reduces crop returns.To spinners, procuring cotton with an expected quality but low in price enables mills to produce the targeted quality in a more cost-controlled manner so that they can run a profitable and competitive business, as cotton remains the major cost component (Bourland et al. 2010;Constable et al. 2015;Yang and Gordon 2016).
New technologies have continuously increased the speed and efficiency of textile processes.Consumers also have an increased interest in cheap cotton brand products.Both factors have driven the ongoing changes in demands and specifications for high-quality cotton (Bradow and Davidonis 2000;Liu et al. 2010;Kelly et al. 2015).Breeding is a major contributor to improved fibre quality (Campbell et al. 2011;Clement et al. 2012;Kuraparthy and Bowman 2013;Constable et al. 2015), and continuing this effort becomes more important than ever while improving cotton productivity.
Cotton fibre quality is defined according to physical fibre properties, namely, fibre length (LEN, mm), strength (STR, g•tex −1 ) [0.980 is a conversion factor to express as cN•tex −1 ], micronaire (MIC), uniformity (UNI, %), short fibre index (SFI, %) and elongation (EL, %).Each of them has an influence on spun yarn and the finished products (Meredith et al. 1991;May and Taylor 1998;Kelly et al. 2015).High Volume Instrument (HVI) lines are widely used to measure these properties, from which breeders use for selection, merchants for determining sale prices, and spinners for deciding types of yarns to spin (Bourland et al. 2010;Kelly et al. 2015;Yang and Gordon 2016).Fibre properties are both heritable and influenced by the environment and the extent varies with traits and breeding populations (Coyle and Smith 1997;Meredith et al. 2012;Campbell and Myers 2015;Koebernick et al. 2019).Breeding for improved fibre quality is to ensure fibre properties meet the evolving demands of textile mills (Bourland et al. 2010).
Routine practice for fibre quality improvement is to assess and identify individual plants or lines with acceptable fibre properties.To warrant the translation of improved fibre properties into yarn quality, two approaches are commonly used in cotton breeding programs.One is to convert fibre property measurements into various composite indices based on their importance for spinning and use those numbers in the breeding selection.For example, the fibre quality index (FQI) or its modifications have commonly been used to assess the quality of cotton bales in textile mills (Lord 1961;Majumdar et al. 2005;Liu et al. 2011).More recently the Q-score based on normalised fibre LEN, STR, and MIC values has been developed and applied (Bourland et al. 2010).Using indices also simplifies selection and saves time (Bourland et al. 2010).Given the negative associations between fibre yield and quality, particularly via STR and LEN (Culp and Harrell 1973;Green and Culp 1990;Clement et al. 2012), the question remains on how index-based selections can affect yield improvement in cotton.The other is to directly conduct spinning tests, either using the industrial scale spinning setting (May and Taylor 1998;Long et al. 2010;Faulkner et al. 2012) or miniature lab-based spinning systems (Bradow and Davidonis 2000;Foulk et al. 2009;Liu et al. 2019).The approach is appreciated as it can verify the impact of improved fibre properties on resultant yarn quality.However, extensive use is difficult, as it is slow and expensive and is unable to assess the large number of individuals tested every season in breeding programs.However, there is an interest in resolving such challenges by developing predictive models from large textile mill datasets.
HVI measurements of cotton bales in textile mills are commonly used to assess the impact of fibre properties on yarn quality and to develop in-house predictive tools to facilitate decision-making mainly in cotton procurement and the making up of laydown for processing in the mills (Majumdar et al. 2005;Yang and Gordon 2016).Two types of predictive tools from fibre to yarn quality have been developed.One is based on the simple causational relation of fibre properties or their index numbers (e.g., FQI) with yarn properties (Ramey et al. 1977;Majumdar et al. 2005;Ureyen and Kadoglu 2006;Foulk et al. 2009;Liu et al. 2011).However, the effect of fibre properties on yarn properties can vary with the properties and spinning systems (Ramey et al. 1977;Ghosh et al. 2005;Ureyen and Kadoglu 2006).Yarn quality is also influenced by the intrinsic relation of different fibre properties as well as yarn structural factors (Faulkner et al. 2012;Yang and Gordon 2016).Therefore, these factors should be taken into account when developing predictive models (Hafez 1987;Ramesh et al. 2008;Yang and Gordon 2016).
Cottonspec is a software developed by the Commonwealth Scientific and Industrial Research Organisation (CSIRO) for textile mills to predict the yarn properties of medium to fine-count (30 tex (Ne 20) to 12 tex (Ne 50)) carded or combed ring-spun yarns based on HVImeasured fibre properties (Yang and Gordon 2016).The models were initially built and validated with large datasets from leading ring-spinning textile mills in China while taking account of the fibre properties, as well as the yarn specification (count and twist) and intrinsic yarn factors.The software has been used in textile mills in South-East Asian regions for around 10 years.However, how the models can be used to assist breeding for better fibre quality and higher-yielding cotton has not been investigated.
Nowadays, cotton bales are traded by the price schedule based on HVI-measured LEN, STR, MIC, colour, and trash content (Bradow and Davidonis 2000;Bourland et al. 2010).The schedule is updated yearly to reflect the changes for quality cotton in the market and premiums and discounts are applied to tell the potential of raw cotton processing into the end products.Using the market value of fibre quality was attempted in breeding selection but further research is required (Bourland et al. 2010).One of the limits is the current price schedule does not account for HVI-measured EL, but increased evidence demonstrates its importance for fibre resilience, consistency, and quality of resultant yarns (Benzina et al. 2007;Kelly et al. 2019;Mathangadeera et al. 2020).Thus, taking EL into account in the market value estimation of fibre quality is also important, even though its calibration standard for HVI remains to be adopted (Delhom et al. 2020).
Australia is one of the global leading cotton exporters not only producing the highest yield but supplying some of the best grades of raw cotton.Most of the Australian cotton is exported for ring spinning into medium to fine count yarns (30-12 tex) (Long et al. 2010;Yang and Gordon 2016).In this study, we examine how FQI, Cottonspec predicted yarn quality, and estimated premium and discount (PD) points are useful to assist our breeding efforts for higher-yielding and desired quality cotton when compared with our routine breeding selection based on the HVI-measured fibre quality.It is hypothesised that breeding effectiveness and progress can be further increased via streamlining and integrating the use of tools developed from the post-harvest value chain.

Test population and field experiments
All test lines in this study are recombinant inbred lines (RILs) from a multi-parent advanced generation intercrossing among four released conventional cultivars, namely, Sicot 71 (Reid 2003), Sicot 75 (Stiller 2008), Sicot F-1 (Reid 2005a), and Siokra 24 (Reid 2005b).The parents were initially paired to make two crosses in 2011, and their F 1 's were crossed to derive a segregated population with over 1 000 F 1 seeds.The F 1 to F 5 generations were advanced according to the single seed descent method (SSD) to increase the homozygosity of single plants in the population from which RILs were derived.From the F 7 generation, 256 F 6 -derived RILs were randomly selected and tested in single-row plot experiments with their four parents over the summer of 2016/2017, 2017/2018, and 2019/2020 in fields at Australian Cotton Research Institute (ACRI) near Narrabri, NSW,Australia (S30° 11' ,E149° 35').The soil is a self-mulching Vertosol classified as a fine, thermic, montmorillonitic Typic Haplustert with high clay content (Soil Survey Staff 1996;Ward et al. 1999).The experiments were planted according to a row-column design with two replicates generated with DiGGer V. 1 software.Single row plots of 12 m length and 1 m row spacing represented units in the experiment.
The experiments were managed for fertilisers, irrigation, control of pests, and weeds according to the commercial farm practices (https:// www.cotto ninfo.com.au/ publi catio ns/ austr alian-cotton-produ ction-manual), and chemically defoliated at the end of the season with thidazuron and ethephon when at least 60% of bolls were open.A single-row plot picker was used to harvest and record plot yield.About 250 g seed cotton sample was grabbed and kept at harvest of each plot and ginned with a 20-saw gin to separate lint and fuzzy seeds.The weight of the lint fraction was recorded and used to calculate lint percentage (LP) and the subsamples were kept for measuring fibre properties using a Uster 1000 High Volume Instrument following Guideline for Standardized Instrument Testing of Cotton (International Cotton Advisory Committee and International Textile Manufacturers Federation 2018).

Lint yield and fibre properties
Plot yield recorded at harvest was multiplied with LP to convert lint yield into lint per hectare (LY, kg•ha -1 ) which was used in the data analysis.The HVI testing of lint samples recorded five fibre properties: fibre length (LEN), uniformity (UNI), short fibre index (SFI), fibre strength (STR), elongation (EL), and micronaire (MIC).In cotton quality assessment, SFI is defined as the proportion of fibre in a test fibre bundle by weight shorter than 12.7 mm in length, and in HVI test, it is estimated indirectly based on fibre length measurements (Kelly et al. 2015).

Fibre quality index (FQI)
FQI is a multiplicative composite variable to transform measured fibre properties of LEN, STR, and fineness, i.e., MIC into a single value.Despite of different forms reported and used in textile mills (Lord 1961;Majumdar et al. 2005;Liu et al. 2011), in this study, a basic form according to the formula, LEN × UNI × STR/MIC, was used to calculate FQI numbers.

Predicted yarn properties
In this study, two yarn properties, yarn evenness (YCV) and tenacity (YT), for English cotton fine count (Ne) 50 (11.8tex) combed yarn were predicted by applying Cottonspec software and used to select for improved fibre properties.This count and type of yarn are typically used in the manufacture of high-quality light-weight knits and woven fabrics, e.g., for shirting.Our breeding effort is aimed at providing improved fibre properties that continuously meet the requirements of top-end yarn textile mills.Due to intellectual property protection, we are unable to disclose the algebra behind the predictions here.However, in the previous validation studies using textile mill datasets, yarn predicted values were highly correlated with measured YCV and YT in the mills where good quality control was practiced.The coefficients of determination (R 2 ) for the models ranged from 0.82 to 0.84 for YCV and 0.85 to 0.98 for YT (Yang and Gordon 2016).

Premium and discount (PD) points
In Australia, cotton bales are classed and sold according to a yearly updated Premiums and Discounts Schedule Sheet (see an example via the link: A basic guide to cotto n prici ng and quali ty-Jan 2017.pdf (cottoninfo.com.au).In the Schedule, premium and discount points are tabulated in terms of the grades of HVI-measured fibre properties including LEN, STR, and MIC combined with fibre colour and trash content.When total points are determined for cotton bales, they are converted into US cents by the rate of 100 points to one and added onto the cotton price in the ICE Futures U.S. Cottons No. 2 Futures as the bale price for sale.The system is similar to the loan scheme used in cotton sales under USA Commodity Credit Corporation (CCC) Loan Premiums and Discounts (Bourland et al. 2010).
In this study, harvest cotton samples were ginned using experimental gin after pre-cleaning and the fibre properties were tested with HVI.In the process, however, trash and colour were not assessed or recorded, therefore, we can only assume these properties were in acceptable ranges with zero premium and discount.On the other hand, there is an increased recognition of the importance of EL in both yarn spinning and fabric manufacturing, despite no calibration HVI standard available hampering its inclusion and assessment in cotton classing and trading.In this study, all samples were measured by the same single HVI in our laboratory so that different experiments could be compared without calibration (Benzina et al. 2007).In order to reflect its importance in breeding better fibre quality cotton, we developed a premium and discount points range for EL, using the range of fibre elongation values observed in G. hirsutum germplasm and also the minimal elongation (4%) required for fibre spinning and weaving (Benzina et al. 2007;Kelly et al. 2019;Mathangadeera et al. 2020) (Additional file 1, Table S1).According to the HVI measured fibre properties, we derived the correspondent premium and discount points for each trait, and the sum represents premium and discount (PD) points for data analysis and breeding selection of fibre property package.

Statistical analysis
The dataset from three season experiments were pooled together and analysed as a multi-environmental trial with a mixed model (Smith et al. 2005).Test lines, season, and their interaction were fitted as fixed and spatial variations, i.e., global and extra effects associated with the block structure of the experiment, were taken into account following our previous study (Liu et al. 2015).The empirical best linear unbiased estimates (E-BLUEs) obtained for RILs and four parents were used in retrospect selection as described below.

Comparison of selection effectiveness
Improved overall fibre quality is achieved by the selection of HVI-measured fibre properties, FQI number, Cottonspec predicted yarn properties and PD points.For each selection practice, truncating points for the traits were determined when considering the means of parents and population as well as the breeding objective of fibre quality and selection was applied to keep approximately 33% of the entire population.The resultant populations were compared for their means and variation in lint yield, fibre properties, and the proportion of RILs commonly retained, and whether yield variations were maintained as such from which further selection resulted in the improved combinations for both yield and fibre quality.For all the comparisons, the resultant population from selection based on the HVI-measured fibre properties were used as references, as that is a common practice in current cotton breeding programs.
All analyses were carried out using ASreml-R (Butler et al. 2018) and R (R Core Team 2022).

Phenotypic variation for line yield, fibre properties, predicted yarn properties, and the composite quality traits
Four parents used for crossing showed larger difference in LY, SFI, FQI, and PD points but smaller for UNI, STR, and EL with the remaining properties being within these differences (Table 1).Sicot 71 was the highest yielder but had short and coarse fibre resulting in the lowest predicted yarn quality and minus PD points (Additional file1, Table S2).In contrast, Sicot 75 and Siokra 24, despite not being the best yielders, possessed good fibre properties, especially LEN and STR, which translated to the best-predicted yarn quality as well as positive PD points.Sicot F-1 was the lowest yielder and had shorter fibre, but high STR made it a base grade fibre quality with small PD points and intermediate predicted yarn quality.
The phenotypic variations for all traits in the population are shown in Table 1.When compared with the parent means, the population exhibited higher for LY, LP, and EL, but lower for SFI, FQI, and PD points and comparable for the other traits.Variations followed normal distributions with larger variability in LY, SFI, EL, MIC, FQI, and PD points (absolute CV > 6%) and the smallest in UNI (CV < 1%).Distributions skewed to small values for LY, LP, LEN, MIC, and PD points and with a large number of individuals around the means for LY, LP, MIC, and PD points (Kurtosis > 4.5, P < 0.001) (Additional file 1, Table S3).However, transgressive segregants, measured as either higher or lower than the best or worst parent, existed for all the traits.There was a proportion of 20% to 63% of individuals better than the best parent for LY, LP, UNI, STR, EL, and PD points, 12% more for the predicted YCV, and YT and 20% with the desirable range of MIC value (3.8 to 4.9).Therefore, the intercrossing of diverse commercial varieties resulted in large segregation and recombination in both lint yield and fibre quality traits.

Interrelation of yield, lint percentage, fibre quality properties, and composite quality traits
The direction and magnitude of trait interrelations are given in Table 2. LY exhibited a moderately positive relation with LP, but only occasionally with fibre properties and their converted index number or predicted quality traits.Even when the relations existed, they were generally weak and/or may vary in direction with experimentation seasons.LP was positively correlated with EL, MIC, and YCV but negatively with LEN, FQI, and YT, however, the relation was weak (r = -0.312 to 0.309).LEN was positively correlated with UNI, FQI, STR, YT, and PD points, but negatively with SFI, EL, MIC, and YCV, in most cases being moderate (0.33<|r|< 0.66) to strong (|r|≥ 0.66) but weak (|r|≤ 0.33) with STR.UNI showed positively weak to moderate relations with STR, FQI, YT, and PD points but negative ones with SFI and YCV with magnitude ranged from strong to weak.SFI had a positive correlation with YCV but negative with STR, FQI, YT, and PD points, either strong or weak.STR was moderately positively correlated with FQI and YT but weak negatively with EL.EL was positively correlated with MIC and YCV but negatively with FQI, YT, and PD points, ranging from weak to moderate.MIC, FQI, YCV, and YT were correlated one another over the experiments.The relations were positive and moderate in magnitude, when FQI and YT paired with PD points, respectively, but exhibited strong when MIC paired with YCV and FQI with YT.The rest of the trait pairs had negative relations, ranging from moderate to strong.Strong and consistent relations among MIC, FQI, and YCV but often in opposite directions support the great influence of fibre fineness on spinning performance and yarn quality.

Comparison of selections for better fibre quality and high yield
Keeping the best 33% of the population when selecting desired fibre quality combinations resulted in the retention of 82 to 86 individual lines, regardless of which selection was used (Table 3).The selected population based on HVI-measured fibre properties was less consistent with those by others, as only 52.3% to 61.6% of individuals in the resultant populations were identical with the highest consistency being with the result of selection based on Cottonspec predicted yarn properties.The consistency increased to 67.1% or more in the resultant populations  The mean and variation of fibre properties of the selected populations were given in Table 4. Overall, selection shifted the mean of fibre properties to the target direction and range.When comparing the population selected based on HVI-measured fibre properties, the populations retained based on the other methods maintained the comparable means but extended the variability for LEN, STR and MIC, suggesting that the lines being outside the truncation ranges based on their HVI-measured fibre properties were kept when the other methods were used for selection.For UNI and EL, selection with the other methods reduced the variation in the resultant populations, particularly for EL, resulting in more lines with low EL retained, which is undesirable.
The number of lines in the selected population belonging to relative lint yield quartiles of the entire population is given in Table 5.Interestingly, regardless of the traits used in selection for better fibre quality, there were more than 54% of retained lines with yields ranking in the top two quartiles (46 to 49 lines).Furthermore, selection based on Cottonspec predicted yarn properties resulted in retaining the highest proportion of high-yielding lines commonly kept based on HVI-measured fibre properties (62%), followed by FQI numbers (58%) and the least for PD points (50%).
Relative yield distributions of the unselected and selected populations are given in Fig. 1.The mean of the resultant population based on selection using FQI numbers and PD points very much unchanged, when compared with those kept by selection based on HVImeasured fibre properties and Cottonspec predicted yarn properties.Consequently, there were fewer lines with yield potential better than the mean of the population being captured after selection with FQI and PD points.The key reason behind this was selection resulted in reducing or discarding the lines in the topyielding bins, which was more obvious in the resultant population based on FQI and best yielders (relative lint yields ≥ 109%) were dismissed.In contrast, the selected population based on Cottonspec predicted yarn properties captured a greater number of test lines in the highest-yielding bin, for examples, increased number of test lines kept for relative lint yield of 105-107% and 115%.This is consistent with the results aforementioned and suggests the merits of predicted yarn quality selection  EL /% 5.8 (4.5-7.3)5.7 (5.0-7.1)5.5 (4.7-6.5)5.6 (4.8-6.5)5.6 (4.8-6.7)MIC 4.6 (3.2-5.5)4.5 (4.0-4.9)4.4 (3.9-5.0)4.4 (3.9-4.7)4.4 (3.9-4.9) for breeding high yielding and better-quality cotton, although the routine selection practice based on HVImeasured fibre properties performed comparably well.

Discussion
A grower's cotton will be eventually transformed into the end-use products, yarn and fabrics that matter to our daily life.There is a long and strong advocacy to link yarn quality with its farm-gate fibre quality to assist and direct the on-farm research to enhance crop performance, resilience, and profitability (Bradow and Davidonis 2000;Bourland et al. 2010;Liu et al. 2019).Using direct spinning tests is challenging, particularly in terms of the considerable input of time and labour to prepare and process sufficient lint before conducting the spinning tests, regardless of what scale settings, e.g., industrial, pilot or miniature spinning systems (Meredith et al. 1991;Foulk et al. 2009;Long et al. 2010;Faulkner et al. 2012).
For breeding programs, where selections are made on hundreds or thousands of single plants and test lines, it is prohibitive, given the short time available before breeders can use the information to make decisions on selection.This study tested the alternatives using the tools developed in textile mills to quantify the end-use value of fibre bales for breeding selection.In contrast to previous research aiming to develop in-house prediction models using limited field experiments and spinning tests conducted for breeding selections (Meredith et al. 1991;Benzina et al. 2007;Faulkner et al. 2012;Liu et al. 2019), we assessed the value of three widely used tools existing in the post-harvest value chain research for cotton breeding selection.These tools are FQI numbers (Majumdar et al. 2005;Liu et al. 2011), PD points from Premiums and Discounts Schedule used in today's cotton classing and trading, and predicted yarn properties from Cottonspec, which is a CSIRO-developed software based on a large ring spinning mill dataset (Yang and Gordon 2016).When they were used for sequential selection for simultaneously improving fibre quality and yield in the test population (Table 1), selection based on Cottonspec predicted yarn properties for fibre quality resulted in a population with better yield when compared with those retained by the other tools examined (Table 5 and Fig. 1).When compared with the resultant population by the routine selection based on HVI fibre properties, a large proportion of individuals in the retained populations was identical, and more interestingly, a greater number of the highest yielding RILs were kept under selection based on Cottonspec predicted yarn properties (Fig. 1).Thus, Cottonspec is a good and rapid tool for selecting improved fibre quality while maintaining good yield variation and potential in the resultant population.
Australian cotton is mainly supplied to ring spinning mills that produce high-end quality yarns (Long et al. 2010;Yang and Gordon 2016), under which HVImeasured LEN, UNI, STR, and MIC, strictly speaking, fibre fineness, are known to play a great influential role on ring-spun yarn quality (Ramey et al. 1977;Meredith et al. 1991;May and Taylor 1998;Bourland et al. 2010; Table 5 The number of lines in the selected populations for better fibre quality belonging to different quartiles of relative lint yield of the unselected population calculated against a parent, Sicot 71 under different selection methods and the proportion of their commonality with those kept by the selection based on HVI-measured fibre properties The values in brackets represent the proportion /%  Kelly et al. 2015).Better agreed outcomes in selections between HVI-measured fibre quality and Cottonspec predicted yarn properties suggest that the latter can fairly describe and reflect the core breeding objective for fibre quality in the CSIRO cotton breeding program.However, that is not the case for selection based on FQI and PD points, as their resultant populations after selection contained a reduced number of individuals common to those selected based on HVI-measured fibre properties, and they kept more lower yields while excluding individuals with higher yields (Tables 3 and 5; Fig. 1).Therefore, both FQI and PD points in the current forms are not useful in breeding selection for the improved combinations of yield and fibre quality.The shortfall in the ability for FQI and PD points to be used in selection can be explained as: 1) The relative importance of HVI fibre properties.In ring-spun yarns, fibre strength is the number one influential factor on yarn tenacity, but fibre length and fineness also played an important role either independently or by their interactions with other properties in the yarn structure (Ramey et al. 1977;Meredith et al. 1991).Unfortunately, this complex is not fully captured in the FQI formula.Majumdar et al. (2005) demonstrated the use of a multiplicative analytic hierarchy process to find and assign the different power for involved fibre properties in the FQI formula.The power could reflect their relative importance in spinning processes.When modified FQI was used in the calculation, the resultant numbers showed a much stronger association with ring-spun yarn tenacity.PD points in this study represent the net sum premium and discount of four fibre properties, i.e., LEN, STR, EL, and MIC.The points of each trait can reflect their relative importance.However, PD points were not distributed normally where more than 41% population concentrated in the zero to positive ranges with small variability (Table 1; Additional file 2, Fig. S1), which is similar to the distribution of fibre premium and discount values calculated from USDA Commodity Credit Corporation Loan Value (Bourland et al. 2010).This characteristic can limit the discrimination of breeding lines for selection and may also unintentionally exclude high yielding RILs.2) Discount the impact of important traits and conditions on estimation.SFI is one of the major HVImeasured fibre properties associated with YCV.Short fibres are increased during ginning, cleaning, and prior-spinning preparations, with the accumulated broken long fibres over mechanical processes (Yang and Gordon 2016).In both the FQI number and PD point estimations used in this study, SFI is not included.However, there are studies that different FQI formulas counting for the SFI effect could improve their predictiveness for both yarn tenancy and evenness (Majumdar et al. 2005;Liu et al. 2011Liu et al. , 2019)).In Cottonspec, the multiple linear regression model for yarn evenness predictions counts on SFI contribution, and via yarn evenness, it influences yarn tenacity based on yarn weak link theory (Yang and Gordon 2016).The current PD points version does not consider both UNI and SFI, as they are not considered or valued in current cotton classing and trading systems.3) Both FQI and PD points are essentially transformed composite variables based on solely HVI measured fibre properties or price points, thus they completely ignore the effect of spinning setting, for example, twist, yarn count, and yarn structure on yarn properties (Hafez 1987;Ureyen and Kadoglu 2006;Yang and Gordon 2016).It is not surprising that they were not as competitive and effective as Cottonspec predicted yarn properties in breeding selection.
Moderate relation between lint yield and LP suggests higher LP is one of the key contributors to higher yielding in the population.This agrees with a trend in cotton breeding that genetic progress made for cotton yield continuously relies on selecting for high lint fraction of harvest yield (Campbell et al. 2011;Conaty and Constable 2020).In the population, the proportion of RILs with LP higher than the highest parent, Sicot 75, was dominant (30%), demonstrating how single parents can drive and skew trait segregation and distribution.Sicot 75 is known for its small seeds with weak seed and seedling vigour.Similar phenomena were also observed in the USA breeding programs (Snider et al. 2014;Dowd et al. 2018).Therefore, for those selected lines with improved yielding and better fibre quality, further work requires to see whether they possess seeds meeting the standards of planting seed quality for growing commercial crops.Our future paper will focus on how to achieve the right balances between seed size, lint fraction, and fibre quality at seed levels based on a good understanding of seed yield trait interrelations while breeding for high lint yielding and better-quality cotton.
That absence of a strong inverse relation between lint yield, HVI-measured fibre properties and other composite traits in the population was interesting.The finding explains why the resultant populations after truncation selection for improved fibre properties still contained sufficient variation for lint yield from which high lint yields can be captured and why a large proportion of the selected populations remain common despite of different trait forms, combinations or nature used in the initial screening for better fibre quality.However, this disagrees with previous studies in which undesired relations were common, for example, yield with LEN and STR (Culp and Harrell 1973;Clement et al. 2012;Koebernick et al. 2019).Australian cotton germplasm is known for its superior fibre quality, especially LEN and STR.The exploitation and utilisation of Australian germplasm in USA cotton breeding programs in recent decades were reported to one of the key elements behind the improved LEN and STR of the released cultivars (Kuraparthy and Bowman 2013).Clement et al. (2012) further demonstrated the negative relations between yield and STR and LEN were weaker in high-quality advanced breeding materials bred and tested in Australia than USA.Taken together, we expect future simultaneous breakthroughs for high productivity and better fibre quality will come more often from Australian breeding efforts and subsequently, that will continuously help global cotton genetic improvement.Furthermore, LEN and MIC are essentially two components of LP, as one represents fibre volume, and the other is fibre weight.The negative and unfavourable relations of these traits suggest LEN and MIC need to be particularly watched in breeding for high yield with improved fibre quality.Using alternative traits, such as fibre density (Clement et al. 2014), may avoid the selection pitfall of short and coarse fibre quality in combination with high yielding.HVI-reported MIC is known to be the outcome of both fibre fineness and maturity (Kelly et al. 2015), breeding programs are encouraged to use the instruments, for examples, AFIS and FMT, to measure and screen fibre fineness and maturity.
The nature and extent of interrelations between HVImeasured fibre properties suggest selection for the combination of long, higher length uniformity, strong, and fine fibres was feasible in the population, but there were challenges for combining long, strong, and fine fibres with increased fibre extensibility (EL) (Table 2).It is clear when EL is ignored (FQI and Cottonspec predicted yarn properties) or given little emphasis (PD points) (Table 4), selections resulted in discarding the RILs with best EL while retaining the RIL with low EL.In contrast, that was ultimately avoided when using the routine selection practice, as individual HVI-measured fibre property was checked and scrutinised.Better EL has been emphasised recently for its important roles in textile processes to reduce stoppages and increase processing efficiency as well as improve yarn evenness and tenacity (Ureyen and Kadoglu 2006;Faulkner et al. 2012;Kelly et al. 2019;Mathangadeera et al. 2020).Therefore, to avoid the reductions for EL when using the other tools in selection for fibre properties, their new versions need to appropriately account for the EL effect.
In this study, EL was one of the uncalibrated HVI measures from our laboratory but was done by a single HVI.Despite subjecting to the influence of experiment seasons, EL was less affected by genotype × season interaction (broad sense heritability (H 2 ) = 0.77), which is very similar to LEN (H 2 = 0.88) and STR (H 2 = 0.82).Thus, EL can be measured as reliable as other important fibre quality traits, in agreement with previous studies for fibre properties (Meredith et al. 1991(Meredith et al. , 2012)).More interestingly, despite no large difference existed among the four parents, the variation for EL in the population was substantial where there was a proportion of 63% of the population with EL exceeding the best parental value (5.6).Therefore, our evidence supports Delhom et al. (2020)'s conclusions that the uncalibrated EL from a single HVI can be used for the selection of improved fibre extensibility in the population.The negative weak relation between fibre strength and elongation in this study is consistent with many previous studies (May and Taylor 1998;Benzina et al. 2007;Faulkner et al. 2012).
Interrelations between HVI-measured fibre properties with FQI, Cottonspec predicted YCV and YT and PD points suggest that low EL and high MIC and SFI results in generally low FQI numbers, and shorter, poor uniform and coarse fibres with high SFI would lead to more uneven and weak predicted yarn properties.Longer and more uniform fibres would reward higher PD points, i.e., premiums, while lower EL and higher MIC reduced PD points, i.e., discounts.All these evidences describe the relative importance of appraising fibre quality using estimated composite or predicted values.A very strong interrelation between FQI and predicted YCV and YT explains why they can pick up a large proportion of the same individuals in the populations when used for fibre quality selections (Table 3).
Cottonspec is a yarn quality prediction software developed using large datasets from leading ring-spinning mills in China.It is particularly relevant to ring-spun medium to fine-count yarns that are in high demand in the market.Cautions should be taken when used in breeding cotton aiming at the other spinning systems.The current software version does not count for the effect of fibre EL on yarn properties; further research and upgrade are needed.Following the same token, to ensure the timely reflection of the ongoing evolutions in market demands for quality cotton fibres, future research should be continued to develop more robust predicted models through exploiting big datasets from spinning mills using more advanced data analytic tools, for example, artificial intelligence and machine learning.Overall, this study does exemplify how the tools developed from the post-harvest value chain can be potentially used in cotton improvement and how they can support and bridge the connectedness and impacts of the field to post-harvest research and development efforts in the cotton value chain.

Conclusions
In this study, three tools existed in cotton post-harvest value chain were examined for their effectiveness on selecting desired fibre properties and subsequently impact on capturing higher yielders in a four-way cross derived RIL population.Evidently, intercrossing multiple and diverse parents can lead to the emerge of recombinants with improved yield and fibre quality in the population.Selection first for desired fibre properties would truncate the variation of yield of the selected population, however, the extent varied with the tools of selection.Cottonspec predicted yarn quality exceeded both FQI and PD points in terms of retaining high yielders, and also showed good competitive to the routine selection practice based on HVI-measured fibre quality.Therefore, Cottonspec predicted yarn quality can be used in preceding selection as a replacement of, or in a complementary manner with the routine selection practice used HVI-measured fibre properties.The significance of this change in the practice should help breeders much quicker glance at how the improvement made at fibre levels impact the quality of end-use products.They can use that information to define fibre quality breeding objectives and in a broader sense, fibre quality goals for the entire cotton value chain from field researchers, growers, traders, and processors so that they can work and invest together in the directions from which the resultant development in technologies can advance global cotton production and utilisation.

Fig. 1
Fig. 1 Distribution of relative lint yield against the highest yield parent, Sicot 71, in the populations prior-and post-selection for fibre quality based on HVI-measured fibre properties (top), fibre quality index (FQI, 2 nd ), Cottonspec predicted yarn properties (3 rd ) and premium and discount (PD) points (bottom).Down triangles represent relative yields of the lowest and highest yielding parents; Squares represent the means of the unselected and selected populations

Table 1
Descriptive statistics of yield, lint percentage, fibre properties, fibre quality index, premium and discount points, and predicted yarn properties for the parents and population CV Coefficient of Variance, LY lint yield, LP lint percentage, LEN fibre length, UNI fibre uniformity, SFI short fibre index, STR fibre strength, EL fibre elongation, MIC fibre micronaire, FQI fibre quality index, YCV yarn evenness, YT yarn tenacity, PD points Premium and Discount points

Table 2
Correlation coefficient and significance of yield, lint percentage, fibre quality properties, predicted yarn properties and other composite fibre quality values from mean estimates of analysis of individual and combined season experiments LY lint yield, LP lint percentage, LEN fibre length, UNI fibre uniformity, SFI short fibre index, STR fibre strength, EL fibre elongation, MIC fibre micronaire, FQI fibre quality index, YCV yarn evenness, YT yarn tenacity, PD points Premium and Discount points, ns non-significant * , ** , *** represent significance at P < 0.05, 0.01, and 0.001, respectively

Table 3
The number and proportion of common lines retained by different selection methods for desirable fibre quality under selection intention of keeping 33 percentage of the population

Table 4
Comparison of six fibre properties between the unselected and selected populations from four different methods for best fibre quality combinations