Differential strength of sex-biased hybrid inferiority in impeding gene flow may be a cause of Haldane's rule

Abstract In animals, if one sex of the F₁ hybrid between two species is sterile or inviable, it is usually the heterogametic (XY or WZ) sex. This phenomenon, known as Haldane's rule, is currently thought to be coincidentally caused by different mechanisms in separate entities. The following questions have never been asked: Are heterogametic and homogametic inferiority (sterility or inviability) equivalent as isolating mechanisms? Could discrepancies between them, if existing, produce Haldane's rule? Here I consider sex‐biased hybrid inferiority strictly as an isolating mechanism, and quantitatively evaluate its strength in impeding gene flow. The comparison reveals that the ability of sex‐biased inferiority to impede gene flow varies according to the sex and chromosome involved. Heterogametic inferiority is a weaker barrier when unidirectional and a much stronger one when in compound reciprocal directions, compared with homogametic inferiority. Such differential strength may affect divergence in speciation and produce Haldane's rule.     

Endosperm-based incompatibilities in hybrid monkeyflowers

Endosperm is an angiosperm innovation central to their reproduction whose development, and thus seed viability, is controlled by genomic imprinting, where expression from certain genes is parent-specific. Unsuccessful imprinting has been linked to failed inter-specific and inter-ploidy hybridization. Despite their importance in plant speciation, the underlying mechanisms behind these endosperm-based barriers remain poorly understood. Here, we describe one such barrier between diploid Mimulus guttatus and tetraploid Mimulus luteus. The two parents differ in endosperm DNA methylation, expression dynamics, and imprinted genes. Hybrid seeds suffer from underdeveloped endosperm, reducing viability, or arrested endosperm and seed abortion when M. guttatus or M. luteus is seed parent, respectively, and transgressive methylation and expression patterns emerge. The two inherited M. luteus subgenomes, genetically distinct but epigenetically similar, are expressionally dominant over the M. guttatus genome in hybrid embryos and especially their endosperm, where paternal imprints are perturbed. In aborted seeds, de novo methylation is inhibited, potentially owing to incompatible paternal instructions of imbalanced dosage from M. guttatus imprints. We suggest that diverged epigenetic/regulatory landscapes between parental genomes induce epigenetic repatterning and global shifts in expression, which, in endosperm, may uniquely facilitate incompatible interactions between divergent imprinting schemes, potentially driving rapid barriers.

Diverged epigenetic/regulatory landscapes between parental genomes result in epigenetic repatterning in hybrids that drive global shifts in endosperm gene expression patterns.     

Reinforcement and the genetics of hybrid incompatibilities

Recent empirical studies suggest that genes involved in speciation are often sex-linked. We derive a general analytic model of reinforcement to study the effects of sex linkage on reinforcement under three forms of selection against hybrids: one-locus, two-locus, and ecological incompatibilities. We show that the pattern of sex linkage can have a large effect on the amount of reinforcement due to hybrid incompatibility. Sex linkage of genes involved in postzygotic isolation generally increases the strength of reinforcement, but only if genes involved in prezygotic isolation are also sex linked. We use exact simulations to test the accuracy of the approximation and find that qualitative predictions made assuming weak selection can hold when selection is strong. Our simulations also show that incompatibilities that evolve in allopatry by drift or weak selection are likely to be lost by swamping during secondary contact, even when selection against hybrids is strong.     

A fine-scale genetic analysis of hybrid incompatibilities in Drosophila

The sterility and inviability of species hybrids is thought to evolve by the accumulation of genes that cause generally recessive, incompatible epistatic interactions between species. Most analyses of the loci involved in such hybrid incompatibilities have suffered from low genetic resolution. Here I present a fine-resolution genetic screen that allows systematic counting, mapping, and characterizing of a large number of hybrid incompatibility loci in a model genetic system. Using small autosomal deletions from D. melanogaster and a hybrid rescue mutation from D. simulans, I measured the viability of hybrid males that are simultaneously hemizygous for a small region of the D. simulans autosomal genome and hemizygous for the D. melanogaster X chromosome. These hybrid males are exposed to the full effects of any recessive-recessive epistatic incompatibilities present in these regions. A screen of approximately 70% of the D. simulans autosomal genome reveals 20 hybrid-lethal and 20 hybrid-semilethal regions that are incompatible with the D. melanogaster X. In further crosses, I confirm the epistatic nature of hybrid lethality by showing that all of the incompatibilities are rescued when the D. melanogaster X is replaced with a D. simulans X. Combined with information from previous studies, these results show that the number of recessive incompatibilities is approximately eightfold larger than the number of dominant ones. Finally, I estimate that a total of approximately 191 hybrid-lethal incompatibilities separate D. melanogaster and D. simulans, indicating extensive functional divergence between these species' genomes.     

