Cis-by-Trans regulatory divergence causes the asymmetric lethal effects of an ancestral hybrid incompatibility gene

Expression pattern of resynthesized allotetraploid Capsella is determined by hybridization, not whole-genome duplication

Polyploidization, the process leading to the increase in chromosome sets, is a major evolutionary transition in plants. Whole-genome duplication (WGD) within the same species gives rise to autopolyploids, whereas allopolyploids result from a compound process with two distinct components: WGD and interspecific hybridization. To dissect the instant effects of WGD and hybridization on gene expression and phenotype, we created a series of synthetic hybrid and polyploid Capsella plants, including diploid hybrids, autotetraploids of both parental species, and two kinds of resynthesized allotetraploids with different orders of WGD and hybridization. Hybridization played a major role in shaping the relative expression pattern of the neo-allopolyploids, whereas WGD had almost no immediate effect on relative gene expression pattern but, nonetheless, still affected phenotypes. No transposable element-mediated genomic shock scenario was observed in either neo-hybrids or neo-polyploids. Finally, WGD and hybridization interacted and the distorting effects of WGD were less strong in hybrids. Whole-genome duplication may even improve hybrid fertility. In summary, while the initial relative gene expression pattern in neo-allotetraploids was almost entirely determined by hybridization, WGD only had trivial effects on relative expression patterns, both processes interacted and had a strong impact on physical attributes and meiotic behaviors.

System drift and speciation

Even if a species' phenotype does not change over evolutionary time, the underlying mechanism may change, as distinct molecular pathways can realize identical phenotypes. Here we use linear system theory to explore the consequences of this idea, describing how a gene network underlying a conserved phenotype evolves, as the genetic drift of small changes to these molecular pathways causes a population to explore the set of mechanisms with identical phenotypes. To do this, we model an organism's internal state as a linear system of differential equations for which the environment provides input and the phenotype is the output, in which context there exists an exact characterization of the set of all mechanisms that give the same input-output relationship. This characterization implies that selectively neutral directions in genotype space should be common and that the evolutionary exploration of these distinct but equivalent mechanisms can lead to the reproductive incompatibility of independently evolving populations. This evolutionary exploration, or system drift, is expected to proceed at a rate proportional to the amount of intrapopulation genetic variation divided by the effective population size ( Ne$N_e$ ). At biologically reasonable parameter values this could lead to substantial interpopulation incompatibility, and thus speciation, on a time scale of Ne$N_e$ generations. This model also naturally predicts Haldane's rule, thus providing a concrete explanation of why heterogametic hybrids tend to be disrupted more often than homogametes during the early stages of speciation.

The Integrity of the HMR complex is necessary for centromeric binding and reproductive isolation in Drosophila

Postzygotic isolation by genomic conflict is a major cause for the formation of species. Despite its importance, the molecular mechanisms that result in the lethality of interspecies hybrids are still largely unclear. The genus Drosophila, which contains over 1600 different species, is one of the best characterized model systems to study these questions. We showed in the past that the expression levels of the two hybrid incompatibility factors Hmr and Lhr diverged in the two closely related Drosophila species, D. melanogaster and D. simulans, resulting in an increased level of both proteins in interspecies hybrids. The overexpression of the two proteins also leads to mitotic defects, a misregulation in the expression of transposable elements and decreased fertility in pure species. In this work, we describe a distinct six subunit protein complex containing HMR and LHR and analyse the effect of Hmr mutations on complex integrity and function. Our experiments suggest that HMR needs to bring together components of centromeric and pericentromeric chromatin to fulfil its physiological function and to cause hybrid male lethality.

A major cause of biological speciation is the sterility and/or lethality of hybrids. This hybrid lethality is thought to be the consequence of two incompatible genomes of the two different species. We used the fruit fly Drosophila melanogaster as a model system to isolate a defined protein complex, which mediates this hybrid lethality. Our data suggest that this complex containing six subunits has evolved in one Drosophila species (Drosophila melanogaster) to bring together components of centromeric and pericentromeric chromatin. We show that the integrity of the complex is necessary for its genomic binding patterns and its ability to maintain fertility in female Drosophila melanogaster flies. Hybrid males between Drosophila melanogaster and the very closely related species Drosophila simulans die because they contain elevated levels of this complex. These high levels result in mitotic defects and a misregulation in the expression of transposable elements in those hybrids. Our results show that mutations that interfere with the complex’s function in Drosophila melanogaster also fail to induce lethality in hybrids suggesting that its evolutionary acquired functions in one species induce lethality in interspecies hybrids.

Drosophila Satellite Repeats at the Intersection of Chromatin, Gene Regulation and Evolution

Satellite repeats make up a large fraction of the genomes of many higher eukaryotes. Until recently these sequences were viewed as molecular parasites with few functions. Drosophila melanogaster and related species have a wealth of diverse satellite repeats. Comparative studies of Drosophilids have been instrumental in understanding how these rapidly evolving sequences change and move. Remarkably, satellite repeats have been found to modulate gene expression and mediate genetic conflicts between chromosomes and between closely related fly species. This suggests that satellites play a key role in speciation. We have taken advantage of the depth of research on satellite repeats in flies to review the known functions of these sequences and consider their central role in evolution and gene expression.

