Population Trends Research Articles

Political Demography: The Political Consequences of Structural Population Change

This article surveys the growing field of political demography, which explores the political consequences of structural population change. It underscores the importance of integrating demography and political sociology research to better understand the complex and nuanced relationship between demography and political dynamics. The existing research demonstrates profound and multifaceted impacts of demographic shifts on the political landscape, with different demographic factors having distinct political consequences. Notably, population composition and distribution tend to hold greater political significance than sheer population size and growth. Furthermore, while more research is needed, the existing work suggests that the effect of structural demographic factors is neither inevitable nor without limit; rather, the political consequences of demographic change often exhibit nonlinear patterns and interact with prevailing socioeconomic and institutional contexts. As demographic shifts continue to unfold globally, political demography stands as a promising and enlightening area of research that merits further inquiry.

Molecular genetic diversity analysis of the Nigerian laughing dove (Streptopelia senegalensis) and related species using selected mitochondrial genes

Nigerian laughing doves (Streptopelia senegalensis) are small birds with long tail and living in bushes of the Sub-Saharan regions of African continent, the middle East and Asia, especially India. They are used for food, medicinal and religious purposes in Nigeria. Despite their usage, there is a lack of information on the genetic diversity of laughing doves in Nigeria. This study investigates taxonomic order and diversity of Nigerian laughing doves based on the mitochondrial cytochrome oxidase subunit I (COI) and cytochrome B (CYTB). The results showed 20 haplotypes within the 28 Nigerian coupled with Global Streptopelia genus using concatenated sequences. The Nigerian laughing dove constitute 16 distinct haplotypes. The haplotype diversity was 0.743 ± 0.070, and nucleotide diversity 0.154 ± 0.101 within Nigerian population using COI sequences. Phylogenetic tree showed that Nigerian laughing doves were in the same monophyletic clade with other Streptopela orientalis, S. decocto and S. chinensis; and this confirmed that Nigerian laughing doves might have shared descendant. The median-joining network further grouped Nigerian laughing doves into two: the first group consisting of Nigerian populations only, while the second group are with Saudi Arabian and Djiboutian populations. Population expansion was revealed in Nigerian dove individuals. This study revealed 16 unique haplotypes among Nigerian laughing dove population using concatenated sequences. Interestingly, CYTB showed clustering in African laughing doves (For instance, Nigerian individuals shared haplotypes with Sao Tome and Principe, an island country in the Gulf of Guinea, the western equatorial coast of Central Africa). The current data is the first report on genetic diversity of Nigerian laughing dove using mitochondrial COI and CYTB genes.

Reconstructing Impact of the 1867 Ionian Sea (Western Greece) Earthquake by Focusing on New Contemporary and Modern Sources for Building Damage, Environmental and Health Effects

The 4 February 1867 Cephalonia (Western Greece) earthquake is the largest in the Ionian Islands and one of the largest in the Eastern Mediterranean. However, it remained one of the least studied historical events. For reconstructing this earthquake, we reevaluated existing knowledge and used new contemporary and modern sources, including scientific and local writers’ reports and books, local and national journals, newspapers, and ecclesiastical chronicles. The extracted information covered the earthquake parameters, population impact, building damage, and earthquake environmental effects (EEEs). The earthquake parameters included the origin time and duration of the main shock, epicenter location, precursors, aftershocks, and characteristics of the earthquake ground motion. The population impact involved direct and indirect health effects and population change. Building data highlighted the dominant building types and the types, grades, and distribution of damage. The EEEs included ground cracks, landslides, liquefaction, hydrological anomalies, and mild sea disturbances. Field surveys were also conducted for validation. The quantitative and qualitative information enabled the application of seismic intensity scales (EMS-98, ESI-07). The study concluded that since the affected areas were mainly composed of post-alpine deposits and secondarily of clay–clastic alpine formations with poor geotechnical properties, they were highly susceptible to failure. Effects and maximum intensities occurred in highly susceptible areas with a rich inventory.

Modifying Factors of the Incidence of Type 1 Diabetes Mellitus in Children under 18 Years of Age in Northern Spain in the Last 3 Decades