Differential introgression across newt hybrid zones: Evidence from replicated transects

Genomic heterogeneity of divergence between hybridizing species may reflect heterogeneity of introgression, but also processes unrelated to hybridization. Heterogeneous introgression and its repeatability can be directly tested in natural hybrid zones by examining multiple transects. Here, we studied hybrid zones between the European newts Lissotriton montandoni and two lineages of Lissotriton vulgaris, with replicate transects within each zone. Over 1,000 nuclear genes located on a linkage map and mitochondrial DNA were investigated using geographical and genomic clines. Overall, the five transects were all similar, showing hallmarks of strong reproductive isolation: bimodal distribution of genotypes in central populations and narrow allele frequency clines. However, the extent of introgression differed between the zones, possibly as a consequence of their different ages, as suggested by the analysis of heterozygosity runs in diagnostic markers. In three transects genomic signatures of small‐scale (~2 km) zone movements were detected. We found limited overlap of cline outliers between transects, and only weak evidence of stronger differentiation of introgression between zones than between transects within zones. Introgression was heterogeneous across linkage groups, with patterns of heterogeneity similar between transects and zones. Predefined candidates for increased or reduced introgression exhibited only a subtle tendency in the expected direction, suggesting that interspecific differentiation is not a reliable indicator for the strength of introgression. These hierarchically sampled hybrid zones of apparently different ages show how introgression unfolds with time and offer an excellent opportunity to dissect the dynamics of hybridization and architecture of reproductive isolation at advanced stages of speciation.     

SNP genotyping for detecting the ‘rare allele phenomenon’ in hybrid zones

Hybrid zones are regions where genetically distinct populations meet, mate and produce offspring. In such zones, genetically less compatible gene combinations are usually generated, resulting in reduced fitness, and hybrid zones are often maintained because of continuous removal of unfit genotypes, balanced by gene flow into the zone from the parental populations (and are then referred to as ‘tension zones’). Tension zones often display unexpectedly high frequencies of gene variants that are rare outside the zone. Previous work has shown that this ‘rare allele phenomenon’ is not the result of intragenic recombination or increased mutation rates. Further understanding of the population genetics of the phenomenon requires an approach in which both the numbers of individuals and the numbers of loci is increased. Here, we report an approach using a combination of Illumina next‐generation sequencing and mass spectrophotometer genotyping to identify markers that may be used for genome‐wide investigations of the rare allele phenomenon. We test this approach on a hybrid zone in the land snail Albinaria hippolyti from Greece.     

Determining SNP allele frequencies in DNA pools

Single nucleotide polymorphisms (SNPs) are among the most common types of polymorphism used for genetic association studies. A method to allow the accurate quantitation of their allele frequencies from DNA pools would both increase throughput and decrease costs for large-scale genotyping. However, to date, most DNA pooling studies have concentrated on the use of microsatellite polymorphisms. In the case of SNPs that are restriction fragment length polymorphisms (RFLPs), studies have tended to use methods for the quantitation of allele frequency from pools that rely on densitometric evaluation of bands on an autoradiograph. Radiation-based methods have well-known drawbacks, and we present two alternative methods for the determination of SNP allele frequencies. For RFLPs, we used agarose gel electrophoresis of digested PCR products with ethidium bromide staining combined with densitometric analysis of gel images on a PC. For all types of SNP, we used allele-specific fluorescent probes in the Taqman assay to determine the relative frequencies of two different alleles. Both methods gave accurate and reproducible results, suggesting they are suitable for use in DNA pooling experiments.     

Estimating mutant microsatellite allele frequencies in somatic cells by small-pool PCR

Identifying microsatellite instability (MSI) by partitioning DNA into multiple small pools containing only single genome amounts of DNA results in trapping both progenitor and low-frequency mutant alleles into pools where they can be identified and counted following PCR. Statistical approaches determining both the frequencies and the significant differences between frequencies of these Poisson-distributed alleles are presented. Results indicate a level of sensitivity and quantification not possible by standard PCR methods. Using material from colon cancer patients with high levels of MSI in their tumors, we also present the molecular and robotic methods for carrying out such studies. Validation experiments indicated mutants detectable at frequencies >0.03 above background. Frequencies obtained in tumor tissue (>0.25) met the expectations of the approach. Significant levels of MSI were detected in the constitutive tissue of the patient carrying a germ-line mutation for mismatch repair, suggesting both mechanistic and clinical applications of the procedure.     

