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

Multiallelic Positions in the Human Genome: Challenges for Genetic Analyses

As the amount of human genomic sequence available from personal genomes and exomes has increased, so too has the observation of genomic positions having two or more alternative alleles, so-called multiallelic sites. For portions of the haploid genome that are present in more than one copy, including segmental duplications, variation at such multisite variant positions becomes even more complex. Despite the frequency of multiallelic variants, a number of commonly used resources and tools in genomic research and diagnostics do not support these multiallelic variants all together or require special modifications. Here, we explore the frequency of multiallelic sites in large samples with whole exome sequencing and discuss potential outcomes of failing to account for multiple variant alleles. We also briefly discuss some commonly utilized resources that fully support multiallelic sites.

Association between immunoglobulin GM and KM genotypes and placental malaria in HIV-1 negative and positive women in western Kenya

Immunoglobulin (Ig) GM and KM allotypes, genetic markers of γ and κ chains, are associated with humoral immune responsiveness. Previous studies have shown the relationships between GM6-carrying haplotypes and susceptibility to malaria infection in children and adults; however, the role of the genetic markers in placental malaria (PM) infection and PM with HIV co-infection during pregnancy has not been investigated. We examined the relationship between the gene polymorphisms of Ig GM6 and KM allotypes and the risk of PM infection in pregnant women with known HIV status. DNA samples from 728 pregnant women were genotyped for GM6 and KM alleles using polymerase chain reaction-restriction fragment length polymorphism method. Individual GM6 and KM genotypes and the combined GM6 and KM genotypes were assessed in relation to PM in HIV-1 negative and positive women, respectively. There was no significant effect of individual GM6 and KM genotypes on the risk of PM infection in HIV-1 negative and positive women. However, the combination of homozygosity for GM6(+) and KM3 was associated with decreased risk of PM (adjusted OR, 0.25; 95% CI, 0.08-0.8; P = 0.019) in HIV-1 negative women while in HIV-1 positive women the combination of GM6(+/-) with either KM1-3 or KM1 was associated with increased risk of PM infection (adjusted OR, 2.10; 95% CI, 1.18-3.73; P = 0.011). Hardy-Weinberg Equilibrium (HWE) tests further showed an overall significant positive F(is) (indication of deficit in heterozygotes) for GM6 while there was no deviation for KM genotype frequency from HWE in the same population. These findings suggest that the combination of homozygous GM6(+) and KM3 may protect against PM in HIV-1 negative women while the HIV-1 positive women with heterozygous GM6(+/-) combined with KM1-3 or KM1 may be more susceptible to PM infection. The deficit in heterozygotes for GM6 further suggests that GM6 could be under selection likely by malaria infection.

Large-scale MHC class II genotyping of a wild lemur population by next generation sequencing

The critical role of major histocompatibility complex (MHC) genes in disease resistance, along with their putative function in sexual selection, reproduction and chemical ecology, make them an important genetic system in evolutionary ecology. Studying selective pressures acting on MHC genes in the wild nevertheless requires population-wide genotyping, which has long been challenging because of their extensive polymorphism. Here, we report on large-scale genotyping of the MHC class II loci of the grey mouse lemur (Microcebus murinus) from a wild population in western Madagascar. The second exons from MHC-DRB and -DQB of 772 and 672 individuals were sequenced, respectively, using a 454 sequencing platform, generating more than 800,000 reads. Sequence analysis, through a stepwise variant validation procedure, allowed reliable typing of more than 600 individuals. The quality of our genotyping was evaluated through three independent methods, namely genotyping the same individuals by both cloning and 454 sequencing, running duplicates, and comparing parent-offspring dyads; each displaying very high accuracy. A total of 61 (including 20 new) and 60 (including 53 new) alleles were detected at DRB and DQB genes, respectively. Both loci were non-duplicated, in tight linkage disequilibrium and in Hardy-Weinberg equilibrium, despite the fact that sequence analysis revealed clear evidence of historical selection. Our results highlight the potential of 454 sequencing technology in attempts to investigate patterns of selection shaping MHC variation in contemporary populations. The power of this approach will nevertheless be conditional upon strict quality control of the genotyping data.

Multiallelic models of genetic effects and variance decomposition in non-equilibrium populations

Quantitative genetics stems from the theoretical models of genetic effects, which are re-parameterizations of the genotypic values into parameters of biological (genetic) relevance. Different formulations of genetic effects are adequate to address different subjects. We thus need to generalize and unify them under a common framework for enabling researchers to easily transform genetic effects between different biological meanings. The Natural and Orthogonal Interactions (NOIA) model of genetic effects has been developed to achieve this aim. Here, we further implement the statistical formulation of NOIA with multiple alleles under Hardy-Weinberg departures (HWD). We show that our developments are straightforwardly connected to the decomposition of the genetic variance and we point out several emergent properties of multiallelic quantitative genetic models, as compared to the biallelic ones. Further, NOIA entails a natural extension of one-locus developments to multiple epistatic loci under linkage equilibrium. Therefore, we present an extension of the orthogonal decomposition of the genetic variance to multiple epistatic, multiallelic loci under HWD. We illustrate this theory with a graphical interpretation and an analysis of published data on the human acid phosphatase (ACP1) polymorphism.

Minireview: the SRC family of coactivators: an entrée to understanding a subset of polygenic diseases?

In this perspective, we present the idea that SRC family coactivators are likely agents in human polygenic disease states based upon a number of interlocking aspects of their biology. We argue that their role as key integrators of environmental signals and their ability to regulate the expression of myriad downstream genes makes them likely candidates for strong positive evolutionary selection pressures. Based on the fact that they work as part of multiprotein coactivator complexes, we predict that individual coactivator alleles exist as weakly penetrant disease alleles, each contributing only a fraction of transcriptional activity to the whole coactivator complex. In this way, individual coactivator alleles are free to evolve in the absence of strong negative selection. Emerging genomic and proteomic approaches promise to advance the characterization of coactivator proteins and their physiological functions, allowing us to have a greater appreciation of their roles as master regulators at the nexus between genetics, reproduction, metabolism, cancer, other human diseases, and our environment.

Detecting selection-induced departures from Hardy-Weinberg proportions

Viability selection influences the genotypic contexts of alleles and leads to quantifiable departures from Hardy-Weinberg proportions. One measure of these departures is Wright's inbreeding coefficient (F), where observed heterozygosity is compared with expected heterozygosity. Here, I extend population genetics theory to describe post-selection genotype frequencies in terms of post-selection allele frequencies and fitness dominance. The resulting equations correspond to non-equilibrium populations, allowing the following questions to be addressed: When selection is present, how large a sample size is needed to detect significant departures from Hardy-Weinberg? How do selection-induced departures from Hardy-Weinberg vary with allele frequencies and levels of fitness dominance? For realistic selection coefficients, large sample sizes are required and departures from Hardy-Weinberg proportions are small.

A fundamental relationship between genotype frequencies and fitnesses

The set of possible postselection genotype frequencies in an infinite, randomly mating population is found. Geometric mean heterozygote frequency divided by geometric mean homozygote frequency equals two times the geometric mean heterozygote fitness divided by geometric mean homozygote fitness. The ratio of genotype frequencies provides a measure of genetic variation that is independent of allele frequencies. When this ratio does not equal two, either selection or population structure is present. Within-population HapMap data show population-specific patterns, while pooled data show an excess of homozygotes.