The magnitude and origin of European-American admixture in the Gila River Indian Community of Arizona: a union of genetics and demography.

Complementary genetic and demographic analyses estimate the total proportion of European-American admixture in the Gila River Indian Community and trace its mode of entry. Among the 9,616 residents in the sample, 2,015 persons claim only partial Native American heritage. A procedure employing 23 alleles or haplotypes at eight loci was used to estimate the proportion of European-American admixture, m(a), for the entire sample and within six categories of Caucasian admixture calculated from demographic data, md. The genetic analysis gave an estimate of total European-American admixture in the community of 0.054 (95% confidence interval [CI] .044-.063), while an estimate from demographic records was similar, .059. Regression of m(a) on md yielded a fitted line m(a) = .922md, r = .959 (P = .0001). When total European-American admixture is partitioned between the contributing populations, Mexican-Americans have provided .671, European-Americans .305, and African-Americans .023. These results are discussed within the context of the ethnic composition of the Gila River Indian Community, the assumptions underlying the methods, and the potential that demographic data have for enriching genetic measurements of human admixture. It is concluded that, despite the severe assumptions of the mathematical methods, accurate, reliable estimates of genetic admixture are possible from allele and haplotype frequencies, even when there is little demographic information for the population.     

On the formula for admixture linkage disequilibrium

Admixture linkage disequilibrium (ALD), a phenomenon created by gene flow between genetically distinct populations, has for some time been used as a tool in gene mapping. It is therefore important to analyze the pattern of ALD over generations. In this study we explore two models of admixture: the gradual admixture (GA) model, in which admixture occurs at a variable rate in every generation; and the immediate admixture (IA) model, a special case of the GA model, in which admixture occurs in a single generation. In the case of ALD, the well-known formula of linkage disequilibrium (Delta(t)=(1-r)t Delta(0)) is not applicable under these two models. We note the effect of a random mating population (RMP) to the gametic frequencies from the parental population to the offspring population, and provide the correct formula for ALD.     

Source populations of Quercus glauca in the last glacial age in Taiwan revealed by nuclear microsatellite markers

In this work, we attempted to study genetic differentiation between populations of Quercus glauca in Taiwan using nuclear microsatellite markers to infer the potential refugium in the last glaciation stage. Four microsatellite loci for 20 individuals each in 10 populations of Taiwan were analyzed. We found that Q. glauca has relatively high within-population diversity (H(E) = 0.741) and low population differentiation (F(ST) = 0.042) but shows isolation by distance. The most divergent populations, according to the average F(ST) for individual populations in comparison with every other population, were found in populations Cy, Sa, and Hy in southern Taiwan and Pa in north-central Taiwan. Moreover, populations Cy, Sa, and Pa were recognized as being the source populations for gene recolonization after the last glaciation stage. In addition, the three sites of Wu, Ym, and Cy exhibited the highest gene diversities that coincided with populations with the highest chloroplast DNA variations. This may have resulted from an admixture of colonization routes. In conclusion, observations of the most divergent populations and source populations suggest that southern and probably north-central Taiwan may have potentially been refugia for Q. glauca in the last glaciation. This agrees with the possible refugium in southern Taiwan revealed by a previous study using chloroplast DNA markers.     

Admixture analysis of a rural population of the state of Guerrero, Mexico

We studied 156 individuals of Native American descent from the city of Tlapa in the state of Guerrero in western Mexico. Most individuals' ethnicity was either Nahua, Mixtec, or Tlapanec, but self-identified Mestizos and individuals of mixed ethnicities were also included in the sample. We typed 24 autosomal, one Y-chromosome, and four mitochondrial ancestry-informative markers (AIMs) to estimate group and individual admixture proportions, and determine whether the admixture process involved directional gene flow between parental groups. When genetically defined (GD) Mestizos were excluded from the analysis, Native American ancestry represented approximately 98% of the population's gene pool, while European and West African ancestry represented approximately 1% each. Maternally inherited markers also showed an exceptionally high Native American contribution (98.5%), as did the paternally inherited marker, DYS199 (90.7%). We did not detect genetic structure in this population using these AIMs, which appears consistent with the homogeneity of the sample in terms of admixture proportions. The addition of GD Mestizos to the sample did not produce a considerable change in admixture estimates, but it had a major effect on population structure. These results show that the population of Tlapa in Guerrero, Mexico, has experienced little admixture with Europeans and/or West Africans. They also show that the impact of a small number of admixed individuals on an otherwise homogeneous population might have profound implications on subsequent ancestry/phenotype analysis and mapping strategies. We suggest that heterogeneity is a major characteristic of Mexican populations and, as a consequence, should not be disregarded when designing epidemiological studies of Mexican and Mexican American populations.     

