How the Leopard Gets Its Spots, Lecture notes of Geometry

He then employed a mathematical model to show that if morphogens react and diffuse in an appropriate way, spatial patterns of morphogen con- centrations can ...

Typology: Lecture notes

2022/2023

Uploaded on 02/28/2023

hollyb
hollyb 🇺🇸

4.8

(44)

431 documents

1 / 8

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
How the Leopard Gets Its Spots
A single pattern-formation mechanism could underlie the wide
variety of animal coat markings found in nature. Results from
the mathematical model open lines of inquiry for the biologist
by James D. Murray
ammals exhibit a remarkable
variety of coat patterns; the
variety has elicited a compa-
rable variety of explanations—many
of them at the level of cogency that
prevails in Rudyard Kipling's delight-
ful "How the Leopard Got Its Spots."
Although genes control the proces-
ses involved in coat pattern forma-
tion, the actual mechanisms that cre-
ate the patterns are still not known. It
would be attractive from the view-
point of both evolutionary and devel-
opmental biology if a single mech-
anism were found to produce the
enormous assortment of coat pat-
terns found in nature.
I should like to suggest that a single
pattern-formation mechanism could
in fact be responsible for most if not
all of the observed coat markings. In
this article I shall briefly describe a
simple mathematical model for how
these patterns may be generated in
the course of embryonic develop-
ment. An important feature of the
model is that the patterns it generates
bear a striking resemblance to the
patterns found on a wide variety of
animals such as the leopard, the
cheetah, the jaguar, the zebra and
the giraffe. The simple model is also
consistent with the observation that
although the distribution of spots on
members of the cat family and of
stripes on zebras varies widely and
is unique to an individual, each kind
of distribution adheres to a general
theme. Moreover, the model also pre-
dicts that the patterns can take only
certain forms, which in turn implies
the existence of developmental con-
straints and begins to suggest how
coat patterns may have evolved.
It is not clear as to precisely what
happens during embryonic develop-
ment to cause the patterns. There are
now several possible mechanisms
that are capable of generating such
patterns. The appeal of the simple
80
model comes from its mathematical
richness and its astonishing ability
to create patterns that correspond to
what is seen. I hope the model will
stimulate experimenters to pose rele-
vant questions that ultimately will
help to unravel the nature of the bio-
logical mechanism itself.
ome facts, of course, are known
about coat patterns. Physically,
spots correspond to regions of differ-
ently colored hair. Hair color is deter-
mined by specialized pigment cells
called melanocytes, which are found
in the basal, or innermost, layer of
the epidermis. The melanocytes gen-
erate a pigment called melanin that
then passes into the hair. In mam-
mals there are essentially only two
kinds of melanin: eumelanin, from
the Greek words eu (good) and
melas (black), which results in black
or brown hairs, and phaeomelanin,
from phaeos (dusty), which makes
hairs yellow or reddish orange.
It is believed that whether or not
melanocytes produce melanin de-
pends on the presence or absence of
chemical activators and inhibitors.
Although it is not yet known what
those chemicals are, each observed
coat pattern is thought to reflect
an underlying chemical prepattern.
The prepattern, if it exists, should re-
side somewhere in or just under the
epidermis. The melanocytes are
thought to have the role of "reading
out" the pattern. The model I shall
describe could generate such a
prepattern.
My work is based on a model de-
veloped by Alan M. Turing (the in-
ventor of the Turing machine and the
founder of modern computing sci-
ence). In 1952, in one of the most im-
portant papers in theoretical biology,
Turing postulated a chemical mecha-
nism for generating coat patterns. He
suggested that biological form fol-
lows a prepattern in the concentra-
tion of chemicals he called
morpho-gens. The existence of
morphogens is still largely
speculative, except for
circumstantial evidence, but Turing's
model remains attractive because it
appears to explain a large number of
experimental results with one or two
simple ideas.
Turing began with the assumption
that morphogens can react with one
another and diffuse through cells. He
then employed a mathematical model
to show that if morphogens react and
diffuse in an appropriate way, spatial
patterns of morphogen con-
centrations can arise from an initial
uniform distribution in an assem-
blage of cells. Turing's model has
spawned an entire class of models
that are now referred to as
reaction-diffusion models. These
models are applicable if the scale of
the pattern is large compared with the
diameter of an individual cell. The
models are applicable to the leopard's
coat, for instance, because the
number of cells in a leopard spot at
the time the pattern is laid down is
probably on the order of 100.
Turing's initial work has been
developed by a number of investi-
gators, including me, into a more
complete mathematical theory. In a
typical reaction-diffusion model one
starts with two morphogens that can
react with each other and diffuse at
varying rates. In the absence of dif-
fusion—in a well-stirred reaction,
for example—the two morphogens
would react and reach a steady uni-
form state. If the morphogens are
now allowed to diffuse at equal rates,
any spatial variation from that steady
state will be smoothed out. If, however,
the diffusion rates are not equal,
LEOPARD reposes. Do mathematical as
well as genetic rules produce its spots?
M
S
pf3
pf4
pf5
pf8

