Capturing the fusion of two ancestries and kinship structures in Merovingian Flanders.

The Merovingian period (5th to 8th cc AD) was a time of demographic, socioeconomic, cultural, and political realignment in Western Europe. Here, we report the whole-genome shotgun sequence data of 30 human skeletal remains from a coastal Late Merovingian site of Koksijde (675 to 750 AD), alongside 18 remains from two Early to Late Medieval sites in present-day Flanders, Belgium. We find two distinct ancestries, one shared with Early Medieval England and the Netherlands, while the other, minor component, reflecting likely continental Gaulish ancestry. Kinship analyses identified no large pedigrees characteristic to elite burials revealing instead a high modularity of distant relationships among individuals of the main ancestry group. In contrast, individuals with >90% Gaulish ancestry had no kinship links among sampled individuals. Evidence for population structure and major differences in the extent of Gaulish ancestry in the main group, including in a mother-daughter pair, suggests ongoing admixture in the community at the time of their burial. The isotopic and genetic evidence combined supports a model by which the burials, representing an established coastal nonelite community, had incorporated migrants from inland populations. The main group of burials at Koksijde shows an abundance of >5 cM long shared allelic intervals with the High Medieval site nearby, implying long-term continuity and suggesting that similarly to Britain, the Early Medieval ancestry shifts left a significant and long-lasting impact on the genetic makeup of the Flemish population. We find substantial allele frequency differences between the two ancestry groups in pigmentation and diet-associated variants, including those linked with lactase persistence, likely reflecting ancestry change rather than local adaptation.

Formation of human kinship structures depending on population size and cultural mutation rate.

How does social complexity depend on population size and cultural transmission? Kinship structures in traditional societies provide a fundamental illustration, where cultural rules between clans determine people's marriage possibilities. Here, we propose a simple model of kinship interactions that considers kin and in-law cooperation and sexual rivalry. In this model, multiple societies compete. Societies consist of multiple families with different cultural traits and mating preferences. These values determine interactions and hence the growth rate of families and are transmitted to offspring with mutations. Through a multilevel evolutionary simulation, family traits and preferences are grouped into multiple clans with interclan mating preferences. It illustrates the emergence of kinship structures as the spontaneous formation of interdependent cultural associations. Emergent kinship structures are characterized by the cycle length of marriage exchange and the number of cycles in society. We numerically and analytically clarify their parameter dependence. The relative importance of cooperation versus rivalry determines whether attraction or repulsion exists between families. Different structures evolve as locally stable attractors. The probabilities of formation and collapse of complex structures depend on the number of families and the mutation rate, showing characteristic scaling relationships. It is now possible to explore macroscopic kinship structures based on microscopic interactions, together with their environmental dependence and the historical causality of their evolution. We propose the basic causal mechanism of the formation of typical human social structures by referring to ethnographic observations and concepts from statistical physics and multilevel evolution. Such interdisciplinary collaboration will unveil universal features in human societies.

Mobility and kinship in the world's first village societies.

Around 10,000 y ago in southwest Asia, the cessation of a mobile lifestyle and the emergence of the first village communities during the Neolithic marked a fundamental change in human history. The first communities were small (tens to hundreds of individuals) but remained semisedentary. So-called megasites appeared soon after, occupied by thousands of more sedentary inhabitants. Accompanying this shift, the material culture and ancient ecological data indicate profound changes in economic and social behavior. A shift from residential to logistical mobility and increasing population size are clear and can be explained by either changes in fertility and/or aggregation of local groups. However, as sedentism increased, small early communities likely risked inbreeding without maintaining or establishing exogamous relationships typical of hunter-gatherers. Megasites, where large populations would have made endogamy sustainable, could have avoided this risk. To examine the role of kinship practices in the rise of megasites, we measured strontium and oxygen isotopes in tooth enamel from 99 individuals buried at Pinarbasi, Boncuklu, and Çatalhöyük (Turkey) over 7,000 y. These sites are geographically proximate and, critically, span both early sedentary behaviors (Pinarbasi and Boncuklu) and the rise of a local megasite (Çatalhöyük). Our data are consistent with the presence of only local individuals at Pinarbasi and Boncuklu, whereas at Çatalhöyük, several nonlocals are present. The Çatalhöyük data stand in contrast to other megasites where bioarchaeological evidence has pointed to strict endogamy. These different kinship behaviors suggest that megasites may have arisen by employing unique, community-specific kinship practices.