Role of the altitude level on cerebral autoregulation in residents at high altitude.

Cerebral autoregulation is impaired in Himalayan high-altitude residents who live above 4,200 m. This study was undertaken to determine the altitude at which this impairment of autoregulation occurs. A second aim of the study was to test the hypothesis that administration of oxygen can reverse this impairment in autoregulation at high altitudes. In four groups of 10 Himalayan high-altitude dwellers residing at 1,330, 2,650, 3,440, and 4,243 m, arterial oxygen saturation (Sa(O(2))), blood pressure, and middle cerebral artery blood velocity were monitored during infusion of phenylephrine to determine static cerebral autoregulation. On the basis of these measurements, the cerebral autoregulation index (AI) was calculated. Normally, AI is between zero and 1. AI of 0 implies absent autoregulation, and AI of 1 implies intact autoregulation. At 1,330 m (Sa(O(2)) = 97%), 2,650 m (Sa(O(2)) = 96%), and 3,440 m (Sa(O(2)) = 93%), AI values (mean +/- SD) were, respectively, 0.63 +/- 0.27, 0.57 +/- 0.22, and 0.57 +/- 0.15. At 4,243 m (Sa(O(2)) = 88%), AI was 0.22 +/- 0.18 (P < 0.0005, compared with AI at the lower altitudes) and increased to 0.49 +/- 0.23 (P = 0.008, paired t-test) when oxygen was administered (Sa(O(2)) = 98%). In conclusion, high-altitude residents living at 4,243 m have almost total loss of cerebral autoregulation, which improved during oxygen administration. Those people living at 3,440 m and lower have still functioning cerebral autoregulation. This study showed that the altitude region between 3,440 and 4,243 m, marked by Sa(O(2)) in the high-altitude dwellers of 93% and 88%, is a transitional zone, above which cerebral autoregulation becomes critically impaired.     

Validity of two-component models for estimating body fat of black men.

Commonly used two-component model conversion formulas that estimate relative body fat (%BF) from body density (Db) were cross-validated on a heterogeneous sample of black men (n = 30; age = 19--45 yr). A four-component model was used to obtain criterion measures of %BF, and linear regression and analysis of individual residual scores were conducted to assess the predictive accuracy of the formulas under investigation. The two-component formula commonly used to estimate %BF of black men (Schutte JE, Townsend EJ, Hugg J, Shoup RF, Malina RM, and Blomqvist CG. J Appl Physiol 56: 1647-1649, 1984) significantly (P < or = 0.01) and systematically (87% of sample) overestimated %BF (-1.28%); thus we developed the following two-component Db conversion formula: %BF = [(4.858/Db) - 4.394] x 100. Because our formula was derived from a four-component model and a larger, more heterogeneous sample than the commonly used two-component formula, we recommend using it to convert Db to %BF for black men. Additionally, there was good agreement between dual-energy X-ray absorptiometry and the four-component model, making this a suitable alternative for estimating the %BF of black men.     

A classical likelihood based approach for admixture mapping using EM algorithm

