A high-resolution map of human evolutionary constraint using 29 mammals

The comparison of related genomes has emerged as a powerful lens for genome interpretation. Here we report the sequencing and comparative analysis of 29 eutherian genomes. We confirm that at least 5.5% of the human genome has undergone purifying selection, and locate constrained elements covering ～4.2% of the genome. We use evolutionary signatures and comparisons with experimental data sets to suggest candidate functions for ～60% of constrained bases. These elements reveal a small number of new coding exons, candidate stop codon readthrough events and over 10,000 regions of overlapping synonymous constraint within protein-coding exons. We find 220 candidate RNA structural families, and nearly a million elements overlapping potential promoter, enhancer and insulator regions. We report specific amino acid residues that have undergone positive selection, 280,000 non-coding elements exapted from mobile elements and more than 1,000 primate- and human-accelerated elements. Overlap with disease-associated variants indicates that our findings will be relevant for studies of human biology, health and disease.     

A haplotype map of the human genome

Inherited genetic variation has a critical but as yet largely uncharacterized role in human disease. Here we report a public database of common variation in the human genome: more than one million single nucleotide polymorphisms (SNPs) for which accurate and complete genotypes have been obtained in 269 DNA samples from four populations, including ten 500-kilobase regions in which essentially all information about common DNA variation has been extracted. These data document the generality of recombination hotspots, a block-like structure of linkage disequilibrium and low haplotype diversity, leading to substantial correlations of SNPs with many of their neighbours. We show how the HapMap resource can guide the design and analysis of genetic association studies, shed light on structural variation and recombination, and identify loci that may have been subject to natural selection during human evolution.     

Scanning for sisnatures of geosraphically restricted selection based on population genomics analysis

Natural selection, as the driving force of human evolution, has direct impact on population differentiation. However, it is still unclear to what extent the genetic differentiation has been caused by natural selection. To explore this question, we performed a genome-wide scan with single nucleotide polymorphism （SNP） data from the International HapMap Project. Single locus FsTanalysis was applied to assess the frequency difference among populations in autosomes. Based on the empirical distribution of FsT, we identified 12669 SNPs correlating to population differentiation and 1853 candidate genes subjected to geographic restricted natural selection. Further interpretation of gene ontogeny revealed 121 categories of biological process with the enrichments of candidate genes. Our results suggest that natural selection may play an important role in human population differentiation. In addition, our analysis provides new clues as well as research methods for our understanding of population differentiation and natural selection.     

Mapping determinants of human gene expression by regional and genome-wide association

To study the genetic basis of natural variation in gene expression, we previously carried out genome-wide linkage analysis and mapped the determinants of approximately 1,000 expression phenotypes. In the present study, we carried out association analysis with dense sets of single-nucleotide polymorphism (SNP) markers from the International HapMap Project. For 374 phenotypes, the association study was performed with markers only from regions with strong linkage evidence; these regions all mapped close to the expressed gene. For a subset of 27 phenotypes, analysis of genome-wide association was performed with >770,000 markers. The association analysis with markers under the linkage peaks confirmed the linkage results and narrowed the candidate regulatory regions for many phenotypes with strong linkage evidence. The genome-wide association analysis yielded highly significant results that point to the same locations as the genome scans for about 50% of the phenotypes. For one candidate determinant, we carried out functional analyses and confirmed the variation in cis-acting regulatory activity. Our findings suggest that association studies with dense SNP maps will identify susceptibility loci or other determinants for some complex traits or diseases.     

Detecting recent positive selection in the human genome from haplotype structure

The ability to detect recent natural selection in the human population would have profound implications for the study of human history and for medicine. Here, we introduce a framework for detecting the genetic imprint of recent positive selection by analysing long-range haplotypes in human populations. We first identify haplotypes at a locus of interest (core haplotypes). We then assess the age of each core haplotype by the decay of its association to alleles at various distances from the locus, as measured by extended haplotype homozygosity (EHH). Core haplotypes that have unusually high EHH and a high population frequency indicate the presence of a mutation that rose to prominence in the human gene pool faster than expected under neutral evolution. We applied this approach to investigate selection at two genes carrying common variants implicated in resistance to malaria: G6PD and CD40 ligand. At both loci, the core haplotypes carrying the proposed protective mutation stand out and show significant evidence of selection. More generally, the method could be used to scan the entire genome for evidence of recent positive selection.     

