Derived Allele Frequency Research Articles

Strong Heterogeneity in Mutation Rate Causes Misleading Hallmarks of Natural Selection on Indel Mutations in the Human Genome

Elucidating the mechanisms of mutation accumulation and fixation is critical to understand the nature of genetic variation and its contribution to genome evolution. Of particular interest is the effect of insertions and deletions (indels) on the evolution of genome landscapes. Recent population-scaled sequencing efforts provide unprecedented data for analyzing the relative impact of selection versus nonadaptive forces operating on indels. Here, we combined McDonald–Kreitman tests with the analysis of derived allele frequency spectra to investigate the dynamics of allele fixation of short (1–50 bp) indels in the human genome. Our analyses revealed apparently higher fixation probabilities for insertions than deletions. However, this fixation bias is not consistent with either selection or biased gene conversion and varies with local mutation rate, being particularly pronounced at indel hotspots. Furthermore, we identified an unprecedented number of loci with evidence for multiple indel events in the primate phylogeny. Even in nonrepetitive sequence contexts (a priori not prone to indel mutations), such loci are 60-fold more frequent than expected according to a model of uniform indel mutation rate. This provides evidence of as yet unidentified cryptic indel hotspots. We propose that indel homoplasy, at known and cryptic hotspots, produces systematic errors in determination of ancestral alleles via parsimony and advise caution interpreting classic selection tests given the strong heterogeneity in indel rates across the genome. These results will have great impact on studies seeking to infer evolutionary forces operating on indels observed in closely related species, because such mutations are traditionally presumed homoplasy-free.

The influence of genomic context on mutation patterns in the human genome inferred from rare variants.

Understanding patterns of spontaneous mutations is of fundamental interest in studies of human genome evolution and genetic disease. Here, we used extremely rare variants in humans to model the molecular spectrum of single-nucleotide mutations. Compared to common variants in humans and human–chimpanzee fixed differences (substitutions), rare variants, on average, arose more recently in the human lineage and are less affected by the potentially confounding effects of natural selection, population demographic history, and biased gene conversion. We analyzed variants obtained from a population-based sequencing study of 202 genes in >14,000 individuals. We observed considerable variability in the per-gene mutation rate, which was correlated with local GC content, but not recombination rate. Using >20,000 variants with a derived allele frequency ≤10−4, we examined the effect of local GC content and recombination rate on individual variant subtypes and performed comparisons with common variants and substitutions. The influence of local GC content on rare variants differed from that on common variants or substitutions, and the differences varied by variant subtype. Furthermore, recombination rate and recombination hotspots have little effect on rare variants of any subtype, yet both have a relatively strong impact on multiple variant subtypes in common variants and substitutions. This observation is consistent with the effect of biased gene conversion or selection-dependent processes. Our results highlight the distinct biases inherent in the initial mutation patterns and subsequent evolutionary processes that affect segregating variants.

Modulation of PF10_0355 (MSPDBL2) Alters Plasmodium falciparum Response to Antimalarial Drugs

Malaria's ability to rapidly adapt to new drugs has allowed it to remain one of the most devastating infectious diseases of humans. Understanding and tracking the genetic basis of these adaptations are critical to the success of treatment and intervention strategies. The novel antimalarial resistance locus PF10_0355 (Pfmspdbl2) was previously associated with the parasite response to halofantrine, and functional validation confirmed that overexpression of this gene lowered parasite sensitivity to both halofantrine and the structurally related antimalarials mefloquine and lumefantrine, predominantly through copy number variation. Here we further characterize the role of Pfmspdbl2 in mediating the antimalarial drug response of Plasmodium falciparum. Knockout of Pfmspdbl2 increased parasite sensitivity to halofantrine, mefloquine, and lumefantrine but not to unrelated antimalarials, further suggesting that this gene mediates the parasite response to a specific class of antimalarial drugs. A single nucleotide polymorphism encoding a C591S mutation within Pfmspdbl2 had the strongest association with halofantrine sensitivity and showed a high derived allele frequency among Senegalese parasites. Transgenic parasites expressing the ancestral Pfmspdbl2 allele were more sensitive to halofantrine and structurally related antimalarials than were parasites expressing the derived allele, revealing an allele-specific effect on drug sensitivity in the absence of copy number effects. Finally, growth competition experiments showed that under drug pressure, parasites expressing the derived allele of Pfmspdbl2 outcompeted parasites expressing the ancestral allele within a few generations. Together, these experiments demonstrate that modulation of Pfmspdbl2 affects malaria parasite responses to antimalarial drugs.

