Hostname: page-component-78c5997874-t5tsf Total loading time: 0 Render date: 2024-11-04T21:38:53.125Z Has data issue: false hasContentIssue false

Cross-cousin marriage among Tsimane forager–horticulturalists during demographic transition and market integration

Published online by Cambridge University Press:  26 March 2024

Arianna Dalzero*
Affiliation:
Department of Human Behavior, Ecology and Culture, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
Bret A. Beheim
Affiliation:
Department of Human Behavior, Ecology and Culture, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
Hillard Kaplan
Affiliation:
Economic Science Institute, Chapman University, Orange, California, USA
Jonathan Stieglitz
Affiliation:
Toulouse School of Economics and Institute for Advanced Study, University of Toulouse 1 Capitole, Toulouse, France
Paul L. Hooper
Affiliation:
Economic Science Institute, Chapman University, Orange, California, USA Department of Anthropology, University of New Mexico, Albuquerque, New Mexico, USA
Cody T. Ross
Affiliation:
Department of Human Behavior, Ecology and Culture, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
Michael Gurven
Affiliation:
Department of Anthropology, University of California, Santa Barbara, California, USA
Dieter Lukas
Affiliation:
Department of Human Behavior, Ecology and Culture, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
*
Corresponding author: Arianna Dalzero; E-mail: [email protected]

Abstract

Although still prevalent in many human societies, the practice of cousin marriage has precipitously declined in populations undergoing rapid demographic and socioeconomic change. However, it is still unclear whether changes in the structure of the marriage pool or changes in the fitness-relevant consequences of cousin marriage more strongly influence the frequency of cousin marriage. Here, we use genealogical data collected by the Tsimane Health and Life History Project to show that there is a small but measurable decline in the frequency of first cross-cousin marriage since the mid-twentieth century. Such changes are linked to concomitant changes in the pool of potential spouses in recent decades. We find only very modest differences in fitness-relevant demographic measures between first cousin and non-cousin marriages. These differences have been diminishing as the Tsimane have become more market integrated. The factors that influence preferences for cousin marriage appear to be less prevalent now than in the past, but cultural inertia might slow the pace of change in marriage norms. Overall, our findings suggest that cultural changes in marriage practices reflect underlying societal changes that shape the pool of potential spouses.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press

Social media summary: Traditions in a changing world: partner availability, fitness effects and the decline in cousin marriage among Tsimane.

1. Introduction

Many human populations are experiencing rapid demographic, socioeconomic and cultural changes, including demographic transition (Borgerhoff Mulder, Reference Borgerhoff Mulder1998), industrialisation, market integration and urbanisation (Mattison & Sear, Reference Mattison and Sear2016). These changes affect – and are affected by – family composition and norms governing spousal partner choice. Cousin marriage – a spousal union between individuals who share one or more grandparents – used to be practised at high rates in many human societies (Shaw & Raz, Reference Shaw and Raz2015). It is still estimated to account for more than 10% of spousal unions worldwide (Bittles & Black, Reference Bittles and Black2010), but its frequency has been declining in many societies undergoing demographic and socioeconomic transitions (Sahoo et al., Reference Sahoo, Debnath, Mandal, Nagarajan and Appunni2022; Brennan, Reference Brennan1981; Kalam et al., Reference Kalam, Sharma, Ghosh and Roy2021). This has led some to forecast a global decline in cousin marriages (Bittles & Black, Reference Bittles and Black2010). Changes in the frequency of cousin marriage are often associated with changes in other cultural norms – e.g. owing to particular religious proscriptions (Cazes, Reference Cazes1996; Schulz et al., Reference Schulz, Bahrami-Rad, Beauchamp and Henrich2019) – but can also reflect broader changes in demographic and socioeconomic circumstances that might influence marriage decisions indirectly (Shenk et al., Reference Shenk, Towner, Voss and Alam2016). More specifically, changes in demographic and socioeconomic circumstances have often been shown to co-occur with changes in the frequency of cousin marriage (Peña et al., Reference Peña, Alfonso-Sánchez and Calderón2002; Givens & Hirschman, Reference Givens and Hirschman1994; Assaf & Khawaja, Reference Assaf and Khawaja2009; Kumari et al., Reference Kumari, Bittles and Saxena2020; Shenk et al., Reference Shenk, Towner, Voss and Alam2016). Demographic circumstances, for example, are expected to shape the pool of potential marriage partners (Hajnal, Reference Hajnal1963). Socioeconomic circumstances might also influence the potential fitness returns to particular marriage decisions (Ross et al., Reference Ross, Borgerhoff Mulder, Oh, Bowles, Beheim, Bunce and Ziker2018), and thus affect people's preferences for them (Borgerhoff Mulder & Ross, Reference Borgerhoff Mulder and Ross2019; Micheletti, Reference Micheletti2020). In this paper, we investigate whether – and how – changes in the occurrence of first cousin marriage are associated with changes in the demographic composition of the pool of eligible partners and the potential fitness returns to marrying a first cousin. We focus our case study on the Tsimane forager–horticulturalists of the Bolivian Amazon. While we consider socioeconomic factors like wealth, social network structure and residence patterns in interpreting these dynamics, our analysis does not directly delve into modelling such factors, owing to a lack of data (see Study Limitations).

A central question raised by researchers studying cousin marriage is whether individuals choose their spouses based on different characteristics and/or simply search for potential partners from a broader pool of individuals during periods of demographic and socioeconomic changes (e.g. owing to globalising forces) (Goode, Reference Goode1963; Givens & Hirschman, Reference Givens and Hirschman1994). The prevalence of cousin marriage can be high in societies where the number of suitable partners is limited – e.g. in geographically isolated populations – and studies have shown that the availability of cousins as partners influences partner choice (Mobarak et al., Reference Mobarak, Chaudhry, Brown, Zelenska, Khan, Chaudry and Li2019; Barakat & Basten, Reference Barakat and Basten2014). In such circumstances, individuals may choose cousins as partners in order to avoid the costs of migration to find a mate and to increase their marriage opportunities (see below for details on the potential fitness benefits of cousin marriage) (Hoben et al., Reference Hoben, Buunk and Fisher2016; Cavalli-Sforza et al., Reference Cavalli-Sforza, Moroni and Zei2013). An increase in social connections between isolated populations and larger population centres is often associated with a decrease in the frequency of cousin marriage (Pettener, Reference Pettener1985). Past studies have thus hypothesised a causal connection between changes in the frequency of cousin marriage over time and shifts in the composition of the pool of potential mates (Givens & Hirschman, Reference Givens and Hirschman1994; Helgason et al., Reference Helgason, Palsson, Guðbjartsson, Kristjansson and Stefansson2008). In particular, past work has linked early periods of demographic transition to an increase in fertility and a decrease in childhood mortality rates – circumstances that may lead to greater availability of kin among potential partners (Givens & Hirschman, Reference Givens and Hirschman1994; Peña et al., Reference Peña, Alfonso-Sánchez and Calderón2002; Pettener, Reference Pettener1985; Murphy, Reference Murphy2011; Brennan, Reference Brennan1981). These studies, however, are generally focused on closed populations or on short-term outcomes, and it is unclear how the pool of eligible partners might change in the long term, and in more interconnected, growing populations, where changes in the number of unrelated potential partners may exceed changes in the number of kin (Lukas et al., Reference Lukas, Reynolds, Boesch and Vigilant2005).

Socioeconomic shifts – such as market integration – might have differential impacts on the fitness outcomes of cousin vs. non-cousin marriages – e.g. through differences in moderation of the direct biological effects of marriage practice or through impacts on the socioeconomic factors that influence marriage partner choice (Mobarak et al., Reference Mobarak, Chaudhry, Brown, Zelenska, Khan, Chaudry and Li2019; Dalzero et al., Reference Dalzero, Ross and Lukas2023). Changes in fitness returns to cousin marriage might in turn cause shifts in individuals’ preferences to marry a cousin. Fitness outcomes and the frequency of cousin marriage are closely tied to the particular socioeconomic conditions of a population. Cousin marriage rates are notably high in specific societies – especially agricultural and horticultural ones – where it offers advantages like increased access to kin-based wealth (Johow et al., Reference Johow, Willführ and Voland2019) and alliances (Lévi-Strauss, Reference Lévi-Strauss1971), greater opportunities for marriage, and a smoother and faster marriage process through matrimonial exchanges and agreements (Hoben et al., Reference Hoben, Buunk and Fisher2016; Cavalli-Sforza et al., Reference Cavalli-Sforza, Moroni and Zei2013; Chagnon et al., Reference Chagnon, Lynch, Shenk, Hames and Flinn2017). All of these factors can contribute to fitness gains – e.g. through increased fertility and offspring survival, and earlier age at reproduction (Bailey et al., Reference Bailey, Hill and Walker2014; Chagnon et al., Reference Chagnon, Lynch, Shenk, Hames and Flinn2017; Mobarak et al., Reference Mobarak, Chaudhry, Brown, Zelenska, Khan, Chaudry and Li2019; Dalzero et al., Reference Dalzero, Ross and Lukas2023). Fitness gains may offset the potential cost arising from the biological consequences of mating with close relatives (Clark et al., Reference Clark, Okada, Moore, Mason, Pirastu, Gandin and Wilson2019). Preferences for cousin marriage have indeed been shown to increase in association with dowry and land ownership systems that prioritise family assets (Peña et al., Reference Peña, Alfonso-Sánchez and Calderón2002). In contrast, declines in cousin marriage have also been linked to economic changes that facilitate more resource independence for women – e.g. increased levels of education, availability of wage labour and engagement in a market economy appear to attenuate the frequency of cousin marriage, rendering intensive kinship systems and kin-based alliances less beneficial (Assaf & Khawaja, Reference Assaf and Khawaja2009; Bras et al., Reference Bras, Van Poppel and Mandemakers2009; Islam, Reference Islam2018; Brennan, Reference Brennan1981; Koc, Reference Koc2008; Shenk et al., Reference Shenk, Towner, Voss and Alam2016; Schaffnit et al., Reference Schaffnit, Page, Lynch, Spake, Sear, Sosis and Shenk2022).

