Genome-wide dissection of hybrid sterility in Drosophila confirms a polygenic threshold architecture

On the impermanence of species: The collapse of genetic incompatibilities in hybridizing populations

Species pairs often become genetically incompatible during divergence, which is an important source of reproductive isolation. An idealized picture is often painted where incompatibility alleles accumulate and fix between diverging species. However, recent studies have shown both that incompatibilities can collapse with ongoing hybridization, and that incompatibility loci can be polymorphic within species. This paper suggests some general rules for the behavior of incompatibilities under hybridization. In particular, we argue that redundancy of genetic pathways can strongly affect the dynamics of intrinsic incompatibilities. Since fitness in genetically redundant systems is unaffected by introducing a few foreign alleles, higher redundancy decreases the stability of incompatibilities during hybridization, but also increases tolerance of incompatibility polymorphism within species. We use simulations and theories to show that this principle leads to two types of collapse: in redundant systems, exemplified by classical Dobzhansky-Muller incompatibilities, collapse is continuous and approaches a quasi-neutral polymorphism between broadly sympatric species, often as a result of isolation-by-distance. In nonredundant systems, exemplified by co-evolution among genetic elements, incompatibilities are often stable, but can collapse abruptly with spatial traveling waves. As both types are common, the proposed principle may be useful in understanding the abundance of genetic incompatibilities in natural populations.

Divergence of X-linked trans regulatory proteins and the misexpression of gene targets in sterile Drosophila pseudoobscura hybrids

BACKGROUND

The genetic basis of hybrid incompatibilities is characterized by pervasive cases of gene interactions. Sex chromosomes play a major role in speciation and X-linked hybrid male sterility (HMS) genes have been identified. Interestingly, some of these genes code for proteins with DNA binding domains, suggesting a capability to act as trans-regulatory elements and disturb the expression of a large number of gene targets. To understand how interactions between trans- and cis-regulatory elements contribute to speciation, we aimed to map putative X-linked trans-regulatory elements and to identify gene targets with disrupted gene expression in sterile hybrids between the subspecies Drosophila pseudoobscura pseudoobscura and D. p. bogotana.

RESULTS

We find six putative trans-regulatory proteins within previously mapped X chromosome HMS loci with sequence changes that differentiate the two subspecies. Among them, the previously characterized HMS gene Overdrive (Ovd) had the largest number of amino acid changes between subspecies, with some substitutions localized within the protein's DNA binding domain. Using an introgression approach, we detected transcriptional responses associated with a sterility/fertility Ovd allele swap. We found a network of 52 targets of Ovd and identified cis-regulatory effects among target genes with disrupted expression in sterile hybrids. However, a combined analysis of polymorphism and divergence in non-coding sequences immediately upstream of target genes found no evidence of changes in candidate regulatory proximal cis-elements. Finally, peptidases were over-represented among target genes.

CONCLUSIONS

We provide evidence of divergence between subspecies within the DNA binding domain of the HMS protein Ovd and identify trans effects on the expression of 52 gene targets. Our results identify a network of trans-cis interactions with possible effects on HMS. This network provides molecular evidence of gene × gene incompatibilities as contributors to hybrid dysfunction.

Hybrid Incompatibilities and Transgressive Gene Expression Between Two Closely Related Subspecies of Drosophila

Genome-wide assays of expression between species and their hybrids have identified genes that become either over- or underexpressed relative to the parental species (i.e., transgressive). Transgressive expression in hybrids is of interest because it highlights possible changes in gene regulation linked to hybrid dysfunction. Previous studies in Drosophila that used long-diverged species pairs with complete or nearly complete isolation (i.e., full sterility and partial inviability of hybrids) and high-levels of genome misregulation have found correlations between expression and coding sequence divergence. The work highlighted the possible effects of directional selection driving sequence divergence and transgressive expression. Whether the same is true for taxa at early stages of divergence that have only achieved partial isolation remains untested. Here, we reanalyze previously published genome expression data and available genome sequence reads from a pair of partially isolated subspecies of Drosophila to compare expression and sequence divergence. We find a significant correlation in rates of expression and sequence evolution, but no support for directional selection driving transgressive expression in hybrids. We find that most transgressive genes in hybrids show no differential expression between parental subspecies and used SNP data to explore the role of stabilizing selection through compensatory mutations. We also examine possible misregulation through cascade effects that could be driven by interacting gene networks or co-option of off-target cis-regulatory elements.