Ecological divergence in sympatry causes gene misexpression in hybrids

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.

Adaptive evolution among cytoplasmic piRNA proteins leads to decreased genomic auto-immunity

In metazoan germlines, the piRNA pathway acts as a genomic immune system, employing small RNA-mediated silencing to defend host DNA from the harmful effects of transposable elements (TEs). Expression of genomic TEs is proposed to initiate self regulation by increasing the production of repressive piRNAs, thereby “adapting” piRNA-mediated control to the most active TE families. Surprisingly, however, piRNA pathway proteins, which execute piRNA biogenesis and enforce silencing of targeted sequences, evolve rapidly and adaptively in animals. If TE silencing is ensured through piRNA biogenesis, what necessitates changes in piRNA pathway proteins? Here we used interspecific complementation to test for functional differences between Drosophila melanogaster and D. simulans alleles of three adaptively evolving piRNA pathway proteins: Armitage, Aubergine and Spindle-E. In contrast to piRNA-mediated transcriptional regulators examined in previous studies, these three proteins have cytoplasmic functions in piRNA maturation and post-transcriptional silencing. Across all three proteins we observed interspecific divergence in the regulation of only a handful of TE families, which were more robustly silenced by the heterospecific piRNA pathway protein. This unexpected result suggests that unlike transcriptional regulators, positive selection has not acted on cytoplasmic piRNA effector proteins to enhance their function in TE repression. Rather, TEs may evolve to “escape” silencing by host proteins. We further discovered that D. simulans alleles of aub and armi exhibit enhanced off-target effects on host transcripts in a D. melanogaster background, as well as modest reductions in the efficiency of piRNA biogenesis, suggesting that promiscuous binding of D. simulans Aub and Armi proteins to host transcripts reduces their participation in piRNA production. Avoidance of genomic auto-immunity may therefore be a critical target of selection. Our observations suggest that piRNA effector proteins are subject to an evolutionary trade-off between defending the host genome from the harmful effect of TEs while also minimizing collateral damage to host genes.

Transposable elements are mobile fragments of selfish DNA that burden host genomes with deleterious mutations and incite genome instability. Host cells employ a specialized small-RNA mediated silencing pathway, the piRNA pathway, to act as a genomic immune system suppressing the mobilization of TEs. Changes in genomic TE content are met with rapid changes in the piRNA pool, thereby maintaining host control over transposition. However, piRNA pathway proteins—which enact piRNA biogenesis and silence target TEs—also evolve adaptively. To isolate forces that underlie this adaptive evolution, we examined functional divergence between two Drosophila species for three adaptively evolving piRNA pathway proteins. To our surprise, we found very few differences in TE regulation, suggesting that evolution has not generally acted to enhance control of TE parasites. Rather, we discovered interspecific differences in the regulation of host mRNAs for two proteins, which suggested that proteins evolve to avoid off-target silencing of host transcripts. We propose that the avoidance of such “genomic autoimmunity” is an important and underappreciated force driving the adaptive evolution of piRNA proteins.

A multi-layered structure of the interphase chromocenter revealed by proximity-based biotinylation

During interphase centromeres often coalesce into a small number of chromocenters, which can be visualized as distinct, DAPI dense nuclear domains. Intact chromocenters play a major role in maintaining genome stability as they stabilize the transcriptionally silent state of repetitive DNA while ensuring centromere function. Despite its biological importance, relatively little is known about the molecular composition of the chromocenter or the processes that mediate chromocenter formation and maintenance. To provide a deeper molecular insight into the composition of the chromocenter and to demonstrate the usefulness of proximity-based biotinylation as a tool to investigate those questions, we performed super resolution microscopy and proximity-based biotinylation experiments of three distinct proteins associated with the chromocenter in Drosophila. Our work revealed an intricate internal architecture of the chromocenter suggesting a complex multilayered structure of this intranuclear domain.

The subgenomes show asymmetric expression of alleles in hybrid lineages of Megalobrama amblycephala × Culter alburnus

Hybridization drives rapid speciation by shaping novel genotypic and phenotypic profiles. Genomic incompatibility and transcriptome shock have been observed in hybrids, although this is rarer in animals than in plants. Using the newly sequenced genomes of the blunt snout bream (Megalobrama amblycephala [BSB]) and the topmouth culter (Culter alburnus [TC]), we focused on the sequence variation and gene expression changes in the reciprocal intergeneric hybrid lineages (F₁-F₃) of BSB × TC. A genome-wide transcriptional analysis identified 145-974 expressed recombinant genes in the successive generations of hybrid fish, suggesting the rapid emergence of allelic variation following hybridization. Some gradual changes of gene expression with additive and dominance effects and various cis and trans regulations were observed from F₁ to F₃ in the two hybrid lineages. These asymmetric patterns of gene expression represent the alternative strategies for counteracting deleterious effects of the subgenomes and improving adaptability of novel hybrids. Furthermore, we identified positive selection and additive expression patterns in transforming growth factor, beta 1b (tgfb1b), which may account for the morphological variations of the pharyngeal jaw in the two hybrid lineages. Our current findings provide insights into the evolution of vertebrate genomes immediately following hybridization.

Support for the Dominance Theory in Drosophila Transcriptomes