Importantly, people practising cousin marriage might also be in different socioeconomic conditions than people not practising it – e.g. in terms of residence, social and kin networks, and wealth – and market integration might affect these conditions (Ghosh et al., Reference Ghosh, Hwang and Squires2023). Socioeconomic changes are also expected to alter the trade-offs between different components of fitness in cousin and non-cousin marriages. For example, if a high fertility rate and early age at reproduction trade off against child or maternal survival, this relationship is likely to be influenced by shifts in healthcare or reproductive behaviours (Davison & Gurven, Reference Davison and Gurven2021). Recent literature on cross-cultural patterns in cousin marriage practices documents the correlates (e.g. industrialisation, dowry, and cultural traditions) of changes in the frequency of cousin marriage, especially as they relate to economic development (Shenk et al., Reference Shenk, Towner, Voss and Alam2016; Bittles & Black, Reference Bittles and Black2010; Caldwell et al., Reference Caldwell, Reddy and Caldwell1983). However, most studies in the literature were not able to determine how changes in partner choice relate to realised differences in fitness returns, as such analyses require data on the costs and benefits of different marriage strategies in terms of both parent reproduction and offspring survival. Likewise, studies that have examined the relationship between kinship and fertility – or offspring survival – in human populations (Mobarak et al., Reference Mobarak, Chaudhry, Brown, Zelenska, Khan, Chaudry and Li2019; Helgason et al., Reference Helgason, Palsson, Guðbjartsson, Kristjansson and Stefansson2008; Husain & Bunyan, Reference Husain and Bunyan1997; Hussain & Bittles, Reference Hussain and Bittles2004; Bittles et al., Reference Bittles, Grant, Sullivan and Hussain2002) did not test for temporal dynamics in these fitness components for cousin and non-cousin marriages. We hope to fill in these gaps in understanding using data from a large, longitudinal sample of individuals from the Tsimane forager–horticulturalists of the Bolivian Amazon.

The Tsimane are a high-fertility, isolated population with limited use of contraception. Tsimane cultural practices feature cross-cousin marriage – where a focal man marries the daughter of his mother's brother or his father's sister. Traditionally, Tsimane stigmatise those who do not marry a cross-cousin, who – according to local legend – are believed to become jaguars after they die (Godoy et al., Reference Godoy, Eisenberg, Reyes-García, Huanca, Leonard, McDade and Tanner2008). While marrying a first cross-cousin (fom cheyas) is considered ideal, in cases when first cousins are not available, the second option is to marry a more distant cousin in the cross lineage (fom moch) (Ellis, Reference Ellis1997). In demographic interviews conducted between 2002 and 2005 in 18 villages across the Tsimane territory (M. Gurven, unpublished), 64 couples (31%) reported being first cross-cousins, 92 couples (45%) reported being distant cross-cousins and 48 couples (24%) reported being unrelated. Similarly, other estimates suggest that the frequency of cross-cousin marriage (including first and distant cousins) is 75% in some areas (unpublished pilot survey reported in Patel et al., Reference Patel, Godoy, Seyfried, Reyes-García, Huanca, Leonard and Tanner2007). Opposite-sex siblings often arrange marriages between their respective children (Winking et al., Reference Winking, Gurven and Kaplan2011), or otherwise exert informal pressure on their children's choice of partners (Godoy et al., Reference Godoy, Eisenberg, Reyes-García, Huanca, Leonard, McDade and Tanner2008). However, marriages are often decided by the couple, in accordance with the custom of preferential cross-cousin marriage. Tsimane marriage partners are generally similar to each other in age and wealth. Tsimane do not commemorate marriages with formal ceremonies, but rather consider a pair to be married when they sleep together in the same house. Newlyweds reside with the wife's family until the birth of their first child, at which point they typically transition to a neolocal residence (Godoy et al., Reference Godoy, Eisenberg, Reyes-García, Huanca, Leonard, McDade and Tanner2008). The median age of marriage among Tsimane is 21 for males and 16.5 for females. Women, on average, become mothers by the age of 18 (Martin & Gurven, Reference Martin and Gurven2022). The Tsimane are primarily endogamous, occasionally intermarrying with lowland groups who have settled in Tsimane territory (Pisor & Gurven, Reference Pisor and Gurven2018). Monogamy is the norm among Tsimane (95% of married men, with 5% practising sororal polygyny), and divorce is rare.

The Tsimane are undergoing a socioeconomic shift from a horticultural subsistence economy to a more market-integrated economy, and simultaneously experiencing a demographic transition (Gurven et al., Reference Gurven, Stieglitz, Trumble, Blackwell, Beheim, Davis and Kaplan2017). Both of these factors may influence the availability of potential partners, the choices individuals make in selecting partners, and individuals’ survival and reproduction. Most Tsimane still subsist on slash-and-burn horticulture, fishing and hunting. However, owing to an agrarian reform in 2006 that distributed forests to highland colonist farmers, they are becoming increasingly dependent on cash cropping, selective logging and/or on employment in logging camps, cattle ranches and farms (Reyes-García et al., Reference Reyes-García, Paneque-Gálvez, Luz, Gueze, Macía, Orta-Martínez and Pino2014). The Tsimane are in the first phase of a demographic transition, with an annual population growth rate of over 3.5% (Gurven et al., Reference Gurven, Stieglitz, Trumble, Blackwell, Beheim, Davis and Kaplan2017). This population growth results from a combination of reduced child and adult mortality, owing to increased access to healthcare facilities, and individuals maintaining a preference for high fertility, with women conceiving an average of nine children over their life course (Kaplan et al., Reference Kaplan, Hooper, Stieglitz and Gurven2015; Mcallister et al., Reference Mcallister, Gurven, Kaplan and Stieglitz2012). While high fertility and the fall of death rates may lead to increases in the number of cross-cousins available as potential partners (Givens & Hirschman, Reference Givens and Hirschman1994), this increase in number of kin may be outpaced by increases in the number of unrelated individuals available as potential partners, owing to the overall growth of the population (Lukas et al., Reference Lukas, Reynolds, Boesch and Vigilant2005), especially in more connected areas (Pettener, Reference Pettener1985; Brennan, Reference Brennan1981).

In fact, while Tsimane have traditionally lived in small dispersed settlements, often connected by rivers, and mainly composed of extended family clusters, they now tend to aggregate in larger settlements, closer to market towns, roads and logging camps (Gurven et al., Reference Gurven, Stieglitz, Trumble, Blackwell, Beheim, Davis and Kaplan2017). Changes in residence and movement patterns may lead to an increase in contact with unrelated potential partners. Traditionally, beginning in adolescence, Tsimane men often travel to visit relatives in other villages. Many young men use these frequent visitations to relatives as an opportunity to seek their future spouses among their cross-cousins and to obtain permission for the marriage from future in-laws (Miner et al., Reference Miner, Gurven, Kaplan and Gaulin2014). According to Godoy et al. (Reference Godoy, Reyes-García, Seyfried, Huanca, Leonard and McDade2007), 79% of individuals find their spouse in a different village. In contrast, members of more connected and larger communities travel more often to market towns (Miner et al., Reference Miner, Gurven, Kaplan and Gaulin2014), possibly expanding their social networks and increasing their chances of establishing relationships with non-kin. As a consequence of migration and re-settlement, changes in kin-based and broader sharing networks might also affect the potential fitness benefits of marrying a cross-cousin (Hooper et al., Reference Hooper, Gurven, Winking and Kaplan2015; Sear & Coall, Reference Sear and Coall2011). Tsimane sharing networks are strongly shaped by kinship (Hooper et al., Reference Hooper, Gurven, Winking and Kaplan2015). However, the access to family wealth and alliances that might have traditionally increased the reproductive success of people practising cousin marriage might be becoming less important (Borgerhoff Mulder et al., Reference Borgerhoff Mulder, Bowles, Hertz, Bell, Beise, Clark and Wiessner2009; Johow et al., Reference Johow, Willführ and Voland2019). With increasing reliance on a cash and wage-based economy, factors other than kin ties might be more relevant for partner choice – for example, individuals who speak Spanish might be becoming more attractive partners because of their ability to communicate with individuals in the wider market economy, thus providing better access to resources important for family well-being and reproductive success (Godoy et al., Reference Godoy, Reyes-García, Seyfried, Huanca, Leonard and McDade2007). Access to modern medicine, and the resulting reduction in mortality (Gurven & Kaplan, Reference Gurven and Kaplan2007), might also reduce differences in fitness between families practising different marriage strategies. In these novel circumstances, individuals might be expected to choose their partners based on new, more personal criteria, rather than adhering to cultural norms or family preferences. Such changes in choice might be further facilitated by a breakdown of traditional and parental authority (Reyes-García et al., Reference Reyes-García, Paneque-Gálvez, Luz, Gueze, Macía, Orta-Martínez and Pino2014; Godoy et al., Reference Godoy, Eisenberg, Reyes-García, Huanca, Leonard, McDade and Tanner2008).

Our focus here is to understand whether normative first cross-cousin marriage might have become less attractive over time in this population, potentially because of changes in the composition of the pool of eligible partners or because of reduced fitness benefits of kin-intensive networks. The detailed long-term data collected by the Tsimane Health and Life History Project (THLHP; Gurven et al., Reference Gurven, Stieglitz, Trumble, Blackwell, Beheim, Davis and Kaplan2017) afford a unique opportunity to unravel the dynamics of cousin marriage in a changing society. We rely solely on the THLHP genealogical data, because these data span the period of time in which the Tsimane livelihood began to transition from horticultural subsistence to market integration, and include individuals born from the end of the nineteenth century to the beginning of the twenty-first century. We begin by quantitatively assessing the extent to which the frequency of first cross-cousin marriage has changed over this several decade period. We then investigate the role of partner availability on the occurrence of first cross-cousin marriage. Specifically, we investigate whether the composition of the pool of potential partners has changed over time, and whether such changes can explain changes in the frequency of first cross-cousin marriage. Next, we explore the age-specific fitness costs and benefits of first cross-cousin marriage for both women and men, by comparing life history traits (i.e. offspring child survival rates, fertility rates, and ages at first reproduction) that might be expected to differ between families practising cousin marriage and families not practising it (Dalzero et al., Reference Dalzero, Ross and Lukas2023; Bailey et al., Reference Bailey, Hill and Walker2014). The joint estimation of these life history traits in demographic models can potentially identify trade-offs that might occur. Finally, we examine whether any patterns in these life-history traits have changed over time. Table 1 provides an overview of our questions, and inferential estimands.

Table 1. Research questions and inferential estimands

2. Methods

2.1. Data source

The THLHP has been conducting genealogical and reproductive history interviews in the Tsimane territory since 2002. The study encompasses 90 villages that are almost exclusively occupied by Tsimane people. The Tsimane population register contains data at the individual level, spanning a period from 1870 to 2014, and currently describes 28,093 individuals. The register includes information on the unique ID, gender, birth date and death date for each individual. For 30% of the register, there is additional data on the biological and socially recognised parents of each focal individual.

2.2. Data preparation and relatedness classification

