Invasion genetics of the Mediterranean fruit fly: variation in multiple nuclear introns

Biological invasions generally start from low initial population sizes, leading to reduced genetic variation in nuclear and especially mitochondrial DNA. Consequently, genetic approaches for the study of invasion history and population structure are difficult. An extreme example is the Mediterranean fruit fly, Ceratitis capitata (Medfly), for which successive invasions during this century have resulted in a loss of 60% of ancestral genetic variation in isozymes and 75% of variation in mitochondrial DNA. Using Medflies as an example, we present a new approach to invasion genetics that measures DNA sequence variation within introns from multiple nuclear loci. These loci are so variable that even relatively recently founded Medfly populations within California and Hawaii retain ample genetic diversity. Invading populations have only lost 35% of the ancestral genetic variation. Intron variation will allow high-resolution genetic characterization of invading populations in both natural and managed systems, although non-equilibrium methods of analysis may be necessary if the genetic diversity represents sorting ancestral polymorphism.     

Comparative molecular phylogeography of two Xenopus species, X. gilli and X. laevis, in the south-western Cape Province, South Africa

Xenopus gilli is a vulnerable anuran with a patchy distribution along the south-western coast of the Cape Province, South Africa. This species is sympatric with Xenopus laevis laevis, a widespread relative found over much of southern Africa. We examined the molecular phylogeography and population structure of the contact zone between these species to obtain information about historical biogeography and conservation management of this region. Analyses of the distribution, frequency, and cladistic and phenetic relationships among mitochondrial DNA haplotypes indicate that population subdivision is present in both taxa but that long-term isolation of sets of populations has occurred in X. gilli only. Haplotype and nucleotide diversity are also considerably higher within and among X. gilli ponds than X. l. laevis ponds in this region. We attribute the genetic segregation of X. gilli populations to ancient habitat fragmentation by ocean transgression into X. gilli habitat and to continued habitat alteration by human activity. The lower level of genetic diversity in X. L. laevis in this region is likely a result of a recent arrival of this taxon to the south-western Cape region relative to X. gilli. Population structure in X. l. laevis may be a result of isolation by distance. Clear evidence exists for at least two management units within X. gilli and strongly supports the establishment of protective measures east of False Bay in order to conserve a substantial portion of this species' extant genetic diversity.     

Phylogeography and pleistocene evolution in the North American black bear

To determine the extent of phylogeographic structuring in North American black bear (Ursus americanus) populations, we examined mitochondrial DNA sequences (n = 118) and restriction fragment length polymorphism profiles (n = 258) in individuals from 16 localities. Among the bears examined, 19 lineages falling into two highly divergent clades were identified. The clades differ at 5.0% of nucleotide positions, a distance consistent with an origin 1.8 MYA, and have different but overlapping geographical distributions. Areas of clade cooccurrence show that eastern and western populations are currently mixing, but regional differences in lineage distribution suggest that mixing has begun only recently. The long-term population history of black bears appears to be characterized predominantly by long-term regional isolation followed by recent contact and hybridization. Congruence between the pattern of diversity observed in black bears and patterns of forest refuge formation during the Pleistocene supports earlier speculation that Pleistocene forest fragmentations underlie a common pattern in the phylogeography of North American forest taxa.     

Genic homozygosity in an ancient reptile (Alligator mississippiensis)

Proteins assayed electrophoretically showed variation at only three of 49 presumed genetic loci in alligators from southwestern Louisiana. Average heterozygosity per individual was 0.021+/-0.012; proportion of polymorphic loci was 0.06. Data on the history, structure, and ecology of this alligator population are consistent with natural selection as the primary factor accounting for this low genetic variability. However, neither a historic population bottleneck nor some genetic mechanism limiting variability can be dismissed as a possible factor.     

A selective difference between human Y-chromosomal DNA haplotypes

