Relationship between single-cross performance and molecular marker heterozygosity

Summary In studies involving isozymes or restriction fragment length polymorphisms (RFLPs), correlations of parental molecular marker diversity with grain yield of maize (Zea mays L.) single-crosses have been too low to be of any predictive value. The relationship of molecular marker heterozygosity (D_ij) with hybrid performance ( _ij) and combining ability was examined. For a simple genetic model involving uncorrelated parental allele frequencies and complete coverage of quantitative trait loci (QTL) by molecular markers, the correlations between _ij and D_ij were 0.25. _ij and D_ij were partitioned into general and specific effects. The expected correlation between specific combining ability and specific molecular marker heterozygosity is high. Expected correlations between general combining ability and general molecular marker heterozygosity are either positive or negative, depending on allele frequencies in the tester lines. Computer simulation was used to investigate a more complex but more realistic genetic model involving incomplete coverage of QTL by molecular markers. All of the following conditions are necessary for effective prediction of hybrid performance based on molecular marker heterozygosity: (1) dominance effects are strong; (2) allele frequencies at individual loci in the parental inbreds are negatively correlated; (3) trait heritability is high; (4) average parental allele frequencies vary only within a narrow range; (5) at least 30–50% of the QTL are linked to molecular markers; and (6) not more than 20–30% of the molecular markers are randomly dispersed or unlinked to QTL.A contribution from Limagrain Genetics, a Groupe Limagrain company     

CYTONUCLEAR THEORY FOR HAPLODIPLOID SPECIES AND X-LINKED GENES. II. STEPPING-STONE MODELS OF GENE FLOW AND APPLICATION TO A FIRE ANT HYBRID ZONE

We develop cytonuclear, hybrid zone models for haplodiploid species or X-linked genes in diploid species using a stepping-stone framework of migration, in which migration rates vary with both direction and sex. The equilibrium clines for the allele frequencies, cytonuclear disequilibria, and frequencies of pure parental types are examined for species with diagnostic markers, under four important migration schemes: uniform migration of both sexes in both directions, greater migration of both sexes from one direction, greater migration of females, and greater migration of males. Of the three cytonuclear variables examined, the allele frequency clines are the most informative in differentiating among the various migration patterns. The cytonuclear disequilibria and the frequency of the pure parental types tend to be useful only in revealing directional asymmetries in migration. The extent of hybrid zone subdivision has quantitative but not qualitative effects on the distribution of cytonuclear variables, in that the allele frequency clines become more gradual, the cytonuclear disequilibria decrease in magnitude, and the frequencies of pure parentals decline with increasing subpopulation number. Also, the only major difference between the X-linked and haplodiploid frameworks is that a higher frequency of pure parentals is found when considering haplodiploids, in which male production does not require mating. The final important theoretical result is that censusing after migration yields greater disequilibria and parental frequencies than censusing after mating. We analyzed cytonuclear data from two transects from a naturally occurring hybrid zone between two haplodiploid fire ant species, Solenopsis invicta and S. richteri, using our stepping-stone framework. The frequency of S. invicta mtDNA exceeds the frequency of the S. invicta nuclear markers through much of this hybrid zone, indicating that sex differences in migration or selection may be occurring. Maximum-likelihood estimates for the migration rates are very high, due to an unexpectedly large number of pure parental types in the hybrid zone, and differ substantially between the two transects. Overall, our model does not provide a good fit, in part because the S. invicta–S. richteri hybrid zone has not yet reached equilibrium.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Populations composed entirely of hybrid colonies: bidirectional hybridization and polyandry in harvester ants

