Statistical power of an exact test of Hardy-Weinberg proportions of genotypic data at a multiallelic locus

Are we underestimating the occurrence of sympatric populations?

Sympatric populations are conspecific populations that coexist spatially. They are of interest in evolutionary biology by representing the potential first steps of sympatric speciation and are important to identify and monitor in conservation management. Reviewing the literature pertaining to sympatric populations, we find that most cases of sympatry appear coupled to phenotypic divergence, implying ease of detection. In comparison, phenotypically cryptic, sympatric populations seem rarely documented. We explore the statistical power for detecting population mixtures from genetic marker data, using commonly applied tests for heterozygote deficiency (i.e., Wahlund effect) and the structure software, through computer simulations. We find that both tests are efficient at detecting population mixture only when genetic differentiation is high, sample size and number of genetic markers are reasonable and the sympatric populations happen to occur in similar proportions in the sample. We present an approximate expression based on these experimental factors for the lower limit of F_ST , beyond which power for structure collapses and only the heterozygote-deficiency tests retain some, although low, power. The findings suggest that cases of cryptic sympatry may have passed unnoticed in population genetic screenings using number of loci typical of the pre-genomics era. Hence, cryptic sympatric populations may be more common than hitherto thought, and we urge more attention being diverted to their detection and characterization.

Assessing the validity of asthma associations for eight candidate genes and age at diagnosis effects

BACKGROUND

Before the advent of genome-wide association studies (GWAS), ADAM33, ADRB2, CD14, MS4A2 (alias FCER1B), IL13, IL4, IL4R, and TNF constituted the most replicated non-HLA candidate genes with asthma and related traits. However, except for the IL13-IL4 region, none of these genes have been found in close proximity of genome-wide significant hits among GWAS for asthma or related traits. Here we aimed to assess the reproducibility of these asthma associations and to test if associations were more evident considering the effect of age at diagnosis.

METHODOLOGY/PRINCIPAL FINDINGS

We systematically evaluated 286 common single nucleotide polymorphisms (SNPs) of these 8 genes in a sample of 1,865 unrelated Spanish individuals (606 asthmatics and 1,259 controls). We found that variants at MS4A2, IL4R and ADAM33 genes demonstrated varying association effects with the age at diagnosis of asthma, with 10 SNPs showing study-wise significance after the multiple comparison adjustment. In addition, in silico replication with GWAS data supported the association of IL4R.

CONCLUSIONS/SIGNIFICANCE

Our results support the important role of MS4A2, IL4R and ADAM33 genes in asthma and/or atopy susceptibility. However, additional studies in larger samples sets are needed to firmly implicate these genes in asthma susceptibility, and also to identify the causal variation underlying the associations found.

A Monte Carlo permutation test for random mating using genome sequences

Testing for random mating of a population is important in population genetics, because deviations from randomness of mating may indicate inbreeding, population stratification, natural selection, or sampling bias. However, current methods use only observed numbers of genotypes and alleles, and do not take advantage of the fact that the advent of sequencing technology provides an opportunity to investigate this topic in unprecedented detail. To address this opportunity, a novel statistical test for random mating is required in population genomics studies for which large sequencing datasets are generally available. Here, we propose a Monte-Carlo-based-permutation test (MCP) as an approach to detect random mating. Computer simulations used to evaluate the performance of the permutation test indicate that its type I error is well controlled and that its statistical power is greater than that of the commonly used chi-square test (CHI). Our simulation study shows the power of our test is greater for datasets characterized by lower levels of migration between subpopulations. In addition, test power increases with increasing recombination rate, sample size, and divergence time of subpopulations. For populations exhibiting limited migration and having average levels of population divergence, the statistical power approaches 1 for sequences longer than 1Mbp and for samples of 400 individuals or more. Taken together, our results suggest that our permutation test is a valuable tool to detect random mating of populations, especially in population genomics studies.

A note on exact conditional and unconditional tests for Hardy-Weinberg equilibrium

The exact conditional approach is frequently used for testing Hardy-Weinberg equilibrium in population genetics. This approach respects the test size as compared to the traditionally used asymptotic approaches. It is a full-enumeration method and very computational. Many efficient algorithms have been successfully developed to implement this exact approach. An alternative to the conditional approach is the unconditional approach, which relaxes the restriction of the fixed number of allelic counts as in the conditional approach. The first unconditional test considered in this study is the one based on maximization, which has been shown to be more powerful than the conditional test to loci with two alleles for small sample sizes. By using the p value of the conditional approach as a test statistic in the following maximization step, the second unconditional test is uniformly more powerful than the conditional approach. We compared these exact tests based on three commonly used test statistics with regards to type I error rate and power. It is recommended to use the second unconditional approach in practice due to the power gain in the case with two alleles.