DNA analysis is making a valuable contribution to the understanding of human evolution [1]. Much attention has focused on mitochondrial DNA (mtDNA) [2] and the Y chromosome [3] [4], both of which escape recombination and so provide information on maternal and paternal lineages, respectively. It is often assumed that the polymorphisms observed at loci on mtDNA and the Y chromosome are selectively neutral and, therefore, that existing patterns of molecular variation can be used to deduce the histories of populations in terms of drift, population movements, and cultural practices. The coalescence of the molecular phylogenies of mtDNA and the Y chromosome to recent common ancestors in Africa [5] [6], for example, has been taken to reflect a recent origin of modern human populations in Africa. An alternative explanation, though, could be the recent selective spread of mtDNA and Y chromosome haplotypes from Africa in a population with a more complex history [7]. It is therefore important to establish whether there are selective differences between classes (haplotypes) of mtDNA and Y chromosomes and, if so, whether these differences could have been sufficient to influence the distributions of haplotypes in existing populations. A precedent for this hypothesis has been established for mtDNA in that one mtDNA background increases susceptibility to Leber hereditary optic neuropathy [8]. Although studies of nucleotide diversity in global samples of Y chromosomes have suggested an absence of recent selective sweeps or bottlenecks [9], selection may, in principle, be very important for the Y chromosome because it carries several loci affecting male fertility [10] [11] and as many as 5% of males are infertile [11] [12]. Here, we show that one class of infertile males, PRKX/PRKY translocation XX males, arises predominantly on a particular Y haplotypic background. Selection is, therefore, acting on Y haplotype distributions in the population.     

Contrasting population structure from nuclear intron sequences and mtDNA of humpback whales

Powerful analyses of population structure require information from multiple genetic loci. To help develop a molecular toolbox for obtaining this information, we have designed universal oligonucleotide primers that span conserved intron-exon junctions in a wide variety of animal phyla. We test the utility of exon-primed, intron-crossing amplifications by analyzing the variability of actin intron sequences from humpback, blue, and bowhead whales and comparing the results with mitochondrial DNA (mtDNA) haplotype data. Humpback actin introns fall into two major clades that exist in different frequencies in different oceanic populations. It is surprising that Hawaii and California populations, which are very distinct in mtDNAs, are similar in actin intron alleles. This discrepancy between mtDNA and nuclear DNA results may be due either to differences in genetic drift in mitochondrial and nuclear genes or to preferential movement of males, which do not transmit mtDNA to offspring, between separate breeding grounds. Opposing mtDNA and nuclear DNA results can help clarify otherwise hidden patterns of structure in natural populations.     

Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data

When a population experiences a reduction of its effective size, it generally develops a heterozygosity excess at selectively neutral loci, i.e., the heterozygosity computed from a sample of genes is larger than the heterozygosity expected from the number of alleles found in the sample if the population were at mutation drift equilibrium. The heterozygosity excess persists only a certain number of generations until a new equilibrium is established. Two statistical tests for detecting a heterozygosity excess are described. They require measurements of the number of alleles and heterozygosity at each of several loci from a population sample. The first test determines if the proportion of loci with heterozygosity excess is significantly larger than expected at equilibrium. The second test establishes if the average of standardized differences between observed and expected heterozygosities is significantly different from zero. Type I and II errors have been evaluated by computer simulations, varying sample size, number of loci, bottleneck size, time elapsed since the beginning of the bottleneck and level of variability of loci. These analyses show that the most useful markers for bottleneck detection are those evolving under the infinite allele model (IAM) and they provide guidelines for selecting sample sizes of individuals and loci. The usefulness of these tests for conservation biology is discussed.     

Mitochondrial lineages in Ladin-speaking communities of the eastern Alps

European mitochondrial alleles cluster into five haplogroups. Haplogroup 2 is rare in general, but represents more than half of the few known sequences among Ladin speakers of the Alps. Here we describe DNA diversity in control region I of the hypervariable D-loop in 43 Ladins, and in 25 Italian speakers. Analysis of these data, and of previously published sequences, confirms a high degree of differentiation among Ladins and their geographical neighbours. This cannot be regarded as a simple effect of isolating factors, geographic or linguistic, as diversity is high within Ladin communities too. Rather, allele genealogies, population trees, and principal component analysis suggest a relationship between Ladin and Near Eastern samples. Two evolutionary hypotheses seem compatible with these findings. The view whereby Ladins could be descended from Palaeolithic inhabitants of the Alps is supported by the identification, in this study, of the probable ancestral haplotype of group 2, never previously observed in central Europe. Alternatively, a comparatively recent, Neolithic immigration of the ancestors of current Ladin speakers seems consistent with recent linguistic theories. In both cases, the number of lineages present, and their extensive diversity, are not compatible with a serious bottleneck in the Ladin population's history.     

