Genetic continuity across a deeply divergent linguistic contact zone in North Maluku, Indonesia.

BACKGROUND: The islands of North Maluku, Indonesia occupy a central position in the major prehistoric dispersal streams that shaped the peoples of Island Southeast Asia and the Pacific. Within this region a linguistic contact zone exists where speakers of Papuan and Austronesian languages reside in close proximity. Here we use population genetic data to assess the extent to which North Maluku populations experienced admixture of Asian genetic material, and whether linguistic boundaries reflect genetic differentiation today. RESULTS: Autosomal and X-linked markers reveal overall Asian admixture of 67% in North Maluku, demonstrating a substantial contribution of genetic material into the region from Asia. We observe no evidence of population structure associated with ethnicity or language affiliation. CONCLUSIONS: Our data support a model of widespread Asian admixture in North Maluku, likely mediated by the expansion of Austronesian-speaking peoples into the region during the mid Holocene. In North Maluku there is no genetic differentiation in terms of Austronesian- versus Papuan-speakers, suggesting extensive gene flow across linguistic boundaries. In a regional context, our results illuminate a major genetic divide at the Molucca Sea, between the islands of Sulawesi and North Maluku. West of this divide, populations exhibit predominantly Asian ancestry, with very little contribution of Papuan genetic material. East of the Molucca Sea, populations show diminished rates of Asian admixture and substantial persistence of Papuan genetic diversity.

The "Frankenplasmid" lab: an investigative exercise for teaching recombinant DNA methods.

We describe an investigative laboratory module designed to give college undergraduates strong practical and theoretical experience with recombinant DNA methods within 3 weeks. After deducing restriction enzyme maps for two different plasmids, students ligate the plasmids together in the same reaction, transform E. coli with this mixture of ligated DNA, and plate the cells on media that specifically select for hybrid plasmids. The main goal of the assignment is for students to deduce the gene map of one hybrid "Frankenplasmid" using the LacZ phenotype of its transformants, PCR, and restriction mapping. Our protocol results in a number of possible outcomes, meaning that students are mapping truly unknown plasmids. The open-ended nature of this assignment results in an effective module that teaches recombinant DNA procedures while engaging students with its investigative approach, increasing complexity, and puzzle-like quality. Moreover, the modular design of the activity allows it to be adapted to a more limited schedule, introductory courses, or more advanced courses.

Distribution and mechanism of α-amanitin tolerance in mycophagous Drosophila (Diptera: Drosophilidae).

Many mycophagous species of Drosophila can tolerate the mushroom poison α-amanitin in wild mushrooms and in artificial diet. We conducted feeding assays with sixteen Drosophila species and α-amanitin in artificial diet to better determine the phylogenetic distribution of this tolerance. For eight tolerant and one related susceptible species, we sequenced the gene encoding the large subunit of RNA Polymerase II, which is the target site of α-amanitin. We found no differences in the gene that could account for differences in susceptibility to the toxin. We also conducted feeding assays in which α-amanitin was combined with chemical inhibitors of cytochrome P450s or glutathione S-transferases (GSTs) in artificial diet to determine if either of these enzyme families is involved in tolerance to α-amanitin. We found that an inhibitor of GSTs did not reduce tolerance to α-amanitin, but that an inhibitor of cytochrome P450s reduced tolerance in several species. It is possible that the same cytochrome P450 activity that produces tolerance of α-amanitin might produce tolerance of other mushroom toxins as well. If so, a general detoxification mechanism based on cytochrome P450s might answer the question of how tolerance to α-amanitin arose in mycophagous Drosophila when this toxin is found in relatively few mushrooms.

Molecular evolution of GYPC: evidence for recent structural innovation and positive selection in humans.