In eusocial Hymenoptera, haplodiploid life cycles, obligate sterile castes, and polyandry may facilitate selection for hybridization. We analyzed a broad hybrid zone between the ecologically distinct seed‐harvester ants Pogonomyrmex occidentalis (Cresson) and Pogonomyrmex maricopa (Wheeler) using mitochondrial (mt)DNA sequence data, eight morphological markers, and 14 random amplified polymorphic DNA (RAPD) markers. Average mtDNA sequence divergence among parental species was 11.34%, indicating secondary contact. RAPD markers were significantly correlated with morphological variation, confirming the interspecific hybrid origin of all morphologically putative hybrid colonies. A morphological hybrid index indicates an abundance of both F₁ hybrids and parental morphotypes within colonies. Individual character frequencies plotted against distance show coincident and concordant clines, suggesting little to no introgression. The structure of the hybrid zone is two‐fold. Within the western region, stark reversals in character frequencies coincide with overt soil differences, indicating a mosaic hybrid zone structure. The eastern region is a riparian habitat where four adjacent populations were composed entirely of hybrid colonies. These habitat associations suggest that hybrid worker genomes permit dispersal into intermediate environments that select against one or both parental species. The present study suggests that, in addition to retaining reproductive compatibility, ecologically distinct species of ants may generate hybrid colonies maintained by environmental selection. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 320–336.     

Genetic and morphological evidence of a geographically widespread hybrid zone between two crocodile species,Crocodylus acutus and Crocodylus moreletii

Hybrid zones represent natural laboratories to study gene flow, divergence and the nature of species boundaries between closely related taxa. We evaluated the level and extent of hybridization between Crocodylus moreletii and Crocodylus acutus using genetic and morphological data on 300 crocodiles from 65 localities. To our knowledge, this is the first genetic study that includes the entire historic range and sympatric zone of the two species. Contrary to expectations, Bayesian admixture proportions and maximum‐likelihood estimates of hybrid indexes revealed that most sampled crocodiles were admixed and that the hybrid zone is geographically extensive, extending well beyond their historical region of sympatry. We identified a few geographically isolated, nonadmixed populations of both parental species. Hybrids do not appear to be F₁s or recent backcrosses, but rather are more likely later‐generation hybrids, suggesting that hybridization has been going on for several to many generations and is mostly the result of natural processes. Crocodylus moreletii is not the sister species of C. acutus, suggesting that the hybrid zone formed from secondary contact rather than primary divergence. Nonadmixed individuals from the two species were distinguishable based on morphological characters, whereas hybrids had a complex mosaic of morphological characters that hinders identification in the wild. Very few nonadmixed C. acutus and C. moreletii populations exist in the wild. Consequently, the last nonadmixed C. moreletii populations have become critically endangered. Indeed, not only the parental species but also the naturally occurring hybrids should be considered for their potential conservation value.     

Tandem repeat markers for population genetic studies of the protected California tiger salamander (<Emphasis Type="Italic">Ambystoma californiense</Emphasis>)

Habitat loss is the single greatest threat to persistence of the critically threatened California tiger salamander (Ambystoma californiense). To aid management plans that designate critical habitat for this species, I developed and characterized 21 tetranucleotide microsatellite markers using two native populations in Santa Barbara and Alameda Counties. Allelic variation and average heterozygosities were lower in the endangered Santa Barbara population (allele range 1–4, mean 2.4; H _O = 0.308 H _E = 0.288) compared with the threatened Alameda population (allele range 2–10, mean 6.7; H _O = 0.712, H _E = 0.722). In-depth population studies using these markers will provide vital information for plans to assign critical habitat that optimize gene flow among breeding populations, as well as for identifying non-native hybrid genotypes that threaten native A. californiense stocks. Beyond the conservation goals for A. californiense, the close phylogenetic relationships within the tiger salamander complex also suggest a broad utility for population studies using these markers.     

Widespread hybridization within mound‐building wood ants in Southern Finland results in cytonuclear mismatches and potential for sex‐specific hybrid breakdown