We use genealogical information to classify marriage types. We test if individuals are in cousin marriages by assessing if individuals who are partners have at least one grandparent in common. In order to distinguish (first degree) cross-cousins (hereafter referred to as ‘cousins’) from parallel-cousins, we checked for any shared maternal and paternal grandparents in each couple. We excluded a small number of parallel cousins (four marriages) from our empirical analyses, and we excluded parallel cousins from the set of eligible partners, because parallel cousin marriage is considered incestuous among Tsimane. Since information on great-grandparents is unavailable for the majority of couples, we classify as cross-cousins only those who are first-degree cross-cousins – i.e. the individuals referred to as ‘fom cheyas’ in the Tsimane language. Distant cousins (‘fom moch’) are categorised as non-cousins (see Study Limitations). As such, our inferences here are specifically limited to patterns in first degree cross-cousin marriage.

In order to identify couples, we assume that two individuals are married if they share parentage for a child, as we lack specific information on marriages per se (see the section ‘Study limitations’). We test if marriages are between first cousins using two different relatedness classification techniques. Our ‘relaxed’ condition classifies all marriages in which at least one of the eight grandparents are known for each partner. Using this classification technique, we find 546 individuals married to cousins, 4217 individuals married to non-cousins, 1230 individuals born to cousins and 12081 individuals born to non-cousins. Because some grandparents are unknown, however, this approach necessarily under-estimates the proportion of cousin marriages (see the section ‘Study limitations’). Our ‘strict’ relatedness classification, in contrast, only includes marriages in which all eight grand-parents are known. Using this technique, we find 319 individuals married to cousins, 808 individuals married to non-cousins, 607 individuals born to cousins and 1920 individuals born to non-cousins. Thus, between 11 and 28% of Tsimane marriages – for which assessment is possible – are between first degree cross-cousins. We use the ‘strict’ relatedness classification for the analyses on the frequency of cousin marriage and eligible partners (questions 1–3), to avoid classifying people as not having cousins available as potential partners just because their grandparents are unknown. We use the ‘relaxed’ relatedness classification for the analyses assessing age-specific fitness and change in fitness over time (questions 4 and 5), in order to increase sample size, especially among individuals born longer ago. The absence of information about grandparents is less problematic in the fitness analyses because we do not need to identify potential cousin partners here (see the section ‘Study limitations’). Repeating the fitness analyses with the sample of individuals classified using the ‘strict’ relatedness criterion leads to qualitatively similar results.

In order to estimate the demographic composition of the pool of potential partners for each focal person, we first calculate the age difference between partners who are in unions – i.e. people who are classified as mothers or fathers of the same children – and summarise this range using the 90% HPDI (highest posterior density interval). Next, we calculate the number of potential eligible cousin and non-cousin partners (excluding the parallel cousins) for each individual at the age when they have their first child. Using the age at the first child as a proxy for when an individual forms a partnership, we rely on the observation that there is a relatively short time gap between the formation of marital union and the birth of the first child among Tsimane. The ‘potential partners’ of each individual are defined by assessing the number of unmarried opposite sex individuals who are within an appropriate age range, based on the empirical distribution. Male potential partners can be up to 9 years older or 1 year younger than a focal woman, and female potential partners can be up to 9 years younger or 1 year older than a focal man. We include individuals who are already married (either to cousins or not) in the analyses of the frequency of cousin marriage and eligible partners (questions 1–3) and in the analyses of fertility and age at first reproduction (questions 4.2, 4.3, 5.2 and 5.3). To estimate fertility for each marriage type, we count the number of offspring of each individual and note their age at each of their births (AOB). Similarly, to estimate offspring survival for each marriage type, we note the age of death (AOD) or age of censoring (AOC) for each child – the AOC represents the age of individuals when they were last sampled, if they are still living.

2.3. Data analyses

For the first question (change in the frequency of cousin marriage over time), we estimate the association between marriage type and annual birth cohort. For the second group of questions (frequency of cousin marriage in relation to the set of eligible partners), we estimate the association between marriage type and the number of cousins and non-cousins available as partners, as well as the association between marriage type and the frequency of cousins in the pool of potential partners (i.e. the number of cousins divided by the number of cousins and non-cousins available as potential partners). For the third group of questions (change in the pool of eligible partners over time), we estimate the association between the frequency of cousins in the pool of potential partners and the annual birth cohort of individuals. Next, we estimate the association between marriage type and the frequency of cousins in the pool of potential partners, controlling for secular trends owing to annual birth cohort. We use annual birth cohort (i.e. birth year) to measure change in the frequency of cousin marriage and change in the pool of potential eligible partners over time (questions 1–3). For the analyses including annual birth cohort, we restricted inference to the years in which at least 10 individuals with known marriage relation were recorded – i.e. the period between 1964 and 1993.

For the fourth group of questions (age-specific fitness), we estimate age-specific probabilities of offspring survival by marriage type. We also estimate age-specific probabilities of reproduction by marriage type. For these analyses, we include individuals born from 1870 to 2014, the full time-span of the database. For the fifth group of questions (change in fitness over time), we estimate the probability of survival from birth to age 5 as a function of decadal birth cohort and parental marriage type. Additionally, we estimate fertility (i.e. number of children) and age at first reproduction, as a function of decadal birth cohort and marriage type. We use decadal birth cohort (i.e. birth decade) to measure change in fitness over time (question 5). We summarise the data by decade in order to capture potential non-linear shifts in key fitness parameters over time. We consider decades of birth in which at least 10 births were recorded for each sub-population (i.e. each group of individuals married or born to cousins or to non-cousins). In the offspring survival model, we include data from 1950 to 2000, and in the fertility and age at first reproduction models, we include data from 1940 to 2000 (see the section ‘Study limitations’).

We perform analysis using R, version 4.3.1 (R Core Team, Reference Core Team2021). We use the rethinking package (version 2.21) to construct all statistical models (McElreath, Reference McElreath2020), apart from the age-specific survival and fertility models, for which we used CmdStan, version 2.33.1, directly (Stan Development Team, Reference Team2021; Carpenter et al., Reference Carpenter, Gelman, Hoffman, Lee, Goodrich, Betancourt and Riddell2017). We use the rethinking package in R (McElreath, Reference McElreath2020) to summarise the results. All the code needed to reproduce our analyses will be maintained at: https://doi.org/10.17605/OSF.IO/PK8F7. Individual-level data are stored in the THLHP Data Repository, and are available through protected access protocols (see the ‘Research transparency, reproducibility and data availability’ section). Data were drawn from the THLHP database as of 2 November 2021. For more information, please see: https://tsimane.anth.ucsb.edu/data.html.

2.4. Analyses of the frequency of cousin marriage over time

In order to see if the frequency of cousin marriage has changed over time (question 1), we estimate the association between marriage type (i.e. cousin vs. non-cousin marriage) and the annual birth cohort of individuals. We build a logistic regression that is specified by the following set of equations:

$$ \hskip 16ptC_{[ i ] }\sim {\rm Bernoulli}( {\theta_{[ i] }} ) $$
$$ \eqalign{ \matrix{ { {\rm logit}( {\theta_{[ i] }} ) = \beta _1 } + {\beta _2{\rm B}_{[ i ] } } \cr \beta \sim {{\rm Normal}( {0, \;2.5} ) } \cr }} $$

where the response variable, C, is a binary indicator of whether an individual ever participated in a cousin marriage, and the predictor variable, B, is the annual birth year of each individual (i.e. of individual i). The parameter vector β gives the estimated coefficients – i.e. the intercept and the slope on annual birth cohort. The slope coefficient on annual birth cohort (β 2) reflects how the frequency of cousin marriage has changed as a logit–linear function of time.

2.5. Analyses on frequency of cousin marriage in relation to eligible partners

In order to understand whether cousin marriage among the Tsimane is limited by the number of cousins that are available as potential partners (question 2.1), we estimate the association between marriage type and the number of cousins available as potential partners for each individual. We build a binomial logistic regression specified by the following set of equations:

$$ \eqalign{ \matrix{ {C_{[ i ] }}\sim {{\rm Bernoulli}( {\theta_{[ i] }} ) } \cr {{\rm logit}( {\theta_{[ i] }} ) = \beta _1 } + {\beta _2{\rm Q}_{[ i ] } } \cr \beta \sim {{\rm Normal}( {0, \;2.5} ) }}} $$

where the response variable, C, is a binary indicator of cousin marriage as above, and the predictor variable, Q, is the number of cousins available as potential partners for each individual (i.e. for individual i). The slope coefficient on number of cousins (β 2) reflects how the frequency of cousin marriage covaries with number of cousins.

In order to understand if cousin marriage among the Tsimane is determined by the (median-centred) frequency of cousins in the pool of potential partners of each individual (question 2.2), and whether Tsimane people marry cousins more than expected by chance in relation to their available pool of potential partners (question 2.3), we estimate the association between the individuals’ marriage type and the frequency of cousins in their pool of potential partners. We build a logistic regression specified by the following set of equations:

$$\matrix{ {C_{[ i ] }}\sim {{\rm Bernoulli}( {\theta_{[ i] }} ) } \cr {{\rm logit}( {\theta_{[ i] }} ) = \beta _1 } + {\beta _2F_{[ i ] } } \cr \beta \sim {{\rm Normal}( {0, \;2.5} ) } \cr } $$

where the predictor variable, F, is the frequency of cousins in the pool of potential partners of each individual (i.e. of individual i). From this estimation, we extract the slope coefficient – i.e. the β 2 coefficient – to see the effect of the frequency of cousins in the pool of potential partners on the chance to marry a cousin (question 2.2). We extract the intercept – more specifically, the inverse logit of β 1 – to test if Tsimane marry cousins more than expected by chance (question 2.3).

2.6. Analyses on the change in the pool of potential eligible partners over time

In order to see if the frequency of cousins in the pool of potential partners has changed over time (question 3.1), we estimate the linear association between the frequency of cousins in the pool of potential partners of each individual and the (median-centred) year of birth of participating individuals. We build a linear regression specified by the following set of equations:

$$\matrix{ {F_{[ i ] }}\sim {{\rm Normal}( \mu _{[ i] }, \;\sigma ) } \cr {\mu_{[ i ] } = \beta_1 } + {\beta_2 B_{[ i ] } } \cr } $$
$$\beta \sim {\rm Normal}( 0, \;2.5) $$
$$\sigma \sim {\rm Exponential}( 1 ) $$

where the response variable F is the observed frequency of cousins eligible as partners and the predictor variable, B, is the annual birth cohort of each individual. Here, σ is a standard deviation parameter. We also calculate the difference in μ for the maximum and minimum value of B.

