Developmentally Unique Cerebellar Processing Prioritizes Self- over Other-Generated Movements.

Animals must distinguish the sensory consequences of self-generated movements (reafference) from those of other-generated movements (exafference). Only self-generated movements entail the production of motor copies (i.e., corollary discharges), which are compared with reafference in the cerebellum to compute predictive or internal models of movement. Internal models emerge gradually over the first three postnatal weeks in rats through a process that is not yet fully understood. Previously, we demonstrated in postnatal day (P) 8 and P12 rats that precerebellar nuclei convey corollary discharge and reafference to the cerebellum during active (REM) sleep when pups produce limb twitches. Here, recording from a deep cerebellar nucleus (interpositus, IP) in P12 rats of both sexes, we compared reafferent and exafferent responses with twitches and limb stimulations, respectively. As expected, most IP units showed robust responses to twitches. However, in contrast with other sensory structures throughout the brain, relatively few IP units showed exafferent responses. Upon finding that exafferent responses occurred in pups under urethane anesthesia, we hypothesized that urethane inhibits cerebellar cortical cells, thereby disinhibiting exafferent responses in IP. In support of this hypothesis, ablating cortical tissue dorsal to IP mimicked the effects of urethane on exafference. Finally, the results suggest that twitch-related corollary discharge and reafference are conveyed simultaneously and in parallel to cerebellar cortex and IP. Based on these results, we propose that twitches provide opportunities for the nascent cerebellum to integrate somatotopically organized corollary discharge and reafference, thereby enabling the development of closed-loop circuits and, subsequently, internal models.

Parallel and Serial Sensory Processing in Developing Primary Somatosensory and Motor Cortex.

It is generally supposed that primary motor cortex (M1) receives somatosensory input predominantly via primary somatosensory cortex (S1). However, a growing body of evidence indicates that M1 also receives direct sensory input from the thalamus, independent of S1; such direct input is particularly evident at early ages before M1 contributes to motor control. Here, recording extracellularly from the forelimb regions of S1 and M1 in unanesthetized rats at postnatal day (P)8 and P12, we compared S1 and M1 responses to self-generated (i.e., reafferent) forelimb movements during active sleep and wake, and to other-generated (i.e., exafferent) forelimb movements. At both ages, reafferent responses were processed in parallel by S1 and M1; in contrast, exafferent responses were processed in parallel at P8 but serially, from S1 to M1, at P12. To further assess this developmental difference in processing, we compared exafferent responses to proprioceptive and tactile stimulation. At both P8 and P12, proprioceptive stimulation evoked parallel responses in S1 and M1, whereas tactile stimulation evoked parallel responses at P8 and serial responses at P12. Independent of the submodality of exafferent stimulation, pairs of S1-M1 units exhibited greater coactivation during active sleep than wake. These results indicate that S1 and M1 independently develop somatotopy before establishing the interactive relationship that typifies their functionality in adults.SIGNIFICANCE STATEMENT Learning any new motor task depends on the ability to use sensory information to update motor outflow. Thus, to understand motor learning, we must also understand how animals process sensory input. Primary somatosensory cortex (S1) and primary motor cortex (M1) are two interdependent structures that process sensory input throughout life. In adults, the functional relationship between S1 and M1 is well established; however, little is known about how S1 and M1 begin to transmit or process sensory information in early life. In this study, we investigate the early development of S1 and M1 as a sensory processing unit. Our findings provide new insights into the fundamental principles of sensory processing and the development of functional connectivity between these important sensorimotor structures.

Sensory Coding of Limb Kinematics in Motor Cortex across a Key Developmental Transition.

Primary motor cortex (M1) undergoes protracted development in mammals, functioning initially as a sensory structure. Throughout the first postnatal week in rats, M1 is strongly activated by self-generated forelimb movements-especially by the twitches that occur during active sleep. Here, we quantify the kinematic features of forelimb movements to reveal receptive-field properties of individual units within the forelimb region of M1. At postnatal day 8 (P8), nearly all units were strongly modulated by movement amplitude, especially during active sleep. By P12, only a minority of units continued to exhibit amplitude tuning, regardless of behavioral state. At both ages, movement direction also modulated M1 activity, though to a lesser extent. Finally, at P12, M1 population-level activity became more sparse and decorrelated, along with a substantial alteration in the statistical distribution of M1 responses to limb movements. These findings reveal a transition toward a more complex and informationally rich representation of movement long before M1 develops its motor functionality.SIGNIFICANCE STATEMENT Primary motor cortex (M1) plays a fundamental role in the generation of voluntary movements and motor learning in adults. In early development, however, M1 functions as a prototypical sensory structure. Here, we demonstrate in infant rats that M1 codes for the kinematics of self-generated limb movements long before M1 develops its capacity to drive movements themselves. Moreover, we identify a key transition during the second postnatal week in which M1 activity becomes more informationally complex. Together, these findings further delineate the complex developmental path by which M1 develops its sensory functions in support of its later-emerging motor capacities.