Similarity in recombination rate and linkage disequilibrium at CYP2C and CYP2D cytochrome P450 gene regions among Europeans indicates signs of selection and no advantage of using tagSNPs in population isolates

Linkage disequilibrium (LD) and recombination rate variations are known to vary considerably between human genome regions and populations mostly because of the combined effects of mutation, recombination, and demographic history. Thus, the pattern of LD is a key issue to disentangle variants associated with complex traits. Here, we aim to describe the haplotype structure and LD variation at the pharmacogenetically relevant cytochrome P450 CYP2C and CYP2D gene regions among European populations. To assess the haplotype structure, LD pattern, and recombination rate variations in the clinically significant CYP2C and CYP2D regions, we genotyped 143 single-nucleotide polymorphisms (SNPs) across these two genome regions in a diverse set of 11 European population samples and one sub-Saharan African sample. Our results showed extended patterns of LD and in general a low rate of recombination at these loci, and a low degree of allele differentiation for the two cytochrome P450 regions among Europeans, with the exception of the Sami and the Finns as European outliers. The Sami sample showed reduced haplotype diversity and higher LD for the two cytochrome P450 regions than the other Europeans, a feature that is proposed to enhance the LD mapping of underlying common complex traits. However, recombination hotspots and LD blocks at these two regions showed highly consistent structures across Europeans including Finns and Sami. Moreover, we showed that the CEPH sample has significantly higher tag transferability among Europeans and a more efficient tagging of both the rare CYP2C9 and the common CYP2C19 functional variants than the Sami. Our data set included CYP2C9*3 (rs1057910) and CYP2C19*2 (rs4244285) enzyme activity-altering variants associated in a recent genome-wide study with acenocoumarol-induced and warfarin-induced anticoagulation or to the antiplatelet effect of clopidogrel, respectively. Including these known activity-altering variants, we showed the haplotype variation and high derived allele frequencies of novel recently identified acenocoumarol genome-wide associated SNPs at CYP2C9 (rs4086116) and CYP2C18 (rs12772169, rs1998591, rs2104543, rs1042194) loci in a comprehensive set of 11 European populations. Furthermore, a significant frequency difference of a CYP2C19*2 gene mutation causing variable drug reactions was observed among Europeans. The CEPH sample representing the general European population as such in the HapMap project seems to be the optimal population sample for the LD mapping of common complex traits among Europeans. Nevertheless, it is still argued that the unique pattern of LD in the Sami may offer an advantage for further association mapping, especially if multiple rare variants play a role in disease etiology. However, besides the activity-altering CYP2C9*3 (rs1057910) and CYP2C19*2 (rs4244285) variants, the high derived allele frequencies of novel recently identified acenocoumarol genome-wide associated SNPs at CYP2C9 (rs4086116) and CYP2C18 (rs12772169, rs1998591, rs2104543, rs1042194) loci variants indicated that the CYP2C region may have been influenced by selection. Thus, this fine-scale haplotype map of the CYP2C and CYP2D regions may help to choose markers for further association mapping of complex pharmacogenetic traits at these loci.

Mutational Bias and Gene Conversion Affect the Intraspecific Nitrogen Stoichiometry of the Arabidopsis thaliana Transcriptome

The transcriptome and proteome of Arabidopsis thaliana are reduced in nitrogen content when compared with other taxa, which may result from ecological nitrogen limitation. We hypothesized that if the A. thaliana transcriptome is selected for a low nitrogen content, nitrogen-reducing derived alleles of single nucleotide polymorphisms (SNPs) should segregate at higher frequencies than nitrogen-increasing alleles. This pattern should be stronger in populations with a larger effective population size (N(e)) if natural selection is more efficient in large than in small populations. We analyzed variation in the nitrogen content in the transcriptome of 80 natural accessions of A. thaliana. In contrast to our expectations, derived alleles increase the nitrogen content in all accessions, and there is a positive correlation between nitrogen difference and derived allele frequency, which is strongest with nonsynonymous SNPs (nsSNPs). Also, there is a positive correlation between nitrogen difference and N(e) that was mainly caused by nsSNPs. These observations led us to reject the hypothesis that the transcriptome of A. thaliana is currently under selection to reduce nitrogen content. Instead, we show that a change in nitrogen content is a side effect of interacting evolutionary factors that influence base composition and include mutational bias, purifying selection of functionally deleterious alleles, and GC-biased gene conversion. We provide strong evidence that GC-biased gene conversion may play an important role for base composition in the highly selfing plant A. thaliana.