GYPC encodes two erythrocyte surface sialoglycoproteins in humans, glycophorin C and glycophorin D (GPC and GPD), via initiation of translation at two start codons on a single transcript. The malaria-causing parasite Plasmodium falciparum uses GPC as a means of invasion into the human red blood cell. Here, we examine the molecular evolution of GYPC among the Hominoidea (Greater and Lesser Apes) and also the pattern of polymorphism at the locus in a global human sample. We find an excess of nonsynonymous divergence among species that appears to be caused solely by accelerated evolution of GYPC in the human lineage. Moreover, we find that the ability of GYPC to encode both GPC and GPD is a uniquely human trait, caused by the evolution of the GPC start codon in the human lineage. The pattern of polymorphism among humans is consistent with a hitchhiking event at the locus, suggesting that positive natural selection affected GYPC in the relatively recent past. Because GPC is exploited by P. falciparum for invasion of the red blood cell, we hypothesize that selection for evasion of P. falciparum has caused accelerated evolution of GYPC in humans (relative to other primates) and that this positive selection has continued to act in the recent evolution of our species. These data suggest that malaria has played a powerful role in shaping molecules on the surface of the human red blood cell. In addition, our examination of GYPC reveals a novel mechanism of protein evolution: co-option of untranslated region (UTR) sequence following the formation of a new start codon. In the case of human GYPC, the ancestral protein (GPD) continues to be produced through leaky translation. Because leaky translation is a widespread phenomenon among genes and organisms, we suggest that co-option of UTR sequence may be an important source of protein innovation.

Molecular population genetics of SLC4A1 and Southeast Asian ovalocytosis.

Southeast Asian ovalocytosis (SAO) is an erythrocyte abnormality that protects affected individuals from cerebral malaria. This trait is caused by a 27-bp deletion in the SLC4A1 gene, which is lethal when homozygous. We reseqeunced approximately 5 kb of SLC4A1 in an Indonesian population where SAO is prevalent to better understand the evolution of this clinically important trait. The four SAO chromosomes we resequenced share a single haplotype that differs from a sampled non-SAO haplotype only by the 27-bp deletion. Comparison of Indonesian sequence data to that from two other Asian populations (aboriginal Taiwanese and Japanese) shows Indonesian SLC4A1 to be strongly differentiated from the Taiwanese, but not the Japanese. Indeed, the Taiwanese sample contains only chromosomes that are highly divergent from all sampled SAO chromosomes. Because earlier studies have found an association between Austronesian-speakers (who most likely originated in Taiwan) and SAO, our failure to find SAO-like chromosomes in Taiwan is unexpected. Finally, our data find a strong excess of high-frequency derived alleles in all three populations. These alleles include the non-synonymous 'Memphis' variant, which is known to affect anion transport across the erythrocyte membrane. Our data suggest a role for recent natural selection acting on Memphis or a linked variant.

Contrasting signatures of population growth for mitochondrial DNA and Y chromosomes among human populations in Africa.

A history of Pleistocene population expansion has been inferred from the frequency spectrum of polymorphism in the mitochondrial DNA (mtDNA) of many human populations. Similar patterns are not typically observed for autosomal and X-linked loci. One explanation for this discrepancy is a recent population bottleneck, with different rates of recovery for haploid and autosomal loci as a result of their different effective population sizes. This hypothesis predicts that mitochondrial and Y chromosomal DNA will show a similar skew in the frequency spectrum in populations that have experienced a recent increase in effective population size. We test this hypothesis by resequencing 6.6 kb of noncoding Y chromosomal DNA and 780 basepairs of the mtDNA cytochrome c oxidase subunit III (COIII) gene in 172 males from 5 African populations. Four tests of population expansion are employed for each locus in each population: Fu's Fs statistic, the R(2) statistic, coalescent simulations, and the mismatch distribution. Consistent with previous results, patterns of mtDNA polymorphism better fit a model of constant population size for food-gathering populations and a model of population expansion for food-producing populations. In contrast, none of the tests reveal evidence of Y chromosome growth for either food-gatherers or food-producers. The distinct mtDNA and Y chromosome polymorphism patterns most likely reflect sex-biased demographic processes in the recent history of African populations. We hypothesize that males experienced smaller effective population sizes and/or lower rates of migration during the Bantu expansion, which occurred over the last 5,000 years.

