Daniel A. Barbash

Eukaryotic genomes are replete with repeated sequences in the form of transposable elements (TEs) dispersed across the genome or as satellite arrays, large stretches of tandemly repeated sequences. Many satellites clearly originated as TEs, but it is unclear how mobile genetic parasites can transform into megabase-sized tandem arrays. Comprehensive population genomic sampling is needed to determine the frequency and generative mechanisms of tandem TEs, at all stages from their initial formation to their subsequent expansion and maintenance as satellites. The best available population resources, short-read DNA sequences, are often considered to be of limited utility for analyzing repetitive DNA due to the challenge of mapping individual repeats to unique genomic locations. Here we develop a new pipeline called ConTExt that demonstrates that paired-end Illumina data can be successfully leveraged to identify a wide range of structural variation within repetitive sequence, including tandem elements. By analyzing 85 genomes from five populations ofDrosophila melanogaster, we discover that TEs commonly form tandem dimers. Our results further suggest that insertion site preference is the major mechanism by which dimers arise and that, consequently, dimers form rapidly during periods of active transposition. This abundance of TE dimers has the potential to provide source material for future expansion into satellite arrays, and we discover one such copy number expansion of the DNA transposonhoboto approximately 16 tandem copies in a single line. The very process that defines TEs—transposition—thus regularly generates sequences from which new satellites can arise.

The question of how new species evolve has been examined at every level, from macroevolutionary patterns of diversification to molecular population genetic analyses of specific genomic regions between species pairs. Drosophila has been at the center of many of these research efforts. Though our understanding of the speciation process has grown considerably over the past few decades, very few genes have been identified that contribute to barriers to reproduction. The development of advanced molecular genetic and genomic methods provides promising avenues for the rapid discovery of more genes that contribute to speciation, particularly those involving prezygotic isolation. The continued expansion of tools and resources, especially for species other than Drosophila melanogaster, will be most effective when coupled with comparative approaches that reveal the genetic basis of reproductive isolation across a range of divergence times. Future research programs in Drosophila have high potential to answer long-standing questions in speciation. These include identifying the selective forces that contribute to divergence between populations and the genetic basis of traits that cause reproductive isolation. The latter can be expanded upon to understand how the genetic basis of reproductive isolation changes over time and whether certain pathways and genes are more commonly involved.

Eukaryotic genomes are replete with repeated sequences, in the form of transposable elements (TEs) dispersed across the genome or as satellite arrays, large stretches of tandemly repeated sequence. Many satellites clearly originated as TEs, but it is unclear how mobile genetic parasites can transform into megabase-sized tandem arrays. Comprehensive population genomic sampling is needed to determine the frequency and generative mechanisms of tandem TEs, at all stages from their initial formation to their subsequent expansion and maintenance as satellites. The best available population resources, short-read DNA sequences, are often considered to be of limited utility for analyzing repetitive DNA due to the challenge of mapping individual repeats to unique genomic locations. Here we develop a new pipeline called ConTExt which demonstrates that paired-end Illumina data can be successfully leveraged to identify a wide range of structural variation within repetitive sequence, including tandem elements. Analyzing 85 genomes from five populations of Drosophila melanogaster we discover that TEs commonly form tandem dimers. Our results further suggest that insertion site preference is the major mechanism by which dimers arise and that, consequently, dimers form rapidly during periods of active transposition. This abundance of TE dimers has the potential to provide source material for future expansion into satellite arrays, and we discover one such copy number expansion of the DNA transposon Hobo to ~16 tandem copies in a single line. The very process that defines TEs —transposition— thus regularly generates sequences from which new satellites can arise.

Mendel’s Law of equal segregation states that each of the two alleles in a diploid organism has an equal probability of being transmitted into...Meiotic drive occurs when a selfish element increases its transmission frequency above the Mendelian ratio by hijacking the asymmetric divisions of female meiosis. Meiotic drive causes genomic conflict and potentially has a major impact on genome evolution, but only a few drive loci of large effect have been described. New methods to reliably detect meiotic drive are therefore needed, particularly for discovering moderate-strength drivers that are likely to be more prevalent in natural populations than strong drivers. Here, we report an efficient method that uses sequencing of large pools of backcross (BC1) progeny to test for deviations from Mendelian segregation genome-wide with single-nucleotide polymorphisms (SNPs) that distinguish the parental strains. We show that meiotic drive can be detected by a characteristic pattern of decay in distortion of SNP frequencies, caused by recombination unlinking the driver from distal loci. We further show that control crosses allow allele-frequency distortion caused by meiotic drive to be distinguished from distortion resulting from developmental effects. We used this approach to test whether chromosomes with extreme telomere-length differences segregate at Mendelian ratios, as telomeric regions are a potential hotspot for meiotic drive due to their roles in meiotic segregation and multiple observations of high rates of telomere sequence evolution. Using four different pairings of long and short telomere strains, we find no evidence that extreme telomere-length variation causes meiotic drive in Drosophila. However, we identify one candidate meiotic driver in a centromere-linked region that shows an ∼8% increase in transmission frequency, corresponding to a ∼54:46 segregation ratio. Our results show that candidate meiotic drivers of moderate strength can be readily detected and localized in pools of BC1 progeny.

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. Animals need to make gametes–sperm or eggs–in order to reproduce. Gametes are produced from a specialized tissue called the germline that is found within the testes or ovaries. These organs contain a small population of stem cells that are able to both self-renew and differentiate to generate gametes and are thus essential for maintaining gamete production throughout the reproductive lifespan of most animals. Surprisingly, some of the genes that control this process evolve rapidly between Drosophila species. We find for a key germline stem cell regulatory gene, bag of marbles (bam), that its rapid evolution affects only female but not male functions. We further report that the endosymbiont bacterium Wolbachia that infects insects and other species interacts with bam and may be contributing to the wider pattern of rapid evolution of germline stem cell regulatory genes.

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.