Active Sleep Promotes Coherent Oscillatory Activity in the Cortico-Hippocampal System of Infant Rats.

Active sleep (AS) provides a unique developmental context for synchronizing neural activity within and between cortical and subcortical structures. In week-old rats, sensory feedback from myoclonic twitches, the phasic motor activity that characterizes AS, promotes coherent theta oscillations (4-8 Hz) in the hippocampus and red nucleus, a midbrain motor structure. Sensory feedback from twitches also triggers rhythmic activity in sensorimotor cortex in the form of spindle bursts, which are brief oscillatory events composed of rhythmic components in the theta, alpha/beta (8-20 Hz), and beta2 (20-30 Hz) bands. Here we ask whether one or more of these spindle-burst components are communicated from sensorimotor cortex to hippocampus. By recording simultaneously from whisker barrel cortex and dorsal hippocampus in 8-day-old rats, we show that AS, but not other behavioral states, promotes cortico-hippocampal coherence specifically in the beta2 band. By cutting the infraorbital nerve to prevent the conveyance of sensory feedback from whisker twitches, cortical-hippocampal beta2 coherence during AS was substantially reduced. These results demonstrate the necessity of sensory input, particularly during AS, for coordinating rhythmic activity between these two developing forebrain structures.

Spatiotemporal organization of myoclonic twitching in sleeping human infants.

During the perinatal period in mammals when active sleep predominates, skeletal muscles twitch throughout the body. We have hypothesized that myoclonic twitches provide unique insight into the functional status of the human infant's nervous system. However, assessments of the rate and patterning of twitching have largely been restricted to infant rodents. Thus, here we analyze twitching in human infants over the first seven postnatal months. Using videography and behavioral measures of twitching during bouts of daytime sleep, we find at all ages that twitching across the body occurs predominantly in bursts at intervals of 10 s or less. We also find that twitching is expressed differentially across the body and with age. For example, twitching of the face and head is most prevalent shortly after birth and decreases over the first several months. In addition, twitching of the hands and feet occurs at a consistently higher rate than does twitching elsewhere in the body. Finally, the patterning of twitching becomes more structured with age, with twitches of the left and right hands and feet exhibiting the strongest coupling. Altogether, these findings support the notion that twitches can provide a unique source of information about typical and atypical sensorimotor development.

Hnrnph1 Is A Quantitative Trait Gene for Methamphetamine Sensitivity.

Psychostimulant addiction is a heritable substance use disorder; however its genetic basis is almost entirely unknown. Quantitative trait locus (QTL) mapping in mice offers a complementary approach to human genome-wide association studies and can facilitate environment control, statistical power, novel gene discovery, and neurobiological mechanisms. We used interval-specific congenic mouse lines carrying various segments of chromosome 11 from the DBA/2J strain on an isogenic C57BL/6J background to positionally clone a 206 kb QTL (50,185,512-50,391,845 bp) that was causally associated with a reduction in the locomotor stimulant response to methamphetamine (2 mg/kg, i.p.; DBA/2J < C57BL/6J)-a non-contingent, drug-induced behavior that is associated with stimulation of the dopaminergic reward circuitry. This chromosomal region contained only two protein coding genes-heterogeneous nuclear ribonucleoprotein, H1 (Hnrnph1) and RUN and FYVE domain-containing 1 (Rufy1). Transcriptome analysis via mRNA sequencing in the striatum implicated a neurobiological mechanism involving a reduction in mesolimbic innervation and striatal neurotransmission. For instance, Nr4a2 (nuclear receptor subfamily 4, group A, member 2), a transcription factor crucial for midbrain dopaminergic neuron development, exhibited a 2.1-fold decrease in expression (DBA/2J < C57BL/6J; p 4.2 x 10-15). Transcription activator-like effector nucleases (TALENs)-mediated introduction of frameshift deletions in the first coding exon of Hnrnph1, but not Rufy1, recapitulated the reduced methamphetamine behavioral response, thus identifying Hnrnph1 as a quantitative trait gene for methamphetamine sensitivity. These results define a novel contribution of Hnrnph1 to neurobehavioral dysfunction associated with dopaminergic neurotransmission. These findings could have implications for understanding the genetic basis of methamphetamine addiction in humans and the development of novel therapeutics for prevention and treatment of substance abuse and possibly other psychiatric disorders.

