Homo sapiens reached the higher latitudes of Europe by 45,000 years ago.

The Middle to Upper Palaeolithic transition in Europe is associated with the regional disappearance of Neanderthals and the spread of Homo sapiens. Late Neanderthals persisted in western Europe several millennia after the occurrence of H. sapiens in eastern Europe1. Local hybridization between the two groups occurred2, but not on all occasions3. Archaeological evidence also indicates the presence of several technocomplexes during this transition, complicating our understanding and the association of behavioural adaptations with specific hominin groups4. One such technocomplex for which the makers are unknown is the Lincombian-Ranisian-Jerzmanowician (LRJ), which has been described in northwestern and central Europe5-8. Here we present the morphological and proteomic taxonomic identification, mitochondrial DNA analysis and direct radiocarbon dating of human remains directly associated with an LRJ assemblage at the site Ilsenhöhle in Ranis (Germany). These human remains are among the earliest directly dated Upper Palaeolithic H. sapiens remains in Eurasia. We show that early H. sapiens associated with the LRJ were present in central and northwestern Europe long before the extinction of late Neanderthals in southwestern Europe. Our results strengthen the notion of a patchwork of distinct human populations and technocomplexes present in Europe during this transitional period.

A computational framework for ultra-high resolution cortical segmentation at 7Tesla.

This paper presents a computational framework for whole brain segmentation of 7Tesla magnetic resonance images able to handle ultra-high resolution data. The approach combines multi-object topology-preserving deformable models with shape and intensity atlases to encode prior anatomical knowledge in a computationally efficient algorithm. Experimental validation on simulated and real brain images shows accuracy and robustness of the method and demonstrates the benefits of an increased processing resolution.

A gradual increase of iron toward the medial-inferior tip of the subthalamic nucleus.

The subthalamic nucleus (STN) is an important node of the cortico-basal ganglia network and the main target of deep brain stimulation (DBS) in Parkinson's disease. Histological studies have revealed an inhomogeneous iron distribution within the STN, which has been related to putative subdivisions within this nucleus. Here, we investigate the iron distribution in more detail using quantitative susceptibility mapping (QSM), a novel magnetic resonance imaging (MRI) contrast mechanism. QSM allows for detailed assessment of iron content in both in vivo and postmortem tissue. Twelve human participants and 7 postmortem brain samples containing the STN were scanned using ultra-high field 7 Tesla (T) MRI. Iron concentrations were found to be higher in the medial-inferior tip of the STN. Using quantitative methods we show that the increase of iron concentration towards the medial-inferior tip is of a gradual rather than a discrete nature.

High-resolution MRI and diffusion-weighted imaging of the human habenula at 7 tesla.

PURPOSE: To investigate the feasibility of discriminating the habenula in human brain using high-resolution structural MRI and diffusion-weighted imaging at 7 Tesla (T). MATERIALS AND METHODS: MRI experiments included a MP2RAGE and GRE sequence to acquire quantitative parameter maps of T1, T2*, and a calculated proton density map and the combined approach of zoomed and parallel imaging (ZOOPPA) to obtain dw images. Probabilistic tractography algorithms were used to identify multiple fiber orientations in submillimetre voxels, and constrained spherical deconvolution to resolve orientations in regions where fibers cross. RESULTS: Maps of T1, T2*, and proton density showed high contrast of the human habenula. The lateral habenula and its commissure can be distinguished from medial habenula and adjacent tissue. DWI data with 0.7 mm isotropic resolution revealed that fiber populations differ in medial and lateral habenula and two major fiber bundles that connect habenular nuclei with forebrain structures and brainstem. CONCLUSION: High resolution 7T MR imaging of the human habenula provides sufficient signal-to-noise and contrast to enable identification of the lateral and medial nuclei. In vivo high resolution DWI at 7T is able to distinguish between lateral and medial habenula, and to detect major fiber tracts that connect the habenula with other brain areas.

The ecology, subsistence and diet of ~45,000-year-old Homo sapiens at Ilsenhöhle in Ranis, Germany.

