Joseph C. Watkins

Genetic evidence for archaic admixture in Africa.

A long-debated question concerns the fate of archaic forms of the genus Homo: did they go extinct without interbreeding with anatomically modern humans, or are their genes present in contemporary populations? This question is typically focused on the genetic contribution of archaic forms outside of Africa. Here we use DNA sequence data gathered from 61 noncoding autosomal regions in a sample of three sub-Saharan African populations (Mandenka, Biaka, and San) to test models of African archaic admixture. We use two complementary approximate-likelihood approaches and a model of human evolution that involves recent population structure, with and without gene flow from an archaic population. Extensive simulation results reject the null model of no admixture and allow us to infer that contemporary African populations contain a small proportion of genetic material (≈2%) that introgressed ≈35 kya from an archaic population that split from the ancestors of anatomically modern humans ≈700 kya. Three candidate regions showing deep haplotype divergence, unusual patterns of linkage disequilibrium, and small basal clade size are identified and the distributions of introgressive haplotypes surveyed in a sample of populations from across sub-Saharan Africa. One candidate locus with an unusual segment of DNA that extends for >31 kb on chromosome 4 seems to have introgressed into modern Africans from a now-extinct taxon that may have lived in central Africa. Taken together our results suggest that polymorphisms present in extant populations introgressed via relatively recent interbreeding with hominin forms that diverged from the ancestors of modern humans in the Lower-Middle Pleistocene.

A Haplotype at STAT2 Introgressed from Neanderthals and Serves as a Candidate of Positive Selection in Papua New Guinea

Signals of archaic admixture have been identified through comparisons of the draft Neanderthal and Denisova genomes with those of living humans. Studies of individual loci contributing to these genome-wide average signals are required for characterization of the introgression process and investigation of whether archaic variants conferred an adaptive advantage to the ancestors of contemporary human populations. However, no definitive case of adaptive introgression has yet been described. Here we provide a DNA sequence analysis of the innate immune gene STAT2 and show that a haplotype carried by many Eurasians (but not sub-Saharan Africans) has a sequence that closely matches that of the Neanderthal STAT2. This haplotype, referred to as N, was discovered through a resequencing survey of the entire coding region of STAT2 in a global sample of 90 individuals. Analyses of publicly available complete genome sequence data show that haplotype N shares a recent common ancestor with the Neanderthal sequence (~80 thousand years ago) and is found throughout Eurasia at an average frequency of ~5%. Interestingly, N is found in Melanesian populations at ~10-fold higher frequency (~54%) than in Eurasian populations. A neutrality test that controls for demography rejects the hypothesis that a variant of N rose to high frequency in Melanesia by genetic drift alone. Although we are not able to pinpoint the precise target of positive selection, we identify nonsynonymous mutations in ERBB3, ESYT1, and STAT2-all of which are part of the same 250 kb introgressive haplotype-as good candidates.

The ratio of human X chromosome to autosome diversity is positively correlated with genetic distance from genes

The ratio of X-linked to autosomal diversity was estimated from an analysis of six human genome sequences and found to deviate from the expected value of 0.75. However, the direction of this deviation depends on whether a particular sequence is close to or far from the nearest gene. This pattern may be explained by stronger locally acting selection on X-linked genes compared with autosomal genes, combined with larger effective population sizes for females than for males.

Balinese Y-Chromosome Perspective on the Peopling of Indonesia: Genetic Contributions from Pre-Neolithic Hunter-Gatherers, Austronesian Farmers, and Indian Traders

