Genetic structure of the Saguenay, 1852-1911: evidence from migration and isonymy matrices

The Saguenay is a region in northeastern Québec populated in the second half of the 19th century through migration from other parts of Québec. The present-day population of nearly 300,000 is the result of both immigration and high rates of intrinsic growth. This population has been of interest to geneticists because of the high incidence of certain hereditary diseases, notably spastic ataxia, tyrosinemia, agenesis of the corpus callosum, vitamin D-dependent rickets, and myotonic dystrophy. Parent-offspring migration and isonymy matrices were used to estimate random kinship using the Malécot model for six 10-year time periods from 1852-1911. Comparisons between two estimates of kinship--one from parent-offspring migration matrices (phi) and the other from isonymy (R)--and geographic distance were made using both product-moment and Mantel correlation. Comparisons of within- and between-subdivision kinship were made using nonparametric and Mantel correlation. Within-subdivision kinship from the phi matrix was also compared with kinship estimated from marriage dispensations for endogamous marriages. The estimates of random kinship from the parent-offspring matrices showed a good fit with geography. However, isonymy did not correlate well with geographic distance; and phi and R showed no correlation until the last two time periods, and the diagonal of phi did not correlate with the marriage dispensations. Examination of scatterplots of phi vs. R suggests that nonrandom migration during the process of settlement formation is responsible for the lack of correlation. While movement across space seems to be highly dependent on distance, nonrandom selection of migrants means that between-subdivision estimates of kinship based on migration are not congruent with those obtained by other methods. On the whole, genetic differentiation seems to have been low due to the high levels of movement between subdivisions and immigration. The weak dependence of genetic structure on geographic distances in the present population is demonstrated by mapping the geographic distribution of cases of three recessively inherited diseases.     

Inbreeding patterns in La Cabrera, Spain: dispensations, multiple consanguinity analysis, and isonymy

Marital structure and inbreeding coefficients were analyzed in La Cabrera, an isolated mountain region in northwestern Spain. A total of 5,714 marriages were celebrated from 1880 to 1989 in the 37 parishes of the area. The total frequency of consanguineous marriages (up to the fourth degree) is 23.05%; multiple consanguineous marriages are remarkably common, reaching 5.43% of the total. The first cousin/second cousin ratio (referred to as kinship-type frequencies) is 0.43. The inbreeding values are the highest recorded in Spain and in Europe: alpha3 is 4.82 x 10(-3) for the whole period and alpha4 is 6.78 x 10(-3) for 1880--1919. The temporal trend of inbreeding shows high values (alpha3 > 4.5 x 10(-3)) for a particularly long period (1900--1959) and a rapid decline from 1960 onward. This historical inbreeding trend is clearly related to changes in population size. The frequencies of multiple consanguineous marriages and the analysis of isonymy show that the inbreeding structure is related to geographic and demographic factors. Comparing the results at two hierarchical levels (La Cabrera as a whole and the 37 parishes individually), we conclude that the inbreeding values are affected by internal geographic subdivision of the population (Wahlund effect). Social and cultural factors, such as avoidance of or preference for consanguineous marriages, are less important but depend on the kinship type involved.     

Genetic structure of the Ards Peninsula,Northern Ireland: evidence from civil registers of marriage 1840-1911

This paper uses marital migration data transcribed from the Civil Registers of Marriage 1840-1911 to estimate kinship from migration matrices and isonymy in the Ards Peninsula, Northern Ireland. The distribution of religious denominations (Presbyterian, Episcopalian, and Roman Catholic) varies systematically throughout the region, with up to 77% Roman Catholic in the south and 81% Presbyterian in the north. Portavogie, a fishing village on the east coast, is exclusively Protestant, with a population 93% Presbyterian. Comparison of migration and isonymy with geographical distance by multidimensional scaling and the MATFIT procedure show Portavogie to be an outlier, more distantly related to other areas than its geographical position would predict. We suggest that this discrepancy is due to settlement history and occupational and religious isolation. Mantel tests show that marital migration is significantly related to geographical distance (rMG = 0.4257), as is the distribution of religious denominations (rRG = 0.5548) through settlement history. Migration is dependent on religion (rMR = 0.3674), and isonymy is dependent on migration (rIM= 0.2531) but not on geography or religion. With Portavogie omitted from the analysis, the dependence of migration on geography and on religion increases (rMG = 0.5583, rMR = 0.5646), as does the correlation between religion and geography (rRG = 0.7213). The dependence of isonymy on migration increases (rIM = 0.5103), and significant correlations between isonymy and religion (rIR = 0.4135) and isonymy and geography (rIG = 0.4660) appear. We argue that a full explanation of population structure requires geographical distance, settlement history, and the influence of religion and occupation to be taken into account.     