Introduction: Type 1 Diabetes Mellitus is a chronic disease of autoimmune etiology that is highly prevalent in childhood. In the last decade, the possibility of the influence of the Bacillus Calmette-Guérin vaccine, different autoimmune diseases, modification of the immune response, and COVID19 itself of a simple population change due to immigration. Objectives: The objective of this work is to study the incidence of Type 1 Diabetes Mellitus in a cohort of children from Euskadi, north of Spain in the last 15 years. Methods: For the epidemiological study in Euskadi, a comparison of accumulated incidences between vaccinated and unvaccinated population groups was made. Results: The results of our study were not conclusive, but we believe that it is necessary to continue with long-term studies to eliminate bias.

Broken expectations: Population decline of a key grouper species within a 30-year-old no-take MPA in the Southwestern Atlantic

Understanding the impacts of human actions on marine species and developing mitigation strategies are crucial conservation challenges. Many fish species, especially marine top predators, are already highly depleted and rely on population monitoring and conservation strategies for their maintenance. In this context, Marine Protected Areas can help recover target species and provide crucial information about population trends when human activities are restricted or managed. Here, we used the longest and most complete temporal series of reef fish data in the Southwestern Atlantic, comprising 16 years of underwater visual census sampling. We evaluated temporal trends in the density and biomass of the dusky grouper (Epinephelus marginatus), both inside and outside a no-take and no-entry Marine Reserve. Our findings reveal a consistent decline in density and biomass inside and outside the MPA over time, reaching alarmingly low densities and experiencing a steeper biomass decline in the last five years. Our results contradict global expectations of positive population recovery trends within MPAs. We highlight the adverse impacts of human activities on the dusky grouper population, including illegal fishing and habitat degradation. Additionally, we address the limited understanding of species' life stages, underscoring that their occurrence outside protected boundaries may render the entire population vulnerable to human pressures. Thus, to guarantee population recovery, stakeholders and researchers must take immediate action to increase knowledge of all life stages and improve ongoing strategies.

Modeling carbon dioxide emissions reduction

The paper proposes a new scientific approach to modeling carbon dioxide (CO2) emissions in individual countries, based on the tools of Kohonen self-organizing maps and the provisions of economic theory. This work is based on conducted a comprehensive analysis of wide set of indicators influencing carbon dioxide emissions and decarbonization processes in 40 selected countries for 10 years time span from 2013 till 2022. During clustering these nations using self-organizing maps, 14 key indicators were chosen that cover economic and demographic growth, energy consumption, CO2 emission data, and includes trade of energy recourses for wider analysis and identification of countries with similar decarbonization potential, economic development and trade energy recourses possibilities. This approach exposed distinct clusters of nations with varying decarbonization capabilities and tracked the progress of analyzed countries in terms of CO2 emission changes in dynamic.The analysis of clustering was focused on identification leaders and followers in decarbonization activities. The leaders cluster includes Sweden, Norway, New Zealand and other countries exhibited leadership in developing a climate-resilient economy among 40 countries. Hence, all countries positioned closer to the leader cluster on the map demonstrated higher efficiency in their decarbonization pathways. Consequently, these properties of Kohonen maps provide a foundation for formulating general recommendations to achieve an efficient and effective low-carbon economy for followers including Ukraine, Kazakhstan and other energy intensive countries.Distinct countries decarbonization profiles resulting from clustering were used to develop practical recommendations towards estimating and establishing targeted CO2 emission levels for analyzed countries, conducting continued promoting renewable energy utilization, enhancing cross-border energy trade, diversifying energy portfolios, etc. Also, the proposed methodology can be utilized to project future emissions trends for each country based on cluster-specific models and facilitate policy decisions making process related to mitigating carbon emissions.

STING coordinates resolution of inflammation during wound repair by modulating macrophage trafficking through STAT3.

Efficient cutaneous wound healing requires a coordinated transition between inflammatory phases mediated by dynamic changes in leukocyte subset populations. Here, we identify STING as a key innate immune mediator governing timely resolution of inflammation by regulating macrophage dynamics during skin repair. Using a mouse model, we show STING deficiency caused delayed wound closure associated with abnormal persistence of TNF-α+ leukocytes. This resulted from the impaired macrophage recruitment. STING controlled the trafficking of bone marrow myeloid cells into blood and wounds, intrinsically enhancing macrophage migratory capacity through STAT3 activation. Specifically, STING modulated the production of monocyte chemokines and their receptors CCR2/CCR5 to enable efficient egress and wound infiltration. Consequently, disrupted systemic and local STING-STAT3-chemokine signaling combine to delay macrophage influx. This study elucidates STING as a critical rheostat tuning macrophage responses through STAT3 to orchestrate inflammatory resolution necessary for efficient wound healing. Our findings have broad implications for targeting STING therapeutically in both regenerative medicine and inflammatory disease contexts. STING regulates the macrophage trafficking through STAT3 in wound healing.

