Genetic architecture of hybrid male sterility in Drosophila: analysis of intraspecies variation for interspecies isolation

Identification and genetic analysis of a pervasive 'needle-eye' sperm phenotype in Drosophila sterile hybrid males

Interspecies hybrid sterility has been extensively studied, especially in the genus Drosophila. Hybrid sterility is more often found in the heterogametic (XY or ZW) sex, a trend called Haldane's rule. Although this phenomenon is pervasive, identification of a common genetic mechanism remains elusive, with modest support found for a range of potential theories. Here, we identify a single precise morphological phenotype, which we call 'needle-eye sperm', that is associated with hybrid sterility in three separate species pairs that span the Drosophila genus. The nature of the phenotype indicates a common point of meiotic failure in sterile hybrid males. We used 10 generations of backcross selection paired with whole-genome pooled sequencing to genetically map the regions underlying the needle-eye (NE) sperm phenotype. Surprisingly, the sterility phenotype was present in ~50% of males even after 10 generations of backcrossing, and only a single region of the X chromosome was associated with sterility in one direction of backcross. Owing to the common phenotype among sterile male hybrids, and the strong effect of individual loci, further exploration of these findings may identify a universal mechanism for the evolution of hybrid sterility.

Hybrid breakdown in male reproduction between recently diverged Drosophila melanogaster populations has a complex and variable genetic architecture

Secondary contact between formerly isolated populations may result in hybrid breakdown, in which untested allelic combinations in hybrids are maladaptive and limit genetic exchange. Studying early-stage reproductive isolation may yield key insights into the genetic architectures and evolutionary forces underlying the first steps toward speciation. Here, we leverage the recent worldwide expansion of Drosophila melanogaster to test for hybrid breakdown between populations that diverged within the last 13,000 years. We found clear evidence for hybrid breakdown in male reproduction, but not female reproduction or viability, supporting the prediction that hybrid breakdown affects the heterogametic sex first. The frequency of non-reproducing F2 males varied among different crosses involving the same southern African and European populations, as did the qualitative effect of cross direction, implying a genetically variable basis of hybrid breakdown and a role for uniparentally inherited factors. The levels of breakdown observed in F2 males were not recapitulated in backcrossed individuals, consistent with the existence of incompatibilities with at least three partners. Thus, some of the very first steps toward reproductive isolation could involve incompatibilities with complex and variable genetic architectures. Collectively, our findings emphasize this system's potential for future studies on the genetic and organismal basis of early-stage reproductive isolation.

Transposable Element Expression and Regulation Profile in Gonads of Interspecific Hybrids of Drosophila arizonae and Drosophila mojavensis wrigleyi

Interspecific hybridization may lead to sterility and/or inviability through differential expression of genes and transposable elements (TEs). In Drosophila, studies have reported massive TE mobilization in hybrids from interspecific crosses of species presenting high divergence times. However, few studies have examined the consequences of TE mobilization upon hybridization in recently diverged species, such as Drosophila arizonae and D. mojavensis. We have sequenced transcriptomes of D. arizonae and the subspecies D. m. wrigleyi and their reciprocal hybrids, as well as piRNAs, to analyze the impact of genomic stress on TE regulation. Our results revealed that the differential expression in both gonadal tissues of parental species was similar. Globally, ovaries and testes showed few deregulated TEs compared with both parental lines. Analyses of small RNA data showed that in ovaries, the TE upregulation is likely due to divergence of copies inherited from parental genomes and lack of piRNAs mapping to them. Nevertheless, in testes, the divergent expression of genes associated with chromatin state and piRNA pathway potentially indicates that TE differential expression is related to the divergence of regulatory genes that play a role in modulating transcriptional and post-transcriptional mechanisms.

Comparative transcriptomics between Drosophila mojavensis and D. arizonae reveals transgressive gene expression and underexpression of spermatogenesis-related genes in hybrid testes