Adaptive evolution drives divergence of a hybrid inviability gene between two species of Drosophila

Speciation--the splitting of one species into two--occurs by the evolution of any of several forms of reproductive isolation between taxa, including the intrinsic sterility and inviability of hybrids. Abundant evidence shows that these hybrid fitness problems are caused by incompatible interactions between loci: new alleles that become established in one species are sometimes functionally incompatible with alleles at interacting loci from another species. However, almost nothing is known about the genes involved in such hybrid incompatibilities or the evolutionary forces that drive their divergence. Here we identify a gene that causes epistatic inviability in hybrids between two fruitfly species, Drosophila melanogaster and D. simulans. Our population genetic analysis reveals that this gene--which encodes a nuclear pore protein--evolved by positive natural selection in both species' lineages. These results show that a lethal hybrid incompatibility has evolved as a by-product of adaptive protein evolution.     

Fine-scale recombination rate differences between sexes, populations and individuals

Meiotic recombinations contribute to genetic diversity by yielding new combinations of alleles. Recently, high-resolution recombination maps were inferred from high-density single-nucleotide polymorphism (SNP) data using linkage disequilibrium (LD) patterns that capture historical recombination events. The use of these maps has been demonstrated by the identification of recombination hotspots and associated motifs, and the discovery that the PRDM9 gene affects the proportion of recombinations occurring at hotspots. However, these maps provide no information about individual or sex differences. Moreover, locus-specific demographic factors like natural selection can bias LD-based estimates of recombination rate. Existing genetic maps based on family data avoid these shortcomings, but their resolution is limited by relatively few meioses and a low density of markers. Here we used genome-wide SNP data from 15,257 parent-offspring pairs to construct the first recombination maps based on directly observed recombinations with a resolution that is effective down to 10 kilobases (kb). Comparing male and female maps reveals that about 15% of hotspots in one sex are specific to that sex. Although male recombinations result in more shuffling of exons within genes, female recombinations generate more new combinations of nearby genes. We discover novel associations between recombination characteristics of individuals and variants in the PRDM9 gene and we identify new recombination hotspots. Comparisons of our maps with two LD-based maps inferred from data of HapMap populations of Utah residents with ancestry from northern and western Europe (CEU) and Yoruba in Ibadan, Nigeria (YRI) reveal population differences previously masked by noise and map differences at regions previously described as targets of natural selection.     

Maintaining a behaviour polymorphism by frequency-dependent selection on a single gene

Accounting for the abundance of genetic variation in the face of natural selection remains a central problem of evolutionary biology. Genetic polymorphisms are constantly arising through mutation, and although most are promptly eliminated, polymorphisms in functionally important traits are common. One mechanism that can maintain polymorphisms is negative frequency-dependent selection on alternative alleles, whereby the fitness of each decreases as its frequency increases. Examples of frequency-dependent selection are rare, especially when attempting to describe the genetic basis of the phenotype under selection. Here we show frequency-dependent selection in a well-known natural genetic polymorphism affecting fruitfly foraging behaviour. When raised in low nutrient conditions, both of the naturally occurring alleles of the foraging gene (for(s) and for(R)) have their highest fitness when rare-the hallmark of negative frequency-dependent selection. This effect disappears at higher resources levels, demonstrating the role of larval competition. We are able to confirm the involvement of the foraging gene by showing that a sitter-like mutant allele on a rover background has similar frequency-dependent fitness as the natural sitter allele. Our study represents a clear demonstration of frequency-dependent selection, and we are able to attribute this effect to a single, naturally polymorphic gene known to affect behaviour.     

应用DHPLC检测多巴胺D2受体基因在中国人群中的多态性