Multi-locus interactions and the build-up of reproductive isolation

All genes interact with other genes, and their additive effects and epistatic interactions affect an organism's phenotype and fitness. Recent theoretical and empirical work has advanced our understanding of the role of multi-locus interactions in speciation. However, relating different models to one another and to empirical observations is challenging. This review focuses on multi-locus interactions that lead to reproductive isolation (RI) through reduced hybrid fitness. We first review theoretical approaches and show how recent work incorporating a mechanistic understanding of multi-locus interactions recapitulates earlier models, but also makes novel predictions concerning the build-up of RI. These include high variance in the build-up rate of RI among taxa, the emergence of strong incompatibilities producing localized barriers to introgression, and an effect of population size on the build-up of RI. We then review recent experimental approaches to detect multi-locus interactions underlying RI using genomic data. We argue that future studies would benefit from overlapping methods like ancestry disequilibrium scans, genome scans of differentiation and analyses of hybrid gene expression. Finally, we highlight a need for further overlap between theoretical and empirical work, and approaches that predict what kind of patterns multi-locus interactions resulting in incompatibilities will leave in genome-wide polymorphism data. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

The genetic basis of hybrid male sterility in sympatric Primulina species

BACKGROUND

Sympatric sister species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain species boundaries. The persistence of morphologically and genetically distinct populations in sympatry can only occur if some degree of reproductive isolation exists. A pair of sympatric sister species of Primulina (P. depressa and P. danxiaensis) was used to explore the genetic architecture of hybrid male sterility.

RESULTS

We mapped one major- and seven minor-effect quantitative trait loci (QTLs) that underlie pollen fertility rate (PFR). These loci jointly explained 55.4% of the phenotypic variation in the F2 population. A Bateson-Dobzhansky-Muller (BDM) model involving three loci was observed in this system. We found genotypic correlations between hybrid male sterility and flower morphology, consistent with the weak but significant phenotypic correlations between PFR and floral traits.

CONCLUSIONS

Hybrid male sterility in Primulina is controlled by a polygenic genetic basis with a complex pattern. The genetic incompatibility involves a three-locus BDM model. Hybrid male sterility is genetically correlated with floral morphology and divergence hitchhiking may occur between them.

Altered expression of cell adhesion genes and hybrid male sterility between subspecies ofDrosophila pseudoobscura

Drosophila pseudoobscura pseudoobscura and Drosophila pseudoobscura bogotana are two closely related subspecies with incomplete reproductive isolation. A genome-wide comparison of expression in hybrids relative to parental subspecies has been previously used to identify genes with significant changes in expression uniquely associated with the sterile condition. The misexpression (i.e., gene expression beyond levels found in parentals) of such genes could be directly linked to the onset of sterility or could alternatively be caused by incompatibilities in a hybrid genome without a direct link to sterility. Cell adhesion was previously found to be one of the largest gene ontologies with changes in expression linked to sterility. Here we used gene expression assays in fertile backcross male progeny, along with introgression progeny in which we swap a major hybrid male sterility (HMS) allele, to generate fertile and sterile males genotypically similar to F1sterile hybrids. We identify a cell adhesion gene (GA10921) whose change in expression is directly linked to sterility and modulated by a previously characterized HMS protein. GA10921 adds to our rather limited knowledge of changes in gene expression associated with HMS, and to the identification of gene interacting partners linked to HMS.

Coadapted genomes and selection on hybrids: Fisher's geometric model explains a variety of empirical patterns

Natural selection plays a variety of roles in hybridization, speciation, and admixture. Most research has focused on two extreme cases: crosses between closely related inbred lines, where hybrids are fitter than their parents, or crosses between effectively isolated species, where hybrids suffer severe breakdown. But many natural populations must fall into intermediate regimes, with multiple types of gene interaction, and these are more difficult to study. Here, we develop a simple fitness landscape model, and show that it naturally interpolates between previous modeling approaches, which were designed for the extreme cases, and invoke either mildly deleterious recessives, or discrete hybrid incompatibilities. Our model yields several new predictions, which we test with genomic data from Mytilus mussels, and published data from plants (Zea, Populus, and Senecio) and animals (Mus, Teleogryllus, and Drosophila). The predictions are generally supported, and the model explains a number of surprising empirical patterns. Our approach enables novel and complementary uses of genome-wide datasets, which do not depend on identifying outlier loci, or "speciation genes" with anomalous effects. Given its simplicity and flexibility, and its predictive successes with a wide range of data, the approach should be readily extendable to other outstanding questions in the study of hybridization.

