Selection, Hitchhiking and Disequilibrium Analysis at Three Linked Loci with Application to Hla Data

The HLA system has been extensively studied from an evolutionary perspective. Although it is clear that selection has acted on the genes in the HLA complex, the nature of this selection has yet to be fully clarified. A study of constrained disequilibrium values is presented that is applicable to HLA and other less polymorphic systems with three or more linked loci, with the purpose of identifying selection events. The method uses the fact that three locus systems impose additional constraints on the range of possible disequilibrium values for any pair of loci. We have thus examined the behavior of the normalized pairwise disequilibrium measures using two locus (D'), and also three locus (D"), constraints on pairwise disequilibria in a three locus system when one of the three loci is under positive selection. The difference between these measures, delta = magnitude of D' - magnitude of D", has a distribution for the two unselected loci differing from that for the selected locus with either of the unselected loci (the hallmark is a high positive value of delta for the two unselected loci). An examination of genetic drift indicates that positive delta values are unlikely to be found in human populations in the absence of selection when recombination is greater than about 0.1%. This measure can thus provide insight into which allele of several linked loci might have been subject to selection. Application of this method to HLA haplotypes from a large French population study (Provinces Francaise) identifies selected alleles on particular haplotypes. Application of a complementary method, disequilibrium pattern analysis also confirms the action of selection on these haplotypes.     

High allelic polymorphism,moderate sequence diversity and diversifying selection for B-NK but not B-lec,the pair of lectin-like receptor genes in the chicken MHC

We previously characterised the C-type lectin-like receptor genes B-NK and B-lec, located next to each other in opposite orientations in the chicken major histocompatibility complex (MHC). We showed that B-NK is an inhibitory receptor expressed on natural killer cells, whereas B-lec is an activation-induced receptor with a broader expression pattern. It is interesting to note that the chicken MHC has been linked with resistance or susceptibility to Marek's disease virus (MDV), an oncogenic herpes virus. Recent reports show that the C-type lectin-like receptors in mouse and rat (Ly49H, NKR-P1 and Clr) are associated with resistance to another herpesvirus, cytomegalovirus (CMV). Therefore, B-NK and B-lec are potential candidate genes for the MHC-mediated resistance to MDV. In this paper, we report that both genes encode glycosylated type II membrane proteins that form disulphide-linked homodimers. The gene sequences from nine lines of domestic chicken representing seven haplotypes show that B-lec is well conserved between the different haplotypes, apparently under purifying selection. In contrast, B-NK has high allelic polymorphism and moderate sequence diversity, with 21 nucleotide changes in the complementary deoxyribonucleic acids (cDNAs) resulting in 20 amino acid substitutions. The allelic variations include substitutions, an indel and loss/gain of three predicted N-linked glycosylation sites. Strikingly, there is as much as 7% divergence between protein sequences of B-NK from different haplotypes, greater than the difference observed between the highly polymorphic human KIR NK receptors. Analysis of ds and dn reveal evidence of strong positive selection for B-NK to be polymorphic at the protein level, and modelling demonstrates significant variation between haplotypes in the predicted ligand binding face of B-NK.     

Population history and natural selection shape patterns of genetic variation in 132 genes

Identifying regions of the human genome that have been targets of natural selection will provide important insights into human evolutionary history and may facilitate the identification of complex disease genes. Although the signature that natural selection imparts on DNA sequence variation is difficult to disentangle from the effects of neutral processes such as population demographic history, selective and demographic forces can be distinguished by analyzing multiple loci dispersed throughout the genome. We studied the molecular evolution of 132 genes by comprehensively resequencing them in 24 African-Americans and 23 European-Americans. We developed a rigorous computational approach for taking into account multiple hypothesis tests and demographic history and found that while many apparent selective events can instead be explained by demography, there is also strong evidence for positive or balancing selection at eight genes in the European-American population, but none in the African-American population. Our results suggest that the migration of modern humans out of Africa into new environments was accompanied by genetic adaptations to emergent selective forces. In addition, a region containing four contiguous genes on Chromosome 7 showed striking evidence of a recent selective sweep in European-Americans. More generally, our results have important implications for mapping genes underlying complex human diseases.     

Characterization of a QTL region affecting clinical mastitis and protein yield on BTA6

Quantitative trait loci affecting clinical mastitis were detected and fine mapped to a narrow region on bovine chromosome 6 in the Norwegian Red cattle population. The region includes the casein gene cluster and several candidate genes thought to influence clinical mastitis. The most significant results were found for SNPs within the Mucin 7 gene. This gene encodes an antimicrobial peptide and constitutes part of the first line of defence for the mucosal immune system. Detection of long haplotypes extending several Mb may indicate that artificial selection has influenced the haplotype structures in the region. A search for selection sweeps supports this observation and coincides with association results found both by single SNP and haplotype analyses. Our analyses identified haplotypes carrying quantitative trait loci alleles associated with high protein yield and simultaneously fewer incidences of clinical mastitis. The fact that such haplotypes are found in relative high frequencies in Norwegian Red may reflect the combined breeding goal that is characterized by selection for both milk production and disease resistance. The identification of these haplotypes raises the possibility of overcoming the unfavourable genetic correlation between these traits through haplotype-assisted selection.     

A hierarchical Bayesian model for next-generation population genomics

The demography of populations and natural selection shape genetic variation across the genome and understanding the genomic consequences of these evolutionary processes is a fundamental aim of population genetics. We have developed a hierarchical Bayesian model to quantify genome-wide population structure and identify candidate genetic regions affected by selection. This model improves on existing methods by accounting for stochastic sampling of sequences inherent in next-generation sequencing (with pooled or indexed individual samples) and by incorporating genetic distances among haplotypes in measures of genetic differentiation. Using simulations we demonstrate that this model has a low false-positive rate for classifying neutral genetic regions as selected genes (i.e., Φ(ST) outliers), but can detect recent selective sweeps, particularly when genetic regions in multiple populations are affected by selection. Nonetheless, selection affecting just a single population was difficult to detect and resulted in a high false-negative rate under certain conditions. We applied the Bayesian model to two large sets of human population genetic data. We found evidence of widespread positive and balancing selection among worldwide human populations, including many genetic regions previously thought to be under selection. Additionally, we identified novel candidate genes for selection, several of which have been linked to human diseases. This model will facilitate the population genetic analysis of a wide range of organisms on the basis of next-generation sequence data.     

Polymorphism,haplotype composition,and selection in the <Emphasis Type="Italic">Mhc-DRB</Emphasis> of wild baboons

General patterns of organization in the major histocompatibility complex (MHC) have been successfully explained by the model of birth-and-death evolution, but understanding why certain MHC genes are maintained together into specific haplotypes remains challenging. The haplotype configurations of the functionally important class II DR region have been described in few primates and display important interspecific variability with respect to the extent of allelic variation, the number of loci and/or combinations of loci present. Understanding the evolutionary mechanisms driving such variation is conditional upon characterizing haplotypes in new species and identifying the selective pressures acting on haplotypes. This study explores the variability of haplotype configurations in the Mhc-DRB region (exon 2) for the first time in wild non-human primates, chacma baboons (Papio ursinus). Paur-DRB haplotypes were characterized through segregation studies and linkage disequilibrium. 23 Paur-DRB sequences and 15 haplotype configurations were identified in 199 animals. The Paur-DRB exon 2 is shown to be subjected to intense positive selection and frequent recombination. An approach recently developed for human vaccine studies was used to classify Paur-DRB sequences into supertypes, based on the physico-chemical properties of amino acids that are positively selected, thus most probably involved in antigen recognition. Sequences grouped into the same supertype (thus presumably sharing antigen-binding affinities) are non-randomly distributed within haplotypes, leading to an increased individual diversity of supertypes. Our results suggest that selection favoring haplotypes with complementary sets of DRB supertypes shapes functionally tuned haplotypes in this natural baboon population.     

Molecular evolution of transferrin: evidence for positive selection in salmonids

Transferrins are iron-binding proteins that are involved in iron storage and resistance to bacterial disease. Previous work has shown that nonsynonymous-to-synonymous-site substitution ratios (d(n)/d(s) ratios) between transferrin genes from some salmonid species were significantly greater than 1.0, providing evidence for positive selection at the transferrin gene. The purpose of the current study was to put these earlier results in a broader evolutionary context by examining variation among 25 previously published transferrin sequences from fish, amphibians, and mammals. The results of the study show that evidence for positive selection at transferrin is limited to salmonids-d(n)/d(s) ratios estimated for nonsalmonid lineages were generally less than 1.0. Within the salmonids, approximately 13% of the transferrin codon sites are estimated to be subject to positive selection, with an estimated d(n)/d(s) ratio of approximately 7. The three- dimensional locations of some of the selected sites were inferred by comparing these sites to homologous sites in the bovine lactoferrin crystallographic structure. The selected sites generally fall on the outside of the molecule, within and near areas that are bound by transferrin-binding proteins from human pathogenic bacteria. The physical locations of sites estimated to be subject to positive selection support previous speculation that competition for iron from pathogenic bacteria could be the source of positive selection.     

Patterns of Population Differentiation and Natural Selection on the Celiac Disease Background Risk Network

Celiac disease is a common small intestinal inflammatory condition induced by wheat gluten and related proteins from rye and barley. Left untreated, the clinical presentation of CD can include failure to thrive, malnutrition, and distension in juveniles. The disease can additionally lead to vitamin deficiencies, anemia, and osteoporosis. Therefore, CD potentially negatively affected fitness in past populations utilizing wheat, barley, and rye. Previous analyses of CD risk variants have uncovered evidence for positive selection on some of these loci. These studies also suggest the possibility that risk for common autoimmune conditions such as CD may be the result of positive selection on immune related loci in the genome to fight infection. Under this evolutionary scenario, disease phenotypes may be a trade-off from positive selection on immunity. If this hypothesis is generally true, we can expect to find a signal of natural selection when we survey across the network of loci known to influence CD risk. This study examines the non-HLA autosomal network of gene loci associated with CD risk in Europe. We reject the null hypothesis of neutrality on this network of CD risk loci. Additionally, we can localize evidence of selection in time and space by adding information from the genome of the Tyrolean Iceman. While we can show significant differentiation between continental regions across the CD network, the pattern of evidence is not consistent with primarily recent (Holocene) selection across this network in Europe. Further localization of ancient selection on this network may illuminate the ecological pressures acting on the immune system during this critically interesting phase of our evolution.     

Network-Based Multiple Sclerosis Pathway Analysis with GWAS Data from 15,000 Cases and 30,000 Controls

Multiple sclerosis (MS) is an inflammatory CNS disease with a substantial genetic component, originally mapped to only the human leukocyte antigen (HLA) region. In the last 5 years, a total of seven genome-wide association studies and one meta-analysis successfully identified 57 non-HLA susceptibility loci. Here, we merged nominal statistical evidence of association and physical evidence of interaction to conduct a protein-interaction-network-based pathway analysis (PINBPA) on two large genetic MS studies comprising a total of 15,317 cases and 29,529 controls. The distribution of nominally significant loci at the gene level matched the patterns of extended linkage disequilibrium in regions of interest. We found that products of genome-wide significantly associated genes are more likely to interact physically and belong to the same or related pathways. We next searched for subnetworks (modules) of genes (and their encoded proteins) enriched with nominally associated loci within each study and identified those modules in common between the two studies. We demonstrate that these modules are more likely to contain genes with bona fide susceptibility variants and, in addition, identify several high-confidence candidates (including BCL10, CD48, REL, TRAF3, and TEC). PINBPA is a powerful approach to gaining further insights into the biology of associated genes and to prioritizing candidates for subsequent genetic studies of complex traits.     

Selecting SNPs for association studies based on population frequencies: a novel interactive tool and its application to polygenic diseases

Common complex polygenic diseases as autoimmune diseases have not been completely understood on a molecular level. While many genes are known to be involved in the pathways responsible for the phenotype, explicit causes for the susceptibility of the disease remain to be elucidated. The susceptibility to disease is thought to be the result of genetic epistatic interactions between common polymorphic genes. This polymorphism is mostly caused by single nucleotide polymorphisms (SNPs). Human subpopulations are known to differ in the susceptibility to the diseases and generally in the distribution of single nucleotide polymorphisms. The here presented approach retrieves SNPs with the most divergent frequencies for selected human subpopulations to help defining properties for the experimental verification of SNPs within defined regions. A web-accessible program implementing this approach was evaluated for multiple sclerosis (MS), a common human polygenic disease. A link to a summary of data from "The SNP Consortium" (TSC) with sex-dependencies of SNPs is available. Associations of SNPs to genes, genetic markers and chromosomal loci are retrieved from the Ensembl project. This tool is recommended to be used in conjunction with microarray analyses or marker association studies that link genes or chromosomal loci to particular diseases.     

Reconstruction of a functional human gene network, with an application for prioritizing positional candidate genes

Most common genetic disorders have a complex inheritance and may result from variants in many genes, each contributing only weak effects to the disease. Pinpointing these disease genes within the myriad of susceptibility loci identified in linkage studies is difficult because these loci may contain hundreds of genes. However, in any disorder, most of the disease genes will be involved in only a few different molecular pathways. If we know something about the relationships between the genes, we can assess whether some genes (which may reside in different loci) functionally interact with each other, indicating a joint basis for the disease etiology. There are various repositories of information on pathway relationships. To consolidate this information, we developed a functional human gene network that integrates information on genes and the functional relationships between genes, based on data from the Kyoto Encyclopedia of Genes and Genomes, the Biomolecular Interaction Network Database, Reactome, the Human Protein Reference Database, the Gene Ontology database, predicted protein-protein interactions, human yeast two-hybrid interactions, and microarray co-expressions. We applied this network to interrelate positional candidate genes from different disease loci and then tested 96 heritable disorders for which the Online Mendelian Inheritance in Man database reported at least three disease genes. Artificial susceptibility loci, each containing 100 genes, were constructed around each disease gene, and we used the network to rank these genes on the basis of their functional interactions. By following up the top five genes per artificial locus, we were able to detect at least one known disease gene in 54% of the loci studied, representing a 2.8-fold increase over random selection. This suggests that our method can significantly reduce the cost and effort of pinpointing true disease genes in analyses of disorders for which numerous loci have been reported but for which most of the genes are unknown.     

Genetic polymorphism and natural selection in the malaria parasite Plasmodium falciparum.

We have studied the genetic polymorphism at 10 Plasmodium falciparum loci that are considered potential targets for specific antimalarial vaccines. The polymorphism is unevenly distributed among the loci; loci encoding proteins expressed on the surface of the sporozoite or the merozoite (AMA-1, CSP, LSA-1, MSP-1, MSP-2, and MSP-3) are more polymorphic than those expressed during the sexual stages or inside the parasite (EBA-175, Pfs25, PF48/45, and RAP-1). Comparison of synonymous and nonsynonymous substitutions indicates that natural selection may account for the polymorphism observed at seven of the 10 loci studied. This inference depends on the assumption that synonymous substitutions are neutral, which we test by analyzing codon bias and G+C content in a set of 92 gene loci. We find evidence for an overall trend towards increasing A+T richness, but no evidence for mutation bias. Although the neutrality of synonymous substitutions is not definitely established, this trend towards an A+T rich genome cannot explain the accumulation of substitutions at least in the case of four genes (AMA-1, CSP, LSA-1, and PF48/45) because the Gleft and right arrow C transversions are more frequent than expected. Moreover, the Tajima test manifests positive natural selection for the MSP-1 and, less strongly, MSP-3 polymorphisms; the McDonald-Kreitman test manifests natural selection at LSA-1 and PF48/45. We conclude that there is definite evidence for positive natural selection in the genes encoding AMA-1, CSP, LSA-1, MSP-1, and Pfs48/45. For four other loci, EBA-175, MSP-2, MSP-3, and RAP-1, the evidence is limited. No evidence for natural selection is found for Pfs25.     

Can Genetic Analysis of Putative Blood Alzheimer’s Disease Biomarkers Lead to Identification of Susceptibility Loci?

Although 24 Alzheimer’s disease (AD) risk loci have been reliably identified, a large portion of the predicted heritability for AD remains unexplained. It is expected that additional loci of small effect will be identified with an increased sample size. However, the cost of a significant increase in Case-Control sample size is prohibitive. The current study tests whether exploring the genetic basis of endophenotypes, in this case based on putative blood biomarkers for AD, can accelerate the identification of susceptibility loci using modest sample sizes. Each endophenotype was used as the outcome variable in an independent GWAS. Endophenotypes were based on circulating concentrations of proteins that contributed significantly to a published blood-based predictive algorithm for AD. Endophenotypes included Monocyte Chemoattractant Protein 1 (MCP1), Vascular Cell Adhesion Molecule 1 (VCAM1), Pancreatic Polypeptide (PP), Beta2 Microglobulin (B2M), Factor VII (F7), Adiponectin (ADN) and Tenascin C (TN-C). Across the seven endophenotypes, 47 SNPs were associated with outcome with a p-value ≤1x10^-7. Each signal was further characterized with respect to known genetic loci associated with AD. Signals for several endophenotypes were observed in the vicinity of CR1, MS4A6A/MS4A4E, PICALM, CLU, and PTK2B. The strongest signal was observed in association with Factor VII levels and was located within the F7 gene. Additional signals were observed in MAP3K13, ZNF320, ATP9B and TREM1. Conditional regression analyses suggested that the SNPs contributed to variation in protein concentration independent of AD status. The identification of two putatively novel AD loci (in the Factor VII and ATP9B genes), which have not been located in previous studies despite massive sample sizes, highlights the benefits of an endophenotypic approach for resolving the genetic basis for complex diseases. The coincidence of several of the endophenotypic signals with known AD loci may point to novel genetic interactions and should be further investigated.     

Genome Patterns of Selection and Introgression of Haplotypes in Natural Populations of the House Mouse (Mus musculus)

General parameters of selection, such as the frequency and strength of positive selection in natural populations or the role of introgression, are still insufficiently understood. The house mouse (Mus musculus) is a particularly well-suited model system to approach such questions, since it has a defined history of splits into subspecies and populations and since extensive genome information is available. We have used high-density single-nucleotide polymorphism (SNP) typing arrays to assess genomic patterns of positive selection and introgression of alleles in two natural populations of each of the subspecies M. m. domesticus and M. m. musculus. Applying different statistical procedures, we find a large number of regions subject to apparent selective sweeps, indicating frequent positive selection on rare alleles or novel mutations. Genes in the regions include well-studied imprinted loci (e.g. Plagl1/Zac1), homologues of human genes involved in adaptations (e.g. alpha-amylase genes) or in genetic diseases (e.g. Huntingtin and Parkin). Haplotype matching between the two subspecies reveals a large number of haplotypes that show patterns of introgression from specific populations of the respective other subspecies, with at least 10% of the genome being affected by partial or full introgression. Using neutral simulations for comparison, we find that the size and the fraction of introgressed haplotypes are not compatible with a pure migration or incomplete lineage sorting model. Hence, it appears that introgressed haplotypes can rise in frequency due to positive selection and thus can contribute to the adaptive genomic landscape of natural populations. Our data support the notion that natural genomes are subject to complex adaptive processes, including the introgression of haplotypes from other differentiated populations or species at a larger scale than previously assumed for animals. This implies that some of the admixture found in inbred strains of mice may also have a natural origin.     

Pervasive adaptive evolution in primate seminal proteins

Seminal fluid proteins show striking effects on reproduction, involving manipulation of female behavior and physiology, mechanisms of sperm competition, and pathogen defense. Strong adaptive pressures are expected for such manifestations of sexual selection and host defense, but the extent of positive selection in seminal fluid proteins from divergent taxa is unknown. We identified adaptive evolution in primate seminal proteins using genomic resources in a tissue-specific study. We found extensive signatures of positive selection when comparing 161 human seminal fluid proteins and 2,858 prostate-expressed genes to those in chimpanzee. Seven of eight outstanding genes yielded statistically significant evidence of positive selection when analyzed in divergent primates. Functional clues were gained through divergent analysis, including several cases of species-specific loss of function in copulatory plug genes, and statistically significant spatial clustering of positively selected sites near the active site of kallikrein 2. This study reveals previously unidentified positive selection in seven primate seminal proteins, and when considered with findings in Drosophila, indicates that extensive positive selection is found in seminal fluid across divergent taxonomic groups.     

Genetic susceptibility to Alzheimer disease

Alzheimer disease (AD) is the leading cause of dementia in the elderly. Less than a decade ago, it was questioned as to whether or not genes were even involved in anything but rare early onset AD. Since that time, using a variety of genetic epidemiological and molecular biological techniques, four loci have been identified that play a role in the genetic susceptibility of AD, AD presents as a prototype of the power of genetic techniques in defining the etiology of a complex disease.     

<Emphasis Type="Italic">KIR</Emphasis> haplotypes defined by segregation analysis in 59 Centre d’Etude Polymorphisme Humain (CEPH) families

The killer cell immunoglobulin-like receptor (KIR) gene cluster exhibits extensive allelic and haplotypic diversity. Variation at the locus is associated with an increasing number of human diseases, reminiscent of the HLA loci. Characterization of diversity at the KIR locus has progressed over the past several years, particularly since the sequence of entire KIR haplotypes have become available. To determine the extent of KIR haplotypic variability among individuals of northern European descent, we genotyped 59 CEPH families for presence/absence of all KIR genes and performed limited allelic subtyping at several KIR loci. A total of 20 unique haplotypes differing in gene content were identified, the most common of which was the previously defined A haplotype (f = 0.52). Several unusual haplotypes that probably arose as a consequence of unequal crossing over events were also identified. Linkage disequilibrium (LD) analysis indicated strong negative and positive LD between several pairs of genes, values that may be useful in determining haplotypic structure when family data are not available. These data provide a resource to aid in the interpretation of disease association data involving individuals of European descent. An erratum to this article can be found at     

Single Nucleotide Polymorphisms in the Wnt and BMP Pathways and Colorectal Cancer Risk in a Spanish Cohort

Background

Colorectal cancer (CRC) is considered a complex disease, and thus the majority of the genetic susceptibility is thought to lie in the form of low-penetrance variants following a polygenic model of inheritance. Candidate-gene studies have so far been one of the basic approaches taken to identify these susceptibility variants. The consistent involvement of some signaling routes in carcinogenesis provided support for pathway-based studies as a natural strategy to select genes that could potentially harbour new susceptibility loci.

Methodology/Principal Findings

We selected two main carcinogenesis-related pathways: Wnt and BMP, in order to screen the implicated genes for new risk variants. We then conducted a case-control association study in 933 CRC cases and 969 controls based on coding and regulatory SNPs. We also included rs4444235 and rs9929218, which did not fulfill our selection criteria but belonged to two genes in the BMP pathway and had consistently been linked to CRC in previous studies. Neither allelic, nor genotypic or haplotypic analyses showed any signs of association between the 37 screened variants and CRC risk. Adjustments for sex and age, and stratified analysis between sporadic and control groups did not yield any positive results either.

Conclusions/Significance

Despite the relevance of both pathways in the pathogenesis of the disease, and the fact that this is indeed the first study that considers these pathways as a candidate-gene selection approach, our study does not present any evidence of the presence of low-penetrance variants for the selected markers in any of the considered genes in our cohort.