Skin Color Variation in Orang Asli Tribes of Peninsular Malaysia

Pigmentation is a readily scorable and quantitative human phenotype, making it an excellent model for studying multifactorial traits and diseases. Convergent human evolution from the ancestral state, darker skin, towards lighter skin colors involved divergent genetic mechanisms in people of European vs. East Asian ancestry. It is striking that the European mechanisms result in a 10–20-fold increase in skin cancer susceptibility while the East Asian mechanisms do not. Towards the mapping of genes that contribute to East Asian pigmentation there is need for one or more populations that are admixed for ancestral and East Asian ancestry, but with minimal European contribution. This requirement is fulfilled by the Senoi, one of three indigenous tribes of Peninsular Malaysia collectively known as the Orang Asli. The Senoi are thought to be an admixture of the Negrito, an ancestral dark-skinned population representing the second of three Orang Asli tribes, and regional Mongoloid populations of Indo-China such as the Proto-Malay, the third Orang Asli tribe. We have calculated skin reflectance-based melanin indices in 492 Orang Asli, which ranged from 28 (lightest) to 75 (darkest); both extremes were represented in the Senoi. Population averages were 56 for Negrito, 42 for Proto-Malay, and 46 for Senoi. The derived allele frequencies for SLC24A5 and SLC45A2 in the Senoi were 0.04 and 0.02, respectively, consistent with greater South Asian than European admixture. Females and individuals with the A111T mutation had significantly lighter skin (p = 0.001 and 0.0039, respectively). Individuals with these derived alleles were found across the spectrum of skin color, indicating an overriding effect of strong skin lightening alleles of East Asian origin. These results suggest that the Senoi are suitable for mapping East Asian skin color genes.

Inferring the History of Population Size Change from Genome-Wide SNP Data

Dense, genome-wide single-nucleotide polymorphism (SNP) data can be used to reconstruct the demographic history of human populations. However, demographic inferences from such data are complicated by recombination and ascertainment bias. We introduce two new statistics, allele frequency-identity by descent (AF-IBD) and allele frequency-identity by state (AF-IBS), that make use of linkage disequilibrium information and show defined relationships to the time of coalescence. These statistics, when conditioned on the derived allele frequency, are able to infer complex population size changes. Moreover, the AF-IBS statistic, which is based on genome-wide SNP data, is robust to varying ascertainment conditions. We constructed an efficient approximate Bayesian computation (ABC) pipeline based on AF-IBD and AF-IBS that can accurately estimate demographic parameters, even for fairly complex models. Finally, we applied this ABC approach to genome-wide SNP data and inferred the demographic histories of two human populations, Yoruba and French. Our results suggest a rather stable ancestral population size with a mild recent expansion for Yoruba, whereas the French seemingly experienced a long-lasting severe bottleneck followed by a drastic population growth. This approach should prove useful for new insights into populations, especially those with complex demographic histories.

Association of melanocortin 1 receptor gene (MC1R) polymorphisms with skin reflectance and freckles in Japanese

Most studies on the genetic basis of human skin pigmentation have focused on people of European ancestry and only a few studies have focused on Asian populations. We investigated the association of skin reflectance and freckling with genetic variants of melanocortin 1 receptor (MC1R) gene in Japanese. DNA samples were obtained from a total of 653 Japanese individuals (ages 19-40 years) residing in Okinawa; skin reflectance was measured using a spectrophotometer and freckling status was determined for each individual. Lightness index (L*) and freckling status were not correlated with age, body mass index or ancestry (Ryukyuan or Main Islanders of Japan). Among the 10 nonsynonymous variants that were identified by direct sequencing of the coding region of MC1R, two variants--R163Q and V92M--with the derived allele frequencies of 78.6 and 5.5%, respectively, were most common. Multiple regression analysis showed that the 163Q allele and the presence of nonsynonymous rare variants (allele frequencies <5%) were significantly associated with an increase in sex-standardized skin lightness (L* of CIELAB (CIE 1976 (L*a*b*) color space)) of the inner upper arm. Relative to the 92V allele, the 92M allele was significantly associated with increased odds of freckling. This is the first study to show an association between the 163Q allele and skin reflectance values; this association indicated that light-toned skin may have been subjected to positive selection in East Asian people.