Descent, marriage, and residence practices of a 3,800-year-old pastoral community in Central Eurasia.

Our understanding of prehistoric societal organization at the family level is still limited. Here, we generated genome data from 32 individuals from an approximately 3,800-y-old burial mound attributed to the Bronze Age Srubnaya-Alakul cultural tradition at the site of Nepluyevsky, located in the Southern Ural region of Central Eurasia. We found that life expectancy was generally very low, with adult males living on average 8 y longer than females. A total of 35 first-degree, 40 second-degree, and 48 third-degree biological relationships connected 23 of the studied individuals, allowing us to propose a family tree spanning three generations with six brothers at its center. The oldest of these brothers had eight children with two women and the most children overall, whereas the other relationships were monogamous. Notably, related female children above the age of five were completely absent from the site, and adult females were more genetically diverse than males. These results suggest that biological relationships between male siblings played a structural role in society and that descent group membership was based on patrilineality. Women originated from a larger mating network and moved to join the men, with whom they were buried. Finally, the oldest brother likely held a higher social position, which was expressed in terms of fertility.

Projections of human kinship for all countries.

Demographers have long attempted to project future changes in the size and composition of populations, but have ignored what these processes will mean for the size, composition, and age distribution of family networks. Kinship structures matter because family solidarity-a crucial source of informal care for millions of people around the world-is conditional on kin being alive. Here, we present innovative projections of biological kin for the 1950 to 2100 period and discuss what they imply for the availability of informal care. Overall, we project that the number of living kin for individuals will decline dramatically worldwide. While a 65-yo woman in 1950 could expect to have 41 living kin, a 65-yo woman in 2095 is projected to have just 25 [18.8 to 34.7] relatives (lower and upper 80% projection intervals). This represents a 38% [15 to 54] global decline. The composition of family networks is also expected to change, with the numbers of living grandparents and great-grandparents markedly increasing, and the numbers of cousins, nieces and nephews, and grandchildren declining. Family networks will age considerably, as we project a widening age gap between individuals and their kin due to lower and later fertility and longer lifespans. In Italy, for example, the average age of a grandmother of a 35-yo woman is expected to increase from 77.9 y in 1950 to 87.7 y [87.1 to 88.5] in 2095. The projected changes in kin supply will put pressure on the already stretched institutional systems of social support, as more individuals age with smaller and older family networks.

Low Genetic Impact of the Roman Occupation of Britain in Rural Communities.

The Roman period saw the empire expand across Europe and the Mediterranean, including much of what is today Great Britain. While there is written evidence of high mobility into and out of Britain for administrators, traders, and the military, the impact of imperialism on local, rural population structure, kinship, and mobility is invisible in the textual record. The extent of genetic change that occurred in Britain during the Roman military occupation remains underexplored. Here, using genome-wide data from 52 ancient individuals from eight sites in Cambridgeshire covering the period of Roman occupation, we show low levels of genetic ancestry differentiation between Romano-British sites and indications of larger populations than in the Bronze Age and Neolithic. We find no evidence of long-distance migration from elsewhere in the Empire, though we do find one case of possible temporary mobility within a family unit during the Late Romano-British period. We also show that the present-day patterns of genetic ancestry composition in Britain emerged after the Roman period.

Measuring the Efficiency of Purging by non-random Mating in Human Populations.

Human populations harbor a high concentration of deleterious genetic variants. Here, we tested the hypothesis that non-random mating practices affect the distribution of these variants, through exposure in the homozygous state, leading to their purging from the population gene pool. To do so, we produced whole-genome sequencing data for two pairs of Asian populations exhibiting different alliance rules and rates of inbreeding, but with similar effective population sizes. The results show that populations with higher rates of inbred matings do not purge deleterious variants more efficiently. Purging therefore has a low efficiency in human populations, and different mating practices lead to a similar mutational load.

A fast linkage method for population GWAS cohorts with related individuals.

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

Molecular genetic diversity and linkage disequilibrium structure of the Egyptian faba bean using Single Primer Enrichment Technology (SPET).