Different patterns of variation at the X- and Y-chromosome-linked microsatellite loci DXYS156X and DXYS156Y in human populations

We compared the global pattern of variation at two homologous microsatellites mapping to the long arm of the X chromosome (DXYS156X) and to the short arm of the Y chromosome (DXYS156Y) in humans. A single pair of oligonucleotide primers amplifies these two nonallelic loci, each of which contains polymorphism in the number of pentanucleotide units. We observed 11 alleles in a sample of 2290 X chromosomes and 2006 Y chromosomes from 50 populations representing 6 major geographic regions. The overlapping size range of the X- and Y-chromosome alleles indicated a more complex distribution of alleles at these two loci than previously reported. Contrasting patterns of X-chromosome-linked and Y-chromosome-linked variation were reflected in statistically significant differences in genetic diversity values among geographic regions and between X and Y chromosomes. Higher levels of diversity characterized the DXYS156X locus in Africa (0.799 +/- 0.004) and the DXYS156Y locus in East Asia (0.700 +/- 0.006) compared with populations from other regions. These different patterns of variation can be explained by a combination of processes at both the molecular and population levels, including variable mutation rates, different effective population sizes, and genetic drift.     

Extended gene diversity at the FMR1 locus and neighbouring CA repeats in a sub-Saharan population

We report on the allele distributions in a normal black African population at two microsatellite loci neighbouring the FRAXA locus and at the CGG repeat in the 5' end of the FMR-1 gene, which causes the fragile X syndrome. The CGG repeat distribution was found to be similar to that of other ethnic groups, as well as to that of other nonhuman primates, possibly predicting a comparable prevalence of fragile X in Africa. Significant linkage disequilibrium has been observed between fragile X mutations and alleles of the DXS548 and FRAXAC1 loci in European and Asian populations, and some founder chromosomes may be extremely old. Those associated with FRAXAC1-A and DXS548-2 alleles are not present in the Asian fragile X samples. We searched for these alleles and their frequency in the well defined Bamileke population of Cameroon. All previously described alleles and some new ones were found in this sample, supporting the hypothesis of their pre-existence and subsequent loss in Asian populations. Finally, the heterozygosity of the Bamileke sample was significantly higher at both marker loci and comparable to that of Europeans at the CGG repeat, confirming the notion that genetic diversity is greater in Africans than in other groups and supporting the view that evolution of modern man started in Africa.     

Mitochondrial DNA diversity in the Kuna Amerinds of Panamá

Mitochondrial DNA (mtDNA) haplotype diversity was determined for 63 Chibcha-speaking Kuna Amerinds sampled widely across their geographic range in eastern Panamá. The Kuna data were compared with mtDNA control region I sequences from two neighboring Chibchan groups, the Ng?bé and the Huetar; two Amerind groups located at the northern and southern extremes of Amerind distribution, the Nuu-Chah-Nulth of the Pacific Northwest and the Chilean Mapuche; and with a single Na-Dene group, the Haida of the Pacific Northwest. The Kuna exhibited low levels of mitochondrial diversity as had been reported for the other two Chibchan groups and, furthermore, carried only two of the four Amerind founding lineages first reported by Schurr and coworkers (Am. J. Hum. Genet. 1990; 46: 613-623). We posit that speakers of modern Chibchan languages (henceforth referred to as the Chibcha) passed through a population bottleneck caused either by ethnogenesis from a small founding population and/or subsequent European conquest and colonization. Using the approach of Harpending et al. (Curr. Anthropol. 1993; 34: 483-496), we estimated a Chibchan population bottleneck and subsequent expansion approximately 10,000 years before present, a date consistent with a bottleneck at the time of Chibchan ethnogenesis. The low mtDNA diversity of Kuna Amerinds, as opposed to the generally high levels of mtDNA variation detected in other Amerind groups, demonstrates the need for adequate sampling of cultural or racial groups when attempting to genetically characterize human populations.     

Habitat islands, genetic diversity, and gene flow in a Patagonian rodent