Partial preview of the text

Download How the Leopard Gets Its Spots and more Lecture notes Geometry in PDF only on Docsity!

How the Leopard Gets Its Spots

A single pattern-formation mechanism could underlie the wide

variety of animal coat markings found in nature. Results from

the mathematical model open lines of inquiry for the biologist

by James D. Murray

ammals exhibit a remarkable variety of coat patterns; the variety has elicited a compa- rable variety of explanations—many of them at the level of cogency that prevails in Rudyard Kipling's delight- ful "How the Leopard Got Its Spots." Although genes control the proces- ses involved in coat pattern forma- tion, the actual mechanisms that cre- ate the patterns are still not known. It would be attractive from the view- point of both evolutionary and devel- opmental biology if a single mech- anism were found to produce the enormous assortment of coat pat- terns found in nature. I should like to suggest that a single pattern-formation mechanism could in fact be responsible for most if not all of the observed coat markings. In this article I shall briefly describe a simple mathematical model for how these patterns may be generated in the course of embryonic develop- ment. An important feature of the model is that the patterns it generates bear a striking resemblance to the patterns found on a wide variety of animals such as the leopard, the cheetah, the jaguar, the zebra and the giraffe. The simple model is also consistent with the observation that although the distribution of spots on members of the cat family and of stripes on zebras varies widely and is unique to an individual, each kind of distribution adheres to a general theme. Moreover, the model also pre- dicts that the patterns can take only certain forms, which in turn implies the existence of developmental con- straints and begins to suggest how coat patterns may have evolved. It is not clear as to precisely what happens during embryonic develop- ment to cause the patterns. There are now several possible mechanisms that are capable of generating such patterns. The appeal of the simple model comes from its mathematical richness and its astonishing ability to create patterns that correspond to what is seen. I hope the model will stimulate experimenters to pose rele- vant questions that ultimately will help to unravel the nature of the bio- logical mechanism itself. ome facts, of course, are known about coat patterns. Physically, spots correspond to regions of differ- ently colored hair. Hair color is deter- mined by specialized pigment cells called melanocytes, which are found in the basal, or innermost, layer of the epidermis. The melanocytes gen- erate a pigment called melanin that then passes into the hair. In mam- mals there are essentially only two kinds of melanin: eumelanin, from the Greek words eu (good) and melas (black), which results in black or brown hairs, and phaeomelanin, from phaeos (dusty), which makes hairs yellow or reddish orange. It is believed that whether or not melanocytes produce melanin de- pends on the presence or absence of chemical activators and inhibitors. Although it is not yet known what those chemicals are, each observed coat pattern is thought to reflect an underlying chemical prepattern. The prepattern, if it exists, should re- side somewhere in or just under the epidermis. The melanocytes are thought to have the role of "reading out" the pattern. The model I shall describe could generate such a prepattern. My work is based on a model de- veloped by Alan M. Turing (the in- ventor of the Turing machine and the founder of modern computing sci- ence). In 1952, in one of the most im- portant papers in theoretical biology, Turing postulated a chemical mecha- nism for generating coat patterns. He suggested that biological form fol- lows a prepattern in the concentra- tion of chemicals he called morpho-gens. The existence of morphogens is still largely speculative, except for circumstantial evidence, but Turing's model remains attractive because it appears to explain a large number of experimental results with one or two simple ideas. Turing began with the assumption that morphogens can react with one another and diffuse through cells. He then employed a mathematical model to show that if morphogens react and diffuse in an appropriate way, spatial patterns of morphogen con- centrations can arise from an initial uniform distribution in an assem- blage of cells. Turing's model has spawned an entire class of models that are now referred to as reaction-diffusion models. These models are applicable if the scale of the pattern is large compared with the diameter of an individual cell. The models are applicable to the leopard's coat, for instance, because the number of cells in a leopard spot at the time the pattern is laid down is probably on the order of 100. Turing's initial work has been developed by a number of investi- gators, including me, into a more complete mathematical theory. In a typical reaction-diffusion model one starts with two morphogens that can react with each other and diffuse at varying rates. In the absence of dif- fusion—in a well-stirred reaction, for example—the two morphogens would react and reach a steady uni- form state. If the morphogens are now allowed to diffuse at equal rates, any spatial variation from that steady state will be smoothed out. If, however, the diffusion rates are not equal, LEOPARD reposes. Do mathematical as well as genetic rules produce its spots?

