Recombination and selection against introgressed DNA.

Introgressed DNA is often deleterious at many loci in the recipient species' genome, and is therefore purged by selection. Here, we use mathematical modeling and whole-genome simulations to study the influence of recombination on this process. We find that aggregate recombination controls the genome-wide rate of purging in the early generations after admixture, when purging is most rapid. Aggregate recombination is influenced by the number of chromosomes and heterogeneity in their size, and by the number of crossovers and their locations along chromosomes. A comparative prediction is that species with fewer chromosomes should purge introgressed ancestry more profoundly, and should therefore exhibit weaker genomic signals of historical introgression. Turning to within-genome patterns, we show that, in species with autosomal recombination in both sexes, more purging is expected on sex chromosomes than autosomes, all else equal. The opposite prediction holds for species without autosomal recombination in the heterogametic sex. Finally, positive correlations between recombination rate and introgressed ancestry have recently been observed within the genomes of several species. We show that these correlations are likely driven not by recombination's effect in unlinking neutral from deleterious introgressed alleles, but by recombination's effect on the rate of purging of deleterious introgressed alleles themselves.

Complex basis of hybrid female sterility and Haldane's rule in Heliconius butterflies: Z-linkage and epistasis.

Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.

Prevalence and Adaptive Impact of Introgression.

Alleles that introgress between species can influence the evolutionary and ecological fate of species exposed to novel environments. Hybrid offspring of different species are often unfit, and yet it has long been argued that introgression can be a potent force in evolution, especially in plants. Over the last two decades, genomic data have increasingly provided evidence that introgression is a critically important source of genetic variation and that this additional variation can be useful in adaptive evolution of both animals and plants. Here, we review factors that influence the probability that foreign genetic variants provide long-term benefits (so-called adaptive introgression) and discuss their potential benefits. We find that introgression plays an important role in adaptive evolution, particularly when a species is far from its fitness optimum, such as when they expand their range or are subject to changing environments.

Synteny-Based Genome Assembly for 16 Species of Heliconius Butterflies, and an Assessment of Structural Variation across the Genus.

Heliconius butterflies (Lepidoptera: Nymphalidae) are a group of 48 neotropical species widely studied in evolutionary research. Despite the wealth of genomic data generated in past years, chromosomal level genome assemblies currently exist for only two species, Heliconius melpomene and Heliconius erato, each a representative of one of the two major clades of the genus. Here, we use these reference genomes to improve the contiguity of previously published draft genome assemblies of 16 Heliconius species. Using a reference-assisted scaffolding approach, we place and order the scaffolds of these genomes onto chromosomes, resulting in 95.7-99.9% of their genomes anchored to chromosomes. Genome sizes are somewhat variable among species (270-422 Mb) and in one small group of species (Heliconius hecale, Heliconius elevatus, and Heliconius pardalinus) expansions in genome size are driven mainly by repetitive sequences that map to four small regions in the H. melpomene reference genome. Genes from these repeat regions show an increase in exon copy number, an absence of internal stop codons, evidence of constraint on nonsynonymous changes, and increased expression, all of which suggest that at least some of the extra copies are functional. Finally, we conducted a systematic search for inversions and identified five moderately large inversions fixed between the two major Heliconius clades. We infer that one of these inversions was transferred by introgression between the lineages leading to the erato/sara and burneyi/doris clades. These reference-guided assemblies represent a major improvement in Heliconius genomic resources that enable further genetic and evolutionary discoveries in this genus.

De novo chromosome-length assembly of the mule deer (Odocoileus hemionus) genome.

The mule deer (Odocoileus hemionus) is an ungulate species that is distributed in a range from western Canada to central Mexico. Mule deer are an essential source of food for many predators, are relatively abundant, and commonly make broad migration movements. A clearer understanding of the mule deer genome can improve our knowledge of its population genetics, movements, and demographic history, aiding in conservation efforts. Their large population size, continuous distribution, and diversity of habitat make mule deer excellent candidates for population genomics studies; however, few genomic resources are currently available for this species. Here, we sequence and assemble the mule deer genome into a highly contiguous chromosome-length assembly for use in future research using long-read sequencing and Hi-C technologies. We also provide a genome annotation and compare demographic histories of the mule deer and white-tailed deer using the pairwise sequentially Markovian coalescent model. We expect this assembly to be a valuable resource in the continued study and conservation of mule deer.

Variation in Genetic Relatedness Is Determined by the Aggregate Recombination Process.