Faba bean is an important legume crop. The genetic diversity among faba bean genotypes is very important for the genetic improvement of target traits. A set of 128 fab bean genotypes that are originally from Egypt were used in this study to investigate the genetic diversity and population structure. The 128 genotypes were genotyped using the Single Primer Enrichment Technology (SPET) by which a set of 6759 SNP markers were generated after filtration. The SNP markers were distributed on all chromosomes with a range extending from 822 (Chr. 6) to 1872 (Chr.1). The SNP markers had wide ranges of polymorphic information content (PIC), gene diversity (GD), and minor allele frequency. The analysis of population structure divided the Egyptian faba bean population into five subpopulations. Considerable genetic distance was found among all genotypes, ranging from 0.1 to 0.4. The highly divergent genotype was highlighted in this study and the genetic distance among genotypes ranged from 0.1 and 0.6. Moreover, the structure of linkage disequilibrium was studied, and the analysis revealed a low level of LD in the Egyptian faba bean population. A slow LD decay at the genomic and chromosomal levels was observed. Interestingly, the distribution of haplotype blocks was presented in each chromosome and the number of haplotype block ranged from 65 (Chr. 4) to 156 (Chr. 1). Migration and genetic drift are the main reasons for the low LD in the Egyptian faba bean population. The results of this study shed light on the possibility of the genetic improvement of faba bean crop in Egypt and conducting genetic association analyses to identify candidate genes associated with target traits (e.g. protein content, grain yield, etc.) in this panel.

Genetic diversity analysis of Inner Mongolia cashmere goats (Erlangshan subtype) based on whole genome re-sequencing.

BACKGROUND: Inner Mongolia cashmere goat (IMCG), renowned for its superior cashmere quality, is a Chinese indigenous goat breed that has been developed through natural and artificial selection over a long period. However, recently, the genetic resources of IMCGs have been significantly threatened by the introduction of cosmopolitan goat breeds and the absence of adequate breed protection systems. RESULTS: In order to assess the conservation effectiveness of IMCGs and efficiently preserve and utilize the purebred germplasm resources, this study analyzed the genetic diversity, kinship, family structure, and inbreeding of IMCGs utilizing resequencing data from 225 randomly selected individuals analyzed using the Plink (v.1.90), GCTA (v.1.94.1), and R (v.4.2.1) software. A total of 12,700,178 high-quality SNPs were selected through quality control from 34,248,064 SNP sites obtained from 225 individuals. The average minor allele frequency (MAF), polymorphic information content (PIC), and Shannon information index (SHI) were 0.253, 0.284, and 0.530, respectively. The average observed heterozygosity (Ho) and the average expected heterozygosity (He) were 0.355 and 0.351, respectively. The analysis of the identity by state distance matrix and genomic relationship matrix has shown that most individuals' genetic distance and genetic relationship are far away, and the inbreeding coefficient is low. The family structure analysis identified 10 families among the 23 rams. A total of 14,109 runs of homozygosity (ROH) were identified in the 225 individuals, with an average ROH length of 1014.547 kb. The average inbreeding coefficient, calculated from ROH, was 0.026 for the overall population and 0.027 specifically among the 23 rams, indicating a low level of inbreeding within the conserved population. CONCLUSIONS: The IMCGs exhibited moderate polymorphism and a low level of kinship with inbreeding occurring among a limited number of individuals. Simultaneously, it is necessary to prevent the loss of bloodline to guarantee the perpetuation of the IMCGs' germplasm resources.

Privacy-aware estimation of relatedness in admixed populations.

BACKGROUND: Estimation of genetic relatedness, or kinship, is used occasionally for recreational purposes and in forensic applications. While numerous methods were developed to estimate kinship, they suffer from high computational requirements and often make an untenable assumption of homogeneous population ancestry of the samples. Moreover, genetic privacy is generally overlooked in the usage of kinship estimation methods. There can be ethical concerns about finding unknown familial relationships in third-party databases. Similar ethical concerns may arise while estimating and reporting sensitive population-level statistics such as inbreeding coefficients for the concerns around marginalization and stigmatization. RESULTS: Here, we present SIGFRIED, which makes use of existing reference panels with a projection-based approach that simplifies kinship estimation in the admixed populations. We use simulated and real datasets to demonstrate the accuracy and efficiency of kinship estimation. We present a secure federated kinship estimation framework and implement a secure kinship estimator using homomorphic encryption-based primitives for computing relatedness between samples in two different sites while genotype data are kept confidential. Source code and documentation for our methods can be found at https://doi.org/10.5281/zenodo.7053352. CONCLUSIONS: Analysis of relatedness is fundamentally important for identifying relatives, in association studies, and for estimation of population-level estimates of inbreeding. As the awareness of individual and group genomic privacy is growing, privacy-preserving methods for the estimation of relatedness are needed. Presented methods alleviate the ethical and privacy concerns in the analysis of relatedness in admixed, historically isolated and underrepresented populations. SHORT ABSTRACT: Genetic relatedness is a central quantity used for finding relatives in databases, correcting biases in genome wide association studies and for estimating population-level statistics. Methods for estimating genetic relatedness have high computational requirements, and occasionally do not consider individuals from admixed ancestries. Furthermore, the ethical concerns around using genetic data and calculating relatedness are not considered. We present a projection-based approach that can efficiently and accurately estimate kinship. We implement our method using encryption-based techniques that provide provable security guarantees to protect genetic data while kinship statistics are computed among multiple sites.