Exact tests for Hardy-Weinberg proportions

Exact conditional tests are often required to evaluate statistically whether a sample of diploids comes from a population with Hardy-Weinberg proportions or to confirm the accuracy of genotype assignments. This requirement is especially common when the sample includes multiple alleles and sparse data, thus rendering asymptotic methods, such as the common chi(2)-test, unreliable. Such an exact test can be performed using the likelihood ratio as its test statistic rather than the more commonly used probability test. Conceptual advantages in using the likelihood ratio are discussed. A substantially improved algorithm is described to permit the performance of a full-enumeration exact test on sample sizes that are too large for previous methods. An improved Monte Carlo algorithm is also proposed for samples that preclude full enumeration. These algorithms are about two orders of magnitude faster than those currently in use. Finally, methods are derived to compute the number of possible samples with a given set of allele counts, a useful quantity for evaluating the feasibility of the full enumeration procedure. Software implementing these methods, ExactoHW, is provided.

Deviations from Hardy-Weinberg proportions for multiple alleles under viability selection

Departures of genotype frequencies from Hardy-Weinberg proportions (HWP) for a single autosomal locus due to viability selection in a random mating population have been studied only for the two-allele case. In this article, the analysis of deviations from HWP due to constant viability selection is extended to multiple alleles. The deviations for an autosomal locus with k alleles are measured by means of k fii fixation indices for homozygotes and k(k-1)/2 fij fixation indices for heterozygotes, and expressions are obtained for these indices (FIS statistics) under the multiallele viability model. Furthermore, expressions for fii and fij when the multiallele polymorphism is at stable equilibrium are also derived and it is demonstrated that the pattern of multiallele Hardy-Weinberg deviations at equilibrium is characterized by a global heterozygote excess and a deficiency of each of the homozygotes. This pattern may be useful for detecting whether a given multiallelic polymorphism is at stable equilibrium in the population due to viability selection. An analysis of Hardy-Weinberg deviations from published data for the three-allele polymorphism at the beta-globin locus in human populations from West Africa is presented for illustration.

A microsatellite guided insight into the genetic status of adi, an isolated hunting-gathering tribe of northeast India

Tibeto-Burman populations of India provide an insight into the peopling of India and aid in understanding their genetic relationship with populations of East, South and Southeast Asia. The study investigates the genetic status of one such Tibeto-Burman group, Adi of Arunachal Pradesh based on 15 autosomal microsatellite markers. Further the study examines, based on 9 common microsatellite loci, the genetic relationship of Adi with 16 other Tibeto-Burman speakers of India and 28 neighboring populations of East and Southeast Asia. Overall, the results support the recent formation of the Adi sub-tribes from a putative ancestral group and reveal that geographic contiguity is a major influencing factor of the genetic affinity among the Tibeto-Burman populations of India.

Allele frequency distribution at 15 autosomal STR loci in Panggi, Komkar and Padam sub tribes of Adi, a Tibeto-Burman speaking population of Arunachal Pradesh, India

Fifteen autosomal STR loci were analyzed in 223 healthy individuals belonging to three remote, isolated Tibeto-Burman speaking sub tribes namely, Panggi, Komkar and Padam of Adi tribe of Arunachal Pradesh, India. The analyzed markers exhibited a high degree of polymorphism in the studied populations. Statistical parameters of forensic interest; observed heterozygosity, probability of homozygosity, exact test, likelihood ratio test, power of discrimination, power of exclusion, match probability and typical paternity index were determined for all loci. The average heterozygosity values were found to be low in the three populations (Panggi: 0.7747; Komkar: 0.7742 and Padam: 0.7663). The combined power of discrimination and power of exclusion were 0.9999 in the studied populations thereby revealing the high forensic significance of the chosen markers. The study indicates the utility of the tested microsatellite markers in forensic human identification, paternity testing and human population genetic studies.

Concordance of allozyme and microsatellite differentiation in a marine fish, but evidence of selection at a microsatellite locus