Recent excavations at Ranis (Germany) identified an early dispersal of Homo sapiens into the higher latitudes of Europe by 45,000 years ago. Here we integrate results from zooarchaeology, palaeoproteomics, sediment DNA and stable isotopes to characterize the ecology, subsistence and diet of these early H. sapiens. We assessed all bone remains (n = 1,754) from the 2016-2022 excavations through morphology (n = 1,218) or palaeoproteomics (zooarchaeology by mass spectrometry (n = 536) and species by proteome investigation (n = 212)). Dominant taxa include reindeer, cave bear, woolly rhinoceros and horse, indicating cold climatic conditions. Numerous carnivore modifications, alongside sparse cut-marked and burnt bones, illustrate a predominant use of the site by hibernating cave bears and denning hyaenas, coupled with a fluctuating human presence. Faunal diversity and high carnivore input were further supported by ancient mammalian DNA recovered from 26 sediment samples. Bulk collagen carbon and nitrogen stable isotope data from 52 animal and 10 human remains confirm a cold steppe/tundra setting and indicate a homogenous human diet based on large terrestrial mammals. This lower-density archaeological signature matches other Lincombian-Ranisian-Jerzmanowician sites and is best explained by expedient visits of short duration by small, mobile groups of pioneer H. sapiens.

Stable isotopes show Homo sapiens dispersed into cold steppes ~45,000 years ago at Ilsenhöhle in Ranis, Germany.

The spread of Homo sapiens into new habitats across Eurasia ~45,000 years ago and the concurrent disappearance of Neanderthals represents a critical evolutionary turnover in our species' history. 'Transitional' technocomplexes, such as the Lincombian-Ranisian-Jerzmanowician (LRJ), characterize the European record during this period but their makers and evolutionary significance have long remained unclear. New evidence from Ilsenhöhle in Ranis, Germany, now provides a secure connection of the LRJ to H. sapiens remains dated to ~45,000 years ago, making it one of the earliest forays of our species to central Europe. Using many stable isotope records of climate produced from 16 serially sampled equid teeth spanning ~12,500 years of LRJ and Upper Palaeolithic human occupation at Ranis, we review the ability of early humans to adapt to different climate and habitat conditions. Results show that cold climates prevailed across LRJ occupations, with a temperature decrease culminating in a pronounced cold excursion at ~45,000-43,000 cal BP. Directly dated H. sapiens remains confirm that humans used the site even during this very cold phase. Together with recent evidence from the Initial Upper Palaeolithic, this demonstrates that humans operated in severe cold conditions during many distinct early dispersals into Europe and suggests pronounced adaptability.

Optimizing T1-weighted imaging of cortical myelin content at 3.0 T.

With increases in the sensitivity and resolution of anatomical MRI for the brain, methods for mapping the organization of the cerebral cortex by imaging its myelin content have emerged. This identifies major sensory and motor regions and could be used in studies of cortical organization, particularly if patterns of myelination can be visualized over the cortical surface robustly in individual subjects. The imaging problem is difficult, however, because of the relative thinness of the cerebral cortex and the low intracortical tissue contrast. In this paper, we optimize the contrast of T(1)-weighted MRI to help better visualize patterns of myelination. We measure a small but statistically significant difference in T(1) of 171 ± 40 ms between cortical regions with low and high myelin contents in the human cortex at 3T, and then perform simulations to choose parameters for an inversion-recovery pulse sequence that utilizes this T(1) difference to increase contrast within the cortex. We show that lengthening the delay between signal acquisition and the next inversion pulse in the sequence increases intracortical contrast more effectively than does image averaging. Using the optimized sequence, we show that major myelinated regions that are relatively thick, such as the primary motor and auditory regions, can be visualized well in individuals at 3T using whole-cortex 3D images made at 1mm isotropic resolution, while thinner regions, such as the primary visual cortex, can be visualized using targeted 3D images made at 0.5mm isotropic resolution. Our findings demonstrate that patterns of myelination can be better visualized in individual subjects when the imaging is optimized to highlight intracortical contrast and can help to pave the way for the creation of matched maps of microanatomy and function in the cortex of living individual humans.

Comparing extraction method efficiency for high-throughput palaeoproteomic bone species identification.