Genetic diversity of Dotilla myctiroides (Crustacea: Brachyura: Dotillidae) based on mitochondrial cox1 gene marker in Peninsular Malaysia

ABSTRACT The present study examined the genetic diversity of two populations of Dotilla myctiroides (Milne-Edwards) from Peninsular Malaysia. Based on 640-bp COX1 (cytochrome c oxidase subunit 1) gene sequences, the D. myctiroides specimens from Peninsular Malaysia, Singapore and India formed a lineage which was distinctly separated from Dotilla wichmanni and D. fenestrata. Sixteen COX1 haplotypes were found in Peninsular Malaysia samples, with 11 haplotypes in the Morib sample, eight in the Tanjung Tuan sample, and three common haplotypes in these samples. The Singapore specimen had the commonest haplotype of Peninsular Malaysia, while the Indian specimen had a distinct haplotype. The proportions of the two most common COX1 haplotypes were not significantly different between the Morib and Tanjung Tuan samples. The difference in the mean haplotype diversity between these samples was statistically significant. However, the nucleotide diversity of the two samples was not significantly different and there was no genetic subdivision/differentiation between the two populations. The high number of closely related haplotypes indicates the possibility of a recent population expansion. The data provide valuable information for future studies on the genetic diversity, population structure and phylogeography of Dotilla species.

The movement pattern changes of population following a disaster: Example of the Aegean Sea earthquake of October 2020

Unexpected seismic events can have a significant impact on our daily lives. Examining post-disaster population movement patterns contributes to disaster response and recovery operations. While most studies analyze population movement or displacement data, only a few examine multi-level movement patterns and consider public transportation demand. This paper aims to analyze changes in movement patterns following a major earthquake that severely affected İzmir in Türkiye, mainly using movement data from the Data for Good platform of Meta. According to the findings, on the first day following the earthquake, those traveling within İzmir considered the areas least damaged to be safer, despite the higher seismic intensity. The timeframe alignment between aftershock activity and the number of people traveling within İzmir suggests the routine movement pattern within the city was regained as the aftershocks decreased day by day. Unlike intra-city travel routines, there was an increase in the number of people moving to summer houses, located far from the center of İzmir. By November 2, 2020, three days after the earthquake, intra-city movements reached a new normal, as confirmed by public transportation use. However, intra-city travel patterns were impacted for approximately one month. This study reveals key findings from the evaluation of post-earthquake movements to improve disaster management decision making and aid crisis response and recovery operations.

Frequency-Dependent Vibronic Effects in Steady State Energy Transport.

The interplay between electronic and intramolecular high-frequency vibrational degrees of freedom is ubiquitous in natural light-harvesting systems. Recent studies have indicated that an intramolecular vibrational donor-acceptor frequency difference can enhance energy transport. Here, we analyze the extent to which different intramolecular donor-acceptor vibrational frequencies affect excitation energy transport in the natural nonequilibrium steady state configuration. Comments are included on the less physical equilibrium case for comparison with the literature. It is found that for constant Huang-Rhys factors, whereas the acceptor population increases in the equilibrium case when the intramolecular vibrational frequency of the acceptor exceeds that of the donor, this increase is negligible for the nonequilibrium steady state. Therefore, these changes in acceptor population do not significantly enhance energy transport in the nonequilibrium steady state for the natural scenario of incoherent light excitation with biologically relevant parameters of typical photosynthetic complexes. Insight about a potential mechanism to optimize energy transfer in the nonequilibrium steady state based on increasing the harvesting time at the reaction center is analyzed.

Rapid phenotypic and genotypic change in a laboratory schistosome population.