Misregulation of gene expression is often observed in interspecific hybrids and is generally attributed to regulatory incompatibilities caused by divergence between the two genomes. However, it has been challenging to distinguish effects of regulatory divergence from secondary effects including developmental and physiological defects common to hybrids. Here, we use RNA-Seq to profile gene expression in F1 hybrid male larvae from crosses of Drosophila melanogaster to its sibling species D. simulans. We analyze lethal and viable hybrid males, the latter produced using a mutation in the X-linked D. melanogaster Hybrid male rescue (Hmr) gene and compare them with their parental species and to public data sets of gene expression across development. We find that Hmr has drastically different effects on the parental and hybrid genomes, demonstrating that hybrid incompatibility genes can exhibit novel properties in the hybrid genetic background. Additionally, we find that D. melanogaster alleles are preferentially affected between lethal and viable hybrids. We further determine that many of the differences between the hybrids result from developmental delay in the Hmr⁺ hybrids. Finally, we find surprisingly modest expression differences in hybrids when compared with the parents, with only 9% and 4% of genes deviating from additivity or expressed outside of the parental range, respectively. Most of these differences can be attributed to developmental delay and differences in tissue types. Overall, our study suggests that hybrid gene misexpression is prone to overestimation and that even between species separated by approximately 2.5 Ma, regulatory incompatibilities are not widespread in hybrids.

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.

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.

The Dobzhansky and Muller (D-M) model explains the evolution of hybrid incompatibility (HI) through the interaction between lineage-specific derived alleles at two or more loci. In agreement with the expectation that HI results from functional divergence, many protein-coding genes that contribute to incompatibilities between species show signatures of adaptive evolution, including Lhr, which encodes a heterochromatin protein whose amino acid sequence has diverged extensively between Drosophila melanogaster and D. simulans by natural selection. The lethality of D. melanogaster/D. simulans F1 hybrid sons is rescued by removing D. simulans Lhr, but not D. melanogaster Lhr, suggesting that the lethal effect results from adaptive evolution in the D. simulans lineage. It has been proposed that adaptive protein divergence in Lhr reflects antagonistic coevolution with species-specific heterochromatin sequences and that defects in LHR protein localization cause hybrid lethality. Here we present surprising results that are inconsistent with this coding-sequence-based model. Using Lhr transgenes expressed under native conditions, we find no evidence that LHR localization differs between D. melanogaster and D. simulans, nor do we find evidence that it mislocalizes in their interspecific hybrids. Rather, we demonstrate that Lhr orthologs are differentially expressed in the hybrid background, with the levels of D. simulans Lhr double that of D. melanogaster Lhr. We further show that this asymmetric expression is caused by cis-by-trans regulatory divergence of Lhr. Therefore, the non-equivalent hybrid lethal effects of Lhr orthologs can be explained by asymmetric expression of a molecular function that is shared by both orthologs and thus was presumably inherited from the ancestral allele of Lhr. We present a model whereby hybrid lethality occurs by the interaction between evolutionarily ancestral and derived alleles. When two different species mate, the hybrid progeny are often sterile or lethal. Such hybrid incompatibilities cause reproductive isolation between species and are an important mechanism for maintaining species as separate units. A gene called Lethal hybrid rescue (Lhr) is part of the cause of hybrid lethality between Drosophila species. Like many other hybrid incompatibility genes, Lhr protein sequences in the hybridizing species have diverged from one another by natural selection. This and other findings led to the hypotheses that the function of Lhr has changed between the two species, and this is what makes Lhr a hybrid lethality gene. Using a series of genetic, molecular, and cytological assays, we report evidence contrary to these hypotheses, that hybrid lethal activity is instead a function shared by both species and inherited from their common ancestor. This result is particularly surprising because the Lhr genes from the two species have different effects on hybrid viability. We discovered that these differential effects are caused by differences in expression levels of Lhr in hybrids rather than by changes in its protein-coding sequence. Our results demonstrate that, while natural selection may be important in evolving hybrid incompatibilities, how it does so in this case remains mysterious.