Interactions among divergent elements of transcriptional networks from different species can lead to misexpression in hybrids through regulatory incompatibilities, some with the potential to generate sterility. While the possible contribution of faster-male evolution to this misexpression has been explored, the role of the hemizygous X chromosome (i.e., the dominance theory for transcriptomes) remains yet to be determined. Here, we study genome-wide patterns of gene expression in females and males of Drosophila yakuba, Drosophila santomea and their hybrids. We used attached-X stocks to specifically test the dominance theory, and we uncovered a significant contribution of recessive alleles on the X chromosome to hybrid misexpression. Our analyses also suggest a contribution of weakly deleterious regulatory mutations to gene expression divergence in genes with sex-biased expression, but only in the sex toward which the expression is biased (e.g., genes with female-biased expression when analyzed in females). In the opposite sex, we found stronger selective constraints on gene expression divergence. Although genes with a high degree of male-biased expression show a clear signal of faster-X evolution of gene expression, we also detected slower-X evolution in other gene classes (e.g., female-biased genes). This slower-X effect is mediated by significant decreases in cis- and trans-regulatory divergence. The distinct behavior of X-linked genes with a high degree of male-biased expression is consistent with these genes experiencing a higher incidence of positively selected regulatory mutations than their autosomal counterparts.

The Hybrid Incompatibility Genes Lhr and Hmr Are Required for Sister Chromatid Detachment During Anaphase but Not for Centromere Function

Crosses between Drosophila melanogaster females and Drosophila simulans males produce hybrid sons that die at the larval stage. This hybrid lethality is suppressed by loss-of-function mutations in the D. melanogaster Hybrid male rescue (Hmr) or in the D. simulans Lethal hybrid rescue (Lhr) genes. Previous studies have shown that Hmr and Lhr interact with heterochromatin proteins and suppress expression of transposable elements within D. melanogaster It also has been proposed that Hmr and Lhr function at the centromere. We examined mitotic divisions in larval brains from Hmr and Lhr single mutants and Hmr; Lhr double mutants in D. melanogaster In none of the mutants did we observe defects in metaphase chromosome alignment or hyperploid cells, which are hallmarks of centromere or kinetochore dysfunction. In addition, we found that Hmr-HA and Lhr-HA do not colocalize with centromeres either during interphase or mitotic division. However, all mutants displayed anaphase bridges and chromosome aberrations resulting from the breakage of these bridges, predominantly at the euchromatin-heterochromatin junction. The few dividing cells present in hybrid males showed fuzzy and irregularly condensed chromosomes with unresolved sister chromatids. Despite this defect in condensation, chromosomes in hybrids managed to align on the metaphase plate and undergo anaphase. We conclude that there is no evidence for a centromeric function of Hmr and Lhr within D. melanogaster nor for a centromere defect causing hybrid lethality. Instead, we find that Hmr and Lhr are required in D. melanogaster for detachment of sister chromatids during anaphase.

Ectopic application of the repressive histone modification H3K9me2 establishes post-zygotic reproductive isolation in Arabidopsis thaliana

Hybrid seed lethality as a consequence of interspecies or interploidy hybridizations is a major mechanism of reproductive isolation in plants. This mechanism is manifested in the endosperm, a dosage-sensitive tissue supporting embryo growth. Deregulated expression of imprinted genes such as ADMETOS (ADM) underpin the interploidy hybridization barrier in Arabidopsis thaliana; however, the mechanisms of their action remained unknown. In this study, we show that ADM interacts with the AT hook domain protein AHL10 and the SET domain-containing SU(VAR)3-9 homolog SUVH9 and ectopically recruits the heterochromatic mark H3K9me2 to AT-rich transposable elements (TEs), causing deregulated expression of neighboring genes. Several hybrid incompatibility genes identified in Drosophila encode for dosage-sensitive heterochromatin-interacting proteins, which has led to the suggestion that hybrid incompatibilities evolve as a consequence of interspecies divergence of selfish DNA elements and their regulation. Our data show that imbalance of dosage-sensitive chromatin regulators underpins the barrier to interploidy hybridization in Arabidopsis, suggesting that reproductive isolation as a consequence of epigenetic regulation of TEs is a conserved feature in animals and plants.

The Drosophila speciation factor HMR localizes to genomic insulator sites

Hybrid incompatibility between Drosophila melanogaster and D. simulans is caused by a lethal interaction of the proteins encoded by the Hmr and Lhr genes. In D. melanogaster the loss of HMR results in mitotic defects, an increase in transcription of transposable elements and a deregulation of heterochromatic genes. To better understand the molecular mechanisms that mediate HMR’s function, we measured genome-wide localization of HMR in D. melanogaster tissue culture cells by chromatin immunoprecipitation. Interestingly, we find HMR localizing to genomic insulator sites that can be classified into two groups. One group belongs to gypsy insulators and another one borders HP1a bound regions at active genes. The transcription of the latter group genes is strongly affected in larvae and ovaries of Hmr mutant flies. Our data suggest a novel link between HMR and insulator proteins, a finding that implicates a potential role for genome organization in the formation of species.

Gene Regulation and Speciation