Several disease-mapping methods have been proposed recently, which use the information generated by recent admixture of populations from historically distinct geographic origins. These methods include both classic likelihood and Bayesian approaches. In this study we directly maximize the likelihood function from the hidden Markov Model for admixture mapping using the EM algorithm, allowing for uncertainty in model parameters, such as the allele frequencies in the parental populations. We determined the robustness of the proposed method by examining the ancestral allele frequency estimate and individual marker-location specific ancestry when the data were generated by different population admixture models and no learning sample was used. The proposed method outperforms a widely used Bayesian MCMC strategy for data generated from various population admixture models. The multipoint information content for ancestry was derived based on the map provided by Smith et al. (2004) and the associated statistical power was calculated. We examined the distribution of admixture LD across the genome for both real and simulated data and established a threshold for genome wide significance applicable to admixture mapping studies. The software ADMIXPROGRAM for performing admixture mapping is available from authors.     

Estimating genetic correlations in natural populations in the absence of pedigree information: accuracy and precision of the Lynch method

Usually, genetic correlations are estimated from breeding designs in the laboratory or greenhouse. However, estimates of the genetic correlation for natural populations are lacking, mostly because pedigrees of wild individuals are rarely known. Recently Lynch (1999) proposed a formula to estimate the genetic correlation in the absence of data on pedigree. This method has been shown to be particularly accurate provided a large sample size and a minimum (20%) proportion of relatives. Lynch (1999) proposed the use of the bootstrap to estimate standard errors associated with genetic correlations, but did not test the reliability of such a method. We tested the bootstrap and showed the jackknife can provide valid estimates of the genetic correlation calculated with the Lynch formula. The occurrence of undefined estimates, combined with the high number of replicates involved in the bootstrap, means there is a high probability of obtaining a biased upward, incomplete bootstrap, even when there is a high fraction of related pairs in a sample. It is easier to obtain complete jackknife estimates for which all the pseudovalues have been defined. We therefore recommend the use of the jackknife to estimate the genetic correlation with the Lynch formula. Provided data can be collected for more than two individuals at each location, we propose a group sampling method that produces low standard errors associated with the jackknife, even when there is a low fraction of relatives in a sample.     

Admixture estimates for Caracas, Venezuela, based on autosomal, Y-chromosome, and mtDNA markers

The present Venezuelan population is the product of admixture of Amerindians, Europeans, and Africans, a process that was not homogeneous throughout the country. Blood groups, short tandem repeats (STRs), mtDNA, and Y-chromosome markers have been used successfully in admixture studies, but few such studies have been conducted in Venezuela. In this study we aim to estimate the admixture components of samples from two different socioeconomic levels from Caracas, Venezuela's capital city, compare their differences, and infer sexual asymmetry in the European Amerindian union patterns. Gene frequencies for blood groups ABO and Rh (CDE) and for the STRs VWA, F13A01, and FES/FPS and mtDNA and Y-chromosome haplogroups were studied in a sample of 60 individuals living in Caracas, taken from a private clinic (high socioeconomic level), and 50 individuals, also living in Caracas, drawn from a public maternity clinic (low socioeconomic level). The admixture analysis for the five autosomal markers gives a high European component (0.78) and an almost negligible African sub-Saharan component (0.06) for the high socioeconomic level, whereas for the low socioeconomic level the sub-Saharan, European, and Amerindian components were 0.21, 0.42, and 0.36, respectively. Estimates of admixture based on mtDNA and Y-chromosome markers reveal that the Amerindian contribution to these Caracas samples is almost entirely through females, because the Y-chromosome Amerindian and African sub-Saharan chromosomes found in this study were scarce. Our study reveals that the identification of the grandparents' geographic origin is an important methodological aspect to take into account in genetic studies related to the reconstruction of historical events.     

Effect of repeated normobaric hypoxia exposures during sleep on acute mountain sickness, exercise performance, and sleep during exposure to terrestrial altitude

