Interethnic and intraethnic variability of CYP2C8 and CYP2C9 polymorphisms in healthy individuals

Cytochrome P450 (CYP) superfamily members CYP2C8 and CYP2C9 are polymorphically expressed enzymes that are involved in the metabolic inactivation of several drugs, including, among others, antiepileptics, NSAIDs, oral hypoglycemics, and anticoagulants. Many of these drugs have a narrow therapeutic index, and growing evidence indicates a prominent role of CYP2C8 and CYP2C9 polymorphisms in the therapeutic efficacy and in the development of adverse effects among patients treated with drugs that are CYP2C8 or CYP2C9 substrates. In this review, we summarize present knowledge on human variability in the frequency of variant CYP2C8 and CYP2C9 alleles. Besides an expected interethnic variability in allele frequencies, a large intraethnic variability exists. Among Asian subjects, for example, statistically significant differences (p < 0.0001) in CYP2C9*3 allele frequencies between Chinese and Japanese individuals have been reported. In addition, individuals from East Asia present different allele frequencies for CYP2C9*2 and CYP2C9*3 compared with South Asian subjects (p < 0.0001). Among Caucasian Europeans, statistically significant differences for the frequency of CYP2C8*3, CYP2C9*2, and CYP2C9*3 exist (p < 0.0001). This indicates that Asian individuals or Caucasian European individuals cannot be considered as homogeneous groups regarding CYP2C8 or CYP2C9 allele frequencies. Caucasian American subjects also show a large variability in allele frequencies, which is likely to be related to ethnic ancestry. A higher frequency of variant CYP2C8 and CYP2C9 alleles is expected among Caucasian Americans with South European ancestry than in individuals with North European ancestry. The findings summarized in this review suggest that among individuals with Asian or European ancestry, intraethnic differences in the risk of developing adverse effects with drugs that are CYP2C8 or CYP2C9 substrates are to be expected. In addition, the observed intraethnic variability reinforces the need for proper selection of control subjects and points against the use of surrogate control groups for studies involving association of CYP2C8 or CYP2C9 alleles with adverse drug reactions or spontaneous diseases.     

D4 dopamine-receptor (DRD4) alleles and novelty seeking in substance-dependent, personality-disorder, and control subjects.

Two reports have been published suggesting an association between the personality trait of novelty seeking and the DRD4*7R allele at the D4 dopamine-receptor locus (with heterozygotes or homozygotes for DRD4*7R having higher novelty seeking). We studied novelty seeking and four coding-sequence polymorphisms affecting protein structure in the D4 dopamine-receptor gene (DRD4) in a sample of 341 American subjects, of whom 224 are of primarily European ancestry and 117 are of primarily African ancestry. These subjects had diagnoses of substance dependence or personality disorder (PD) or were screened to exclude major psychiatric diagnosis. We found that, although the substance-dependent subjects had significantly higher novelty seeking than the control and PD subjects, they did not differ in DRD4*7R allele frequency. There was no association between any DRD4 polymorphism and novelty seeking in any population or diagnostic group, except for a significant association between the DRD4*7R allele and lower novelty seeking among European American females and African American substance abusers. The novelty seeking of subjects heterozygous for a null mutation did not differ from that of subjects with two functional alleles. We conclude that the most likely explanation of these results is that the DRD4 VNTR does not influence directly the trait of novelty seeking, in these samples.     

High frequency of CYP1A1*2C allele in Brazilian populations

The genetic variability of the CYP1A1 I462V polymorphism (CYP1A1*2C) was investigated in four Brazilian populations: three groups of African descent and one group of European descent. The CYP1A1 polymorphism was analyzed by two different procedures, first by the allele-specific polymerase chain reaction (PCR) method and then by the PCR-restricted fragment length polymorphism (PCR-RFLP) method before digestion with BsrDI. The frequency of CYP1A1 *2C was 11% in Brazilians of European descent, a frequency that is slightly higher but not statistically different from that observed in European populations. In Brazilians of African ancestry this value was very high (12% to 15%). This allele was not observed in the only two African populations investigated thus far. By themselves, the two factors of interethnic admixture (with populations of European descent and/or Amerindian populations) and genetic drift cannot explain the high values observed here. Our findings suggest that the CYP1A1 *2C allele may possibly be present in Africa, but restricted to some ethnic groups not yet investigated. Environmental factors in South America might also have acted as selective factors increasing the CYP1A1 *2C gene frequency. Our data also suggest that the CYP1A1 *2C allele might possibly have originated in Africa.     