Understanding the genetic architecture of speciation is a major goal in evolutionary biology. Hybrid dysfunction is thought to arise most commonly through negative interactions between alleles at two or more loci. Divergence between interacting regulatory elements that affect gene expression (i.e., regulatory divergence) may be a common route for these negative interactions to arise. We review here how regulatory divergence between species can result in hybrid dysfunction, including recent theoretical support for this model. We then discuss the empirical evidence for regulatory divergence between species and evaluate evidence for misregulation as a source of hybrid dysfunction. Finally, we review unresolved questions in gene regulation as it pertains to speciation and point to areas that could benefit from future research.

The Drosophila bag of marbles Gene Interacts Genetically with Wolbachia and Shows Female-Specific Effects of Divergence

Many reproductive proteins from diverse taxa evolve rapidly and adaptively. These proteins are typically involved in late stages of reproduction such as sperm development and fertilization, and are more often functional in males than females. Surprisingly, many germline stem cell (GSC) regulatory genes, which are essential for the earliest stages of reproduction, also evolve adaptively in Drosophila. One example is the bag of marbles (bam) gene, which is required for GSC differentiation and germline cyst development in females and for regulating mitotic divisions and entry to spermatocyte differentiation in males. Here we show that the extensive divergence of bam between Drosophila melanogaster and D. simulans affects bam function in females but has no apparent effect in males. We further find that infection with Wolbachia pipientis, an endosymbiotic bacterium that can affect host reproduction through various mechanisms, partially suppresses female sterility caused by bam mutations in D. melanogaster and interacts differentially with bam orthologs from D. melanogaster and D. simulans. We propose that the adaptive evolution of bam has been driven at least in part by the long-term interactions between Drosophila species and Wolbachia. More generally, we suggest that microbial infections of the germline may explain the unexpected pattern of evolution of several GSC regulatory genes.

A screen for F1 hybrid male rescue reveals no major-effect hybrid lethality loci in the Drosophila melanogaster autosomal genome

Hybrid sons between Drosophila melanogaster females and D. simulans males die as 3rd instar larvae. Two genes, D. melanogaster Hybrid male rescue (Hmr) on the X chromosome, and D. simulans Lethal hybrid rescue (Lhr) on chromosome II, interact to cause this lethality. Loss-of-function mutations in either gene suppress lethality, but several pieces of evidence suggest that additional factors are required for hybrid lethality. Here we screen the D. melanogaster autosomal genome by using the Bloomington Stock Center Deficiency kit to search for additional regions that can rescue hybrid male lethality. Our screen is designed to identify putative hybrid incompatibility (HI) genes similar to Hmr and Lhr which, when removed, are dominant suppressors of lethality. After screening 89% of the autosomal genome, we found no regions that rescue males to the adult stage. We did, however, identify several regions that rescue up to 13% of males to the pharate adult stage. This weak rescue suggests the presence of multiple minor-effect HI loci, but we were unable to map these loci to high resolution, presumably because weak rescue can be masked by genetic background effects. We attempted to test one candidate, the dosage compensation gene male specific lethal-3 (msl-3), by using RNA interference with short hairpin microRNA constructs targeted specifically against D. simulans msl-3 but failed to achieve knockdown, in part due to off-target effects. We conclude that the D. melanogaster autosomal genome likely does not contain additional major-effect HI loci. We also show that Hmr is insufficient to fully account for the lethality associated with the D. melanogaster X chromosome, suggesting that additional X-linked genes contribute to hybrid lethality.

Hybrid incompatibility despite pleiotropic constraint in a sequence-based bioenergetic model of transcription factor binding

Hybrid incompatibility can result from gene misregulation produced by divergence in trans-acting regulatory factors and their cis-regulatory targets. However, change in trans-acting factors may be constrained by pleiotropy, which would in turn limit the evolution of incompatibility. We employed a mechanistically explicit bioenergetic model of gene expression wherein parameter combinations (number of transcription factor molecules, energetic properties of binding to the regulatory site, and genomic background size) determine the shape of the genotype-phenotype (G-P) map, and interacting allelic variants of mutable cis and trans sites determine the phenotype along that map. Misregulation occurs when the phenotype differs from its optimal value. We simulated a pleiotropic regulatory pathway involving a positively selected and a conserved trait regulated by a shared transcription factor (TF), with two populations evolving in parallel. Pleiotropic constraints shifted evolution in the positively selected trait to its cis-regulatory locus. We nevertheless found that the TF genotypes often evolved, accompanied by compensatory evolution in the conserved trait, and both traits contributed to hybrid misregulation. Compensatory evolution resulted in "developmental system drift," whereby the regulatory basis of the conserved phenotype changed although the phenotype itself did not. Pleiotropic constraints became stronger and in some cases prohibitive when the bioenergetic properties of the molecular interaction produced a G-P map that was too steep. Likewise, compensatory evolution slowed and hybrid misregulation was not evident when the G-P map was too shallow. A broad pleiotropic "sweet spot" nevertheless existed where evolutionary constraints were moderate to weak, permitting substantial hybrid misregulation in both traits. None of these pleiotropic constraints manifested when the TF contained nonrecombining domains independently regulating the respective traits.

Hybrid incompatibility arises in a sequence-based bioenergetic model of transcription factor binding

