Estimates of natural selection due to protein tertiary structure inform the ancestry of biallelic loci

We consider the inference of which of two alleles is ancestral when the alleles have a single nonsynonymous difference and when natural selection acts via protein tertiary structure. Whereas the probability that an allele is ancestral under neutrality is equal to its frequency, under selection this probability depends on allele frequency and on the magnitude and direction of selection pressure. Although allele frequencies can be well estimated from intraspecific data, small fitness differences have a large evolutionary impact but can be difficult to estimate with only intraspecific data. Methods for predicting aspects of phenotype from genotype can supplement intraspecific sequence data. Recently developed statistical techniques can assess effects of phenotypes, such as protein tertiary structure on molecular evolution. While these techniques were initially designed for comparing protein-coding genes from different species, the resulting interspecific inferences can be assigned population genetic interpretations to assess the effect of selection pressure, and we use them here along with intraspecific allele frequency data to estimate the probability that an allele is ancestral. We focus on 140 nonsynonymous single nucleotide polymorphisms of humans that are in proteins with known tertiary structures. We find that our technique for employing protein tertiary structure information yields some biologically plausible results but that it does not substantially improve the inference of ancestral human allele types.     

Novel (CA)n marker DXYS241 on the nonrecombinant part of the human Y chromosome

The origin of modern humans can be traced by comparing polymorphic sites in either mitochondria or genomic sequences between humans and other primates. The human Y chromosome has both a non-recombining region and X-Y homologous pseudo-autosomal regions. In the nonrecombining region events during evolution can be directly detected. At least a part of homology between Xq21 and Yp11 is a result of rather recent translocations from the X chromosome to the Y chromosome. DNA markers residing in the nonrecombining region of the human Y chromosome are potentially useful in tracing male-specific gene flow in human evolution. However, the number of available markers in the region is limited. Here, we report a novel X-Y homologous (CA)n repeat locus in the nonrecombining region of the Y chromosome. This marker, DXYS241, has several interesting features. Y- and X-chromosome alleles are distinguishable because the Y-chromosome alleles are shorter than the X-chromosome alleles most of the time. We developed 2 primer sets for specific examination of Y- and X-chromosome alleles. The marker should be useful in establishing relationships between populations based on patrilineal gene flow. Sequences homologous to DXYS241 are also found on the X chromosome of primates. Four events during primate evolution that led to the modern human Y chromosome were identified.     

Signatures of selection in the human olfactory receptor OR5I1 gene

The human olfactory receptor (OR) repertoire is reduced in comparison to other mammals and to other nonhuman primates. Nonetheless, this olfactory decline opens an opportunity for evolutionary innovation and improvement. In the present study, we focus on an OR gene, OR5I1, which had previously been shown to present an excess of amino acid replacement substitutions between humans and chimpanzees. We analyze the genetic variation in OR5I1 in a large worldwide human panel and find an excess of derived alleles segregating at relatively high frequencies in all populations. Additional evidence for selection includes departures from neutrality in allele frequency spectra tests but no unusually extended haplotype structure. Moreover, molecular structural inference suggests that one of the nonsynonymous polymorphisms defining the presumably adaptive protein form of OR5I1 may alter the functional binding properties of the OR. These results are compatible with positive selection having modeled the pattern of variation found in the OR5I1 gene and with a relatively ancient, mild selective sweep predating the "Out of Africa" expansion of modern humans.     

Estimation of allele frequencies in polyploids under certain patterns of inheritance

Allele frequencies have long been studied by biologists interested in evolution and speciation. More recently, with the application of molecular markers in human DNA profiling we have also seen the need for reliable population allele frequency estimates for making probabilistic inferences. There is now interest in applying the same DNA profiling technology to identification of plant varieties. HortResearch maintains a large germplasm of horticultural plant species. It is becoming evident that accurate identification of these accessions through DNA fingerprinting is essential for effective utilisation and maintenance of this germplasm. Microsatellites are the markers of choice for this fingerprinting. However, such markers do not reveal the dosage of alleles in a polyploid. Polyploidy is common amongst horticultural plants. Estimating allele frequencies in a polyploid population is, therefore, complicated because of some marker genotypes being phenotypically indistinguishable. For example, in a tetraploid, with four alleles at a locus showing polysomic inheritance, although 35 genotypes are possible, these will fall into only 15 marker phenotypic classes. Furthermore 'null' individuals are rarely detected in polyploids. Furthermore, some polyploids can be cryptic exhibiting disomy, instead of the polysomic inheritance. We will discuss the implications of these factors and present an EM-type algorithm for estimating allele frequencies of a polyploid population under certain patterns of inheritance. The method will be demonstrated on simulated data. We also discuss the nature of some of the additional problems that may be encountered with estimating allele frequencies in polyploids for which other solutions still need to be developed.     