Population structure in the Connecticut Valley: II. A comparison of multidimensional scaling solutions of migration matrices and isonymy

In a previous paper (Swedlund et al., 1984) we have described the population structure of the historical Connecticut River Valley of Massachusetts in terms of matrimonial migration matrices. Using procedures described by Morton (1973), Harpending and Jenkins (1974), Jorde (1980), and others the exchanges between subdivisions which make up the matrices are made column stochastic and analyzed to predict genetic kinship. Subsequently the kinship estimates within and between subdivisions can be interpreted as genetic covariance and compared to the actual geographic distances between the respective subdivisions using a principal components analysis. In the present paper we extend these results by applying nonmetric multidimensional scaling to the migration matrices, and to isonymy matrices based on the same communities. We demonstrate that the multidimensional scaling configurations of marital migration represent the actual geographic relationships between the communities quite effectively for this particular case study from historical Massachusetts. Moreover, we argue that while these migration data may provide good estimates of social and genetic exchange between the subdivisions, surname analysis may also be informative of processes not revealed in the migration matrices alone.     

Population relationships by isonymy in frontier Pennsylvania.

Data from the 1800 US census were used to study relationships by isonymy among 7 civil subdivisions of Bedford County, Pennsylvania. Two analyses of the data were conducted. In the first analysis heads of household served as the unit of analysis. In the second analysis the total number of individuals in each household was used to correct for varying family sizes. All measures of internal differentiation were approximately doubled when the complete population numbers were used. The head-of-household analysis produced FST and RST values of 0.0012 and 0.0007, respectively; the complete population analysis yielded 0.0021 for FST and 0.0015 for RST. Interpopulation a priori kinship estimations were similar using both methods. Conditional kinship estimations varied more, with almost all values negative, but the head-of-household estimates were less negative. Multidimensional scaling of isonymy values coincided fairly well with actual geographic relationships, but a Mantel test revealed no significant relationship between geographic distances and isonymy, and isolation by distance values indicated a low relationship between the 2 measures. The population of the county was heterogeneous, with low kinship between its constituent communities. This appears to be a result of kin-structured long-distance migration rather than of local processes. The head-of-household values are more comparable with other studies and more representative of population relationships; complete population values exaggerate heterogeneity because of random fluctuations in household size.     

Genetic structure of the Utah Mormons: isonymy analysis

Isonymy analysis is reported for a sample of 188,895 marriages extracted from the Utah Genealogical Database. Inbreeding rates estimated by isonymy are low, ranging from 0.005 for the earliest marriage cohort (1800-1809) to 0.0008 in the most recent cohort (1950-1959). The inbreeding values decrease considerably through time, but they are consistently higher than inbreeding values estimated from pedigrees. Several explanations are offered for this, including polyphyletism of surnames and the presence of Scandinavian patronyms in this population. Random isonymy between subdivisions is also compared with random kinship estimated from migration matrices. In terms of within-subdivision kinship, the two approaches yield similar results. However, the results are quite dissimilar for between-subdivision kinship. This reflects the recent and nonrandom settlement of Utah by different ethnic groups with different surname distributions. In later time periods, the correlations between the two types of kinship estimates increase, showing that migration patterns (which are strongly determined by geographic distance) exert an increasing influence on the distribution of surnames. Logistic regression is performed on a subset of marriages (n = 88,202), using isonymous vs. nonisonymous marriage as the dependent variable. The independent variables are year of marriage, geographic distance between husband's and wife's birthplaces, endogamous vs. exogamous marriage, and population sizes of husband's and wife's birthplaces. Year of marriage and geographic distance are shown to be significant independent predictors of isonymous marriage.     

Surnames in the Azores: analysis of the isonymy structure