In order to understand whether the change in frequency of cousin marriage over time is determined by the change in the frequency of cousins available as partners for individuals born in different years (question 3.2), we estimate the association between marriage type and the frequency of cousins in the pool of potential partners of each individual, controlling for annual birth cohort. We build a binomial logistic regression specified by the following set of equations:

$$\matrix{ {C_{[ i ] }}\sim {{\rm Bernoulli}( {\theta_{[ i] }} ) } \cr {{\rm logit}( {\theta_{[ i] }} ) = \beta _1 + \beta _2F_{[ i] }} + { \beta _3B_{[ i ] }} \cr } $$
$$\beta \sim \;{\rm Normal}( 0, \;\;2.5) $$

where the response variable, C, is a binary indicator of cousin marriage. The parameter vector β gives the intercept and slopes indicating the effect of the frequency of cousins as available partners and the effect of annual birth cohort on the response variable.

2.7. Analyses on fitness costs and benefits of cousin marriage

In order to assess if individuals married to cousins experience costs and/or benefits in terms of offspring survival (question 4.1), we estimate age-specific survival probabilities (i.e. for each age x, the probability of survival to age x + 1, given that individuals survived to age x), for individuals born to parents who are cousins and for individuals born to parents who are non-cousins. For these analyses, we include individuals born from 1870 to 2014. We build aggregated binomial regressions with varying effects on age categories and type of marriage, specified by the following set of equations:

$$S_{[ t, \ a ] } \sim \;{\rm Binomial}( N_{[ t, \ a ] } , \;\theta _{[ t, \ a ] } ) $$

where S [t,a ] is the number of individuals surviving for each parental marriage type, t, and age category, a, N [t,a] is the number of individuals at risk and θ [t,a] is the survival probability for each parental marriage type and age class. We then model:

$${\rm logit}( \theta _{[ t, \ a ] } ) = \mu_{[ t ] } + \kappa_{[ t, \ a ] } $$

where μ [t ] is an intercept for given parental marriage type and κ [t,a ] is an offset for each age category.

$$\mu _{[ t ] }\sim {\rm Normal}( \bar\mu , \;\sigma ) $$

The hyperparameters μ and σ partially pool information between parental marriage types.

$$\bar\mu \ \sim {\rm Normal}( 0, \;2.5) $$
$$\sigma \sim {\rm Exponential}( 1 ) $$

To define the age-category offsets, we use a Gaussian Process approach:

$$\kappa_{[ t, \ a ] } = L_{[ t ] } \ast \;\hat{\kappa }_{[ t, \ a ] } $$

where $\hat{\kappa }_{[ t, a ]}$ are unit normal random effects and L [t] is a Cholesky factor from the decomposition of the covariance matrix, Ω[t ], which describes variation and covariation among posterior random effects for each age category, for marriage type t (McElreath, Reference McElreath2020). The covariance matrix is defined using an exponential decay function in Stan. This model structure allows for the probability of survival to take on arbitrary functional forms with respect to age, while still sharing information across neighbouring age categories. More specifically, the Gaussian process approach: (1) allows for partially pooled estimation of parameters; (2) imposes no a priori functional form (e.g. linear, quadratic, etc.) on random effects; and (3) allows for a reduction in parameter complexity of the model by reducing the extent to which neighbouring random effects on age category can vary independently. The Cholesky decomposition formulation is simply a numerically efficient way to parameterise the Gaussian Process model (interested readers should refer to Barnard et al., Reference Barnard, McCulloch and Meng2000, for technical details).

We extract the posterior distributions of age-specific survival probabilities, for individuals born to cousins and born to non-cousins, and we use them to reconstruct survival curves. We compute child survival as the product of age-specific survival coefficients up to age 5. We compute the contrasts in child survival between those individuals born to cousin and those individuals born to non-cousins.

We use similar Stan models to assess if individuals married to cousins have costs and/or benefits in terms of fertility (question 4.2) and age at first reproduction (question 4.3). In these cases, the Stan model is used to estimate the probability of reproducing at each age of life for individuals married to partners who are cousins and for individuals married to partners who are non-cousins. The variable S [t,a ] is therefore replaced with the variable R [t,a ], which defines the number of individuals reproducing for each marriage type, t, and age category, a.

To estimate age-specific fertility, and age at first reproduction, we extract the posterior distributions of the probability of reproducing at each age, for individuals married to cousins and individuals married to non-cousins, and we use them to reconstruct the relevant fertility curves. Overall fertility here is defined as the sum of age-specific fertility across reproductive ages. The age at first reproduction here is defined as the first age when the expected number of births exceeds 1. We compute contrasts in fertility and age at first reproduction between each marriage type.

2.8. Analyses of temporal trends in fitness costs and benefits of cousin marriage

Here, we measure longitudinal changes in offspring survival, fertility, and age at first reproduction. In order to assess how the chance of surviving up to age 5 has changed over time for offspring of cousin and non-cousin parents (question 5.1), we estimate the chance to survive up to age 5 by parental marriage type for individuals born in each decade, D, between 1950 and 2000. We build a fixed-effect logistic regression specified by the following set of equations:

$$\matrix{ {S_{[ i ] } } \sim {{\rm Bernoulli}( {\theta_{[ i] }} ) } \cr {{\rm logit}( \theta _{[ i] } ) = \alpha_{[ {D( i ) } ] } } + {\beta_{[ {D( i ) } ] } \hat{C}_{[ i ] } } \cr } $$
$$\alpha _{[ d ] }\sim {\rm Normal}( 0, \;2.5) $$
$$\beta _{[ d ] }\sim {\rm Normal}( 0, \;2.5) $$

where the response variable, S, is a binary indicator for survival up to age 5, ${\hat C}_{[ i ] }$ is an indicator for being born to cousins, α is a vector of intercept parameters for each decadal birth cohort and β is a vector of slopes giving the estimated difference in log-odds of survival up to age 5 between individuals born to cousins and individuals born to non-cousins, for each decadal cohort. Here, and below, the function D(i) returns the birth decade of individual i. We fit this model independently for focal male and female individuals.

In order to assess how the fertility of individuals married to cousins and individuals married to non-cousins changed over time (question 5.2), we estimate longitudinal changes in fertility (for individuals born in decades between 1940 and 2000), using Poisson regression specified by the following set of equations:

$$R_{[ i ] }\sim {\rm Poisson}( {\lambda_{[ i] }} ) $$
$$\matrix{ {{\rm log}( \lambda_{[ i ] } ) = \alpha_{[ {D( i ) } ] } } + {\beta_{[ 1 , D( i ) ]} \tilde{C}_{[ i ] } + \beta_{[ 2 , D( i ) ]} E_{[ i ] } } \cr {\alpha _{[ d ] }\sim {\rm Normal}( 0, \;2.5) } \cr } $$
$$\beta _{[ d ] }\sim {\rm Normal}( 0, \;2.5) $$

where the response variable R is the observed number of children of each individual, α is a vector of intercept parameters for each decadal birth cohort, β 1 is a vector of slope parameters for the effect of spousal cousin marriage, C˜, on fertility for each decadal birth cohort, and β 2 is a vector of slope parameters for the effect of exposure time, E, on fertility for each decadal birth cohort. We control for exposure time (i.e. number of reproductive years) to adjust for individuals who have not completed their reproductive life. We fit this model independently for male and female individuals.

In order to assess how the age at first reproduction of individuals married to cousins and individuals married to non-cousins changed over time (question 5.3), we estimate age at first birth (for individuals born in decades between 1940 and 2000) as a function of decadal birth cohort and marriage type, using a Poisson regression specified by the following set of equations:

$$\matrix{ {A_{[ i ] } } \sim {{\rm Poisson}( {\lambda_{[ i ] } } ) } \cr {{\rm log}( \lambda_{[ i ] } ) = \alpha_{[ {D( i ) } ] } } + { \beta_{[ {D( i ) } ] } \tilde{C}_{[ i ] } } } $$
$$\alpha _{[ d ] }\sim {\rm Normal}( 0, \;2.5) $$
$$\beta _{[ d ] }\sim {\rm Normal}( 0, \;2.5) $$

where the response variable A is the observed age at first birth of each individual and other terms are as defined above. We fit this model independently for male and female individuals.

2.9. Results computation

For each model of interest, we used the rethinking package in R (McElreath, Reference McElreath2020) to compute and analyse the mean and the compatibility intervals (89% PI, percentile interval) of the posterior distribution of each model parameter. We also computed the Type S error, ps – the proportion of the posterior distribution on the far side of 0 – which gives the probability that an estimate's sign is reversed (Gelman & Carlin, Reference Gelman and Carlin2014). We interpret values of ps close to 0 (i.e. ps < 0.05) as indicative of a reliable non-zero effect, while values closer to 1 indicate a less reliable effect.

3. Results

In the following sub-sections, we describe our main results. Table 2 (questions 1–3) and Figure 5 (questions 4 and 5) present a summary of these results. Additional results are provided in the Supplementary Materials, Tables S1–S20.

Table 2. Summary of results on frequency of cousin marriage and pool of potential partners. Total sample size of married adults for all these analyses n = 1310. The dashed lines separate analyses conducted with different statistical models. The parameters extracted for questions 2.2 and 2.3 are, respectively, the slope (β) and intercept (θ) coefficients related to the same predictor from the same model. The parameters extracted for question 3.2 pertain to the effects of two predictors in the same model. The estimates for questions 2.2, 2.3, 3.1 and 3.2 are from analyses with standardised predictors

3.1. Frequency of cousin marriage over time

The frequency of cousin marriage among Tsimane has slightly decreased over time, from 29% of individuals born in 1964 to 22% of individuals born in 1993, β: −0.01 (89% PI −0.03, 0.00; ps= 0.091). See Figure 1 for a visualisation of the trend.

Figure 1. Observed (circles) and predicted (shaded region) frequency of cousin marriage (y axis) for birth-years between 1964 and 1993 (x axis). The shaded region represents the 89% compatibility interval of the predicted distribution. The predicted mean estimate slightly declined from 29% in 1964 to 22% in 1993. The circles represent annual birth-year ‘bins’, within which the empirical frequency of cousin marriage is displayed for illustration purposes; all analyses are based on the individual level data, and the sample size of individuals combined in each bin is indicated by the scale of the circles (ranging from 13 to 88 individuals). The total sample here includes n = 1331 individuals.

3.2. Frequency of cousin marriage in relation to the pool of eligible partners

