The endangered Sonoran topminnow: Examination of species and ESUs using three mtDNA genes

There has been controversy over the species status of Sonoran topminnows and debate about the presence of ESUs in the Gila topminnow. From examination of sequence variation at 2626 base pairs over three mtDNA genes, we found a 29 (1.1%) nucleotide genetic difference between Gila and Yaqui topminnows. This provides strong support that these two taxa are separate species, Poeciliopsis occidentalis (Gila topminnow) and P. sonoriensis (Yaqui topminnow) and have been separated for approximately one million years. All the Gila topminnows within Arizona have the same sequence for the three mtDNA genes, that is, there is not reciprocal monophyly for mtDNA sequence data for the two previously designated ESUs. However, evidence of the unique habitat for Monkey Spring, its long-term isolation from other Gila topminnow habitats, and the presence of unique fish and invertebrate taxa in Monkey Spring support the designation of the Monkey Spring topminnows as an ESU. Finally, theoretical considerations using molecular data and estimates of heterozygosity and genetic distance for nuclear genes between populations of the Gila topminnow show that the lack of mtDNA variation is not inconsistent with the level and pattern of nuclear genetic variation observed.     

MHC VARIATION IN THE ENDANGERED GILA TOPMINNOW

Sequence variation at a major histocompatibility complex (MHC) gene, assumed to be involved in parasite and pathogen resistance, was examined in the endangered Gila topminnow (Poeciliopis o. occidentalis), from the four watersheds where they remain in the United States. This is the first estimate of variation in MHC genes in natural populations of an endangered species. The population that has experienced the most bottlenecks historically was monomorphic for MHC variation. Another population, which earlier had been found to be the only population polymorphic for allozymes, had five MHC alleles, four different from those found in the other populations. Overall, nine different alleles were found. The four populations were highly divergent at MHC with four of the six population pairs not sharing any alleles. However, the magnitude of differentiation between populations on the amino-acid level varied fivefold for the populations that shared no alleles. Using single-stranded conformational polymorphism (SSCP), these alleles segregated consistently with Mendelian expectations in families. Because of the high genetic differentiation between these populations for a potentially adaptive gene, we recommend that the four watersheds be examined further for separate conservation and management.     

MHC variation and tissue transplantation in fish

Major histocompatibility complex (MHC) genes were originally discovered because of their role in tissue rejection in mammals and have subsequently been implicated in the incidence of autoimmune diseases and resistance to infectious diseases. Here we present the first demonstration that a gene defined by molecular sequence in the fish MHC, specifically a class II locus, plays an important role in tissue rejection. This effect in the endangered Gila topminnows appears to be additive and depends on the number of MHC alleles shared between the host and the recipient fish of the scale transplants. In addition, there was lower success of scale transplants in MHC-matched individuals in a population with high microsatellite variation than in a population with low variation. This suggests that other loci, presumably other MHC loci, play a significant role in transplantation success in fishes, as they do in mammals.     

Genetic variation of the major histocompatibility complex (MHC class II β gene) in the threatened Gila trout, <Emphasis Type="Italic">Oncorhynchus gilae gilae</Emphasis>

Gila trout (Oncorhynchus gilae gilae) was federally protected in 1973 because of severe declines in abundance and geographic range size. At present, four relict genetic lineages of the species remain in mountain streams of New Mexico and Arizona, USA. Management actions aimed at species recovery, including hatchery production and restocking of formerly occupied streams, have been guided by information from non-functional genetic markers. In this study, we investigated genetic variation at exon 2 of the major histocompatibility complex (MHC) class II β gene that is involved in pathogen resistance and thus presumably under natural selection. Phylogenetic analysis revealed trans-species polymorphism and a significantly high ratio of non-synonymous to synonymous amino acid changes consistent with the action of historical balancing selection that maintained diversity at this locus in the past. However, Gila trout exhibited low allelic diversity (five alleles from 142 individuals assayed) compared to some other salmonid fishes, and populations that originated exclusively from hatcheries possessed three or fewer MHC alleles. Comparative analysis of genetic variation at MHC and six presumably neutrally evolving microsatellite loci revealed that genetic drift cannot be rejected as a primary force governing evolution of MHC in contemporary populations of Gila trout. Maintenance of diversity at MHC will require careful implementation of hatchery breeding protocols and continued protection of wild populations to prevent loss of allelic diversity due to drift.     

Discriminating coho salmon (Oncorhynchus kisutch) populations within the Fraser River, British Columbia, using microsatellite DNA markers