Allelic dimorphism in the human tissue-type plasminogen activator (TPA) gene as a result of an Alu insertion/deletion event

Summary Polymerase chain reaction and direct sequencing were used to investigate an amplified DNA fragment containing the suspected polymorphic site of all known intragenic restriction fragment length polymorphisms (RFLPs) within the human tissue-type plasminogen activator (TPA) gene. Sequence data obtained showed that these RFLPs were all generated by the presence or absence of one of the two Alu sequences located in intron h of the human TPA gene. Furthermore, one of the direct repeats flanking this Alu sequence was absent in the minor allele. In addition to indicating the presence of an Alu insertion in an ancestral human TPA gene, these findings suggest a slip-replication mechanism for the deletion of this Alu repeat, once inserted into the gene. As both alleles have been observed in similar frequencies among different ethnic groups, the insertion or subsequent deletion of this Alu sequence in the human TPA gene must have occurred early in human evolution.     

Context-dependent mutation rates may cause spurious signatures of a fixation bias favoring higher GC-content in humans

Understanding the proximate and ultimate causes underlying the evolution of nucleotide composition in mammalian genomes is of fundamental interest to the study of molecular evolution. Comparative genomics studies have revealed that many more substitutions occur from G and C nucleotides to A and T nucleotides than the reverse, suggesting that mammalian genomes are not at equilibrium for base composition. Analysis of human polymorphism data suggests that mutations that increase GC-content tend to be at much higher frequencies than those that decrease or preserve GC-content when the ancestral allele is inferred via parsimony using the chimpanzee genome. These observations have been interpreted as evidence for a fixation bias in favor of G and C alleles due to either positive natural selection or biased gene conversion. Here, we test the robustness of this interpretation to violations of the parsimony assumption using a data set of 21,488 noncoding single nucleotide polymorphisms (SNPs) discovered by the National Institute of Environmental Health Sciences (NIEHS) SNPs project via direct resequencing of n = 95 individuals. Applying standard nonparametric and parametric population genetic approaches, we replicate the signatures of a fixation bias in favor of G and C alleles when the ancestral base is assumed to be the base found in the chimpanzee outgroup. However, upon taking into account the probability of misidentifying the ancestral state of each SNP using a context-dependent mutation model, the corrected distribution of SNP frequencies for GC-content increasing SNPs are nearly indistinguishable from the patterns observed for other types of mutations, suggesting that the signature of fixation bias is a spurious artifact of the parsimony assumption.     

Ancestral alleles and population origins: inferences depend on mutation rate.

Previous studies have found that at most human loci, ancestral alleles are "African," in the sense that they reach their highest frequency there. Conventional wisdom holds that this reflects a recent African origin of modern humans. This paper challenges that view by showing that the empirical pattern (of elevated allele frequencies within Africa) is not as pervasive as has been thought. We confirm this African bias in a set of mainly protein-coding loci, but find a smaller bias in Alu insertion polymorphisms, and an even smaller bias in noncoding loci. Thus, the strong bias that was originally observed must reflect some factor that varies among data sets--something other than population history. This factor may be the per-locus mutation rate: the African bias is most pronounced in loci where this rate is high. The distribution of ancestral alleles among populations has been studied using 2 methods. One of these involves comparing the fractions of loci that reach maximal frequency in each population. The other compares the average frequencies of ancestral alleles. The first of these methods reflects history in a manner that depends on the mutation rate. When that rate is high, ancestral alleles at most loci reach their highest frequency in the ancestral population. When that rate is low, the reverse is true. The other method--comparing averages--is unresponsive. Average ancestral allele frequencies are affected neither by mutation rate nor by the history of population size and migration. In the absence of selection and ascertainment bias, they should be the same everywhere. This is true of one data set, but not of 2 others. This also suggests the action of some factor, such as selection or ascertainment bias, that varies among data sets.     

Type 2 diabetes risk alleles demonstrate extreme directional differentiation among human populations, compared to other diseases

