MHC-dependent mate choice in humans: why genomic patterns from the HapMap European American dataset support the hypothesis

The role of the major histocompatibility complex (MHC) in mate choice in humans is controversial. Nowadays, the availability of genetic variation data at genomic scales allows for a careful assessment of this question. In 2008, Chaix et al. reported evidence for MHC-dependent mate choice among European American spouses from the HapMap 2 dataset. Recently, Derti et al. suggested that this observation was not robust. Furthermore, when Derti et al. applied similar analyses to the HapMap 3 European American samples, they did not see a significant effect. Although some of the points raised by Derti et al. are relevant, we disagree with the reported absence of evidence for MHC-dependent mate choice within the HapMap samples. More precisely, we show here that the MHC dissimilarity among HapMap 3 European American spouses is still extreme in comparison to the rest of the genome, even after multiple testing correction. This finding supports the hypothesis of MHC-dependent mate choice in some human populations.     

Navigating the HapMap

With the availability of the HapMap--a resource which describes common patterns of linkage disequilibrium (LD) in four different human population samples, we now have a powerful tool to help dissect the role of genetic variation in the biology of the genome. HapMap is entirely complimentary to the human genome map and so it is particularly fitting that it should be viewed in a full genomic context. However, characterization of high resolution LD across the genome can be a challenging task, owing in part to the sheer volume of data and the inherent dimensionality that its analysis entails. However, a number of tools are now available to make this task easier for researchers. This review will examine tools for viewing and analysing haplotype and LD data, enabling a number of tasks; including identification of optimal sets of haplotype tagging single nucleotide polymorphisms (SNPs); drawing links between associated SNPs and putative causal alleles; or simply viewing LD and haplotypes across a gene or region of interest. The data generated by the HapMap also has other important applications, informing, for example, on the demographic history and evidence of selection in human populations and on previously undetected regulatory relationships and gene networks. All of these properties make the HapMap no less an important resource than the human genome sequence itself and so this makes it essential viewing for all in the field of human biology.     

HapMap filter 1.0: a tool to preprocess the HapMap genotypic data for association studies