Three microsatellite loci were used to examine genetic variation among 16 coho salmon ( Oncorhynchus kisutch ) populations within the Fraser River drainage system, in British Columbia, Canada. Each locus was highly polymorphic with 30 alleles at the Ots 101 locus, 15 alleles at the Ots 3 locus and 38 alleles at the Ots 103 locus. Average observed heterozygosities were 86.1%, 70%, and 56.1%, respectively. With the exception of the Dunn and Lemieux River populations, Chi-square tests and F _ST values indicated that all populations had significantly different allele frequencies. Two distinct population groups within the Fraser River drainage were observed. Lower Fraser River populations were strongly differentiated from populations spawning in the upper Fraser River, which includes the Thompson River (a tributary flowing into the upper Fraser) and the portion of the Fraser River beyond the precipitous Fraser River canyon. This regional population structure may have resulted from colonization of the upper and lower Fraser River regions by different founder populations following Pleistocene glaciation, and be maintained by adaptive differences between the two groups of coho salmon. Coho salmon populations in the upper Fraser and Thompson River drainages form an evolutionarily significant unit (ESU) of importance for conservation of biodiversity in coho salmon. Microsatellite DNA loci show promise as technically simple and highly informative genetic markers for coho salmon population management.     

Low diversity in the major histocompatibility complex class II DRB1 gene of the Spanish ibex, Capra pyrenaica

During the last two centuries, the Spanish ibex (Capra pyrenaica) has shown a significant demographic decline as a result of the progressive destruction of its natural habitat, disease epidemics, and uncontrolled hunting. Partial sequencing of the class II MHC DRB1 gene revealed that the Spanish ibex has remarkably low levels of genetic variation at this locus, with only six different DRB1 alleles and an observed heterozygosity of 0.429-0.579. The rates of nonsynonymous vs synonymous substitutions were significantly different in the peptide-binding region (dN/dS=5.347, P=0.002), a feature that indicates that the DRB1 gene is under positive selection. A phylogenetic analysis of the Spanish ibex and a set of domestic goat DRB1 alleles revealed that the reported sequences represent four major allelic lineages. The limited allelic repertoire of the DRB1 gene in the Spanish ibex is likely the direct result of the recent history of population bottlenecks and marked demographic decline of this species. A genetic survey of 13 microsatellite loci was consistent with this idea. The Spanish ibex subspecies C. p. hispanica and C. p. victoriae consistently showed considerably lower levels of microsatellite heterozygosity (Ho=0.184-0.231) and allelic diversity (mean number of alleles per locus=2-2.4) than those reported in other wild ruminants. This study demonstrates the significance of both natural selection and the demographic history of populations in determining patterns of genetic variation at MHC loci. In addition, our results emphasize the importance of locally adapted populations for the preservation of genetic diversity.     

Phylogeographical population structure of tiger quolls Dasyurus maculatus (Dasyuridae: Marsupialia), an endangered carnivorous marsupial

Tiger quolls, Dasyurus maculatus, are the largest carnivorous marsupials still extant on the mainland of Australia, and occupy an important ecological niche as top predators and scavengers. Two allopatric subspecies are recognized, D.m. gracilis in north Queensland, and D.m. maculatus in the southeast of the mainland and Tasmania. D.m. gracilis is considered endangered while D.m. maculatus is listed as vulnerable to extinction; both subspecies are still in decline. Phylogeographical subdivision was examined to determine evolutionarily significant units (ESUs) and management units (MUs) among populations of tiger quolls to assist in the conservation of these taxa. Ninety-three tiger quolls from nine representative populations were sampled from throughout the species range. Six nuclear microsatellite loci and the mitochondrial DNA (mtDNA) control region (471 bp) were used to examine ESUs and MUs in this species. We demonstrated that Tasmanian tiger quolls are reciprocally monophyletic to those from the mainland using mtDNA analysis, but D.m. gracilis was not monophyletic with respect to mainland D.m. maculatus. Analysis of microsatellite loci also revealed significant differences between the Tasmanian and mainland tiger quolls, and between D.m. gracilis and mainland D.m. maculatus. These results indicate that Tasmanian and mainland tiger quolls form two distinct evolutionary units but that D.m. gracilis and mainland D.m. maculatus are different MUs within the same ESU. The two marker types used in this study revealed different male and female dispersal patterns and indicate that the most appropriate units for short-term management are local populations. A revised classification and management plan are needed for tiger quolls, particularly in relation to conservation of the Tasmanian and Queensland populations.     

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake,the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long‐term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.     

PCR primers for microsatellite loci in the desert tortoise (Gopherus agassizii,Testudinidae)