M

S

ZEBRA STRIPES at the junction of the foreleg and body (left) can be produced by a reaction-diffusion mechanism (above). sion models concerns the outcome of beginning with a uniform steady state and holding all the parameters fixed except one, which is varied. To be specific, suppose the scale of the tissue is increased. Then eventually a critical point called a bifurcation value is reached at which the uni- form steady state of the morphogens becomes unstable and spatial pat- terns begin to grow. The most visually dramatic exam- ple of reaction-diffusion pattern for- mation is the colorful class of chemi- cal reactions discovered by the Soviet investigators B. P. Belousov and A. M. Zhabotinsky in the late 1950's [see "Rotating Chemical Reactions," by Arthur T. Winfree; S CIENTIFIC AMERICAN, June, 1974]. The reactions visibly organize themselves in space and time, for example as spiral waves. Such reactions can oscillate with clocklike precision, changing from, say, blue to orange and back to blue again twice a minute. Another example of reaction-dif- fusion patterns in nature was dis- covered and studied by the French chemist Daniel Thomas in 1975. The patterns are produced during reac- tions between uric acid and oxygen on a thin membrane within which the chemicals can diffuse. Although the membrane contains an immobil- ized enzyme that catalyzes the reac- tion, the empirical model for describ- ing the mechanism involves only the two chemicals and ignores the en- | zyme. In addition, since the mem- brane is thin, one can assume cor- rectly that the mechanism takes I place in a two-dimensional space. I should like to suggest that a good candidate for the universal mecha- nism that generates the prepattern for mammalian coat patterns is a re- action-diffusion system that exhibits diffusion-driven spatial patterns. Such patterns depend strongly on the geometry and scale of the do- main where the chemical reaction takes place. Consequently the size and shape of the embryo at the time the reactions are activated should determine the ensuing spatial pat- terns. (Later growth may distort the initial pattern.) ny reaction-diffusion mechanism Lcapable of generating diffusion-driven spatial patterns would provide a plausible model for animal coat markings. The numerical and mathematical results I present in this article are based on the model that grew out of Thomas' work. Employing typical values for the parameters, the time to form coat patterns during embryogenesis would be on the order of a day or so. Interestingly, the mathematical problem of describing the initial stag- es of spatial pattern formation by re- action-diffusion mechanisms (when departures from uniformity are mi- nute) is similar to the mathematical problem of describing the vibration of thin plates or drum surfaces. The ways in which pattern growth de- pends on geometry and scale can therefore be seen by considering analogous vibrating drum surfaces. If a surface is very small, it simply will not sustain vibrations; the distur- bances die out quickly. A minimum size is therefore needed to drive any sustainable vibration. Suppose the drum surface, which corresponds to the reaction-diffusion domain, is a rectangle. As the size of the rectangle is increased, a set of increasingly complicated modes of possible vi- bration emerge. An important example of how the geometry constrains the possible modes of vibration is found when the domain is so narrow that only simple—essentially one-dimension- al—modes can exist. Genuine two-dimensional patterns require the domain to have enough breadth as well as length. The analogous requirement for vibrations on the surface of a cylinder is that the radius cannot be too small, otherwise only quasi-one-dimensional modes can exist; only ringlike patterns can form, in other words. If the radius is large enough, however, two-dimensional patterns can exist on the surface. As a consequence, a tapering cylinder can exhibit a gradation from a two-di- mensional pattern to simple stripes [see illustration on opposite page]. Returning to the actual two-mor-phogen reaction-diffusion mechanism I considered, I chose a set of reaction and diffusion parameters that could produce a diffusion-driven instability and kept them fixed for all the calculations. I varied only the scale and geometry of the domain. As initial conditions for my calculations, which I did on a computer, I chose random perturbations about the uni- form steady state. The resulting pat- terns are colored dark and light in re- gions where the concentration of one of the morphogens is greater than or less than the concentration in the ho- mogeneous steady state. Even with such limitations on the parameters and the initial conditions the wealth of possible patterns is remarkable. A