Many disease-susceptible SNPs exhibit significant disparity in ancestral and derived allele frequencies across worldwide populations. While previous studies have examined population differentiation of alleles at specific SNPs, global ethnic patterns of ensembles of disease risk alleles across human diseases are unexamined. To examine these patterns, we manually curated ethnic disease association data from 5,065 papers on human genetic studies representing 1,495 diseases, recording the precise risk alleles and their measured population frequencies and estimated effect sizes. We systematically compared the population frequencies of cross-ethnic risk alleles for each disease across 1,397 individuals from 11 HapMap populations, 1,064 individuals from 53 HGDP populations, and 49 individuals with whole-genome sequences from 10 populations. Type 2 diabetes (T2D) demonstrated extreme directional differentiation of risk allele frequencies across human populations, compared with null distributions of European-frequency matched control genomic alleles and risk alleles for other diseases. Most T2D risk alleles share a consistent pattern of decreasing frequencies along human migration into East Asia. Furthermore, we show that these patterns contribute to disparities in predicted genetic risk across 1,397 HapMap individuals, T2D genetic risk being consistently higher for individuals in the African populations and lower in the Asian populations, irrespective of the ethnicity considered in the initial discovery of risk alleles. We observed a similar pattern in the distribution of T2D Genetic Risk Scores, which are associated with an increased risk of developing diabetes in the Diabetes Prevention Program cohort, for the same individuals. This disparity may be attributable to the promotion of energy storage and usage appropriate to environments and inconsistent energy intake. Our results indicate that the differential frequencies of T2D risk alleles may contribute to the observed disparity in T2D incidence rates across ethnic populations.     

Polymorphic NumtS trace human population relationships

The human genome is constantly subjected to evolutionary forces which shape its architecture. Insertions of mitochondrial DNA sequences into nuclear genome (NumtS) have been described in several eukaryotic species, including Homo sapiens and other primates. The ongoing process of the generation of NumtS has made them valuable markers in primate phylogenetic studies, as well as potentially informative loci for reconstructing the genetic history of modern humans. Here, we report the identification of 53 human-specific NumtS by inspection of the UCSC genome browser, showing that they may be direct insertions of mitochondrial DNA into the human nuclear DNA after the human-chimpanzee split. In silico analyses allowed us to identify 14 NumtS which are polymorphic in terms of their presence/absence within the human genome in individuals of different ancestry. The allele frequencies of these polymorphic NumtS were calculated for 1000 Genomes Project sequence data from 13 populations worldwide, and principal components analysis and hierarchical clustering methods allowed the detection of strong signals of geographical structure related to the genetic diversity of these loci. All identified polymorphic human-specific NumtS together with a tandemly duplicated NumtS have also been validated by PCR amplification on a panel of 60 samples belonging to five native populations worldwide, confirming the expected NumtS variability. On the basis of these findings, we have succeeded in depicting the landscape of variation of a series of NumtS in several ethnic groups, making an advance in their identification as useful markers in the study on human population genetics.     

Molecular Views of Human Origins

Over the last half century, comparative genomics has increasingly contributed to the definition, resolution and interpretation of human evolution. Early comparisons demonstrated that African apes and humans were more closely related and diverged later than commonly thought. However, it was difficult to determine the branching between humans, chimpanzees and gorillas. By the 1990s, sufficient biomolecular data had accumulated to demonstrate that chimpanzees and humans shared a common ancestor after the divergence of the gorilla. Current reconstructions place the divergence of humans and chimpanzees at 6–8 million years. Comparative genomics from complete genome sequencing to chromosome painting provide a scenario for the origin of the human genome. Starting form the ancestral mammalian karyotype, we can determine the major steps over the last 90 million years leading to the formation of each human chromosome. Despite considerable technical problems, studies of ancient DNA now provide a direct genetic witness of human evolution and add a temporal dimension to reconstructions of our evolutionary history and phylogeny. Ancient DNA has shown that Neanderthals probably did not interbreed with anatomically modern humans and did not make a significant contribution to the gene pool of our species. Ancient DNA has also contributed to the studies of the colonization of the Americas and the Pacific Island, and the domestication of plants and animals. Understanding the genetic basis of the physical and behavioral traits that distinguish humans from other primates presents one of the great future challenges of science.     

Human genome diversity

Human genome diversity studies analyse genetic variation among individuals and between populations in order to understand the origins and evolution of anatomically modern humans (Homo sapiens sapiens). The availability of thousands of DNA polymorphisms (genetic markers) brings analytic power to these studies. Human genome diversity studies have clearly shown that the large part of genetic variability is due to differences among individuals within populations rather than to differences between populations, effectively discrediting a genetic basis of the concept of ‘race’. Evidence from paleontology, archaeology and genetic diversity studies is quite consistent with an African origin of modern humans more than 100 000 years ago. The evidence favors migrations out of African as the source of the original peopling of Asia, Australia, Europe and Oceania. An international program for the scientific analysis of human genome diversity and of human evolution has been developed. The Human Genome Diversity Project (HGDP) aims to collect and preserve biologic samples from hundreds of populations throughout the world, make DNA from these samples available to scientists and distribute to the scientific community the results of DNA typing with hundreds of genetic markers.     

