Environmental selection during the last ice age on the mother-to-infant transmission of vitamin D and fatty acids through breast milk

Leslea J. Hlusko, Joshua P. Carlson, George Chaplin, Scott A. Elias, John F. Hoffecker, Michaela Huffman, Nina G. Jablonski, Tesla A. Monson, Dennis H. O’Rourke, Marin A. Pilloud, and G. Richard Scott
PNAS May 8, 2018 115 (19) E4426-E4432; published ahead of print April 23, 2018 https://doi.org/10.1073/pnas.1711788115

  1. Edited by C. Owen Lovejoy, Kent State University, Kent, OH, and approved March 23, 2018 (received for review June 30, 2017)

Significance

The frequency of the human-specific EDAR V370A isoform is highly elevated in North and East Asian populations. The gene is known to have several pleiotropic effects, among which are sweat gland density and ductal branching in the mammary gland. The former has led some geneticists to argue that the near-fixation of this allele was caused by selection for modulation of thermoregulatory sweating. We provide an alternative hypothesis, that selection instead acted on the allele’s effect of increasing ductal branching in the mammary gland, thereby amplifying the transfer of critical nutrients to infants via mother’s milk. This is likely to have occurred during the Last Glacial Maximum when a human population was genetically isolated in the high-latitude environment of the Beringia.

Abstract

Because of the ubiquitous adaptability of our material culture, some human populations have occupied extreme environments that intensified selection on existing genomic variation. By 32,000 years ago, people were living in Arctic Beringia, and during the Last Glacial Maximum (LGM; 28,000–18,000 y ago), they likely persisted in the Beringian refugium. Such high latitudes provide only very low levels of UV radiation, and can thereby lead to dangerously low levels of biosynthesized vitamin D. The physiological effects of vitamin D deficiency range from reduced dietary absorption of calcium to a compromised immune system and modified adipose tissue function. The ectodysplasin A receptor (EDAR) gene has a range of pleiotropic effects, including sweat gland density, incisor shoveling, and mammary gland ductal branching. The frequency of the human-specific EDAR V370A allele appears to be uniquely elevated in North and East Asian and New World populations due to a bout of positive selection likely to have occurred circa 20,000 y ago. The dental pleiotropic effects of this allele suggest an even higher occurrence among indigenous people in the Western Hemisphere before European colonization. We hypothesize that selection on EDAR V370A occurred in the Beringian refugium because it increases mammary ductal branching, and thereby may amplify the transfer of critical nutrients in vitamin D-deficient conditions to infants via mothers’ milk. This hypothesized selective context for EDAR V370A was likely intertwined with selection on the fatty acid desaturase (FADS) gene cluster because it is known to modulate lipid profiles transmitted to milk from a vitamin D-rich diet high in omega-3 fatty acids.

From Thomas Jefferson’s archaeological excavations (1) to modern genomics (2⇓–4), scientists have been fascinated by the first migration of humans into the Americas. A myriad of new evidence reveals that the earliest people in the Western Hemisphere dispersed from a population that lived in genetic isolation for thousands of years on the exposed Beringian platform in the Arctic during the Last Glacial Maximum [LGM; 28,000–18,000 y ago (2⇓⇓⇓⇓–7)]. The Arctic is an extreme environment because of the very low UV radiation (UV) reaching the earth’s surface at such high latitude. UV is essential to almost all life forms because it catalyzes biochemical processes, especially the synthesis of vitamin D (8).

Extreme environments can impact genetic variation and provide opportunities to elucidate relationships between genotype and phenotype (9). The classic human example is the range of our physiological adaptations to the hypoxic conditions of high altitude, such as are experienced in the highlands of Ethiopia, the Tibetan Plateau, and the Andean Mountains (10⇓–12). Evidence shows that the populations which have long-occupied these high-elevation regions have an increased frequency of red blood cell polymorphisms that likely underlie these physiologies (10⇓–12).

Here, we investigate whether the population occupying Beringia during the LGM represents another example of human adaptation to an extreme environment, this time adapting to very low UV exposure (Fig. 1). There are two lines of genetic evidence for this: variation in the fatty acid desaturase (FADS) gene cluster that modulates the manufacture of polyunsaturated fatty acids and variation in the ectodysplasin A receptor (EDAR) gene that influences ectodermally derived structures, such as teeth, hair, and mammary gland ductal branching. A study on selection on the FADS gene cluster in the ancestral population of Native Americans has been published previously (13), but, here, we shift the emphasis from phenotypic effects on older adults to focus on those that influence fertility via breast milk. We then present evidence that EDAR may have undergone an episode of selection in the same population, likely due to its influence on mammary ductal branching. We hypothesize that the genetically isolated population which occupied Beringia during the LGM experienced selection for an increase in vitamin D in breast milk in response to the low UV environment. Traces of this previously intense selection appear to still be present in the genetic variation of Native American and North and East Asian populations today.

Fig. 1.
Fig. 1.

Geography of Beringia and levels of UV radiation. (A) Map of Beringia today. Cross-hatching indicates the region in which levels of UVMED (defined as the amount of UV radiation that will produce minimal erythema) that reach the Earth’s surface are too low to promote cutaneous synthesis of vitamin D in humans on a year-by-year basis, requiring dietary supplementation (modified from ref. 84; projected to show an equal distance map of Beringia). The black and white region marks the Arctic Circle, for which there are no Total Ozone Mapping Satellite Data version 7. Other data show that this region has even less UV-B exposure, as would be expected from the increased latitude. The areas below the Arctic Circle in white and light blue are shallow seas as discerned from modern bathymetry using Etopo2 data, indicating land that would have been exposed during the LGM. (B) Map of Beringia during the LGM showing the exposure of land at 117 m below current sea level and the reconstructed terrestrial environments. The shrub tundra is the only area biologically productive enough to support a human population of the size estimated by molecular data. This population was genetically isolated for ∼2,500–9,000 y during the LGM because of the ice to the east and extensive mesic tundra to the west.

The FADS Gene Cluster

Allelic variation in the FADS gene complex corresponds strongly to geographic ancestry (14, 15), and is the most pronounced allele frequency difference between the Greenlandic Inuit and other human populations (16). Because modern Western cultures place much attention on the cardiovascular disease risks associated with diets rich in omega-6 versus omega-3 long-chain polyunsaturated fatty acids (LC-PUFAs) (17, 18), the high frequencies of these alleles in the Inuit have been primarily interpreted as an adaptation to their traditionally meat-rich (and omega-3–rich) diet (16, 19). However, and surprisingly, 95% of native Central and South Americans also show evidence of selection on the same FADS polymorphisms, yet these cultural groups do not traditionally consume diets nearly as rich in omega-3 fatty acids as do the Inuit, suggesting that their common ancestral population is characterized by one which did (13).