The island of Bali lies near the center of the southern chain of islands in the Indonesian archipelago, which served as a stepping-stone for early migrations of hunter-gatherers to Melanesia and Australia and for more recent migrations of Austronesian farmers from mainland Southeast Asia to the Pacific. Bali is the only Indonesian island with a population that currently practices the Hindu religion and preserves various other Indian cultural, linguistic, and artistic traditions (Lansing 1983). Here, we examine genetic variation on the Y chromosomes of 551 Balinese men to investigate the relative contributions of Austronesian farmers and pre-Neolithic hunter-gatherers to the contemporary Balinese paternal gene pool and to test the hypothesis of recent paternal gene flow from the Indian subcontinent. Seventy-one Y-chromosome binary polymorphisms (single nucleotide polymorphisms, SNPs) and 10 Y-chromosome-linked short tandem repeats (STRs) were genotyped on a sample of 1,989 Y chromosomes from 20 populations representing Indonesia (including Bali), southern China, Southeast Asia, South Asia, the Near East, and Oceania. SNP genotyping revealed 22 Balinese lineages, 3 of which (O-M95, O-M119, and O-M122) account for nearly 83.7% of Balinese Y chromosomes. Phylogeographic analyses suggest that all three major Y-chromosome haplogroups migrated to Bali with the arrival of Austronesian speakers; however, STR diversity patterns associated with these haplogroups are complex and may be explained by multiple waves of Austronesian expansion to Indonesia by different routes. Approximately 2.2% of contemporary Balinese Y chromosomes (i.e., K-M9*, K-M230, and M lineages) may represent the pre-Neolithic component of the Indonesian paternal gene pool. In contrast, eight other haplogroups (e.g., within H, J, L, and R), making up approximately 12% of the Balinese paternal gene pool, appear to have migrated to Bali from India. These results indicate that the Austronesian expansion had a profound effect on the composition of the Balinese paternal gene pool and that cultural transmission from India to Bali was accompanied by substantial levels of gene flow.

Coevolution of languages and genes on the island of Sumba, eastern Indonesia

Numerous studies indicate strong associations between languages and genes among human populations at the global scale, but all broader scale genetic and linguistic patterns must arise from processes originating at the community level. We examine linguistic and genetic variation in a contact zone on the eastern Indonesian island of Sumba, where Neolithic Austronesian farming communities settled and began interacting with aboriginal foraging societies ≈3,500 years ago. Phylogenetic reconstruction based on a 200-word Swadesh list sampled from 29 localities supports the hypothesis that Sumbanese languages derive from a single ancestral Austronesian language. However, the proportion of cognates (words with a common origin) traceable to Proto-Austronesian (PAn) varies among language subgroups distributed across the island. Interestingly, a positive correlation was found between the percentage of Y chromosome lineages that derive from Austronesian (as opposed to aboriginal) ancestors and the retention of PAn cognates. We also find a striking correlation between the percentage of PAn cognates and geographic distance from the site where many Sumbanese believe their ancestors arrived on the island. These language–gene–geography correlations, unprecedented at such a fine scale, imply that historical patterns of social interaction between expanding farmers and resident hunter-gatherers largely explain community-level language evolution on Sumba. We propose a model to explain linguistic and demographic coevolution at fine spatial and temporal scales.

Neandertal Origin of Genetic Variation at the Cluster of OAS Immunity Genes

Analyses of ancient DNA from extinct humans reveal signals of at least two independent hybridization events in the history of non-African populations. To date, there are very few examples of specific genetic variants that have been rigorously identified as introgressive. Here, we survey DNA sequence variation in the OAS gene cluster on chromosome 12 and provide strong evidence that a haplotype extending for ~185 kb introgressed from Neandertals. This haplotype is nearly restricted to Eurasians and is estimated to have diverged from the Neandertal sequence ~125 kya. Despite the potential for novel functional variation, the observed frequency of this haplotype is consistent with neutral introgression. This is the second locus in the human genome, after STAT2, carrying distinct haplotypes that appear to have introgressed separately from both Neandertals and Denisova.

Evidence for Increased Levels of Positive and Negative Selection on the X Chromosome versus Autosomes in Humans