The genetic architecture of selection at the human dopamine receptor D4 (DRD4) gene locus

Associations of the seven-repeat (7R) allele of the human dopamine receptor D4 (DRD4) gene with both the personality trait of novelty seeking and attention deficit/hyperactivity disorder have been reported. Recently, on the basis of the unusual DNA sequence organization of the DRD4 7R 48-bp tandem repeat (VNTR), we proposed that the 7R allele originated as a rare mutational event that increased to high frequency by positive selection. We now have resequenced the entire DRD4 locus from 103 individuals homozygous for 2R, 4R, or 7R variants of the VNTR, a method developed to directly estimate haplotype diversity. DNA from individuals of African, European, Asian, North and South American, and Pacific Island ancestry were used. 4R/4R homozygotes exhibit little linkage disequilibrium (LD) over the region examined, with more polymorphisms observed in DNA samples from African individuals. In contrast, the evidence for strong LD surrounding the 7R allele is dramatic, with all 7R/7R individuals (including those from Africa) exhibiting the same alleles at most polymorphic sites. By intra-allelic comparison at 18 high-heterozygosity sites spanning the locus, we estimate that the 7R allele arose prior to the upper Paleolithic era (approximately 40000-50000 years ago). Further, the pattern of recombination at these polymorphic sites is the pattern expected for selection acting at the 7R VNTR itself, rather than at an adjacent site. We propose a model for selection at the DRD4 locus consistent with these observed LD patterns and with the known biochemical and physiological differences between receptor variants.     

Acid phosphatase, adenosine deaminase and esterase D polymorphisms in the Spanish Basque population.

The 3 red-cell polymorphic systems acid phosphatase (ACP), adenosine deaminase (ADA) and esterase D (ESD) have been studied in a random sample of 1,112 individuals from the Basque country: The allelic frequencies obtained were ACP*A = 0.275, ACP*B = 0.718 and ACP*C = 0.007; ADA*2 = 0.021, and, ESD*2 = 0.066. The allelic frequencies have been compared with those of other Basque and other European populations. In comparison with Basques, significant differences were detected only for ACP, whereas as regards other Europeans significant differences were obtained with practically all the populations compared for the 3 genetic systems studied. The low values of the less frequent alleles, especially that for the ACP*C allele which is the lowest reported in Europe, are noteworthy.     

HLA-DRB1 Alleles Are Associated with the Susceptibility to Sporadic Parkinson’s Disease in Chinese Han Population

Immune disorders may play an important role in the pathogenesis of Parkinson''s disease (PD). Recently, polymorphisms in the HLA-DR region have been found to be associated with sporadic PD in European ancestry populations. However, polymorphisms in the HLA complex are highly variable with ethnic and geographic origin. To explore the relationships between polymorphisms of the HLA-DR region and sporadic PD in Chinese Han population, we genotyped 567 sporadic PD patients and 746 healthy controls in two independent series for the HLA-DRB1 locus with Polymerase chain reaction-sequence based typing(PCR-SBT). The χ² test was used to evaluate the distribution of allele frequencies between the patients and healthy controls. The impact of HLA-DRB1 alleles on PD risk was estimated by unconditional logistic regression. We found a significant higher frequency of HLA-DRB1*0301 in sporadic PD patients than in healthy controls and a positive association, which was independent of onset age, between HLA-DRB1*0301 and PD risk. Conversely, a lower frequency of HLA-DRB1*0406 was found in sporadic PD patients than in healthy controls, with a negative association between HLA-DRB1*0406 and PD risk. Furthermore, a meta-analysis involving 195205 individuals was conducted to summarize the frequencies of these two alleles in populations from various ethnic regions, we found a higher frequency of HLA-DRB1*0301, but a lower frequency of HLA-DRB1*0406 in European ancestry populations than that in Asians, this was consistent with the higher prevalence of sporadic PD in European ancestry populations. Based on these results, we speculate that HLA-DRB1 alleles are associated with the susceptibility to sporadic PD in Chinese Han population, among them HLA-DRB1*0301 is a risk allele while the effect of HLA-DRB1*0406 deserves debate.     