From an evolutionary perspective, the possibility that the selective force favoring these FADS alleles was a healthier ratio of omega-3 to omega-6 fatty acids in adults would indicate that the effect must have been very strong, as selective pressure is relatively low on genetic variation that influences later and post-reproductive years of life (20). In contrast, even very small positive effects on fertility and childhood survival have far greater selective power (20). Considering this, the selective benefits of these FADS polymorphisms are more likely connected to their significant role in modulating the relative levels of omega-3 and omega-6 fatty acids during gestation and postnatal growth of infants.

To elaborate, the fatty acids arachidonic acid, docosahexaenoic acid (DHA), and eicosapentaenoic acid (EPA) are essential to cognitive and visual development (21). Maternal polymorphisms in the FADS gene cluster have strong effects on the levels of EPA and DHA in breast milk [in humans (22⇓⇓⇓⇓–27) and in other mammals (27, 28)]. This suggests that lipid levels in a mother’s diet are optimized by specific polymorphisms in her FADS gene cluster during milk synthesis (Fig. 2A). These polymorphisms are therefore a potentially strong target of selection because fetuses and infants have a very limited ability to synthesize LC-PUFAs on their own and are dependent on their mother’s genotype to modulate their relative proportions (22⇓⇓⇓⇓–27). We propose that the phenotypic effects of FADS polymorphisms on the LC-PUFA content of breast milk was the primary target of selection in past human populations that consumed diets with compromised proportions of omega-3 and omega-6 fatty acids, such as is seen in the traditional diets of people living in the Arctic.

Fig. 2.
Fig. 2.

Overview of breast anatomy, development, milk production, and histology. (A) Genetic influences on milk production summarized from the main text. The rodent, human, and cow figures indicate the systems from which the results are known. (B) Three stages of mammary development alongside the genetic mutations known to increase ductal branching. Hormones induce the pubertal and gestation/lactation stages (denoted in pink). In vitamin D-deficient adipose tissue, ductal branching increases during these two stages of development. The embryonic stage (in green) is not hormonally induced; branching density is influenced by the ectodysplasin pathway (108) and, specifically, increased by EDAR V370A. (C) Cross-section of a ductal lumen showing the bilayer of luminal and myoepithelial cells. (D) Schematic of a desmosome, one of the main adhesive structures between mammary luminal cells. The orange lines represent desmosomal cadherins in the extracellular space. The gray lines in the purple area represent the intracellular filaments. (E) Close-up view of the desmocolin proteins that comprise the desmo-adhesome. An increase in the activity of the ectodyplasin pathway alters the relative proportions of Dsc2 and Dsc3, leading to a reduction in the adhesive strength.

EDAR.

The human EDAR V370A variant encodes a change in amino acid sequence in the highly conserved death domain of EDAR (29⇓⇓–32). EDAR functions as a protein receptor on a cell’s surface that activates the transcription factor NF-ĸB in the ectodysplasin pathway (29). Comparative studies show this pathway is functionally conserved across virtually all vertebrates, playing essential roles in the development of ectodermal structures from bird feathers to fish scales (33, 34).

EDAR is one of four genes implicated in hypohidrotic ectodermal dysplasia (HED), a set of ∼150 syndromes characterized by mild to severe defects in ectodermally derived structures, such as hair, teeth, breasts, and sweat glands (33, 35⇓–37). In contrast to HED, the EDAR V370A allele has the opposite effect. Genome-wide association studies (GWASs) of multiple Asian populations show that EDAR V370A is correlated with hair shaft caliber (38, 39), earlobe and chin shape (40), and a suite of morphological variants on teeth (40⇓⇓–43). Knock-in mouse studies reveal that the V370A allele leads to a twofold increase in NF-ĸB activation (30, 31, 44). Just like humans, the EDAR V370A mouse has thicker hair shafts, an increase in the branching density of mammary gland ducts, and (inconsistently) an increase in the number of eccrine glands on the footpads (30, 31).

Several studies conclude that EDAR V370A experienced a bout of intense positive selection ∼20,000 y ago in northern China (29, 31, 45). Previous interpretations as to why this selection occurred focused on the associated increase in eccrine sweat gland density on the fingertips, concluding that the selection was for improved thermoregulatory sweating during a warm spell during the LGM (31) or that the increase in sebaceous glands would offer protection from the cold, dry air of the LGM (30). These interpretations rely primarily on the current distribution of the allele in living populations (Fig. 3A), which shows high frequencies of EDAR V370A in North and East Asian and Native American populations but virtual absence in other populations around the world (29⇓–31, 45).

Fig. 3.
Fig. 3.

Relationship between EDAR V370A and degree of incisor shoveling. (A) Current world-wide allele frequencies of EDAR V370A (in yellow) and other EDAR haplotypes (in blue) (data from ref. 29). Because these data are from living people, the Native American data include significant admixture from European colonization, a population with essentially no occurrence of EDAR V370A. These modern data likely vastly underrepresent the occurrence of EDAR V370A among indigenous people before European contact. Current Asian allele frequencies may be higher than pre-LGM frequencies due to back-migration from the Beringian population after the end of the LGM, a migration supported by linguistic data (46). (B) Histogram showing EDAR V370A genotype and degree of incisor shoveling, demonstrating the imperfect but strongly additive nature of EDAR V370A’s influence on incisor shoveling (adapted from ref. 43). (C) Frequencies of incisor shoveling scores observed in archaeological populations from Africa, Europe, South Asia, East Asia, North America, and South America (Table S1). Purple and blue represent a lack of incisor shoveling, and as such, an individual who is EDAR V370A−/−. Note the lack of shoveling scores of 0 and 1, and very low occurrence of 2, in the indigenous people of the Western Hemisphere, indicating a very high frequency of EDAR V370A before European contact.

The post-LGM diaspora into the Americas and the subsequent centuries of European colonization (and genetic admixture) have dramatically overwritten the allelic variation of ∼20,000 y ago when EDAR V370A experienced a bout of selective pressure. To best interpret the evolutionary history of EDAR V370A, we need information on past allele frequencies. This task is difficult from a genome sequence approach because of the paucity of ancient DNA available. However, the pleiotropic effects of EDAR V370A provide a unique opportunity to more precisely reconstruct its evolutionary history without DNA sequence data.

Temporal and Geographic Variation of EDAR V370A

As noted previously, EDAR V370A is correlated with a suite of morphological variants on teeth, the most notable of which is incisor shoveling (40⇓⇓–43). GWASs report a close relationship between genotype (EDAR V370A) and this dental phenotype in living Chinese, Japanese, and Koreans (40⇓⇓–43). In these populations, incisor shoveling of any degree correlates with the presence of one or two copies of the EDAR V370A allele in an imperfect but clearly additive manner (41⇓–43) (Fig. 3B).

We scored the degree of expression in maxillary incisor shoveling for 5,333 people from >54 archaeological populations from across Europe, Asia, and North and South America (Supporting Information and Tables S1 and S2). Incisor shoveling is absent or expressed only to a slight degree outside of North and East Asian populations (where EDAR V370A frequencies are close to zero) but close to ubiquity in Native American populations, a pattern long interpreted as the result of an unknown selective pressure (47, 48) (Fig. 3C).