Postzygotic isolation between incipient species results from the accumulation of incompatibilities that arise as a consequence of genetic divergence. When phenotypes are determined by regulatory interactions, hybrid incompatibility can evolve even as a consequence of parallel adaptation in parental populations because interacting genes can produce the same phenotype through incompatible allelic combinations. We explore the evolutionary conditions that promote and constrain hybrid incompatibility in regulatory networks using a bioenergetic model (combining thermodynamics and kinetics) of transcriptional regulation, considering the bioenergetic basis of molecular interactions between transcription factors (TFs) and their binding sites. The bioenergetic parameters consider the free energy of formation of the bond between the TF and its binding site and the availability of TFs in the intracellular environment. Together these determine fractional occupancy of the TF on the promoter site, the degree of subsequent gene expression and in diploids, and the degree of dominance among allelic interactions. This results in a sigmoid genotype-phenotype map and fitness landscape, with the details of the shape determining the degree of bioenergetic evolutionary constraint on hybrid incompatibility. Using individual-based simulations, we subjected two allopatric populations to parallel directional or stabilizing selection. Misregulation of hybrid gene expression occurred under either type of selection, although it evolved faster under directional selection. Under directional selection, the extent of hybrid incompatibility increased with the slope of the genotype-phenotype map near the derived parental expression level. Under stabilizing selection, hybrid incompatibility arose from compensatory mutations and was greater when the bioenergetic properties of the interaction caused the space of nearly neutral genotypes around the stable expression level to be wide. F2's showed higher hybrid incompatibility than F1's to the extent that the bioenergetic properties favored dominant regulatory interactions. The present model is a mechanistically explicit case of the Bateson-Dobzhansky-Muller model, connecting environmental selective pressure to hybrid incompatibility through the molecular mechanism of regulatory divergence. The bioenergetic parameters that determine expression represent measurable properties of transcriptional regulation, providing a predictive framework for empirical studies of how phenotypic evolution results in epistatic incompatibility at the molecular level in hybrids.

The Hmr and Lhr hybrid incompatibility genes suppress a broad range of heterochromatic repeats

Hybrid incompatibilities (HIs) cause reproductive isolation between species and thus contribute to speciation. Several HI genes encode adaptively evolving proteins that localize to or interact with heterochromatin, suggesting that HIs may result from co-evolution with rapidly evolving heterochromatic DNA. Little is known, however, about the intraspecific function of these HI genes, the specific sequences they interact with, or the evolutionary forces that drive their divergence. The genes Hmr and Lhr genetically interact to cause hybrid lethality between Drosophila melanogaster and D. simulans, yet mutations in both genes are viable. Here, we report that Hmr and Lhr encode proteins that form a heterochromatic complex with Heterochromatin Protein 1 (HP1a). Using RNA-Seq analyses we discovered that Hmr and Lhr are required to repress transcripts from satellite DNAs and many families of transposable elements (TEs). By comparing Hmr and Lhr function between D. melanogaster and D. simulans we identify several satellite DNAs and TEs that are differentially regulated between the species. Hmr and Lhr mutations also cause massive overexpression of telomeric TEs and significant telomere lengthening. Hmr and Lhr therefore regulate three types of heterochromatic sequences that are responsible for the significant differences in genome size and structure between D. melanogaster and D. simulans and have high potential to cause genetic conflicts with host fitness. We further find that many TEs are overexpressed in hybrids but that those specifically mis-expressed in lethal hybrids do not closely correlate with Hmr function. Our results therefore argue that adaptive divergence of heterochromatin proteins in response to repetitive DNAs is an important underlying force driving the evolution of hybrid incompatibility genes, but that hybrid lethality likely results from novel epistatic genetic interactions that are distinct to the hybrid background.

Sister species capable of mating often produce hybrids that are sterile or die during development. This reproductive isolation is caused by incompatibilities between the two sister species' genomes. Some hybrid incompatibilities involve genes that encode rapidly evolving proteins that localize to heterochromatin. Heterochromatin is largely made up of highly repetitive transposable elements and satellite DNAs. It has been hypothesized that rapid changes in heterochromatic DNA drives the changes in these HI genes and thus the evolution of reproductive isolation. In support of this model, we show that two rapidly evolving HI proteins, Lhr and Hmr, which reproductively isolate the fruit fly sister species D. melanogaster and D. simulans, repress transposable elements and satellite DNAs. These proteins also help regulate the length of the atypical Drosophila telomeres, which are themselves made of domesticated transposable elements. Our data suggest that these proteins are part of the adaptive machinery that allows the host to respond to changes and increases in heterochromatin and to maintain the activity of genes located within or adjacent to heterochromatin.

Heterochronic meiotic misexpression in an interspecific yeast hybrid