The desert tortoise, Gopherus agassizii, is a threatened species native to the North American desert southwest and is recognized as having distinct Mojave and Sonoran populations. We identified six polymorphic microsatellite loci in the desert tortoise. All six loci were polymorphic in Sonoran samples. Five of the loci were variable in Mojave samples with varying degrees of amplification success. Two of the loci exhibited low allelic variation (2–3 alleles) while four were highly variable (8–27 alleles).     

Genetic Variation of Microsatellite Loci in the Major Histocompatibility Complex (MHC) Region in the Southeast Asian House Mouse (<Emphasis Type="Italic">Mus musculus castaneus</Emphasis>)

Major histocompatibility complex (MHC) genes are the most polymorphic loci known for vertebrates. Here we employed five microsatellite loci closely linked to the MHC region in an attempt to study the amount of genetic variation in 19 populations of the southeast Asian house mouse (Mus musculus castaneus) in Taiwan. The overall polymorphism at the five loci was high (H_e = 0.713), and the level of polymorphism varied from locus to locus. Furthermore, in order to investigate if selection is operating on MHC genes in natural mouse populations, we compared the extent and pattern of genetic variation for the MHC-linked microsatellite loci (the MHC loci) with those for the microsatellite loci located outside the MHC region (the non-MHC loci). The number of alleles and the logarithm of variance in repeat number were significantly higher for the MHC loci than for the non-MHC loci, presumably reflecting linkage to a locus under balancing selection. Although three statistical tests used do not provide support for selection, their lack of support may be due to low statistical power of the tests, to weakness of selection, or to a profound effect of genetic drift reducing the signature of balancing selection. Our results also suggested that the populations in the central and the southwestern regions of Taiwan might be one part of a metapopulation structure.     

Initial stages of reproductive isolation in two species of the endangered Sonoran topminnow

Long-term geographic isolation can result in reproductive incompatibilities due to forces such as mutation, genetic drift, and differential selection. In the Sonoran topminnow, molecular genetic studies of mtDNA, microsatellites, and MHC genes have shown that the endangered Gila and Yaqui topminnows are substantially different, suggesting that divergence took place approximately two million years ago. Here we examined hybrid crosses and backcrosses between these two allopatric taxa to evaluate the accumulation of postmating barriers to reproduction. These results are then compared with results from a previous study where male topminnows were shown to mate assortatively with conspecific females. Despite their preference for conspecific mates, both types of interspecific crosses successfully produced offspring. There was evidence of reduced hybrid fitness, including smaller mean brood size and male-biased sex ratio, for some classes of backcrosses. Brood sizes and interbrood intervals varied significantly when hybrids were subdivided into different cross categories. Our results illustrate the importance of distinctly defining hybrid classes in studies of reproductive isolation. To our knowledge, this is the first such detailed evolutionary analysis in endangered fish taxa.     

Conservation units based on mitochondrial and nuclear DNA variation among European bullhead populations (Cottus gobio L., 1758) from Flanders, Belgium

In this study, we aimed to delineateevolutionarily significant units (ESUs) andmanagement units (MUs) for the Europeanbullhead in Flanders (Belgium). Therefore, wedetermined the genetic interrelationshipsbetween 11 bullhead populations, using lengthvariation at 7 polymorphic microsatellite lociand sequence variation in the d-loop of themitochondrial DNA (mtDNA). Despite therelatively small geographical scale of ourstudy, the analysis of the d-loop sequencesshows that the Flemish bullhead populationscontain 3 haplotype groups, which can beassigned to 3 previously described EuropeanmtDNA clades. Because of the importantdifferences between these clades, they may bedefined as evolutionarily significant units,which should be managed separately. Analysis ofmicrosatellite data reveals very high degreesof isolation between populations, with theexception of 3 pairwise comparisons whichinvolved adjacent populations. Our data suggestthat the 3 haplotype groups probably qualify asESUs, as they show phylogeographicdifferentiation for mtDNA variants as well assignificant divergence of allele frequencies atnuclear loci. However, one of these units,limited to a single population, may be ofCentral European origin. All populations of theScheldt basin meet the criteria for MUrecognition, since significantly differentmicrosatellite allele frequencies as well asprivate alleles are found. In contrast, geneticdifferentiation among the 3 populations of theMeuse basin is very low.     

Contrasting patterns of variation in MHC loci in the Alpine newt