Genomic analysis has revealed extensive contamination among laboratory-maintained microbes including malaria parasites, Mycobacterium tuberculosis and Salmonella spp. Here, we provide direct evidence for recent contamination of a laboratory schistosome parasite population, and we investigate its genomic consequences. The Brazilian Schistosoma mansoni population SmBRE has several distinctive phenotypes, showing poor infectivity, reduced sporocysts number, low levels of cercarial shedding and low virulence in the intermediate snail host, and low worm burden and low fecundity in the vertebrate rodent host. In 2021 we observed a rapid change in SmBRE parasite phenotypes, with a ~10x increase in cercarial production and ~4x increase in worm burden. To determine the underlying genomic cause of these changes, we sequenced pools of SmBRE adults collected during parasite maintenance between 2015 and 2023. We also sequenced another parasite population (SmLE) maintained alongside SmBRE without phenotypic changes. While SmLE allele frequencies remained stable over the eight-year period, we observed sudden changes in allele frequency across the genome in SmBRE between July 2021 and February 2023, consistent with expectations of laboratory contamination. (i) SmLE-specific alleles rose in the SmBRE population from 0 to 41-46% across the genome between September and October 2021, documenting the timing and magnitude of the contamination event. (ii) After contamination, strong selection (s = ~0.23) drove replacement of low fitness SmBRE with high fitness SmLE alleles. (iii) Allele frequency changed rapidly across the whole genome, except for a region on chromosome 4 where SmBRE alleles remained at high frequency. We were able to detect contamination in this case because SmBRE shows distinctive phenotypes. However, this would likely have been missed with phenotypically similar parasites. These results provide a cautionary tale about the importance of tracking the identity of parasite populations, but also showcase a simple approach to monitor changes within populations using molecular profiling of pooled population samples to characterize fixed single nucleotide polymorphisms. We also show that genetic drift results in continuous change even in the absence of contamination, causing parasites maintained in different labs (or sampled from the same lab at different times) to diverge.

Treating gaps and biases in biodiversity data as a missing data problem.

Big biodiversity data sets have great potential for monitoring and research because of their large taxonomic, geographic and temporal scope. Such data sets have become especially important for assessing temporal changes in species' populations and distributions. Gaps in the available data, especially spatial and temporal gaps, often mean that the data are not representative of the target population. This hinders drawing large-scale inferences, such as about species' trends, and may lead to misplaced conservation action. Here, we conceptualise gaps in biodiversity monitoring data as a missing data problem, which provides a unifying framework for the challenges and potential solutions across different types of biodiversity data sets. We characterise the typical types of data gaps as different classes of missing data and then use missing data theory to explore the implications for questions about species' trends and factors affecting occurrences/abundances. By using this framework, we show that bias due to data gaps can arise when the factors affecting sampling and/or data availability overlap with those affecting species. But a data set per se is not biased. The outcome depends on the ecological question and statistical approach, which determine choices around which sources of variation are taken into account. We argue that typical approaches to long-term species trend modelling using monitoring data are especially susceptible to data gaps since such models do not tend to account for the factors driving missingness. To identify general solutions to this problem, we review empirical studies and use simulation studies to compare some of the most frequently employed approaches to deal with data gaps, including subsampling, weighting and imputation. All these methods have the potential to reduce bias but may come at the cost of increased uncertainty of parameter estimates. Weighting techniques are arguably the least used so far in ecology and have the potential to reduce both the bias and variance of parameter estimates. Regardless of the method, the ability to reduce bias critically depends on knowledge of, and the availability of data on, the factors creating data gaps. We use this review to outline the necessary considerations when dealing with data gaps at different stages of the data collection and analysis workflow.

20-Year trends in the social participation of the oldest old.

To investigate 20-year trends in social participation among the oldest old (77+ years) in Sweden and assess the extent to which changes in educational attainment and functional abilities explain these trends. Seven waves of the Swedish Panel Study of Living Conditions of the Oldest Old (SWEOLD) spanning 2002-2021 were used with a repeated cross-sectional design. To analyse the association between time and social participation we employed the Karlson-Holm-Breen method of decomposition. The study focused on informal social participation (contact with friends and family), leisure participation (public or semi-public gatherings), and formal participation (organisational engagement and study circle attendance). Both leisure and informal participation peaked in 2014 and declined in 2021, influenced by the COVID-19 pandemic, whereas formal participation showed a slight increase in 2021. Total participation increased at least until 2011. Overall, older adults have increased their levels of social participation in recent decades, disregarding the influence of the pandemic. Decomposition analysis revealed that population-level changes in educational attainment and functional abilities explained a substantial portion of the observed trends. As the proportion of older adults continues to rise, it becomes increasingly important to understand the developments and drivers of behavioural change in the older population. As more people are socially active, there may be increasing differences between those participating and those not - which could lead to increased inequalities. The observed trend in increasing participation, influenced by changes in education and health, emphasises the importance of fostering age-friendly environments and addressing potential social inequalities among older adults.