Regulatory changes rapidly accumulate between species, and interspecific hybrids often misexpress genes. Hybrid misexpression, expression levels outside the range of both parental species, can result from cis- and trans-acting regulatory changes that interact abnormally in hybrids. Thus, misexpressed genes may contribute to hybrid sterility. However, in the context of a whole organism, misexpression may not result directly from cis-trans interactions but rather indirectly from differences between hybrid and parental abundance of cell types. Here we eliminate the confounding effects of cell types by examining gene expression in a sterile interspecific yeast hybrid during meiosis. We investigated gene expression of the yeasts Saccharomyces cerevisiae, S. paradoxus, and their hybrid at multiple meiotic stages. Although the hybrid and parents exhibit similar changes in expression levels across meiosis, the hybrid meiotic program occurs earlier than either parent. The timing change produces a heterochronic pattern of misexpression during midmeiosis. Coincident with the timing of misexpression, we find a transition from predominantly trans-acting to cis-acting expression divergence and an increase in the number of opposing cis-trans changes. However, we find no direct relationship between opposing cis-trans changes and misexpression. Contrary to the notion that cis-trans interactions cause misexpression, a heterochronic shift in the normal meiotic gene expression program produces patterns of misexpression in an yeast hybrid. Our results imply that temporal dynamics of single cell types is important to understanding hybrid misexpression and its relationship to cis-trans interactions.

Tempo and mode of regulatory evolution in Drosophila

Genetic changes affecting gene expression contribute to phenotypic divergence; thus, understanding how regulatory networks controlling gene expression change over time is critical for understanding evolution. Prior studies of expression differences within and between species have identified properties of regulatory divergence, but technical and biological differences among these studies make it difficult to assess the generality of these properties or to understand how regulatory changes accumulate with divergence time. Here, we address these issues by comparing gene expression among strains and species of Drosophila with a range of divergence times and use F₁ hybrids to examine inheritance patterns and disentangle cis- and trans-regulatory changes. We find that the fixation of compensatory changes has caused the regulation of gene expression to diverge more rapidly than gene expression itself. Specifically, we observed that the proportion of genes with evidence of cis-regulatory divergence has increased more rapidly with divergence time than the proportion of genes with evidence of expression differences. Surprisingly, the amount of expression divergence explained by cis-regulatory changes did not increase steadily with divergence time, as was previously proposed. Rather, one species (Drosophila sechellia) showed an excess of cis-regulatory divergence that we argue most likely resulted from positive selection in this lineage. Taken together, this work reveals not only the rate at which gene expression evolves, but also the molecular and evolutionary mechanisms responsible for this evolution.

Allelic asymmetry of the Lethal hybrid rescue (Lhr) gene expression in the hybrid between Drosophila melanogaster and D. simulans: confirmation by using genetic variations of D. melanogaster

In the cross between Drosophila melanogaster females and D. simulans males, hybrid males die at the late larval stage, and the sibling females also die at later stages at high temperatures. Removing the D. simulans allele of the Lethal hybrid rescue gene (Lhr (sim) ) improves the hybrid incompatibility phenotypes. However, the loss-of-function mutation of Lhr (sim) (Lhr (sim0) ) does not rescue the hybrid males in crosses with several D. melanogaster strains. We first describe the genetic factor possessed by the D. melanogaster strains. It has been suggested that removing the D. melanogaster allele of Lhr (Lhr (mel) ), that is Lhr (mel0) , does not have the hybrid male rescue effect, contrasting to Lhr (sim0) . Because the expression level of the Lhr gene is known to be Lhr (sim) > Lhr (mel) in the hybrid, Lhr (mel0) may not lead to enough of a reduction in total Lhr expression. Then, there is a possibility that the D. melanogaster factor changes the expression level to Lhr (sim) < Lhr (mel) . But in fact, the expression level was Lhr (sim) > Lhr (mel) in the hybrid irrespectively of the presence of the factor. At last, we showed that Lhr (mel0) slightly improves the viability of hybrid females, which was not realized previously. All of the present results are consistent with the allelic asymmetry model of the Lhr gene expression in the hybrid.

A pair of centromeric proteins mediates reproductive isolation in Drosophila species

Speciation involves the reproductive isolation of natural populations due to the sterility or lethality of their hybrids. However, the molecular basis of hybrid lethality and the evolutionary driving forces that provoke it remain largely elusive. The hybrid male rescue (Hmr) and the lethal hybrid rescue (Lhr) genes serve as a model to study speciation in Drosophilids because their interaction causes lethality in male hybrid offspring. Here, we show that HMR and LHR form a centromeric complex necessary for proper chromosome segregation. We find that the Hmr expression level is substantially higher in Drosophila melanogaster, whereas Lhr expression levels are increased in Drosophila simulans. The resulting elevated amount of HMR/LHR complex in hybrids results in an extensive mislocalization of the complex, an interference with the regulation of transposable elements, and an impairment of cell proliferation. Our findings provide evidence for a major role of centromere divergence in the generation of biodiversity.

Research progress in allele-specific expression and its regulatory mechanisms

Although the majority of genes are expressed equally from both alleles, some genes are differentially expressed. Organisms possess characteristics to preferentially express a particular allele under regulatory factors, which is termed allele-specific expression (ASE). It is one of the important genetic factors that lead to phenotypic variation and can be used to identify the variance of gene regulation factors. ASE indicates mechanisms such as DNA methylation, histone modifications, and non-coding RNAs function. Here, we review a broad survey of progress in ASE studies, and what this simple yet very effective approach can offer in functional genomics, and possible implications toward our better understanding of the underlying mechanisms of complex traits.

Genome-wide patterns of natural variation reveal strong selective sweeps and ongoing genomic conflict in Drosophila mauritiana