We find a positive association between the frequency of cousin marriage and the availability of cousins as eligible partners, both in terms of the absolute number of cousins – β: 0.26 (89% PI 0.22, 0.30; ps < 0.001, Figure 2a) – and the frequency of cousins in the pool of eligible partners of each individual – β: 0.98 (89% PI 0.85, 1.12; ps < 0.001, Figure 2b). The frequency of cousin marriage among Tsimane is also much higher than the frequency of cousins available at any realised level – e.g. we estimate a 23% chance of cousin marriage (θ = 0.23, 89% PI 0.21, 0.25; ps < 0.001) even among individuals for whom cousins constitute only 1% of the pool of eligible partners. See Table 2.

Figure 2. Observed (circles) and predicted (shaded region) frequency of cousin marriage (y axes) as a function of number of cousins available as potential partners (x axis, panel a) and as a function of the frequency of cousins in the pool of potential partners of each individual (x axis, panel b). The shaded regions represent the 89% compatibility intervals of the predicted frequencies. All analyses are based on the individual level data. For illustration, each circle represents a ‘bin’, an x-axis range in which we calculate frequency estimates for illustration purposes. The size of the circles is scaled by sample size (from 2 individuals, for the highest x-axis ranges, to 329 individuals, for the lowest). The dashed line—which has a slope of 1—represents a linear relationship between the frequency of cousin marriage and the frequency of cousins available as partners; it represents the expected level of cousin marriage under random mating.

3.3. Change in the pool of eligible partners over time

Individuals born more recently have a higher number of cousins and – to a greater extent – non-cousins available as eligible partners than individuals born in earlier cohorts. However, we find a negative association between the frequency of cousins in the pool of potential partners of a given individual and annual birth cohort, β: −0.21 (89% PI −0.26, −0.17; ps < 0.001). The expected difference in the frequency of cousins eligible as partners for individuals in the 1964 birth cohort and individuals in the 1993 birth cohort is −0.01 (89% PI −0.01, −0.01; ps < 0.001, Figure 3). In the model estimating the association between the frequency of cousin marriage and the frequency of cousins in the pool of potential partners of individuals, controlling for annual birth cohort, we find a positive association between the frequency of cousin marriage and the frequency of cousins as eligible partners, β: 1.02 (89% PI 0.89, 1.16; ps < 0.001). In this model, the association between the frequency of cousin marriage and annual birth cohort is lower and less reliable, β: 0.16 (89% PI 0.03, 0.28; ps: 0.027; see Table 2).

Figure 3. Observed (circles) and predicted (shaded region) frequency of cousins in the pool of potential partners of each individual (y axis) as a function of annual birth cohort (x axis). The shaded region represents the 89% compatibility interval. The predicted mean estimate of frequency of cousins declined from 1% in 1964 to 0.3% in 1993. Each circle represents a ‘bin’, an x-axis averaged frequency estimate for each birth year calculated for illustration purposes. The size of the circles is scaled by sample size (ranging from 13 to 88 individuals).

3.4. Age-specific offspring survival and fertility

Both male and female children born to cousins have similar probability of survival to age 5 compared with children born to non-cousin parents. The difference in survival probability as a function of parental cousin marriage is: for male children, 0.01 (89% PI −0.01, 0.04; ps = 0.161); for female children, −0.01 (89% PI −0.03, 0.02; ps = 0.328) (Figure 4b). Fertility across the life course is similar for men married to cousins and for men married to non-cousins. The difference in total predicted number of offspring for a man of age 50 is: −0.20 (89% PI −0.54, 0.15; ps = 0.178). Fertility is similar for women married to cousins and for women married to non-cousins. The difference in total predicted number of offspring for a woman of age 50 is −0.17 (89% PI −0.50, 0.19; ps = 0.224; Figure 4a). Age at first reproduction within cousin marriages is reliably earlier than within non-cousin marriages, both for men – expected difference in estimated years −0.75 (89% PI −1.00, 0.00; ps < 0.001) – and for women – expected difference: −0.35 (89% PI −1.00, 0.00; ps < 0.001).

Figure 4. Fertility and offspring survival as a function of cousin marriage. In panel (a) we plot expected number of children (y axis) for each age of life (x axis). In panel (b) we plot the posterior probabilities of cumulative survival (y axis) for each age of life from 0 to 5 (x axis). We plot estimates for female individuals in the left frames and male individuals in the right frames. Estimates are colored by marriage type categories. The bands represent 89% compatibility intervals and reflect the precision of the estimated mean.

3.5. Temporal trends in survival and fertility

We find no reliable evidence of differences in survival to age 5 as a function of parental marriage type, for any of the decadal birth cohorts between 1950 and 2000 (see Figure 5a). There are no strong differences in fertility between individuals married to cousins and individuals married to non-cousins in most decades. Only men married to cousins born between 1940 and 1950 and men married to cousins born between 1970 and 1980 have reliably higher fertility than men married to non-cousins born in the same decades. The expected difference in number of children for men born between 1940 and 1950 is 0.99 (89% PI 0.20, 1.87; ps = 0.025). The expected difference in number of children for men born between 1970 and 1980 is 0.56 (89% PI 0.13, 0.99; ps = 0.018). Also, women married to cousins and born between 1960 and 1970 have higher fertility than women married to non-cousins born in the same decades, with an expected difference in number of children of 0.68 (89% PI 0.13, 1.24; ps = 0.024). Slight differences in fertility as a function of cousin marriage apparent in early decades disappear in more recent decades, where we have more observations and thus statistical power (see Figure 5b). We find no strong evidence of temporal trends in age at first reproduction as a function of marriage type for either sex (see Figure 5c).

Figure 5. Temporal trends in the association between cousin marriage and fertility/survival. In panel (a), we plot offspring survival outcomes. In panel (b), we plot fertility outcomes. And, in panel (c), we plot age at first reproduction outcomes. Each horizontal bar represents the posterior distribution of the contrasts (i.e., differences) between families practicing cousin marriage and families not practicing it, in each decadal birth cohort (questions 5). Estimates for males are plotted with black circles, estimates for females are plotted with white circles. In the first rows, estimates are averaged across the whole sample (questions 4). The numbers represent the sample sizes of each decadal birth cohort. The first number is the sample size of individuals in families practicing cousin marriage, and the second number is the sample size of individuals in families not practicing cousin marriage.

4. Discussion

4.1. Overview

The frequency of cousin marriage has decreased only slightly over time among Tsimane people born between 1964 and 1993 (question 1). The probability of cousin marriage relates to the availability of eligible cousins, both in absolute numbers and relative to non-cousins in the marriage market. Consistent with stated Tsimane cultural preferences, cousin marriage occurs at higher rates than expected by chance (questions 2). While market integration and increasing exposure to Bolivian cultural norms may affect Tsimane marriage practices, we find that changes in the frequency of cousin marriage over time remain tightly connected to the frequency of cousins in the pool of eligible partners (questions 3). Individuals married to cousins exhibit only very slight differences in life histories, with earlier ages of first reproduction, but similar completed fertility and offspring child survival compared with those married to non-cousins (question 4). Analyses on temporal trends in fitness (questions 5) did not reveal clear differences across decadal birth cohorts as a function of marriage type in offspring survival, fertility and age at first reproduction. Minimal differences in fitness visible for individuals born in the early decades tend to disappear for individuals born in more recent decades. Overall, our findings support the prediction that young Tsimane people choose their spouses from a broader pool of potential partners compared to their ancestors. Demographic and socioeconomic shifts – e.g. those that shape kin and non-kin networks, and therefore fitness consequences of different marriage strategies – may lead individuals to use different partner selection criteria. Nevertheless, entrenched cultural traditions may counteract this change and keep the frequency of cousin marriage relatively stable in this population.

4.2. Demographic factors shaping the frequency of cousin marriage

The decline in the frequency of cousin marriage – from 29 to 22% – over birth cohorts spanning a period of 30 years has been remarkably modest. This finding recapitulates ethnographic reports describing Tsimane as particularly successful in preserving their cultural identity, for instance by continuing to speak and teach their language, and by maintaining traditional fishing and hunting practices. The slight decline in the frequency of Tsimane cousin marriage, however, accompanies a demographic shift that is reshaping social networks. Individuals born more recently may have more cousins available as potential partners than individuals born earlier, as a result of the increase in fertility among Tsimane people born from the 1930s through the 1970s (see supplement of: Costa et al., Reference Costa, Trumble, Kaplan and Gurven2018). In particular, a drop in mortality rates (Gurven et al., Reference Gurven, Kaplan and Supa2007) may have increased the number of cousins in each age class. However, the number of potential partners who are non-cousins increased at a much higher rate in recent years, as a consequence of the population growth rate (Lukas et al., Reference Lukas, Reynolds, Boesch and Vigilant2005; Gurven et al., Reference Gurven, Stieglitz, Trumble, Blackwell, Beheim, Davis and Kaplan2017). Accordingly, the frequency of cousins in the pool of eligible partners of a given individual has decreased over time (question 3.1). The slight change in the frequency of cousin marriage over time is in fact almost completely associated with the change in the frequency of cousins as eligible partners over time (question 3.2). The remaining effect of birth year, after we account for the changes in frequency of cousins, is negligible. Tsimane people make marriage decisions on the basis of partner availability. The frequency of cousin marriage is constrained by the availability of cousins as eligible partners, However, Tsimane marry cousins at rates higher than what would be expected if individuals chose partners randomly, purely according to the demographic composition of the pool of the eligible partners. The frequency of cousin marriage is in fact higher than what would be predicted from the frequency of marriageable cousins (Hajnal, Reference Hajnal1963) (question 2.3). According to our results, when individuals have more than 12 cousins as potential partners, more than 80% of them will marry a cousin. Another potential, although untested, explanation for the slight decrease in the frequency of cousin marriage could be heightened competition among same-sex siblings for marrying a cousin (Chagnon et al., Reference Chagnon, Lynch, Shenk, Hames and Flinn2017). This competition may arise from the increased number of siblings resulting from population growth. It is important to note, however, that the number of siblings is expected to increase linearly, much like the absolute number of cousins, and therefore any potential effect of this competition should be minimal, because the exponential increase in the number of first cousins is even larger in a growing population.

4.3. Fitness consequences of cousin marriage

