Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis).

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.

Genetic variation across trophic levels: A test of the correlation between population size and genetic diversity in sympatric desert lizards.

Understanding the causes of genetic variation in real populations has been elusive. Competing theories claim that neutral vs. selective processes have a greater influence on the genetic variation within a population. A key difference among theories is the relationship between population size and genetic diversity. Our study tests this empirically by sampling two species of herbivorous lizards (Dipsosaurus dorsalis and Sauromalus ater) and two species of carnivorous lizards (Crotaphytus bicinctores and Gambelia wislizenii) that vary in population size at the same locality, and comparing metrics of genetic diversity. Contrary to neutral expectations, results from four independent loci showed levels of diversity were usually higher for species with smaller population sizes. This suggests that selective processes may be having an important impact on intraspecific diversity in this reptile community, although tests showed little evidence for selection on the loci sequenced for this study. It is also possible that idiosyncratic histories of the focal species may be overriding predictions from simple neutral models. If future studies show that lack of correlation between population size and genetic diversity is common, methods using genetic diversity to estimate population parameters like population size or time to common ancestor should be used with caution, as these estimates are based on neutral theory predictions.

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice.

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.

Analysis of quantitative trait locus effects on the size and shape of mandibular molars in mice.

While >50 genes have been found to influence the development of teeth in mice, we still know very little about the genetic basis for the adaptive characteristics of teeth, such as size and shape. We applied interval mapping procedures to Procrustes size and shape data obtained from 10 morphological landmarks on the mandibular molar row of the F(2) progeny from a cross between the LG/J and SM/J strains of mice. This revealed many more QTL for molar shape (18) than for molar centroid size (3), although levels of dominance effects were comparable among QTL for size and shape. Comparisons of patterns of Procrustes additive and dominance shape effects and ordination of QTL effects by principal components analysis suggested that the effects of the shape QTL were dispersed among the three molars and thus that none of these molars represents a genetically distinct developmental structure. The results of an analysis of co-occurrence of QTL for molar shape, mandible shape, and cranial dimensions in these mice suggested that many of the QTL for molar shape may be the same as those affecting these other sets of characters, although in some cases this could be due to effects of closely linked genes.

An epistatic genetic basis for fluctuating asymmetry of mandible size in mice.

The genetic basis of fluctuating asymmetry (FA), or nondirectional variation in the subtle differences between left and right sides of bilateral characters, continues to be of considerable theoretical interest. FA generally has been thought to arise from random noise during development and therefore to have a largely or entirely environmental origin. Whereas additive genetic variation for FA generally has been small and often insignificant, a number of investigators have hypothesized that interactions between loci, or epistasis, significantly influence FA. We tested this hypothesis by conducting a whole-genome scan to detect any epistasis in FA of centroid size in the mandibles of more than 400 mice from an F2 intercross population formed from crossing the Large (LG/J) and Small (SM/J) inbred strains. Genotypic deviations were imputed at each site 2 cM apart on all 19 autosomes, and these and centroid size asymmetry values were used in canonical correlation analyses for each of the 171 possible pairs of 19 autosomes to identify the most probable sites for epistasis. Epistasis for centroid size asymmetry was abundant, occurring far more often than was expected by chance alone (there were 30 separate instances of epistasis at the 0.001 significance level, when only two were expected by chance alone). The contributions of epistasis from 30 pairwise combinations of loci tended to suppress the additive and dominance genetic variance, but greatly increased the epistatic genetic variance for FA in centroid size given the intermediate allele frequencies of an F2 intercross population.

Quantitative Trait Loci for Early- and Late-Developing Skull Characters in Mice: A Test of the Genetic Independence Model of Morphological Integration.