Interspecific hybridization is a stressful condition that can lead to sterility and/or inviability through improper gene regulation in Drosophila species with a high divergence time. However, the extent of these abnormalities in hybrids of recently diverging species is not well known. Some studies have shown that in Drosophila, the mechanisms of postzygotic isolation may evolve more rapidly in males than in females and that the degree of viability and sterility is associated with the genetic distance between species. Here, we used transcriptomic comparisons between two Drosophila mojavensis subspecies and D. arizonae (repleta group, Drosophila) and identified greater differential gene expression in testes than in ovaries. We tested the hypothesis that the severity of the interspecies hybrid phenotype is associated with the degree of gene misregulation. We showed limited gene misregulation in fertile females and an increase in the amount of misregulation in males with more severe sterile phenotypes (motile vs. amotile sperm). In addition, for these hybrids, we identified candidate genes that were mostly associated with spermatogenesis dysfunction.

A common suite of cellular abnormalities and spermatogenetic errors in sterile hybrid males in Drosophila

When two species interbreed, the resulting hybrid offspring are often sterile, with the heterogametic (e.g. XY) hybrid usually being more severely affected. The prevailing theory for this pattern of sterility evokes divergent changes in separate lineages having maladaptive interactions when placed together in a hybrid individual, with recessive factors on the sex chromosome interacting with dominant factors on the autosomes. The effect of these interactions on gametogenesis should not be uniform across species pairs unless genetic divergence follows the same paths in different lineages or if a specific stage of gametogenesis is more susceptible to detrimental genetic interactions. Here, we perform a detailed cellular characterization of hybrid male sterility across three recently diverged species pairs of Drosophila. Across all three pairs, sterile hybrid sperm are alive but exhibit rapid nuclear de-condensation with age, with active, but non-differentiated, mitochondria. Surprisingly, all three sets of interspecies hybrids produce half of the number of sperm per round of spermatogenesis, with each sperm cell containing two tails. We identify non-disjunction failures during meiosis I as the likely cause. Thus, errors during meiosis I may be a general phenomenon underlying Drosophila male sterility, indicating either a heightened sensitivity of this spermatogenic stage to failure, or a basis to sterility other than the prevailing model.

Identification of misexpressed genetic elements in hybrids between Drosophila-related species

Crosses between close species can lead to genomic disorders, often considered to be the cause of hybrid incompatibility, one of the initial steps in the speciation process. How these incompatibilities are established and what are their causes remain unclear. To understand the initiation of hybrid incompatibility, we performed reciprocal crosses between two species of Drosophila (D. mojavensis and D. arizonae) that diverged less than 1 Mya. We performed a genome-wide transcriptomic analysis on ovaries from parental lines and on hybrids from reciprocal crosses. Using an innovative procedure of co-assembling transcriptomes, we show that parental lines differ in the expression of their genes and transposable elements. Reciprocal hybrids presented specific gene categories and few transposable element families misexpressed relative to the parental lines. Because TEs are mainly silenced by piwi-interacting RNAs (piRNAs), we hypothesize that in hybrids the deregulation of specific TE families is due to the absence of such small RNAs. Small RNA sequencing confirmed our hypothesis and we therefore propose that TEs can indeed be major players of genome differentiation and be implicated in the first steps of genomic incompatibilities through small RNA regulation.

Environmentally induced changes in correlated responses to selection reveal variable pleiotropy across a complex genetic network

Selection in novel environments can lead to a coordinated evolutionary response across a suite of characters. Environmental conditions can also potentially induce changes in the genetic architecture of complex traits, which in turn could alter the pattern of the multivariate response to selection. We describe a factorial selection experiment using the nematode Caenorhabditis remanei in which two different stress-related phenotypes (heat and oxidative stress resistance) were selected under three different environmental conditions. The pattern of covariation in the evolutionary response between phenotypes or across environments differed depending on the environment in which selection occurred, including asymmetrical responses to selection in some cases. These results indicate that variation in pleiotropy across the stress response network is highly sensitive to the external environment. Our findings highlight the complexity of the interaction between genes and environment that influences the ability of organisms to acclimate to novel environments. They also make clear the need to identify the underlying genetic basis of genetic correlations in order understand how patterns of pleiotropy are distributed across complex genetic networks.

X-y interactions underlie sperm head abnormality in hybrid male house mice