Low frequency of the Mx allele for viral resistance predates recent intensive selection in domestic chickens

Avian influenza is a serious threat to the poultry industry and, as the potential source of a human pandemic virus, to public health. Different Mx alleles have been reported to confer resistance or susceptibility to influenza virus replication, and so knowledge of their frequencies is important when considering the potential for improvement of modern commercial flocks. We analysed a range of chicken lines and ancestral breeds for the relevant Mx codon that confers resistance or susceptibility to influenza virus replication. We confirmed the high frequency of the susceptibility allele in contemporary meat-type (broiler) birds compared to egg-laying strains and found this difference is present already in ancestral breeds. We sequenced full-length complementary DNA (cDNA) and noted additional substitutions, which may be associated with the resistance haplotypes. High frequencies of the susceptibility allele could be readily reduced by modern breeding techniques.     

Spatial patterns of human gene frequencies in Europe

The aims of this study of spatial patterns of human gene frequencies in Europe are twofold. One is to present new methodology developed for the analysis of such data. The other is to report on the diversity of spatial patterns observed in Europe and their interpretation as evidence of population processes. Spatial variation in 59 allele and haplotype frequencies (26 genetic systems) for polymorphisms in blood antigens, enzymes, and proteins is analyzed for an aggregate of 3,384 localities, using homogeneity tests, one-dimensional and directional spatial correlograms, and SYMAP interpolated surfaces. The data matrices are reduced to reveal the principal patterns by clustering techniques. The findings of this study can be summarized as follows: 1) There is significant heterogeneity in allele frequencies among the localities for all but one genetic system. 2) There are significant spatial patterns for most allele frequencies. 3) There is a substantial minority of clinal patterns in these populations. Clinal trends are found more frequently in HLA alleles than for other variables. North-south and northwest-southwest gradients predominate. 4) There is a strong decline in overall genetic similarity with geographic distance for most variables. 5) There are few, if any, appreciable correlations in pairs of allele frequencies over the continent, and there is little interesting correlation structure in the resulting correlation matrix. 6) Few spatial correlograms are markedly similar to each other, yet they form well-defined clusters. Spatial variation patterns, therefore, differ among allele frequencies. Patterns of human gene frequencies in modern Europe are diverse and complex. No single model suffices for interpretation of the observed genetic structure. Some clinal patterns reported here support the Neolithic demic-expansion hypothesis, others suggest latitudinal selection. Most of the clinal patterns are in HLA alleles, but there is also evidence from ABO for east-west migration diffusion. The majority of patterns are patchy, consistent with hypotheses of isolation by distance or of settlement of genetically differing, subsequently expanding ethnic groups. While undoubtedly there has been an ongoing stochastic process of differentiation consistent with the isolation-by-distance model, this has not obscured the directional patterns caused by migration (demic diffusion), and has perhaps only reinforced the contribution from settlement of ethnic units to patterns of genetic variation. However, the impact of the latter is most difficult to discern and requires further methodological developments.     

Genetic aspects in hominid evolution

Genomic comparison between apes and humans have made important contributions to our understanding of human evolution. The modern period of karyological comparisons between humans and other primates began about forty years ago and has been marked by a series of technical revolutions. In the 1960s pioneering genetic and chromosomal comparisons of human and great apes suggested, as had Darwin a century before, that our closest relative were the African apes. Early immunological analyses placed human/apes divergence at about five million year ago. Acceptance of man’s late divergence from the African apes was delayed by the scarcity of paleontological evidence coupled with a fallacious Asiatic origin hypothesis of the hominoids. Chromosome banding techniques in the seventies and high resolution methods in the eighties allowed a detailed comparison of the chromosomes between closely related primates and reinforced the hypothesis of an African origin for humans. It was clearly shown that humans were more closely related to African apes than to the orang-utan. The last decade has seen a vigorous integration of molecular and cytogenetic. This powerful combination promises to be quite fruitful because chromosomes can be compared directly at the DNA level. Fluorescentin situ hybridisation (FISH), chromosome painting, is a colourful technique for establishing chromosomal homology between species. Results obtained by FISH over the last ten years have resolved the cytogenetic problem of the homology between humans, apes, hylobates and Old World monkeys and defined the chromosomal syntenies and major translocations involved in the genome evolution of higher primates.     

Evolution of Haplotypes at the DRD2 Locus