Hybridization and gene flow between diverging lineages are increasingly recognized as common evolutionary processes, and their consequences can vary from hybrid breakdown to adaptive introgression. We have previously found a population of wood ant hybrids between Formica aquilonia and F. polyctena that shows antagonistic effects of hybridization: females with introgressed alleles show hybrid vigour, whereas males with the same alleles show hybrid breakdown. Here, we investigate whether hybridization is a general phenomenon in this species pair and analyse 647 worker samples from 16 localities in Finland using microsatellite markers and a 1200‐bp mitochondrial sequence. Our results show that 27 sampled nests contained parental‐like gene pools (six putative F. polyctena and 21 putative F. aquilonia) and all remaining nests (69), from nine localities, contained hybrids of varying degrees. Patterns of genetic variation suggest these hybrids arise from several hybridization events or, instead, have backcrossed to the parental gene pools to varying extents. In contrast to expectations, the mitochondrial haplotypes of the parental species were not randomly distributed among the hybrids. Instead, nests that were closer to parental‐like F. aquilonia for nuclear markers preferentially had F. polyctena's mitochondria and vice versa. This systematic pattern suggests there may be underlying selection favouring cytonuclear mismatch and hybridization. We also found a new hybrid locality with strong genetic differences between the sexes similar to those predicted under antagonistic selection on male and female hybrids. Further studies are needed to determine the selective forces that act on male and female genomes in these newly discovered hybrids.     

Phenotypic and genetic introgression across a moving woodpecker hybrid zone

In hybrid zones in which two divergent taxa come into secondary contact and interbreed, selection can maintain phenotypic diversity despite widespread genetic introgression. Red‐breasted (Sphyrapicus ruber) and red‐naped (S. nuchalis) sapsuckers meet and hybridize along a narrow contact zone that stretches from northern California to southern British Columbia. We found strong evidence for changes in the structure of this hybrid zone across time, with significant temporal shifts in allele frequencies and in the proportions of parental phenotypes across the landscape. In addition to these shifts, we found that differences in plumage predict genetic differences (R² = 0.80), suggesting that plumage is a useful proxy for assessing ancestry. We also found a significant bimodal distribution of hybrids across the contact zone, suggesting that premating barriers may be driving reproductive isolation, perhaps as a result of assortative mating based on plumage differences. However, despite evidence of selection and strong patterns of population structure between parental samples, we found only weak patterns of genetic divergence. Using museum specimens and genomic data, this study of sapsuckers provides insight into the ways in which phenotypic and genetic structure have changed over a 40‐year period, as well as insight into the mechanisms that may contribute to the maintenance of the hybrid zone over time.     

Population estimators or progeny tests: what is the best method to assess null allele frequencies at SSR loci?

Nuclear SSRs are notorious for having relatively high frequencies of null alleles, i.e. alleles that fail to amplify and are thus recessive and undetected in heterozygotes. In this paper, we compare two kinds of approaches for estimating null allele frequencies at seven nuclear microsatellite markers in three French Fagus sylvatica populations: (1) maximum likelihood methods that compare observed and expected homozygote frequencies in the population under the assumption of Hardy-Weinberg equilibrium and (2) direct null allele frequency estimates from progeny where parent genotypes are known. We show that null allele frequencies are high in F. sylvatica (7.0% on average with the population method, 5.1% with the progeny method), and that estimates are consistent between the two approaches, especially when the number of sampled maternal half-sib progeny arrays is large. With null allele frequencies ranging between 5% and 8% on average across loci, population genetic parameters such as genetic differentiation (F _ST) may be mostly unbiased. However, using markers with such average prevalence of null alleles (up to 15% for some loci) can be seriously misleading in fine scale population studies and parentage analysis.     

Identification of cyprinid hybrids by using geometric morphometrics and microsatellites