Partially recessive variants under positive selection are expected to go to fixation more quickly on the X chromosome as a result of hemizygosity, an effect known as faster-X. Conversely, purifying selection is expected to reduce substitution rates more effectively on the X chromosome. Previous work in humans contrasted divergence on the autosomes and X chromosome, with results tending to support the faster-X effect. However, no study has yet incorporated both divergence and polymorphism to quantify the effects of both purifying and positive selection, which are opposing forces with respect to divergence. In this study, we develop a framework that integrates previously developed theory addressing differential rates of X and autosomal evolution with methods that jointly estimate the level of purifying and positive selection via modeling of the distribution of fitness effects (DFE). We then utilize this framework to estimate the proportion of nonsynonymous substitutions fixed by positive selection (α) using exome sequence data from a West African population. We find that varying the female to male breeding ratio (β) has minimal impact on the DFE for the X chromosome, especially when compared with the effect of varying the dominance coefficient of deleterious alleles (h). Estimates of α range from 46% to 51% and from 4% to 24% for the X chromosome and autosomes, respectively. While dependent on h, the magnitude of the difference between α values estimated for these two systems is highly statistically significant over a range of biologically realistic parameter values, suggesting faster-X has been operating in humans.

Male dominance rarely skews the frequency distribution of Y chromosome haplotypes in human populations.

A central tenet of evolutionary social science holds that behaviors, such as those associated with social dominance, produce fitness effects that are subject to cultural selection. However, evidence for such selection is inconclusive because it is based on short-term statistical associations between behavior and fertility. Here, we show that the evolutionary effects of dominance at the population level can be detected using noncoding regions of DNA. Highly variable polymorphisms on the nonrecombining portion of the Y chromosome can be used to trace lines of descent from a common male ancestor. Thus, it is possible to test for the persistence of differential fertility among patrilines. We examine haplotype distributions defined by 12 short tandem repeats in a sample of 1269 men from 41 Indonesian communities and test for departures from neutral mutation-drift equilibrium based on the Ewens sampling formula. Our tests reject the neutral model in only 5 communities. Analysis and simulations show that we have sufficient power to detect such departures under varying demographic conditions, including founder effects, bottlenecks, and migration, and at varying levels of social dominance. We conclude that patrilines seldom are dominant for more than a few generations, and thus traits or behaviors that are strictly paternally inherited are unlikely to be under strong cultural selection.

A de novo missense mutation of GABRB2 causes early myoclonic encephalopathy

Background Early myoclonic encephalopathy (EME), a disease with a devastating prognosis, is characterised by neonatal onset of seizures and massive myoclonus accompanied by a continuous suppression-burst EEG pattern. Three genes are associated with EMEs that have metabolic features. Here, we report a pathogenic mutation of an ion channel as a cause of EME for the first time. Methods Sequencing was performed for 214 patients with epileptic seizures using a gene panel with 109 genes that are known or suspected to cause epileptic seizures. Functional assessments were demonstrated by using electrophysiological experiments and immunostaining for mutant γ-aminobutyric acid-A (GABAA) receptor subunits in HEK293T cells. Results We discovered a de novo heterozygous missense mutation (c.859A>C [p.Thr287Pro]) in the GABRB2-encoded β2 subunit of the GABAA receptor in an infant with EME. No GABRB2 mutations were found in three other EME cases or in 166 patients with infantile spasms. GABAA receptors bearing the mutant β2 subunit were poorly trafficked to the cell membrane and prevented γ2 subunits from trafficking to the cell surface. The peak amplitudes of currents from GABAA receptors containing only mutant β2 subunits were smaller than that of those from receptors containing only wild-type β2 subunits. The decrease in peak current amplitude (96.4% reduction) associated with the mutant GABAA receptor was greater than expected, based on the degree to which cell surface expression was reduced (66% reduction). Conclusion This mutation has complex functional effects on GABAA receptors, including reduction of cell surface expression and attenuation of channel function, which would significantly perturb GABAergic inhibition in the brain.

Clinical implications of SCN1A missense and truncation variants in a large Japanese cohort with Dravet syndrome

Two major classes of SCN1A variants are associated with Dravet syndrome (DS): those that result in haploinsufficiency (truncating) and those that result in an amino acid substitution (missense). The aim of this retrospective study was to describe the first large cohort of Japanese patients with SCN1A mutation–positive DS (n = 285), and investigate the relationship between variant (type and position) and clinical expression and response to treatment.