Quantitative genetical theory suggests that characters existing in developmentally or functionally integrated groups are expected to be genetically correlated because they share a common inheritance. The genetic independence model for the cause of this genetic integration predicts that pleiotropic effects of single genes are mostly restricted to the characters in these integrated groups. We tested this model by estimating the additive and dominance effects of quantitative trait loci (QTLs) affecting early- (cranial vault) and late-developing (face) skeletal characters in F2 house mice originally derived from a cross of the Large and Small inbred strains. Interval mapping procedures were used that resulted in the identification of 26 QTLs on 17 of the 19 autosomes that significantly affected these characters. Additive, but not dominance, genotypic effects of many of these QTLs predominantly affected either the cranial vault or face characters, which supports the genetic independence model. Only two QTLs had positive pleiotropic effects on one group of characters but negative pleiotropic effects on the other (antagonistic pleiotropy).

Adaptive color polymorphism and unusually high local genetic diversity in the side-blotched lizard, Uta stansburiana.

Recently, studies of adaptive color variation have become popular as models for examining the genetics of natural selection. We examined color pattern polymorphism and genetic variation in a population of side-blotched lizards (Uta stansburiana) that is found in habitats with both dark (lava) and light colored (granite) substrates. We conducted a limited experiment for adult phenotypic plasticity in laboratory conditions. We recorded both substrate and lizard color patterns in the field to determine whether lizards tended to match their substrate. Finally we examined genetic variation in a gene (melanocortin 1 receptor) that has been shown to affect lizard color in other species and in a presumably neutral gene (mitochondrial cytochrome b). Populations were sampled in the immediate area of the lava flows as well as from a more distant site to examine the role of population structure. Our captive Uta did not change color to match their background. We show that side-blotched lizards tend to match the substrate on which it was caught in the field and that variation in the melanocortin 1 receptor gene does not correlate well with color pattern in this population. Perhaps the most remarkable result is that this population of side-blotched lizards shows extremely high levels of variation at both genetic markers, in the sense of allele numbers, with relatively low levels of between-allele sequence variation. Genetic variation across this small region was as great or greater than that seen in samples of pelagic fish species collected worldwide. Statistical analysis of genetic variation suggests rapid population expansion may be responsible for the high levels of variation.

Sperm development and motility are regulated by PP1 phosphatases in Caenorhabditis elegans.

Sperm from different species have evolved distinctive motility structures, including tubulin-based flagella in mammals and major sperm protein (MSP)-based pseudopods in nematodes. Despite such divergence, we show that sperm-specific PP1 phosphatases, which are required for male fertility in mouse, function in multiple processes in the development and motility of Caenorhabditis elegans amoeboid sperm. We used live-imaging analysis to show the PP1 phosphatases GSP-3 and GSP-4 (GSP-3/4) are required to partition chromosomes during sperm meiosis. Postmeiosis, tracking fluorescently labeled sperm revealed that both male and hermaphrodite sperm lacking GSP-3/4 are immotile. Genetic and in vitro activation assays show lack of GSP-3/4 causes defects in pseudopod development and the rate of pseudopodial treadmilling. Further, GSP-3/4 are required for the localization dynamics of MSP. GSP-3/4 shift localization in concert with MSP from fibrous bodies that sequester MSP at the base of the pseudopod, where directed MSP disassembly facilitates pseudopod contraction. Consistent with a role for GSP-3/4 as a spatial regulator of MSP disassembly, MSP is mislocalized in sperm lacking GSP-3/4. Although a requirement for PP1 phosphatases in nematode and mammalian sperm suggests evolutionary conservation, we show PP1s have independently evolved sperm-specific paralogs in separate lineages. Thus PP1 phosphatases are highly adaptable and employed across a broad range of sexually reproducing species to regulate male fertility.

Phylogeography of the diamond turbot (Hypsopsetta guttulata) across the Baja California Peninsula.