There is an expectation that repeated daily exposures to normobaric hypoxia (NH) will induce ventilatory acclimatization and lessen acute mountain sickness (AMS) and the exercise performance decrement during subsequent hypobaric hypoxia (HH) exposure. However, this notion has not been tested objectively. Healthy, unacclimatized sea-level (SL) residents slept for 7.5 h each night for 7 consecutive nights in hypoxia rooms under NH [n = 14, 24 ± 5 (SD) yr] or "sham" (n = 9, 25 ± 6 yr) conditions. The ambient percent O(2) for the NH group was progressively reduced by 0.3% [150 m equivalent (equiv)] each night from 16.2% (2,200 m equiv) on night 1 to 14.4% (3,100 m equiv) on night 7, while that for the ventilatory- and exercise-matched sham group remained at 20.9%. Beginning at 25 h after sham or NH treatment, all subjects ascended and lived for 5 days at HH (4,300 m). End-tidal Pco(2), O(2) saturation (Sa(O(2))), AMS, and heart rate were measured repeatedly during daytime rest, sleep, or exercise (11.3-km treadmill time trial). From pre- to posttreatment at SL, resting end-tidal Pco(2) decreased (P < 0.01) for the NH (from 39 ± 3 to 35 ± 3 mmHg), but not for the sham (from 39 ± 2 to 38 ± 3 mmHg), group. Throughout HH, only sleep Sa(O(2)) was higher (80 ± 1 vs. 76 ± 1%, P < 0.05) and only AMS upon awakening was lower (0.34 ± 0.12 vs. 0.83 ± 0.14, P < 0.02) in the NH than the sham group; no other between-group rest, sleep, or exercise differences were observed at HH. These results indicate that the ventilatory acclimatization induced by NH sleep was primarily expressed during HH sleep. Under HH conditions, the higher sleep Sa(O(2)) may have contributed to a lessening of AMS upon awakening but had no impact on AMS or exercise performance for the remainder of each day.     

Estimating the timing of multiple admixture events using 3-locus linkage disequilibrium

Estimating admixture histories is crucial for understanding the genetic diversity we see in present-day populations. Allele frequency or phylogeny-based methods are excellent for inferring the existence of admixture or its proportions. However, to estimate admixture times, spatial information from admixed chromosomes of local ancestry or the decay of admixture linkage disequilibrium (ALD) is used. One popular method, implemented in the programs ALDER and ROLLOFF, uses two-locus ALD to infer the time of a single admixture event, but is only able to estimate the time of the most recent admixture event based on this summary statistic. To address this limitation, we derive analytical expressions for the expected ALD in a three-locus system and provide a new statistical method based on these results that is able to resolve more complicated admixture histories. Using simulations, we evaluate the performance of this method on a range of different admixture histories. As an example, we apply the method to the Colombian and Mexican samples from the 1000 Genomes project. The implementation of our method is available at https://github.com/Genomics-HSE/LaNeta.     

Ancestral proportions and admixture dynamics in geographically defined African Americans living in South Carolina

We analyzed admixture in samples of six different African-American populations from South Carolina: Gullah-speaking Sea Islanders in coastal South Carolina, residents of four counties in the "Low Country" (Berkeley, Charleston, Colleton, and Dorchester), and persons living in the city of Columbia, located in central South Carolina. We used a battery of highly informative autosomal, mtDNA, and Y-chromosome markers. Two of the autosomal markers (FY and AT3) are linked and lie 22 cM apart on chromosome 1. The results of this study indicate, in accordance with previous historical, cultural, and anthropological evidence, a very low level of European admixture in the Gullah Sea Islanders (m = 3.5 +/- 0.9%). The proportion of European admixture is higher in the Low Country (m ranging between 9. 9 +/- 1.8% and 14.0 +/- 1.9%), and is highest in Columbia (m = 17.7 +/- 3.1%). A sex-biased European gene flow and a small Native American contribution to the African-American gene pool are also evident in these data. We studied the pattern of pairwise allelic associations between the FY locus and the nine other autosomal markers in our samples. In the combined sample from the Low Country (N = 548), a high level of linkage disequilibrium was observed between the linked markers, FY and AT3. Additionally, significant associations were also detected between FY and 4 of the 8 unlinked markers, suggesting the existence of significant genetic structure in this population. A continuous gene flow model of admixture could explain the observed pattern of genetic structure. A test conditioning on the overall admixture of each individual showed association of ancestry between the two linked markers (FY and AT3), but not between any of the unlinked markers, as theory predicts. Thus, even in the presence of genetic structure due to continuous gene flow or some other factor, it is possible to differentiate associations due to linkage from spurious associations due to genetic structure.     