通过变性高效液相色谱 (DHPLC)和DNA测序在 4 1个中国汉族人DNA样本中检测DRD2基因编码区和拼接区的单核苷酸多态性 (SNP) ,结果发现 3个SNP :Intron5的 2 77G/A、Exon7的 4 2C/T和Exon7的 1 2 9T/C .Intron5 2 77G/A是在中国汉族人群中发现的新SNP ,Exon7 4 2C/T和 1 2 9T/C在NCBIdbSNP中已有相应记录 (分别为rs4 986 92 1和rs6 2 75) ,它们均导致DRD2基因的同义突变 ,其中Exon7 1 2 9C等位基因频率在研究样本中高达 4 3.9% .这些结果为在中国汉族人群中开展DRD2基因相关的群体遗传学研究提供了遗传标记 .另外 ,还探讨了DHPLC检测突变的干扰因素及控制措施 ,为国内同行开展类似工作提供参考     

Initial sequence of the chimpanzee genome and comparison with the human genome

Here we present a draft genome sequence of the common chimpanzee (Pan troglodytes). Through comparison with the human genome, we have generated a largely complete catalogue of the genetic differences that have accumulated since the human and chimpanzee species diverged from our common ancestor, constituting approximately thirty-five million single-nucleotide changes, five million insertion/deletion events, and various chromosomal rearrangements. We use this catalogue to explore the magnitude and regional variation of mutational forces shaping these two genomes, and the strength of positive and negative selection acting on their genes. In particular, we find that the patterns of evolution in human and chimpanzee protein-coding genes are highly correlated and dominated by the fixation of neutral and slightly deleterious alleles. We also use the chimpanzee genome as an outgroup to investigate human population genetics and identify signatures of selective sweeps in recent human evolution.     

A second generation human haplotype map of over 3.1 million SNPs

We describe the Phase II HapMap, which characterizes over 3.1 million human single nucleotide polymorphisms (SNPs) genotyped in 270 individuals from four geographically diverse populations and includes 25-35% of common SNP variation in the populations surveyed. The map is estimated to capture untyped common variation with an average maximum r2 of between 0.9 and 0.96 depending on population. We demonstrate that the current generation of commercial genome-wide genotyping products captures common Phase II SNPs with an average maximum r2 of up to 0.8 in African and up to 0.95 in non-African populations, and that potential gains in power in association studies can be obtained through imputation. These data also reveal novel aspects of the structure of linkage disequilibrium. We show that 10-30% of pairs of individuals within a population share at least one region of extended genetic identity arising from recent ancestry and that up to 1% of all common variants are untaggable, primarily because they lie within recombination hotspots. We show that recombination rates vary systematically around genes and between genes of different function. Finally, we demonstrate increased differentiation at non-synonymous, compared to synonymous, SNPs, resulting from systematic differences in the strength or efficacy of natural selection between populations.     

Tracking the evolution of insecticide resistance in the mosquito Culex pipiens.

The evolution of pesticide resistance provides some of the most striking examples of darwinian evolution occurring over a human life span. Identification of resistance alleles opens an outstanding framework in which to study the evolution of adaptive mutations from the beginning of pesticide application, the evolution of interactions between alleles (dominance) or between loci (epistasis). Here we show that resistance alleles can also be used as markers to dissect population processes at a microevolutionary scale. We have focused on the antagonistic roles of selection and migration involved in the dynamics of local adaptation with reference to allelic frequencies at two resistance loci in the mosquito Culex pipiens. We find that their frequencies follow an annual cycle of large amplitude (25%), and we precisely unravel the seasonal variation of migration and selection underlying this cycle. Our results provide a firm basis on which to devise an insecticide treatment strategy that will better control the evolution of resistance genes and the growth of mosquito populations.     

Contribution of Distal-less to quantitative variation in butterfly eyespots.

The colour patterns decorating butterfly wings provide ideal material to study the reciprocal interactions between evolution and development. They are visually compelling products of selection, often with a clear adaptive value, and are amenable to a detailed developmental characterization. Research on wing-pattern evolution and development has focused on the eyespots of the tropical butterfly Bicyclus anynana. There is quantitative variation for several features of eyespot morphology but the actual genes contributing to such variation are unknown. On the other hand, studies of gene expression patterns in wing primordia have implicated different developmental pathways in eyespot formation. To link these two sets of information we need to identify which genes within the implicated pathways contribute to the quantitative variation accessible to natural selection. Here we begin to bridge this gap by demonstrating linkage between DNA polymorphisms in the candidate gene Distal-less (Dll) and eyespot size in B. anynana.     