The genomic proportion that two relatives share identically by descent-their genetic relatedness-can vary depending on the history of recombination and segregation in their pedigree. Previous calculations of the variance of genetic relatedness have defined genetic relatedness as the proportion of total genetic map length (cM) shared by relatives, and have neglected crossover interference and sex differences in recombination. Here, we consider genetic relatedness as the proportion of the total physical genome (bp) shared by relatives, and calculate its variance for general pedigree relationships, making no assumptions about the recombination process. For the relationships of grandparent-grandoffspring and siblings, the variance of genetic relatedness is a simple decreasing function of [Formula: see text], the average proportion of locus pairs that recombine in meiosis. For general pedigree relationships, the variance of genetic relatedness is a function of metrics analogous to [Formula: see text] Therefore, features of the aggregate recombination process that affect [Formula: see text] and analogs also affect variance in genetic relatedness. Such features include the number of chromosomes and heterogeneity in their size, the number of crossovers and their spatial organization along chromosomes, and sex differences in recombination. Our calculations help to explain several recent observations about variance in genetic relatedness, including that it is reduced by crossover interference (which is known to increase [Formula: see text]). Our methods further allow us to calculate the neutral variance of ancestry among F2s in a hybrid cross, enabling precise statistical inference in F2-based tests for various kinds of selection.

Genomic architecture and introgression shape a butterfly radiation.

We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.

Geographic contrasts between pre- and postzygotic barriers are consistent with reinforcement in Heliconius butterflies.

Identifying the traits causing reproductive isolation and the order in which they evolve is fundamental to understanding speciation. Here, we quantify prezygotic and intrinsic postzygotic isolation among allopatric, parapatric, and sympatric populations of the butterflies Heliconius elevatus and Heliconius pardalinus. Sympatric populations from the Amazon (H. elevatus and H. p. butleri) exhibit strong prezygotic isolation and rarely mate in captivity; however, hybrids are fertile. Allopatric populations from the Amazon (H. p. butleri) and Andes (H. p. sergestus) mate freely when brought together in captivity, but the female F1 hybrids are sterile. Parapatric populations (H. elevatus and H. p. sergestus) exhibit both assortative mating and sterility of female F1s. Assortative mating in sympatric populations is consistent with reinforcement in the face of gene flow, where the driving force, selection against hybrids, is due to disruption of mimicry and other ecological traits rather than hybrid sterility. In contrast, the lack of assortative mating and hybrid sterility observed in allopatric populations suggests that geographic isolation enables the evolution of intrinsic postzygotic reproductive isolation. Our results show how the types of reproductive barriers that evolve between species may depend on geography.

Simultaneous TE Analysis of 19 Heliconiine Butterflies Yields Novel Insights into Rapid TE-Based Genome Diversification and Multiple SINE Births and Deaths.

Transposable elements (TEs) play major roles in the evolution of genome structure and function. However, because of their repetitive nature, they are difficult to annotate and discovering the specific roles they may play in a lineage can be a daunting task. Heliconiine butterflies are models for the study of multiple evolutionary processes including phenotype evolution and hybridization. We attempted to determine how TEs may play a role in the diversification of genomes within this clade by performing a detailed examination of TE content and accumulation in 19 species whose genomes were recently sequenced. We found that TE content has diverged substantially and rapidly in the time since several subclades shared a common ancestor with each lineage harboring a unique TE repertoire. Several novel SINE lineages have been established that are restricted to a subset of species. Furthermore, the previously described SINE, Metulj, appears to have gone extinct in two subclades while expanding to significant numbers in others. This diversity in TE content and activity has the potential to impact how heliconiine butterflies continue to evolve and diverge.

No evidence for intracellular magnetite in putative vertebrate magnetoreceptors identified by magnetic screening.

The cellular basis of the magnetic sense remains an unsolved scientific mystery. One theory that aims to explain how animals detect the magnetic field is the magnetite hypothesis. It argues that intracellular crystals of the iron oxide magnetite (Fe3O4) are coupled to mechanosensitive channels that elicit neuronal activity in specialized sensory cells. Attempts to find these primary sensors have largely relied on the Prussian Blue stain that labels cells rich in ferric iron. This method has proved problematic as it has led investigators to conflate iron-rich macrophages with magnetoreceptors. An alternative approach developed by Eder et al. [Eder SH, et al. (2012) Proc Natl Acad Sci USA 109(30):12022-12027] is to identify candidate magnetoreceptive cells based on their magnetic moment. Here, we explore the utility of this method by undertaking a screen for magnetic cells in the pigeon. We report the identification of a small number of cells (1 in 476,000) with large magnetic moments (8-106 fAm(2)) from various tissues. The development of single-cell correlative light and electron microscopy (CLEM) coupled with electron energy loss spectroscopy (EELS) and energy-filtered transmission electron microscopy (EFTEM) permitted subcellular analysis of magnetic cells. This revealed the presence of extracellular structures composed of iron, titanium, and chromium accounting for the magnetic properties of these cells. Application of single-cell CLEM to magnetic cells from the trout failed to identify any intracellular structures consistent with biogenically derived magnetite. Our work illustrates the need for new methods to test the magnetite hypothesis of magnetosensation.

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.

Drosophila interspecific hybrids phenocopy piRNA-pathway mutants.

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.