High-throughput proteomic analysis of archaeological skeletal remains provides information about past fauna community compositions and species dispersals in time and space. Archaeological skeletal remains are a finite resource, however, and therefore it becomes relevant to optimize methods of skeletal proteome extraction. Ancient proteins in bone specimens can be highly degraded and consequently, extraction methods for well-preserved or modern bone might be unsuitable for the processing of highly degraded skeletal proteomes. In this study, we compared six proteomic extraction methods on Late Pleistocene remains with variable levels of proteome preservation. We tested the accuracy of species identification, protein sequence coverage, deamidation, and the number of post-translational modifications per method. We find striking differences in obtained proteome complexity and sequence coverage, highlighting that simple acid-insoluble proteome extraction methods perform better in highly degraded contexts. For well-preserved specimens, the approach using EDTA demineralization and protease-mix proteolysis yielded a higher number of identified peptides. The protocols presented here allowed protein extraction from ancient bone with a minimum number of working steps and equipment and yielded protein extracts within three working days. We expect further development along this route to benefit large-scale screening applications of relevance to archaeological and human evolution research.

From genotypes to organisms: State-of-the-art and perspectives of a cornerstone in evolutionary dynamics.

Understanding how genotypes map onto phenotypes, fitness, and eventually organisms is arguably the next major missing piece in a fully predictive theory of evolution. We refer to this generally as the problem of the genotype-phenotype map. Though we are still far from achieving a complete picture of these relationships, our current understanding of simpler questions, such as the structure induced in the space of genotypes by sequences mapped to molecular structures, has revealed important facts that deeply affect the dynamical description of evolutionary processes. Empirical evidence supporting the fundamental relevance of features such as phenotypic bias is mounting as well, while the synthesis of conceptual and experimental progress leads to questioning current assumptions on the nature of evolutionary dynamics-cancer progression models or synthetic biology approaches being notable examples. This work delves with a critical and constructive attitude into our current knowledge of how genotypes map onto molecular phenotypes and organismal functions, and discusses theoretical and empirical avenues to broaden and improve this comprehension. As a final goal, this community should aim at deriving an updated picture of evolutionary processes soundly relying on the structural properties of genotype spaces, as revealed by modern techniques of molecular and functional analysis.

Diffusion tensor imaging segments the human amygdala in vivo.

The amygdala plays an important role in emotion, learning, and memory. It would be highly advantageous to understand more precisely its internal structure and connectivity for individual human subjects in vivo. Earlier cytoarchitectural research in post-mortem human and animal brains has revealed multiple subdivisions and connectivity patterns, probably related to different functions. With standard magnetic resonance imaging (MRI) techniques, however, the amygdala appears as an undifferentiated area of grey matter. Using high-quality diffusion tensor imaging (DTI) at 3 Tesla, we show diffusion anisotropy in this grey matter area. Such data allowed us to subdivide the amygdala for the first time in vivo. In 15 living subjects, we applied a spectral clustering algorithm to the principal diffusion direction in each amygdala voxel and found a consistent subdivision of the amygdala into a medial and a lateral region. The topography of these regions is in good agreement with the fibre architecture visible in myelin-stained sections through the amygdala of a human post-mortem brain. From these in vivo results we derived a probabilistic map of amygdalar fibre orientations. This segmentation technique has important implications for functional studies in the processing of emotions, cognitive function, and psychiatric disorders and in studying morphometry and volumetry of amygdala subdivisions.

The Lichtenberg Keilmesser - it's all about the angle.

The presence of the 'Keilmesser-concept' in late Middle Paleolithic assemblages of Central and Eastern Europe defines the eponymous 'Keilmessergruppen'. The site of Lichtenberg (Lower Saxony, Germany) was discovered in 1987 and yielded one of the most important Keilmessergruppen assemblages of the northwestern European Plain. At that time, researchers used the bifacial backed knives to define a new type, the 'Lichtenberger Keilmesser', which they characterized by an aesthetic form-function concept with a specific range of morphological variability on the one hand, and a standardized convex cutting edge one the other hand. Thereby, a shape continuum was observed between different form-function concepts in the Lichtenberg assemblage, from Keilmesser through to Faustkeilblätter and handaxes. In a contrasting view, it was recently suggested that the morphology of Keilmesser, including what is defined here as type Lichtenberg, is the result of solutions to establish and maintain edge angles during resharpening. With the intention to evaluate these contrasting hypotheses, I conducted a re-analysis of the Keilmesser from Lichtenberg and their relationship to central German late Middle Paleolithic knives, using 3D geometric morphometric analyses and an automatized approach to measure edge angles on 3D models. Despite a morphological overlap of the tools from both regions, I could show that the Lichtenberg Keilmesser concept refers to one solution to create a tool with specific functionalities, like potentially cutting, prehension, and reusability. To establish and maintain its functionality, certain angles where created by the knappers along the active edges. This behavior resulted in specific shapes and positions of the active parts and created what looks like a standardized or template morphology of this Keilmesser type.