Wakefulness suppresses retinal wave-related neural activity in visual cortex.

In the developing visual system before eye opening, spontaneous retinal waves trigger bursts of neural activity in downstream structures, including visual cortex. At the same ages when retinal waves provide the predominant input to the visual system, sleep is the predominant behavioral state. However, the interactions between behavioral state and retinal wave-driven activity have never been explicitly examined. Here we characterized unit activity in visual cortex during spontaneous sleep-wake cycles in 9- and 12-day-old rats. At both ages, cortical activity occurred in discrete rhythmic bursts, ~30-60 s apart, mirroring the timing of retinal waves. Interestingly, when pups spontaneously woke up and moved their limbs in the midst of a cortical burst, the activity was suppressed. Finally, experimentally evoked arousals also suppressed intraburst cortical activity. All together, these results indicate that active wake interferes with the activation of the developing visual cortex by retinal waves. They also suggest that sleep-wake processes can modulate visual cortical plasticity at earlier ages than has been previously considered.NEW & NOTEWORTHY By recording in visual cortex in unanesthetized infant rats, we show that neural activity attributable to retinal waves is specifically suppressed when pups spontaneously awaken or are experimentally aroused. These findings suggest that the relatively abundant sleep of early development plays a permissive functional role for the visual system. It follows, then, that biological or environmental factors that disrupt sleep may interfere with the development of these neural networks.

Sensorimotor processing in the newborn rat red nucleus during active sleep.

Sensory feedback from sleep-related myoclonic twitches is thought to drive activity-dependent development in spinal cord and brain. However, little is known about the neural pathways involved in the generation of twitches early in development. The red nucleus (RN), source of the rubrospinal tract, has been implicated in the production of phasic motor activity during active sleep in adults. Here we hypothesized that the RN is also a major source of motor output for twitching in early infancy, a period when twitching is an especially abundant motor behavior. We recorded extracellular neural activity in the RN during sleep and wakefulness in 1-week-old unanesthetized rats. Neurons in the RN fired phasically before twitching and wake movements of the contralateral forelimb. A subpopulation of neurons in the RN exhibited a significant peak of activity after forelimb movement onset, suggesting reafferent sensory processing. Consistent with this observation, manual stimulation of the forelimb evoked RN responses. Unilateral inactivation of the RN using a mixture comprising GABAA, GABAB, and glycine receptor agonists caused an immediate and temporary increase in motor activity followed by a marked and prolonged decrease in twitching and wake movements. Altogether, these data support a causal role for the RN in infant motor behavior. Furthermore, they indicate that twitching, which is characterized by discrete motor output and reafferent input, provides an opportunity for sensorimotor integration and activity-dependent development of topography within the newborn RN.

Spontaneous activity and functional connectivity in the developing cerebellorubral system.

The development of the cerebellar system depends in part on the emergence of functional connectivity in its input and output pathways. Characterization of spontaneous activity within these pathways provides insight into their functional status in early development. In the present study we recorded extracellular activity from the interpositus nucleus (IP) and its primary downstream target, the red nucleus (RN), in unanesthetized rats at postnatal days 8 (P8) and P12, a period of dramatic change in cerebellar circuitry. The two structures exhibited state-dependent activity patterns and age-related changes in rhythmicity and overall firing rate. Importantly, sensory feedback (i.e., reafference) from myoclonic twitches (spontaneous, self-generated movements that are produced exclusively during active sleep) drove neural activity in the IP and RN at both ages. Additionally, anatomic tracing confirmed the presence of cerebellorubral connections as early as P8. Finally, inactivation of the IP and adjacent nuclei using the GABAA receptor agonist muscimol caused a substantial decrease in neural activity in the contralateral RN at both ages, as well as the disappearance of rhythmicity; twitch-related activity in the RN, however, was preserved after IP inactivation, indicating that twitch-related reafference activates the two structures in parallel. Overall, the present findings point to the contributions of sleep-related spontaneous activity to the development of cerebellar networks.