EXAMPLES OF DRAMATIC PATTERNS occurring naturally are found in the anteater (left) and the Valais goat, Capra aegagrus hircus (right). Such patterns can be accounted for by the author's reaction-diffusion mechanism (see bottom illustration on these

How do the results of the model cheetah (Acinonyx jubatus), the jag- relatively short, and so one would

compare with typical coat markings uar (Panthem onca) and the genet expect that it could support spots to

and general features found on ani- {Genetta genetta) provide good exam- the very tip. (The adult leopard tail is

mals? I started by employing taper- pies of such pattern behavior. The long but has the same number of ver-

ing cylinders to model the patterns spots of the leopard reach almost to tebrae.)The tail of the genet embryo,

on the tails and legs of animals. The the tip of the tail. The tails of the at the other extreme, has a remark-

results are mimicked by the results cheetah and the jaguar have distinct- ably uniform diameter that is quite

from the vibrating-plate analogue, ly striped parts, and the genet has thin. The genet tail should therefore

namely, if a two-dimensional region a totally striped tail. These obser- not be able to support spots,

marked by spots is made sufficiently vations are consistent with what is The model also provides an in-

thin, the spots will eventually change known about the embryonic struc- stance of a developmental con-

to stripes, ture of the four animals. The prenatal straint, documented examples of

The leopard (Panthera pardus), the leopard tail is sharply tapered and which are exceedingly rare. If the

SCALE AFFECTS PATTERNS generated within the constraints of a generic animal shape in the author's model. Increasing the scale and holding all other parameters fixed produces a remark- able variety of patterns. The model agrees with the fact that

on the initial conditions, the geome- try and the scale. An important aspect of the mechanism is that, for a given geometry and scale, the patterns generated for a variety of random initial conditions are qualitatively similar. In the case of a spotted pattern, for example, only the distribution of spots varies. The finding is consistent with the individuality of an animal's markings within a species. Such individuality allows for kin recognition and also for general group recognition. The patterns generated by the model mechanism are thought to correspond to spatial patterns of morphogen concentrations. If the concentration is high enough, mela-nocytes will produce the melanin pigments. For simplicity we assumed that the uniform steady state is the threshold concentration, and we rea- soned that melanin will be generated if the value is equal to or greater than that concentration. The assumption is somewhat arbitrary, however. It is reasonable to expect that the thresh- old concentration may vary, even within species. To investigate such DIFFERENT GIRAFFES have different kinds of markings. The subspecies Gi-raffa camelopardalis tippelskirchi is characterized by rather small spots separated by wide spaces (top left); G. camelopardalis reticulata, in contrast, is covered by large, closely spaced spots (top right). Both kinds of pattern can be accounted for by the author's reaction-diffusion model (bottom left and bottom right). The assumption is that at the time the pattern is laid down the embryo is between 35 and 45 days old and has a length of roughly eight to 10 centimeters. (The gestation period of the giraffe is about 457 days.)