Neutral components show a hierarchical community structure in the genotype-phenotype map of RNA secondary structure.

Genotype-phenotype (GP) maps describe the relationship between biological sequences and structural or functional outcomes. They can be represented as networks in which genotypes are the nodes, and one-point mutations between them are the edges. The genotypes that map to the same phenotype form subnetworks consisting of one or multiple disjoint connected components-so-called neutral components (NCs). For the GP map of RNA secondary structure, the NCs have been found to exhibit distinctive network features that can affect the dynamical processes taking place on them. Here, we focus on the community structure of RNA secondary structure NCs. Building on previous findings, we introduce a method to reveal the hierarchical community structure solely from the sequence constraints and composition of the genotypes that form a given NC. Thereby, we obtain modularity values similar to common community detection algorithms, which are much more complex. From this knowledge, we endorse a sampling method that allows a fast exploration of the different communities of a given NC. Furthermore, we introduce a way to estimate the community structure from genotype samples, which is useful when an exhaustive analysis of the NC is not feasible, as is the case for longer sequence lengths.

Using small samples to estimate neutral component size and robustness in the genotype-phenotype map of RNA secondary structure.

In genotype-phenotype (GP) maps, the genotypes that map to the same phenotype are usually not randomly distributed across the space of genotypes, but instead are predominantly connected through one-point mutations, forming network components that are commonly referred to as neutral components (NCs). Because of their impact on evolutionary processes, the characteristics of these NCs, like their size or robustness, have been studied extensively. Here, we introduce a framework that allows the estimation of NC size and robustness in the GP map of RNA secondary structure. The advantage of this framework is that it only requires small samples of genotypes and their local environment, which also allows experimental realizations. We verify our framework by applying it to the exhaustively analysable GP map of RNA sequence length L = 15, and benchmark it against an existing method by applying it to longer, naturally occurring functional non-coding RNA sequences. Although it is specific to the RNA secondary structure GP map in the first place, our framework can probably be transferred and adapted to other sequence-to-structure GP maps.

Introducing platform surface interior angle (PSIA) and its role in flake formation, size and shape.

Four ways archaeologists have tried to gain insights into how flintknapping creates lithic variability are fracture mechanics, controlled experimentation, replication and attribute studies of lithic assemblages. Fracture mechanics has the advantage of drawing more directly on first principles derived from physics and material sciences, but its relevance to controlled experimentation, replication and lithic studies more generally has been limited. Controlled experiments have the advantage of being able to isolate and quantify the contribution of individual variables to knapping outcomes, and the results of these experiments have provided models of flake formation that when applied to the archaeological record of flintknapping have provided insights into past behavior. Here we develop a linkage between fracture mechanics and the results of previous controlled experiments to increase their combined explanatory and predictive power. We do this by documenting the influence of Herztian cone formation, a constant in fracture mechanics, on flake platforms. We find that the platform width is a function of the Hertzian cone constant angle and the geometry of the platform edge. This finding strengthens the foundation of one of the more influential models emerging from the controlled experiments. With additional work, this should make it possible to merge more of the experimental results into a more comprehensive model of flake formation.

Timing of the Saalian- and Elsterian glacial cycles and the implications for Middle - Pleistocene hominin presence in central Europe.

By establishing a luminescence-based chronology for fluvial deposits preserved between the Elsterian- and Saalian tills in central Germany, we obtained information on the timing of both the Middle Pleistocene glacial cycles and early human appearance in central Europe. The luminescence ages illustrate different climatic driven fluvial aggradation periods during the Saalian glacial cycle spanning from 400-150 ka. The ages of sediments directly overlying the Elsterian till are approximately 400 ka and prove that the first extensive Fennoscandian ice sheet extension during the Quaternary correlates with MIS 12 and not with MIS 10. Furthermore, the 400 ka old fluvial units contain Lower Paleolithic stone artefacts that document the first human appearance in the region. In addition, we demonstrate that early MIS 8 is a potential date for the onset of the Middle Paleolithic in central Germany, as Middle Paleolithic stone artefacts are correlated with fluvial units deposited between 300 ka and 200 ka. However, the bulk of Middle Paleolithic sites date to MIS 7 in the region. The fluvial units preserved directly under the till of the southernmost Saalian ice yield an age of about 150 ka, and enable a correlation of the Drenthe stage to late MIS 6.