Immunodetection of selected pancreatic hormones under intragastric administration of apelin-13, a novel endogenous ligand for an angiotensin-like orphan G-protein coupled receptor, in unweaned rats.

This study investigated the effects of intragastric administration of apelin-13 on the secretion of critical pancreatic hormones in a cohort of three-week-old Wistar rats. The research aimed to uncover apelin's modulatory roles in endocrine interactions dictating metabolic homeostasis during early life. Rats were randomly assigned to control or experimental groups, receiving apelin-13 or saline for 14 days. The study population consisted of three-week-old Wistar rats of both sexes, weighing between 20 and 25 grams. Histological examination, analysis of variance and t-tests were employed to assess significant differences. Distinctive alterations in large islet morphology were observed, indicating a notable reduction in size. Additionally, an increase in alpha- and beta-cell density within specific islet sizes was noted, suggesting significant changes in cell populations. The study found a substantial increase in mitotic activity and a decrease in apoptosis in small and medium-sized islets post apelin-13 administration, indicating its potential role in regulating cell survival and proliferation. The notable reduction in large islet size coupled with increased alpha and beta cell density implies a targeted impact of apelin-13 on pancreatic cell dynamics. Also, the observed increase in mitotic activity and decrease in apoptosis in small and medium-sized islets suggest its potential regulatory role in cell survival and proliferation within the pancreatic microenvironment.

Place-based drivers and effective management of population growth and change in regional Australia

Place-based drivers and effective management of population growth and change in regional Australia

Trends in incidence of self-harm, neurodevelopmental and mental health conditions among university students compared with the general population: nationwide electronic data linkage study in Wales.

Concern that self-harm and mental health conditions are increasing in university students may reflect widening access to higher education, existing population trends and/or stressors associated with this setting. To compare population-level data on self-harm, neurodevelopmental and mental health conditions between university students and non-students with similar characteristics before and during enrolment. This cohort study linked electronic records from the Higher Education Statistics Agency for 2012-2018 to primary and secondary healthcare records. Students were undergraduates aged 18 to 24 years at university entry. Non-students were pseudo-randomly selected based on an equivalent age distribution. Logistic regressions were used to calculate odds ratios. Poisson regressions were used to calculate incidence rate ratios (IRR). The study included 96 760 students and 151 795 non-students. Being male, self-harm and mental health conditions recorded before university entry, and higher deprivation levels, resulted in lower odds of becoming a student and higher odds of drop-out from university. IRRs for self-harm, depression, anxiety, autism spectrum disorder (ASD), drug use and schizophrenia were lower for students. IRRs for self-harm, depression, attention-deficit hyperactivity disorder, ASD, alcohol use and schizophrenia increased more in students than in non-students over time. Older students experienced greater risk of self-harm and mental health conditions, whereas younger students were more at risk of alcohol use than non-student counterparts. Mental health conditions in students are common and diverse. While at university, students require person-centred stepped care, integrated with local third-sector and healthcare services to address specific conditions.

Loss of genetic variation and ancestral sex determination system in North American northern pike characterized by whole-genome resequencing.