Major histocompatibility complex (MHC) genes are essential in pathogen recognition and triggering an adaptive immune response. Although they are the most polymorphic genes in vertebrates, very little information on MHC variation and patterns of evolution are available for amphibians, a group known to be declining rapidly worldwide. As infectious diseases are invoked in the declines, information on MHC variation should contribute to devising appropriate conservation strategies. In this study, we examined MHC variation in 149 Alpine newts ( Mesotriton alpestris ) from three allopatric population groups in Poland at the northeastern margin of the distribution of this species. The genetic distinctiveness of the population groups has previously been shown by studies of skin graft rejection, allozymes and microsatellites. Two putative expressed MHC II loci with contrasting levels of variation and clear evidence of gene conversion/recombination between them were detected. The Meal-DAB locus is highly polymorphic (37 alleles), and shows evidence of historical positive selection for amino acid replacements and substantial geographical differentiation in allelic richness. On the contrary, the Meal-DBB locus exhibits low polymorphism (three alleles differing by up to two synonymous substitutions) and a uniform distribution of three alleles among geographical regions. The uniform frequencies of the presumptively neutral Meal-DBB alleles may be explained by linkage to Meal-DAB . We found differences in allelic richness in Meal-DAB between regions, consistent with the hypothesis that genetic drift prevails with increasing distance from glacial refugia. Pseudogene loci appear to have evolved neutrally. The level of DAB variation correlated with variation in microsatellite loci, implying that selection and drift interplayed to produce the pattern of MHC variation observed in marginal populations of the Alpine newt.     

Genetic signatures of translocations and habitat fragmentation for two evolutionarily significant units of a protected fish species

Genetic population structure was evaluated for the White Sands pupfish (Cyprinodon tularosa), a protected fish species comprised of two Evolutionarily Significant Units (ESUs); the Malpais Spring ESU and the Salt Creek ESU. The Malpais Spring ESU is restricted to Malpais Spring, whereas the Salt Creek ESU includes the native Salt Creek population and two Salt Creek-derived populations at Mound Spring and Lost River; all three of these habitats are physically fragmented. We sampled the upper and lower reaches of the four populations, examining 13 DNA microsatellite loci from 40 individuals per population. As expected, significant genetic structure was observed between the two ESUs; Malpais Spring and Salt Creek. Substantial genetic drift was observed for the introduced Lost River population, with modest genetic drift for the introduced Mound Spring population. Taken together with ecological data, neither of the introduced populations successfully replicates the Salt Creek population. We also report significant reductions in genetic diversity for the upper reaches of both Salt Creek and Lost River, indicating that recent habitat changes have altered the genetic structure of these two populations. We consider these findings along with previously reported ecological data to develop guidelines for managing C. tularosa.     

Size homoplasy and mutational processes of interrupted microsatellites in two bee species, Apis mellifera and Bombus terrestris (Apidae)

Similar microsatellite electromorphs (PCR products of the same size) can arise from independent mutational events. Such alleles are not identical by descent. This phenomenon, termed size homoplasy, was studied by sequencing electromorphs of two microsatellite loci in which the stretch of basic repeats is interrupted by different short (1-2 bp) DNA motifs. The number and position of these interruptions were established for electromorphs from closely and distantly related populations of honeybees and bumblebees. No sequence difference was found when electromorphs came from the same subspecies or from closely related subspecies, suggesting that they were probably identical by descent. In contrast, sequence differences were often detected in distantly related subspecies, showing that size homoplasy frequently occurs at this level of population differentiation. Size homoplasy is increased by limits to free length variation of alleles, a phenomenon that seems to act on interrupted microsatellites when comparing distantly related taxa, that is, honeybee subspecies from different evolutionary lineages. Electromorph sequences suggest that, within the scope of these limits, large mutation events have occurred frequently at both interrupted loci studied. In good agreement with the molecular data, computations based on the observed heterozygosity and number of electromorphs and simulation studies showed that neither locus fits the one-step stepwise mutant model (SMM). We speculate that interrupted microsatellites in general could be characterized by a higher variance in repeat number and consequently a lower homoplasy rate than pure ones. Hence, interrupted microsatellites should be most appropriate for investigating population differentiation and evolutionary relationship between relatively distant populations.      

Founder effect in an island population of bighorn sheep

The Tiburon Island population of desert bighorn sheep has increased in size from 20 founders in 1975 to approximately 650 in 1999. This population is now the only population being used as the source stock for transplantations throughout northern Mexico. To evaluate the genetic variation in this population, we examined 10 microsatellite loci and a major histocompatibility complex (MHC) locus. The genetic variation was significantly less than found in other populations of the same subspecies in Arizona. Using a model that takes into account the effects of genetic drift on genetic distance, most of the genetic distance observed between the Tiburon population and Arizona samples could be explained. Because of the low genetic variation found in the Tiburon population, it is suggested that the Tiburon population should be supplemented with additional unrelated animals or that the transplant populations should be supplemented with unrelated animals.     