Insulin-IGF2 region on chromosome 11p encodes a gene implicated in HLA-DR4-dependent diabetes susceptibility

A class of alleles at the VNTR (variable number of tandem repeat) locus in the 5' region of the insulin gene (INS) on chromosome 11p is associated with increased risk of insulin-dependent diabetes mellitus (IDDM), but family studies have failed to demonstrate linkage. INS is thought to contribute to IDDM susceptibility but this view has been difficult to reconcile with the lack of linkage evidence. We thus investigated polymorphisms of INS and neighbouring loci in random diabetics, IDDM multiplex families and controls. HLA-DR4-positive diabetics showed an increased risk associated with common variants at polymorphic sites in a 19-kilobase segment spanned by the 5' INS VNTR and the third intron of the gene for insulin-like growth factor II (IGF2). As INS is the major candidate gene from this region, diabetic and control sequence were compared to identify all INS polymorphisms that could contribute to disease susceptibility. In multiplex families the IDDM-associated alleles were transmitted preferentially to HLA-DR4-positive diabetic offspring from heterozygous parents. The effect was strongest in paternal meioses, suggesting a possible role for maternal imprinting. Our results strongly support the existence of a gene or genes affecting HLA-DR4 IDDM susceptibility which is located in a 19-kilobase region of INS-IGF2. Our results also suggest new ways to map susceptibility loci in other common diseases.     

Linkage disequilibrium in the human genome

With the availability of a dense genome-wide map of single nucleotide polymorphisms (SNPs), a central issue in human genetics is whether it is now possible to use linkage disequilibrium (LD) to map genes that cause disease. LD refers to correlations among neighbouring alleles, reflecting 'haplotypes' descended from single, ancestral chromosomes. The size of LD blocks has been the subject of considerable debate. Computer simulations and empirical data have suggested that LD extends only a few kilobases (kb) around common SNPs, whereas other data have suggested that it can extend much further, in some cases greater than 100 kb. It has been difficult to obtain a systematic picture of LD because past studies have been based on only a few (1-3) loci and different populations. Here, we report a large-scale experiment using a uniform protocol to examine 19 randomly selected genomic regions. LD in a United States population of north-European descent typically extends 60 kb from common alleles, implying that LD mapping is likely to be practical in this population. By contrast, LD in a Nigerian population extends markedly less far. The results illuminate human history, suggesting that LD in northern Europeans is shaped by a marked demographic event about 27,000-53,000 years ago.     

Paternal inheritance of a female moth's mating preference

Females of the arctiid moth Utetheisa ornatrix mate preferentially with larger males, receiving both direct phenotypic and indirect genetic benefits. Here we demonstrate that the female's mating preference is inherited through the father rather than the mother, indicating that the preference gene or genes lie mostly or exclusively on the Z sex chromosome, which is strictly paternally inherited by daughters. Furthermore, we show that the preferred male trait and the female preference for that trait are correlated, as females with larger fathers have a stronger preference for larger males. These findings are predicted by the protected invasion theory, which asserts that male homogametic sex chromosome systems (ZZ/ZW) found in lepidopterans and birds promote the evolution of exaggerated male traits through sexual selection. Specifically, the theory predicts that, because female preference alleles arising on the Z chromosome are transmitted to all sons that have the father's attractive trait rather than to only a fraction of the sons, such alleles will experience stronger positive selection and be less vulnerable to chance loss than would autosomal alleles.     

Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection

The major histocompatibility complex (MHC) loci are known to be highly polymorphic in humans, mice and certain other mammals, with heterozygosity as high as 80-90% (ref. 1). Four different hypotheses have been proposed to explain this high degree of polymorphism: (1) a high mutation rate, (2) gene conversion or interlocus genetic exchange, (3) over dominant selection and (4) frequency-dependent selection. In an attempt to establish which of these hypotheses is correct, we examined the pattern of nucleotide substitution between polymorphic alleles in the region of the antigen recognition site (ARS) and other regions of human and mouse class I MHC genes. The results indicate that in ARS the rate of nonsynonymous (amino acid altering) substitution is significantly higher than that of synonymous substitution in both humans and mice, whereas in other regions the reverse is true. This observation, together with a theoretical study and other considerations, supports the hypothesis of overdominant selection (heterozygote advantage).     