The genetic basis of hybrid male sterility in house mice is complex, highly polygenic, and strongly X linked. Previous work suggested that there might be interactions between the Mus musculus musculus X and the M. m. domesticus Y with a large negative effect on sperm head morphology in hybrid males with an F1 autosomal background. To test this, we introgressed the M. m. domesticus Y onto a M. m. musculus background and measured the change in sperm morphology, testis weight, and sperm count across early backcross generations and in 11th generation backcross males in which the opportunity for X-autosome incompatibilities is effectively eliminated. We found that abnormality in sperm morphology persists in M. m. domesticus Y introgression males, and that this phenotype is rescued by M. m. domesticus introgressions on the X chromosome. In contrast, the severe reductions in testis weight and sperm count that characterize F1 males were eliminated after one generation of backcrossing. These results indicate that X-Y incompatibilities contribute specifically to sperm morphology. In contrast, X-autosome incompatibilities contribute to low testis weight, low sperm count, and sperm morphology. Restoration of normal testis weight and sperm count in first generation backcross males suggests that a small number of complex incompatibilities between loci on the M. m. musculus X and the M. m. domesticus autosomes underlie F1 male sterility. Together, these results provide insight into the genetic architecture of F1 male sterility and help to explain genome-wide patterns of introgression across the house mouse hybrid zone.

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

To date, different studies about the genetic basis of hybrid male sterility (HMS), a postzygotic reproductive barrier thoroughly investigated using Drosophila species, have demonstrated that no single major gene can produce hybrid sterility without the cooperation of several genetic factors. Early work using hybrids between Drosophila koepferae (Dk) and Drosophila buzzatii (Db) was consistent with the idea that HMS requires the cooperation of several genetic factors, supporting a polygenic threshold (PT) model. Here we present a genome-wide mapping strategy to test the PT model, analyzing serially backcrossed fertile and sterile males in which the Dk genome was introgressed into the Db background. We identified 32 Dk-specific markers significantly associated with hybrid sterility. Our results demonstrate 1) a strong correlation between the number of segregated sterility markers and males' degree of sterility, 2) the exchangeability among markers, 3) their tendency to cluster into low-recombining chromosomal regions, and 4) the requirement for a minimum number (threshold) of markers to elicit sterility. Although our findings do not contradict a role for occasional major hybrid-sterility genes, they conform more to the view that HMS primarily evolves by the cumulative action of many interacting genes of minor effect in a complex PT architecture.

Insect speciation rules: unifying concepts in speciation research

The study of speciation is concerned with understanding the connection between causes of divergent evolution and the origin and maintenance of barriers to gene exchange between incipient species. Although the field has historically focused either on examples of recent divergence and its causes or on the genetic basis of reproductive isolation between already divergent species, current efforts seek to unify these two approaches. Here we integrate these perspectives through a discussion of recent progress in several insect speciation model systems. We focus on the evolution of speciation phenotypes in each system (i.e., those phenotypes causally involved in reducing gene flow between incipient species), drawing an explicit connection between cause and effect (process and pattern). We emphasize emerging insights into the genomic architecture of speciation as well as timely areas for future research.

A novel approach identifying hybrid sterility QTL on the autosomes of Drosophila simulans and D. mauritiana

When species interbreed, the hybrid offspring that are produced are often sterile. If only one hybrid sex is sterile, it is almost always the heterogametic (XY or ZW) sex. Taking this trend into account, the predominant model used to explain the genetic basis of F₁ sterility involves a deleterious interaction between recessive sex-linked loci from one species and dominant autosomal loci from the other species. This model is difficult to evaluate, however, as only a handful of loci influencing interspecies hybrid sterility have been identified, and their autosomal genetic interactors have remained elusive. One hindrance to their identification has been the overwhelming effect of the sex chromosome in mapping studies, which could ‘mask’ the ability to accurately map autosomal factors. Here, we use a novel approach employing attached-X chromosomes to create reciprocal backcross interspecies hybrid males that have a non-recombinant sex chromosome and recombinant autosomes. The heritable variation in phenotype is thus solely caused by differences in the autosomes, thereby allowing us to accurately identify the number and location of autosomal sterility loci. In one direction of backcross, all males were sterile, indicating that sterility could be entirely induced by the sex chromosome complement in these males. In the other direction, we identified nine quantitative trait loci that account for a surprisingly large amount (56%) of the autosome-induced phenotypic variance in sterility, with a large contribution of autosome-autosome epistatic interactions. These loci are capable of acting dominantly, and thus could contribute to F₁ hybrid sterility.