Drosophila Females Undergo Genome Expansion after Interspecific Hybridization

Genome size (or C-value) can present a wide range of values among eukaryotes. This variation has been attributed to differences in the amplification and deletion of different noncoding repetitive sequences, particularly transposable elements (TEs). TEs can be activated under different stress conditions such as interspecific hybridization events, as described for several species of animals and plants. These massive transposition episodes can lead to considerable genome expansions that could ultimately be involved in hybrid speciation processes. Here, we describe the effects of hybridization and introgression on genome size of Drosophila hybrids. We measured the genome size of two close Drosophila species, Drosophila buzzatii and Drosophila koepferae, their F1 offspring and the offspring from three generations of backcrossed hybrids; where mobilization of up to 28 different TEs was previously detected. We show that hybrid females indeed present a genome expansion, especially in the first backcross, which could likely be explained by transposition events. Hybrid males, which exhibit more variable C-values among individuals of the same generation, do not present an increased genome size. Thus, we demonstrate that the impact of hybridization on genome size can be detected through flow cytometry and is sex-dependent.

Expression of the Retrotransposon Helena Reveals a Complex Pattern of TE Deregulation in Drosophila Hybrids

Transposable elements (TEs), repeated mobile sequences, are ubiquitous in the eukaryotic kingdom. Their mobilizing capacity confers on them a high mutagenic potential, which must be strongly regulated to guarantee genome stability. In the Drosophila germline, a small RNA-mediated silencing system, the piRNA (Piwi-interacting RNA) pathway, is the main responsible TE regulating mechanism, but some stressful conditions can destabilize it. For instance, during interspecific hybridization, genomic stress caused by the shock of two different genomes can lead, in both animals and plants, to higher transposition rates. A recent study in D. buzatii-D. koepferae hybrids detected mobilization of 28 TEs, yet little is known about the molecular mechanisms explaining this transposition release. We have characterized one of the mobilized TEs, the retrotransposon Helena, and used quantitative expression to assess whether its high transposition rates in hybrids are preceded by increased expression. We have also localized Helena expression in the gonads to see if cellular expression patterns have changed in the hybrids. To give more insight into changes in TE regulation in hybrids, we analysed Helena-specific piRNA populations of hybrids and parental species. Helena expression is not globally altered in somatic tissues, but male and female gonads have different patterns of deregulation. In testes, Helena is repressed in F1, increasing then its expression up to parental values. This is linked with a mislocation of Helena transcripts along with an increase of their specific piRNA levels. Ovaries have additive levels of Helena expression, but the ping-pong cycle efficiency seems to be reduced in F1 hybrids. This could be at the origin of new Helena insertions in hybrids, which would be transmitted to F1 hybrid female progeny.

Genetic Linkage Map Construction and QTL Analysis of Two Interspecific Reproductive Isolation Traits in Sponge Gourd

The hybrids between Luffa acutangula (L.) Roxb. and L.cylindrica (L.) Roem. have strong heterosis effects. However, some reproductive isolation traits hindered their normal hybridization and fructification, which was mainly caused by the flowering time and hybrid pollen sterility. In order to study the genetic basis of two interspecific reproductive isolation traits, we constructed a genetic linkage map using an F2 population derived from a cross between S1174 [L. acutangula (L.) Roxb.] and 93075 [L. cylindrica (L.) Roem.]. The map spans 1436.12 CentiMorgans (cM), with an average of 8.11 cM among markers, and consists of 177 EST-SSR markers distributed in 14 linkage groups (LG) with an average of 102.58 cM per LG. Meanwhile, we conducted colinearity analysis between the sequences of EST-SSR markers and the genomic sequences of cucumber, melon and watermelon. On the basis of genetic linkage map, we conducted QTL mapping of two reproductive isolation traits in sponge gourd, which were the flowering time and hybrid male sterility. Two putative QTLs associated with flowering time (FT) were both detected on LG 1. The accumulated contribution of these two QTLs explained 38.07% of the total phenotypic variance (PV), and each QTL explained 15.36 and 22.71% of the PV respectively. Four QTLs for pollen fertility (PF) were identified on LG 1 (qPF1.1 and qPF1.2), LG 3 (qPF3) and LG 7 (qPF7), respectively. The percentage of PF explained by these QTLs varied from 2.91 to 16.79%, and all together the four QTLs accounted for 39.98% of the total PV. Our newly developed EST-SSR markers and linkage map are very useful for gene mapping, comparative genomics and molecular marker-assisted breeding. These QTLs for interspecific reproductive isolation will also contribute to the cloning of genes relating to interspecific reproductive isolation and the utilization of interspecific heterosis in sponge gourd in further studies.