The effects of terrestrial habitat islands on gene flow and genetic diversity in animal populations have been predicted and discussed in theoretical terms, but empirical data are needed to test these predictions and provide an understanding of the relationships of life-history characteristics to genetics of insular species. We studied saxicolous mice (Phyllotis xanthopygus) in Patagonia to explore genetic structure, phylogeography, and gene flow in a species inhabiting natural habitat islands. Phylogeographic analyses based on mtDNA sequences revealed two haplotype clades, which presumably reflect early Pleistocene factors that temporarily separated the mice into two geographically isolated groups. The Río Chubut, which lies within a glacial drainage basin bisecting northern Patagonia, might have affected gene flow in the species. Although we anticipated isolation by distance and founder phenomena associated with habitat islands, in some habitat patches we found evidence of high local genetic diversity. The amount of divergence in the mitochondrial cytochrome b gene (approximately 3.4%) in animals at a single locality could best be explained through a combination of historical factors and metapopulation source-sink theory. Demographic shifts, dispersal, and episodic recolonization are important in the life history and genetic population structure of P. xanthopygus.     

Population substructure and isolation by distance in three continental regions

Isolation by distance and divergence from a shared population history are two sources of population substructure. Isolation by distance erases population history as populations approach migration-drift equilibrium, while diverging populations descended from a single ancestral population will accumulate genetic differences with time. Here I investigate how much of the worldwide genetic diversity from Jorde et al.'s ([1997] Proc. Natl. Acad. Sci. USA 94:3100-3103) 60 tetranucleotide short tandem repeat (STR) data can be explained by isolation by distance. I use Slatkin's measure of population substructure, R(ST), principal components analyses, and Mantel tests to investigate the pattern of genetic diversity at both the intercontinental and intracontinental levels. Geographic distance accounts for almost 60% of world-wide interpopulation genetic relationships. Within continents, the correlations are less, although not significantly so because of wide confidence intervals. These results suggest that populations have not reached migration-drift equilibrium and that there is information in STR data to reconstruct population history. The level of population substructure worldwide is consistent with previous observations, but at the intracontinental level substructure is less. When one examines diversity against distance from the centroid, one sees excess heterozygosity in Africa, a pattern also noted by Stoneking et al. ([1998] Genome Research 7:1061-1071). A larger effective population size in Africa could explain the excess diversity. Greater gene flow in Africa is an unlikely explanation because the African R(ST) value is slightly larger than the Asian and European R(ST)s, pointing to less gene flow and greater substructure among African populations. Furthermore, there are differences in patterns between heterozygosity and allele size variance. Heterozygosity has a higher correlation with distance from the centroid than does allele size variance, and this may reflect demographic history. Kimmel et al. ([1998] Genetics 148:1921-1930) have shown that after a population expansion heterozygosity returns to equilibrium more quickly than does allele size variance. The contrasting patterns between heterozygosity and allele size variance may reflect different times after an expansion. However, simulations and further work need to be done to more thoroughly investigate the possibility that these data reflect population expansions.     

Field trial to evaluate the immunogenicity of pseudorabies virus vaccines with deletions for glycoproteins G and E

This study is a geographically systematic genetic survey of the easternmost subspecies of chimpanzee, Pan troglodytes schweinfurthii. DNA was noninvasively collected in the form of shed hair from chimpanzees of known origin in Uganda, Rwanda, Tanzania, and Za?re. Two hundred sixty-two DNA sequences from hypervariable region 1 of which of the mitochondrial control region were generated. Eastern chimpanzees display levels of mitochondrial genetic variation which are low and which are similar to levels observed in humans (Homo sapiens). Also like humans, between 80% and 90% of the genetic variability within the eastern chimpanzees is apportioned within populations. Spatial autocorrelation analysis shows that genetic similarity between eastern chimpanzees decreases clinically with distance, in a pattern remarkably similar to one seen for humans separated by equivalent geographic distances. Eastern chimpanzee mismatch distributions (frequency distributions of pairwise genetic differences between individuals) are similar in shape to those for humans, implying similar population histories of recent demographic expansion. The overall pattern of genetic variability in eastern chimpanzees is consistent with the hypothesis that the subject has responded demographically to paleoclimatically driven changes in the distribution of eastern African forests during the recent Pleistocene.     

Mitochondrial DNA haplotype frequencies in natural and experimental populations of Drosophila subobscura