Considering the GWAS results that indicate an additive effect of EDAR V370A’s contribution to incisor shoveling, we interpreted a score of 1 or greater as evidence of an individual likely carrying at least one copy of EDAR V370A. There are no reported instances of a score of 0 among the 3,183 individuals assessed from North and South America, and as such, no evidence of any Native Americans before European contact without at least one copy of the EDAR V370A allele (Fig. 3C and Supporting Information). A series of ANOVAs comparing the frequency of each shoveling grade (range: 0–7) across geographic regions demonstrates that a score of 1 has the greatest variation between geographic groups, further suggesting that a score of 1 or higher is reflective of the presence of at least one copy of EDAR V370A (Supporting Information and Table S3).

These phenotypic data strongly suggest that incisor shoveling (and, concomitantly, the EDAR V370A allele) reached near-fixation in the population ancestral to all indigenous people of the Western Hemisphere. Given the timing of the dispersal into the Americas ∼17,000 y ago (3, 49), these dental data support the conclusion that EDAR V370A underwent positive selection in the ancestral Native American population during the LGM, similar to what we see with the FADS gene cluster (13), rather than in a population in what is now China (our evolutionary quantitative genetic analysis further supports this; Supporting Information).

The Ancestral Population of Native Americans and the LGM

Indigenous people in the Western Hemisphere derive from people who occupied latitudes above 55°N in Asia before 40,000 y ago (50). Well-dated archaeological sites near the mouth of the Yana River indicate a year-round adaptation to the Beringian Arctic Zone (latitude of 70°N) by 32,000 calibrated radiocarbon years (cal) B.P. (51).

Between 28,000 and 18,000 y ago, during the LGM, as aridity increased and biological productivity reduced dramatically (52), plants and animals in many parts of the world retreated to refugia, leaving signals of genetic isolation and hybridization that are documented widely by molecular ecologists (53, 54). Human populations similarly abandoned arid regions in Africa, Eurasia, and Australia, as indicated by genetic bottlenecks and evidence of local settlement hiatus (e.g., refs. 55, 56).

Various lines of evidence suggest that the Native American lineages during the LGM were genetically isolated in one or more refugia within Beringia (2, 5⇓–7), and at least one archaeological site appears to confirm a human presence in northeastern Beringia (latitude of 67°N) during that time frame (57). Analyses of Native American genetic variation indicate that the ancestral population was isolated for up to 9,000 y before dispersal in North and South America (after 15,000 cal B.P.), with an effective population size of at least a few thousand individuals (3, 49, 58).

The “Beringian Standstill” model places the ancestral Native American population in the Beringian refugium during much or all of the LGM (2, 5⇓–7), where it was geographically isolated from the Siberian population by uninhabitable areas in northeastern Asia and unable to expand into the Americas because of the coalesced Laurentide and Cordilleran ice sheets (Fig. 1B). The genetic isolation is a key point, as evolution occurs more readily in the absence of gene flow (59). The FADS polymorphisms and EDAR V370A likely existed in the pre-LGM ancestral populations [although the allele is not present in Denisovans (60)]. However, there was a significant opportunity for selection to shift the frequencies of functional alleles dramatically when this population became genetically isolated in the extreme UV environment of Beringia, including Arctic Beringia (above 66°N).

Biological Consequences of Low UV

Most vitamin D necessary for human health and reproductive success is produced through biosynthesis initiated by skin exposure to UV-B photons (61, 62). The lower rates of cutaneous synthesis of vitamin D that result from low UV environments have a wide range of health consequences. Vitamin D is among a large family of fat-soluble micronutrients that accumulate in adipose tissue proportional to circulating serum levels, and is released slowly as serum levels reduce (63). The bioactive form of vitamin D plays a secosteroid hormonal role in permitting absorption of dietary calcium through the lining of the gut (64). Through the Vitamin D receptor (Vdr), vitamin D also regulates expression of more than 220 genes in the human genome, with significant involvement in immune function and autoimmune disorders (65, 66). It is therefore not surprising that variation in immune function and disease risk correlates with latitude and UV-B exposure (67). Vitamin D also plays an important immunomodulatory role at the maternal/fetal interface (the placenta) (68) and is associated with sex-specific variation in birth weight (69).

However important the immunological role of vitamin D, adipose tissue appears to be the main phenotypic target of Vdr (63, 65, 66, 70). Vitamin D deficiency significantly compromises the metabolic function of adipose tissue, with a wide range of deleterious health effects (63). Body fat has a remarkable degree of phenotypic plasticity that enables a wide range of essential physiological functions, from maintaining body heat, to responding to cold, to forming the pink adipocytes in the mammary glands during pregnancy in anticipation of lactation (71⇓–73) (Fig. 2A).

Nutritional studies have long shown interrelationships between vitamin D and LC-PUFAs (74⇓–76). For example, the vitamin D synthesized through skin exposure to UV-B (vitamin D3) has positive effects on lipid profiles that cannot be replicated with dietary vitamin D2 supplementation (77). These interrelationships result from the effects of vitamin D on adipose tissue and immunological function, which are modulated by environmental exposure to UV-B.

Human populations that dispersed into geographic regions above 30° latitude must have experienced selection for reduced cutaneous eumelanin pigmentation to facilitate vitamin D biosynthesis (78, 79). Selection for depigmented skin phenotypes has occurred at least twice in modern human evolution and once in Neanderthals (78, 80), with selective pressure readily detected in genomic analyses (81⇓–83).

At latitudes above 48° (e.g., 450 km north of Hokkaido, Japan; Fig. 1A), sunlight reaching the Earth’s surface contains almost no UV-B except for low levels at or near the summer solstice (84, 85). Year-round human habitation of even more extreme high latitudes leads to vitamin D deficiency that cannot be counterbalanced by depigmentation alone (84) (Fig. 1A). In these geographic regions, diets centered on vitamin D-rich foods, such as marine mammals, oily fish, reindeer, or caribou, have been adopted across a range of traditional Arctic cultures (85). Archaeological evidence shows that these dietary innovations correlate with the earliest occupation of these latitudes (55, 86, 87) and shifts in allele frequencies of the genes involved in fatty acid synthesis (16).

Dietary practices can mitigate vitamin D deficiency for older children and adults, but pregnant and nursing mothers and their breast-feeding infants are still at significant risk because they need particularly high levels of vitamin D to avoid the litany of deleterious effects of vitamin D insufficiency (88, 89).

Fat and Vitamin D Content in Milk

Human breast milk varies greatly in fat content, with the relative proportions of omega-3 and omega-6 fatty acids varying by ethnicity independent of maternal diet (90), most likely because of polymorphisms in the FADS gene cluster as discussed earlier (22⇓⇓⇓⇓–27). The LC-PUFAs in milk are drawn from the specialized pink adipocytes found only in mammary glands, which are characterized by a higher number of intracellular lipid droplets (71⇓–73) and are dependent on healthy vitamin D levels to function properly (63).