Traditional morphological studies need to be complemented with modern genetic methods to facilitate the identification of hybrids. By using a combination of landmark‐based techniques and microsatellite markers, natural hybridization was investigated between the cyprinids Abramis brama (L.) and Blicca bjoerkna (L.) and their hybrids. Geometric morphometrics revealed significant differences in body shape between A. brama, B. bjoerkna, and hybrids. Hybrids were of intermediate body shape with a tendency of being more like A. brama. Genetic differentiation was found between both parental species and their hybrids. However, hybrids revealed a higher genetical similarity with A. brama. Based on sequencing of the mitochondrial ATP synthase subunit 6 and 8 region a clear split was found between the two sibling species. Seventeen out of 19 hybrid specimens clustered within the A. brama clade. Data indicate that hybridization between A. brama and B. bjoerkna is mainly unidirectional and has not yet resulted in fusion of the two parental gene pools. Genetic integrity is maintained in B. bjoerkna, but F₁ hybrid backcrosses might lead to introgression into the genepool of A. brama.     

Molecular patterns of introgression in a classic hybrid zone between the Australian tree frogs,Litoria ewingii and L. paraewingi: evidence of a tension zone

Hybrid zones provide a rare opportunity to explore the processes involved in reproductive isolation and speciation. The southern hybrid zone between the southeastern Australian tree frogs Litoria ewingii and L. paraewingi has been comprehensively studied over the last 40 years, primarily using reproductive compatibility experiments and male advertisement calls. We used mitochondrial DNA (mtDNA) and eight nuclear microsatellite markers to characterize this hybrid zone along a historically studied transect and to test various dispersal‐dependent and dispersal‐independent hybrid zone models. The species are genetically distinct and the level of hybridization within the contact zone is low, with the majority of admixed individuals representing later‐generation hybrids. Based on previous experimental genetic compatibility studies, we predicted that hybrids with L. paraewingi mtDNA would be more frequent than hybrids with L. ewingii mtDNA. Surprisingly, a greater proportion of the identified hybrids had L. ewingii mtDNA. Geographical cline analyses showed a sharp transition in allele frequencies across the transect, and both the mtDNA and microsatellite data showed concordant cline centres, but were best supported by a model that allowed width to vary. Overall, the L. ewingii–L. paraewingi hybrid zone is best characterized as a tension zone, due to the narrow cline width, concordant genetic clines and low levels of hybridization.     

Identification of contemporary selection signatures using composite log likelihood and their associations with marbling score in Korean cattle

Positive selection not only increases beneficial allele frequency but also causes augmentation of allele frequencies of sequence variants in close proximity. Signals for positive selection were detected by the statistical differences in subsequent allele frequencies. To identify selection signatures in Korean cattle, we applied a composite log‐likelihood (CLL)‐based method, which calculates a composite likelihood of the allelic frequencies observed across sliding windows of five adjacent loci and compares the value with the critical statistic estimated by 50 000 permutations. Data for a total of 11 799 nucleotide polymorphisms were used with 71 Korean cattle and 209 foreign beef cattle. As a result, 147 signals were identified for Korean cattle based on CLL estimates (P < 0.01). The signals might be candidate genetic factors for meat quality by which the Korean cattle have been selected. Further genetic association analysis with 41 intragenic variants in the selection signatures with the greatest CLL for each chromosome revealed that marbling score was associated with five variants. Intensive association studies with all the selection signatures identified in this study are required to exclude signals associated with other phenotypes or signals falsely detected and thus to identify genetic markers for meat quality.     

ECOLOGICAL AND GENETIC ASSOCIATIONS IN AN IRIS HYBRID ZONE