The International HapMap Project provides a resource of genotypic data on single nucleotide polymorphisms (SNPs), which can be used in various association studies to identify the genetic determinants for phenotypic variations. Prior to the association studies, the HapMap dataset should be preprocessed in order to reduce the computation time and control the multiple testing problem. The less informative SNPs including those with very low genotyping rate and SNPs with rare minor allele frequencies to some extent in one or more population are removed. Some research designs only use SNPs in a subset of HapMap cell lines. Although the HapMap website and other association software packages have provided some basic tools for optimizing these datasets, a fast and user-friendly program to generate the output for filtered genotypic data would be beneficial for association studies. Here, we present a flexible, straight-forward bioinformatics program that can be useful in preparing the HapMap genotypic data for association studies by specifying cell lines and two common filtering criteria: minor allele frequencies and genotyping rate. The software was developed for Microsoft Windows and written in C++. AVAILABILITY: The Windows executable and source code in Microsoft Visual C++ are available at Google Code (http://hapmap-filter-v1.googlecode.com/) or upon request. Their distribution is subject to GNU General Public License v3.     

Leveraging the HapMap correlation structure in association studies

Recent high-throughput genotyping technologies, such as the Affymetrix 500k array and the Illumina HumanHap 550 beadchip, have driven down the costs of association studies and have enabled the measurement of single-nucleotide polymorphism (SNP) allele frequency differences between case and control populations on a genomewide scale. A key aspect in the efficiency of association studies is the notion of "indirect association," where only a subset of SNPs are collected to serve as proxies for the uncollected SNPs, taking advantage of the correlation structure between SNPs. Recently, a new class of methods for indirect association, multimarker methods, has been proposed. Although the multimarker methods are a considerable advancement, current methods do not fully take advantage of the correlation structure between SNPs and their multimarker proxies. In this article, we propose a novel multimarker indirect-association method, WHAP, that is based on a weighted sum of the haplotype frequency differences. In contrast to traditional indirect-association methods, we show analytically that there is a considerable gain in power achieved by our method compared with both single-marker and multimarker tests, as well as traditional haplotype-based tests. Our results are supported by empirical evaluation across the HapMap reference panel data sets, and a software implementation for the Affymetrix 500k and Illumina HumanHap 550 chips is available for download.     

On the challenges of the HapMap resource

The International HapMap Project provides a key resource of genotypic data on human lymphoblastoid cell lines derived from four major world populations of European, African, Chinese and Japanese ancestry for researchers to associate with various phenotypic data to find genes affecting health, disease and response to drugs. Recently, the HapMap resource has significantly benefited research areas such as gene expression variation studies. Besides some intrinsic limitations, there are a few challenges that should be considered in the next wave of research using this tremendous resource. We suggest that overcoming these challenges or considering the confounding variables in the interpretation of results can provide more insights into the current views of the human genome as well as complex traits such as drug response variation and susceptibility to common diseases.     

Mining the HapMap to dissect complex traits

A report on the Keystone Symposium 'Genome Sequence Variation and the Inherited Basis of Common Disease and Complex Traits', Big Sky, USA, 8-13 January 2006.     

Sequencing complex diseases With HapMap

Determining the patterns of DNA sequence variation in the human genome is a useful first step toward identifying the genetic basis of a common disease. A haplotype map (HapMap), aimed at describing these variation patterns across the entire genome, has been recently developed by the International HapMap Consortium. In this article, we present a novel statistical model for directly characterizing specific sequence variants that are responsible for disease risk based on the haplotype structure provided by HapMap. Our model is developed in the maximum-likelihood context, implemented with the EM algorithm. We perform simulation studies to investigate the statistical properties of this disease-sequencing model. A worked example from a human obesity study with 155 patients was used to validate this model. In this example, we found that patients carrying a haplotype constituted by allele Gly16 at codon 16 and allele Gln27 at codon 27 genotyped within the beta2AR candidate gene display significantly lower body mass index than patients carrying the other haplotypes. The implications and extensions of our model are discussed.     

Functional analysis of HapMap SNPs

Genome-wide association studies (GWAS) have successfully identified many genetic variants associated with complex diseases and traits. However, functional consequence of genetic variants studied in GWAS is not yet fully investigated, which would hinder the application of GWAS. We therefore performed a systematic functional analysis of HapMap SNPs, which have been most commonly used as the reference panel for GWAS. Our study highlights several characteristics of HapMap SNPs and identifies subsets of genetic variants with interesting functional implication. The results show that HapMap SNPs have good coverage within RefSeq genes, especially within known disease-related genes. On the other hand, only a small percentage of SNPs are non-synonymous SNPs while many SNPs are actually located at gene deserts. Moreover, many functionally important variants are not yet still interrogated. A redesigned SNP reference panel with additional functionally important variants would be useful to identify disease-causal variants in the future genome-wide studies.     

Perspectives on human genetic variation from the HapMap Project

The completion of the International HapMap Project marks the start of a new phase in human genetics. The aim of the project was to provide a resource that facilitates the design of efficient genome-wide association studies, through characterising patterns of genetic variation and linkage disequilibrium in a sample of 270 individuals across four geographical populations. In total, over one million SNPs have been typed across these genomes, providing an unprecedented view of human genetic diversity. In this review we focus on what the HapMap project has taught us about the structure of human genetic variation and the fundamental molecular and evolutionary processes that shape it.     

Unexpected Relationships and Inbreeding in HapMap Phase III Populations

Correct annotation of the genetic relationships between samples is essential for population genomic studies, which could be biased by errors or omissions. To this end, we used identity-by-state (IBS) and identity-by-descent (IBD) methods to assess genetic relatedness of individuals within HapMap phase III data. We analyzed data from 1,397 individuals across 11 ethnic populations. Our results support previous studies (Pemberton et al., 2010; Kyriazopoulou-Panagiotopoulou et al., 2011) assessing unknown relatedness present within this population. Additionally, we present evidence for 1,657 novel pairwise relationships across 9 populations. Surprisingly, significant Cotterman''s coefficients of relatedness K1 (IBD1) values were detected between pairs of known parents. Furthermore, significant K2 (IBD2) values were detected in 32 previously annotated parent-child relationships. Consistent with a hypothesis of inbreeding, regions of homozygosity (ROH) were identified in the offspring of related parents, of which a subset overlapped those reported in previous studies (Gibson et al. 2010; Johnson et al. 2011). In total, we inferred 28 inbred individuals with ROH that overlapped areas of relatedness between the parents and/or IBD2 sharing at a different genomic locus between a child and a parent. Finally, 8 previously annotated parent-child relationships had unexpected K0 (IBD0) values (resulting from a chromosomal abnormality or genotype error), and 10 previously annotated second-degree relationships along with 38 other novel pairwise relationships had unexpected IBD2 (indicating two separate paths of recent ancestry). These newly described types of relatedness may impact the outcome of previous studies and should inform the design of future studies relying on the HapMap Phase III resource.     

Evaluating coverage of exons by HapMap SNPs

Genome-wide association (GWA) studies are currently one of the most powerful tools in identifying disease-associated genes or variants. In typical GWA studies, single-nucleotide polymorphisms (SNPs) are often used as genetic makers. Therefore, it is critical to estimate the percentage of genetic variations which can be covered by SNPs through linkage disequilibrium (LD). In this study, we use the concept of haplotype blocks to evaluate the coverage of five SNP sets including the HapMap and four commercial arrays, for every exon in the human genome. We show that although some Chips can reach similar coverage as the HapMap, only about 50% of exons are completely covered by haplotype blocks of HapMap SNPs. We suggest further high-resolution genotyping methods are required, to provide adequate genome-wide power for identifying variants.     

Risky mate search and mate preference in the golden orb-web spider (Nephila plumipes)

Mate searching is a risky behavior that decreases survival byincreasing predation risk and the risk of energy depletion.However, few studies have quantified actual mortality duringmate search, making it difficult to predict mate searching andmating strategies. Using a mark and recapture study, we examinedmate-searching success in a highly sexually dimorphic species,the golden orb-web spider (Nephila plumipes). We show that despitethe high-density aggregations of this species, male survivalduring mate searching is extremely low (36%) and is phenotypeindependent. Surprisingly, males that survived mate search werein better condition after recapture than prior to release, mostlikely due to kleptoparasitism on females' webs. In a complementaryrelease experiment in a field enclosure, we show that malesare choosy and adjust their choice of female depending on theirown condition and weight. Thus, the high mortality rate of searchingmales in the field may be a cost of choosiness because releasedmales traveled further than necessary to settle on females.Although males were choosy about female phenotypes, they didnot avoid webs with rival males already present. This suggeststhat the cost of continued searching outweighs the cost of competitionbut not the cost of mating with certain females. Further examinationsof mate-searching risk in other species in reference to theirmating system and environmental conditions are necessary todetermine the occurrence and effects of high mortality ratesduring searching.     

Population differences in the rate of proliferation of international HapMap cell lines

The International HapMap Project is a resource for researchers containing genotype, sequencing, and expression information for EBV-transformed lymphoblastoid cell lines derived from populations across the world. The expansion of the HapMap beyond the four initial populations of Phase 2, referred to as Phase 3, has increased the sample number and ethnic diversity available for investigation. However, differences in the rate of cellular proliferation between the populations can serve as confounders in phenotype-genotype studies using these cell lines. Within the Phase 2 populations, the JPT and CHB cell lines grow faster (p < 0.0001) than the CEU or YRI cell lines. Phase 3 YRI cell lines grow significantly slower than Phase 2 YRI lines (p < 0.0001), with no widespread genetic differences based on common SNPs. In addition, we found significant growth differences between the cell lines in the Phase 2 ASN populations and the Han Chinese from the Denver metropolitan area panel in Phase 3 (p < 0.0001). Therefore, studies that separate HapMap panels into discovery and replication sets must take this into consideration.     

Evidence that sex‐specific signals may support mate finding and limit aggregation in the dung beetle Aphodius fossor

1. In temperate climates, dung is often colonised by several species of endocoprid (dwelling) dung beetles which use pats for feeding, shelter, and reproduction. 2. Endocoprid beetles aggregate even when offered patches (dung pats) of consistent age, size, and origin, suggesting that beetles themselves might influence the attractiveness of patches to members of their own species. Both pheromones, and physical changes to the structure of dung pats caused by colonising beetles have been suggested as mechanisms facilitating intraspecific aggregation, but neither of these hypotheses have been empirically tested. 3. Using a common European dung beetle (Aphodius fossor L.), we conducted a simple choice experiment designed to test whether (i) earlier colonisation by conspecifics could alter dung attractiveness and (ii) whether attraction was influenced by sex‐specific signals. 4. We found that female beetles are repelled by dung colonised by conspecific females and are attracted to dung colonised by conspecific males. Male beetles show no evidence of attraction or repellence for dung colonised by either sex. Neither in females nor males was uncolonised dung found to be significantly more or less attractive than predicted by non‐preference. 5. Our results suggest that for A. fossor male‐produced signals may support mate finding in patchy environments, and that female‐produced signals may serve to discourage subsequent colonisation by additional females.     

Reproductive consequences of mate quantity versus mate diversity in a wind-pollinated plant

Since most pollen travels limited distances in wind-pollinated plants, both the local quantity and diversity of mates may limit female reproductive success. Yet little evidence exists on their relative contribution, despite the importance of viable seed production to population dynamics.To study how variation in female reproductive success is affected by the quantity versus the diversity of surrounding mates contributing pollen, we integrated pollination experiments, data on natural seed set and seed viability, and AFLP genetic marker data in the wind-pollinated dioecious clonal forest herb Mercurialis perennis.Pollination experiments indicated weak quantitative pollen limitation effects on seed set. Among-population crosses showed reduced seed viability, suggesting outbreeding depression due to genetic divergence. Pollination with pollen from a single source did not negatively affect reproductive success. These findings were consistent with results of the survey of natural female reproductive success. Seed set decreased with the distance to males in a female plants’ local neighborhood, suggesting a shortage of pollen in isolated female plants, and increased with the degree of local genetic diversity. Spatial isolation to other populations and population size did not affect seed set. None of these variables were related to seed viability.We conclude that pollen movement in M. perennis is likely very limited. Both male proximity and the local degree of genetic diversity influenced female reproductive success.