Although it is well understood that selection shapes the polymorphism pattern in Drosophila, signatures of classic selective sweeps are scarce. Here, we focus on Drosophila mauritiana, an island endemic, which is closely related to Drosophila melanogaster. Based on a new, annotated genome sequence, we characterized the genome-wide polymorphism by sequencing pooled individuals (Pool-seq). We show that the interplay between selection and recombination results in a genome-wide polymorphism pattern characteristic for D. mauritiana. Two large genomic regions (>500 kb) showed the signature of almost complete selective sweeps. We propose that the absence of population structure and limited geographic distribution could explain why such pronounced sweep patterns are restricted to D. mauritiana. Further evidence for strong adaptive evolution was detected for several nucleoporin genes, some of which were not previously identified as genes involved in genomic conflict. Since this adaptive evolution is continuing after the split of D. mauritiana and Drosophila simulans, we conclude that genomic conflict is not restricted to short episodes, but rather an ongoing process in Drosophila.

Drosophila interspecific hybrids phenocopy piRNA-pathway mutants

Hybrids of two Drosophila species show transposable element derepression and piRNA pathway malfunction, revealing adaptive evolution of piRNA pathway components.

The Piwi-interacting RNA (piRNA) pathway defends the germline of animals from the deleterious activity of selfish transposable elements (TEs) through small-RNA mediated silencing. Adaptation to novel invasive TEs is proposed to occur by incorporating their sequences into the piRNA pool that females produce and deposit into their eggs, which then propagates immunity against specific TEs to future generations. In support of this model, the F1 offspring of crosses between strains of the same Drosophila species sometimes suffer from germline derepression of paternally inherited TE families, caused by a failure of the maternal strain to produce the piRNAs necessary for their regulation. However, many protein components of the Drosophila piRNA pathway exhibit signatures of positive selection, suggesting that they also contribute to the evolution of host genome defense. Here we investigate piRNA pathway function and TE regulation in the F1 hybrids of interspecific crosses between D. melanogaster and D. simulans and compare them with intraspecific control crosses of D. melanogaster. We confirm previous reports showing that intraspecific crosses are characterized by derepression of paternally inherited TE families that are rare or absent from the maternal genome and piRNA pool, consistent with the role of maternally deposited piRNAs in shaping TE silencing. In contrast to the intraspecific cross, we discover that interspecific hybrids are characterized by widespread derepression of both maternally and paternally inherited TE families. Furthermore, the pattern of derepression of TE families in interspecific hybrids cannot be attributed to their paucity or absence from the piRNA pool of the maternal species. Rather, we demonstrate that interspecific hybrids closely resemble piRNA effector-protein mutants in both TE misregulation and aberrant piRNA production. We suggest that TE derepression in interspecific hybrids largely reflects adaptive divergence of piRNA pathway genes rather than species-specific differences in TE-derived piRNAs.

Eukaryotic genomes contain large quantities of transposable elements (TEs), short self-replicating DNA sequences that can move within the genome. The selfish replication of TEs has potentially drastic consequences for the host, such as disruption of gene function, induction of sterility, and initiation or exacerbation of some cancers. Like the adaptive immune system that defends our bodies against pathogens, the Piwi-interacting RNA (piRNA) pathway defends animal genomes against the harmful effects of TEs. Fundamental to piRNA-mediated defense is the production of small noncoding RNAs that act like antibodies to target replicating TEs for destruction by piRNA-effector proteins. piRNAs are expected to diverge rapidly between species in response to genome infection by increasingly disparate TEs. Here, we tested this hypothesis by examining how differences in piRNAs between two species of fruit fly relate to TE “immunity” in their hybrid offspring. Because piRNAs are maternally deposited, we expected excessive replication of paternal TEs in hybrids. Surprisingly, we observe increased activity of both maternal and paternal TEs, together with defects in piRNA production that are reminiscent of piRNA effector-protein mutants. Our observations reveal that piRNA effector-proteins do not function properly in hybrids, and we propose that adaptive evolution among piRNA effector-proteins contributes to host genome defense and leads to the functional incompatibilities that we observe in hybrids.

Interallelic and intergenic incompatibilities of the Prdm9 (Hst1) gene in mouse hybrid sterility

The Dobzhansky-Muller model of incompatibilities explains reproductive isolation between species by incorrect epistatic interactions. Although the mechanisms of speciation are of great interest, no incompatibility has been characterized at the gene level in mammals. The Hybrid sterility 1 gene (Hst1) participates in the arrest of meiosis in F₁ males of certain strains from two Mus musculus subspecies, e.g., PWD from M. m. musculus and C57BL/6J (henceforth B6) from M. m. domesticus. Hst1 has been identified as a meiotic PR-domain gene (Prdm9) encoding histone 3 methyltransferase in the male offspring of PWD females and B6 males, (PWD×B6)F₁. To characterize the incompatibilities underlying hybrid sterility, we phenotyped reproductive and meiotic markers in males with altered copy numbers of Prdm9. A partial rescue of fertility was observed upon removal of the B6 allele of Prdm9 from the azoospermic (PWD×B6)F₁ hybrids, whereas removing one of the two Prdm9 copies in PWD or B6 background had no effect on male reproduction. Incompatibility(ies) not involving Prdm9^B6 also acts in the (PWD×B6)F₁ hybrids, since the correction of hybrid sterility by Prdm9^B6 deletion was not complete. Additions and subtractions of Prdm9 copies, as well as allelic replacements, improved meiotic progression and fecundity also in the progeny-producing reciprocal (B6×PWD)F₁ males. Moreover, an increased dosage of Prdm9 and reciprocal cross enhanced fertility of other sperm-carrying male hybrids, (PWD×B6-C3H.Prdm9)F₁, harboring another Prdm9 allele of M. m. domesticus origin. The levels of Prdm9 mRNA isoforms were similar in the prepubertal testes of all types of F₁ hybrids of PWD with B6 and B6-C3H.Prdm9 despite their different prospective fertility, but decreased to 53% after removal of Prdm9^B6. Therefore, the Prdm9^B6 allele probably takes part in posttranscriptional dominant-negative hybrid interaction(s) absent in the parental strains.