Inferring human population sizes, divergence times and rates of gene flow from mitochondrial, X and Y chromosome resequencing data.

We estimate parameters of a general isolation-with-migration model using resequence data from mitochondrial DNA (mtDNA), the Y chromosome, and two loci on the X chromosome in samples of 25-50 individuals from each of 10 human populations. Application of a coalescent-based Markov chain Monte Carlo technique allows simultaneous inference of divergence times, rates of gene flow, as well as changes in effective population size. Results from comparisons between sub-Saharan African and Eurasian populations estimate that 1500 individuals founded the ancestral Eurasian population approximately 40 thousand years ago (KYA). Furthermore, these small Eurasian founding populations appear to have grown much more dramatically than either African or Oceanian populations. Analyses of sub-Saharan African populations provide little evidence for a history of population bottlenecks and suggest that they began diverging from one another upward of 50 KYA. We surmise that ancestral African populations had already been geographically structured prior to the founding of ancestral Eurasian populations. African populations are shown to experience low levels of mitochondrial DNA gene flow, but high levels of Y chromosome gene flow. In particular, Y chromosome gene flow appears to be asymmetric, i.e., from the Bantu-speaking population into other African populations. Conversely, mitochondrial gene flow is more extensive between non-African populations, but appears to be absent between European and Asian populations.

Ascertainment bias and the pattern of nucleotide diversity at the human ALDH2 locus in a Japanese population.

Many East Asian human populations harbor a high-frequency deficiency allele for the aldehyde dehydrogenase 2 (ALDH2) enzyme, a critical protein involved in the metabolism of ethanol. Here we use resequencing and long-range SNP haplotype data from a Japanese sample to test whether patterns of nucleotide diversity and linkage disequilibrium at this locus are compatible with a standard neutral model of evolution. Examination of the pattern of polymorphism at a locus such as this, where the frequency of a common allele is known a priori, introduces an ascertainment bias that must be corrected for in analyses of the frequency spectrum of polymorphisms. We apply a flexible and generally applicable simulation approach to correct for this bias in our ALDH2 data and, also, to explore the effect of bias on the commonly used summary statistics Tajima's D, Fu and Li's D, and Fay and Wu's H. Our study finds no evidence that the pattern of genetic variation at ALDH2 differs from that expected under a standard neutral model. However, our general examination of ascertainment bias indicates that a priori knowledge of segregating alleles greatly affects the expected distributions of summary statistics. Under many parameter combinations we find that ascertainment bias introduces an elevated rate of false positives when summary statistics are used to test for deviations from a standard neutral model. However, we also show that over a wide range of conditions the power of all summary statistics can be greatly increased by incorporating prior knowledge of segregating alleles.

Deep haplotype divergence and long-range linkage disequilibrium at xp21.1 provide evidence that humans descend from a structured ancestral population.

Fossil evidence links human ancestry with populations that evolved from modern gracile morphology in Africa 130,000-160,000 years ago. Yet fossils alone do not provide clear answers to the question of whether the ancestors of all modern Homo sapiens comprised a single African population or an amalgamation of distinct archaic populations. DNA sequence data have consistently supported a single-origin model in which anatomically modern Africans expanded and completely replaced all other archaic hominin populations. Aided by a novel experimental design, we present the first genetic evidence that statistically rejects the null hypothesis that our species descends from a single, historically panmictic population. In a global sample of 42 X chromosomes, two African individuals carry a lineage of noncoding 17.5-kb sequence that has survived for >1 million years without any clear traces of ongoing recombination with other lineages at this locus. These patterns of deep haplotype divergence and long-range linkage disequilibrium are best explained by a prolonged period of ancestral population subdivision followed by relatively recent interbreeding. This inference supports human evolution models that incorporate admixture between divergent African branches of the genus Homo.