Effects and considerations of multiplexing ForenSeq Kintelligence libraries with a negative control.

The negative template control or negative amplification control has been an essential component of the forensic DNA analysis workflow that helps monitor contamination. As such, the inclusion of a negative control in forensic DNA analysis has been a requirement for all laboratories audited under the FBI's Quality Assurance Standards. As massively parallel sequencing (MPS) becomes more conventional in forensic laboratories, considerations for the inclusion of a negative control in every sequencing run can be evaluated. Although the inclusion of a negative control in library preparation and the first sequencing run has a practical function, there is less utility for its inclusion in all subsequent sequencing runs for that library preparation. Although this is universal to all MPS assays, it is most relevant for an assay that has a low sample multiplexing capacity, such as the ForenSeq Kintelligence Kit (Qiagen/Verogen, Inc.). The ForenSeq Kintelligence Kit is an investigative genetic genealogy (IGG) sequencing-based assay that targets 10,230 forensically relevant single-nucleotide polymorphisms. The manufacturer recommends multiplexing 3 libraries per sequencing run, which includes controls. The purpose of this study was to investigate the effect of the inclusion of a negative control in every Kintelligence sequencing run. We observed that the library generated from a negative amplification control will take 7%-14% of the run output. The loss of sequencing space taken by a negative control decreased the available output for DNA-containing samples, leading in some cases to allele or locus dropout and accompanying higher numbers of sixth to seventh order unknown associations in GEDmatch PRO.

Unveiling facial kinship: The BioKinVis dataset for facial kinship verification and genetic association studies.

Facial image-based kinship verification represents a burgeoning frontier within the realms of computer vision and biomedicine. Recent genome-wide association studies have underscored the heritability of human facial morphology, revealing its predictability based on genetic information. These revelations form a robust foundation for advancing facial image-based kinship verification. Despite strides in computer vision, there remains a discernible gap between the biomedical and computer vision domains. Notably, the absence of family photo datasets established through biological paternity testing methods poses a significant challenge. This study addresses this gap by introducing the biological kinship visualization dataset, encompassing 5773 individuals from 2412 families with biologically confirmed kinship. Our analysis delves into the distribution and influencing factors of facial similarity among parent-child pairs, probing the potential association between forensic short tandem repeat polymorphisms and facial similarity. Additionally, we have developed a machine learning model for facial image-based kinship verification, achieving an accuracy of 0.80 in the dataset. To facilitate further exploration, we have established an online tool and database, accessible at http://120.55.161.230:88/.

From grouping and cooperation to menstruation: Spiny mice (Acomys cahirinus) are an emerging mammalian model for sociality and beyond.

While spiny mice are primarily used as a model for Type II diabetes and for studying complex tissue regeneration, they are also an emerging model for a variety of studies examining hormones, behavior, and the brain. We began studying the spiny mouse to take advantage of their highly gregarious phenotype to examine how the brain facilitates large group-living. However, this unique rodent can be readily bred and maintained in the lab and can be used to ask a wide variety of scientific questions. In this brief communication we provide an overview of studies that have used spiny mice for exploring physiology and behavior. Additionally, we describe how the spiny mouse can serve as a useful model for researchers interested in studying precocial development, menstruation, cooperation, and various grouping behaviors. With increasingly available technological advancements for non-traditional organisms, spiny mice are well-positioned to become a valuable organism in the behavioral neuroscience community.

Fine-scale genetic structure in the orchid Gymnadenia conopsea is not associated with local density of flowering plants.