We compared morphology and sequenced nuclear and mitochondrial genes from 11 populations of a previously genetically unstudied "Baja California disjunct" species, the diamond turbot (Hypsopsetta guttulata). This species exhibits very limited adult movement and restriction to soft-bottom habitats but has a moderately long pelagic larval duration. Therefore, if pelagic larval duration is correlated with gene flow between Gulf of California and Pacific populations, we expect a reduced level of genetic and morphological differentiation. However, if adult habitat and ecology have more effect on gene flow, we expect the populations in the two bodies of water to be more highly differentiated. We used logistic regression to compare morphological features and phylogenetic and population genetic analyses to compare nucleotide sequence data. Gulf of California H. guttulata are different from Pacific populations in morphology and both mitochondrial and nuclear gene sequences. MtDNA shows reciprocal monophyly, and nuclear sequences from the Gulf of California formed a monophyletic group. Population genetic analyses also suggest further population subdivision within the Pacific and within the Gulf of California. We argue that adult ecology has a significant effect on migration rates among populations in the Pacific Ocean and the Gulf of California.

Cytonuclear discordance across a leopard frog contact zone.

We conducted a genetic analysis of the extent of hybridization within and outside a contact zone between two North American leopard frogs, Rana blairi and Rana pipiens by comparing distribution patterns of mitochondrial and nuclear haplotypes. The contact zone, located between South Dakota, Nebraska, and Iowa, USA, was previously examined using morphological data in the 1970s, leading to the conclusion that hybridization was rare between R. pipiens and R. blairi. Our genetic analysis of 51 populations (611 samples) shows strong cytonuclear discordance. Mitochondrial-haplotype distribution matches the same pattern as the documented species spatial distributions based on morphology. However, the geographic distribution of the nuclear haplotypes reveals asymmetrical swamping of the R. pipiens nuclear haplotypes by R. blairi haplotypes. Phylogenetic analyses of both mitochondrial and nuclear markers provide strong evidence for the presence of R. blairi-R. pipiens introgression for the nuclear marker. A pattern of mitochondrial isolation with nuclear introgression is extremely unusual, and predicted to occur much less often than the reverse.

GENE EFFECTS ON A QUANTITATIVE TRAIT: TWO-LOCUS EPISTATIC EFFECTS MEASURED AT MICROSATELLITE MARKERS AND AT ESTIMATED QTL.

Most evolutionarily and agriculturally important traits are affected by many genes (quantitative trait loci, or QTL) of relatively small effect. Usually the genetics of these traits are examined by indirect statistical analysis of the covariance among relatives, rather than by direct analyses. We use new analytical and molecular techniques to examine nonadditive interactions of microsatellite markers and estimated QTL that influence adult body weight in mice. Offspring of a cross between a large inbred mouse strain (LG/J) and a small inbred strain (SM/J) were intercrossed to form a segregating F2 generation. Using 76 microsatellite markers and 19 estimated QTL, we estimate gene-level epistasis and population-level epistasis for body weight at 10 weeks for 534 F2 mice. Significant epistasis was found for large numbers of the two locus comparisons using both markers and previously detected QTL. There are many genes segregating for adult body weight in this cross and many of these genes appear to interact epistatically. The discovery of potentially extensive epistasis has important implications for evolutionary models.

EPISTASIS AS A SOURCE OF INCREASED ADDITIVE GENETIC VARIANCE AT POPULATION BOTTLENECKS.

The role of epistasis in evolution and speciation has remained controversial. We use a new parameterization of physiological epistasis to examine the effects of epistasis on levels of additive genetic variance during a population bottleneck. We found that all forms of epistasis increase average additive genetic variance in finite populations derived from initial populations with intermediate allele frequencies. Average additive variance continues to increase over many generations, especially at larger population sizes (N = 32 to 64). Additive-by-additive epistasis is the most potent source of additive genetic variance in this situation, whereas dominance-by-dominance epistasis contributes smaller amounts of additive genetic variance. With additive-by-dominance epistasis, additive genetic variance decreases at a relatively high rate immediately after a population bottleneck, rebounding to higher levels after several generations. Empirical examples of epistasis for murine adult body weight based on measured genotypes are provided illustrating the varying effects of epistasis on additive genetic variance during population bottlenecks.

PLEIOTROPIC EFFECTS OF INDIVIDUAL GENE LOCI ON MANDIBULAR MORPHOLOGY.