Vitamin D content in human milk also varies between individuals and across populations (90⇓⇓–93). Although the effects of dietary supplementation, sun exposure, and genetic variation affect the vitamin’s concentration within an individual, there is nevertheless a strong correlation between maternal and infant circulating vitamin D levels (94, 95), and both mother and infant levels are related to each other’s sun exposure (93, 96). Importantly, vitamin D levels in milk are about threefold lower than the circulating levels in the mother (97). Endogenously produced vitamin D3 from UV-B exposure appears to be preferentially transferred into milk compared with the vitamin D2 consumed through diet (93). It is thus reasonable to suggest that infants in extremely low UV environments would be particularly dependent on the mother’s potential capacity to transfer vitamin D from her adipose stores across the placenta and also into her breast milk.

Fat, Vitamin D, and the Mammary Gland

Given that vitamin D is intricately intertwined with adipose physiology, it is not surprising that vitamin D plays a role in virtually all aspects of mammary gland function, as evinced by the presence and activation of the Vdr in mesenchymal and epithelial breast tissue (98, 99). Vitamin D also appears to influence ductal branching and, with it, milk content.

There are three distinct phases of breast development (Fig. 2B): (i) hormone-independent embryonic development during which buds form, sprout, and start to branch into the precursor of the fatty stroma to form a ductal tree; (ii) the onset of pubertal hormonal signaling when terminal end buds form at the tips of the mammary ducts and start to invade the fat pad from which they will extract milk fats, vitamins, and other nutrients (100⇓–102); and (iii) the repeatable hormone-induced cycles of pregnancy, lactation, and weaning (101, 102).

Vitamin D-deficient adipose tissue has been modeled by an adipose-specific Vdr knockout (Vdr KO) mouse, revealing significant and, from an evolutionary perspective, fascinating sex-by-diet effects (99). Female mice that are unable to bioactivate vitamin D in adipocytes exhibit increased visceral adipose tissue overall, but they have a normal mammary fat pad mass. Although the mass of the mammary fat pad is unaffected by the loss of Vdr, there is an increase in the ductal branching density within it (99). Surprisingly, however, this only occurred when the mice are fed a high-fat diet (reminiscent of traditional human diets in Artic cultures) (99). These results suggest that ductal branching is stimulated by adjacent vitamin D-deficient adipose tissue in response to high-fat diets during the hormonally induced stages of mammary development (98, 99, 103). In sharp contrast to the results for female mice, there are no effects in males with the adipose-specific Vdr deletion in either the normal or high-fat diet (99).

The results of the Vdr KO study, bolstered by additional research (reviewed in ref. 73), reveal a sex-specific relationship between vitamin D deficiency and adipose tissue that includes an increase in mammary ductal branching during the hormone-induced stages of breast development (Fig. 2B). This physiological increase in ductal branching in response to very low levels of vitamin D indicates that other genetic mechanisms that increase ductal branching may be positively selected for in populations experiencing chronic vitamin D deficiency, such as in UV-poor environments at high latitudes.

A completely independent line of evidence suggests that variation in ductal branching influences milk content (Fig. 2A). Variation in milk vitamin D levels is heritable in cows and shows evidence of positive selection over the short time frames associated with animal domestication, demonstrating that the underlying genetic mechanisms are highly responsive to selective pressure (104, 105). Milk vitamin D content varies by cattle breed, with some producing notably higher concentrations irrespective of UV exposure (reviewed in ref. 106). Functional associations of 31 genes expressed in cows producing different milk protein and fat concentrations include transcriptomes of genes involved in mammary gland bud elongation (107), suggesting that mammary duct branching contributes to the genetic underpinnings of variation in milk composition.

EDAR V370A’s Influence on Ductal Branching

The relationship between vitamin D and ductal branching calls for a renewed consideration of EDAR V370A’s influence on mammary duct branching density as the phenotypic target of selection. As with other ectodermally derived structures (e.g., teeth), the earliest phase of mammary development results from epithelial–mesenchymal tissue interactions at embryonic day 11 in the mouse (100⇓–102). Ductal branching is established via NF-ĸB signaling, a downstream target of the ectodysplasin pathway (108). The increase in branching observed in mouse models of EDAR V370A (30, 31, 44) is likely the result of the EDA/EDAR/NF-ĸB pathway’s influence on the proteins that modulate the strength of the junctions that adhere the myoepithelial cells of the mammary ducts (109⇓–111) (Fig. 2 C–E), facilitating or hindering branching during morphogenesis (112). The more elaborate ductal branching induced by EDAR V370A during the embryonic stage of mammary development likely enhances the physiological effect of the increased branching induced by vitamin D deficiency in the later stages of breast development (during puberty and gestation/lactation).

Conclusions

The Arctic is an extreme environment for humans because of the almost complete lack of UV-B exposure that is required for the biosynthesis of vitamin D. This lack of vitamin D would have reduced immunological function and bone development and hindered the healthy function of adipose tissue. Here, we have presented a diverse array of evidence supporting the hypothesis that a genetically isolated population living in this environment during the LGM ∼20,000 y ago experienced selection for polymorphisms in the FADS gene cluster and for EDAR V370A because of the advantage these genetic variants likely confer in transmitting nutrients from mother to infant through breast milk under conditions of extremely low UV.

Acknowledgments

The synthesis presented here is a product of the American Association for the Advancement of Science’s 2017 annual meeting, as the ideas came together in the session entitled “Beringia and the Dispersal of Modern Humans to the Americas.” We thank Leslie Aiello, Marianne Brasil, Sarah Greenlee, Ophir Klein, Peter Kloess, Owen Lovejoy, Kunxin Luo, Whitney Reiner, Cat Taylor, Tim White, and three anonymous reviewers for critical feedback on earlier versions of this manuscript. Development of the ideas presented herein also benefited from conversations with Christopher Barrett, Michael Bell, Angelique Corthals, Liliana Davalos-Alvarez, Robert Dudley, Jennifer Raff, Ellen Simms, and Justin Tackney. L.J.H. also thanks the Department of Ecology and Evolution at the State University of New York at Stony Brook for feedback when she shared this research in their departmental colloquium series. Many thanks also to Christy G. Turner II, as his extensive study of incisor shoveling variation provided the foundation for the dental anthropology presented herein.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: hlusko{at}berkeley.edu.

  • Author contributions: L.J.H. designed research; L.J.H. and G.R.S. performed research; M.H. and M.A.P. contributed new reagents/analytic tools; L.J.H., G.C., T.A.M., and G.R.S. analyzed data; and L.J.H., J.P.C., G.C., S.A.E., J.F.H., N.G.J., T.A.M., D.H.O., M.A.P., and G.R.S. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1711788115/-/DCSupplemental.