Evidence for archaic Asian ancestry on the human X chromosome.

The human RRM2P4 pseudogene has a pattern of nucleotide polymorphism that is unlike any locus published to date. A gene tree constructed from a 2.4-kb fragment of the RRM2P4 locus sequenced in a sample of 41 worldwide humans clearly roots in East Asia and has a most-recent common ancestor approximately 2 Myr before present. The presence of this basal lineage exclusively in Asia results in higher nucleotide diversity among non-Africans than among Africans. A global survey of a single-nucleotide polymorphism that is diagnostic for the basal, Asian lineage in 570 individuals shows that it occurs at frequencies up to 53% in south China, whereas only one of 177 surveyed Africans carries this archaic lineage. We suggest that this ancient lineage is a remnant of introgressive hybridization between expanding anatomically modern humans emerging from Africa and archaic populations in Eurasia.

Global patterns of human mitochondrial DNA and Y-chromosome structure are not influenced by higher migration rates of females versus males.

Global-scale patterns of human population structure may be influenced by the rate of migration among populations that is nearly eight times higher for females than for males. This difference is attributed mainly to the widespread practice of patrilocality, in which women move into their mates' residences after marriage. Here we directly test this hypothesis by comparing global patterns of DNA sequence variation on the Y chromosome and mitochondrial DNA (mtDNA) in the same panel of 389 individuals from ten populations (four from Africa and two each from Europe, Asia and Oceania). We introduce a new strategy to assay Y-chromosome variation that identifies a high density of single-nucleotide polymorphisms, allows complete sequencing of all individuals rather than relying on predetermined markers and provides direct sequence comparisons with mtDNA. We found the overall proportion of between-group variation (Phi(ST)) to be 0.334 for the Y chromosome and 0.382 for mtDNA. Genetic differentiation between populations was similar for the Y chromosome and mtDNA at all geographic scales that we tested. Although patrilocality may be important at the local scale, patterns of genetic structure on the continental and global scales are not shaped by the higher rate of migration among females than among males.

Heterogeneous patterns of variation among multiple human x-linked Loci: the possible role of diversity-reducing selection in non-africans.

Studies of human DNA sequence polymorphism reveal a range of diversity patterns throughout the genome. This variation among loci may be due to natural selection, demographic influences, and/or different sampling strategies. Here we build on a continuing study of noncoding regions on the X chromosome in a panel of 41 globally sampled humans representing African and non-African populations by examining patterns of DNA sequence variation at four loci (APXL, AMELX, TNFSF5, and RRM2P4) and comparing these patterns with those previously reported at six loci in the same panel of 41 individuals. We also include comparisons with patterns of noncoding variation seen at five additional X-linked loci that were sequenced in similar global panels. We find that, while almost all loci show a reduction in non-African diversity, the magnitude of the reduction varies substantially across loci. The large observed variance in non-African levels of diversity results in the rejection of a neutral model of molecular evolution with a multi-locus HKA test under both a constant size and a bottleneck model. In non-Africans, some loci harbor an excess of rare mutations over neutral equilibrium predictions, while other loci show no such deviation in the distribution of mutation frequencies. We also observe a positive relationship between recombination rate and frequency spectra in our non-African, but not in our African, sample. These results indicate that a simple out-of-Africa bottleneck model is not sufficient to explain the observed patterns of sequence variation and that diversity-reducing selection acting at a subset of loci and/or a more complex neutral model must be invoked.

Genetic evidence for unequal effective population sizes of human females and males.