Genotypic proportions in hybrid zones

We study a hybrid zone between two populations of a diploid organism. The populations differ at one locus. Homozygotes have equal fitnesses and the heterozygote fitness is reduced by + ( is the birth rate deviation and is the death rate deviation). The populations extend along a one dimensional continuous habitat, and migration occurs by diffusion of individuals. The model is formulated as a set of simple continuous time demographic models without age structure for the three genotypes, and the system is transformed into three new variables, the total population size N, the gene frequency p, and the deviation from Hardy-Weinberg proportions F. The gene frequency in a steady state cline always follows a hyperbolic tangent closely. Analysis of the asymptotic behavior of the cline far from the hybrid zone suggests a qualitative prediction of the shape of N, p and F over the zone. For weak selection the shape is determined by a central steepness of , as observed by Bazykin in 1969, where D is the diffusion coefficient. For strong selection the cline is less steep than the Bazykin cline, and the form is dominated by the migration process. The steepness at the center of the cline is close to where b is the birth rate of homozygotes.Supported in part by grant 11-9639-1 from the Danish Natural Science Research Council and grant 5.18.2023 from the Danish Biotechnological Research and Development Programme     

Bimodal hybrid zones and speciation

Contact zones exemplify a series of stages in speciation. In unimodal hybrid zones intermediates predominate; in bimodal zones hybrids are rare and parental forms predominate; and finally, species might overlap, but never hybridize. Recent studies show bimodality to be associated strongly with assortative mating or fertilization, and only weakly with overall levels of genetic divergence or intrinsic genomic incompatibility. Ecological divergence across most bimodal hybrid zones suggests that ecology contributes more to speciation than genomic incompatibility. This continuum of stable contact zones provides empirical evidence for a route to speciation, which does not require allopatry.     

Genetic mapping in hybrid zones

Recent advances in genetic mapping methodologies make it feasible to localize quantitative trait loci (QTL) that contribute to adaptation and speciation. However, it has not been possible to employ these methods in many wild species because of difficulties associated with creating and propagating recombinant populations of sufficient size for QTL mapping. Natural hybrid zones contain recombinant individuals resulting from many generations of hybridization and thus offer a potential solution to these problems. For studies of speciation, hybrid zones offer the possibility of mapping QTL simultaneously with assessments of their effects on assortative mating, hybrid fitness, and interspecific gene flow. Here, we explore the problems and prospects associated with genetic map building and QTL analyses in natural hybrid zones by analyzing correlations among markers of known genomic location in four hybrid zones between the wild sunflower species Helianthus annuus and Helianthus petiolaris. Results indicate that mapping in hybrid zones presents many challenges. These include overlap in the strength of marker correlations between linked and unlinked markers, unevenness in marker frequencies along linkages, and heterogeneity in the relationship between marker distances and correlations. All make it difficult to accurately group and order markers or to estimate the distances between them. These problems can be ameliorated by sampling strategies that maximize the difference in linkage disequilibria between linked and unlinked markers and that minimize differences in frequencies among markers or QTL. In addition, studies that employ a previously determined molecular marker map for gene localization have a greater likelihood of success than those that rely on the hybrid zone data for both map construction and QTL analyses.     

Mutant allele frequencies in domestic cats of Taiwan

Preliminary data for Taiwanese cats generally agree with previous findings in far eastern populations, especially Vladivostok. However, surveys are still too few in number to achieve any detailed perspective for this vast region.     

Estimation of allele frequencies in ethnically heterogeneous populations