PREMISE: Density-dependent pollinator visitation can lead to density-dependent mating patterns and within-population genetic structure. In Gymnadenia conopsea, individuals in low-density patches receive more self pollen than individuals in high-density patches, suggesting higher relatedness at low density. Ongoing fragmentation is also expected to cause more local matings, potentially leading to biparental inbreeding depression. METHODS: To evaluate whether relatedness decreases with local density, we analyzed 1315 SNP loci in 113 individuals within two large populations. We quantified within-population genetic structure in one of the populations, recorded potential habitat barriers, and visualized gene flow using estimated effective migration surfaces (EEMS). We further estimated the magnitude of biparental inbreeding depression that would result from matings restricted to within 5 m. RESULTS: There was no significant relationship between local density and relatedness in any population. We detected significant fine-scale genetic structure consistent with isolation by distance, with positive kinship coefficients at distances below 10 m. Kinship coefficients were low, and predicted biparental inbreeding depression resulting from matings within the closest 5 m was a modest 1-3%. The EEMS suggested that rocks and bushes may act as barriers to gene flow within a population. CONCLUSIONS: The results suggest that increased self-pollen deposition in sparse patches does not necessarily cause higher selfing rates or that inbreeding depression results in low establishment success of inbred individuals. The modest relatedness suggests that biparental inbreeding depression is unlikely to be an immediate problem following fragmentation of large populations. The results further indicate that habitat structure may contribute to governing fine-scale genetic structure in G. conopsea.

The influence of biological relatedness on sexual dimorphism and sex classification based on external morphology of the frontal bone.

The most significant sexual differences in the human skull are located in the upper third of the face (the frontal bone), which is a useful research object, mainly in combination with virtual anthropology methods. However, the influence of biological relatedness on sexual dimorphism and frontal bone variability remains unknown. This study was directed at sexual difference description and sex classification using the form and shape of the external surface of the frontal bones from a genealogically documented Central European osteological sample (nineteenth to twentieth centuries). The study sample consisted of 47 cranial CT images of the adult members of several branches of one family group over 4 generations. Three-dimensional virtual models of the frontal bones were analyzed using geometric morphometrics and multidimensional statistics. Almost the entire external frontal surface was significantly different between males and females, especially in form. Significant differences were also found between this related sample and an unrelated one. Sex estimation of the biologically related individuals was performed using the classification models developed on a sample of unrelated individuals from the recent Czech population (Cechová et al. in Int J Legal Med 133: 1285 1294, 2019), with a result of 74.46% and 63.83% in form and shape, respectively. Failure of this classifier was caused by the existence of typical traits found in the biologically related sample different from the usual manifestation of sexual dimorphism. This can be explained as due to the increased degree of similarity and the reduction of variability in biologically related individuals. The results show the importance of testing previously published methods on genealogical data.

Genotyping SNPs in lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2) exhibits association with wood density in teak (Tectona grandis L.f.).

BACKGROUND: Teak (Tectona grandis L.f.), an important source of tropical timber with immense economic value, is a highly outcrossing forest tree species. 150 unrelated accessions of teak (Tectona grandis L.f.) plus trees assembled as clones at National Teak Germplasm Bank, Chandrapur, Maharashtra, India was investigated for association mapping of candidate lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2). METHODS AND RESULTS: The CAD1, MYB1 and MYB2 were amplified using specifically designed primers. The amplified sequences were then sequenced and genotyped for 112 SNPs/11 indels. We evaluated the association between SNPs and wood density in teak accessions using GLM and MLM statistical models, with Bonferroni correction applied. The teak accessions recorded an average wood density of 416.69 kg.m-3 (CV 4.97%) and comprised of three loosely structured admixed sub-populations (K = 3), containing 72.05% genetic variation within sub-populations with low intragenic LD (0-21% SNP pairs) at P < 0.05 and high LD decay (33-934 bp) at R2 = 0.1. GLM and MLM models discounting systematic biases (Q and K matrices) to avoid false discovery revealed five loci at rare variants (MAF 0.003) and three loci at common variants (MAF 0.05) to be significantly (P < 0.05) associated with the wood density. However, the stringent Bonferroni correction (4.06-7.04 × 10-4) yielded only a single associated locus (B1485C/A) from exon of MYB1 transcription factor, contributing to about 10.35% phenotypic variation in wood density trait. CONCLUSION: Scored SNP locus (B1485C/A) can be developed as a molecular probe for selection of improved planting stock with proven wood density trait for a large-scale teak plantation.

Implications of methodologies for integrating empirical kinships into ex situ population management using PMx: A case study of Baer's Pochard (Aythya baeri) in North America.