Mapping of within-species segregation distortion in Drosophila persimilis and hybrid sterility between D. persimilis and D. pseudoobscura

In contrast to the prevailing dogma in the 1990s, recent studies have suggested that an evolutionary history of segregation distortion within species may contribute to sterility in species hybrids. However, this recent work identified segregation distortion exclusively in species hybrids that may never have had an evolutionary history of segregation distortion in either parent species. We expand on previous work using a strain of Drosophila persimilis exhibiting segregation distortion within species to generate QTL maps for segregation distortion and hybrid sterility in crosses between D. persimilis and D. pseudoobscura. The maps localize regions along the XR contributing to both phenotypes, and they indicate one region of overlap between the two maps. This overlap could provide preliminary evidence for an association between segregation distortion within species and hybrid sterility, but the localizations are currently too broad to have confidence in this conclusion. This work is a first step towards possibly supporting a genetic conflict model of speciation in this system.

Dissecting the genetic architecture of F1 hybrid sterility in house mice

Hybrid sterility as a postzygotic reproductive isolation mechanism has been studied for over 80 years, yet the first identifications of hybrid sterility genes in Drosophila and mouse are quite recent. To study the genetic architecture of F(1) hybrid sterility between young subspecies of house mouse Mus m. domesticus and M. m. musculus, we conducted QTL analysis of a backcross between inbred strains representing these two subspecies and probed the role of individual chromosomes in hybrid sterility using the intersubspecific chromosome substitution strains. We provide direct evidence that the asymmetry in male infertility between reciprocal crosses is conferred by the middle region of M. m. musculus Chr X, thus excluding other potential candidates such as Y, imprinted genes, and mitochondrial DNA. QTL analysis identified strong hybrid sterility loci on Chr 17 and Chr X and predicted a set of interchangeable autosomal loci, a subset of which is sufficient to activate the Dobzhansky-Muller incompatibility of the strong loci. Overall, our results indicate the oligogenic nature of F(1) hybrid sterility, which should be amenable to reconstruction by proper combination of chromosome substitution strains. Such a prefabricated model system should help to uncover the gene networks and molecular mechanisms underlying hybrid sterility.

Incompatibility between X chromosome factor and pericentric heterochromatic region causes lethality in hybrids between Drosophila melanogaster and its sibling species

The Dobzhansky-Muller model posits that postzygotic reproductive isolation results from the evolution of incompatible epistatic interactions between species: alleles that function in the genetic background of one species can cause sterility or lethality in the genetic background of another species. Progress in identifying and characterizing factors involved in postzygotic isolation in Drosophila has remained slow, mainly because Drosophila melanogaster, with all of its genetic tools, forms dead or sterile hybrids when crossed to its sister species, D. simulans, D. sechellia, and D. mauritiana. To circumvent this problem, we used chromosome deletions and duplications from D. melanogaster to map two hybrid incompatibility loci in F(1) hybrids with its sister species. We mapped a recessive factor to the pericentromeric heterochromatin of the X chromosome in D. simulans and D. mauritiana, which we call heterochromatin hybrid lethal (hhl), which causes lethality in F(1) hybrid females with D. melanogaster. As F(1) hybrid males hemizygous for a D. mauritiana (or D. simulans) X chromosome are viable, the lethality of deficiency hybrid females implies that a dominant incompatible partner locus exists on the D. melanogaster X. Using small segments of the D. melanogaster X chromosome duplicated onto the Y chromosome, we mapped a dominant factor that causes hybrid lethality to a small 24-gene region of the D. melanogaster X. We provide evidence suggesting that it interacts with hhl(mau). The location of hhl is consistent with the emerging theme that hybrid incompatibilities in Drosophila involve heterochromatic regions and factors that interact with the heterochromatin.

Genetic variation for postzygotic reproductive isolation between Caenorhabditis briggsae and Caenorhabditis sp. 9