Positive selection at sites of multiple amino acid replacements since rat-mouse divergence

New alleles become fixed owing to random drift of nearly neutral mutations or to positive selection of substantially advantageous mutations. After decades of debate, the fraction of fixations driven by selection remains uncertain. Within 9,390 genes, we analysed 28,196 codons at which rat and mouse differ from each other at two nucleotide sites and 1,982 codons with three differences. At codons where rat-mouse divergence involved two non-synonymous substitutions, both of them occurred in the same lineage, either rat or mouse, in 64% of cases; however, independent substitutions would occur in the same lineage with a probability of only 50%. All three non-synonymous substitutions occurred in the same lineage for 46% of codons, instead of the 25% expected. Furthermore, comparison of 12 pairs of prokaryotic genomes also shows clumping of multiple non-synonymous substitutions in the same lineage. This pattern cannot be explained by correlated mutation or episodes of relaxed negative selection, but instead indicates that positive selection acts at many sites of rapid, successive amino acid replacement.     

美洲黑杨种质材料倍性鉴定

【目的】对美洲黑杨种质资源库保存的种质材料进行倍性评估,并对候选多倍体进行倍性鉴定。【方法】利用9对用于多倍体鉴定的SSR引物检测448份美洲黑杨种质材料,根据每个引物扩增位点的等位基因数目对种质材料的倍性进行评估,基于初步鉴定结果筛选出候选多倍体材料。在此基础上,使用流式细胞仪(FCM)对候选多倍体种质材料的倍性进行鉴定。【结果】分子标记基因型分型统计结果显示,9个SSR位点共扩增出129条多态性条带,每个位点的等位基因数为5(Ploidp-10)~25(Ploidp-07)个,平均每个位点的等位基因数为14.33个。9个SSR位点的多态性信息含量(PIC)变动范围为0.46~0.93,平均PIC值为0.76。在检测的448份美洲黑杨种质材料中,440份美洲黑杨种质材料在9个引物位点上扩增出的等位基因数不多于2个;有8份种质材料在2个引物位点上扩增出了3个不同的等位基因,其中编号为1346、16-42、17-49、18-25的4份种质材料在Ploidp-02引物位点扩增出了3个等位基因,而编号为5-43、5-45、15-14、19-37的另外4份种质材料在Ploidp-04引物位点上也扩增出了3个等位基因。分子检测结果表明:该8份种质材料在对应的染色体区域发生了遗传物质增加现象,为候选三倍体种质材料。用FCM对8份三倍体候选种质材料进行倍性鉴定发现,编号为1346的种质材料为三倍体,其余7份种质材料疑似为发生局部染色片段增加的候选非整倍体。【结论】9对SSR分子标记均具有较高的多态性信息含量,可用于美洲黑杨种质资源多倍体初步筛选,结合FCM检测,可在大量样品中快速筛选出多倍体种质材料。     

Calcineurin sets the bandwidth for discrimination of signals during thymocyte development

At critical times in development, cells are able to convert graded signals into discrete developmental outcomes; however, the mechanisms involved are poorly understood. During thymocyte development, cell fate is determined by signals originating from the alphabeta T-cell receptor. Low-affinity/avidity interactions between the T-cell receptor and peptide-MHC complexes direct differentiation to the single-positive stage (positive selection), whereas high-affinity/avidity interactions induce death by apoptosis (negative selection). Here we show that mice deficient in both calcineurin and nuclear factor of activated T cells (NFAT)c2/c3 lack a population of preselection thymocytes with enhanced ability to activate the mitogen-activated protein kinase (Raf-MEK-ERK) pathway, and fail to undergo positive selection. This defect can be partially rescued with constitutively active Raf, indicating that calcineurin controls MAPK signalling. Analysis of mice deficient in both Bim (which is required for negative selection) and calcineurin revealed that calcineurin-induced ERK (extracellular signal-regulated kinase) sensitization is required for differentiation in response to 'weak' positive selecting signals but not in response to 'strong' negative selecting signals (which normally induce apoptosis). These results indicate that early calcineurin/NFAT signalling produces a developmental period of ERK hypersensitivity, allowing very weak signals to induce positive selection. This mechanism might be generally useful in the discrimination of graded signals that induce different cell fates.