Effect of ACP1*C on early life viability

Associations with past malarial morbidity, season of conception, and common diseases such as obesity, type 2 diabetes, and allergy argue against neutrality of the ACP1 genetic polymorphism. Comparison of ACP1 distribution in mothers and their newborns and analysis of the joint wife-husband ACP1 phenotype distribution in couples with repeated spontaneous abortion suggest a negative effect of the ACP1*C allele on early life viability. Analysis of the polymorphism of the ACP1 gene suggests that, unlike the ACP1*A and ACP1*B alleles, the ACP1*C allele is independent of sequences in the 5' flanking region, resulting in an inverted F/S isoform ratio.     

Analysis of CYP2C9*2, CYP2C9*3 and VKORC1 ‐1639 G>A polymorphisms in a population from South‐Eastern Europe

The CYP2C9 enzyme metabolizes a wide range of relevant drugs, among which are oral anticoagulants. VKORC1 is the pharmacodynamic target of the oral anticoagulants. The genetic polymorphisms CYP2C9*2, CYP2C9*3 and VKORC1 ‐1639 G>A are the major determinants of the inter‐individual variability in the dosage requirements of oral anticoagulants. This study provides a first evaluation of these 3 polymorphisms in a Romanian population. A total of 332 Romanian individuals were genotyped for the CYP2C9*2, CYP2C9*3 and VKORC1 ‐1639 G>A polymorphisms using the PCR‐RFLP technique. Sixty‐two individuals (18.7%) were heterozygous for CYP2C9*2, whereas 47 individuals (14.1%) were heterozygous for CYP2C9*3. Fourteen individuals (4.2%) had a CYP2C9*2 homozygous, CYP2C9*3 homozygous or CYP2C9*2/CYP2C9*3 compound heterozygous genotype. These individuals are predicted to have the lowest CYP2C9 enzymatic activity. The allele frequencies of the CYP2C9*2 and CYP2C9*3 polymorphisms were 11.3% and 9.3% respectively. For the VKORC1 ‐1639 G>A polymorphism, there were 170 heterozygotes (51.2%) and 55 (16.6%) homozygotes for the A allele. The frequency of the A allele was 42.2%. Overall, the distribution of the CYP2C9*2, CYP2C9*3 and VKORC1 ‐1639 G>A polymorphisms observed in our cohort is in accordance with other Caucasian populations. A large number of Romanians are expected to harbour at least one CYP2C9 variant allele and/or one VKORC1 ‐1639 G>A allele. This frequency has major implications in the pharmacogenomics of oral anticoagulants in Romanians.     

Silencing of duplicate genes: a null allele polymorphism for lactate dehydrogenase in brown trout (Salmo trutta)

A previously described isozyme polymorphism at one of two skeletal muscle LdhA loci in brown trout is due to a null allele, Ldh1(n), producing no detectable catalytic activity. Homozygotes for this allele have approximately only 56% of the LDH activity in skeletal muscle relative to homozygotes for the active allele. The remaining activity results from enzyme subunits produced by other LDH loci. The Ldh1(n) allele is common and widespread throughout brown trout populations in Sweden and is also found in populations from Ireland. The persistence of duplicate gene expression for the LdhA loci in almost all salmonid species is best explained by natural selection against individuals containing null alleles. However, there is no indication of natural selection against brown trout with the Ldh1(n) allele: We suggest that the selection against individuals containing null alleles that is apparently responsible for the persistence of duplicate LdhA loci in salmonids occurs only under certain environmental conditions.      

High domestic pig contribution to the local gene pool of free-living European wild boar: a case study in Poland

Rates of hybridization between wild and domesticated animals appear to be increasing worldwide. Recent results suggest that genetic introgression from domestic swine into European wild boar is much more common in local populations than expected, based on pan-European studies. Thus, we screened the genetic purity of 265 free-living wild boars from two hunting areas in Poland by genotyping the melanocortin receptor 1 gene (MC1R) for polymorphism. Unexpectedly, high numbers of individuals with domestic genes (24%) were identified. This suggests that mixed ancestry may be common in Polish wild boar. Among admixed individuals, backcrosses with domestic pig and/or introgressed wild boars were detected (2%). Multiple commercial domestic pig breeds are possibly involved in the introgression observed in the study populations. In addition, the absence of significant differences in the frequency of wild-type allele among two hunting areas suggests high dispersal of individuals and gene flow among populations. We conclude that further study is needed to better understand the mechanisms and sources of introgression in wild boars in Poland.     