References

  1. ↵
    1. Trigger BG
    (1989) A History of Archaeological Thought (Cambridge Univ Press, New York).
    .
  2. ↵
    1. Tamm E, et al.
    (2007) Beringian standstill and spread of Native American founders. PLoS One 2:e829.
    .
    OpenUrlCrossRefPubMed
  3. ↵
    1. Raghavan M, et al.
    (2015) POPULATION GENETICS. Genomic evidence for the Pleistocene and recent population history of Native Americans. Science 349:aab3884.
    .
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Reich D, et al.
    (2012) Reconstructing Native American population history. Nature 488:370–374.
    .
    OpenUrlCrossRefPubMed
  5. ↵
    1. Tackney JC, et al.
    (2015) Two contemporaneous mitogenomes from terminal Pleistocene burials in eastern Beringia. Proc Natl Acad Sci USA 112:13833–13838.
    .
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Moreno-Mayar JV, et al.
    (2018) Terminal Pleistocene Alaskan genome reveals first founding population of Native Americans. Nature 553:203–207.
    .
    OpenUrl
  7. ↵
    1. Hoffecker JF,
    2. Elias SA,
    3. O’Rourke DH,
    4. Scott GR,
    5. Bigelow NH
    (2016) Beringia and the global dispersal of modern humans. Evol Anthropol 25:64–78.
    .
    OpenUrl
  8. ↵
    1. Björn LO
    (2015) Photobiology: The Science of Light and Life (Springer, New York).
    .
  9. ↵
    1. Pibernat RA,
    2. Ellis-Evans C,
    3. Hinghofer-Szalkay HG
    (2007) Life in Extreme Environments (Springer Science & Business Media, Dordrecht, The Netherlands).
    .
  10. ↵
    1. Beall CM
    (2013) Human adaptability studies at high altitude: Research designs and major concepts during fifty years of discovery. Am J Hum Biol 25:141–147.
    .
    OpenUrlCrossRefPubMed
  11. ↵
    1. Azad P, et al.
    (2017) High-altitude adaptation in humans: From genomics to integrative physiology. J Mol Med (Berl) 95:1269–1282.
    .
    OpenUrl
  12. ↵
    1. Cho JI, et al.
    (2017) Ethnically Tibetan women in Nepal with low hemoglobin concentration have better reproductive outcomes. Evol Med Public Health 2017:82–96.
    .
    OpenUrl
  13. ↵
    1. Amorim CEG, et al.
    (2017) Genetic signature of natural selection in first Americans. Proc Natl Acad Sci USA 114:2195–2199.
    .
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Mathias RA,
    2. Pani V,
    3. Chilton FH
    (2014) Genetic variants in the FADS gene: Implications for dietary recommendations for fatty acid intake. Curr Nutr Rep 3:139–148.
    .
    OpenUrlCrossRefPubMed
  15. ↵
    1. Buckley MT, et al.
    (2017) Selection in Europeans on fatty acid desaturases associated with dietary changes. Mol Biol Evol 34:1307–1318.
    .
    OpenUrl
  16. ↵
    1. Fumagalli M, et al.
    (2015) Greenlandic Inuit show genetic signatures of diet and climate adaptation. Science 349:1343–1347.
    .
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Chilton FH, et al.
    (2014) Diet-gene interactions and PUFA metabolism: A potential contributor to health disparities and human diseases. Nutrients 6:1993–2022.
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Oscarsson J,
    2. Hurt-Camejo E
    (2017) Omega-3 fatty acids eicosapentaenoic acid and docosahexaenoic acid and their mechanisms of action on apolipoprotein B-containing lipoproteins in humans: A review. Lipids Health Dis 16:149.
    .
    OpenUrl
  19. ↵
    1. Zhang JY,
    2. Kothapalli KSD,
    3. Brenna JT
    (2016) Desaturase and elongase-limiting endogenous long-chain polyunsaturated fatty acid biosynthesis. Curr Opin Clin Nutr Metab Care 19:103–110.
    .
    OpenUrl
  20. ↵
    1. Williams GC
    (1957) Pleiotropy, natural selection, and the evolution of senescence. Evolution 11:398–411.
    .
    OpenUrlCrossRef
  21. ↵
    1. Haggarty P
    (2010) Fatty acid supply to the human fetus. Annu Rev Nutr 30:237–255.
    .
    OpenUrlCrossRefPubMed
  22. ↵
    1. Moltó-Puigmartí C, et al.
    (2010) FADS1 FADS2 gene variants modify the association between fish intake and the docosahexaenoic acid proportions in human milk. Am J Clin Nutr 91:1368–1376.
    .
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Lattka E, et al.
    (2011) Genetic variants in the FADS gene cluster are associated with arachidonic acid concentrations of human breast milk at 1.5 and 6 mo postpartum and influence the course of milk dodecanoic, tetracosenoic, and trans-9-octadecenoic acid concentrations over the duration of lactation. Am J Clin Nutr 93:382–391.
    .
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Moltó-Puigmartí C, et al.
    (2014) Maternal but not fetal FADS gene variants modify the association between maternal long-chain PUFA intake in pregnancy and birth weight. J Nutr 144:1430–1437.
    .
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Mennitti LV, et al.
    (2015) Type of fatty acids in maternal diets during pregnancy and/or lactation and metabolic consequences of the offspring. J Nutr Biochem 26:99–111.
    .
    OpenUrlCrossRefPubMed
  26. ↵
    1. Ding Z, et al.
    (2016) Association of polyunsaturated fatty acids in breast milk with fatty acid desaturase gene polymorphisms among Chinese lactating mothers. Prostaglandins Leukot Essent Fatty Acids 109:66–71.
    .
    OpenUrl
  27. ↵
    1. Sosa-Castillo E,
    2. Rodríguez-Cruz M,
    3. Moltó-Puigmartí C
    (2017) Genomics of lactation: Role of nutrigenomics and nutrigenetics in the fatty acid composition of human milk. Br J Nutr 118:161–168.
    .
    OpenUrl
  28. ↵
    1. Ibeagha-Awemu EM,
    2. Akwanji KA,
    3. Beaudoin F,
    4. Zhao X
    (2014) Associations between variants of FADS genes and omega-3 and omega-6 milk fatty acids of Canadian Holstein cows. BMC Genet 15:25.
    .
    OpenUrlPubMed
  29. ↵
    1. Bryk J, et al.
    (2008) Positive selection in East Asians for an EDAR allele that enhances NF-kappaB activation. PLoS One 3:e2209.
    .
    OpenUrlCrossRefPubMed
  30. ↵
    1. Chang SH,
    2. Jobling S,
    3. Brennan K,
    4. Headon DJ
    (2009) Enhanced Edar signalling has pleiotropic effects on craniofacial and cutaneous glands. PLoS One 4:e7591.
    .
    OpenUrlCrossRefPubMed
  31. ↵
    1. Kamberov YG, et al.
    (2013) Modeling recent human evolution in mice by expression of a selected EDAR variant. Cell 152:691–702.
    .
    OpenUrlCrossRefPubMed
  32. ↵
    1. Kowalczyk-Quintas C,
    2. Schneider P
    (2014) Ectodysplasin A (EDA) - EDA receptor signalling and its pharmacological modulation. Cytokine Growth Factor Rev 25:195–203.
    .
    OpenUrlCrossRefPubMed
  33. ↵
    1. Sadier A,
    2. Viriot L,
    3. Pantalacci S,
    4. Laudet V
    (2014) The ectodysplasin pathway: From diseases to adaptations. Trends Genet 30:24–31.
    .
    OpenUrlCrossRefPubMed
  34. ↵
    1. Glazer AM,
    2. Cleves PA,
    3. Erickson PA,
    4. Lam AY,
    5. Miller CT
    (2014) Parallel developmental genetic features underlie stickleback gill raker evolution. Evodevo 5:19.
    .
    OpenUrlCrossRefPubMed
  35. ↵
    1. Lefebvre S,
    2. Mikkola ML
    (2014) Ectodysplasin research–Where to next? Semin Immunol 26:220–228.
    .
    OpenUrlCrossRef
  36. ↵
    1. Goodwin AF, et al.
    (2014) Craniofacial morphometric analysis of individuals with X-linked hypohidrotic ectodermal dysplasia. Mol Genet Genomic Med 2:422–429.
    .
    OpenUrlCrossRefPubMed
  37. ↵
    1. Lindfors PH,
    2. Voutilainen M,
    3. Mikkola ML
    (2013) Ectodysplasin/NF-κB signaling in embryonic mammary gland development. J Mammary Gland Biol Neoplasia 18:165–169.
    .
    OpenUrlCrossRefPubMed
  38. ↵
    1. Fujimoto A, et al.
    (2008) A replication study confirmed the EDAR gene to be a major contributor to population differentiation regarding head hair thickness in Asia. Hum Genet 124:179–185.
    .
    OpenUrlCrossRefPubMed
  39. ↵
    1. Tan J, et al.
    (2013) The adaptive variant EDARV370A is associated with straight hair in East Asians. Hum Genet 132:1187–1191.
    .
    OpenUrlCrossRef
  40. ↵
    1. Peng Q, et al.
    (2016) EDARV370A associated facial characteristics in Uyghur population revealing further pleiotropic effects. Hum Genet 135:99–108.
    .
    OpenUrl
  41. ↵
    1. Park JH, et al.
    (2012) Effects of an Asian-specific nonsynonymous EDAR variant on multiple dental traits. J Hum Genet 57:508–514.
    .
    OpenUrlCrossRefPubMed
  42. ↵
    1. Tan J, et al.
    (2014) Characteristics of dental morphology in the Xinjiang Uyghurs and correlation with the EDARV370A variant. Sci China Life Sci 57:510–518.
    .
    OpenUrl
  43. ↵
    1. Kimura R, et al.
    (2009) A common variation in EDAR is a genetic determinant of shovel-shaped incisors. Am J Hum Genet 85:528–535.
    .
    OpenUrlCrossRefPubMed
  44. ↵
    1. Mou C, et al.
    (2008) Enhanced ectodysplasin-A receptor (EDAR) signaling alters multiple fiber characteristics to produce the East Asian hair form. Hum Mutat 29:1405–1411.
    .
    OpenUrlCrossRefPubMed
  45. ↵
    1. Sabeti PC, et al., International HapMap Consortium
    (2007) Genome-wide detection and characterization of positive selection in human populations. Nature 449:913–918.
    .
    OpenUrlCrossRefPubMed
  46. ↵
    1. Sicoli MA,
    2. Holton G
    (2014) Linguistic phylogenies support back-migration from Beringia to Asia. PLoS One 9:e91722.
    .
    OpenUrl
  47. ↵
    1. Mizoguchi Y
    (1985) Shovelling: A statistical analysis of its morphology. Univ Tokyo Bull 26:1–176.
    .
    OpenUrl
  48. ↵
    1. Scott GR,
    2. Turner CG II
    (1997) The Anthropology of Modern Human Teeth: Dental Morphology and Its Variation in Recent Human Populations (Cambridge Univ Press, New York).
    .
  49. ↵
    1. Llamas B, et al.
    (2016) Ancient mitochondrial DNA provides high-resolution time scale of the peopling of the Americas. Sci Adv 2:e1501385.
    .
    OpenUrlFREE Full Text
  50. ↵
    1. Fu Q, et al.
    (2014) Genome sequence of a 45,000-year-old modern human from western Siberia. Nature 514:445–449.
    .
    OpenUrlCrossRefPubMed
  51. ↵
    1. Pitulko V, et al.
    (2013) Human habitation in arctic western Beringia prior to the LGM. Paleoamerican Odyssey, eds Graf KE, Ketron CV, Waters MR (Texas A&M Univ Press, College Station, TX), pp 13–44.
    .
  52. ↵
    1. Claussen M,
    2. Selent K,
    3. Brovkin V,
    4. Raddatz T,
    5. Gayler V
    (2013) Impact of CO2 and climate on Last Glacial Maximum vegetation–A factor separation. Biogeosciences 10:3593–3604.
    .
    OpenUrl
  53. ↵
    1. Hewitt GM
    (2004) Genetic consequences of climatic oscillations in the Quaternary. Philos Trans R Soc Lond B Biol Sci 359:183–195, discussion 195.
    .
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Provan J,
    2. Bennett KD
    (2008) Phylogeographic insights into cryptic glacial refugia. Trends Ecol Evol 23:564–571.
    .
    OpenUrlCrossRefPubMed
  55. ↵
    1. Cai X, et al., Genographic Consortium
    (2011) Human migration through bottlenecks from Southeast Asia into East Asia during Last Glacial Maximum revealed by Y chromosomes. PLoS One 6:e24282.
    .
    OpenUrlCrossRefPubMed
  56. ↵
    1. Williams AN,
    2. Ulm S,
    3. Cook AR,
    4. Langley MC,
    5. Collard M
    (2013) Human refugia in Australia during the Last Glacial Maximum and terminal Pleistocene: A geospatial analysis of the 25–12 ka Australian archaeological record. J Archaeol Sci 40:4612–4625.
    .
    OpenUrlCrossRef
  57. ↵
    1. Hoffecker JF,
    2. Elias SA
    (2007) Human Ecology of Beringia (Columbia Univ Press, New York).
    .
  58. ↵
    1. Mulligan CJ,
    2. Kitchen A
    (2013) Three-stage colonization model for the peopling of the Americas. Paleoamerican Odyssey, eds Graf KE, Ketron CV, Waters MR (Texas A&M Univ Press, College Station, TX), pp 171–181.
    .
  59. ↵
    1. Hartl DL
    (1988) A Primer of Population Genetics (Sinauer, Sunderland, MA), 2nd Ed.
    .
  60. ↵
    1. Zanolli C,
    2. Hourset M,
    3. Esclassan R,
    4. Mollereau C
    (2017) Neanderthal and Denisova tooth protein variants in present-day humans. PLoS One 12:e0183802.
    .
    OpenUrl
  61. ↵
    1. Holick MF, et al.
    (1980) Photosynthesis of previtamin D3 in human skin and the physiologic consequences. Science 210:203–205.
    .
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Webb AR,
    2. Kline L,
    3. Holick MF
    (1988) Influence of season and latitude on the cutaneous synthesis of vitamin D3: Exposure to winter sunlight in Boston and Edmonton will not promote vitamin D3 synthesis in human skin. J Clin Endocrinol Metab 67:373–378.
    .
    OpenUrlCrossRefPubMed
  63. ↵
    1. Abbas MA
    (2017) Physiological functions of vitamin D in adipose tissue. J Steroid Biochem Mol Biol 165:369–381.
    .
    OpenUrlCrossRef
  64. ↵
    1. Lips P
    (2006) Vitamin D physiology. Prog Biophys Mol Biol 92:4–8.
    .
    OpenUrlCrossRefPubMed
  65. ↵
    1. Carlberg C,
    2. Seuter S,
    3. Heikkinen S
    (2012) The first genome-wide view of vitamin D receptor locations and their mechanistic implications. Anticancer Res 32:271–282.
    .
    OpenUrlAbstract/FREE Full Text
  66. ↵
    1. Ramagopalan SV, et al.
    (2010) A ChIP-seq defined genome-wide map of vitamin D receptor binding: Associations with disease and evolution. Genome Res 20:1352–1360.
    .
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Hart PH,
    2. Gorman S,
    3. Finlay-Jones JJ
    (2011) Modulation of the immune system by UV radiation: More than just the effects of vitamin D? Nat Rev Immunol 11:584–596.
    .
    OpenUrlCrossRefPubMed
  68. ↵
    1. Liu NQ,
    2. Hewison M
    (2012) Vitamin D, the placenta and pregnancy. Arch Biochem Biophys 523:37–47.
    .
    OpenUrlCrossRefPubMed
  69. ↵
    1. Workalemahu T, et al.
    (2017) Placental genetic variations in vitamin D metabolism and birthweight. Placenta 50:78–83.
    .
    OpenUrl
  70. ↵
    1. Mutt SJ,
    2. Hyppönen E,
    3. Saarnio J,
    4. Järvelin MR,
    5. Herzig KH
    (2014) Vitamin D and adipose tissue-more than storage. Front Physiol 5:228.
    .
    OpenUrlPubMed
  71. ↵
    1. Giordano A,
    2. Smorlesi A,
    3. Frontini A,
    4. Barbatelli G,
    5. Cinti S
    (2014) White, brown and pink adipocytes: The extraordinary plasticity of the adipose organ. Eur J Endocrinol 170:R159–R171.
    .
    OpenUrlAbstract/FREE Full Text
  72. ↵
    1. Cinti S
    (2017) UCP1 protein: The molecular hub of adipose organ plasticity. Biochimie 134:71–76.
    .
    OpenUrl
  73. ↵
    1. Valencak TG,
    2. Osterrieder A,
    3. Schulz TJ
    (2017) Sex matters: The effects of biological sex on adipose tissue biology and energy metabolism. Redox Biol 12:806–813.
    .
    OpenUrl
  74. ↵
    1. Harris RS,
    2. Bunker JW
    (1939) Vitamin D potency of human breast milk. Am J Public Health Nations Health 29:744–747.
    .
    OpenUrlCrossRefPubMed
  75. ↵
    1. Lee JH,
    2. O’Keefe JH,
    3. Lavie CJ,
    4. Harris WS
    (2009) Omega-3 fatty acids: Cardiovascular benefits, sources and sustainability. Nat Rev Cardiol 6:753–758.
    .
    OpenUrlCrossRefPubMed
  76. ↵
    1. Linday LA
    (2010) Cod liver oil, young children, and upper respiratory tract infections. J Am Coll Nutr 29:559–562.
    .
    OpenUrlPubMed
  77. ↵
    1. Challoumas D
    (2014) Vitamin D supplementation and lipid profile: What does the best available evidence show? Atherosclerosis 235:130–139.
    .
    OpenUrlCrossRefPubMed
  78. ↵
    1. Norton HL, et al.
    (2007) Genetic evidence for the convergent evolution of light skin in Europeans and East Asians. Mol Biol Evol 24:710–722.
    .
    OpenUrlCrossRefPubMed
  79. ↵
    1. Lamason RL, et al.
    (2005) SLC24A5, a putative cation exchanger, affects pigmentation in zebrafish and humans. Science 310:1782–1786.
    .
    OpenUrlAbstract/FREE Full Text
  80. ↵
    1. Lalueza-Fox C, et al.
    (2007) A melanocortin 1 receptor allele suggests varying pigmentation among Neanderthals. Science 318:1453–1455.
    .
    OpenUrlAbstract/FREE Full Text
  81. ↵
    1. Nielsen R, et al.
    (2017) Tracing the peopling of the world through genomics. Nature 541:302–310.
    .
    OpenUrlCrossRefPubMed
  82. ↵
    1. Fan S,
    2. Hansen ME,
    3. Lo Y,
    4. Tishkoff SA
    (2016) Going global by adapting local: A review of recent human adaptation. Science 354:54–59.
    .
    OpenUrlAbstract/FREE Full Text
  83. ↵
    1. Lachance J,
    2. Tishkoff SA
    (2013) Population genomics of human adaptation. Annu Rev Ecol Evol Syst 44:123–143.
    .
    OpenUrlCrossRefPubMed
  84. ↵
    1. Jablonski NG,
    2. Chaplin G
    (2010) Colloquium paper: Human skin pigmentation as an adaptation to UV radiation. Proc Natl Acad Sci USA 107:8962–8968.
    .
    OpenUrlAbstract/FREE Full Text
  85. ↵
    1. Chaplin G,
    2. Jablonski NG
    (2013) The human environment and the vitamin D compromise: Scotland as a case study in human biocultural adaptation and disease susceptibility. Hum Biol 85:529–552.
    .
    OpenUrlCrossRefPubMed
  86. ↵
    1. Chen TC, et al.
    (2007) Factors that influence the cutaneous synthesis and dietary sources of vitamin D. Arch Biochem Biophys 460:213–217.
    .
    OpenUrlCrossRefPubMed
  87. ↵
    1. Richards MP,
    2. Pettitt PB,
    3. Stiner MC,
    4. Trinkaus E
    (2001) Stable isotope evidence for increasing dietary breadth in the European mid-Upper Paleolithic. Proc Natl Acad Sci USA 98:6528–6532.
    .
    OpenUrlAbstract/FREE Full Text
  88. ↵
    1. Thacher TD,
    2. Fischer PR,
    3. Strand MA,
    4. Pettifor JM
    (2006) Nutritional rickets around the world: Causes and future directions. Ann Trop Paediatr 26:1–16.
    .
    OpenUrlCrossRefPubMed
  89. ↵
    1. Valentine CJ,
    2. Wagner CL
    (2013) Nutritional management of the breastfeeding dyad. Pediatr Clin North Am 60:261–274.
    .
    OpenUrlCrossRefPubMed
  90. ↵
    1. Picciano MF
    (2001) Nutrient composition of human milk. Pediatr Clin North Am 48:53–67.
    .
    OpenUrlCrossRefPubMed
  91. ↵
    1. Andreas NJ,
    2. Kampmann B,
    3. Mehring Le-Doare K
    (2015) Human breast milk: A review on its composition and bioactivity. Early Hum Dev 91:629–635.
    .
    OpenUrlCrossRefPubMed
  92. ↵
    1. Kasalová E,
    2. Aufartová J,
    3. Krčmová LK,
    4. Solichová D,
    5. Solich P
    (2015) Recent trends in the analysis of vitamin D and its metabolites in milk–A review. Food Chem 171:177–190.
    .
    OpenUrl
  93. ↵
    1. Specker BL,
    2. Tsang RC,
    3. Hollis BW
    (1985) Effect of race and diet on human-milk vitamin D and 25-hydroxyvitamin D. Am J Dis Child 139:1134–1137.
    .
    OpenUrlCrossRefPubMed
  94. ↵
    1. Lee JM, et al.
    (2007) Vitamin D deficiency in a healthy group of mothers and newborn infants. Clin Pediatr (Phila) 46:42–44.
    .
    OpenUrlCrossRefPubMed
  95. ↵
    1. Saraf R,
    2. Morton SM,
    3. Camargo CA Jr,
    4. Grant CC
    (2016) Global summary of maternal and newborn vitamin D status–A systematic review. Matern Child Nutr 12:647–668.
    .
    OpenUrl
  96. ↵
    1. Dawodu A, et al.
    (2014) Heightened attention to supplementation is needed to improve the vitamin D status of breastfeeding mothers and infants when sunshine exposure is restricted. Matern Child Nutr 10:383–397.
    .
    OpenUrlCrossRefPubMed
  97. ↵
    1. Greer FR,
    2. Hollis BW,
    3. Napoli JL
    (1984) High concentrations of vitamin D2 in human milk associated with pharmacologic doses of vitamin D2. J Pediatr 105:61–64.
    .
    OpenUrlCrossRefPubMed
  98. ↵
    1. Ching S,
    2. Kashinkunti S,
    3. Niehaus MD,
    4. Zinser GM
    (2011) Mammary adipocytes bioactivate 25-hydroxyvitamin D3 and signal via vitamin D3 receptor, modulating mammary epithelial cell growth. J Cell Biochem 112:3393–3405.
    .
    OpenUrlCrossRefPubMed
  99. ↵
    1. Matthews DG,
    2. D’Angelo J,
    3. Drelich J,
    4. Welsh J
    (2016) Adipose-specific Vdr deletion alters body fat and enhances mammary epithelial density. J Steroid Biochem Mol Biol 164:299–308.
    .
    OpenUrl
  100. ↵
    1. Paine IS,
    2. Lewis MT
    (2017) The terminal end bud: The little engine that could. J Mammary Gland Biol Neoplasia 22:93–108.
    .
    OpenUrl
  101. ↵
    1. Inman JL,
    2. Robertson C,
    3. Mott JD,
    4. Bissell MJ
    (2015) Mammary gland development: Cell fate specification, stem cells and the microenvironment. Development 142:1028–1042.
    .
    OpenUrlAbstract/FREE Full Text
  102. ↵
    1. Robinson GW
    (2004) Identification of signaling pathways in early mammary gland development by mouse genetics. Breast Cancer Res 6:105–108.
    .
    OpenUrlCrossRefPubMed
  103. ↵
    1. Zinser G,
    2. Packman K,
    3. Welsh J
    (2002) Vitamin D(3) receptor ablation alters mammary gland morphogenesis. Development 129:3067–3076.
    .
    OpenUrlAbstract/FREE Full Text
  104. ↵
    1. Qanbari S, et al.
    (2010) A genome-wide scan for signatures of recent selection in Holstein cattle. Anim Genet 41:377–389.
    .
    OpenUrlPubMed
  105. ↵
    1. Pan D, et al.
    (2013) Genome-wide detection of selective signature in Chinese Holstein. PLoS One 8:e60440.
    .
    OpenUrlCrossRefPubMed
  106. ↵
    1. Weir RR, et al.
    (2017) Environmental and genetic factors influence the vitamin D content of cows’ milk. Proc Nutr Soc 76:76–82.
    .
    OpenUrl
  107. ↵
    1. Cui X, et al.
    (2014) Transcriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing. BMC Genomics 15:226.
    .
    OpenUrl
  108. ↵
    1. Voutilainen M, et al.
    (2012) Ectodysplasin regulates hormone-independent mammary ductal morphogenesis via NF-κB. Proc Natl Acad Sci USA 109:5744–5749.
    .
    OpenUrlAbstract/FREE Full Text
  109. ↵
    1. Basham KJ, et al.
    (2013) Chemical genetic screen reveals a role for desmosomal adhesion in mammary branching morphogenesis. J Biol Chem 288:2261–2270.
    .
    OpenUrlAbstract/FREE Full Text
  110. ↵
    1. Tokonzaba E, et al.
    (2013) Plakoglobin as a regulator of desmocollin gene expression. J Invest Dermatol 133:2732–2740.
    .
    OpenUrlCrossRefPubMed
  111. ↵
    1. Celentano A,
    2. Mignogna MD,
    3. McCullough M,
    4. Cirillo N
    (2017) Pathophysiology of the desmo-adhesome. J Cell Physiol 232:496–505.
    .
    OpenUrl
  112. ↵
    1. Owens MB,
    2. Hill AD,
    3. Hopkins AM
    (2013) Ductal barriers in mammary epithelium. Tissue Barriers 1:e25933.
    .
    OpenUrl
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Environmental selection during the last ice age on the mother-to-infant transmission of vitamin D and fatty acids through breast milk
Leslea J. Hlusko, Joshua P. Carlson, George Chaplin, Scott A. Elias, John F. Hoffecker, Michaela Huffman, Nina G. Jablonski, Tesla A. Monson, Dennis H. O’Rourke, Marin A. Pilloud, G. Richard Scott
Proceedings of the National Academy of Sciences May 2018, 115 (19) E4426-E4432; DOI: 10.1073/pnas.1711788115
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Proceedings of the National Academy of Sciences: 116 (7)
Current Issue

Submit

Jump to section

You May Also be Interested in

Nonmonogamous strawberry poison frog (Oophaga pumilio). Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Putative signature of monogamy
A study suggests a putative gene-expression hallmark common to monogamous male vertebrates of some species, namely cichlid fishes, dendrobatid frogs, passeroid songbirds, common voles, and deer mice, and identifies 24 candidate genes potentially associated with monogamy.
Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Active lifestyles. Image courtesy of Pixabay/MabelAmber.
Meaningful life tied to healthy aging
Physical and social well-being in old age are linked to self-assessments of life worth, and a spectrum of behavioral, economic, health, and social variables may influence whether aging individuals believe they are leading meaningful lives.
Image courtesy of Pixabay/MabelAmber.