The time to the most recent common ancestor (TMRCA) of the human mitochondria (mtDNA) is estimated to be older than that of the nonrecombining portion of the Y chromosome (NRY). Surveys of variation in globally distributed humans typically result in mtDNA TMRCA values just under 200 thousand years ago (kya), whereas those for the NRY range between 46 and 110 kya. A favored hypothesis for this finding is that natural selection has acted on the NRY, leading to a recent selective sweep. An alternate hypothesis is that sex-biased demographic processes are responsible. Here, we re-examine the disparity between NRY and mtDNA TMRCAs using data collected from individual human populations--a sampling strategy that minimizes the confounding influence of population subdivision in global data sets. We survey variation at 782 bp of the mitochondrial cytochrome c oxidase subunit 3 gene as well as at 26.5 kb of noncoding DNA from the NRY in a sample of 25 Khoisan, 24 Mongolians, and 24 Papua New Guineans. Data from both loci in all populations are best described by a model of constant population size, with the exception of Mongolian mtDNA, which appears to be experiencing rapid population growth. Taking these demographic models into account, we estimate the TMRCAs for each locus in each population. A pattern that is remarkably consistent across all three populations is an approximately twofold deeper coalescence for mtDNA than for the NRY. The oldest TMRCAs are observed for the Khoisan (73.6 kya for the NRY and 176.5 kya for mtDNA), whereas those in the non-African populations are consistently lower (averaging 47.7 kya for the NRY and 92.8 kya for mtDNA). Our data do not suggest that differential natural selection is the cause of this difference in TMRCAs. Rather, these results are most consistent with a higher female effective population size.

European ACP1*C allele has recessive deleterious effects on early life viability.

The acid phosphatase locus (ACP1) is a classical polymorphism that has been surveyed in hundreds of human populations worldwide. Among individuals of European ancestry, the ACP1*C allele occurs with an average frequency of approximately 0.05, whereas it is nearly absent in all other human populations. It has been hypothesized that this allele is maintained by overdominant selection among European populations. Here, we analyze ACP1 protein polymorphism data from more than 50,000 individuals previously surveyed in 67 populations across Europe as well as inheritance data from more than 6,000 European parent-offspring pairs to assess the signature of natural selection currently acting on this allele. Although we see a significant excess of ACP1*C heterozygotes relative to Hardy-Weinberg expectations, we find no evidence that natural selection favors ACP1*C heterozygotes. Instead, ACP1*C appears to have a strongly deleterious and recessive fitness effect. We observed only 48.9% of expected homozygous offspring from heterozygous parents and significantly fewer homozygotes than expected within populations. Because parent-offspring pairs indicate a significant deficiency of ACP1*C homozygotes, we infer that viability selection is acting on ACP1*C homozygotes very early in life, perhaps before birth. We estimate that approximately 1.2% of all couples of European ancestry are composed of individuals who both carry the APC1*C allele. As such, selection against ACP1*C homozygotes may represent a nonnegligible contribution to the overall number of spontaneous abortions among women of European ancestry and may cause substantial fertility reductions among some combinations of parental genotypes.

Human population structure and its effects on sampling Y chromosome sequence variation.

The excess of rare variants in global sequencing studies of the nonrecombining portion of the Y chromosome (NRY) has been interpreted as evidence for the effects of human demographic expansion. However, many NRY polymorphisms are geographically localized and the effect of different geographical sampling on patterns of NRY variation is unknown. We use two sampling designs to detect population structure and its effects on patterns of human NRY polymorphism. First, we sequence 26.5 kb of noncoding Y chromosome DNA from 92 globally distributed males representing 35 populations. We find that the number of polymorphisms with singleton variants is positively correlated with the number of populations sampled and that there is a significant negative correlation of Tajima's D (TD) and Fu and Li's D (FD) statistics with the number of pooled populations. We then sequence the same region in a total of 73 males sampled from 3 distinct populations and find that TD and FD values for the 3 pooled and individual population samples were much less negative than those in the aforementioned global sample. Coalescent simulations show that a simple splitting model of population structure, with no changes in population size, is sufficient to produce the negative values of TD seen in our pooled samples. These empirical and simulation results suggest that observed levels of NRY population structure may lead to an upward bias in the number of singleton variants in global surveys and call into question inferences of population expansion based on global sampling strategies.