Geographic isolation is a significant factor to consider when characterizing human populations. The knowledge of the genetic structure of isolated populations has been of great importance to disease-locus positioning and gene identification. To investigate the genetic structure of the Azorean population, we conducted a survey based on the frequencies of surnames listed in the 2001 telephone book. We calculated the following parameters: isonymy (I), the random component of inbreeding (F(ST)), genetic diversity according to Fisher (alpha), Karlin-McGregor's migration rate (v), and Nei's distance. For the 1,271 subscribers and 163 different surnames, Graciosa island presented the lowest value of abundance of surnames (alpha = 15.75), suggesting great genetic isolation compared to the other eight islands. Migration, calculated on the basis of the diversity of surnames within islands, ranged from 0.2747 (Corvo island) to 0.0026 (S?o Miguel island), indicating that people migrated preferentially toward the economically more developed islands. The value of the random component of inbreeding obtained for the whole population (F(ST) = 0.0039) indicates little genetic differentiation (Wright's F(ST) < 0.05). Moreover, isonymy similarity revealed using the UPGMA method shows three subclusters corresponding to the geographic distribution of the islands.     

Surnames in Argentina: a population study through isonymy

In order to study the isonymic structure of Argentina, the surname distributions of 22.6 million electors registered for the year 2001 were analyzed in the 24 districts (distritos) and 541 municipalities (municipios) of the country. The number of different surnames found was 414,441. Matrices of isonymic distances between districts were constructed and tested for correlation with the geographic distance between the capital towns of the districts. We found that, for the whole of Argentina, Euclidean distance was correlated with the log of geographic distance (r=0.480+/- 0.067). A dendrogram of the 24 regions was built from the matrix of Euclidean distances, using the UPGMA method. The clusters identified by the dendrogram are coincident with conterminous geographical regions of the country. Random inbreeding calculated from isonymy, F(ST), was highest in La Rioja, Corrientes, and Santiago del Estero. It was lowest in the area of Buenos Aires and in the north-central region of Santa Fé. Average Fisher's alpha for municipalities was 358; for districts, it was 422; and for Argentina as a unit, it was 602. The geographical distribution of alpha in 541 municipalities, high in the east and lower in the west of the country, is compatible with the settlement in the 20th century of subsequent waves of immigrants moving from the North Atlantic coast toward the foot of the Andes and toward the south. The present structure of Argentina indicates that migration dominates over drift.     

Interpopulation relationship by isonymy: application to ethnosocial groups and illegitimacy

The application of the isonymic method to establish interpopulation relationships is made difficult by such factors as: (1) a small population size; (2) the subdivision of the population into ethnosocial groups; and (3) the existence of individuals born extramaritally. The present study analyzes the validity of the isonymy method in populations where such difficulties exist. Lasker's R(ij) relationship coefficients were calculated in base to marital records from six Pocho parishes (Argentina) for the period 1766 to 1840. Three endogamous ethnosocial groups were considered-Spanish, American Indians, "Mestizos"--and a fourth group combining mates of the three previous groups. For each ethnosocial category, and taking into account paternal surnames, maternal surnames (both for legitimate and illegitimate mates), and whole surnames merged, R(i) interparish relationship matrices were obtained. All these matrices were correlated by means of the Mantel test. Maternal surnames of illegitimate mates show a similar pattern to the maternal surnames of legitimate mates and to all surnames of each category. Groups with larger sample size in every parish intercorrelate regardless of ethnosocial category. Results suggest the convenience of maximizing the sample size and using maternal surnames in populations with high illegitimacy.     

Microsatellite Polymorphisms in Cassava Landraces from the Cerrado Biome, Mato Grosso do Sul, Brazil

Using nine microsatellite loci, we investigated genetic structure and diversity in 83 Brazilian cassava accessions, including several landraces, in the Cerrado biome in Mato Grosso do Sul, Brazil. All nine loci were polymorphic, averaging 6.00 alleles per locus. Treating each of seven municipalities as a cassava group or population, they averaged 3.5 alleles per locus, with 97% polymorphic loci, high values for observed heterozygosity (0.32) and gene diversity (0.56). Total genetic variability was high (0.668), and most of this genetic variability was concentrated within municipalities (0.577). Cluster and structure analyses divided accessions into two major clusters or populations (K = 2). Also, a significant genetic versus geographic correlation was found (r = 0.4567; P < 0.0260). Migratory routes in the Cerrado are considered main contributors to the region’s high cassava diversity and spatial genetic structure, amplifying interactions between traditional farmers and the evolutionary dynamics of this crop.     

Geographic patterns of (genetic, morphologic, linguistic) variation: how barriers can be detected by using Monmonier's algorithm

When sampling locations are known, the association between genetic and geographic distances can be tested by spatial autocorrelation or regression methods. These tests give some clues to the possible shape of the genetic landscape. Nevertheless, correlation analyses fail when attempting to identify where genetic barriers exist, namely, the areas where a given variable shows an abrupt rate of change. To this end, a computational geometry approach is more suitable because it provides the locations and the directions of barriers and because it can show where geographic patterns of two or more variables are similar. In this frame we have implemented Monmonier's (1973) maximum difference algorithm in a new software package to identify genetic barriers. To provide a more realistic representation of the barriers in a genetic landscape, we implemented in the software a significance test by means of bootstrap matrices analysis. As a result, the noise associated with genetic markers can be visualized on a geographic map and the areas where genetic barriers are more robust can be identified. Moreover, this multiple matrices approach can visualize the patterns of variation associated with different markers in the same overall picture. This improved Monmonier's method is highly reliable and can be applied to nongenetic data whenever sampling locations and a distance matrix between corresponding data are available.     

Consanguinity in two Spanish regions: La Cabrera and Fuentes Carrionas. Dispensations and isonymy

In this work the level and structure of consanguinity is analysed in two Spanish rural regions of similar geographic and orographic characteristics for the period between 1880 and 1979, employing two different methodologies. The estimates according dispensations shows that the total levels (alpha4 = 0.00552 in La Cabrera and 0.00405 in Fuentes Carrionas) and the structures of consanguinity (C22/C33 = 0.43 and 0.34) are similar in both regions, but have evolved differently. Whereas in La Cabrera both parameters have remained stable, in Fuentes Carrionas they fluctuated through the period analysed. On studying the structure of consanguinity more closely using the isonymic method, it can be seen that in La Cabrera total consanguinity (Ft = 0.0206) should be attributed mainly to environmental factors (Fr = 0.0193) and to a lesser extent to socio-cultural factors (Fn = 0.0013), whereas in Fuentes Carrionas it derives, almost exclusively, from the former (Ft = 0.01270; Fr = 0.01589; Fn = -0.00325).     

Surnames in Texas: a population study through isonymy

To study the isonymy structure of Texas, we analyzed the surname distributions of 3.6 million telephone users registered for the year 1996 in 232 towns distributed in the 7 regions of the state. The number of different surnames was 235,740. Matrices of isonymy distances between towns and between geographic regions were constructed and tested for correlation with geographic distance. We found that isonymy distances between the seven regions showed borderline or no correlation with geographic distance, with r = 0.089 +/- 0.232, r = 0.492 +/- 0.232, and r = 0.337 +/- 0.232 for Lasker's, Euclidean, and Nei's distances, respectively. Isonymy distances between towns were significantly correlated with geographic distance, with r = 0.249 +/- 0.006 for Lasker's distance, r = 0.338 +/- 0.006 for the Euclidean distance, and r = 0.418 +/- 0.006 for Nei's distance. Two dendrograms, one for the 7 regions and one for the 232 towns, were built from the matrices of Nei's distances. The dendrogram for regions indicates that a main surname differentiation exists between the East and West areas of Texas, with West Texas being predominantly Hispanic and East Texas being predominantly English-speaking. The dendrogram for the towns confirms in detail the differences identified by the matrix of distances between regions. Random inbreeding calculated from isonymy, F(ST), was highest in the west and in the south of the state. It was lowest in the area of Austin and Houston. Average Fisher's alpha for towns was 734, for regions it was 1,047, and for Texas as a whole it was 1,230. The geographic distribution of alpha in the state shows distinctly lower values in the traditionally Hispanic west and higher values in the east and on the Gulf of Mexico.     

A statistical procedure to assess the significance level of barriers to gene flow

Although several methods are available to study the extent of isolation by distance （IBD） among natural populations, comparatively few exist to detect the presence of sharp genetic breaks in genetic distance datasets. In recent years, Monmonier＇s maximum-difference algorithm has been increasingly used by population geneticists. However, this method does not provide means to measure the statistical significance of such barriers, nor to determine their relative contribution to population differentiation with respect to IBD. Here, we propose an approach to assess the significance of genetic boundaries. The method is based on the calculation of a multiple regression from distance matrices, where binary matrices represent putative genetic barriers to test, in addition to geographic and genetic distances. Simulation results suggest that this method reliably detects the presence of genetic barriers, even in situations where IBD is also significant. We also illustrate the methodology by analyzing previously published datasets. Conclusions about the importance of genetic barriers can be misleading if one does not take into consideration their relative contribution to the overall genetic structure of species.     

Population structure of São Miguel Island, Azores: a surname study

The knowledge of a population structure may constitute a powerful tool for mapping genes underlying susceptibility to Mendelian and complex diseases. To obtain a better understanding of the population structure of S?o Miguel Island (Azorean Archipelago, Portugal), we carried out a surname survey using the surnames listed in the most recent telephone book (2001). We identified 1315 different surnames in a total of 27,621 subscribers. The frequency of the different surnames was used to calculate the following parameters: isonymy (I), random component of inbreeding (FST), genetic diversity according to Fisher (alpha), migration rate according to Karlin-McGregor (v) and Nei's genetic distance. Eleven localities were selected, according to population size and geographic distribution, for analysis using the above parameters. Our results show that 51% of Salga's population and 52% of Sete Cidades's population are represented by six and eight surnames, respectively. These figures demonstrate the effective isolation of these two small places, which are located at opposite ends of S?o Miguel Island. Salga, Achada, and Sete Cidades present the lowest values of Fisher's alpha, indicating less genetic diversity. In contrast, the capital, Ponta Delgada, presents the highest value of alpha (78.13), indicating more genetic diversity. Our data indicate that the clustering of the localities corresponds to the geographic features of the island, where localities close together tend to share similar surnames. In conclusion, the population of S?o Miguel is relatively homogeneous and may constitute an ideal model for genetic mapping studies.     

Ethnicity and biodemographic structure in the Arbëreshe of the province of Cosenza, southern Italy, in the XIX century

Cultural and environmental factors interact in determining the genetic structure of human populations. Bio-demographic investigations of ethnic minorities are able to disentangle the influences that these two components have on the evolution of the genetic structure of a population. The ethnic minority of the Arb?reshe of the province of Cosenza (Calabria, southern Italy) is analyzed in this paper and its bio-demographic structure in the early 1800s is compared with that of neighboring Italian populations. The data derive from surnames recorded in the birth registers of the 19 Arbdreshe municipalities of the province of Cosenza and in 5 non-Arb?reshe municipalities of the same province. Isonymy and repeated pairs of surnames are used to analyze the bio-demographic structure of these populations, while analysis of isonymic relationships is used to investigate the variability between populations. Higher values of marital isonymy and subdivision into subpopulations characterize the Arb?reshe populations with respect to their non-Arb?reshe neighbors. However, the high range of variability of these parameters suggests a strong influence of geographic location on the marriage pattern of each community. At the same time, cultural differences linked to group identity had a strong impact in limiting marriage exchanges between the different ethnic groups living in the province of Cosenza in the early 1800s. In fact, the analysis of isonymic relationships demonstrates that geographic location shaped kinship patterns among the Arbereshe communities, but it also shows that the non-Arb?reshe neighbors formed a clearly separate reproductive cluster.     

Environmental factors affecting population level genetic divergence of the striped field mouse (<Emphasis Type="Italic">Apodemus agrarius</Emphasis>) in South Korea

The cognizing of connectivity among small mammal populations across heterogeneous landscapes is complicated due to complex influences of landscape and anthropogenic factors on gene flow. A landscape genetics approach offers inferences on how landscape features drive population structure. Through a landscape genetics approach, we investigated influences of geographical, environmental, and anthropogenic features on populations of Apodemus agrarius, the striped field mouse, the prime vector of hemorrhagic fever by a landscape genetic approach. We identified landscape features that might affect the population structure of striped field mice by analyzing microsatellite markers of 197 striped field mice from 21 populations throughout South Korea. We developed Maximum-likelihood population effects models based on landscape distances and resistance matrices and pairwise F_ST values for meta-populations of striped field mouse. We also conducted Mantel and partial Mantel tests to investigate geographic patterns of genetic similarities. In Mantel and partial Mantel tests, the F_ST was significantly correlated with all three models of movement; movement cost, Euclidian distance and least-cost distance, although the magnitudes of correlations varied. The 4 top-ranked models included three variables; temperature, precipitation and one human disturbance factor (population). We did not attain a significant effect for anthropogenic factors on genetic similarities among populations in the Korean striped field mouse, but we confirmed a significant association for genetic similarity with climatic features (temperature and precipitation).     

A geographic mosaic of passive dispersal: population structure in the endemic Hawaiian amber snail Succinea caduca (Mighels, 1845)

We used 276 cytochrome c oxidase subunit I (COI, 645 bp) and a subset of 84 16S large ribosomal subunit (16S, 451 bp) sequences to evaluate geographic patterns of genetic variation in 24 populations of the endemic Hawaiian land snail Succinea caduca spanning its range on six islands. Haplotype networks, gene tree topologies, pairwise molecular divergence and F _ST matrices suggest substantial geographic genetic structuring and complex dispersal patterns. Low nucleotide diversity and low pairwise molecular divergence values within populations coupled with higher between population values suggest multiple founder events. High overall haplotype diversity suggests diversification involving rare interpopulation dispersal, fragmentation by historical lava flows and variation in habitat structure. Within-island rather than between-island population comparisons accounted for the majority of molecular variance. Although 98% of 153 COI haplotypes were private by population, a Mantel test showed no evidence for isolation by distance. Mismatch distributions and population partitioning patterns suggest that genetic fragmentation has been driven by punctuated, passive dispersal of groups of closely related haplotypes that subsequently expanded and persisted in isolation for long periods (average > 2 million years ago), and that Pleistocene island connections may have been important in enhancing gene flow. Historical availability of mesic coastal habitat, together with effective dispersal may explain the long-term persistence and unusual multi-island distribution of this species, contrasting with the single-island endemism of much of the Hawaiian biota.     

The names of Spain: a study of the isonymy structure of Spain

In order to estimate the isonymy structure of Spain, we studied surname distribution in 283 Spanish towns based on 3.625 million telephone users selected from 6.328 million users, downloaded from a commercial CD-ROM which contains all 13 million users in the country. Since in Spain the surname is made by the paternal and the maternal surname, it was possible to classify surnames according to parental origin. Two matrices of isonymy distances, one for paternal and one for maternal surnames, were constructed and tested for correlation with geographic distance. For the whole of Spain, Euclidean distance was significantly but weakly correlated with geographic distance both for paternal and maternal surnames, with r = 0.205 +/- 0.013 and r = 0.263 +/- 0.012, respectively. Two dendrograms of the 283 sampled towns were built from the two matrices of Euclidean distance. They are largely colinear. Four main clusters identified by the dendrograms are correlated with geography. Given the surname structure of Spain, we were able to calculate from isonymy and for each town 1). total or expressed inbreeding, 2). random or expected inbreeding, and 3). local inbreeding. Total inbreeding, F(IT), was highest in the North Atlantic regions and lowest along the Mediterranean Coast. The lowest levels were found in Andalusia, Catalunyia, Valencia, and Navarra. Random inbreeding, F(ST), had a similar geographical pattern. Local inbreeding, F(IS), was relatively uniform in the whole of Spain. In towns, random inbreeding dominates over local inbreeding. From the analysis, it emerges that the northwestern area of Spain is the most inbred.     