Identification and characterization of lineage-specific highly conserved noncoding sequences in Mammalian genomes.

Vertebrate genome comparisons revealed that there are highly conserved noncoding sequences (HCNSs) among a wide range of species and many of which contain regulatory elements. However, recently emerged sequences conserved in specific lineages have not been well studied. Toward this end, we identified 8,198 primate and 21,128 specific HCNSs as representative ones among mammals from human–marmoset and mouse–rat comparisons, respectively. Derived allele frequency analysis of primate-specific HCNSs showed that these HCNSs were under purifying selection, indicating that they may harbor important functions. We selected the top 1,000 largest HCNSs and compared the lineage-specific HCNS-flanking genes (LHF genes) with ultraconserved element (UCE)-flanking genes. Interestingly, the majority of LHF genes were different from UCE-flanking genes. This lineage-specific set of LHF genes was more enriched in protein-binding function. Conversely, the number of LHF genes that were also shared by UCEs was small but significantly larger than random expectation, and many of these genes were involved in anatomical development as transcriptional regulators, suggesting that certain groups of genes preferentially recruit new HCNSs in addition to old HCNSs that are conserved among vertebrates. This group of LHF genes might be involved in the various levels of lineage-specific evolution among vertebrates, mammals, primates, and rodents. If so, the emergence of HCNSs in and around these two groups of LHF genes developed lineage-specific characteristics. Our results provide new insight into lineage-specific evolution through interactions between HCNSs and their LHF genes.

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.

Genome-Wide Survey of Natural Selection on Functional, Structural, and Network Properties of Polymorphic Sites in Saccharomyces paradoxus

To characterize the genetic basis of phenotypic evolution, numerous studies have identified individual genes that have likely evolved under natural selection. However, phenotypic changes may represent the cumulative effect of similar evolutionary forces acting on functionally related groups of genes. Phylogenetic analyses of divergent yeast species have identified functional groups of genes that have evolved at significantly different rates, suggestive of differential selection on the functional properties. However, due to environmental heterogeneity over long evolutionary timescales, selection operating within a single lineage may be dramatically different, and it is not detectable via interspecific comparisons alone. Moreover, interspecific studies typically quantify selection on protein-coding regions using the D(n)/D(s) ratio, which cannot be extended easily to study selection on noncoding regions or synonymous sites. The population genetic-based analysis of selection operating within a single lineage ameliorates these limitations. We investigated selection on several properties associated with genes, promoters, or polymorphic sites, by analyzing the derived allele frequency spectrum of single nucleotide polymorphisms (SNPs) in 28 strains of Saccharomyces paradoxus. We found evidence for significant differential selection between many functionally relevant categories of SNPs, underscoring the utility of function-centric approaches for discovering signatures of natural selection. When comparable, our findings are largely consistent with previous studies based on interspecific comparisons, with one notable exception: our study finds that mutations from an ancient amino acid to a relatively new amino acid are selectively disfavored, whereas interspecific comparisons have found selection against ancient amino acids. Several of our findings have not been addressed through prior interspecific studies: we find that synonymous mutations from preferred to unpreferred codons are selected against and that synonymous SNPs in the linker regions of proteins are relatively less constrained than those within protein domains. We present the first global survey of selection acting on various functional properties in S. paradoxus. We found that selection pressures previously detected over long evolutionary timescales have also shaped the evolution of S. paradoxus. Importantly, we also make novel discoveries untenable via conventional interspecific analyses.

Interleukin 6 Promoter Polymorphisms in Normal Han Chinese Population

BackgroundInterleukin 6 (IL-6) is a cytokine involved in different physiologic and pathophysiologic processes, and circulating levels of IL-6 differ greatly between individuals, but the results have not always been concordant...

Predicting functionally important SNP classes based on negative selection.