Structural and bioindicative aspects of fluctuated asymmetry of bilateral organisms

We have proposed and validated a method for quantitative assessment of phenotypic diversity of natural populations. Method is based on the fluctuated asymmetry (FA) indices of bilateral organisms, and it is applicable for biondicative investigations. Convolution of functions was proposed to estimate the mean (population) value of FA complex of features. This function could be written as finit sum [formula: see text] where eta is power of sample invariance (symmetry) for m individuals (i = 1, m). Eta is characterized by n asymmetric characteristics (j = 1, n) for the left (L) and the right (R) sides of the body. We have validated applicability of generalized function of desirability [formula: see text] (where di is partial desirability function [0,1]) for cumulative characterization of environment quality with results of bioindicative investigations. The value of function coincides with the value of symmetry of indicating species in this case.     

African ancestry of the population of Buenos Aires

The population of Argentina today does not have a “visible” black African component. However, censuses conducted during most of the 19th century registered up to 30% of individuals of African origin living in Buenos Aires city. What has happened to this African influence? Have all individuals of African origin died, as lay people believe? Or is it possible that admixture with the European immigrants made the African influence “invisible?” We investigated the African contribution to the genetic pool of the population of Buenos Aires, Argentina, typing 12 unlinked autosomal DNA markers in a sample of 90 individuals. The results of this analysis suggest that 2.2% (SEM = 0.9%) of the genetic ancestry of the Buenos Aires population is derived from Africa. Our analysis of individual admixture shows that those alleles that have a high frequency in populations of African origin tend to concentrate among 8 individuals in our sample. Therefore, although the admixture estimate is relatively low, the actual proportion of individuals with at least some African influence is approximately 10%. The evidence we are presenting of African ancestry is consistent with the known historical events that led to the drastic reduction of the Afro‐Argentine population during the second half of the 19th century. However, as our results suggest, this reduction did not mean a total disappearance of African genes from the genetic pool of the Buenos Aires population. Am J Phys Anthropol 128:164‐170, 2005. © 2005 Wiley‐Liss, Inc.     

Caucasian genes in American blacks: new data.

Data on 15 polymorphic protein-coding loci are used to estimate the proportion of Caucasian genes in U.S. blacks from the greater-metropolitan area of Pittsburgh. Allele frequencies from U.S. whites of the same region and from a sample of Nigerians are used as representatives of the genetic contributions of the source populations between which admixture has occurred. These materials provide 18 unique variants that occur exclusively in the blacks and 5 variants that are restricted to the Caucasians only. As a result, when all segregating alleles (52) at these 15 loci are considered, the proportion (mean +/- SE) of Caucasian genes in U.S. blacks (25.2% +/- 2.7%) is estimated with a precision much better than that of all other previous estimates. The estimate based on the frequencies of these 18 unique variants of African origin (24.8% +/- 6.2%) is also consistent with the pooled estimate obtained from the 15 loci by the weighted least-square method. The homogeneity of locus-specific estimates of admixture indicates that these loci are appropriate for studying the proportion of black genes in any admixed population involving African admixture. The advantages of employing such loci for genetic-epidemiologic studies in U.S. blacks is discussed in the context of the availability of these large number of unique African variants at these protein loci.     

Mitochondrial DNA (mtDNA) haplotypes reveal maternal population genetic affinities of Sea Island Gullah-speaking African Americans