Haptoglobin 1F allele frequency is high among indigenous populations in the state of Durango, Mexico

We studied haptoglobin polymorphism in Mexican populations with high Indian ancestry living in isolated and urban areas in the state of Durango. Analysis with respect to the HP*1F and HP*1S allelic subtypes by isoelectric focusing showed unusually high HP*1F allele frequencies among urban (0.370) and isolated Mexican Indians (0.383). Comparison with other population studies demonstrated a geographical cline of the HP*1F allele increasing in the same direction of the HP*1 allele, while HP*1S frequency does not show racial differences.     

Haplotypes of single cancer driver genes and their local ancestry in a highly admixed long-lived population of Northeast Brazil

Admixed populations have not been examined in detail in cancer genetic studies. Here, we inferred the local ancestry of cancer-associated single nucleotide polymorphisms (SNPs) and haplotypes of a highly admixed Brazilian population. SNP array was used to genotype 73 unrelated individuals aged 80-102 years. Local ancestry inference was performed by merging genotyped regions with phase three data from the 1000 Genomes Project Consortium using RFmix. The average ancestry tract length was 9.12-81.71 megabases. Strong linkage disequilibrium was detected in 48 haplotypes containing 35 SNPs in 10 cancer driver genes. All together, 19 risk and eight protective alleles were identified in 23 out of 48 haplotypes. Homozygous individuals were mainly of European ancestry, whereas heterozygotes had at least one Native American and one African ancestry tract. Native-American ancestry for homozygous individuals with risk alleles for HNF1B, CDH1, and BRCA1 was inferred for the first time. Results indicated that analysis of SNP polymorphism in the present admixed population has a high potential to identify new ancestry-associated alleles and haplotypes that modify cancer susceptibility differentially in distinct human populations. Future case-control studies with populations with a complex history of admixture could help elucidate ancestry-associated biological differences in cancer incidence and therapeutic outcomes.     

Several Different Lactase Persistence Associated Alleles and High Diversity of the Lactase Gene in the Admixed Brazilian Population

Adult-type hypolactasia is a common phenotype caused by the lactase enzyme deficiency. The −13910 C>T polymorphism, located 14 Kb upstream of the lactase gene (LCT) in the MCM6 gene was associated with lactase persistence (LP) in Europeans. This polymorphism is rare in Africa but several other variants associated with lactase persistence were observed in Africans. The aims of this study were to identify polymorphisms in the MCM6 region associated with the lactase persistence phenotype and to determine the distribution of LCT gene haplotypes in 981 individuals from North, Northeast and South Brazil. These polymorphisms were genotyped by PCR based methods and sequencing. The −13779*C,−13910*T, −13937*A, −14010*C, −14011*T LP alleles previously described in the MCM6 gene region that acts as an enhancer for the LCT gene were identified in Brazilians. The most common LP allele was −13910*T. Its frequency was highly correlated with European ancestry in the Brazilian populations investigated. The −13910*T was higher (0.295) in southern Brazilians of European ancestry and lower (0.175) in the Northern admixed population. LCT haplotypes were derived from the 10 LCT SNPs genotyped. Overall twenty six haplotypes previously described were identified in the four Brazilian populations studied. The Multidimensional Scaling analysis showed that Belém, in the north, was closer to Amerindians. Northeastern and southern Afro-descendants were more related with Bantu-speaking South Africans whereas the Southern population with European ancestry grouped with Southern and Northern Europeans. This study shows a high variability considering the number of LCT haplotypes observed. Due to the highly admixed nature of the Brazilian populations, the diagnosis of hypolactasia in Brazil, based only in the investigation of the −13910*T allele is an oversimplification.     

Differential performance among LDH-B genotypes in Rana lessonae tadpoles