Previous studies have reported higher levels of divergence for microsatellites than for allozymes in several species, suggested to reflect stabilizing selection on the allozymes. We compared the differentiation patterns of 11 allozyme and nine microsatellite loci using 679 spawning Atlantic herring (Clupea harengus) collected in the Baltic and North Seas to test for differential natural selection on these markers. Observed distributions of F statistics for the two types of markers are conspicuously dissimilar, but we show that these differences can largely be explained by sampling phenomena caused by different allele frequency distributions and degrees of variability. The results show consistently low levels of differentiation for both marker types, with the exception of one outlier microsatellite locus with a notably high F(ST). The aberrant pattern at this locus is primarily due to two alleles occurring at markedly high frequencies in the Baltic, suggesting selection at this locus, or a closely linked one. When excluding this locus, the two marker types show similar, weak differentiation patterns with F(ST) values between the Baltic and the North Seas of 0.001 and 0.002 for allozymes and microsatellites, respectively. This small heterogeneity, and weak isolation by distance, is easier to distinguish statistically with microsatellites than with allozymes that have fewer alleles and skewed frequency distributions. The allozymes, however, also detect surprisingly low levels of divergence. Our results support suggestions that previously described differences between marker types are primarily caused by a small number of outlier loci.

Genetic diversity at 15 microsatellite loci among the Adi Pasi population of Adi tribal cluster in Arunachal Pradesh, India

Genetic polymorphisms at 15 tetrameric short tandem repeat (STR) loci were studied in 203 healthy individuals of Adi Pasi population from Arunachal Pradesh, India. All the loci analyzed were highly polymorphic and there was no significant deviation from the Hardy-Weinberg equilibrium (HWE) excepting D8S1179 and D18S51. Other forensic useful statistical parameters were also calculated and the 15 microsatellite markers selected for this study were found to be suitable for human identification and population genetic studies.

Genetic Polymorphism at 15 Tetrameric Short Tandem Repeat Loci in Four Aboriginal Tribal Populations of Bengal

POPULATIONS

This study reports the genetic polymorphism observed at 15 short tandem repeat loci D3S1358, TH01, D21S11, D18S51, D5S818, D13S317, D7S820, D16S539, CSF1PO, vWA, D8S1179, TPOX, D2S1338, D19S433, and FGA in four aboriginal populations of Bengal. The analysis was performed to decipher the suitability of CODIS as well as six other highly polymorphic and unlinked markers in Forensic Testing. Studied populations include four tribes: Karmali, Kora, Maheli, and Lodha.

Hardy–Weinberg equilibrium in genetic association studies: an empirical evaluation of reporting, deviations, and power

We evaluated the testing and reporting of Hardy-Weinberg equilibrium (HWE) in recent genetic association studies, detected how frequently HWE was violated and estimated the power for HWE testing in this literature. Genetic association studies published in 2002 in Nature Genetics, American Journal of Human Genetics, and American Journal of Medical Genetics were assessed. Data were analyzed on 239 biallelic associations using 154 distinct genotype distribution data sets where HWE could be tested. Any information on HWE was given only for 150 (62.8%) associations (92 (59.7%) data sets). Reanalysis of the data showed significant deviation from HWE in the disease-free controls of 20 associations (13 data sets), but only four of them (two data sets) were admitted in the published articles. Another four deviations (in two data sets) were observed in the combined sample of cases and controls of studies where both cases and controls were diseased, and none were reported in the papers. In all six tested multiallelic associations (six data sets), there was violation of HWE, but this was not admitted in the published articles. Power calculations showed that most studies conforming to HWE simply were largely underpowered to detect HWE deviation; for example, power to detect an inbreeding of magnitude F=0.10 exceeded 80% in only 11 (7%) of the data sets being tested. This empirical evidence suggests that, even in high profile genetics journals, testing and reporting for HWE is often neglected and deviations are rarely admitted in the published reports. Moreover, power is limited for HWE testing in most current genetic association studies.

Characterization of Pacific golden chanterelle (Cantharellus formosus) genet size using co-dominant microsatellite markers