To better understand the population substructure of African Americans living in coastal South Carolina, we used restriction site polymorphisms and an insertion/deletion in mitochondrial DNA (mtDNA) to construct seven-position haplotypes across 1,395 individuals from Sierra Leone, Africa, from U.S. European Americans, and from the New World African-derived populations of Jamaica, Gullah-speaking African Americans of the South Carolina Sea Islands (Gullahs), African Americans living in Charleston, South Carolina, and West Coast African Americans. Analyses showed a high degree of similarity within the New World African-derived populations, where haplotype frequencies and diversities were similar. Phi-statistics indicated that very little genetic differentiation has occurred within New World African-derived populations, but that there has been significant differentiation of these populations from Sierra Leoneans. Genetic distance estimates indicated a close relationship of Gullahs and Jamaicans with Sierra Leoneans, while African Americans living in Charleston and the West Coast were progressively more distantly related to the Sierra Leoneans. We observed low maternal European American admixture in the Jamaican and Gullah samples (m = 0.020 and 0.064, respectively) that increased sharply in a clinal pattern from Charleston African Americans to West Coast African Americans (m = 0.099 and 0.205, respectively). The appreciably reduced maternal European American admixture noted in the Gullah indicates that the Gullah may be uniquely situated to allow genetic epidemiology studies of complex diseases in African Americans with low European American admixture.     

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.     

Dissection of a continuous distribution: red cell galactokinase activity in blacks.

A significant difference between blacks and whites in the distribution of red cell galactokinase (GALK) has been found by Tedesco et al. [2]. From the shapes of the distributions, it was inferred that whites are essentially all homozygous for one allele (GALKA), but blacks are polymorphic. A second allele (GALKP), for lower GALK activity, is presented at high frequency in blacks but rare or absent in whites. This paper presents a method which, assuming the genetic model presented, estimates the genotype composition of the black sample. We make some reasonable biochemical assumptions and fit a mixture of three normal distributions to the black data to obtain an estimate of p, the frequency of GALKA in blacks. The fit of the model to the data is excellent and the best estimate of p is .217 +/- .025. Since admixture of white genes in blacks from the United States is known to be about 20%, the value of p implies that virtually all GALKA alleles were introduced by admixture, and that the ancestral black population was monomorphic for GALKP. If whites are indeed monomorphic for GALKA, they differ from unmixed blacks by a full gene substitution at the locus for GALK.     

Monovalent and divalent salt correction algorithms for Tm prediction--recommendations for Primer3 usage

Primer3 is a widely used program for selection of oligonucleotide primers for PCR. The websites used for implementation of Primer3 have recently been updated. PCR requires Mg(2+)(,) which has a significant dsDNA stabilizing effect that must be taken into account when designing PCR primers. The data sets and formulas used to correct for salt concentrations have been updated in Primer3 to give better prediction of melting temperature (T(m)). The liberal combination of different formulas for monovalent and divalent salt correction can lead to different results, depending on the formula chosen by the user. Using published T(m) for 475 different oligonucleotides, it is shown that the combination of the implemented conversion of divalent to monovalent cation concentration works well with one salt correction formula but not with an alternative one. Use of a more recently described alternative formula would lead to equally good T(m) predictions if divalent cations are present. The proper selection of compatible primer pairs depends on the choice of a good combination of salt correction formulas. Currently the SantaLucia salt correction formula should be used if Mg(2+) is present. The alternative formula should be updated to its recent form for future releases.     

Evolution of marital structure in 20 Sardinian villages from 1800 to 1974

This paper reports the data for endogamy, exogamy and consanguinity in 20 Sardinian villages from 1800 to 1974, divided into generational periods (25 years). The data are taken from the records of 48,262 marriages celebrated in the parishes of 20 Sardinian villages. The results for Spearman's rank correlation coefficient indicate significant positive correlations of endogamy with consanguinity (r(s) = 0.6551, percentage of consanguineous marriages; r(s) = 0.5477, Bernstein's alpha), with altitude (r(s) = 0.6386), with population size (r(s) = 0.2519), as well as a significant negative correlation of endogamy with time (r(s) = -0.4210). In addition, consanguinity shows a significant positive correlation with altitude (r(s) = 0.5717, percentage of consanguineous marriages; r(s) = 0.6295, Bernstein's alpha) and a significant negative correlation with time (r(s) = -0.2363, percentage of consanguineous marriages); the negative correlation between the mean level of consanguinity (alpha) and time is statistically not significant. In the single parishes, the levels of endogamy and consanguinity are generally higher in the mountain villages than in those situated in the hill and plain areas. Finally, it is hypothesized that, in the time period considered, the level of consanguinity and especially of endogamy were largely determined by isolation.