Archaeologists often argue whether Paleolithic works of art, cave paintings in particular, constitute reflections of the natural environment of humans at the time. They also debate the extent to which these paintings actually contain creative artistic expression, reflect the phenotypic variation of the surrounding environment, or focus on rare phenotypes. The famous paintings “The Dappled Horses of Pech-Merle,” depicting spotted horses on the walls of a cave in Pech-Merle, France, date back ∼25,000 y, but the coat pattern portrayed in these paintings is remarkably similar to a pattern known as “leopard” in modern horses. We have genotyped nine coat-color loci in 31 predomestic horses from Siberia, Eastern and Western Europe, and the Iberian Peninsula. Eighteen horses had bay coat color, seven were black, and six shared an allele associated with the leopard complex spotting ( LP ), representing the only spotted phenotype that has been discovered in wild, predomestic horses thus far. LP was detected in four Pleistocene and two Copper Age samples from Western and Eastern Europe, respectively. In contrast, this phenotype was absent from predomestic Siberian horses. Thus, all horse color phenotypes that seem to be distinguishable in cave paintings have now been found to exist in prehistoric horse populations, suggesting that cave paintings of this species represent remarkably realistic depictions of the animals shown. This finding lends support to hypotheses arguing that cave paintings might have contained less of a symbolic or transcendental connotation than often assumed.

With a range that covered most of northern Eurasia and parts of North America, the cave lion (Panthera spelaea) was one of the most widespread carnivores of the Late Pleistocene. Earlier ancient DNA analyses have shown that it is distinct from modern lions, and have suggested a demographic decline in Beringia during marine isotope stage 3 (MIS 3). Here, we further investigate the Late Pleistocene population dynamics in more detail by combining a powerful algorithm that couples MCMC with coalescent simulations under an approximate Bayesian computation framework. We use an ancient DNA dataset of previously published (n = 34) and new radiocarbon dated specimens (n = 14). Phylogenetic and network analyses based on the mitochondrial control region and the ATP8 gene identified two major haplogroups, one of which appears to vanish around 41,000 cal a BP. The approximate Bayesian computation analysis suggested a decline in effective population size (Ne) in Beringia of at least a 2-fold magnitude that began approximately 47,000 cal a BP, followed by an increase in Ne, most likely around 18,000 cal a BP. The cave lion went through a demographic bottleneck during MIS 3, which may have lasted for several tens of thousands of years, and only recovered shortly before the species' extinction. Several other large mammal species display similar declines in genetic diversity in Beringia during MIS 3, suggesting that major environmental changes might have affected megafaunal populations during this time period.

Leopard complex spotting is inherited by the incompletely dominant locus, LP , which also causes congenital stationary night blindness in homozygous horses. We investigated an associated single nucleotide polymorphism in the TRPM1 gene in 96 archaeological bones from 31 localities from Late Pleistocene (approx. 17 000 YBP) to medieval times. The first genetic evidence of LP spotting in Europe dates back to the Pleistocene. We tested for temporal changes in the LP associated allele frequency and estimated coefficients of selection by means of approximate Bayesian computation analyses. Our results show that at least some of the observed frequency changes are congruent with shifts in artificial selection pressure for the leopard complex spotting phenotype. In early domestic horses from Kirklareli–Kanligecit (Turkey) dating to 2700–2200 BC, a remarkably high number of leopard spotted horses (six of 10 individuals) was detected including one adult homozygote. However, LP seems to have largely disappeared during the late Bronze Age, suggesting selection against this phenotype in early domestic horses. During the Iron Age, LP reappeared, probably by reintroduction into the domestic gene pool from wild animals. This picture of alternating selective regimes might explain how genetic diversity was maintained in domestic animals despite selection for specific traits at different times.

The milpa is a traditional maize-based polyculture in Mexico that is typically practiced as rainfed agriculture. Because milpa cultivation has been practiced over a vast range of environmental and cultural conditions, this agroecosystem is recognized as an important repository of biological and cultural diversity. As for any agroecosystem, the relationship between plant development and the biogeochemical processes of the soil is critical. Although the milpa has been studied from different perspectives, the diversity and structure of microbial communities within milpa soils remain largely unexplored. In this study, we surveyed a milpa system in Central Mexico across cropping season: before planting (dry season; t1), during the early growth of plants (onset of the rainy season; t2), and before harvest (end of the rainy season; t3). In order to examine changes in community structure through time, we characterized bacterial diversity through high-throughput sequencing of 16S rRNA gene amplicons and recorded the nutrient status of multiple (5–10) soil samples from our milpa plots. We estimated microbial diversity from a total of 90 samples and constructed co-occurrence networks. Although we did not find significant changes in diversity or composition of bacterial communities across time, we identified significant rearrangements in their co-occurrence network structure. We found particularly drastic changes between the first and second time points. Co-occurrence analyses showed that the bacterial community changed from a less structured network at (t1) into modules with a non-random composition of taxonomic groups at (t2). We conclude that changes in bacterial communities undetected by standard diversity analyses can become evident when performing co-occurrence network analyses. We also postulate possible functional associations among keystone groups suggested by biogeochemical processes. This study represents the first contribution on soil microbial diversity of a maize-based polyculture and shows its dynamic nature in short-term scales.