We characterized five co-dominant microsatellite markers and used them to study Pacific golden chanterelle (Cantharellus formosus) genet size and its relation to forest age and disturbance. Fruit-bodies were mapped in and collected from nine replicate study plots in old-growth, recently thinned, and unthinned 40-60-year-old second-growth stands dominated by Douglas fir (Pseudotsuga menziesii). Information from microsatellite loci, combined with random fragment length polymorphism analysis of the nuclear DNA internal transcribed spacer indicates that putative 'C. formosus' fruit-body collections may include a cryptic chanterelle species. Small genets were characterized for both genetic types with mean maximum widths of 3.2 +/- 3.6 m for C. formosus and 1.5 +/- 1.7 m for the alternative genetic group. Variance in genet size was high and some multilocus genotypes were observed on multiple plots separated by 0.3 km or more, indicating that genets were not fully resolved by the loci described here. There was no evidence that genet size differed across the three disturbance treatments.

Global genetic variation at nine short tandem repeat loci and implications on forensic genetics

We have studied genetic variation at nine autosomal short tandem repeat loci in 20 globally distributed human populations defined by geographic and ethnic origins, viz., African, Caucasian, Asian, Native American and Oceanic. The purpose of this study is to evaluate the utility and applicability of these nine loci in forensic analysis in worldwide populations. The levels of genetic variation measured by number of alleles, allele size variance and heterozygosity are high in all populations irrespective of their effective sizes. Single- as well as multi-locus genotype frequencies are in conformity with the assumptions of Hardy-Weinberg equilibrium. Further, alleles across the entire set of nine loci are mutually independent in all populations. Gene diversity analysis shows that pooling of population data by major geographic groupings does not introduce substructure effects beyond the levels recommended by the National Research Council, validating the establishment of population databases based on major geographic and ethnic groupings. A network tree based on genetic distances further supports this assertion, in which populations of common ancestry cluster together. With respect to the power of discrimination and exclusion probabilities, even the relatively reduced levels of genetic variation at these nine STR loci in smaller and isolated populations provide an exclusionary power over 99%. However, in paternity testing with unknown genotype of the mother, the power of exclusion could fall below 80% in some isolated populations, and in such cases use of additional loci supplementing the battery of the nine loci is recommended.

Genetic polymorphism at 15 STR loci among three important subpopulation of Bihar, India

Genotype polymorphism studies at 15 highly polymorphic short tandem repeat (STR) loci were carried out in three genetically important minor caste groups (Yadav, Kurmi and Baniya) of Bihar, a eastern state of India to evaluate their significance in human identification and population genetics study. The selected communities practice endogamy. Despite of same geographical area, the physical features of Yadavs and Baniyas resemble North Indian Indo-Caucasoids whereas Kurmis resemble more to Indo-Austroloids. Among the chosen 15 loci, two are penta-nucleotide repeat: Penta-D and Penta-E, and 13 are tetra-nucleotide repeat: vWA D8S1179, TPOX, FGA, D5S818, D13S317, D7S820, D16S539, D3S1358, THO1, CSF1PO, D21S11, D18S51 and are validated for other population of India and world for forensic testing and human population study. Thirteen of these STR loci are present in the combined DNA index system (CODIS) [J. Forensic Sci. 44 (1999) 1277] and world-wide data is available.

Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis

Homoplasy has recently attracted the attention of population geneticists, as a consequence of the popularity of highly variable stepwise mutating markers such as microsatellites. Microsatellite alleles generally refer to DNA fragments of different size (electromorphs). Electromorphs are identical in state (i.e. have identical size), but are not necessarily identical by descent due to convergent mutation(s). Homoplasy occurring at microsatellites is thus referred to as size homoplasy. Using new analytical developments and computer simulations, we first evaluate the effect of the mutation rate, the mutation model, the effective population size and the time of divergence between populations on size homoplasy at the within and between population levels. We then review the few experimental studies that used various molecular techniques to detect size homoplasious events at some microsatellite loci. The relationship between this molecularly accessible size homoplasy size and the actual amount of size homoplasy is not trivial, the former being considerably influenced by the molecular structure of microsatellite core sequences. In a third section, we show that homoplasy at microsatellite electromorphs does not represent a significant problem for many types of population genetics analyses realized by molecular ecologists, the large amount of variability at microsatellite loci often compensating for their homoplasious evolution. The situations where size homoplasy may be more problematic involve high mutation rates and large population sizes together with strong allele size constraints.

Heterozygote deficiencies in small lacustrine populations of brook charr Salvelinus Fontinalis Mitchill (Pisces, Salmonidae): a test of alternative hypotheses