effects, we considered the various

kinds of giraffe. For a given type

of pattern, we varied the parameter

that corresponds to the morphogen

threshold concentration for

melano-cyte activity. By varying the

parameter, we found we could

produce patterns that closely

resemble those of two different kinds

of giraffe [see illustration on opposite

page].

ecently the results of our model

-have been corroborated dramat-

ically by Charles M. Vest and

You-ren Xu of the University of

Michigan. They generated

standing-wave patterns on a

vibrating plate and changed the

nature of the patterns by changing the

frequency of vibration. The patterns

were made visible by a holographic

technique in which the plate was

bathed in laser light. Light reflected

from the plate interfered with a

reference beam, so that crests of

waves added to crests, troughs

added to troughs, and crests and

troughs canceled, and the resulting

pattern was recorded on a piece of

photographic emulsion [see illustra-

tion at right].

Vest and Youren found that low fre-

quencies of vibration produce simple

patterns and high frequencies of vi-

bration produce complex patterns.

The observation is interesting, be-

cause it has been shown that if a pat-

tern forms on a plate vibrating at a

given frequency, the pattern formed

on the same plate vibrated at a higher

frequency is identical with the

pattern formed on a proportionally

larger plate vibrated at the original

frequency. In other words, Vest and

Youren's data support our con-

clusion that more complex patterns

should be generated as the scale of

! the reaction-diffusion domain is in-

creased. The resemblance between

our patterns and the patterns subse-

quently produced by the Michigan

workers is striking.

I should like to stress again that

all the patterns generated were pro-

duced by varying only the scale and

geometry of the reaction domain; all

the other parameters were held fixed

(with the exception of the different

threshold concentrations in the case

of the giraffe). Even so, the diversity

of pattern is remarkable. The mod-

I el also suggests a possible explana-

tion for the various pattern anom-

alies seen in some animals. Under

some circumstances a change in the

[ value of one of the parameters can re-

sult in a marked change in the pat-

tern obtained. The size of the effect

R STANDING-WAVE PATTERNS generated on a thin vibrating plate resemble coat pat- terns and confirm the author's work. More complex patterns correspond to higher fre- quencies of vibration. The experiments were done by Charles M. Vest and Youren Xu.

depends on how close the value of the

parameter is to a bifurcation value: the

bolic rate are among some of them.

value at which a qualitative change in

Although the effects of such factors

probably could be mimicked by ma-

the pattern is generated.

If one of the parameters, say a rate

nipulating various parameters, there

constant in the reaction kinetics, is

is little point in doing so until more is

varied continuously, the mechanism

known about how the patterns re-

passes from a state in which no spatial

flected in the melanin pigments are

pattern can be generated to a patterned

actually produced. In the meantime

state and finally back to a state

one cannot help but note the wide va-

containing no patterns. The fact that

riety of patterns that can be generat-

such small changes in a parameter near

ed with a reaction-diffusion model by

a bifurcation value can result in such

varying only the scale and geometry.

large changes in pattern is consistent

The considerable circumstantial evi-

with the punctuated-equi-librium

dence derived from comparison with

theory of evolution. This theory holds

specific animal-pattern features is en-

that long periods of little evolutionary

couraging. I am confident that most

change are punctuated by short bursts

of the observed coat patterns can

be generated by a reaction-diffusion

of sudden and rapid change. mechanism. The fact that many gen-

eral and specific features of mamma-

any factors, of course, affect an-

imal coloration. Temperature,

humidity, diet, hormones and meta-

lian coat patterns can be explained

M by this simple theory, however, does

not make it right. Only experimental

observation can confirm the theory.