The northern pike Esox lucius is a freshwater fish with low genetic diversity but ecological success throughout the Northern Hemisphere. Here, we generate an annotated chromosome-level genome assembly of 941 Mbp in length with 25 chromosome-length scaffolds. We then genotype 47 northern pike from Alaska through New Jersey at a genome-wide scale and characterize a striking decrease in genetic diversity along the sampling range. Individuals west of the North American Continental Divide have substantially higher diversity than those to the east (e.g. Interior Alaska and St. Lawrence River have on average 181 and 64K heterozygous SNPs per individual, or a heterozygous SNP every 5.2 and 14.6 kbp, respectively). Individuals clustered within each population with strong support, with numerous private alleles observed within each population. Evidence for recent population expansion was observed for a Manitoba hatchery and the St. Lawrence population (Tajima's D = -1.07 and -1.30, respectively). Several chromosomes have large regions with elevated diversity, including LG24, which holds amhby, the ancestral sex determining gene. As expected amhby was largely male-specific in Alaska and the Yukon and absent southeast to these populations, but we document some amhby(-) males in Alaska and amhby(+) males in the Columbia River, providing evidence for a patchwork of presence of this system in the western region. These results support the theory that northern pike recolonized North America from refugia in Alaska and expanded following deglaciation from west to east, with probable founder effects resulting in loss of both neutral and functional diversity (e.g. amhby).

Contrasting demographic history and mutational load in three threatened whitebark pines (Pinus subsect. Gerardianae): implications for conservation.

Demographic history and mutational load are of paramount importance for the adaptation of the endangered species. However, the effects of population evolutionary history and genetic load on the adaptive potential in endangered conifers remain unclear. Here, using population transcriptome sequencing, whole chloroplast genomes and mitochondrial DNA markers, combined with niche analysis, we determined the demographic history and mutational load for three threatened whitebark pines having different endangered statuses, Pinus bungeana, P. gerardiana and P. squamata. Demographic inference indicated that severe bottlenecks occurred in all three pines at different times, coinciding with periods of major climate and geological changes; in contrast, while P. bungeana experienced a recent population expansion, P. gerardiana and P. squamata maintained small population sizes after bottlenecking. Abundant homozygous-derived variants accumulated in the three pines, particularly in P. squamata, while the species with most heterozygous variants was P. gerardiana. Abundant moderately and few highly deleterious variants accumulated in the pine species that have experienced the most severe demographic bottlenecks (P. gerardiana and P. squamata), most likely because of purging effects. Finally, niche modeling showed that the distribution of P. bungeana might experience a significant expansion in the future, and the species' identified genetic clusters are also supported by differences in the ecological niche. The integration of genomic, demographic and niche data has allowed us to prove that the three threatened pines have contrasting patterns of demographic history and mutational load, which may have important implications in their adaptive potential and thus are also key for informing conservation planning.

Translocations, rising populations, and phylogeographic consequences: genomic implications for conservation of introduced Aoudad (Ammotragus lervia) in the southwestern United States

Abstract Over the last century, Texas has been inundated with exotic ungulates, with as many as 67 species represented by captive and free-ranging populations. One of the most successful species to have established free-ranging populations is the Aoudad (Ammotragus lervia). Herein, we expand on previous genetic studies in Aoudad by including new genomic profiles based on double-digest restriction site-associated DNA sequencing methods with the intent of better understanding the effects of >70 years of human-mediated translocations, population expansion, and gene flow among populations. Using a range of 4,338 to 5,529 single-nucleotide polymorphisms obtained from 73 individuals that potentially represented A. l. blainei, A. l. lervia, and A. l. sahariensis, we measured genetic patterns of diversity and connectivity between source populations from California and New Mexico and resultant aoudad populations in Texas. Principal component analyses identified phylogeographically structured populations across Texas. Further, we identified a signature of isolation by distance among Texas populations. To assess geographic regions that facilitated or blocked gene flow between populations, we used estimation of effective migration surfaces, which indicated 3 barriers to movement between populations as well as connectivity among populations in close proximity. Comparison of mtDNA haplogroups and nuclear genotypes among the 3 examined subspecies were discordant, indicating that the current subspecific taxonomic classification needs refinement. Given the vulnerable status of Aoudad in their native range of northern Africa and the exponential population growth and potential competition with native ungulates in Texas, we propose that efforts to manage non-native Aoudad should strive to maintain and preserve unique genetic units while minimizing competition and other detriments to native species in Texas. Conservation efforts are especially important given the vulnerable status of Aoudad within their native range of northern Africa.

A model of long-term population growth with an application to Central West Argentina.

We propose an Ideal Specialization Model to help explain the diversity of population growth trajectories exhibited across archaeological regions over thousands of years. The model provides a general set of expectations useful for guiding empirical research, and we provide a concrete example by conducting a preliminary evaluation of three expectations in Central West Argentina. We use kernel density estimates of archaeological radiocarbon, estimates of paleoclimate, and human bone stable isotopes from archaeological remains to evaluate three expectations drawn from the model's dynamics. Based on our results, we suggest that innovations in the production of food and social organization drove demographic transitions and population expansion in the region. The consistency of population expansion in the region positively associates with changes in diet and, potentially, innovations in settlement and social integration.