Chloroplast DNA (cpDNA) markers and 12 nuclear (random amplified polymorphic DNA, or RAPD) markers were used to examine the distribution of genetic variation among individuals and the genetic and ecological associations in a hybrid iris population. Plants in the population occurred at various distances from the edge of a bayou in a relatively undisturbed mixed hardwood forest and in an adjacent pasture dominated by herbaceous perennials with interspersed oak and cypress trees. The majority of plants sampled possessed combinations of markers from the different Iris species. Genetic markers diagnostic for Iris fulva and I. brevicaulis occurred at high frequencies, whereas markers diagnostic for I. hexagona were infrequent. For the majority of the nuclear markers, significant levels of cytonuclear disequilibria existed because of intraspecific associations among the markers in both the pasture and the forest. The distribution of nuclear markers among individuals was bimodal; intermediate genotypes were absent and the majority of RAPD markers were associated with their intraspecific cpDNA haplotypes. Strong intraspecific associations existed among RAPD markers in the forest, but associations tended to be weaker in the pasture area. Ecological correlations were detected for all but one of the I. fulva and I. brevicaulis RAPD markers. The ecological associations of hybrids similar to I. brevicaulis resembled associations of I. brevicaulis parental genotypes, suggesting that these hybrid genotypes may be relatively fit in the same habitats. The hybrids similar to I. fulva, however, were distributed in habitats that were unique relative to the parental species. The patterns of genetic and environmental associations along with other available data suggest that (1) only advanced generation hybrids were present in the population; (2) formation of F₁ hybrids among Louisiana irises is rare, leading to sporadic formation of hybrid populations; and (3) selection and assortative mating have contributed to the formation of hybrid genotypes that tend to be similar to parental genotypes. The patterns of ecological and genetic associations detected in this population suggest that assortative mating and environmental and viability selection are important in the structuring and maintenance of this hybrid zone.     

Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus)

Local adaptation is often invoked to explain hybrid zone structure, but empirical evidence of this is generally rare. Hybrid zones between two poeciliid fishes, Xiphophorus birchmanni and X. malinche, occur in multiple tributaries with independent replication of upstream‐to‐downstream gradients in morphology and allele frequencies. Ecological niche modelling revealed that temperature is a central predictive factor in the spatial distribution of pure parental species and their hybrids and explains spatial and temporal variation in the frequency of neutral genetic markers in hybrid populations. Among populations of parentals and hybrids, both thermal tolerance and heat‐shock protein expression vary strongly, indicating that spatial and temporal structure is likely driven by adaptation to local thermal environments. Therefore, hybrid zone structure is strongly influenced by interspecific differences in physiological mechanisms for coping with the thermal environment.     

Comparing the van Oosterhout and Chybicki‐Burczyk methods of estimating null allele frequencies for inbred populations

In spite of the usefulness of codominant markers in population genetics, the existence of null alleles raises challenging estimation issues in natural populations that are characterized by positive inbreeding coefficients (F > 0). Disregarding the possibility of F > 0 in a population will generally lead to overestimates of null allele frequencies. Conversely, estimates of inbreeding coefficients (F) may be strongly biased upwards (excess homozygotes), in the presence of nontrivial frequencies of null alleles. An algorithm has been presented for the estimation of null allele frequencies in inbred populations (van Oosterhout method), using external estimates of the F‐statistics. The goal of this study is to introduce a modification of this method and to provide a formal comparison with an alternative likelihood‐based method (Chybicki‐Burczyk). Using simulated data, we illustrate the strengths and limitations of these competing methods. Under most circumstances, the likelihood method is preferable, but for highly inbred organisms, a modified van Oosterhout method offers some advantages.     

Variable post‐zygotic isolation in Drosophila melanogaster/D. simulans hybrids

The study of hybrid inviability reveals cryptic divergence between the genetic interactions that maintain stable phenotypes in the pure species . We characterized the effects of natural variation on the penetrance of hybrid inviability phenotypes in crosses between Drosophila melanogaster and two species of the D. simulans subcomplex, D. simulans and D. sechellia. Using a panel of wild‐caught lines, we studied the levels of genetic variance present in D. simulans and D. sechellia affecting prezygotic and post‐zygotic isolation in hybridizations with D. melanogaster females. We observed extensive variability in the viability of hybrid individuals, dependent on the genotype of the parents, suggesting that intraspecific natural variation manifests directly in hybrid phenotypes. Furthermore, we found that genetic background significantly affects the penetrance of a well‐studied determinant of hybrid inviability: the interaction between Hmr^mel–Lhr^sim. Our results suggest that hybrid inviability – and reproductive isolation generally – can be modified by polymorphisms at multiple loci segregating within the parental species. Just as the penetrance of most mutant phenotypes can be modified by the genetic background within the pure species, the penetrance of hybrid inviability phenotypes is highly influenced by the parental genotypes.     