The genotypic basis of morphological variation is largely unknown. In this study we examine patterns of pleiotropic effects on mandibular morphology at individual gene loci to determine whether the pleiotropic effects of individual genes are restricted to functionally and developmentally related traits. Mandibular measurements were obtained from 480 mice from the F2 generation of an intercross between the LG/J and SM/J mouse strains. DNA was also extracted from these animals, and 76 microsatellite loci covering the autosomes were scored. Interval mapping was used to detect chromosomal locations with significant effects on various mandibular measurements. Sets of traits mapping to a common chromosomal region were considered as being affected by a single quantitative trait locus (QTL) for mandibular morphology. Thirty-seven such chromosomal regions were identified spread throughout the autosomes. Gene effects were small to moderate with the allele derived from the LG/J strain typically leading to larger size. When dominance was present, the LG/J allele was typically dominant to the SM/J allele. Most loci affected restricted functional and developmental regions of the mandible. Of the 26 chromosomal regions affecting more than two traits, 50% affect the muscular processes of the ascending ramus, 27% affect the alveolar processes carrying the teeth, and 23% affect the whole mandible. Four additional locations affecting two traits had effects significantly associated with alveolar regions. Pleiotropic effects are typically restricted to morphologically integrated complexes.

A SEARCH FOR QUANTITATIVE TRAIT LOCI AFFECTING ASYMMETRY OF MANDIBULAR CHARACTERS IN MICE.

An interval mapping procedure was used to search for and describe the effects of any quantitative trait loci (QTLs) for directional asymmetry (DA) and fluctuating asymmetry (FA) of 10 bilateral mandible characters in house mice. It was hypothesized that more QTLs would be found for DA than for FA, but that any discovered for FA should tend to exhibit dominance. All mandible characters were triply measured and 76 microsatellite markers were scored in an average of 471 mice from the F2 intercross of the Large (LG/J) and Small (SM/J) inbred strains. A total of 16 QTLs significantly affected DA in nine of the 10 mandible characters, and this was more than the 9.5 expected by chance alone. These QTLs were found on seven of the 19 chromosomes, often at or near locations of QTLs affecting the mean of the two sides for various dimensions on the mandible. It was concluded that there is genetical variability for DA in these characters, although its level was low (4.4% of the total variation in this particular F2 population). Eleven QTLs were detected for FA, suggesting that there is very little genetic variability for FA, at least as seen in the mandible characters in this particular F2 population. As hypothesized, however, these QTLs did tend to exhibit dominance.

Molecular systematics, hybridization, and phylogeography of the Bufo americanus complex in Eastern North America.

We reconstruct phylogenetic relationships among a well-studied group of toads and find relationships that differ greatly from the current taxonomic understanding. We use mitochondrial sequences encoding ND1, tRNA(Leu(UUR)), and part of 16S to infer relationships among members of the Bufo americanus complex. Focusing on the four taxa that historically have been most problematic due to morphological similarity and hybridization in sympatry, we sample 150 individuals from multiple populations across each species' geographic range. Our evidence conflicts with previous taxonomic hypotheses that were based on ability to hybridize, geographic distribution, and call variation. First, sequences from B. fowleri do not comprise the sister clade to sequences of B. woodhousii; therefore the previous classifications of B. fowleri as sister species to, or eastern subspecies of, B. woodhousii are both called into question. Second, sequences from B. americanus are more closely related to those of B. woodhousii than to those of B. terrestris, indicating that similar advertisement call characteristics evolved independently. Third, sequences of B. fowleri are paraphyletic, with sequences of B. terrestris embedded within. Lastly, sequences from B. fowleri cluster into three distinct mitochondrial clades, with some divergences corresponding to greater than 2mya. These clades are somewhat geographically structured, suggesting divergence in allopatry during the Pleistocene. These mitochondrial divergences are not accompanied by known phenotypic differences, however, suggesting either evolutionary stasis in morphology and behavior, cryptic phenotypic evolution, or that hybridization in secondary contact has homogenized phenotypic differences that may have arisen in allopatry.