Empirical studies of natural populations have commonly reported departures from Hardy-Weinberg expected proportions of heterozygote individuals. Recent advances in statistical population genetics now offer the potential to exploit individual multilocus genotypic information to test more rigorously for possible sources of heterozygote deficiencies. In a previous study in lacustrine brook charr (Salvelinus fontinalis), we reported stronger deficits in small than in large lakes. In the present paper, we propose a methodology for empirically testing alternative hypotheses to identify the cause of the deficits observed in three of the smallest lakes (85, 109 and 182 ha) analysed. First, as in several salmonid species, brook charr may exhibit a trophic polymorphism in north temperate lakes. If morphs are genetically divergent, indiscriminate sampling of both forms would result in less heterozygote individuals than expected in a randomly mating population (Wahlund effect). Using an individual-based method aiming at detecting cryptic population structure, we can reject this explanation as the sole source of deficits for all three lakes. Secondly, mating among relatives could also be frequent in small lakes and lead to heterozygote deficiencies. Significantly more fish than expected at random had low individual multilocus heterozygosity in two of the lakes, suggesting that inbred fish may have been present. Thirdly, sampling of genetically related fish would also lead to departures from Hardy-Weinberg proportions. In the same two lakes, the distribution of pairwise relatedness coefficients departed from its random expectation, suggesting that non-random sampling of kin may have occurred.

Population structure, stepwise mutations, heterozygote deficiency and their implications in DNA forensics

In a substructured population the overall heterozygote deficiency can be predicted from the number of subpopulations (s), their time of divergence (t), and the nature of the mutations. At present the true mutational mechanisms at the hypervariable DNA loci are not known. However, the two existing mutation models (the infinite allele model (IAM) and the stepwise mutation model (SMM)) provide some guides to predictions from which the possible effect of population substructuring may be evaluated, assuming that the subpopulations do not exchange any genes among them during evolution. The theory predicts that the loci with larger mutation rate, and consequently showing greater heterozygosity within subpopulations, should exhibit a smaller proportional heterozygote deficiency (GST) and, hence, the effects of population substructuring should be minimal at the hypervariable DNA loci (an order of magnitude smaller than that at the blood group and protein loci). Applications of this theory to data on six Variable Number of Tandem Repeat (VNTR) loci and five short tandem repeat (STR) loci in the major cosmopolitan populations of the USA show that while the VNTR loci often exhibit a large significant heterozygote deficiency, the STR loci do not show a similar tendency. This discordant finding may be ascribed to the limitations, coalescence and nondetectability of alleles associated with the restriction fragment length polymorphism (RFLP) analysis through which the VNTR loci are scored. Such limitations do not apply to the polymerase chain reaction (PCR) method, through which the STR loci are scored. The implications of these results are discussed in the context of the forensic use of DNA typing data.

Kin selection, social structure, gene flow, and the evolution of chimpanzees

Hypotheses about chimpanzee social behavior, phylogeography, and evolution were evaluated by noninvasive genotyping of free-ranging individuals from 20 African sites. Degrees of relatedness among individuals in one community were inferred from allele-sharing at eight nuclear simple sequence repeat (SSR) loci. Males are related on the order of half-siblings, and homozygosity is significantly increased at several SSR loci compared to Hardy-Weinberg expectations. These data support the kin-selection hypothesis for the evolution of cooperation among males. Sequence variation patterns at two mitochondrial loci indicate historically high long-distance gene flow and clarify the relationships among three allopatric subspecies. The unexpectedly large genetic distance between the western subspecies, Pan troglodytes verus, and the other two subspecies suggests a divergence time of about 1.58 million years. This result, if confirmed at nuclear loci and supported by eco-behavioral data, implies that P. t. verus should be elevated to full species rank.

Likelihood ratios for DNA identification

Likelihood ratio (LR) tests are provided for the three alternatives to DNA identity: exclusion, coincidence, and kinship. The coincidence test uses the radius of coalescence to conserve the observed frequency of single band phenotypes. Genotype probabilities under kinship are derived for mating groups, specified relatives, and structured populations; and unbiased estimates of the genetic parameters are provided. The LR is made robust to gene frequency errors by specifying the mean matching probability, and the tolerable loss of information this entails is determined by LR theory. This straightforward application of the seminal work of Jerzy Neyman and Sewall Wright strongly supports the use of LRs and kinship for presentation of DNA evidence by expert witnesses and committees.