Twitch-related and rhythmic activation of the developing cerebellar cortex.

The cerebellum is a critical sensorimotor structure that exhibits protracted postnatal development in mammals. Many aspects of cerebellar circuit development are activity dependent, but little is known about the nature and sources of the activity. Based on previous findings in 6-day-old rats, we proposed that myoclonic twitches, the spontaneous movements that occur exclusively during active sleep (AS), provide generalized as well as topographically precise activity to the developing cerebellum. Taking advantage of known stages of cerebellar cortical development, we examined the relationship between Purkinje cell activity (including complex and simple spikes), nuchal and hindlimb EMG activity, and behavioral state in unanesthetized 4-, 8-, and 12-day-old rats. AS-dependent increases in complex and simple spike activity peaked at 8 days of age, with 60% of units exhibiting significantly more activity during AS than wakefulness. Also, at all three ages, approximately one-third of complex and simple spikes significantly increased their activity within 100 ms of twitches in one of the two muscles from which we recorded. Finally, we observed rhythmicity of complex and simple spikes that was especially prominent at 8 days of age and was greatly diminished by 12 days of age, likely due to developmental changes in climbing fiber and mossy fiber innervation patterns. All together, these results indicate that the neurophysiological activity of the developing cerebellum can be used to make inferences about changes in its microcircuitry. They also support the hypothesis that sleep-related twitches are a prominent source of discrete climbing and mossy fiber activity that could contribute to the activity-dependent development of this critical sensorimotor structure.

A valuable and promising method for recording brain activity in behaving newborn rodents.

Neurophysiological recording of brain activity has been critically important to the field of neuroscience, but has contributed little to the field of developmental psychobiology. The reasons for this can be traced largely to methodological difficulties associated with recording neural activity in behaving newborn rats and mice. Over the last decade, however, the evolution of methods for recording from head-fixed newborns has heralded a new era in developmental neurophysiology. Here, we review these recent developments and provide a step-by-step primer for those interested in applying the head-fix method to their own research questions. Until now, this method has been used primarily to investigate spontaneous brain activity across sleep and wakefulness, the contributions of the sensory periphery to brain activity, or intrinsic network activity. Now, with some ingenuity, the uses of the head-fix method can be expanded to other domains to benefit our understanding of brain-behavior relations under normal and pathophysiological conditions across early development.

Genetic analysis in the Collaborative Cross breeding population.

Genetic reference populations in model organisms are critical resources for systems genetic analysis of disease related phenotypes. The breeding history of these inbred panels may influence detectable allelic and phenotypic diversity. The existing panel of common inbred strains reflects historical selection biases, and existing recombinant inbred panels have low allelic diversity. All such populations may be subject to consequences of inbreeding depression. The Collaborative Cross (CC) is a mouse reference population with high allelic diversity that is being constructed using a randomized breeding design that systematically outcrosses eight founder strains, followed by inbreeding to obtain new recombinant inbred strains. Five of the eight founders are common laboratory strains, and three are wild-derived. Since its inception, the partially inbred CC has been characterized for physiological, morphological, and behavioral traits. The construction of this population provided a unique opportunity to observe phenotypic variation as new allelic combinations arose through intercrossing and inbreeding to create new stable genetic combinations. Processes including inbreeding depression and its impact on allelic and phenotypic diversity were assessed. Phenotypic variation in the CC breeding population exceeds that of existing mouse genetic reference populations due to both high founder genetic diversity and novel epistatic combinations. However, some focal evidence of allele purging was detected including a suggestive QTL for litter size in a location of changing allele frequency. Despite these inescapable pressures, high diversity and precision for genetic mapping remain. These results demonstrate the potential of the CC population once completed and highlight implications for development of related populations.

Developmentally unique cerebellar processing prioritizes self-over other-generated movements.