A method for reconstructing allele frequencies characteristic of an original ethnically homogeneous population before the start of migration processes is described. Information on both the ethnic group studied and offspring of interethnic marriages is used to estimate the allele frequencies. This makes it possible to increase the informativeness of the sample, which, in the case of ethnic heterogeneity, depends not only on allele frequencies and the total sample size, but also on the ethnic structure of the sample. The problem of estimating allele frequency in an ethnically heterogeneous sample has been solved analytically for diallelic loci. It has been demonstrated that, if offspring of interethnic marriages with the same degree of outbreeding is added to a sample of the ethnic group studied, the sample informativeness does not change. To utilize the information contained in the phenotypes of the offspring of interethnic marriages, representatives of the population from which migration occurs should be included into the sample. The size of the sample ensuring the preassigned accuracy of estimation is minimized at a certain ratio between the numbers of the offspring of interethnic marriages and the "immigrants." To analyze polyallelic loci, a software package has been developed that allows estimating allele frequencies, determining the errors of these estimates, and planning the sample ensuring the preassigned accuracy of estimation. The package is available free at http://mga.bionet.bionet.nsc.ru/PopMixed/PopMixed.html.     

Temperature stress increases hybrid incompatibilities in the parasitic wasp genus Nasonia

Hybrid incompatibilities, measured as mortality and sterility, are caused by the disruption of gene interactions. They are important post-zygotic isolation barriers to species hybridization, and much effort is put into the discovery of the genes underlying these incompatibilities. In hybridization studies of the haplodiploid parasitic wasp genus Nasonia, genic incompatibilities have been shown to affect mortality and sterility. The genomic regions associated with mortality have been found to depend on the cytotype of the hybrids and thus suggest cytonuclear incompatibilities. As environmental conditions can affect gene expression and gene interaction, we here investigate the effect of developmental temperature on sterility and mortality in Nasonia hybrids. Results show that extreme temperatures strongly affect both hybrid sterility (mainly spermatogenic failure) and mortality. Molecular mapping revealed that extreme temperatures increase transmission ratio distortion of parental alleles at incompatible loci, and thus, cryptic incompatible loci surface under temperature stress that remain undiscovered under standard temperatures. Our results underline the sensitivity of hybrid incompatibilities to environmental factors and the effects of unstable epistasis.     

Microsatellite allele frequencies in humans and chimpanzees, with implications for constraints on allele size

The distributions of allele sizes at eight simple-sequence repeat (SSR) or microsatellite loci in chimpanzees are found and compared with the distributions previously obtained from several human populations. At several loci, the differences in average allele size between chimpanzees and humans are sufficiently small that there might be a constraint on the evolution of average allele size. Furthermore, a model that allows for a bias in the mutation process shows that for some loci a weak bias can account for the observations. Several alleles at one of the loci (Mfd 59) were sequenced. Differences between alleles of different lengths were found to be more complex than previously assumed. An 8-base-pair deletion was present in the nonvariable region of the chimpanzee locus. This locus contains a previously unrecognized repeated region, which is imperfect in humans and perfect in chimpanzees. The apparently greater opportunity for mutation conferred by the two perfect repeat regions in chimpanzees is reflected in the higher variance in repeat number at Mfd 59 in chimpanzees than in humans. These data indicate that interspecific differences in allele length are not always attributable to simple changes in the number of repeats.      

Morphological distinctiveness of Ligularia tongolensis and L. cymbulifera is maintained between habitats despite bidirectional and asymmetrical introgression in multiple hybrid zones

Natural hybridization is a crucial evolutionary process and a long-standing topic of study in evolutionary biology. Hybrid zones, where two congeneric species interact, can provide insight into the process of natural hybridization, especially with respect to how taxon diversity is maintained. In this study, we used double digest restriction-site associated DNA sequencing technology (ddRAD-seq) to examine genetic structure and estimate introgression in four hybrid zones of Ligularia tongolensis and Ligularia cymbulifera. Our analysis demonstrated that parental species were highly differentiated, whereas pairwise F_ST between parents and their hybrids was low, indicating that sympatric sites can form hybrid zones. As most F₁ hybrid individuals were observed within these zones, our finding also implied the presence of substantial barriers to interbreeding. Furthermore, some individuals that possessed the typical morphology of the parental species belonged to the F₁ generation. Genomic clines analysis revealed that a large fraction of single nucleotide polymorphisms (SNPs) deviated from a model of neutral introgression in the four hybrid zones, and most SNPs exhibited selection favoring the L. cymbulifera genotype. Bidirectional but asymmetric introgression was revealed as evident in the four hybrid zones. Habitat differences between the four hybrid zones may affect isolation barriers between both species. Taken together, these findings suggest that where incomplete reproductive barriers allow natural hybridization, the introgression between species generates rich genetic recombination that contributes to the fast adaptation and diversification of the widespread Ligularia in the Hengduan Mountains Region (HMR).