Overall, the life-history consequences associated with cousin marriage in Tsimane are quite limited. The average relatedness among individuals in the whole population might be relatively high owing to endogamy and because many Tsimane individuals who are not married to a first cousin are nevertheless married to a more distant cousin (Swinford et al., Reference Swinford, Prall, Williams, Sheehama, Scelza and Henn2023). Socioeconomic factors may complicate the influence of direct biological factors in determining the modest differences between families practising different marriage strategies. This is due to the reported higher prevalence of cousin marriage in more remote areas (unpublished observations reported in: Patel et al., Reference Patel, Godoy, Seyfried, Reyes-García, Huanca, Leonard and Tanner2007). Fitness differences as a function of marriage type are minimal in the analyses for decadal birth cohorts, and disappear, in particular, for individuals born in more recent decades. The diminishing disparities found in fitness indicators suggest that the socioeconomic conditions of Tsimane people might be changing, possibly becoming more homogeneous across locations and environmental contexts, owing to increasing market integration and acculturation, changes in kinship and social networks, changes in residence and travel patterns, and increasing access to healthcare and contraception. Averaging offspring survival across years, we found no difference between sons and daughters born to cousins and sons and daughters born to non-cousins (question 4.1). In cases when children born to cousins live in remote villages where child mortality rates are higher than in acculturated villages closer to market towns (Gurven et al., Reference Gurven, Kaplan and Supa2007), we speculate that potential costs in child mortality from socioeconomic conditions are compensated for by benefits from increased kin support in families practising cousin marriage (Sear & Mace, Reference Sear and Mace2008; Willführ et al., Reference Willführ, Johow and Voland2018). We were unable to detect stable differences in offspring child survival by decadal birth cohort (question 5.1). This finding is in line with previous research reporting little change in child mortality across all the geographical and social contexts of the Tsimane territory over the second half of the twentieth century (Gurven et al., Reference Gurven, Kaplan and Supa2007). The effects of cousin marriage on fertility are very limited (question 4.2). Averaging fertility across the whole sample, we find no differences in fertility as a function of marriage type. As in the survival model, we were unable to detect stable differences in fertility as a function of marriage type by decadal birth cohort (question 5.2). The only differences we detected were higher fertility for cousin-married women born in the 1960s, and higher fertility for cousin-married men born in the 1940s and 1970s. The reproductive advantage found for individuals born in these periods may be explained by improved support from kin in families practising cousin marriage (Hooper et al., Reference Hooper, Gurven, Winking and Kaplan2015; Sear & Coall, Reference Sear and Coall2011) or norms for higher fertility among families practising cousin marriage (Mcallister et al., Reference Mcallister, Gurven, Kaplan and Stieglitz2012). However, differences in number of children as a function of marriage type become less marked for the most recent birth cohorts, for whom the estimates are more reliable because of the larger sample size. Kin support may have become less significant in influencing fertility for the younger generations, owing to shifts in kin-based and sharing networks. The fading of differences for individuals born in the decades from 1980 and onward could be also a consequence of a nascent demographic transition – with fertility reduction – across all the areas of the Tsimane territory (Gurven et al., Reference Gurven, Stieglitz, Trumble, Blackwell, Beheim, Davis and Kaplan2017). We observed that both women and men marrying cousins may gain slight fitness benefits in terms of earlier age at first reproduction (question 4.3). This suggests that cousin marriage might ease the problem of finding a mate in this population, especially in the more remote villages where the number of eligible partners can be limited (Cavalli-Sforza et al., Reference Cavalli-Sforza, Moroni and Zei2013). Increased parental control in marriage choices might also be determining the earlier age at first reproduction for individuals in cousin marriage among the Tsimane. In demographic interviews conducted between 2002 and 2005 (M. Gurven, unpublished observations), arranged marriages were four times more likely than non-arranged marriages if the spouses belonged to the cross-lineage. The opportunity to get married earlier can be an incentive to marry a cousin for Tsimane: prestigious Tsimane men marry wives who give birth at earlier ages (Von Rueden et al., Reference Von Rueden, Gurven and Kaplan2011). Early age at reproduction is in fact the strongest pathway between fitness and status in Tsimane (Von Rueden et al., Reference Von Rueden, Gurven and Kaplan2011). As for the temporal change in trends of age at first reproduction (question 5.3), slight differences found for individuals born in the early decades disappear for individuals born in the most recent decades, probably as a result of the increased availability of non-cousin partners in the more recent period or the reduced risk of adverse outcomes of giving birth at a young age.

4.4. Study limitations

The genealogical data used in the current study provide information only on individuals’ birth and death dates and their parents’ identity (see the ‘Methods’ section), because only these data span across several decades. Data do not contain information on individuals’ village of provenance, health or wealth. These details would have been valuable for investigating the impact of socioeconomic factors on the frequency of cousin marriage and its effects on fitness. Another potential data limitation is the absence of information regarding marriages, such as when the couples were married and whether the unions were arranged or not. Our method of constructing family unions is based on shared parentage for a given child (see the ‘Methods’ section). Consequently, we do not include couples who never had children, which could affect fitness estimates if childless couples were more or less likely to be cousins. Moreover, extra-marital offspring may inflate the number of apparent marriages. Because we lack information on the great-grandparents for most couples, we cannot discern patterns related to second cousin marriage, which is also frequent among Tsimane. Nevertheless, the dynamics we identified in first-cousin marriage practice and kinship structure can provide insights into potential dynamics in other types of kin marriages. Our sample is biased towards young individuals (especially in the offspring groups), because we need to know the grandparents of individuals to infer whether they married a cousin or not. For individuals born recently, we do not have data to infer survival and fertility across the lifespan. This limits our ability to use age-structured population projection models to investigate the potential lineage consequences of trade-offs between survival and reproduction (Dalzero et al., Reference Dalzero, Ross and Lukas2023). As for the analyses on temporal change in fitness, we could not make strong inferences about individuals born in the decades from 1900 to 1940, because of the low sample size. We used a strict relatedness classification (i.e. including only individuals with eight grandparents known) for the analyses on change in frequency of cousin marriage in relation to the available partners, and a relaxed relatedness classification (i.e. discarding only individuals with eight grandparents unknown) for the analyses on fitness and temporal change in fitness over time. The rationale for this decision was to strike a balance between avoiding the possible confound of classifying people as not having cousins as available partners because their grandparents are unknown, while increasing sample size in the fitness analyses. Some individuals in cousin marriages could be incorrectly classified as in non-cousin marriages under the relaxed condition. Nevertheless, we conducted robustness checks and found that our results qualitatively hold using either approach for the comparisons with a sufficient sample size (code and output files are provided in the repository). The sample might also be biased because genealogical data is likely to mainly include individuals who survived and reproduced, and therefore appear as ancestors. However, this bias applies to both families practising cousin marriage and families not practising it and therefore should not influence our estimation of the differences in survival and fertility.

4.5. Conclusion and future directions

Our study of the dynamics of cousin marriage among Tsimane highlights the importance of changes in the composition of the pool of eligible partners. In recent decades, the Tsimane have been choosing their spouses from a broader pool of potential partners. This finding is presumably due to the high population growth rate of Tsimane and the migration of Tsimane into larger and less isolated settlements. However, cousin marriage is practised far more frequently among Tsimane than expected by chance. Cousin marriage among Tsimane is a preferred marriage choice, and this preference, in the examined decades, has been relatively stable. The pace of change may be countered by cultural inertia, as observed in other contexts (e.g. Bashir & Nazir, Reference Bashir and Nazir2022; Shenk et al., Reference Shenk, Towner, Voss and Alam2016). Cousin marriage entails very modest or negligible fitness consequences in terms of offspring survival and fertility. In contrast, cousin marriage seems to have given Tsimane people some advantages on the marriage market, easing the marriage process in this population, as suggested by the earlier age at first reproduction of cousin-married individuals. However, any detectable differences in the examined traits between individuals practising cousin marriage and not practising it have disappeared closer to the present. This reduction in fitness differences by marriage type suggests a homogenization owing to increasing acculturation, healthcare access and the use of contraception, and social connectedness across the Tsimane territory. Accordingly, the factors that have guided marriage choice and led to a strong preference for cousins as spouses (e.g. increased marriage opportunities and kin-based support for those choosing cousins) might also be changing. This change might contribute to a sharper reduction in the frequency of cousin marriages in the future, as has been observed in other societies undergoing market integration (Bittles & Black, Reference Bittles and Black2010).

In order to better understand dynamics in the frequency of cousin marriage among Tsimane, future research could explore the impact of cousin marriage on social networks, economic production and health. One interesting angle could be to investigate how cousin marriage can persist at rates higher than what we could expect considering the demographic and socioeconomic changes. This persistence is also observable in other cultures, including immigrants adapting to entirely new socioeconomic conditions (Reniers, Reference Reniers2001; Shaw, Reference Shaw2001). An analysis of the regional variation in cousin marriage across the Tsimane territory, which quantifies its frequency both near market towns and in the most remote areas, would help to investigate the socioeconomic correlates of this practice. Further research could also explore how cultural traits such as religious beliefs determine marital choices. Our analyses reveal the complexity of the potential connections between cultural changes in marriage practices and the underlying societal changes that shape partner availability and fitness returns to different marriage choices. More broadly, our approach showcases a path to use longitudinal data to gain insights into how sociodemographic changes can lead to behavioural changes in human societies.

Acknowledgements

The Department of Human Behavior, Ecology and Culture at the Max Plank Institute for Evolutionary Anthropology provided valuable feedback and support on this analysis. We thank James Holland Jones for input and inspiration on demographic factors affecting cousin marriage. We thank two anonymous reviewers for feedback and suggestions.

Author contributions

AD and DL designed the structure of the study. AD prepared data and constructed statistical models, with DL and BAB's guidance. BAB designed Stan models, organised data, and edited workflow and manuscript. AD wrote the article, with input from DL. MG and HK coordinated data collection and edited the manuscript. JS and PLH contributed to the data collection and edited the manuscript. CTR edited the manuscript and provided statistical consulting. All the authors have read the article.

Financial support