The European pool frog, Rana lessonae, is widely polymorphic for two common alleles (b,e) at the lactate dehydrogenase-B (LDH-B) locus. We compared fitness-related larval life-history traits among LDH-B genotypes, which originated from segregation in heterozygous parents, in an artificial pond experiment where tadpoles of R. lessonae from a Swiss population were raised together with tadpoles of the hemiclonal hybrid R. esculenta at two densities. In R. lessonae, LDH-B e/e homozygotes at each density had a higher proportion of metamorphs among survivors, reached metamorphosis earlier, and were heavier at metamorphosis than b/b homozygotes; b/e heterozygotes had intermediate values. That e/e individuals were superior to b/b in both time to and mass at metamorphosis is surprising because these two life-history traits are thought to reflect a performance trade-off; e/e genotypes apparently compensated for shorter time to metamorphosis by a higher growth rate. The two alleles showed the same performance ranking when combined in hybrids with a R. ridibunda allele: When R. esculenta from Swiss populations reared in the same ponds had received the e allele rather than the b allele from their R. lessonae parent, they reached metamorphosis earlier, but did not differ in mass at metamorphosis. The degree of linkage disequilibrium in the source population of the eight R. lessonae used as parents of the R. lessonae tadpoles is unknown, so we cannot exclude the possibility that the performance differences are caused by some anonymous tightly linked gene, rather than the LDH-B locus, that constitutes the genomically localized target of natural selection. A causal involvement of LDH-B is plausible, nevertheless, because this enzyme takes part in the central energy-metabolizing processes and has been reported to underlie fitness differences in other animals; also, differential performance of LDH-B genotypes has been observed in R. lessonae larvae from another population. The present results suggest strong directional selection for allele e; the sum of available data, including an independent laboratory experiment, suggests that partial environment-dependent overdominance combined with balancing selection favoring e/e homozygotes under some and b/b homozygotes under other conditions may be partially responsible for the broad maintenance of the LDH-B polymorphism in R. lessonae.     

Effects of a Locus Affecting Floral Pigmentation in Ipomoea Purpurea on Female Fitness Components

A locus influencing floral pigment intensity in the morning glory, Ipomoea purpurea, is polymorphic throughout the southeastern United States. Previous work has suggested that the white allele at this locus has a transmission advantage during mating because of the effect of flower color on pollinator behavior. The experiment described here was designed to determine whether other effects of the W locus may contribute an opposing selective advantage to the dark allele. Dark homozygotes were vegetatively smaller and produced fewer flowers, seed capsules and seeds than either light heterozygotes or white homozygotes. In addition, dark homozygotes produced smaller seeds than heterozygotes, and there is some indication that white homozygotes also produced smaller seeds than heterozygotes. Pleiotropic effects on seed number thus do not seem to contribute to selection opposing the mating advantage associated with the white allele. However, pleiotropic effects on seed size might contribute to overdominance that could stabilize the W locus polymorphism.     

Ancestral components of admixed genomes in a Mexican cohort

For most of the world, human genome structure at a population level is shaped by interplay between ancient geographic isolation and more recent demographic shifts, factors that are captured by the concepts of biogeographic ancestry and admixture, respectively. The ancestry of non-admixed individuals can often be traced to a specific population in a precise region, but current approaches for studying admixed individuals generally yield coarse information in which genome ancestry proportions are identified according to continent of origin. Here we introduce a new analytic strategy for this problem that allows fine-grained characterization of admixed individuals with respect to both geographic and genomic coordinates. Ancestry segments from different continents, identified with a probabilistic model, are used to construct and study "virtual genomes" of admixed individuals. We apply this approach to a cohort of 492 parent-offspring trios from Mexico City. The relative contributions from the three continental-level ancestral populations-Africa, Europe, and America-vary substantially between individuals, and the distribution of haplotype block length suggests an admixing time of 10-15 generations. The European and Indigenous American virtual genomes of each Mexican individual can be traced to precise regions within each continent, and they reveal a gradient of Amerindian ancestry between indigenous people of southwestern Mexico and Mayans of the Yucatan Peninsula. This contrasts sharply with the African roots of African Americans, which have been characterized by a uniform mixing of multiple West African populations. We also use the virtual European and Indigenous American genomes to search for the signatures of selection in the ancestral populations, and we identify previously known targets of selection in other populations, as well as new candidate loci. The ability to infer precise ancestral components of admixed genomes will facilitate studies of disease-related phenotypes and will allow new insight into the adaptive and demographic history of indigenous people.     

Short alleles revealed by PCR demonstrate no heterozygote deficiency at minisatellite loci D1S7, D7S21, and D12S11.