Biodemographic and Genetic Structure of Zamora Province (Spain): Insights from Surname Analysis

Located in the north-western extreme of the Iberian northern plateau with the Portugal border on the west, Zamora province is one of the current divisions of the autonomous region of Castile-Leon (Spain). According to natural boundaries and historical records, the province of Zamora can be subdivided into six different regions: Aliste, Benavente, Bajo Duero, Campos-Pan, Sanabria and Sayago. The geography of the area is configured by two major rivers: the Duero that crosses the province in a west-east direction and its affluent the Esla that crosses the province in a northwest-southeast course. In this work, we analyze both surnames of 166,349 individuals from Zamora province. The main goal of this study is to explore the differential weight of historical, demographic, geographic and sociocultural factors in shaping the bio-demographic and genetic structure of this specific population. The highest value of total consanguinity based on random isonymy is observed in Sanabria, which also has the lowest values of heteronymy and the highest intra-population a priori kinship. By contrast, Sayago exhibits the highest value of heteronymy and the lowest intra-population a priori kinship. The consanguinity values observed in Zamora are, in part, related to the population structure. Notwithstanding, the values of non-random consanguinity observed in some municipalities could be associated with particular mating behaviors that favor isonymic marriages. Finally, the results of analyses indicate that the levels of micro-differentiation in the province cannot be considered high; however, some degree of substructure could be deduced. It appears that the current population structure of Zamora has been determined not only by settlement history, geographic distance and geographic boundaries, but also, by land system distribution and population size.