BackgroundWith the advent of cost-effective genotyping technologies, genome-wide association studies allow researchers to examine hundreds of thousands of single nucleotide polymorphisms (SNPs) for association with human disease. Recently, many researchers applying this strategy have detected strong associations to disease with SNP markers that are either not in linkage disequilibrium with any nonsynonymous SNP or large distances from any annotated gene. In such cases, no well-established standard practice for effective SNP selection for follow-up studies exists. We aim to identify and prioritize groups of SNPs that are more likely to affect phenotypes in order to facilitate efficient SNP selection for follow-up studies.ResultsBased on the annotations available in the Ensembl database, we categorized SNPs in the human genome into classes related to regulatory attributes, such as epigenetic modifications and transcription factor binding sites, in addition to classes related to gene structure and cross-species conservation. Using the distribution of derived allele frequencies (DAF) within each class, we assessed the strength of natural selection for each class relative to the genome as a whole. We applied this DAF analysis to Perlegen resequenced SNPs genome-wide. Regulatory elements annotated by Ensembl such as specific histone methylation sites as well as classes defined by cross-species conservation showed negative selection in comparison to the genome as a whole.ConclusionsThese results highlight which annotated classes are under purifying selection, have putative functional importance, and contain SNPs that are strong candidates for follow-up studies after genome-wide association. Such SNP annotation may also be useful in interpreting results of whole-genome sequencing studies.

Ongoing GC-Biased Evolution Is Widespread in the Human Genome and Enriched Near Recombination Hot Spots

Fast evolving regions of many metazoan genomes show a bias toward substitutions that change weak (A,T) into strong (G,C) base pairs. Single-nucleotide polymorphisms (SNPs) do not share this pattern, suggesting that it results from biased fixation rather than biased mutation. Supporting this hypothesis, analyses of polymorphism in specific regions of the human genome have identified a positive correlation between weak to strong (W→S) SNPs and derived allele frequency (DAF), suggesting that SNPs become increasingly GC biased over time, especially in regions of high recombination. Using polymorphism data generated by the 1000 Genomes Project from 179 individuals from 4 human populations, we evaluated the extent and distribution of ongoing GC-biased evolution in the human genome. We quantified GC fixation bias by comparing the DAFs of W→S mutations and S→W mutations using a Mann–Whitney U test. Genome-wide, W→S SNPs have significantly higher DAFs than S→W SNPs. This pattern is widespread across the human genome but varies in magnitude along the chromosomes. We found extreme GC-biased evolution in neighborhoods of recombination hot spots, a significant correlation between GC bias and recombination rate, and an inverse correlation between GC bias and chromosome arm length. These findings demonstrate the presence of ongoing fixation bias favoring G and C alleles throughout the human genome and suggest that the bias is caused by a recombination-associated process, such as GC-biased gene conversion.

GC-Biased Evolution Near Human Accelerated Regions

Regions of the genome that have been the target of positive selection specifically along the human lineage are of special importance in human biology. We used high throughput sequencing combined with methods to enrich human genomic samples for particular targets to obtain the sequence of 22 chromosomal samples at high depth in 40 kb neighborhoods of 49 previously identified 100–400 bp elements that show evidence for human accelerated evolution. In addition to selection, the pattern of nucleotide substitutions in several of these elements suggested an historical bias favoring the conversion of weak (A or T) alleles into strong (G or C) alleles. Here we found strong evidence in the derived allele frequency spectra of many of these 40 kb regions for ongoing weak-to-strong fixation bias. Comparison of the nucleotide composition at polymorphic loci to the composition at sites of fixed substitutions additionally reveals the signature of historical weak-to-strong fixation bias in a subset of these regions. Most of the regions with evidence for historical bias do not also have signatures of ongoing bias, suggesting that the evolutionary forces generating weak-to-strong bias are not constant over time. To investigate the role of selection in shaping these regions, we analyzed the spatial pattern of polymorphism in our samples. We found no significant evidence for selective sweeps, possibly because the signal of such sweeps has decayed beyond the power of our tests to detect them. Together, these results do not rule out functional roles for the observed changes in these regions—indeed there is good evidence that the first two are functional elements in humans—but they suggest that a fixation process (such as biased gene conversion) that is biased at the nucleotide level, but is otherwise selectively neutral, could be an important evolutionary force at play in them, both historically and at present.

Single-nucleotide evolutionary constraint scores highlight disease-causing mutations

Single-nucleotide evolutionary constraint scores highlight disease-causing mutations

Positive selection drives population differentiation in the skeletal genes in modern humans

During the course of evolution, the human skeletal system has evolved rapidly leading to an incredible array of phenotypic diversity, including variations in height and bone mineral density. However, the genetic basis of this phenotypic diversity and the relatively rapid tempo of evolution have remained largely undocumented. Here, we discover that skeletal genes exhibit a significantly greater level of population differentiation among humans compared with other genes in the genome. The pattern is exceptionally evident at amino acid-altering sites within these genes. Divergence is greater between Africans and both Europeans and East Asians. In contrast, relatively weak differentiation is observed between Europeans and East Asians. SNPs with higher levels of differentiation have correspondingly higher derived allele frequencies in Europeans and East Asians. Thus, it appears that positive selection has operated on skeletal genes in the non-African populations and this may have been initiated with the human colonization of Eurasia. In conclusion, we provide genetic evidence supporting the rapid evolution of the human skeletal system and the associated diversity of phenotypes.

Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken

BackgroundThe chicken (Gallus gallus), like most avian species, has a very distinct karyotype consisting of many micro- and a few macrochromosomes. While it is known that recombination frequencies are much higher for micro- as compared to macrochromosomes, there is limited information on differences in linkage disequilibrium (LD) and haplotype diversity between these two classes of chromosomes. In this study, LD and haplotype diversity were systematically characterized in 371 birds from eight chicken populations (commercial lines, fancy breeds, and red jungle fowl) across macro- and microchromosomes. To this end we sampled four regions of ~1 cM each on macrochromosomes (GGA1 and GGA2), and four 1.5 -2 cM regions on microchromosomes (GGA26 and GGA27) at a high density of 1 SNP every 2 kb (total of 889 SNPs).ResultsAt a similar physical distance, LD, haplotype homozygosity, haploblock structure, and haplotype sharing were all lower for the micro- as compared to the macrochromosomes. These differences were consistent across populations. Heterozygosity, genetic differentiation, and derived allele frequencies were also higher for the microchromosomes. Differences in LD, haplotype variation, and haplotype sharing between populations were largely in line with known demographic history of the commercial chicken. Despite very low levels of LD, as measured by r2 for most populations, some haploblock structure was observed, particularly in the macrochromosomes, but the haploblock sizes were typically less than 10 kb.ConclusionDifferences in LD between micro- and macrochromosomes were almost completely explained by differences in recombination rate. Differences in haplotype diversity and haplotype sharing between micro- and macrochromosomes were explained by differences in recombination rate and genotype variation. Haploblock structure was consistent with demography of the chicken populations, and differences in recombination rates between micro- and macrochromosomes. The limited haploblock structure and LD suggests that future whole-genome marker assays will need 100+K SNPs to exploit haplotype information. Interpretation and transferability of genetic parameters will need to take into account the size of chromosomes in chicken, and, since most birds have microchromosomes, in other avian species as well.

The Neandertal genome and ancient DNA authenticity

Recent advances in high-thoughput DNA sequencing have made genome-scale analyses of genomes of extinct organisms possible. With these new opportunities come new difficulties in assessing the authenticity of the DNA sequences retrieved. We discuss how these difficulties can be addressed, particularly with regard to analyses of the Neandertal genome. We argue that only direct assays of DNA sequence positions in which Neandertals differ from all contemporary humans can serve as a reliable means to estimate human contamination. Indirect measures, such as the extent of DNA fragmentation, nucleotide misincorporations, or comparison of derived allele frequencies in different fragment size classes, are unreliable. Fortunately, interim approaches based on mtDNA differences between Neandertals and current humans, detection of male contamination through Y chromosomal sequences, and repeated sequencing from the same fossil to detect autosomal contamination allow initial large-scale sequencing of Neandertal genomes. This will result in the discovery of fixed differences in the nuclear genome between Neandertals and current humans that can serve as future direct assays for contamination. For analyses of other fossil hominins, which may become possible in the future, we suggest a similar ‘boot-strap' approach in which interim approaches are applied until sufficient data for more definitive direct assays are acquired.

SNP deserts of Asian cultivated rice: genomic regions under domestication

When performing a genome-wide comparison between indica (93-11) and japonica (Nipponbare), we find 8% of the genome, which have an extremely low SNP rate (< 1 SNP/kb). Inside these 'SNP deserts', experimentally confirmed genes show increased K(a)/K(s) that indicate adaptive selection. To further elucidate this connection, we survey the level and pattern of genetic variation in both cultivated and wild rice groups, using 155 noncoding regions located within SNP deserts. The results suggest that cultivated rice has greatly reduced genetic variation within SNP deserts as compared to either the nondesert or corresponding genomic regions in wild rice. Consistent with this reduction in genetic variation, we find a biased distribution of derived allele frequency in the cultivated group, indicative of positive selection. Furthermore, over half of the confirmed, domestication-related genes are found within SNP deserts, also suggesting that SNP deserts are strongly related to domestication, and might be the key sites in the process of domestication.