The evolution of Drosophila subobscura mitochondrial DNA has been studied in experimental populations, founded with flies from a natural population from Esporles (Majorca, Balearic Islands, Spain). This population, like other European ones, is characterized by the presence of two very common (>96%) mitochondrial haplotypes (called I and II) and rare and endemic haplotypes that appear at very low frequencies. There is no statistical evidence of positive Darwinian selection acting on the mitochondrial DNA variants according to Tajima's neutrality test. Two experimental populations, with one replicate each, were established with flies having a heterogeneous nuclear genetic background, which was representative of the composition of the natural population. Both populations were started with the two most frequent mitochondrial haplotypes, but at different initial frequencies. After 13 to 16 generations, haplotype II reached fixation in three cages and its frequency was 0.89 by generation 25 in the fourth cage. Random drift can be rejected as the force responsible for the observed changes in haplotype frequencies. There is not only statistical evidence of a linear trend favoring a mtDNA (haploid) fitness effect, but also of a significant nonlinear deviation that could be due to a nuclear component.     

Finding the genes in genomic DNA

We analyzed two-locus disequilibria for 16 polymorphic loci of seven susceptibility genes for coronary artery disease located in five chromosomal regions distributed across four chromosomes. Included were the genes coding for apolipoprotein B (ApoB, chromosome 2, four marker loci), lipoprotein lipase (LPL, chromosome 8, three marker loci), apolipoproteins AI, CIII AIV (ApoAI-CIII-AIV, chromosome 11, three marker loci), apolipoprotein E (ApoE, chromosome 19, two marker loci), and the low density lipoprotein receptor (LDLR, chromosome 19, four marker loci). Our sample included 540 unrelated individuals from the Rochester, Minn. population. There were no statistically significant deviations of single-locus genotypes from Hardy-Weinberg equilibrium. The strongest associations within genes were for composite diallelic disequilibria; 17/19 were significant (13 at Pr <0.001, 1 at Pr <0.01, 3 at Pr <0.05). These observations suggest marker alleles within genes have a shared evolutionary history reflected by disequilibria that have not been dissipated by recombination. Disequilibrium was not generally concordant with the physical orderings of markers. Only two significant higher-order disequilibria were observed although 12 triallelic disequilibria were at maximum possible values. We observed 19 statistically significant disequilibria (Pr <0.05; 4 composite diallelic, 13 triallelic, and 2 quadriallelic) between 101 pairs of marker loci, where each locus in a pair was from a different unlinked region. These unexpected results are most likely explained by recent historical factors, including worldwide population expansion and amalgamation with continuous admixture, that influence the genetic structure (organization of alleles and non-alleles into genotypes) of a population. We conclude that disequilibria between loci from unlinked regions may be more extensive than is commonly assumed. Our findings also suggest that it is, on average, at least 15 times more likely to not detect significant disequilibrium among unlinked loci when it is really present than to make a false positive inference. Disequilibria between functional loci within or between regions will impact estimates of genetic variance associated with particular functional mutations within a susceptibility gene region.     

The genetic relationship between the Finns and the Finnish Saami (Lapps): analysis of nuclear DNA and mtDNA

The genetic relationships between two Finno-Ugric-speaking populations, the Finns and the Finnish Saami (Lapps), were studied by using PCR for six nuclear-DNA marker loci, mitochondrial restriction-site polymorphism, and sequence variation of a 360-bp segment of the mitochondrial control region. The allele frequencies of each of the nuclear-DNA marker loci and the frequencies of mtDNA restriction haplotypes were significantly different between the populations. The Saami showed exceptionally low variation in their mtDNA restriction sites. The 9-bp deletion common in East Asian populations was not observed, nor did the haplotype data fit into the haplogroup categorization of Torroni et al. The average number of nucleotide substitutions from the mtDNA haplotype data indicated that the Finnish Saami may be closer to the Finns than to the other reference populations, whereas nuclear DNA suggested that the Finns are more closely related to the European reference populations than to the Finnish Saami. The similarity of the Finns to the other Europeans was even more pronounced according to the sequence data. We were unable to distinguish between the Finns and either the Swiss or Sardinian reference populations, whereas the Finnish Saami clearly stood apart. The Finnish Saami are distinct from other Circumarctic populations, although two of the lineages found among the Saami showed closer relationship to the Circumarctic than to the European lineages. The sequence data indicated an exceptionally high divergence for the Saami mtDNA control lineages. The distribution of the pairwise nucleotide differences in the Saami suggested that this population has not experienced an expansion similar to what was indicated for the Finns and the reference populations.     

