Analytical Review of Integument – Research project update May 2022
Tiltwolf; Lathreas; Davies, D. (Athamanatha Kitsune); Roggenbuck, A. (Birdy); Cam Cam; Rakeela; Norsworthy, M. (Zennith)
Please refer to the glossary at the end of this document for an appendix with definitions of terms we use here.
Why are we working on this?
Our ultimate goal is to collect and/or design the required biochemical pathways for developing nonhuman integument types from human skin for the purpose of enabling freedom of form of the integument. There are a great number of methods of gene therapy already in existence (e.g. viral vectors, recombination-based methods, etc.), several of which having been used clinically (e.g. Luxturna™, voretigene neparvovec), meaning that gene delivery is not a major concern for the integument. The key challenge is identifying what genes we need to transfer into cells using these existing methods.
We are combining literature review and novel pathway engineering in our approach here, with the goal of collecting all the information we need in order to make first attempts at generating nonhuman integument from human cells for transspecies procedures. Based on this work, within the next few years we can begin in vitro experimentation using artificially cultured skin, and we will proceed to in vivo studies after experimental validation of methods we develop, along with extensive preclinical safety testing.
This update document is intended to give a glance at our progress on assembling our application-specific review of the integument. Compared to update documents 1 and 2 [link and link], this document is more focused on engineering rather than basic understanding of molecular and cellular biology. We have placed the majority of our emphasis for this update on scales, though we also have some interesting information on feathers and fur to share.
1. Scale development in reptiles and avians
Scales are fundamentally different from the existing types of integument within humans, and as such, there is very little published literature on their morphogenesis in terms of biochemical pathways involved. There are no published studies on scale growth patterning mechanisms to date that we are aware of, after substantial literature review. As such, there is insufficient information at present to create scales from human cells using natural pathways. Therefore, we will opt for a rationally-designed engineering approach to build a scale from components we can find.
In order to do this, we must first list the main properties we want to engineer as an end goal for us to work towards. Scales come in many different types and forms depending on the species. Snakes, for example, have directional and partially overlapping scales akin to roof shingles, whereas iguanas have bump-like non-overlapping scales. This is accompanied by a complex set of underlying epidermal and dermal layers constituting scale microarchitecture, and has been a subject of our review, but will not be discussed in detail here due to the type of system we are proposing to engineer. It will be discussed further in a future update.
As such, in our efforts, we should try to recreate the scale morphology depicted below as closely as possible. While there are also many constituent layers and cell populations in true reptilian scales, these are less important than the overall morphology, so we will not discuss those intermediate layers here.
Figure 1: Simple illustration of scale characteristics
These properties are quite varied, and ideally we would be able to tune all properties individually. The main properties to focus on in scales are the following, based on the scale morphology:
- Scale height (thickness)
- Scale size
- Scale asymmetry/shape
- Overlap between neighbouring scales
- Colouration
- Spatial distribution and variation
Each of these can be tuned by different genetic networks constituting patterns that direct morphogenesis, some of which may be artificially inserted into the stem cell pool of patient skin. However, we are currently limited by technical constraints due to lack of known patterning mechanisms in the literature as described, and limitations on the degree of genetic alteration able to be imposed at once, so we have focused on essential morphological characteristics: development of scale height and the hinge regions (valleys) between scales. We have, at this moment, assigned lower priority to patterning mechanisms for other features such as overlap, to fully focus on the essentials in the near term.
It is likely that these morphological features are governed by mechanisms akin to Turing patterns in nature, although they have not been described in detail in the literature to our knowledge. At minimum, we may assume that such mechanisms can produce said features, but whether they do would require experimental confirmation. For our purposes of engineering, the idea that patterning can produce scales in theory is sufficient for now, and underlies our further work described here.
1.1 Scale growth can be split into pattern-forming and morphogenesis parts
Natural systems often tightly couple pattern formation and morphogenesis, such as how the gradient of Wnt ligands simultaneously establishes a pattern and directs cell proliferation and other behaviours (e.g. Li et al., 2017). However, for engineering, we want more discrete systems.
In order to engineer a proper scale-forming mechanism, we need to separate this into tangible components that are independent of each other as compartmentalised “black boxes”. For our purposes, we will split the growth process up in two sequential steps: the initial pattern formation, and then morphogenesis. In this model, pattern formation would tell each small group of cells how fast to grow (determining scale thickness), in what direction to grow (determining scale directionality and overlap), whether it is in a border region where tissue needs to be more flexible, what colour it should have (in the case of melanocytes), etc. In such an engineered system, after each cell knows what to do, all cells can be given a collective ‘start signal’ to start growing, kickstarting morphogenesis. Using the information they have gained during pattern formation, the cells will start executing their growth and differentiation program. In essence, the morphogenesis step uses the pattern generated in the first step as its input, giving us two separated mechanisms we can study and engineer.
In practice, a growing tissue continuously cycles between pattern formation and morphogenesis, possibly even at the same time, pipelining the two processes. Even so, considering these as two separate mechanisms will allow us to engineer the system more easily.
1.2 Methods of pattern generation
We first need a method to tell specific cells what exactly they need to do, depending on their location. As such, we need to first pattern the skin such that each cell knows whether it is part of the thick part of a scale, whether it is part of the gap between scales, what the orientation is of the scale, how thick the scale should be, etc. Again, we aim to use simplified methods engineered on a de novo basis, as this will be substantially faster and likely yield a viable result within the coming decades. To do this, we propose two separate mechanisms to approach this problem with: a ‘paintbrush’ method, and a Turing pattern method.
1.2.1 The paintbrush method: easy to control/engineer, but less robust over time
The first simple mechanism approach we propose involves ‘painting on’ the pattern of scales manually, hence why we call it the ‘paintbrush method’. The ‘paint’ can be any signal, such as light (optogenetics, such as cryptochromes fused to transcription factors that are activated by a particular wavelength of light), a custom ligand for an inserted receptor, or something else (Hernandez-Candia & Tucker, 2020).
Naïve painting: population extinction causes pattern decay
Before we explain the full system, we will first explore a strawman system that doesn’t yet have the right properties. A naïve approach might be to use gene therapy to create two populations of epithelial stem cells: one population produces the thick parts of scales, and the other population produces the thin membrane that spans between neighbouring scales. One could imagine the latter being the default human skin thickness, with the thick parts of the scales constituting local skin thickening above and beyond what is normal for humans.
This approach sounds rather intuitive: simply inject each population of cells into hexagon shapes on a patient’s skin, akin to applying a tattoo, and hexagonal scales will form. The boundaries will be established by the lack of scale-forming epithelial stem cells, or even the injection of the special ‘boundary’ stem cells that only produce a thin membrane. This alternating pattern between scale-forming and boundary cells would then give rise to scales on skin.
However, this approach has several key problems, as we have hinted. The most important problem relates to pattern decay: without some way to maintain the boundaries between scales, over time, the boundaries will disappear, necessitating another solution which we have developed. In this first case, as epithelial stem cells divide, they can choose either to divide horizontally to maintain the stem cell population, or to divide vertically to differentiate into skin cells. Normal, random cell death or symmetric differentiation will remove a stem cell from the population. Overall, this has the effect that the amount of stem cells of a certain type will fluctuate. Smaller populations have a high chance to go extinct if the number of stem cells in the population ever fluctuates to zero. As such, the population of boundary cells, indicated in figure 2 as red cells, will over time disappear, causing scales to merge into each other over time.
Figure 2: Population extinction causes pattern decay
Some simple simulations have shown that decay occurs quite rapidly. The width of the scale boundaries wildly fluctuates over the span of weeks to months, and in only a couple of years, or even as little as a few months, entire scale boundaries could disappear causing scales to merge or deform.
To verify, we have performed a simple simulation of the stem cell populations in the epidermis. Each population is represented as N stem cells, each of which independently divides at each cell cycle. Each cell cycle is taken to be approximately 1.75 days. At each cell division step, 25% of the cells will divide into two differentiated cells, removing them from the stem cell population (Kuri & Rompolas, 2018). 50% of the cells will undergo asymmetric cell division, meaning one of the daughter cells will differentiate and one of the daughter cells will remain a stem cell. This effectively maintains the stem cell population, since a single stem cell will give a single stem cell after division. The remaining 25% of the cells will undergo symmetric cell division to produce two stem cells. Since each cell will do this maintenance at random, the population size will fluctuate over time. In some of the worst cases, this may lead to the complete extinction of a small population due to pure chance. We have plotted these cases in figure 3 and 4, for different population sizes. Although this is not guaranteed to occur for each population, this happens relatively often, and each scale boundary has a chance of decaying as a result.
Figure 3: A stochastic simulation of the amount of ‘scale boundary’ stem cells over time, starting with a population of 30 cells. Over time, the population size will fluctuate, until randomly populations may go extinct. y-axis: cell population size. x-axis: time in cell cycles. We assume one cell cycle is approximately 1.7 days (calculated from data in Kuri & Rompolas, 2018 and Oakley 2008).
Figure 4: A stochastic simulation of the amount of ‘scale boundary’ stem cells over time, starting with a population of 300 cells. Despite the much larger boundary size, the population still regularly goes extinct after a few years. y-axis: cell population size. x-axis: time in cell cycles. We assume one cell cycle is approximately 1.7 days.
This model does not take random cell death or cell location into account, which we expect to make the onset of pattern distortion or decay somewhat worse.
To prevent this problem, we need to ensure that the population of ‘scale-forming’ stem cells and ‘boundary-forming’ stem cells are maintained over long periods of time.
Although not perfect, due to the fact that pattern decay only becomes severe after a few months to years, one possible strategy would be to manually re-paint the pattern in the clinic periodically (say, annually). This presents a different set of technical challenges and risks, in that while patterning of naïve skin with an optogenetic light array is straightforward, doing the same for scaled skin requires extra steps. There are a number of methods we may use to address this, which would involve recalibration of the original instrumentation used. There are several techniques we are considering for this purpose, but for the sake of brevity, we will discuss them in a future update.
A more advanced and difficult-to-engineer approach is also being considered, where the pattern is maintained without requiring clinical intervention, by using engineered cell-to-cell communication. At this stage, such a method is hypothetical and would not be used in a “generation 1” system.
We will first discuss the ‘compromise’ strategy of clinical re-patterning, since it is technically the easiest to implement. This less-than-perfect, but viable strategy also provides the building blocks for the more advanced self-maintenance mechanism, allowing us to use it as an introduction for both methods.
“Smart paint” method: periodic clinical maintenance
For this method, the entire stem cell population in skin will be homogeneous, and contain a single circuit. Each cell listens to input, and switches to a ‘scale-forming’ state or ‘boundary’ state by flipping an internal bi-stable switch analogous to those that exist in nature (Verdugo et al., 2015). This way, the input signals can induce the right hexagonal pattern we need for scales, and when pattern decay occurs, we simply re-apply the pattern (e.g. an array of light) to reprogram cells that are at the wrong position. This mechanism, in a more advanced form, is actually present in many natural biological systems as well, as we will go into later.
Figure 5: the simplest possible Gene Regulatory Network (GRN) or Chemical Reaction Network (CRN) bi-stable switch. This switch forms a memory module for a cell. A and B are (macro)molecules that promote themselves, but repress each other. As such, when A is present, B is repressed and vice versa. Like a see-saw, this system is bi-stable, meaning it can switch back and forth by nudging the concentrations of A and B using the inputs. The concentrations of A and B can subsequently be read out by cellular circuits, such as by having A and B act as transcription factors, which forms the output of the system. For example, if A is present, this could mean that the cell is a scale-forming cell and must produce the thick layers of scales.
One way in which this can be done is depicted in figure 5. This system uses two inputs to switch a memory module called a bi-stable switch into a certain position. By applying input B, we can flip the cell into state ‘B’. By correctly programming the cell, this circuit could continuously listen to inputs to flip states, thereby allowing us to reprogram each cell that is mispatterned as often as necessary.
The ‘paint’, in this case, is the input signal of the circuit. It can be anything we choose, ranging from small diffusible morphogens to light waves that influence optogenetic receptors. We lean towards the use of optogenetic receptors due to their ease of programming and the ability to avoid using small molecules that themselves would require safety testing. All that being said, regardless of the ‘paint’ being used, we do need to insert this switching system into skin stem cells before they can switch based on morphogens or light waves. There are several ways to accomplish this using gene therapy, such as using a viral vector, virus-like nanoparticles, etc.
Constantly listening means constantly changing
One caveat that needs to be addressed is that if the circuit is constantly listening for input signals, especially if the signal is common (e.g. ambient light), the circuit might repattern scales in uncontrolled ways. As such, we need some mechanism to transiently enable this ‘listening’ behaviour, and then after clinical re-patterning is complete, remove the ability to change patterns again. We can do this by transiently expressing input receptors, such as through inducible expression systems activated by a small molecule (e.g. tetracycline, although many systems exist to select from for our range of use cases) that is administered before the re-patterning treatment begins. As an example, although not one we have committed to yet, we note that tetracycline is well-known in inducible expression systems and is used clinically as an antibiotic, where a one-time low dose for induction of expression is likely to be very well-tolerated.
Since any small molecule used for induction of expression would have a fairly short half-life, it will have cleared out of the system by the time (re-)patterning is complete, and the stem cell state becomes ‘fixed’. Now the patient will be able to expose their skin to many environments without risking accidental re-patterning.
The full circuit would therefore have to look something like what is presented in figure 6:
Figure 6: The proposed gene regulatory circuit that listens to and remembers external patterning information from the optogenetic or chemical ‘paintbrush’.
The full paintbrush system
The full system is summarised in figure 7. In step 1, the novel gene circuit is embedded into the epithelial stem cells of the skin using gene therapy. In step 2, the transient listening circuit is induced by a small-molecule inducer, enabling the cells to switch states. In step 3 & 4, the scale pattern is applied to the skin, causing the cells to switch states appropriately. In step 5, the stem cells have switched to be either a scale-forming or boundary stem cell, and produce the appropriate cell layers comprising scales.
The downstream effects of scale development requires a deep examination of scale-forming pathways in either reptiles or birds, which is an ongoing research effort.
Figure 7: two-state pattern induction
1.2.2 The Turing pattern method: likely harder to engineer, but much more robust
As a parallel method in development, we are attempting to engineer an artificial or semi-artificial Turing pattern which resembles the paintbrush method as described above. However, instead of requiring repeated application of an exogenous pattern for maintenance, it would be self-stabilising, as occurs in nature. The mechanism of scale morphogenesis by Turing patterning is not fully understood, so an artificial pattern could either make use of a natural one if it is discovered, or else be engineered orthogonally on a de novo basis. Currently, a small amount is known about natural patterning, in that it may share components with the first stage of feather patterning (Li et al., 2017).
Such an engineered mechanism would essentially have two steps. Step 1 inserts a genetic control circuit into the relevant stem cell population, which will be dormant until provoked by a manually administered global signal that starts the patterning process. Step 2 then involves cells communicating in a complex network to generate an alternating pattern that induces scale development.
Simple Turing patterns have recently been engineered in both prokaryotic and eukaryotic cell colonies (Duran-Nebreda et al., 2020, Sekine et al., 2018). Even so, bioengineering of controlled, robust or complex Turing patterns in mammalian cell populations remains a serious challenge. One of the major issues is orthogonality, where it is important that the patterning network components do not have off-target effects on other pathways in the cell. Generating orthogonal pathways can take two forms; either it can involve rational protein engineering of custom receptor-ligand pairs, or it can use non-protein methods such as nucleic acids, which we view as more tunable based on unpublished proprietary designs. Our future work will involve assessing methods for construction of patterning mechanisms that are tuneable, robust, and effective in vivo for generating structures that are otherwise unnatural in humans, but this is necessarily a longer term goal which we feel will likely come after simpler methods are established as described previously.
Figure 8: The proposed system of scale maintenance in vivo. Scale boundaries are maintained through long-range inhibition and short-range activation. In nature, this leads to so-called ‘Turing patterns’, which is the process in which hair follicles, scales, and feather follicles are distributed on the skin during embryogenesis. In this system, we would mimic this process, albeit at a larger length scale to maintain the patterns.
1.2.3 The real system may use a combination of both methods
In reality, we will likely use a combination of both methods. For example, Turing patterns are relatively robust, and will tend to produce patterns with the right size and properties. However, without a good nudge, the formed pattern may not look exactly the way we want it to look. We can supply the Turing mechanism with an initial pattern using the paintbrush method to help it along, which in mathematical terms is known as an ‘initial condition’. This has the advantage of giving us a lot of control over the final position and size of scales, while also ensuring that the pattern remains stable over long periods of time.
1.2.4 The genetic network complexity is limited by gene vector size
Regardless of what we design, we are limited in the complexity of artificial gene networks that we can package into clinically viable viral vectors. This is because for most vectors used in gene therapy, 10 kb (kilobases) is the maximum size, with a hard limit of about 12 kb, including regulatory sequence elements. In practice, about 8 kb can be used for transgenes, and we would aim to include the necessary genetic modification in as few vectors as possible due to concerns about transduction efficiency. To achieve this, we may, for example, package multiple genes in the same multicistronic vector, separated by internal ribosome entry sites (IRES).
Very small proteins take as little as 0.6 kb, but some are large – for example, a receptor for Shh, called PTCH1 (Patched) in humans, is 4.3 kb. An average gene is around 1.5 kb, so on average, we could expect to include 5-6 genes, making for a very minimalistic network.
There are more exotic methods for gene transduction that can insert larger sequences, such as hydrodynamic gene therapy (HGD), sonoporation, lipid nanoparticles (LNPs), etc., but they are not currently in clinical use, so research and development would take longer. It is therefore important for us to create a barebones gene network for our purposes that can work with currently clinically usable minimalist vectors. We will keep the option of larger gene networks with the aforementioned future methods open as a possibility for use at a later date.
1.3 Methods of morphogenesis
1.3.1 Top-down approach: 3D bioprinting
The field of bioprinting and engineered tissues has introduced many new possibilities in the way we approach both development in vitro as well as integration in vivo of artificial integument (Weng et al, 2021, Kolesky et al, 2014). Artificial constructs may be constituted by direct application of suspended cells in the form of “bioinks” and/or fabrication of biocompatible scaffolds to be infused with cells in a bioreactor afterwards. Importantly, this approach sidesteps patterning by instead manually constructing the desired morphology.
Bioprinting and biofabrication serve two important purposes in our future research. The first would be to generate scaled skin for use in grafting, although applying grafts to wide areas of skin is very invasive and thus not ideal. More interestingly, the second purpose is to generate artificial skin for testing of the paintbrush method in vitro long before patient use.
1.3.2 Bottom-up approach: Cellular Growth & Asymmetric Cell Division
When scales grow, they grow outwards, in the so-called ‘apical’ direction. The amount of growth is probably modulated using an underlying pattern which controls the proliferation rate, although we are not aware of a known patterning mechanism for this yet. More complex patterns that are rationally designed may be able to introduce more complex scale shapes. However, for overlapping scales, which many forms require, just growing upwards is not enough to cause overlap, no matter how complex the growth pattern. As such, for overlapping scales, the growth direction must also be taken into account. There are two ways in which this directionality might be induced in nature. A generally accepted common method, present in many species, requires interfacing with the polarity axes of the epidermal stem cells. A second method might involve asymmetric patterning, causing the scales to fold and overlap like origami. These two mechanisms might even be intrinsically linked.
In either case, some way to control the direction of growth, whether perpendicular or parallel to the skin, or at some controlled angle, will be necessary to induce overlap in scales through simple growth or complex folding. As such, we will briefly examine the polarity axes within skin cells.
Figure 9: Overlapping scales with simple patterns.
1.3.3 Cell polarity axes in skin
There are two polarity axes of relevance present in these cells: the apical-basal (up-down) polarity axis and the planar (side-to-side) polarity axis. Normally, skin growth is aligned with the apical-basal polarity axis when the epithelial stem cells differentiate. During stem cell maintenance, cell division is aligned along the planar polarity axis (Muroyama & Lechler, 2013). These two processes must remain in balance for the skin to properly maintain itself.
It is possible that scale growth may be oriented at an angle in scales such that it does not grow perpendicular to the skin. Alternatively, there may be a balance between apical-basal and planar growth that yields a net angular direction of growth. The orientation of polarity, whether it be along the scale’s normals or aligned with an overall angled axis, should be assayed in reptilian skin by immunohistochemistry.
In hair follicles, the growth direction is dependent on the basal polarity axis (Chen & Chuong, 2013). However, simply defining the angle with respect to the vector pointing out of the skin is not enough, because the orientation that the hair/scale/feather grows in – the direction of the ‘grain’ of the hair, i.e. the direction in which you should pet a cat – must also be known. That is likely defined by the planar cell polarity of the actively growing cells, i.e. those differentiating from the epidermal cell pool near the base.
For overlapping scale development (snake scales, fish scales, etc), it is important that the growth direction is modulated with respect to the planar cell polarity axis as well as the apical-basal cell polarity axis. Cell division is oriented along the polarity axis by aligning the mitotic spindle with that axis, and this, in turn, is controlled by forces applied to the spindle through microtubules and accessory plasma membrane proteins. It thus becomes important to control orientation of polarity, as this is directly upstream of controlling cell growth direction.
1.4 Section Conclusion
Together, these mechanisms provide a clear research direction for creating scales from scratch. Although these theoretical systems are expected to work, our remaining challenge will be to identify physical targets that abide by the design’s parameters, and this forms the content of some of our current work.
2. Fur and hair development in mammals
As discussed in previous updates, to produce realistic fur, we must take into consideration a number of morphological features. These include:
- Total follicle density (hair follicles per cm2), required for sufficient coverage and coat thickness
- Topcoat/undercoat ratio and morphological differences between the constituent follicles
- Fur length regulation through hair cycle timing (i.e. anagen phase length)
- Hair distribution, kink properties, and curvature
- Hair roughness, texture, and eccentricity (i.e. degree of oval-ness in a cross-section)
2.1 Increasing follicle density for complete fur coverage
We previously covered some of the basic mechanisms for how fur follicle density is naturally controlled. These include major signalling pathways such as Wnt signalling, along with evolutionary differences such as the loss of fur-specific KRT41 in humans. At this stage, we are beginning to design concrete ways to modify follicle density in adult humans.
We have become very interested in methods of culturing whole hair follicles in vitro from co-culture of epidermal and dermal stem cells, as has been previously described (Žnidarič et al., 2021). These can be produced either in skin explants, artificially produced skin, or more interestingly, artificial extracellular matrix-like substances (e.g. Matrigel™). This method offers the flexibility of first modifying the co-cultured cells in ways that alter the morphogenesis of hair, as gene editing is drastically simpler when done in vitro.
As an alternative, simpler approach, with the tradeoff of being harder to modify morphogenesis, it may be possible to take existing hair follicles as explants and divide them into multiple follicles by longitudinal microdissection. To our knowledge, microdissection of biopsied hair follicles followed by re-implantation has not been attempted before, and whether follicles can heal from severe damage (shearing) is not well-established in this context. This would be a relatively simple experiment to conduct in a simple tissue culture lab, although since it involves humans (as laboratory mouse hair follicles differ too much to yield interpretable results), it would need to be approved by a research ethics board first; we may pursue it in the coming years.
2.2 Modifying fur physical properties
We are also further investigating morphological control of fur thickness by genetic means, with the best method appearing to be modulation of the Wnt10b/Dkk1 axis, where Wnt10b upregulates size (Lei et al., 2014). In that study, modulation was done through adenoviral gene therapy in vivo, but this approach is somewhat unrefined and does not target specific cell populations, instead causing transgene overexpression in all cells within a defined radius of a hypodermic injection of the vector. We aim to target the hair follicle dermal papilla directly, such as by using a pseudogenized viral vector or similar method, or by putting transgenes on cell type-specific promoter sequences.
However, it also makes sense to experimentally validate this approach in humans in a less permanent manner first. Wnt signalling has previously been modulated using artificial peptides, such as one homologous to Wnt3a that operates as a competitive inhibitor of Frizzled – this peptide is called UM206 and carries the sequence Ac-CNKTSEGMDGCEL-NH2 (Fu et al., 2019). It is possible to use the same approach to inhibit Dickkopf proteins (Dkk1) (Gregory et al., 2005). We have defined the corresponding peptide for Wnt10b, of a proprietary sequence which has about 60% sequence identity to UM206. Upon ethics approval, our intention is to apply it locally to patches of skin to determine if the peptide causes a change in hair thickness, which would serve as verification of the approach.
We are also in the process of constructing a table of morphological features of fur across different species, with corresponding approaches for modifying human cell-derived hairs to match these features. Some data for this table has proven difficult to find and may require primary research; we intend to research this further in future.
3 Feather production in avians
Feathers are thought to have first developed in pre-dinosaurs, having several complex features which developed over time (Di-Poi & Milinkovitch, 2016). The avian species of today take advantage of bringing those features together, although some birds such as ostriches lack the full set of feather types and are unable to fly (Xu & Guo, 2009; Foth, 2011; Yang et al., 2018 and 2020; Unwin & Martill, 2020).
We have spent a significant amount of time compiling and understanding the extensive and scattered research on different genetic aspects of feathers, with there likely being yet more out there to find. Chickens are most heavily studied, with significant industrial interest in breeding and Mendelian genetics of chicken phenotypes. The wealth of scientific papers and articles out there is testament to the effect such interest has when compared to the relatively few papers on scales.
While we have previously described a nearly complete patterning mechanism in our last update based on existing research, many other determinants of feather morphology exist (Li et al., 2017). As we learn more about these accessory factors, we may become able to incorporate them into a gene network that can operate within human skin.
As we find relevant pathways, we have generated interactome diagrams as shown. We are also marking genes based on sequence identity (100% being identical) between chicken and human versions, to predict which ones require chicken versions to be transposed into the human genome. For this, we have used NCBI BLAST and COBALT.
Our list of data on feathers continues to grow from the previous update and onwards into the future, and also contains much more detail and functional annotation. Interestingly, a significant proportion of the research papers we have found and reviewed since the last update handle feather pigmentation, including genes applicable not only to feathers but to all animal integuments.
3.1 Feather gene homologies based on rational literature review
We have analysed the genes we found in our literature search for their similarities to human equivalent genes, with a view to understanding what the necessary changes would be to enable feather growth from a human baseline.
Figure 10: histogram of sequence identities in feather-related genes
Here we have a histogram of the distribution of sequence identities (percentages of the protein residues which match each other) between chicken genes for feather development and their human equivalents (Figure 10). From this, we can see a statistical distribution (perhaps chi-squared or beta) skewed towards the 80-90% bin. It should be noted that the methods used to calculate these sequence homologies may not be optimal, and in our final version of the review, we will compare multiple methods and seek to identify which is the most suitable and consistent method for such comparisons. However, this histogram and bar chart can be viewed as offering ‘ballpark figures’ in this update.
The bar chart below (Figure 11) is the list of all the genes we have found so far related to feathers (and a few other items), arranged by sequence identity, with a significant portion we have as yet not been able to compare. At sequence identities below a certain percentage, it may be more prudent to replace the human version of a gene with that of an avian species, than to try to copy and/or modify an existing one.
The conservation of these genes between chickens and humans will at least partly inform us of the degree to which these genes are likely to affect the differences between human and avian integument. More diverged genes may actively encode components necessary for the specific integument type, while more conserved genes may be regulators and transcription factors worth adjusting in small ways to effect larger changes in expression (for example, the CTCF pathway, see Figure 12). Other strongly conserved genes (potentially acting as “master regulators”) with roles in chromatin remodelling are revealed in the literature we reviewed, such as BMI1/Polycomb (Hernández-Muñoz et al, 2005), SATB1 (Yasui et al, 2002; Fessing et al, 2011; Liang et al, 2020; Ramanujam et al, 2021) and SATB2 (Liang et al, 2020).
Figure 11: Chicken-to-Human Homology for Currently Analysed Feather-Related Genes
3.2 Known feather follicle patterning mechanism
This refers to how feather follicles are laid out and what the mechanisms are that control the distribution and positions of cells and signalling molecules for feathers to grow to be the shapes we observe at the densities we observe.
CTCF mediated chromatin folding is a transcription adjustment pathway which makes genes more or less available for transcription by unwinding them from histones, and presenting them in loops to the transcription complex. It has a major function in the availability of feather and scale keratin genes (Liang et al, 2020, see Figure 12).
Figure 12: CTCF signalling mediates chromatin folding, including in Chicken Chromosomes 25 and 27 to unravel them from their histones and create loops to expose specific genes to transcription. (Liang et al, 2020).
In our Research Update #2, we discussed the putative Turing pattern proposed by Li et al., which describes all but one of the components of the gene network that leads to feather formation in avians (Li et al., 2017). These included elements of BMP signalling (BMP2, GREM1), Wnt signalling (WNT3A, possibly others), RA signalling (CRABP1, CYP26B1), GDF10, and a mystery activator of BMP2.
We have been considering the role of a mystery activator of feather follicle morphogenesis. Our research has found that the literature is not 100% clear on this activation agent or set of agents, but its properties necessarily must include:
- Activates BMP2 (or inhibits it)
- Is inhibited by BMP2 (or activates it)
- Must activate itself
Potentially, this could be or include SHH (Sonic Hedgehog), which is known to potentiate BMP signalling, is inhibited by BMP2 (Rios et al., 2004), and controls expression of GREM1. SHH has the best characteristics, although there is no evidence that SHH overexpression does not perturb feather morphogenesis, and this was taken as evidence against it being the mystery activator (Li et al., 2017). We aim to identify the mystery activator by omics-level studies of developing feathers in future studies.
In the coming months, we will be focussed on determining the best way to recapitulate this patterning mechanism in humans, whether it be to do so in vitro followed by follicle transplantation, or direct modification of the skin. We have identified several missing genes that must be introduced into human cells, including beta-keratins and their accessory proteins, and several more with poor homology where an avian copy may be superior. Once the entire proposed pathway is laid out, we will be in a good position to see if we can recapitulate feathers from human cells by applying the pathway to cultured hair follicle cells for attempted reprogramming to produce feathers instead.
Where do we go from here?
Within this update, we have extensively focussed on the rational engineering of scale morphogenesis in adult human skin in concept, as well as presented some progress in our other project sections on fur and feathers. With our specific focus on scales, this gives us an extensive set of research directions for taking our proposals from concept to practice. Our next steps should aim to integrate the theoretical designs into practical targets that can be tested in the lab.
There are remaining pieces of information we need to collect in order to develop full sets of hypotheses for inducing the development of our three integument types of interest in humans. Our intention is to finalise that information in the coming months, with the goal of publishing a complete research paper that can then lead to several spinoff projects involving further research and application.
Beside just feathers, scales, or fur, there are of course many more types of integument that are not currently within our scope to review or design. This includes chitin on insects, cetacean skin adapted to aquatic environments, etc. To avoid scope creep, we will only briefly touch upon these integumentary types, and leave a detailed review for a later project.
Accessory integument types that will be reviewed in a subsequent project include dog noses, horns, claws, etc.
We are excited by the progress we have made over the past year! If we are able to obtain sufficient funding, we should be ready to proceed to wet lab experimentation in the very near future. Our team welcomes discussions with future donors for substantial directed donations to fund this next stage in research to take concept to practise.
PS: As always, if there are STEM and/or business professionals out there who read this and are interested in working with us, please do not hesitate to reach out!
Sincerely,
The Project Team
Analytical Review of Integument
Freedom of Form Foundation
501(c)(3) EIN 82-4415111
www.freedomofform.org
Acknowledgements
NG, L, DD, and CC conceptualised the study. L, NG, and AR designed and performed simulations for the paintbrush method for scale development. NG designed Wnt10b-based methods for hair follicle regulation. DD, L, and NG performed analysis of feather developmental pathways. NG, AR, DD, and CC performed comparative gene sequence analysis for multiple sections. NG, L, DD, and MN produced figures. All authors contributed to literature searching throughout. All authors also contributed to writing the manuscript. NG and MN performed final editing and review prior to publication.
We would also like to thank Keiran Stevenson (Syralth) for discussions on intellectual property disclosure in this project.
Glossary
Term | Definition |
Turing pattern | Biomolecular distribution constituting a reaction diffusion mechanism governing development of anatomical structures |
Assay | Laboratory analytical test for the composition or characteristics of a substance, object, or other element |
Protein | Executor molecule in cells composed of a chain of amino acids, usually referred to by the name of the gene that codes for it |
Peptide | Very short protein with only a limited number of amino acids |
Homology | The degree of similarity between proteins or between genes, where ones that are similar are said to be homologous |
BLAST | Basic local alignment search tool for determining significance of homology between two gene or protein sequences |
COBALT | Constraint-based alignment tool for finding the best alignment between two sequences previously determined to have significant homology |
Viral vector | Engineered virus for inserting new genes into cells |
Bistable switch | Gene/protein circuit which has two equilibrium states of different functions that are stable, but which can be forced to switch by an inducer |
Inducible expression | Drug-controlled on/off state of whether a new gene is expressed (in use) at any given time |
Explants | Sections of tissue removed from a body |
Microdissection | Physically cutting small structures into yet-smaller structures along precise boundaries for analysis and engineering |
Gene network | A network defined by chain reactions, activations, up- and down-regulation, etc, between genes, usually via their protein products or via epigenetics. |
Upregulation | Intensification of a function, such as by reading a particular gene more often or otherwise rendering its product more active |
Downregulation | Lessening of a function; opposite to the above |
Chromatin | The normal state of DNA in a non-dividing cell’s nucleus |
Integument | The protective or decorative layers of the external surfaces of an organism’s body, such as hairs, scales, feathers or an exoskeleton |
Transcription Factor | Protein that binds to DNA and helps to promote specific genes to be upregulated by being read more often |
Multicistronic vector | Having multiple genes inside of a single gene vector for transfer into cells all at once |
IRES | Internal ribosome entry site, a sequence separating two genes destined to be transcribed onto the same RNA strand while being synthesized as separate proteins |
Wnt Signalling | A cascade of biochemical reactions that leads to specific outcomes, starting from or involving the gene family we know as ‘Wnt’ (Wingless and Int-1) |