The process of speciation is key to the origins of biodiversity, and yet the Caenorhabditis nematode model system has contributed little to this topic. Genetic studies of speciation in the genus are now feasible, owing to crosses between the recently discovered Caenorhabditis sp. 9 and the well-known C. briggsae producing fertile F(1) hybrid females. We dissected patterns of postzygotic reproductive isolation between these species by crossing eight isogenic strains of C. briggsae reciprocally with six strains of C. sp. 9. We determined that overall patterns of reproductive isolation are robust across these genetic backgrounds. However, we also quantified significant heritable variation within each species for interspecific hybrid incompatibilities for total adult progeny, egg-to-adult viability, and the percentage of male progeny. This demonstrates that intraspecific variation for interspecific hybrid incompatibility occurs despite extensive, albeit incomplete, reproductive isolation. Therefore, this emerging general phenomenon of variable reproductive isolation is not restricted to highly interfertile, early-stage incipient species, but also applies to species in the latest stages of the speciation process. Furthermore, we confirm Haldane's rule and demonstrate strongly asymmetric parent-of-origin effects (Darwin's corollary) that consistently manifest more extremely when hermaphroditic C. briggsae serves as maternal parent. These findings highlight Caenorhabditis as an emerging system for understanding the genetics of general patterns of reproductive isolation.

An unusual barrier to gene flow: perpetually immature larvae from inter-population crosses in the flour beetle, Tribolium castaneum

We genetically characterize an unusual hybrid incompatibility phenotype manifest in F(1) offspring of crosses between two populations of Tribolium castaneum. Hybrid larvae cease development at the third larval instar, persisting as 'perpetually immature larvae' thereafter. Although unable to produce viable adult hybrid offspring with one another, each population produces abundant, fertile hybrids with other populations, indicating a recent origin of the incompatibility and facilitating genetic studies. We mapped the paternal component of the hybrid phenotype to a single region, which exhibits two characteristics common to hybrid incompatibility: marker transmission ratio distortion within crosses and elevated genetic divergence between populations. The incompatible variation and an elevation in between-population genetic divergence is associated with a region containing the T. castaneum ecdysone receptor homologue, a major regulatory switch, controlling larval moults, pupation and metamorphosis. This contributes to understanding the genetics of speciation in the Coleoptera, one of the most speciose of all arthropod taxa.

Genetics of hybrid male sterility among strains and species in the Drosophila pseudoobscura species group

Taxa in the early stages of speciation may bear intraspecific allelic variation at loci conferring barrier traits in hybrids such as hybrid sterility. Additionally, hybridization may spread alleles that confer barrier traits to other taxa. Historically, few studies examine within- and between-species variation at loci conferring reproductive isolation. Here, we test for allelic variation within Drosophila persimilis and within the Bogota subspecies of D. pseudoobscura at regions previously shown to contribute to hybrid male sterility. We also test whether D. persimilis and the USA subspecies of D. pseudoobscura share an allele conferring hybrid sterility in a D. pseudoobscura bogotana genetic background. All loci conferred similar hybrid sterility effects across all strains studied, although we detected some statistically significant quantitative effect variation among D. persimilis alleles of some hybrid incompatibility QTLs. We also detected allelism between D. persimilis and D. pseudoobscura USA at a second chromosome hybrid sterility QTL. We hypothesize that either the QTL is ancestral in D. persimilis and D. pseudoobscura USA and lost in D. pseudoobscura bogotana, or gene flow transferred the QTL from D. persimilis to D. pseudoobscura USA. We discuss our findings in the context of population features that may contribute to variation in hybrid incompatibilities.

Nuclear-cytoplasmic interactions reduce male fertility in hybrids of Arabidopsis lyrata subspecies