This research received support from the Department of Human Behavior, Ecology, and Culture at the Max Planck Institute for Evolutionary Anthropology. THLHP data collection is supported by the National Science Foundation (BCS0136274) and National Institutes of Health (R01AG024119). JS acknowledges IAST (Institute for Advanced Study in Toulouse) funding from the French National Research Agency under the Investments for the Future (Investissements d'Avenir) programme, grant ANR-17-EURE-0010.

Conflicts of interest

The authors declare no competing interests.

Ethics

All data collection procedures were approved by institutional review boards at the University of California – Santa Barbara and the University of New Mexico, and permissions were obtained from the Gran Consejo Tsimane, community leaders and study participants.

Research transparency, reproducibility and data availability

Code for reproducing our methods and results can be accessed on the Open Science Framework: https://doi.org/10.17605/OSF.IO/PK8F7. Individual-level data are stored in the Tsimane Health and Life History Project (THLHP) Data Repository, and are available through controlled access protocols (for ethical reasons related to data sovereignty). The THLHP's highest priority is the safeguarding of human subjects and minimisation of risk to study participants. The THLHP adheres to the CARE Principles for Indigenous Data Governance, which assure that the Tsimane: (1) have sovereignty over how data are shared; (2) are the primary gatekeepers determining ethical use; (3) are actively engaged in the data generation; and (4) derive benefit from data generated and shared whenever possible. The THLHP is also committed to the FAIR Principles to facilitate data use. Requests for individual-level data should take the form of an application that minimally details the exact uses of the data and the research questions to be addressed, procedures that will be employed for data security and individual privacy, potential benefits to the study communities, and procedures for assessing and minimising stigmatising interpretations of the research results (see the following webpage for links to the data sharing policy and data request forms: https://tsimane.anth.ucsb.edu/data.html). Requests for individual-level data will require institutional approval and will be reviewed by an Advisory Council composed of tribal leaders, tribal community members, Bolivian scientists and the THLHP leadership. The study authors and the THLHP leadership are committed to open science and are available to assist interested investigators in preparing data access requests.

Supplementary material

To view supplementary material for this article, please visit https://doi.org/10.1017/ehs.2024.11

References

Assaf, S., & Khawaja, M. (2009). Consanguinity trends and correlates in the Palestinian territories. Journal of Biosocial Science, 41(1), 107124.CrossRefGoogle ScholarPubMed
Bailey, D. H., Hill, K. R., & Walker, R. S. (2014). Fitness consequences of spousal relatedness in 46 small-scale societies. Biology Letters, 10(5), 20140160.CrossRefGoogle ScholarPubMed
Barakat, B., & Basten, S. (2014). Modelling the constraints on consanguineous marriage when fertility declines. Demographic Research, 30, 277312.CrossRefGoogle Scholar
Barnard, J., McCulloch, R., & Meng, X. -L. (2000). Modeling covariance matrices in terms of standard deviations and correlations, with application to shrinkage. Statistica Sinica, 10(4), 12811311.Google Scholar
Bashir, S., & Nazir, S. (2022). Thirty years of consanguineous marriages in Pakistan. NUST Journal of Social Sciences and Humanities, 8(3), 1732.CrossRefGoogle Scholar
Bittles, A. H., & Black, M. L. (2010). Consanguineous marriage and human evolution. Annual Review of Anthropology, 39, 193207.Google Scholar
Bittles, A. H., Grant, J., Sullivan, S., & Hussain, R. (2002). Does inbreeding lead to decreased human fertility? Annals of Human Biology, 29(2), 111130.CrossRefGoogle ScholarPubMed
Borgerhoff Mulder, M. (1998). The demographic transition: Are we any closer to an evolutionary explanation? Trends in Ecology & Evolution, 13(7), 266270.Google Scholar
Borgerhoff Mulder, M., & Ross, C. T. (2019). Unpacking mating success and testing Bateman's principles in a human population. Proceedings of the Royal Society B, 286(1908), 20191516.Google Scholar
Borgerhoff Mulder, M., Bowles, S., Hertz, T., Bell, A., Beise, J., Clark, G., … Wiessner, P. (2009). Intergenerational wealth transmission and the dynamics of inequality in small-scale societies. Science, 326(5953), 682688.CrossRefGoogle ScholarPubMed
Bras, H., Van Poppel, F., & Mandemakers, K. (2009). Relatives as spouses: Preferences and opportunities for kin marriage in a western society. American Journal of Human Biology: The Official Journal of the Human Biology Association, 21(6), 793804.CrossRefGoogle Scholar
Brennan, E. R. (1981). Kinship, demographic, social, and geographic characteristics of mate choice in Sanday, Orkney islands, Scotland. American Journal of Physical Anthropology, 55(1), 129138.CrossRefGoogle Scholar
Caldwell, J. C., Reddy, P. H., & Caldwell, P. (1983). The causes of marriage change in South India. Population studies, 37(3), 343361.Google Scholar
Carpenter, B., Gelman, A., Hoffman, M. D., Lee, D., Goodrich, B., Betancourt, M., …, Riddell, A. (2017). Stan: A probabilistic programming language. Journal of Statistical Software, 76(1), 1–32.CrossRefGoogle Scholar
Cavalli-Sforza, L. L., Moroni, A., & Zei, G. (2013). Consanguinity, inbreeding, and genetic drift in Italy (mpb-39). Princeton University Press.Google Scholar
Cazes, M. (1996). Demographic anthropology in the study of matrimonial exchanges in a Dogon isolate. International Journal of Anthropology, 11(2), 6171.CrossRefGoogle Scholar
Chagnon, N. A., Lynch, R. F., Shenk, M. K., Hames, R., & Flinn, M. V. (2017). Cross-cousin marriage among the Yanomamö shows evidence of parent–offspring conflict and mate competition between brothers. Proceedings of the National Academy of Sciences, 114(13), E2590E2607.CrossRefGoogle ScholarPubMed
Clark, D. W., Okada, Y., Moore, K. H., Mason, D., Pirastu, N., Gandin, I., … Wilson, J. F. (2019). Associations of autozygosity with a broad range of human phenotypes. Nature Communications, 10(1), 117.CrossRefGoogle ScholarPubMed
Costa, M. E., Trumble, B., Kaplan, H., & Gurven, M. D. (2018). Child nutritional status among births exceeding ideal family size in a high fertility population. Maternal & Child Nutrition, 14(4), e12625.CrossRefGoogle Scholar
Dalzero, A., Ross, C. T., & Lukas, D. (2023). Fitness consequences of cousin marriage: A life-history assessment in two populations. Evolutionary Human Sciences, 5, e3.Google Scholar
Davison, R. J., & Gurven, M. D. (2021). Human uniqueness? Life history diversity among small-scale societies and chimpanzees. PLoS ONE, 16(2), e0239170.Google Scholar
Ellis, R. (1997). A taste of movement: An exploration of the social ethics of the Tsimanes of lowland Bolivia. PhD thesis, The University of St Andrews.Google Scholar
Gelman, A., & Carlin, J. (2014). Beyond power calculations: Assessing type s (sign) and type m (magnitude) errors. Perspectives on Psychological Science, 9(6), 641651.CrossRefGoogle ScholarPubMed
Ghosh, A., Hwang, S. I. M., & Squires, M. (2023). Economic consequences of kinship: Evidence from US bans on cousin marriage. The Quarterly Journal of Economics, 138(4), 25592606.CrossRefGoogle Scholar
Givens, B. P., & Hirschman, C. (1994). Modernization and consanguineous marriage in Iran. Journal of Marriage and the Family, 820834.CrossRefGoogle Scholar
Godoy, R., Reyes-García, V., Seyfried, C., Huanca, T., Leonard, W. R., McDade, T., …, Vadez, V. (2007). Language skills and earnings: Evidence from a pre-industrial economy in the Bolivian Amazon. Economics of Education Review, 26(3), 349360.CrossRefGoogle Scholar
Godoy, R., Eisenberg, D. T., Reyes-García, V., Huanca, T., Leonard, W. R., McDade, T. W., … Tanner, S. (2008). Assortative mating and offspring well-being: Theory and empirical findings from a native Amazonian society in Bolivia. Evolution and Human Behavior, 29(3), 201210.Google Scholar
Goode, W. J. (1963). World revolution and family patterns. Free Press Glencoe.Google Scholar
Gurven, M., & Kaplan, H. (2007). Longevity among hunter-gatherers: A cross-cultural examination. Population and Development Review, 33(2), 321365.CrossRefGoogle Scholar
Gurven, M., Kaplan, H., & Supa, A. Z. (2007). Mortality experience of Tsimane Amerindians of Bolivia: Regional variation and temporal trends. American Journal of Human Biology: The Official Journal of the Human Biology Association, 19(3), 376398.CrossRefGoogle ScholarPubMed
Gurven, M., Stieglitz, J., Trumble, B., Blackwell, A. D., Beheim, B., Davis, H., …, Kaplan, H. (2017). The Tsimane health and life history project: Integrating anthropology and biomedicine. Evolutionary Anthropology: Issues, News, and Reviews, 26(2), 5473.Google Scholar
Hajnal, J. (1963). Concepts of random mating and the frequency of consanguineous marriages. Proceedings of the Royal Society of London. Series B. Biological Sciences, 159(974), 125177.Google ScholarPubMed
Helgason, A., Palsson, S., Guðbjartsson, D. F., Kristjansson, þ., & Stefansson, K. (2008). An association between the kinship and fertility of human couples. Science, 319(5864), 813816.CrossRefGoogle ScholarPubMed
Hoben, A. D., Buunk, A. P., & Fisher, M. L. (2016). Factors influencing the allowance of cousin marriages in the standard cross cultural sample. Evolutionary Behavioral Sciences, 10(2), 98.CrossRefGoogle Scholar
Hooper, P. L., Gurven, M., Winking, J., & Kaplan, H. S. (2015). Inclusive fitness and differential productivity across the life course determine intergenerational transfers in a small-scale human society. Proceedings of the Royal Society B: Biological Sciences, 282(1803), 20142808.Google Scholar
Hussain, R., & Bittles, A. H. (2004). Assessment of association between consanguinity and fertility in asian populations. Journal of Health, Population and Nutrition, 112.Google Scholar
Husain, M. A., & Bunyan, M. A. (1997). Consanguineous marriages in a Saudi population and the effect of inbreeding on prenatal and postnatal mortality. Annals of Tropical Paediatrics, 17(2), 155160.CrossRefGoogle Scholar
Islam, M. M. (2018). The changing pattern and determinants of declining consanguinity in Jordan during 1990–2012. Annals of Human Biology, 45(2), 140147.CrossRefGoogle ScholarPubMed
Johow, J., Willführ, K. P., & Voland, E. (2019). High consanguinity promotes intergenerational wealth concentration in socioeconomically privileged Krummhörn families of the 18th and 19th centuries. Evolution and Human Behavior, 40(2), 204213.CrossRefGoogle Scholar
Kalam, M. A., Sharma, S. K., Ghosh, S., & Roy, S. (2021). Change in the prevalence and determinants of consanguineous marriages in India between national family and health surveys of 1992–1993 and 2015–2016. Human Biology, 92(2), 93113.Google Scholar
Kaplan, H., Hooper, P. L., Stieglitz, J., & Gurven, M. (2015). The causal relationship between fertility and infant mortality. Population in the Human Sciences: Concepts, Models, Evidence, 5, 361376.Google Scholar
Koc, I. (2008). Prevalence and sociodemographic correlates of consanguineous marriages in Turkey. Journal of Biosocial Science, 40(1), 137148.CrossRefGoogle ScholarPubMed
Kumari, N., Bittles, A. H., & Saxena, P. (2020). Has the long-predicted decline in consanguineous marriage in India occurred? Journal of Biosocial Science, 52(5), 746755.Google Scholar
Lévi-Strauss, C. (1971). The elementary structures of kinship, no. 340. Beacon Press.Google Scholar
Lukas, D., Reynolds, V., Boesch, C., & Vigilant, L. (2005). To what extent does living in a group mean living with kin? Molecular Ecology, 14(7), 21812196.CrossRefGoogle Scholar
Martin, M. A., & Gurven, M. (2022). Traditional and biomedical maternal and neonatal care practices in a rural indigenous population of the Bolivian Amazon. Global Public Health, 17(6), 971985.CrossRefGoogle Scholar
Mattison, S. M., & Sear, R. (2016). Modernizing evolutionary anthropology: Introduction to the special issue. Human Nature, 27, 335350.CrossRefGoogle ScholarPubMed
Mcallister, L., Gurven, M., Kaplan, H., & Stieglitz, J. (2012). Why do women have more children than they want? Understanding differences in women's ideal and actual family size in a natural fertility population. American Journal of Human Biology, 24(6), 786799.Google Scholar
McElreath, R. (2020). Statistical rethinking: A Bayesian course with examples in R and Stan. Chapman and Hall/CRC.Google Scholar
Micheletti, A. J. (2020). Modelling cultural selection on biological fitness to integrate social transmission and adaptive explanations for human behaviour. Evolutionary Human Sciences, 2.CrossRefGoogle ScholarPubMed
Miner, E. J., Gurven, M., Kaplan, H., & Gaulin, S. J. (2014). Sex difference in travel is concentrated in adolescence and tracks reproductive interests. Proceedings of the Royal Society B: Biological Sciences, 281(1796), 20141476.Google Scholar
Mobarak, A. M., Chaudhry, T., Brown, J., Zelenska, T., Khan, M. N., Chaudry, S., …, Li, S. (2019). Estimating the health and socioeconomic effects of cousin marriage in south asia. Journal of Biosocial Science, 51(3), 418435.Google Scholar
Murphy, M. (2011). Long-term effects of the demographic transition on family and kinship networks in Britain. Population and Development Review, 37, 5580.Google Scholar
Patel, A. M., Godoy, R. A., Seyfried, C., Reyes-García, V., Huanca, T., Leonard, W. R., …, Tanner, S. (2007). On the accuracy of perceived parental height in a native amazonian society. Economics & Human Biology, 5(1), 165178.Google Scholar
Peña, J., Alfonso-Sánchez, M. A., & Calderón, R. (2002). Inbreeding and demographic transition in the Orozco valley (Basque country, Spain). American Journal of Human Biology, 14(6), 713720.CrossRefGoogle ScholarPubMed
Pettener, D. (1985). Consanguineous marriages in the upper Bologna Appennine (1565–1980): Microgeographic variations, pedigree structure and correlation of inbreeding secular trend with changes in population size. Human Biology, 57(2), 267288.Google Scholar
Pisor, A. C., & Gurven, M. (2018). When to diversify, and with whom? Choosing partners among out-group strangers in lowland Bolivia. Evolution and Human Behavior, 39(1), 3039.CrossRefGoogle Scholar
Core Team, R (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing.Google Scholar
Reniers, G. (2001). The post-migration survival of traditional marriage patterns: Consanguineous marriages among Turks and Moroccans in Belgium. Journal of Comparative Family Studies, 32(1), 2145.CrossRefGoogle Scholar
Reyes-García, V., Paneque-Gálvez, J., Luz, A., Gueze, M., Macía, M., Orta-Martínez, M., & Pino, J. (2014). Cultural change and traditional ecological knowledge: An empirical analysis from the Tsimane'in the Bolivian Amazon. Human Organization, 73(2), 162173.Google Scholar
Ross, C. T., Borgerhoff Mulder, M., Oh, S. -Y., Bowles, S., Beheim, B., Bunce, J.Ziker, J. (2018). Greater wealth inequality, less polygyny: Rethinking the polygyny threshold model. Journal of the Royal Society interface, 15(144), 20180035.Google Scholar
Sahoo, H., Debnath, P., Mandal, C., Nagarajan, R., & Appunni, S. (2022). Changing trends of consanguineous marriages in South India. Journal of Asian and African Studies, 57(2), 209225.CrossRefGoogle Scholar
Schaffnit, S., Page, A. E., Lynch, R., Spake, L., Sear, R., Sosis, R.,…Shenk, M. K. (2022). The impact of market integration on arranged marriages in Matlab, Bangladesh . Evolutionary Human Sciences, 5, e5.Google Scholar
Schulz, J. F., Bahrami-Rad, D., Beauchamp, J. P., & Henrich, J. (2019). The church, intensive kinship, and global psychological variation. Science, 366(6466), eaau5141.CrossRefGoogle ScholarPubMed
Sear, R., & Coall, D. (2011). How much does family matter? Cooperative breeding and the demographic transition. Population and Development Review, 37, 81112.Google Scholar
Sear, R., & Mace, R. (2008). Who keeps children alive? A review of the effects of kin on child survival. Evolution and Human Behavior, 29(1), 118.CrossRefGoogle Scholar
Shaw, A. (2001). Kinship, cultural preference and immigration: Consanguineous marriage among British Pakistanis. Journal of the Royal Anthropological Institute, 7(2), 315334.CrossRefGoogle Scholar
Shaw, A., & Raz, A. E. (2015). Cousin marriages: Between tradition, genetic risk and cultural change, vol. 28. Berghahn Books.Google Scholar
Shenk, M. K., Towner, M. C., Voss, E. A., & Alam, N. (2016). Consanguineous marriage, kinship ecology, and market transition. Current Anthropology, 57(S13), S167S180.CrossRefGoogle Scholar
Team, Stan Development (2021). RStan: The R interface to Stan. R package version 2.21.3.Google Scholar
Swinford, N., Prall, S., Williams, C. M., Sheehama, J., Scelza, B., & Henn, B. M. (2023). Increased homozygosity due to endogamy results in fitness consequences in a human population. Proceedings of the National Academy of Sciences, 120(43), e2309552120.CrossRefGoogle Scholar
Von Rueden, C., Gurven, M., & Kaplan, H. (2011). Why do men seek status? Fitness payoffs to dominance and prestige. Proceedings of the Royal Society B: Biological Sciences, 278(1715), 22232232.CrossRefGoogle ScholarPubMed
Willführ, K. P., Johow, J., & Voland, E. (2018). When the mother-in-law is just as good – Differential mortality of reproductive females by family network composition. PLoS One, 13(3), e0193252.CrossRefGoogle ScholarPubMed
Winking, J., Gurven, M., & Kaplan, H. (2011). Father death and adult success among the Tsimane: Implications for marriage and divorce. Evolution and Human Behavior, 32(2), 7989.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Research questions and inferential estimands

Figure 1

Table 2. Summary of results on frequency of cousin marriage and pool of potential partners. Total sample size of married adults for all these analyses n = 1310. The dashed lines separate analyses conducted with different statistical models. The parameters extracted for questions 2.2 and 2.3 are, respectively, the slope (β) and intercept (θ) coefficients related to the same predictor from the same model. The parameters extracted for question 3.2 pertain to the effects of two predictors in the same model. The estimates for questions 2.2, 2.3, 3.1 and 3.2 are from analyses with standardised predictors

Figure 2

Figure 1. Observed (circles) and predicted (shaded region) frequency of cousin marriage (y axis) for birth-years between 1964 and 1993 (x axis). The shaded region represents the 89% compatibility interval of the predicted distribution. The predicted mean estimate slightly declined from 29% in 1964 to 22% in 1993. The circles represent annual birth-year ‘bins’, within which the empirical frequency of cousin marriage is displayed for illustration purposes; all analyses are based on the individual level data, and the sample size of individuals combined in each bin is indicated by the scale of the circles (ranging from 13 to 88 individuals). The total sample here includes n = 1331 individuals.

Figure 3

Figure 2. Observed (circles) and predicted (shaded region) frequency of cousin marriage (y axes) as a function of number of cousins available as potential partners (x axis, panel a) and as a function of the frequency of cousins in the pool of potential partners of each individual (x axis, panel b). The shaded regions represent the 89% compatibility intervals of the predicted frequencies. All analyses are based on the individual level data. For illustration, each circle represents a ‘bin’, an x-axis range in which we calculate frequency estimates for illustration purposes. The size of the circles is scaled by sample size (from 2 individuals, for the highest x-axis ranges, to 329 individuals, for the lowest). The dashed line—which has a slope of 1—represents a linear relationship between the frequency of cousin marriage and the frequency of cousins available as partners; it represents the expected level of cousin marriage under random mating.

Figure 4

Figure 3. Observed (circles) and predicted (shaded region) frequency of cousins in the pool of potential partners of each individual (y axis) as a function of annual birth cohort (x axis). The shaded region represents the 89% compatibility interval. The predicted mean estimate of frequency of cousins declined from 1% in 1964 to 0.3% in 1993. Each circle represents a ‘bin’, an x-axis averaged frequency estimate for each birth year calculated for illustration purposes. The size of the circles is scaled by sample size (ranging from 13 to 88 individuals).

Figure 5

Figure 4. Fertility and offspring survival as a function of cousin marriage. In panel (a) we plot expected number of children (y axis) for each age of life (x axis). In panel (b) we plot the posterior probabilities of cumulative survival (y axis) for each age of life from 0 to 5 (x axis). We plot estimates for female individuals in the left frames and male individuals in the right frames. Estimates are colored by marriage type categories. The bands represent 89% compatibility intervals and reflect the precision of the estimated mean.

Figure 6

Figure 5. Temporal trends in the association between cousin marriage and fertility/survival. In panel (a), we plot offspring survival outcomes. In panel (b), we plot fertility outcomes. And, in panel (c), we plot age at first reproduction outcomes. Each horizontal bar represents the posterior distribution of the contrasts (i.e., differences) between families practicing cousin marriage and families not practicing it, in each decadal birth cohort (questions 5). Estimates for males are plotted with black circles, estimates for females are plotted with white circles. In the first rows, estimates are averaged across the whole sample (questions 4). The numbers represent the sample sizes of each decadal birth cohort. The first number is the sample size of individuals in families practicing cousin marriage, and the second number is the sample size of individuals in families not practicing cousin marriage.

Supplementary material: File

Dalzero et al. supplementary material

Dalzero et al. supplementary material
Download Dalzero et al. supplementary material(File)
File 116.9 KB