Animals must distinguish the sensory consequences of self-generated movements (reafference) from those of other-generated movements (exafference). Only self-generated movements entail the production of motor copies (i.e., corollary discharges), which are compared with reafference in the cerebellum to compute predictive or internal models of movement. Internal models emerge gradually over the first three postnatal weeks in rats through a process that is not yet fully understood. Previously, we demonstrated in postnatal day (P) P8 and P12 rats that precerebellar nuclei convey corollary discharge and reafference to the cerebellum during active (REM) sleep when pups produce limb twitches. Here, recording from a deep cerebellar nucleus (interpositus, IP) in P12 rats of both sexes, we compared reafferent and exafferent responses to twitches and limb stimulations, respectively. As expected, most IP units showed robust responses to twitches. However, in contrast with other sensory structures throughout the brain, relatively few IP units showed exafferent responses. Upon finding that exafferent responses occurred in pups under urethane anesthesia, we hypothesized that urethane inhibits cerebellar cortical cells, thereby disinhibiting exafferent responses in IP. In support of this hypothesis, ablating cortical tissue dorsal to IP mimicked the effects of urethane on exafference. Finally, the results suggest that twitch-related corollary discharge and reafference are conveyed simultaneously and in parallel to cerebellar cortex and IP. Based on these results, we propose that twitches provide opportunities for the nascent cerebellum to integrate somatotopically organized corollary discharge and reafference, thereby enabling the development of closed-loop circuits and, subsequently, internal models.

Coincident development and synchronization of sleep-dependent delta in the cortex and medulla.

In early development, active sleep is the predominant sleep state before it is supplanted by quiet sleep. In rats, the developmental increase in quiet sleep is accompanied by the sudden emergence of the cortical delta rhythm (0.5-4 Hz) around postnatal day 12 (P12). We sought to explain the emergence of the cortical delta by assessing developmental changes in the activity of the parafacial zone (PZ), a medullary structure thought to regulate quiet sleep in adults. We recorded from the PZ in P10 and P12 rats and predicted an age-related increase in neural activity during increasing periods of delta-rich cortical activity. Instead, during quiet sleep, we discovered sleep-dependent rhythmic spiking activity-with intervening periods of total silence-phase locked to a local delta rhythm. Moreover, PZ and cortical delta were coherent at P12 but not at P10. PZ delta was also phase locked to respiration, suggesting sleep-dependent modulation of PZ activity by respiratory pacemakers in the ventral medulla. Disconnecting the main olfactory bulbs from the cortex did not diminish cortical delta, indicating that the influence of respiration on delta at this age is not mediated indirectly through nasal breathing. Finally, we observed an increase in parvalbumin-expressing terminals in the PZ across these ages, supporting a role for local GABAergic inhibition in the PZ's rhythmicity. The unexpected discovery of delta-rhythmic neural activity in the medulla-when cortical delta is also emerging-provides a new perspective on the brainstem's role in regulating sleep and promoting long-range functional connectivity in early development.

DELTA-RHYTHMIC ACTIVITY IN THE MEDULLA DEVELOPS COINCIDENT WITH CORTICAL DELTA IN SLEEPING INFANT RATS.

In early development, active sleep is the predominant sleep state before it is supplanted by quiet sleep. In rats, the developmental increase in quiet sleep is accompanied by the sudden emergence of the cortical delta rhythm (0.5-4 Hz) around postnatal day 12 (P12). We sought to explain the emergence of cortical delta by assessing developmental changes in the activity of the parafacial zone (PZ), a medullary structure thought to regulate quiet sleep in adults. We recorded from PZ in P10 and P12 rats and predicted an age-related increase in neural activity during increasing periods of delta-rich cortical activity. Instead, during quiet sleep we discovered sleep-dependent rhythmic spiking activity-with intervening periods of total silence-phase-locked to a local delta rhythm. Moreover, PZ and cortical delta were coherent at P12, but not at P10. PZ delta was also phase-locked to respiration, suggesting sleep-dependent modulation of PZ activity by respiratory pacemakers in the ventral medulla. Disconnecting the main olfactory bulbs from the cortex did not diminish cortical delta, indicating that the influence of respiration on delta at this age is not mediated indirectly through nasal breathing. Finally, we observed an increase in parvalbumin-expressing terminals in PZ across these ages, supporting a role for GABAergic inhibition in PZ's rhythmicity. The discovery of delta-rhythmic neural activity in the medulla-when cortical delta is also emerging-opens a new path to understanding the brainstem's role in regulating sleep and synchronizing rhythmic activity throughout the brain.

