Papers by Daniel Garrigan
To gain an understanding of the evolution and genomic context of avian major histocompatibility c... more To gain an understanding of the evolution and genomic context of avian major histocompatibility complex (Mhc) genes, we sequenced a 38.8-kb Mhc-bearing cosmid insert from a red-winged blackbird (Agelaius phoeniceus). The DNA sequence, the longest yet retrieved from a bird other than a chicken, provides a detailed view of the process of gene duplication, divergence, and degeneration (''birth and death'') in the avian Mhc, as well as a glimpse into major noncoding features of a songbird genome. The peptide-binding region (PBR) of the single Mhc class II B gene in this region, Agph-DAB2, is almost devoid of polymorphism, and a still-segregating single-base-pair deletion and other features suggest that it is nonfunctional. Agph-DAB2 is estimated to have diverged about 40 MYA from a previously characterized and highly polymorphic blackbird Mhc gene, Aph-DAB1, and is therefore younger than most mammalian Mhc paralogs and arose relatively late in avian evolution. Despite its nonfunctionality, Agph-DAB2 shows very high levels of nonsynonymous divergence from Agph-DAB1 and from reconstructed ancestral sequences in antigen-binding PBR codons-a strong indication of a period of adaptive divergence preceding loss of function. We also found that the region sequenced contains very few other unambiguous genes, a partial Mhcclass II gene fragment, and a paucity of simple-sequence and other repeats. Thus, this sequence exhibits some of the genomic streamlining expected for avian as compared with mammalian genomes, but is not as densely packed with functional genes as is the chicken Mhc.
PloS one, 2015
Secondary contact between divergent populations or incipient species may result in the exchange a... more Secondary contact between divergent populations or incipient species may result in the exchange and introgression of genomic material. We develop a simple DNA sequence measure, called Gmin, which is designed to identify genomic regions experiencing introgression in a secondary contact model. Gmin is defined as the ratio of the minimum between-population number of nucleotide differences in a genomic window to the average number of between-population differences. Although it is conceptually simple, one advantage of Gmin is that it is computationally inexpensive relative to model-based methods for detecting gene flow and it scales easily to the level of whole-genome analysis. We compare the sensitivity and specificity of Gmin to those of the widely used index of population differentiation, FST, and suggest a simple statistical test for identifying genomic outliers. Extensive computer simulations demonstrate that Gmin has both greater sensitivity and specificity for detecting recent int...
Genetics, 2003
The excess of rare variants in global sequencing studies of the nonrecombining portion of the Y c... more 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 p...
European journal of human genetics : EJHG, 2004
The relative roles of natural selection and accentuated genetic drift as explanations for the hig... more The relative roles of natural selection and accentuated genetic drift as explanations for the high frequency of more than 20 Ashkenazi Jewish disease alleles remain controversial. To test for the effects of a maternal bottleneck on the Ashkenazi Jewish population, we performed an extensive analysis of mitochondrial DNA (mtDNA) hypervariable segment 1 (HVS-1) sequence and restriction site polymorphisms in 565 Ashkenazi Jews from different parts of Europe. These patterns of variation were compared with those of five Near Eastern (n=327) and 10 host European (n=849) non-Jewish populations. Only four mtDNA haplogroups (Hgs) (defined on the basis of diagnostic coding region RFLPs and HVS-1 sequence variants) account for approximately 70% of Ashkenazi mtDNA variation. While several Ashkenazi Jewish mtDNA Hgs appear to derive from the Near East, there is also evidence for a low level of introgression from host European non-Jewish populations. HVS-1 sequence analysis revealed increased freq...
Bighorn sheep populations have greatly declined in numbers and distribution since European settle... more Bighorn sheep populations have greatly declined in numbers and distribution since European settlement, primarily because of high susceptibility to infectious diseases transmitted to them from domestic livestock. It has been suggested that low variation at major histocompatibility complex (MHC) genes, the most important genetic aspect of the vertebrate immune system, may result in high susceptibility to infectious disease. Therefore, we examined genetic polymorphism at a MHC gene (Ovca-DRB) in a large sample, both numerically and geographically, of bighorn sheep. Strikingly, there were 21 dierent alleles that showed extensive nucleotide and amino acid sequence divergence. In other words, low MHC variation does not appear to be the basis of the high disease susceptibility and decline in bighorn sheep. On the other hand, analysis of the pattern of the MHC polymorphism suggested that nonsynonymous substitutions predominated, especially at amino acids in the antigenbinding site. The average overall heterozygosity for the 16 amino acid positions that are part of the antigen binding site is 0.389 whereas that for the 67 amino acid positions not involved with antigen binding is 0.076. These ®ndings imply that the diversity present in this gene is functionally signi®cant and is, or has been, maintained by balancing selection. To examine the evolution of DRB alleles in related species, a phylogenetic analysis including other published ruminant (Bovidae and Cervidae) species, was carried out. An intermixture of sequences from bighorn sheep, domestic sheep, goats, cattle, bison, and musk ox was observed supporting trans-species polymorphism for these species. To reconcile the species and gene trees for the 104 sequences examined, 95`deep coalescent' events were necessary, illustrating the importance of balancing selection maintaining variation over speciation events.
Studies of human DNA sequence polymorphism reveal a range of diversity patterns throughout the ge... more 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.
Genome Biology and Evolution, 2014
One common form of reticulate evolution arises as a consequence of secondary contact between prev... more One common form of reticulate evolution arises as a consequence of secondary contact between previously allopatric populations. Using extensive coalescent simulations, we describe the conditions for, and extent of, the introgression of genetic material into the genome of a colonizing population from an endemic population. The simulated coalescent histories are sampled from models that describe the evolution of entire chromosomes, thereby allowing the expected length of introgressed haplotypes to be estimated. The results indicate that our ability to identify reticulate evolution from genetic data is highly variable and depends critically upon the duration of the period of allopatry, the timing of the secondary contact event, as well as the sizes of the populations at the time of contact. One particularly interesting result arises when secondary contact occurs close to the time of a severe founder event, in this case, genetic introgression can be substantially more difficult to detect. However, if secondary contact occurs after such a founding event, when the range of the colonizing population increases, introgression is more readily detectable across the genome. This result may have important implications for our ability to detect introgression between ancestrally bottlenecked modern human populations and archaic hominin species, such as Neanderthals.
Proceedings of the National Academy of Sciences, 2008
Nature Reviews Genetics, 2006
Molecular Biology and Evolution, 2000
To gain an understanding of the evolution and genomic context of avian major histocompatibility c... more To gain an understanding of the evolution and genomic context of avian major histocompatibility complex (Mhc) genes, we sequenced a 38.8-kb Mhc-bearing cosmid insert from a red-winged blackbird (Agelaius phoeniceus). The DNA sequence, the longest yet retrieved from a bird other than a chicken, provides a detailed view of the process of gene duplication, divergence, and degeneration (''birth and death'') in the avian Mhc, as well as a glimpse into major noncoding features of a songbird genome. The peptide-binding region (PBR) of the single Mhc class II B gene in this region, Agph-DAB2, is almost devoid of polymorphism, and a still-segregating single-base-pair deletion and other features suggest that it is nonfunctional. Agph-DAB2 is estimated to have diverged about 40 MYA from a previously characterized and highly polymorphic blackbird Mhc gene, Aph-DAB1, and is therefore younger than most mammalian Mhc paralogs and arose relatively late in avian evolution. Despite its nonfunctionality, Agph-DAB2 shows very high levels of nonsynonymous divergence from Agph-DAB1 and from reconstructed ancestral sequences in antigen-binding PBR codons-a strong indication of a period of adaptive divergence preceding loss of function. We also found that the region sequenced contains very few other unambiguous genes, a partial Mhcclass II gene fragment, and a paucity of simple-sequence and other repeats. Thus, this sequence exhibits some of the genomic streamlining expected for avian as compared with mammalian genomes, but is not as densely packed with functional genes as is the chicken Mhc.
Molecular Biology and Evolution, 1999
Twenty-three sequence haplotypes spanning the boundary of the second exon and intron of a red-win... more Twenty-three sequence haplotypes spanning the boundary of the second exon and intron of a red-winged blackbird Mhc class II B gene, Agph-DAB1, are presented. The polymorphism of the exon segment is distributed in two divergent allelic lineages which appear to be maintained by balancing selection. The silent nucleotide diversity of the exon ( Ï 0.101) is more than five times that of the intron ( Ï 0.018) and decays rapidly across the exonintron boundary. Additionally, genealogical reconstruction indicates that divergence from a common ancestor in the exon sample is over four times that of the intron. The intron sequences reveal a pattern of polymorphism which is characteristic of directional selection, rather than a pattern expected from linkage to a balanced polymorphism. These results suggest that the evolutionary histories of these two adjacent regions have been disassociated by recombination or gene conversion. The estimated population recombination parameter between the exon and the intron is sufficiently high (4N e c Ï 8.545) to explain the homogenization of intron sequences. Compatibility analyses estimate that these events primarily occur from the exon-intron boundary to about 20-30 bases into the intron. Additionally, the observation that divergent exon alleles share identical intron sequence supports the conclusion of disassociation of exon and intron evolutionary histories by recombination.
Molecular Biology and Evolution, 2013
Adaptive mutations that accumulate during species divergence are likely to contribute to reproduc... more Adaptive mutations that accumulate during species divergence are likely to contribute to reproductive incompatibilities and hinder gene flow; however, there may also be a class of mutations that are generally advantageous and can spread across species boundaries. In this study, we characterize a 15 kb region on chromosome 3R that has introgressed from the cosmopolitan generalist species Drosophila simulans into the island endemic D. sechellia, which is an ecological specialist. The introgressed haplotype is fixed in D. sechellia over almost the entirety of the resequenced region, whereas a core region of the introgressed haplotype occurs at high frequency in D. simulans. The observed patterns of nucleotide variation and linkage disequilibrium are consistent with a recently completed selective sweep in D. sechellia and an incomplete sweep in D. simulans. Independent estimates of both the time to the introgression and sweep events are all close to 10,000 years before the present. Interestingly, the most likely target of selection is a highly occupied transcription factor binding region. This work confirms that it is possible for mutations to be globally advantageous, despite their occurrence in divergent genomic and ecological backgrounds.
Molecular Biology and Evolution, 2010
A large number of statistical tests have been proposed to detect natural selection based on a sam... more A large number of statistical tests have been proposed to detect natural selection based on a sample of variation at a single genetic locus. These tests measure the deviation of the allelic frequency distribution observed within populations from the distribution expected under a set of assumptions that includes both neutral evolution and equilibrium population demography. The present study considers a new way to assess the statistical properties of these tests of selection, by their behavior in response to direct perturbations of the steady-state allelic frequency distribution, unconstrained by any particular nonequilibrium demographic scenario. Results from Monte Carlo computer simulations indicate that most tests of selection are more sensitive to perturbations of the allele frequency distribution that increase the variance in allele frequencies than to perturbations that decrease the variance. Simulations also demonstrate that it requires, on average, 4N generations (N is the diploid effective population size) for tests of selection to relax to their theoretical, steady-state distributions following different perturbations of the allele frequency distribution to its extremes. This relatively long relaxation time highlights the fact that these tests are not robust to violations of the other assumptions of the null model besides neutrality. Lastly, genetic variation arising under an example of a regularly cycling demographic scenario is simulated. Tests of selection performed on this last set of simulated data confirm the confounding nature of these tests for the inference of natural selection, under a demographic scenario that likely holds for many species. The utility of using empirical, genomic distributions of test statistics, instead of the theoretical steady-state distribution, is discussed as an alternative for improving the statistical inference of natural selection.
Molecular Biology and Evolution, 2004
The human RRM2P4 pseudogene has a pattern of nucleotide polymorphism that is unlike any locus pub... more 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.
Immunological Reviews, 1999
We review recent developments in the ongoing study of the evolution ofthe Mhc gene family iti bir... more We review recent developments in the ongoing study of the evolution ofthe Mhc gene family iti birds, with emphasis on class IIB genes and results from songbirds obtained in our laboratory. Southern blots suggest a surprising diversity in Mhc class II gene ntunber among various songbird species (Passeriformes), We have sequenced ~30 kb contigs from Mhc-bearing cosmid clones from two species, red-winged blackbirds (Agelaius phoeniceus) and bouse fincbes (Carpodacus mexicanus). whose demography, lifetime reproductive success, epizootics, parasitology and mate choice are among the best studied for natural populations of birds. Of three genes cloned from these species, only one appears strongly polymorphic, and one (from ihe bouse finch) is likely a pseudogene. All are similar in structure to those in chickens, albeit witb introns intermediate in length between chickens and mammals. Phylogenetic analysis of avaiiabie class II B peptide-binding region exons suggests that the overwhelming longterm force operating on avian genes sampled thus far has been post-speciation gene duplication and/or concerted evolution. These and other results suggest that the evolution of class II B genes in birds conforms to a mixture of several models of multigene family evolution proposed for the mamtnalian Mk. incorporating ongoing homogenization. duplication and pseudogene formation. Large-scale sequencing studies in these and other species, though still in their infancy, will prove invaluable for studying the comparative structures of avian Mhcs, as well as patterns of selection, mutation and linkage disequilibrium at several scales.
Immunogenetics, 2001
Twelve MHC class I exon 2 sequences were uncovered in a sample from the endangered Sacramento Riv... more Twelve MHC class I exon 2 sequences were uncovered in a sample from the endangered Sacramento River winter-run Chinook salmon in the central valley of California. Phylogenetic analysis of the 12 sequences indicates that the alleles descend from two of six major allelic lineages found among four Pacific salmon species. Nine of the 12 alleles belong to an allelic lineage that began diversifying 8 million years ago, just prior to the estimated time of Chinook speciation. The most recent common ancestor of all 12 winter-run alleles is estimated to be 15 million years ago, approximately 5 million years before the radiation of the Pacific salmon species. The average nonsynonymous distance among the peptide binding-region codons of exon 2 for the 12 alleles is significantly higher than the average synonymous distance in these codons. We estimate the symmetrical overdominant selection coefficient against homozygotes for this exon to be 0.038. Thus, strong positive and balancing selection has maintained functional diversity in the peptide-binding region of the exon over millions of years and this variation has not yet been substantially eliminated by increased genetic drift due to the recent dramatic decline in abundance of this Chinook salmon population.
Human Genetics, 2004
In order to investigate the genetic history of autochthonous South Siberian populations and to es... more In order to investigate the genetic history of autochthonous South Siberian populations and to estimate the contribution of distinct patrilineages to their gene pools, we have analyzed 17 Y-chromosomal binary markers (YAP, RPS4Y 711 , SRY-8299, M89, M201, M52, M170, 12f2, M9, M20, 92R7, SRY-1532, DYS199, M173, M17, Tat, and LLY22 g) in a total sample of 1,358 males from 14 ethnic groups of Siberia (Altaians-
Uploads
Papers by Daniel Garrigan