We present here the first evolutionary perspective on haplotypes at DRD2, the locus for the dopamine D₂ receptor. The dopamine D₂ receptor plays a critical role in the functioning of many neural circuits in the human brain. If functionally relevant variation at the DRD2 locus exists, understanding the evolution of haplotypes on the basis of polymorphic sites encompassing the gene should provide a powerful framework for identifying that variation. Three DRD2 polymorphisms (TaqI “A” and “B” RFLPs and the (CA)_n short tandem repeat polymorphism) encompassing the coding sequences have been studied in 15 populations; these markers are polymorphic in all the populations studied, and they display strong and significant linkage disequilibria with each other. The common haplotypes for the two TaqI RFLPs are separately derived from the ancestral haplotype but predate the spread of modern humans around the world. The knowledge of how the various haplotypes have evolved, the allele frequencies of the haplotypes in human populations, and the physical relationships of the polymorphisms to each other and to the functional parts of the gene should now allow proper design and interpretation of association studies.     

A Functional MiR-124 Binding-Site Polymorphism in IQGAP1 Affects Human Cognitive Performance

As a product of the unique evolution of the human brain, human cognitive performance is largely a collection of heritable traits. Rather surprisingly, to date there have been no reported cases to highlight genes that underwent adaptive evolution in humans and which carry polymorphisms that have a marked effect on cognitive performance. IQ motif containing GTPase activating protein 1 (IQGAP1), a scaffold protein, affects learning and memory in a dose-dependent manner. Its expression is regulated by miR-124 through the binding sites in the 3′UTR, where a SNP (rs1042538) exists in the core-binding motif. Here we showed that this SNP can influence the miR-target interaction both in vitro and in vivo. Individuals carrying the derived T alleles have higher IQGAP1 expression in the brain as compared to the ancestral A allele carriers. We observed a significant and male-specific association between rs1042538 and tactile performances in two independent cohorts. Males with the derived allele displayed higher tactual performances as compared to those with the ancestral allele. Furthermore, we found a highly diverged allele-frequency distribution of rs1042538 among world human populations, likely caused by natural selection and/or recent population expansion. These results suggest that current human populations still carry sequence variations that affect cognitive performances and that these genetic variants may likely have been subject to comparatively recent natural selection.     

Cytochrome P4501A1 polymorphisms in the Amerindian and Mestizo populations of Mexico

Several polymorphisms in the CYP1A1 locus have been identified and their genotypes appear to exhibit population frequencies that depend on ethnicity. We studied two CYP1A1 polymorphic sites (position 4889 and 6235) in a group of 212 unrelated healthy individuals belonging to three different Mexican populations (106 Mexican Mestizos, 52 Teenek and 54 Mayos). Comparison among Mexican populations showed increased frequency of the *Ile allele (A on position 4889) in Mexican Mestizos when compared to Amerindians (p < 0.05). The analysis of position 6235 showed increased frequencies of *m2 (C in this position) allele in Teenek when compared to Mestizos and Mayos (p < 0.05) and of *m2/*m2 genotype when compared to Mestizos (p < 0.05). Amerindian populations (from Mexico and South America) presented the lowest frequencies of *Ile (position 4889) and *m1 (position 6235) alleles, however these frequencies vary according to the ethnic group studied. Mexican Amerindian groups together with other South Amerindian populations showed the highest frequencies for *Val at position 4889 and the *m2 allele at position 6235. The present study corroborates the high frequencies of*Val and *m2 alleles in the Amerindian populations and detects some differences between Mexican populations that correlate with linguistic differences. Our data could be helpful in understanding the distribution of these polymorphisms and in clarifying their roles as genetic and evolution markers in Amerindian populations.     

Evidence of amino acid diversity-enhancing selection within humans and among primates at the candidate sperm-receptor gene PKDREJ

Sperm-egg interaction is a crucial step in fertilization, yet the identity of most interacting sperm-egg proteins that mediate this process remains elusive. Rapid evolution of some fertilization proteins has been observed in a number of species, including evidence of positive selection in the evolution of components of the mammalian egg coat. The rapid evolution of the egg-coat proteins could strongly select for changes on the sperm receptor, to maintain the interaction. Here, we present evidence that positive selection has driven the evolution of PKDREJ, a candidate sperm receptor of mammalian egg-coat proteins. We sequenced PKDREJ from a panel of 14 primates, including humans, and conducted a comparative maximum-likelihood analysis of nucleotide changes and found evidence of positive selection. An additional panel of 48 humans was surveyed for nucleotide polymorphisms at the PKDREJ locus. The regions predicted to have been subject to adaptive evolution among primates show several amino acid polymorphisms within humans. The distribution of polymorphisms suggests that balancing selection may maintain diverse PKDREJ alleles in some populations. It remains unknown whether there are functional differences associated with these diverse alleles, but their existence could have consequences for human fertility.