Genetic evidence for the proto-Austronesian homeland in Asia: mtDNA and nuclear DNA variation in Taiwanese aboriginal tribes

Previous studies of mtDNA variation in indigenous Taiwanese populations have suggested that they held an ancestral position in the spread of mtDNAs throughout Southeast Asia and Oceania (Melton et al. 1995; Sykes et al. 1995), but the question of an absolute proto-Austronesian homeland remains. To search for Asian roots for indigenous Taiwanese populations, 28 mtDNAs representative of variation in four tribal groups (Ami, Atayal, Bunun, and Paiwan) were sequenced and were compared with each other and with mtDNAs from 25 other populations from Asia and Oceania. In addition, eight polymorphic Alu insertion loci were analyzed, to determine if the pattern of mtDNA variation is concordant with nuclear DNA variation. Tribal groups shared considerable mtDNA sequence identity (P>.90), where gene flow is believed to have been low, arguing for a common source or sources for the tribes. mtDNAs with a 9-bp deletion have considerable mainland-Asian diversity and have spread to Southeast Asia and Oceania through a Taiwanese bottleneck. Only four Taiwanese mtDNA haplotypes without the 9-bp deletion were shared with any other populations, but these shared types were widely dispersed geographically throughout mainland Asia. Phylogenetic and principal-component analyses of Alu loci were concordant with conclusions from the mtDNA analyses; overall, the results suggest that the Taiwanese have temporally deep roots, probably in central or south China, and have been isolated from other Asian populations in recent history.     

Patterns of extensive genetic differentiation and variation among European harbor seals (Phoca vitulina vitulina) revealed using microsatellite DNA polymorphisms

The harbor seal (Phoca vitulina) has the most extensive distribution of any phocid seal species. An analysis of population structure in this species across its European range was made using 7 phocid derived microsatellites in a sample of 1,029 individuals from 12 separate geographic areas. Despite the species potential for long-distance movement, significant genetic differentiation between areas was observed using an unbiased estimator of RST. Six distinct population units were identified: Ireland-Scotland, English east coast, Waddensea, western Scandinavia (Norway-Kattegat-Skagerrak-west Baltic), east Baltic, and Iceland. Little local substructuring is present along coastlines with a continuous distribution of breeding animals, but differentiation does increase with geographic distance. The degree of differentiation is greater over equivalent distances where the distribution is discontinuous, such as along coasts where breeding colonies are separated by large distances or by stretches of open sea. Patterns of population differentiation derived from microsatellites are very similar to those obtained from previous mitochondrial DNA analysis and suggest that philopatry in harbor seals operates over 300-500 km. In Europe, harbor seals have experienced a complex demographic history and patterns of population structure are likely to have been affected by natural environmental influences such as Pleistocene glaciations and epizootics. Comparison of Nm values from an unbiased estimator of RST, GST, and theta are consistent and, in some cases, may indicate populations where conditions deviate from the expectations of the RST model.     

Mutation and evolution of microsatellites in Drosophila melanogaster

Levels of nucleotide polymorphism in the Drosophila melanogaster genome are correlated with rates of recombination. This relationship may be due to hitchhiking of advantageous mutations (selective sweeps) or to continual removal of deleterious mutations from the genome (background selection). One test of the relative contributions of selective sweeps and background selection to the observed levels of variation in the genome of D. melanogaster is to compare levels of nucleotide variability (with a mutation rate on the order of 10(-9) per nucleotide per generation) with more rapidly evolving DNA loci such as microsatellites. This test depends critically on details of the mutational process of microsatellites. In this paper, we summarize our studies of microsatellite characteristics and mutation rates in D. melanogaster. We find that D. melanogaster microsatellites are short and have a mutation rate (6.5 x 10(-6) per locus per generation) several orders of magnitude lower than mammals studied to date. We further show that genetic variation at 18 dinucleotide repeat microsatellites in a population of D. melanogaster from Maryland is correlated with regional rates of recombination. These and other microsatellite data suggest that both background selection and selective sweeps may contribute to the correlation between DNA sequence variation and recombination in Drosophila.