Short VNTR alleles that go undetected after conventional Southern blot hybridization may constitute an alternative explanation for the heterozygosity deficiency observed at some minisatellite loci. To examine this hypothesis, we have employed a screening procedure based on PCR amplification of those individuals classified as homozygotes in our databases for the loci D1S7, D7S21, and D12S11. The results obtained indicate that the frequency of these short alleles is related to the heterozygosity deficiency observed. For the most polymorphic locus, D1S7, approximately 60% of those individuals previously classified as homozygotes were in fact heterozygotes for a short allele. After the inclusion of these new alleles, the agreement between observed and expected heterozygosity, along with other statistical tests employed, provide additional evidence for lack of population substructuring. Comparisons of allele frequency distributions reveal greater differences between racial groups than between closely related populations.     

Investigation of the distribution of haptoglobin genotypes in human populations

Distribution of haptoglobin genotypes in a number of population specimens of the main World regions is studied. New statistics permitting to reveal deviation from the Hardy - Weinberg equilibrium and its direction (towards either deficiency or excess of heterozygotes) are proposed. S-statistics, based on the model of selection against homozygotes, can be used for quantitative estimation of the deviation. In this way, statistically significant excess of heterozygotes was found in European, North Asia and North American populations. In some areas, the coefficient of selection (S-statistics) is about 10%. Additional data which reflect different adaptive values of haptoglobin genotypes permit to suppose that the existence of haptoglobin polymorphism is provided by both higher fertility and lower viability of homozygotes Hp 2-2.     

Frequency of HFE gene mutations C282Y and H63D in Bosnia and Herzegovina

Genetic epidemiology studies of hereditary hemochromatosis (HHC) have shown a high prevalence of the C282Y mutation in individuals of the North Western European origin, whereas lower prevalence of HFE gene mutations was detected in the populations from southern European countries. However, no HFE mutation prevalence data have been provided for the population of Bosnia-Herzegovina so far. Therefore, the aim of this study was to determine the frequency of the C282Y and H63D HFE gene mutations in the population of Bosnia-Herzegovina. Among 200 analysed subjects 8 (4%) were C282Y heterozygotes; no C282Y homozygotes were found. The frequency of the H63D allele was 11.5%. There were 33 (16.5%) heterozygotes and 6 (3%) homozygotes for the H63D mutation. One (0.5%) compound heterozygote C282Y/H63D was identified. The observed C282Y and H63D allele frequency was 2.25% (95% confidence interval: 1.2-4.2) and 11.5% (95% confidence interval: 8.7-14.9), respectively. The prevalence of the C282Y and H63D mutations was estimated in Bosnia-Herzegovina, which fit well in the European northwest-to-southeast gradient of the C282Y mutation distribution. In addition, these results have an important implication for clinical evaluation of HHC in Bosnia-Herzegovina.     

Use of microsatellite loci to classify individuals by relatedness

This study investigates the use of microsatellite loci for estimating relatedness between individuals in wild, outbred, vertebrate populations. We measured allele frequencies at 20 unlinked, dinucleotide-repeat microsatellite loci in a population of wild mice ( Mus musculus ), and used these observed frequencies to generate the expected distributions of pairwise relatedness among full sib, half sib, and unrelated pairs of individuals, as would be estimated from the microsatellite data. In this population one should be able to discriminate between unrelated and full-sib dyads with at least 97% accuracy, and to discriminate half-sib pairs from unrelated pairs or from full-sib pairs with better than 80% accuracy. If one uses the criterion that parent-offspring pairs must share at least one allele per locus, then only 15% of full-sib pairs, 2% of half-sib pairs, and 0% of unrelated pairs in this population would qualify as potential parent-offspring pairs. We verified that the simulation results (which assume a random mating population in Hardy-Weinberg and linkage equilibrium) accurately predict results one would obtain from this population in real life by scoring laboratory-bred full- and half-sib families whose parents were wild-caught mice from the study population. We also investigated the effects of using different numbers of loci, or loci of different average heterozygosities ( He ), on misclassification frequencies. Both variables have strong effects on misclassification rate. For example, it requires almost twice as many loci of He = 0.62 to achieve the same accuracy as a given number of loci of He = 0.75. Finally, we tested the ability of UPGMA clustering to identify family groups in our population. Clustering of allele matching scores among the offspring of four sets of independent maternal half sibships (four females, each mated to two different males) perfectly recovered the true family relationships.