Disturbed gametogenesis in the progeny of two fertile parental forms is called hybrid sterility; it is an important part of reproductive barriers between species. The Dobzhansky-Muller model of incompatibilities explains reproductive isolation between species by incorrect interactions between genes. Hybrid sterility 1 (Hst1) is one of the genes causing meiotic arrest in F₁ male hybrids between certain Mus musculus musculus (e.g., the PWD strain) and M. m. domesticus (C57BL/6J etc.) mice. Hst1, the first mammalian candidate for a speciation gene, was identified as a meiotic PR/SET-domain gene, Prdm9, but the mechanism causing sterility has remained unknown. While the F₁ male offspring of C57BL/6J males and PWD females produce no sperm, the males from the reciprocal cross using PWD males and C57BL/6J females yield progeny. Here we show that the meiotic progress and fertility of hybrid males from both F₁ crosses improved by removal as well as overexpression of the C57BL/6J allele of Prdm9, suggesting that Prdm9 interactions not present in the parental species (incompatibilities) play a role in hybrid sterility. Furthermore, the Prdm9 dosage also controlled fecundity in other F₁ hybrids, indicating that this gene is an important regulator of mouse hybrid fertility.

Coevolution within and between regulatory loci can preserve promoter function despite evolutionary rate acceleration

Phenotypes that appear to be conserved could be maintained not only by strong purifying selection on the underlying genetic systems, but also by stabilizing selection acting via compensatory mutations with balanced effects. Such coevolution has been invoked to explain experimental results, but has rarely been the focus of study. Conserved expression driven by the unc-47 promoters of Caenorhabditis elegans and C. briggsae persists despite divergence within a cis-regulatory element and between this element and the trans-regulatory environment. Compensatory changes in cis and trans are revealed when these promoters are used to drive expression in the other species. Functional changes in the C. briggsae promoter, which has experienced accelerated sequence evolution, did not lead to alteration of gene expression in its endogenous environment. Coevolution among promoter elements suggests that complex epistatic interactions within cis-regulatory elements may facilitate their divergence. Our results offer a detailed picture of regulatory evolution in which subtle, lineage-specific, and compensatory modifications of interacting cis and trans regulators together maintain conserved gene expression patterns.

Some phenotypes, including gene expression patterns, are conserved between distantly related species. However, the molecular bases of those phenotypes are not necessarily conserved. Instead, regulatory DNA sequences and the proteins with which they interact can change over time with balanced effects, preserving expression patterns and concealing regulatory divergence. Coevolution between interacting molecules makes gene regulation highly species-specific, and it can be detected when the cis-regulatory DNA of one species is used to drive expression in another species. In this way, we identified regions of the C. elegans and C. briggsae unc-47 promoters that have coevolved with the lineage-specific trans-regulatory environments of these organisms. The C. briggsae promoter experienced accelerated sequence change relative to related species. All of this evolution occurred without changing the expression pattern driven by the promoter in its endogenous environment.

An indel polymorphism in the hybrid incompatibility gene lethal hybrid rescue of Drosophila is functionally relevant

Hybrid incompatibility (HI) genes are frequently observed to be rapidly evolving under selection. This observation has led to the attractive conjecture that selection-derived protein-sequence divergence is culpable for incompatibilities in hybrids. The Drosophila simulans HI gene Lethal hybrid rescue (Lhr) is an intriguing case, because despite having experienced rapid sequence evolution, its HI properties are a shared function inherited from the ancestral state. Using an unusual D. simulans Lhr hybrid rescue allele, Lhr(2), we here identify a conserved stretch of 10 amino acids in the C terminus of LHR that is critical for causing hybrid incompatibility. Altering these 10 amino acids weakens or abolishes the ability of Lhr to suppress the hybrid rescue alleles Lhr(1) or Hmr(1), respectively. Besides single-amino-acid substitutions, Lhr orthologs differ by a 16-aa indel polymorphism, with the ancestral deletion state fixed in D. melanogaster and the derived insertion state at very high frequency in D. simulans. Lhr(2) is a rare D. simulans allele that has the ancestral deletion state of the 16-aa polymorphism. Through a series of transgenic constructs we demonstrate that the ancestral deletion state contributes to the rescue activity of Lhr(2). This indel is thus a polymorphism that can affect the HI function of Lhr.