Cortical Representation of Movement Across the Developmental Transition to Continuous Neural Activity.

Primary motor cortex (M1) exhibits a protracted period of development that includes the establishment of a somatosensory map long before motor outflow emerges. In rats, the sensory representation is established by postnatal day (P) 8 when cortical activity is still "discontinuous." Here, we ask how the representation survives the sudden transition to noisy "continuous" activity at P12. Using neural decoding to predict forelimb movements based solely on M1 activity, we show that a linear decoder is sufficient to predict limb movements at P8, but not at P12; in contrast, a nonlinear decoder effectively predicts limb movements at P12. The change in decoder performance at P12 reflects an increase in both the complexity and uniqueness of kinematic information available in M1. We next show that the representation at P12 is more susceptible to the deleterious effects of "lesioning" inputs and to "transplanting" M1's encoding scheme from one pup to another. We conclude that the emergence of continuous cortical activity signals the developmental onset in M1 of more complex, informationally sparse, and individualized sensory representations.

Neural decoding reveals specialized kinematic tuning after an abrupt cortical transition.

The primary motor cortex (M1) exhibits a protracted period of development, including the development of a sensory representation long before motor outflow emerges. In rats, this representation is present by postnatal day (P) 8, when M1 activity is "discontinuous." Here, we ask how the representation changes upon the transition to "continuous" activity at P12. We use neural decoding to predict forelimb movements from M1 activity and show that a linear decoder effectively predicts limb movements at P8 but not at P12; instead, a nonlinear decoder better predicts limb movements at P12. The altered decoder performance reflects increased complexity and uniqueness of kinematic information in M1. We next show that M1's representation at P12 is more susceptible to "lesioning" of inputs and "transplanting" of M1's encoding scheme from one pup to another. Thus, the emergence of continuous M1 activity signals the developmental onset of more complex, informationally sparse, and individualized sensory representations.

Genome-wide association for fear conditioning in an advanced intercross mouse line.

Fear conditioning (FC) may provide a useful model for some components of post-traumatic stress disorder (PTSD). We used a C57BL/6J × DBA/2J F(2) intercross (n = 620) and a C57BL/6J × DBA/2J F(8) advanced intercross line (n = 567) to fine-map quantitative trait loci (QTL) associated with FC. We conducted an integrated genome-wide association analysis in QTLRel and identified five highly significant QTL affecting freezing to context as well as four highly significant QTL associated with freezing to cue. The average percent decrease in QTL width between the F(2) and the integrated analysis was 59.2%. Next, we exploited bioinformatic sequence and expression data to identify candidate genes based on the existence of non-synonymous coding polymorphisms and/or expression QTLs. We identified numerous candidate genes that have been previously implicated in either fear learning in animal models (Bcl2, Btg2, Dbi, Gabr1b, Lypd1, Pam and Rgs14) or PTSD in humans (Gabra2, Oprm1 and Trkb); other identified genes may represent novel findings. The integration of F(2) and AIL data maintains the advantages of studying FC in model organisms while significantly improving resolution over previous approaches.

Fine-mapping alleles for body weight in LG/J × SM/J F2 and F(34) advanced intercross lines.

The present study measured variation in body weight using a combined analysis in an F(2) intercross and an F(34) advanced intercross line (AIL). Both crosses were derived from inbred LG/J and SM/J mice, which were selected for large and small body size prior to inbreeding. Body weight was measured at 62 (± 5) days of age. Using an integrated GWAS and forward model selection approach, we identified 11 significant QTLs that affected body weight on ten different chromosomes. With these results we developed a full model that explained over 18% of the phenotypic variance. The median 1.5-LOD support interval was 5.55 Mb, which is a significant improvement over most prior body weight QTLs. We identified nonsynonymous coding SNPs between LG/J and SM/J mice in order to further narrow the list of candidate genes. Three of the genes with nonsynonymous coding SNPs (Rad23b, Stk33, and Anks1b) have been associated with adiposity, waist circumference, and body mass index in human GWAS, thus providing evidence that these genes may underlie our QTLs. Our results demonstrate that a relatively small number of loci contribute significantly to the phenotypic variance in body weight, which is in marked contrast to the situation in humans. This difference is likely to be the result of strong selective pressure and the simplified genetic architecture, both of which are important advantages of our system.