Changes of Osvaldo expression patterns in germline of male hybrids between the species Drosophila buzzatii and Drosophila koepferae

Hybridization between different genomes is a source of genomic instability, sometimes associated with transposable element (TE) mobilization. Previous work showed that hybridization between the species Drosophila buzzatii and Drosophila koepferae induced mobilization of different (TEs), the Osvaldo retrotransposon being the most unstable. However, we ignore the mechanisms involved in this transposition release in interspecific hybrids. In order to disentangle the mechanisms involved in this process, we performed Osvaldo expression studies in somatic and germinal tissues from hybrids and parental species. There was a trend towards increased Osvaldo expression in the somatic tissues of hybrid females and males, which was always significant in males compared to the parental species D. buzzatii but, not in females compared to maternal species D. koepferae. There were massive changes of Osvaldo expression in the testes, which varied depending on the hybrid generation and family. Moreover, Osvaldo hybridization signals, restricted to the apical and primary spermatocyte regions in parents, occupied broader region in the hybrids. In ovaries, there were no significant differences in Osvaldo expression rates between hybrids and the maternal species D. koepferae. The transcript location was restricted to ovarian nurse cells in both parents and hybrids, undetectable in some hybrids. This research highlights first, the existence of putative complex deregulation mechanisms different between sexes and cell types and second, disruption of Osvaldo activity particularly evident in testes from sterile hybrid males. Deeper studies of the total transcriptome in hybrids and parental species are necessary to gain a better knowledge of the TE deregulation pathways in hybrids.

Genome-wide mapping in a house mouse hybrid zone reveals hybrid sterility loci and Dobzhansky-Muller interactions

Mapping hybrid defects in contact zones between incipient species can identify genomic regions contributing to reproductive isolation and reveal genetic mechanisms of speciation. The house mouse features a rare combination of sophisticated genetic tools and natural hybrid zones between subspecies. Male hybrids often show reduced fertility, a common reproductive barrier between incipient species. Laboratory crosses have identified sterility loci, but each encompasses hundreds of genes. We map genetic determinants of testis weight and testis gene expression using offspring of mice captured in a hybrid zone between M. musculus musculus and M. m. domesticus. Many generations of admixture enables high-resolution mapping of loci contributing to these sterility-related phenotypes. We identify complex interactions among sterility loci, suggesting multiple, non-independent genetic incompatibilities contribute to barriers to gene flow in the hybrid zone.

Interspecific hybridization as a genomic stressor inducing mobilization of transposable elements in Drosophila

Transposable elements (TEs) are DNA sequences able to be mobilized in host genomes. They are currently recognized as the major mutation inducers because of their insertion in the target, their effect on neighboring regions, or their ectopic recombination. A large number of factors including chemical and physical factors as well as intraspecific crosses have traditionally been identified as inducers of transposition. Besides environmental factors, interspecific crosses have also been proposed as promoters of transposition of particular TEs in plants and different animals. Our previous published work includes a genome-wide survey with the set of genomic TEs and shows that interspecific hybridization between the species Drosophila buzzatii and Drosophila koepferae induces genomic instability by transposition bursts. A high percentage of this instability corresponds to TEs belonging to classes I and II. The detailed study of three TEs (Osvaldo, Helena, and Galileo), representative of the different TE families, shows an increase of transposition in hybrids compared with parental species, that varies depending on the element. This study suggests ample variation in TE regulation mechanisms and the question is why this variation occurs. Interspecific hybridization is a genomic stressor that disrupts the stability of TEs probably contributing to a relaxation of the mechanisms controlling TEs in the Drosophila genome. In this commentary paper we will discuss these results and the molecular mechanisms that could explain these increases of transposition rates observed in interspecific Drosophila hybrids.