Population Genetic Patterns of Threatened European Mudminnow (Umbra krameri Walbaum, 1792) in a Fragmented Landscape: Implications for Conservation Management

The European mudminnow (Umbra krameri) is a Middle Danubian endemic fish species, which is characterised by isolated populations living mainly in artificial habitats in the centre of its range, in the Carpathian Basin. For their long term preservation, reliable information is needed about the structure of stocks and the level of isolation. The recent distribution pattern, and the population genetic structure within and among regions were investigated to designate the Evolutionary Significant, Conservation and Management Units (ESUs, CUs, MUs) and to explore the conservation biological value of the shrinking populations. In total, eight microsatellite loci were studied in 404 specimens originating from eight regions. The results revealed a pronounced population structure, where strictly limited gene flow was detected among regions, as well as various strengths of connections within regions. Following the results of hierarchical structure analyses, two ESUs were supposed in the Carpathian Basin, corresponding to the Danube and Tisza catchments. Our results recommend designating the borders of CUs in an 80–90km range and 16 clusters should be set up as MUs for the 33 investigated populations. How these genetic findings can be used to better allocate conservation resources for the long term maintenance of the metapopulation structure of this threathened endemic fish is discussed.     

Microsatellite identification of individual sockeye salmon in Barkley Sound, British Columbia

Population structure of sockeye salmon Oncorhynchus nerka in Barkley Sound, British Columbia, Canada was examined by analysis of microsatellite variation at 14 loci in three populations sampled in each of 3 years. The mean F _ST over all 14 loci was 0·063. Differences among populations accounted for 12 times the variation observed among years within populations. The number of alleles present at a locus was related to the power of the locus to provide accurate identification of individuals to population. The more alleles that were present at a locus, the greater was the power of the locus for individual identification. Individuals were correctly classified to one of three lakes of origin at a rate of 89%, and to one of two river drainages at a rate of 96%.     

Microsatellite Variation in Honey Bee (Apis Mellifera L.) Populations: Hierarchical Genetic Structure and Test of the Infinite Allele and Stepwise Mutation Models

Samples from nine populations belonging to three African (intermissa, scutellata and capensis) and four European (mellifera, ligustica, carnica and cecropia) Apis mellifera subspecies were scored for seven microsatellite loci. A large amount of genetic variation (between seven and 30 alleles per locus) was detected. Average heterozygosity and average number of alleles were significantly higher in African than in European subspecies, in agreement with larger effective population sizes in Africa. Microsatellite analyses confirmed that A. mellifera evolved in three distinct and deeply differentiated lineages previously detected by morphological and mitochondrial DNA studies. Dendrogram analysis of workers from a given population indicated that super-sisters cluster together when using a sufficient number of microsatellite data whereas half-sisters do not. An index of classification was derived to summarize the clustering of different taxonomic levels in large phylogenetic trees based on individual genotypes. Finally, individual population X loci data were used to test the adequacy of the two alternative mutation models, the infinite allele model (IAM) and the stepwise mutation models. The better fit overall of the IAM probably results from the majority of the microsatellites used including repeats of two or three different length motifs (compound microsatellites).     

Variation in MHC genotypes in two populations of house sparrow (Passer domesticus) with different population histories

Small populations are likely to have a low genetic ability for disease resistance due to loss of genetic variation through inbreeding and genetic drift. In vertebrates, the highest genetic diversity of the immune system is located at genes within the major histocompatibility complex (MHC). Interestingly, parasite‐mediated selection is thought to potentially maintain variation at MHC loci even in populations that are monomorphic at other loci. Therefore, general loss of genetic variation in the genome may not necessarily be associated with low variation at MHC loci. We evaluated inter‐ and intrapopulation variation in MHC genotypes between an inbred (Aldra) and a relatively outbred population (Hestmannøy) of house sparrows (Passer domesticus) in a metapopulation at Helgeland, Norway. Genomic (gDNA) and transcribed (cDNA) alleles of functional MHC class I and IIB loci, along with neutral noncoding microsatellite markers, were analyzed to obtain relevant estimates of genetic variation. We found lower allelic richness in microsatellites in the inbred population, but high genetic variation in MHC class I and IIB loci in both populations. This suggests that also the inbred population could be under balancing selection to maintain genetic variation for pathogen resistance.