Hair curvature underpins structural diversity and function in mammalian coats, but what causes curl in keratin hair fibres? To obtain structural data to determine one aspect of this question, we used confocal microscopy to provide in situ measurements of the two cell types that make up the cortex of merino wool fibres, which was chosen as a well-characterised model system representative of narrow diameter hairs, such as underhairs. We measured orthocortical and paracortical cross-sectional areas, and cortical cell lengths, within individual fibre snippets of defined uniplanar curvature. This allowed a direct test of two long-standing theories of the mechanism of curvature in hairs. We found evidence contradicting the theory that curvature results from there being more cells on the side of the fibre closest to the outside, or convex edge, of curvature. In all cases, the orthocortical cells close to the outside of curvature were longer than paracortical cells close to the inside of the curvature, which supports the theory that curvature is underpinned by differences in cell type length. However, the latter theory also implies that, for all fibres, curvature should correlate with the proportions of orthocortical and paracortical cells, and we found no evidence for this. In merino wool, it appears that the absolute length of cells of each type and proportion of cells varies from fibre to fibre, and only the difference between the length of the two cell types is important. Implications for curvature in higher diameter hairs, such as guard hairs and those on the human scalp, are discussed.
Hair, along with milk production, is a definitive phenotypic trait for mammals. Hair appears to have evolved from scales of synapsid reptiles around 200-million years ago (Alibardi, 2006; Maderson, 2003), and hair was already a well-established trait at the very dawn of mammalian evolution because most aspects of hair biology are highly conserved across all existing mammals, including general fibre morphology.
Mammalian coats are typically composed of mixtures of differently functional hair types that are arranged in a three-dimensional (3D) structure, or pelage. The emergent properties of the pelage, and its modification over time (typically seasonal), forms the functional phenotype upon which evolutionary selection occurs, and the adaptability of which has enabled mammals to colonise a wide range of habitats, including climatic extremes (Ryder, 1973). In a typical pelage, for example that of a deer, straight high-diameter guard hairs scaffold a mass of narrow diameter curly underhairs (Brunner and Coman, 1974; Woods et al., 2011), with the combined effect of thermoregulation and mechanical protection. Single-fibre properties, in particular length, diameter and curvature, are important in defining the emergent properties of the pelage. The underpinning mechanism that connects structures at the protein and microstructure level to single-fibre curvature is not well resolved. Here, we use merino wool as a model for understanding the basis of mammalian hair curvature. Whereas wild sheep have a typical pelage of high diameter guard hairs and low diameter underhairs, in that of merino domestic sheep, both guard and underhairs have been selected to be of low diameter and high curvature (Harland et al., 2015; Ryder, 1964).
It has long been known that, in wool fibres, there is an approximate correlation between the amounts and distributions of different wool cortical cell types and fleece crimp [the primary parameter of the sheep pelage structure, observed as a wave within a clipped staple (tress/tuft) of wool typically containing thousands of fibres]. Crimp is an emergent effect of single-fibre curvature. Fibre curvature is of two sorts: intrinsic curvature and imposed curvature. Intrinsic curvature is built into a fibre during its development in the follicle, and it is this curvature, sometimes called inherent curvature, to which an undamaged fibre returns when relaxed in water and dried without mechanical constraint (Fish et al., 1999). Imposed curvature is any chemically or physically induced change to intrinsic curvature; for example, from drying a fibre in a constrained state. Here, we investigate intrinsic curvature (Fig. S1).
The observation that high-crimp wools tend to have a cortex in which the cell types are bilaterally arranged, with the paracortical cells always on the inside (or concave) half of the fibre curvature, led to the theory that it is the relative proportion and distribution of the orthocortex and paracortex patches that underpins intrinsic single-fibre curvature (Fraser and Rogers, 1954; Horio and Kondo, 1953). The theory has some experimental support from light- and electron-microscopy studies measuring average cell type distribution against average single-fibre curvature (Orwin et al., 1984; Snyman, 1963). The theory explains that, during keratinisation in the wool follicle, a fibre goes from a wet to a dry state, and that this affects the macrofibrils making up the cells of each cell type differently. Due to a greater number of intermediate filaments (IFs) in orthocortical macrofibrils that are initially tilted away from the fibre axis, longitudinal extension occurs during drying, and this extension in the orthocortex is believed to exceed that of paracortex, which is composed of macrofibrils in which the IFs are aligned along the fibre axis. The theory predicts that, in order to relieve internal strain energies, the wool fibre will curl in the direction of the less-extended paracortical side (Brown and Onions, 1961; Dobozy, 1959; Munro and Carnaby, 1999). We refer to this theory as the ‘cell type distribution model’.
The most convincing and well developed alternative theory is what we will call the ‘cell proliferation model’, and was established from observations of the relationship between follicle morphology and resultant fibre morphology (Auber, 1951). The proposed mechanism of the cell proliferation theory is that there is an asymmetry in the rate of cell division on either side of the follicle's dermal papilla. Consequently, the more prolific side of the follicle produces what will become the cortex closest to the outside (convex side) of the fibre's curvature. The cell proliferation model explains both intrinsic fibre curvature and the observation that there is a possible association between fibre curvature and degree of follicle bending (Nay and Johnson, 1967). In the cell proliferation model, the cell types per se are not directly involved in defining intrinsic curvature; there is just more material on one side of the fibre than the other. More recently, an advanced version of this theory has been developed in which resultant curvature is a trade-off between asymmetric cell production and differential onset of keratinisation on either side of the fibre axis within the follicle (Hynd et al., 2009).
We directly tested the two models using confocal microscopy by measuring cell type proportions and cell numbers and lengths in small pieces of fibre (snippets) of known specific intrinsic curvature but no apparent torsion, and found that, although neither theory fully explains curl, the cell types play a critical role.
MATERIALS AND METHODS
Six adult female merino sheep were sampled from a research flock aged between 3 and 10 years (Table S1). Samples were collected in June 2011 (winter) from sheep housed on a commercial farm in North Canterbury, New Zealand that had been on the farm under uniform feed and environmental conditions for 6 months prior to sampling. All samples were taken from a small mid-side patch using electric animal clippers and cut within a few millimetres of the skin surface. Sheep housing and wool collection procedures complied fully with New Zealand law and were reviewed and approved (AE Application 12188) by the AgResearch Invermay Animal Ethics Committee.
Sample cleaning, relaxation to intrinsic curvature and selection
Staples of wool were cleaned with heptane, ethanol and ethoxylated nonylphenol (Teric GN9) (0.15% in 60°C), and rinsed in ultrapure water. Two-thirds of each staple closest to the tip was then cut away and discarded. Individual fibres were removed by handling with fine forceps at their tip end only and attached by the root end to a strip of double-sided tape along one edge of a microscope slide. Only fibres in good condition with a regular appearance were selected. Each fibre was then cut down further so that only 20 mm closest to the root was retained, tethered at one end to the slide. Slides of fibres were submerged for 18 h in 500 ml of degassed ultrapure water containing 400 µl of 0.15% Teric GN9 (to disrupt surface tension bubbles on fibre surfaces). Individual fibres were removed by breaking them with forceps where they were attached to the tape, transferred to a dry sheet of cardboard surface on the baseplate of a high-vibration shaker (Retsch AS200; Retsch, Germany) and air dried at ambient conditions (20°C, 50% relative humidity) at a frequency and amplitude that kept the fibre mobile on the surface for 60 s. Fibres were stored in small Petri dishes until further cutting into snippets. Snippets were cut under a magnifier with dark-field illumination on a clean ceramic surface with a GEM single-edge razor blade. Individual snippets that had uniform curvature and no torsion were transferred, using static adhesion, onto a microscope slide. Using a dissecting microscope and transmitted light, snippets were carefully anchored at one end into a thin line of chemically inert acetic cure silicone cement (Selleys Glass Silicone Sealant; Selley's Ltd, Australia) to lie flat on the glass surface. Silicone cement was cured for at least 72 h at room temperature before any further processing. Each snippet was examined using a stereo light microscope, and those that were not lying flat against the slide and those in which curvature had become non-uniform were marked and not used in the study (about 50%). A total of 717 snippets were mounted and five from each sheep were selected for imaging. With no previous guide as to appropriate sample size, we chose a manageable target of 30–50 fibres for each experiment, and concentrated on selecting snippets with most uniform curvature and least torsion, thereby reducing two major potential sources of measurement error.
Measurement of intrinsic curvature
Slides were first dried over dry silica gel under a low vacuum for 24 h (measured 5–9% relative humidity) then immediately scanned dry at 2400 dpi in an Epson 3200 flatbed scanner (Fig. 1A). Measurements (deg mm–1) were made from scans by fitting a circle to the region of curvature later used for confocal measurements (AnalySIS Five software; Olympus Imaging Systems, Germany). Curvature measurements were made only after confocal imaging so that the precise region of the snippet from which data were obtained could be identified. Following scanning, slides were stored under vacuum in the dark at room temperature until staining.
Staining for microscopy
Two sets of slides were produced for each sheep. One was stained to highlight orthocortical and paracortical differences, and the other stained to highlight cell boundaries.
For cell-type staining, we used an en bloc modification of the approach used in Bryson et al. (2009). Slides were placed in a staining trough, covered with a solution of 0.002% fluorescein sodium in 4% benzyl alcohol and left in the dark at 45°C for 18 h with gentle rocking. Fluorescein was drained off and slides washed briefly (30 s) twice in ultrapure water and then once in 4% benzyl alcohol before being immersed in a solution of 0.0005% sulforhodamine 101 in 4% benzyl alcohol in a staining trough and left for 1.5 h at 45°C in the dark. Finally, the stain was drained off the slides and they were washed briefly (30 s) three times in ultrapure water, blotted gently with tissue, dried carefully with compressed air from a can and stored under low-vacuum over dry silica gel in darkness.
For cell boundary staining, we followed a method reported by Sideris et al. (1990). Slides were sealed in an upright Coplin jar with a screw-top lid covered in a solution of 0.2 mg ml−1 rhodamine B in 9:1 ethanol:methanol for 18 h at 60°C. Slides were then rinsed for 20 s in perchloroethylene and dried in the dark for 30 min under a stream of nitrogen gas (Fig. S2). Slides were stored under low-vacuum over dry silica gel in darkness until use.
Confocal imaging and analysis
Slides were cover slipped with Dako anti-fade mounting medium and imaged using a laser scanning confocal microscope (Leica TCS SP5; Leica Microsystems, Wetzlar, Germany). Fibre snippets previously selected for imaging were located in bright field mode and then examined using laser confocal mode.
For cell type imaging, excitation wavelengths were tuned to the 488 and 561 nm laser lines (argon laser) for fluorescein and sulforhodamine, respectively. Snippets were first oriented and centred at a point >100 µm from the cut end and where staining was well defined. Then x–z images (optical cross-sections) were made directly at regular 25 µm intervals until a span of 100 µm had been covered. The orientation of the sample was periodically checked and adjusted so that each x–z section was transverse to the fibre longitudinal axis. Each x–z image from each snippet was enhanced by adjusting brightness and contrast. Fibre perimeter was manually traced and the same process was used for images of each cell type to segment the orthocortical and paracortical regions into red and green, respectively (Fig. 1B). Measurements of orthocortical and paracortical regions were made using automatic image analysis tools (AnalySIS Five) (Fig. S3).
For cell boundary samples, a Leica 63× glycerol immersion objective (HCX PL APO 1.30NA) with laser tuned to sulforhodamine excitation spectrum was used to collect stacks through selected snippets at a point of uniform curvature. Measurements of cell length were made from Leica's LAS AF software (V 2.6 2011; Leica Microsystems, Wetzlar, Germany) directly from stacks. Fibre diameter was measured from the projected image of all slices of the stack combined into a single image and is the diameter across the fibre's plane of curvature.
Cell length was measured by identifying the location of cell ends by moving through the image stack and by changing image brightness/contrast. A straight line was then drawn between cell ends (Fig. 1C). It was possible to measure cell lengths and edges throughout the width of the fibre, but we chose to only measure cells within two full cell widths of the fibre cuticle (Fig. S4). In preliminary studies, we established that this ensured that those cells measured on the outside of the curvature were always orthocortex type and those on the inside always paracortex type, and avoided any subjective bias that might have arisen from trying to discriminate cell type based on faint differences in stain colour and intensity.
To examine the hypothesis that curvature is caused by higher numbers of cells end-on-end along the side of the fibre on the outside of curvature, it was necessary develop a metric to enable sample-to-sample comparison measured in a number of cells per unit fibre length. This was difficult because the number of cells end-to-end along one or the other side of the fibre could not easily be directly measured by counting. There were multiple reasons for this: (1) the staining in many fibres did not clearly indicate the edges of every cell; (2) cells were often tilted or not straight; and (3) the outer surfaces of the fibre were often lumpy.
Our solution was to calculate the number of cells on the inside and the outside of the curvature based on geometric data that we could collect accurately from each sample (Fig. 2). Our calculation of cells/length was based on the assumption that the average percentage overlap for neighbouring cells was similar in both the orthocortex and paracortex, something that is supported by earlier transmission electron microscopy (TEM) and current confocal observations. We based our calculation on a 250 µm length of fibre because this was a similar scale to the region over which cell lengths were measured.
Statistical software and data availability
All statistical analyses were carried out on GraphPad Prism V4.03 (San Diego, CA, USA). Measurement data are available in online repository figshare (https://doi.org/10.6084/m9.figshare.5500873.v1).
Cell type measurements
In all cases, cell type regions were clearly differentiated in fibre cross-sections, cell type regions were bilaterally distributed, the paracortex was always within the half of the cortex closest to the inside of the fibre curvature, and the orthocortex was predominantly within the cortex half closest to the outside of the curvature.
The selected snippets (n=5 each from 6 sheep) covered a wide range of diameters and curvatures, and there was a negative trend between curvature and diameter (Fig. S5, dataset 1).
With the exception of one snippet, the orthocortical cell region always had a greater cross-sectional area than the paracortex. Over the range of curvatures in this study, the cross-sectional areas of both the orthocortical and paracortical regions increased with increasing fibre diameter (Fig. 3A). However, there appeared to be no clear relationship between curvature and the area of each cell type (Fig. 3B) (Table S2, datasets 2 and 3).
We considered that it may not be the cell type alone that may relate to curvature but other measures relating to the proportion of cell types. However, none of the following measures appeared to correlate with curvature (Table S2): individual proportion within a snippet (Fig. 4A); the net difference between the two types as a proportion of the overall diameter of the snippet (Fig. 4B); or the distance between the central points (centre of mass/gravity) of each cell type region (Fig. 4C,D).
All virtual (x–z) sections of all snippets in this study were bilaterally symmetric and divided into single orthocortical and paracortical patches. We investigated whether the shape of the patches (as opposed to the size) related to fibre curvature. In transverse section, the region of paracortex was typically observed to be semi-circular or elliptical. In the case of a paracortical region of elliptical shape, ‘arms’ of orthocortex wrapped partly around the paracortex close to the cuticle at its narrow ends. We considered that the placement of orthocortical material alongside paracortical material close to the cuticle on the inside curve side of the cortex could potentially influence snippet curvature. We studied this possible effect by measuring orthocortical convexity (the difference between the area of the shape's convex hull and its actual area, i.e. a measure of how convoluted a shape's profile is) (Fig. S2) and comparing it with curvature (Fig. 5A), but found no relationship (Table S2). In addition, we examined the aspect ratio (Fig. 5B) and the circularity (a measure of roundness, where a perfect circle=1) of the regions (Fig. 5C), but found no correspondence between shape of cell type regions and curvature.
Cell length/number measurements
In all stacks of snippets stained for cell boundaries, we were able to use difference in the intensity of staining across the fibre to differentiate between orthocortex and paracortex. This allowed us to confirm that all samples had cell types that were bilaterally distributed and that the paracortex was always within the half of the cortex closest to the inside of the fibre curvature. The selected snippets (n=5 each from 6 sheep and an additional n=10 each from 2 sheep) covered a wide range of diameters and curvatures (Fig. S6A) and were highly variable between different snippets from single sheep. Again, the relationship between diameter and curvature was found to be noisy (Fig. S6B) but negatively correlated (Table S3).
Measurements of cell length resulted in a sizable data set (n=344 cells) across the 50 wool snippets. The data from the inside and outside of the curve (Fig. 6A) were normally distributed (Kolmogorov–Smirnov test), and cells from close to the inside of the curl were significantly shorter than cells close to the outer side (unpaired two-tailed t-test, P<0.0001, t=12.84, d.f.=342).
The two distributions of cell lengths had considerable overlap (Fig. 6B). A key finding here is that in every fibre we examined the average length of orthocortical cells for that fibre was longer than that of paracortical cells from the same fibre. This is clearly seen when the data are plotted as pairs (Fig. 6C).
We observed that the average length of cells measured for a single fibre increased with fibre diameter irrespective of position with respect to curvature (Fig. 7A) (Table S3). Against curvature, only the paracortical cells showed a statistically significant relationship, in which higher curvature was associated with shorter paracortical cells (Fig. 7B). The difference in cell length (outside minus inside length) for each snippet increased with increasing fibre curvature (Fig. 7C), and this was also the case when the data were normalised against fibre diameter (Fig. 7D).
We estimated the number of cells along 250 µm of fibre along the inside and the outside of each fibre by combining data on the curvature, diameter and cell length measurements from each snippet (see Fig. 2). From these data we determined the ratio of cell numbers on each side of the fibres and analysed this with respect to fibre diameter (Fig. 8A) and curvature (Fig. 8B). Both results indicate that there are generally more cells on the inside of curvature of a fibre despite the shorter arc length on the inside of a curved fibre. Although increases in diameter do not appear to affect this relationship, increasing fibre curvature appears to increase the relative difference in cell numbers on the two sides of the fibre (Table S3). It is, however, a weak statistical relationship and may be associated with the tendency for cells on the inside of the fibre curvature to become shorter with increasing curvature (Fig. 7B).
Our in situ measurements of cell types and cell lengths for regions in which the curvature at the point of measurement was also measured allowed us to directly test current theories of fibre curvature, which were built on indirect data. However, it should be noted that the 3D shape of a single fibre is generally a product of fibre curvature and torsion. Our findings are focused explicitly on the curvature component of single-fibre shape.
One of the current models of fibre curvature, the cell type distribution model, proposes that curvature is caused by the relative amount and distribution of orthocortical and paracortical cells, because orthocortical cells will be uniformly longer than paracortical cells. The hypothesis is that this length difference is caused by orthocortical extension, as predicted by Munro and Carnaby (1999); their work used mathematical modelling to predict how macrofibrils with differently arranged intermediate filaments would change in length because of the drying and shrinkage that occurs during fibre hardening. Munro and Carnaby concluded that macrofibrils with a double-twist arrangement, typical of orthocortical cells, should extend significantly more due to lateral shrinkage than macrofibrils with filaments arranged in parallel, typical of paracortical cells. Microbeam small angle X-ray studies have shown that average intermediate filament angle increases on the outside of curvature of both merino wool (Kajiura et al., 2005) and in the much higher diameter curly human scalp hair from a range of ethnic origins (Kajiura et al., 2006), suggesting that average macrofibril double-twist intensity (Harland et al., 2014) varies across the cortex in a way that is consistent with longer cells on the outside of the curvature.
Our data strongly support this model, but with some key additions. As it stands, the cell distribution model assumes that all orthocortical cells have a similar distribution of length irrespective of fibre or animal, and that the same is true for paracortical cells. This assumed universal difference in cell type length is what should allow curvature to be predicted from proportions and distributions of the cell types within the model (Liu and Bryson, 2002). We would expect this relationship to be particularly clear in fibres that have the combination of especially high curvature, clearly differentiated cell types and simple cell type distributions, such as in merino wool. Although orthocortical cells were significantly longer than paracortical cells (Fig. 6A,B) when averaged across fibres and animals, we found no evidence to support the assumed universal link between the relative amount of each cell type and curvature (Fig. 4). Our data indicate that cell length differences between orthocortical and paracortical cells that correlate with curvature are location specific. That is, two hairs that have the same sized patches of paracortex and orthocortex could have different curvatures because the fibres differ in the relative length of their orthocortical and paracortical cells.
Our results were obtained from merino fibres that had low diameter and high curvature, cell types that were clearly differentiated into orthocortex and paracortex, and that showed clear bilateral symmetry in their distribution. Further study will be required to establish how this variation between fibres and animals extrapolates to fibres with less differentiated morphology. In particular, the relationship between curvature, diameter and cell length will require investigation in fibres with higher diameters. For example, it is possible that variation in cell type length may decrease with increasing diameter, leading to some fibres in which the proportion of cell types does correlate with curvature – higher diameter fibres such as human scalp hairs appear to have cell types that are less clearly differentiated than those found in merino wool (Bryson et al., 2009; Harland et al., 2014).
Further study will also be required to understand the mechanism that controls cell length in developing fibres. Two factors in particular that may affect cell length are the process of cell reshaping and ordering that occurs in the elongation region of the follicle (Orwin and Woods, 1982), and differences in the proportions of macrofibrils with increased double-twist intensity (related to the average helical tilt of the intermediate filaments) within and between cortical cells (Harland et al., 2011, 2014). The two factors are not mutually exclusive, and may also be underpinned by differential protein expression across the bulb (Thibaut et al., 2007) and associated differences in the macrofibril self-assembly processes (McKinnon and Harland, 2011).
Another theory of fibre curvature, which we term the cell proliferation model, explains fibre curvature as a product of differential cell division rates in follicles that results in the side of the fibre on the outside of the curvature having more cells than the inner side (Auber, 1951; Hynd et al., 2009). We found no evidence in the current study that supports this theory, although we recognise that our cell number data are derived from a calculation of cell numbers and not a direct count, and the precision of the data is less than if we were able to count the cells directly. However, there were never more cells on the outside of curve than the inside because the ratio of cells outside:inside was one or less in every fibre (Fig. 8).
Our study measured the length of cortical cells in situ, and the data indicate an average cell length of 73.0 µm for paracortical cells and 89.4 µm for orthocortical cells (Fig. 6). Cortical cells have traditionally been assumed to be of slightly variable dimensions but approximately 100 μm long and 5 μm wide irrespective of cell type (Hynd, 1994; Orwin, 1979). Our results complement earlier TEM studies of merino wool cross-sections, which concluded that orthocortical cells are larger than paracortical cells (Bones and Sikorski, 1967; Kassenbeck and Leveau, 1957; Orwin et al., 1984) based on TEM measurements of the maximum width of cortical cells. Our own unpublished TEM study on 21 µm merino fibres (Harland et al., 2006) also found that orthocortical cells had a significantly higher diameter (mean minimum=3.16 µm, mean maximum=5.89 µm; n=605) than did paracortical cells (mean minimum=2.46 µm, mean maximum=4.56 µm; n=557). The confocal data presented here add to this conclusion by indicating that, irrespective of cell type, cortical cells that are wider are also longer.
The few studies to investigate merino fibre cortical cell length have relied on measuring cells that have been enzymatically or chemo/mechanically isolated, usually dried and typically measured using scanning electron microscopy. Historic measurements include (length×maximum width): 95×5.5 μm (Lockhart, 1960); 102×3.6 μm (Kulkarni et al., 1971); and 109×5.2 µm from our own unpublished studies on 300 merino fibres with average diameter of 21 μm (Harland et al., 2006). This raises the question of why cells measured in this study were considerably shorter than in these earlier studies. One possible reason is that we were not measuring from the same points as for isolated cells because the packing of fibres in situ may constrain them. Also, our in situ cells were in fibres immersed in an anti-fade compound, which probably contains water, which may reduce the length of orthocortical cells due to swelling.
Conclusion: a revised model of intrinsic curvature in hairs
Our data led us to hypothesize that intrinsic fibre curvature at a specific location in merino wool fibres is a function of the relative difference in length between orthocortical and paracortical cells on the outside and inside, respectively, of the cortex along the plane of curvature. What we propose is a modified version of the cell type distribution model (Fig. 9). The relative difference is dynamic because cell length varies between fibres and possibly also along fibres. We hypothesize that intrinsic cortical cell length is controlled by processes in the follicle (including cell shaping), the inclusion of double-twist macrofibrils (i.e. similar to those in the orthocortex) within paracortical cells, or through variation in macrofibril double-twist intensity (see McKinnon and Harland, 2011). Intrinsic curvature is then modified by extrinsic factors such as dry set, disulfide exchange or the introduction of other crosslinks to result in a final curvature. This curvature contributes to the overall 3D shape of a fibre, which is a product of curvature and torsion. We believe that the inclusion of variation in cortical cell length as a parameter of the mechanism underpinning curvature should be universally applicable to mammalian keratin hairs.
We thank Manfred Ingerfeld and David Collings from the School of Biological Sciences and Biomolecular Interaction Centre, University of Canterbury, for assistance with confocal microscopy. We are also grateful for technical assistance and advice during the course of the research from Erin Lee, Ancy Thomas, Denis O'Connell and Chikako van Koten.
Conceptualization: D.P.H., J.L.W., S.N., T.I., K.K.; Methodology: D.P.H., J.A.V., J.L.W.; Validation: J.L.W.; Formal analysis: D.P.H.; Investigation: D.P.H., J.A.V., J.L.W.; Resources: D.A.S., A.J.G.; Data curation: J.A.V.; Writing - original draft: D.P.H., J.L.W.; Writing - review & editing: J.A.V., J.L.W., S.N., T.I., K.K., D.A.S., A.J.G., J.M.D., S.C.; Supervision: T.I., K.K., J.M.D., S.C.; Project administration: J.M.D., S.C.; Funding acquisition: S.N., T.I., K.K.
This work was supported by Kao Corporation and by AgResearch through its Integrated Wool Sciences Programme (New Zealand Government Ministry of Business Innovation & Employment, Science strategic investment fund).
Measurement data are available from figshare (https://doi.org/10.6084/m9.figshare.5500873.v1).
The authors declare no competing or financial interests.