Microsatellite mutation rates in the eastern tiger salamander (<Emphasis Type="Italic">Ambystoma tigrinum tigrinum</Emphasis>) differ 10-fold across loci

Microsatellites are commonly used for mapping and population genetics because of their high heterozygosities and allelic variability (i.e., polymorphism). Microsatellite markers are generally more polymorphic than other types of molecular markers such as allozymes or SNPs because the insertions/deletions that give rise to microsatellite variability are relatively common compared to nucleotide substitutions. Nevertheless, direct evidence of microsatellite mutation rates (MMRs) is lacking in most vertebrate groups despite the importance of such estimates to key population parameters (e.g., genetic differentiation or θ = 4N _e μ). Herein, we present empirical data on MMRs in eastern tiger salamanders (Ambystoma tigrinum tigrinum). We conducted captive breeding trials and genotyped over 1,000 offspring at a suite of microsatellite loci. These data on 7,906 allele transfers provide the first direct estimates of MMRs in amphibians, and they illustrate that MMRs can vary by more than an order of magnitude across loci within a given species (one locus had ten mutations whereas the others had none).     

Ambiguous species boundaries: Hybridization and morphological variation in two closely related Rubus species along altitudinal gradients

Although hybridization frequently occurs among plant species, hybrid zones of divergent lineages formed at species boundaries are less common and may not be apparent in later generations of hybrids with more parental‐like phenotypes, as a consequence of backcrossing. To determine the effects of dispersal and selection on species boundaries, we compared clines in leaf traits and molecular hybrid index along two hybrid zones on Yakushima Island, Japan, in which a temperate (Rubus palmatus) and subtropical (Rubus grayanus) species of wild raspberry are found. Leaf sinus depth in the two hybrid zones had narrower clines at 600 m a.s.l. than the molecular hybrid index and common garden tests confirmed that some leaf traits, including leaf sinus depth that is a major trait used in species identification, are genetically divergent between these closely related species. The sharp transition in leaf phenotypic traits compared to molecular markers indicated divergent selection pressure on the hybrid zone structure. We suggest that species boundaries based on neutral molecular data may differ from those based on observed morphological traits.     

Genetic discrimination of middle Mississippi River <Emphasis Type="Italic">Scaphirhynchus</Emphasis> sturgeon into pallid,shovelnose, and putative hybrids with multiple microsatellite loci

The pallid sturgeon (Scaphirhynchus albus), which is protected under the US endangered species act, and shovelnose sturgeon (S. platorhynchus), which is legally harvested in some locations, are sympatric throughout the range of pallid sturgeon. There is considerable morphological overlap between the species making discrimination problematic. The inability to reliably differentiate between species across all life stages has hampered pallid sturgeon recovery efforts. Furthermore, the two species are believed to hybridize. This study used allele frequency data at multiple microsatellite loci to perform Bayesian and likelihood-based assignment testing and morphological measures and meristics to discriminate pallid, shovelnose, and putative hybrid sturgeons from the middle Mississippi River. Bayesian model-based clustering of the genetic data indicated that two natural genetic units occur in the region. These units correspond to morphologically identified pallid and shovelnose sturgeon. Some individuals were morphologically intermediate and many of these failed to strongly assign genetically as either pallid or shovelnose sturgeon, suggesting they may be hybrids. These data indicate that pallid sturgeon and shovelnose sturgeon are genetically distinct in the middle Mississippi River (F _ST = 0.036, P < 0.0001) and suggest that hybridization between pallid sturgeon and shovelnose sturgeon has occurred in this region with genetic distance estimates indicating the greatest distance is between pallid and shovelnose sturgeon, while hybrid sturgeon are intermediate but closer to shovelnose. This study demonstrates that assignment testing with multiple microsatellite markers can be successful at discriminating pallid sturgeon and shovelnose sturgeon, providing a valuable resource for pallid sturgeon recovery and conservation.