The application of molecular tools to population management can improve the long-term genetic viability of ex situ populations. In this study, we aimed to understand the implications of integrating empirical kinships into the genetic management of an ex situ population of the endangered waterfowl, Baer's pochard (Aythya baeri), in North America. Single nucleotide polymorphism data were generated for 141 Baer's pochard using double digest restriction site-associated DNA sequencing and empirical kinships were derived and integrated into the population management software PMx. Analyses suggested 37.7% of pairwise relationships previously assumed to be unrelated were first, second, or third-order relatives. We determined that most genetic summary statistics were impacted through the calculation of the population's mean kinship, which increased from MK¯=0.0772 to MK¯=0.2074 after empirical kinships were integrated into our analyses. Our results also revealed the importance of understanding how molecular kinships derived from a particular estimator are scaled, if the scale differs significantly from pedigree-based kinships. We describe the theory behind the genetic metrics impacted and provide general guidance on incorporating empirical kinships into ex situ population management as well as provide suggestions for sampling strategies to minimize the biases inherent in merging two types of kinship estimators.

United by Conflict: Convergent Signatures of Parental Conflict in Angiosperms and Placental Mammals.

Endosperm in angiosperms and placenta in eutherians are convergent innovations for efficient embryonic nutrient transfer. Despite advantages, this reproductive strategy incurs metabolic costs that maternal parents disproportionately shoulder, leading to potential inter-parental conflict over optimal offspring investment. Genomic imprinting-parent-of-origin-biased gene expression-is fundamental for endosperm and placenta development and has convergently evolved in angiosperms and mammals, in part, to resolve parental conflict. Here, we review the mechanisms of genomic imprinting in these taxa. Despite differences in the timing and spatial extent of imprinting, these taxa exhibit remarkable convergence in the molecular machinery and genes governing imprinting. We then assess the role of parental conflict in shaping evolution within angiosperms and eutherians using four criteria: (1) Do differences in the extent of sibling relatedness cause differences in the inferred strength of parental conflict? (2) Do reciprocal crosses between taxa with different inferred histories of parental conflict exhibit parent-of-origin growth effects? (3) Are these parent-of-origin growth effects caused by dosage-sensitive mechanisms and do these loci exhibit signals of positive selection? (4) Can normal development be restored by genomic perturbations that restore stoichiometric balance in the endosperm/placenta? Although we find evidence for all criteria in angiosperms and eutherians, suggesting that parental conflict may help shape their evolution, many questions remain. Additionally, myriad differences between the two taxa suggest that their respective biologies may shape how/when/where/to what extent parental conflict manifests. Lastly, we discuss outstanding questions, highlighting the power of comparative work in quantifying the role of parental conflict in evolution.

The Roots of Science, Technology, Engineering, and Mathematics: What Are the Evolutionary and Neural Bases of Human Mathematics and Technology?

INTRODUCTION: Neural exaptations represent descent via transitions to novel neural functions. A primary transition in human cognitive and neural evolution was from a predominantly socially oriented primate brain to a brain that also instantiates and subserves science, technology, and engineering, all of which depend on mathematics. Upon what neural substrates and upon what evolved cognitive mechanisms did human capacities for science, technology, engineering, and mathematics (STEM), and especially its mathematical underpinnings, emerge? Previous theory focuses on roles for tools, language, and arithmetic in the cognitive origins of STEM, but none of these factors appears sufficient to support the transition. METHODS: In this article, I describe and evaluate a novel hypothesis for the neural origins and substrates of STEM-based cognition: that they are based in human kinship systems and human maximizing of inclusive fitness. RESULTS: The main evidence for this hypothesis is threefold. First, as demonstrated by anthropologists, human kinship systems exhibit complex mathematical and geometrical structures that function under sets of explicit rules, and such systems and rules pervade and organize all human cultures. Second, human kinship underlies the core algebraic mechanism of evolution, maximization of inclusive fitness, quantified as personal reproduction plus the sum of all effects on reproduction of others, each multiplied by their coefficient of relatedness to self. This is the only "natural" equation expected to be represented in the human brain. Third, functional imaging studies show that kinship-related cognition activates frontal-parietal regions that are also activated in STEM-related tasks. In turn, the decision-making that integrates kinship levels with costs and benefits from alternative behaviors has recently been shown to recruit the lateral septum, a hub region that combines internal (from the prefrontal cortex, amygdala, and other regions) and external information relevant to social behavior, using a dedicated subsystem of neurons specific to kinship. CONCLUSIONS: Taken together, these lines of evidence suggest that kinship systems and kin-associated behaviors may represent exaptations for the origin of human STEM.