Phenotypes can be robust and evolvable if mutations have non-local effects on sequence constraints.

The mapping between biological genotypes and phenotypes plays an important role in evolution, and understanding the properties of this mapping is crucial to determine the outcome of evolutionary processes. One of the most striking properties observed in several genotype-phenotype (GP) maps is the positive correlation between the robustness and evolvability of phenotypes. This implies that a phenotype can be strongly robust against mutations and at the same time evolvable to a diverse range of alternative phenotypes. Here, we examine the causes for this positive correlation by introducing two analytically tractable GP map models that follow the principles of real biological GP maps. The first model is based on gene-like GP maps, reflecting the way in which genetic sequences are organized into protein-coding genes, and the second one is based on the GP map of RNA secondary structure. For both models, we find that a positive correlation between phenotype robustness and evolvability only emerges if mutations at one sequence position can have non-local effects on the sequence constraints at another position. This highlights that non-local effects of mutations are closely related to the coexistence of robustness and evolvability in phenotypes, and are likely to be an important feature of many biological GP maps.

Spatial normalization of ultrahigh resolution 7 T magnetic resonance imaging data of the postmortem human subthalamic nucleus: a multistage approach.

In this paper, we describe a novel processing strategy for the spatial normalization of ultrahigh resolution magnetic resonance imaging (MRI) data of small ex vivo samples into MNI standard space. We present a multistage scanning and registration method for data of the subthalamic nucleus (STN) obtained using ultrahigh 7 T MRI on four human postmortem brain samples. Four whole brains were obtained and subjected to multistage MRI scanning, corresponding to four different brain dissection stages. Data sets were acquired with an isotropic resolution of 100 µm enabling accurate manual segmentation of the STN. Spatial normalization to MNI reference space was performed, probability maps were calculated, and results were cross-checked with an independent in vivo dataset showing significant overlay. Normalization of results obtained from small tissue samples into MNI standard space will facilitate comparison between individual subjects, as well as between studies. When combining ultrahigh resolution MRI of ex vivo samples with histological studies via blockface imaging, our method enables further insight and inference as multimodal data can be compared within the same reference space. This novel technique may be of value for research purposes using functional MRI techniques, and in the future may be of assistance for anatomical orientation in clinical practice.

Mapping of the internal structure of human habenula with ex vivo MRI at 7T.

The habenula is a small but important nucleus located next to the third ventricle in front of the pineal body. It helps to control the human reward system and is considered to play a key role in emotion, showing increased activation in major depressive disorders. Its dysfunction may underlie several neurological and psychiatric disorders. It is now possible to visualize the habenula and its anatomical subdivisions-medial habenula (MHB) and lateral habenula (LHB)-using MR techniques. The aim of this study was to further differentiate substructures within human lateral habenula (LHB) using ex vivo ultra-high field MR structural imaging, distinguishing between a medial part (m-LHB) and a lateral part (l-LHB). High resolution T1w images with 0.3-mm isotropic resolution and T2(*)w images with 60-micrometer isotropic resolution were acquired on a 7T MR scanner and quantitative maps of T1 and T2(*) were calculated. Cluster analysis of image intensity was performed using the Fuzzy and Noise Tolerant Adaptive Segmentation Method (FANTASM) tool. Ultra-high resolution structural MRI of ex vivo brain tissue at 7T provided sufficient SNR and contrast to discriminate the medial and lateral habenular nuclei. Heterogeneity was observed in the lateral habenula (LHB) nuclei, with clear distinctions between lateral and medial parts (m-LHB, l-LHB) and with the neighboring medial habenula (MHB). Clustering analysis based on the T1 and T2(*) maps strongly showed 4-6 clusters as subcomponents of lateral and medial habenula.