We examined the level of postzygotic reproductive isolation in F(1) and F(2) hybrids of reciprocal crosses between the Arabidopsis lyrata subspecies lyrata (North American) and petraea (European). Our main results are: first, the percentage of fertile pollen was significantly reduced in the F(1) and F(2) compared to the parental populations. Second, mean pollen fertility differed markedly between reciprocal crosses: 84% in the F(2) with ssp. lyrata cytoplasm and 61% in the F(2) with ssp. petraea cytoplasm. Third, 17% of the F(2) with ssp. petraea cytoplasm showed male sterility (produced less than 30 pollen grains in our subsample). The hybrids were female fertile. We used QTL mapping to find the genomic regions that determine pollen fertility and that restore cytoplasmic male sterility (CMS). In the F(2) with ssp. lyrata cytoplasm, an epistatic pair of QTLs was detected. In the reciprocal F(2) progeny, four QTLs demonstrated within-population polymorphism for hybrid male sterility. In addition, in the F(2) with ssp. petraea cytoplasm, there was a strong male fertility restorer locus on chromosome 2 where a cluster of CMS restorer gene-related PPR genes have been found in A. lyrata. Our results underline the importance of cytonuclear interactions in understanding genetics of the early stages of speciation.

Drosophilagenomes and the development of affordable molecular markers for species genotyping

The recent completion of genome sequencing of 12 species of Drosophila has provided a powerful resource for hypothesis testing, as well as the development of technical tools. Here we take advantage of genome sequence data from two closely related species of Drosophila, Drosophila simulans and Drosophila sechellia, to quickly identify candidate molecular markers for genotyping based on expected insertion or deletion (indel) differences between species. Out of 64 candidate molecular markers selected along the second and third chromosome of Drosophila, 51 molecular markers were validated using PCR and gel electrophoresis. We found that the 20% error rate was due to sequencing errors in the genome data, although we cannot rule out possible indel polymorphisms. The approach has the advantage of being affordable and quick, as it only requires the use of bioinformatics tools for predictions and a PCR and agarose gel based assay for validation. Moreover, the approach could be easily extended to a wide variety of taxa with the only limitation being the availability of complete or partial genome sequence data.

Cryptic postzygotic isolation in an eruptive species of bark beetle (Dendroctonus ponderosae)

Studies of postzygotic isolation often involve well-differentiated taxa that show a consistent level of incompatibility, thereby limiting our understanding of the initial stages and development of reproductive barriers. Dendroctonus ponderosae provides an informative system because recent evidence suggests that distant populations produce hybrids with reproductive incompatibilities. Dendroctonus ponderosae shows an isolation-by-distance gene flow pattern allowing us to characterize the evolution of postzygotic isolation (e.g., hybrid inviability, hybrid sterility) by crossing populations along a continuum of geographic/genetic divergence. We found little evidence of hybrid inviability among these crosses. However, crosses between geographically distant populations produced sterile males (consistent with Haldane's rule). This effect was not consistent with the fixation of mutations in an isolation-by-distance pattern, but instead is spatially localized. These reproductive barriers are uncorrelated with a reduction in gene flow suggesting their recent development. Crosses between geographically proximal populations bounding the transition from compatibility to hybrid male sterility showed evidence of unidirectional reduction in hybrid male fecundity. Our study describes significant postzygotic isolation occurring across a narrow and molecularly cryptic geographic zone between the states of Oregon and Idaho. This study provides a view of the early stages of postzygotic isolation in a geographically widespread species.

Epistasis among Drosophila persimilis factors conferring hybrid male sterility with D. pseudoobscura bogotana

The Bateson-Dobzhansky-Muller model posits that hybrid incompatibilities result from genetic changes that accumulate during population divergence. Indeed, much effort in recent years has been devoted to identifying genes associated with hybrid incompatibilities, often with limited success, suggesting that hybrid sterility and inviability are frequently caused by complex interactions between multiple loci and not by single or a small number of gene pairs. Our previous study showed that the nature of epistasis between sterility-conferring QTL in the Drosophila persimilis-D. pseudoobscura bogotana species pair is highly specific. Here, we further dissect one of the three QTL underlying hybrid male sterility between these species and provide evidence for multiple factors within this QTL. This result indicates that the number of loci thought to contribute to hybrid dysfunction may have been underestimated, and we discuss how linkage and complex epistasis may be characteristic of the genetics of hybrid incompatibilities. We further pinpoint the location of one locus that confers hybrid male sterility when homozygous, dubbed "mule-like", to roughly 250 kilobases.