Demographically explicit scans for barriers to gene flow using gIMble.

Identifying regions of the genome that act as barriers to gene flow between recently diverged taxa has remained challenging given the many evolutionary forces that generate variation in genetic diversity and divergence along the genome, and the stochastic nature of this variation. Progress has been impeded by a conceptual and methodological divide between analyses that infer the demographic history of speciation and genome scans aimed at identifying locally maladaptive alleles i.e. genomic barriers to gene flow. Here we implement genomewide IM blockwise likelihood estimation (gIMble), a composite likelihood approach for the quantification of barriers, that bridges this divide. This analytic framework captures background selection and selection against barriers in a model of isolation with migration (IM) as heterogeneity in effective population size (Ne) and effective migration rate (me), respectively. Variation in both effective demographic parameters is estimated in sliding windows via pre-computed likelihood grids. gIMble includes modules for pre-processing/filtering of genomic data and performing parametric bootstraps using coalescent simulations. To demonstrate the new approach, we analyse data from a well-studied pair of sister species of tropical butterflies with a known history of post-divergence gene flow: Heliconius melpomene and H. cydno. Our analyses uncover both large-effect barrier loci (including well-known wing-pattern genes) and a genome-wide signal of a polygenic barrier architecture.

Chromosome Fissions and Fusions Act as Barriers to Gene Flow between Brenthis Fritillary Butterflies.

Chromosome rearrangements are thought to promote reproductive isolation between incipient species. However, it is unclear how often, and under what conditions, fission and fusion rearrangements act as barriers to gene flow. Here we investigate speciation between two largely sympatric fritillary butterflies, Brenthis daphne and Brenthis ino. We use a composite likelihood approach to infer the demographic history of these species from whole-genome sequence data. We then compare chromosome-level genome assemblies of individuals from each species and identify a total of nine chromosome fissions and fusions. Finally, we fit a demographic model where effective population sizes and effective migration rate vary across the genome, allowing us to quantify the effects of chromosome rearrangements on reproductive isolation. We show that chromosomes involved in rearrangements experienced less effective migration since the onset of species divergence and that genomic regions near rearrangement points have a further reduction in effective migration rate. Our results suggest that the evolution of multiple rearrangements in the B. daphne and B. ino populations, including alternative fusions of the same chromosomes, have resulted in a reduction in gene flow. Although fission and fusion of chromosomes are unlikely to be the only processes that have led to speciation between these butterflies, this study shows that these rearrangements can directly promote reproductive isolation and may be involved in speciation when karyotypes evolve quickly.

Experimental introgression in Drosophila: Asymmetric postzygotic isolation associated with chromosomal inversions and an incompatibility locus on the X chromosome.

Interspecific gene flow (introgression) is an important source of new genetic variation, but selection against it can reinforce reproductive barriers between interbreeding species. We used an experimental approach to trace the role of chromosomal inversions and incompatibility genes in preventing introgression between two partly sympatric Drosophila virilis group species, D. flavomontana and D. montana. We backcrossed F1 hybrid females from a cross between D. flavomontana female and D. montana male with the males of the parental species for two generations and sequenced pools of parental strains and their reciprocal second generation backcross (BC2 mon and BC2 fla) females. Contrasting the observed amount of introgression (mean hybrid index, HI) in BC2 female pools along the genome to simulations under different scenarios allowed us to identify chromosomal regions of restricted and increased introgression. We found no deviation from the HI expected under a neutral null model for any chromosome for the BC2 mon pool, suggesting no evidence for genetic incompatibilities in backcrosses towards D. montana. In contrast, the BC2 fla pool showed high variation in the observed HI between different chromosomes, and massive reduction of introgression on the X chromosome (large X-effect). This observation is compatible with reduced recombination combined with at least one dominant incompatibility locus residing within the X inversion(s). Overall, our study suggests that genetic incompatibilities arising within chromosomal inversions can play an important role in speciation.

Males That Silence Their Father's Genes: Genomic Imprinting of a Complete Haploid Genome.

Genetic conflict is considered a key driver in the evolution of reproductive systems with non-Mendelian inheritance, where parents do not contribute equally to the genetic makeup of their offspring. One of the most extraordinary examples of non-Mendelian inheritance is paternal genome elimination (PGE), a form of haplodiploidy which has evolved repeatedly across arthropods. Under PGE, males are diploid but only transmit maternally inherited chromosomes, while the paternally inherited homologues are excluded from sperm. This asymmetric inheritance is thought to have evolved through an evolutionary arms race between the paternal and maternal genomes over transmission to future generations. In several PGE clades, such as the mealybugs (Hemiptera: Pseudococcidae), paternal chromosomes are not only eliminated from sperm, but also heterochromatinized early in development and thought to remain inactive, which could result from genetic conflict between parental genomes. Here, we present a parent-of-origin allele-specific transcriptome analysis in male mealybugs showing that expression is globally biased toward the maternal genome. However, up to 70% of somatically expressed genes are to some degree paternally expressed, while paternal genome expression is much more restricted in the male reproductive tract, with only 20% of genes showing paternal contribution. We also show that parent-of-origin-specific gene expression patterns are remarkably similar across genotypes, and that genes with completely biparental expression show elevated rates of molecular evolution. Our results provide the clearest example yet of genome-wide genomic imprinting in insects and enhance our understanding of PGE, which will aid future empirical tests of evolutionary theory regarding the origin of this unusual reproductive strategy.

Sex-specific expression and DNA methylation in a species with extreme sexual dimorphism and paternal genome elimination.

Phenotypic differences between sexes are often mediated by differential expression and alternative splicing of genes. However, the mechanisms that regulate these expression and splicing patterns remain poorly understood. The mealybug, Planococcus citri, displays extreme sexual dimorphism and exhibits an unusual instance of sex-specific genomic imprinting, paternal genome elimination (PGE), in which the paternal chromosomes in males are highly condensed and eliminated from the sperm. Planococcus citri has no sex chromosomes and both sexual dimorphism and PGE are predicted to be under epigenetic control. We recently showed that P. citri females display a highly unusual DNA methylation profile for an insect species, with the presence of promoter methylation associated with lower levels of gene expression. Here, we therefore decided to explore genome-wide differences in DNA methylation between male and female P. citri using whole-genome bisulphite sequencing. We identified extreme differences in genome-wide levels and patterns between the sexes. Males display overall higher levels of DNA methylation which manifest as more uniform low levels across the genome. Whereas females display more targeted high levels of methylation. We suggest these unique sex-specific differences are due to chromosomal differences caused by PGE and may be linked to possible ploidy compensation. Using RNA-Seq, we identify extensive sex-specific gene expression and alternative splicing, but we find no correlation with cis-acting DNA methylation.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species.

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.

Secretion of an Argonaute protein by a parasitic nematode and the evolution of its siRNA guides.

Extracellular RNA has been proposed to mediate communication between cells and organisms however relatively little is understood regarding how specific sequences are selected for export. Here, we describe a specific Argonaute protein (exWAGO) that is secreted in extracellular vesicles (EVs) released by the gastrointestinal nematode Heligmosomoides bakeri, at multiple copies per EV. Phylogenetic and gene expression analyses demonstrate exWAGO orthologues are highly conserved and abundantly expressed in related parasites but highly diverged in free-living genus Caenorhabditis. We show that the most abundant small RNAs released from the nematode parasite are not microRNAs as previously thought, but rather secondary small interfering RNAs (siRNAs) that are produced by RNA-dependent RNA Polymerases. The siRNAs that are released in EVs have distinct evolutionary properties compared to those resident in free-living or parasitic nematodes. Immunoprecipitation of exWAGO demonstrates that it specifically associates with siRNAs from transposons and newly evolved repetitive elements that are packaged in EVs and released into the host environment. Together this work demonstrates molecular and evolutionary selectivity in the small RNA sequences that are released in EVs into the host environment and identifies a novel Argonaute protein as the mediator of this.

Comparative genomics of the tardigrades Hypsibius dujardini and Ramazzottius varieornatus.

Tardigrada, a phylum of meiofaunal organisms, have been at the center of discussions of the evolution of Metazoa, the biology of survival in extreme environments, and the role of horizontal gene transfer in animal evolution. Tardigrada are placed as sisters to Arthropoda and Onychophora (velvet worms) in the superphylum Panarthropoda by morphological analyses, but many molecular phylogenies fail to recover this relationship. This tension between molecular and morphological understanding may be very revealing of the mode and patterns of evolution of major groups. Limnoterrestrial tardigrades display extreme cryptobiotic abilities, including anhydrobiosis and cryobiosis, as do bdelloid rotifers, nematodes, and other animals of the water film. These extremophile behaviors challenge understanding of normal, aqueous physiology: how does a multicellular organism avoid lethal cellular collapse in the absence of liquid water? Meiofaunal species have been reported to have elevated levels of horizontal gene transfer (HGT) events, but how important this is in evolution, and particularly in the evolution of extremophile physiology, is unclear. To address these questions, we resequenced and reassembled the genome of H. dujardini, a limnoterrestrial tardigrade that can undergo anhydrobiosis only after extensive pre-exposure to drying conditions, and compared it to the genome of R. varieornatus, a related species with tolerance to rapid desiccation. The 2 species had contrasting gene expression responses to anhydrobiosis, with major transcriptional change in H. dujardini but limited regulation in R. varieornatus. We identified few horizontally transferred genes, but some of these were shown to be involved in entry into anhydrobiosis. Whole-genome molecular phylogenies supported a Tardigrada+Nematoda relationship over Tardigrada+Arthropoda, but rare genomic changes tended to support Tardigrada+Arthropoda.

No evidence for extensive horizontal gene transfer in the genome of the tardigrade Hypsibius dujardini.

Tardigrades are meiofaunal ecdysozoans that are key to understanding the origins of Arthropoda. Many species of Tardigrada can survive extreme conditions through cryptobiosis. In a recent paper [Boothby TC, et al. (2015) Proc Natl Acad Sci USA 112(52):15976-15981], the authors concluded that the tardigrade Hypsibius dujardini had an unprecedented proportion (17%) of genes originating through functional horizontal gene transfer (fHGT) and speculated that fHGT was likely formative in the evolution of cryptobiosis. We independently sequenced the genome of H. dujardini As expected from whole-organism DNA sampling, our raw data contained reads from nontarget genomes. Filtering using metagenomics approaches generated a draft H. dujardini genome assembly of 135 Mb with superior assembly metrics to the previously published assembly. Additional microbial contamination likely remains. We found no support for extensive fHGT. Among 23,021 gene predictions we identified 0.2% strong candidates for fHGT from bacteria and 0.2% strong candidates for fHGT from nonmetazoan eukaryotes. Cross-comparison of assemblies showed that the overwhelming majority of HGT candidates in the Boothby et al. genome derived from contaminants. We conclude that fHGT into H. dujardini accounts for at most 1-2% of genes and that the proposal that one-sixth of tardigrade genes originate from functional HGT events is an artifact of undetected contamination.

Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana.

Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.

The genome sequence of the Chalkhill Blue, Lysandra coridon (Poda, 1761).

We present a genome assembly from an individual male Lysandra coridon (the Chalkhill Blue; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequence is 541 megabases in span. Most of the assembly is scaffolded into 90 chromosomal pseudomolecules, including the assembled Z sex chromosome. The mitochondrial genome has also been assembled and is 15.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,334 protein coding genes.

The genome sequence of the Mazarine Blue, Cyaniris semiargus (Rottemburg, 1775).

We present a genome assembly from an individual male Cyaniris semiargus (the Mazarine Blue; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequence is 441.5 megabases in span. Most of the assembly is scaffolded into 24 chromosomal pseudomolecules, including the assembled Z sex chromosome. The mitochondrial genome has also been assembled and is 15.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 16,408 protein coding genes.

The genome sequence of the Brown Argus, Aricia agestis (Denis & Schiffermüller, 1775).

We present genome assemblies from two male Aricia agestis specimens (the Brown Argus; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequences are 435.3 and 437.4 megabases in span. Each assembly is scaffolded into 23 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genomes were assembled and are 15.47 and 15.45 kilobases in length. Gene annotation of these assemblies on Ensembl identified 12,688 and 12,654 protein coding genes.

The phylogenetics of Anguillicolidae (Nematoda: Anguillicoloidea), swimbladder parasites of eels.

BACKGROUND: Anguillicolidae Yamaguti, 1935 is a family of parasitic nematode infecting fresh-water eels of the genus Anguilla, comprising five species in the genera Anguillicola and Anguillicoloides. Anguillicoloides crassus is of particular importance, as it has recently spread from its endemic range in the Eastern Pacific to Europe and North America, where it poses a significant threat to new, naïve hosts such as the economic important eel species Anguilla anguilla and Anguilla rostrata. The Anguillicolidae are therefore all potentially invasive taxa, but the relationships of the described species remain unclear. Anguillicolidae is part of Spirurina, a diverse clade made up of only animal parasites, but placement of the family within Spirurina is based on limited data. RESULTS: We generated an extensive DNA sequence dataset from three loci (the 5' one-third of the nuclear small subunit ribosomal RNA, the D2-D3 region of the nuclear large subunit ribosomal RNA and the 5' half of the mitochondrial cytochrome c oxidase I gene) for the five species of Anguillicolidae and used this to investigate specific and generic boundaries within the family, and the relationship of Anguillicolidae to other spirurine nematodes. Neither nuclear nor mitochondrial sequences supported monophyly of Anguillicoloides. Genetic diversity within the African species Anguillicoloides papernai was suggestive of cryptic taxa, as was the finding of distinct lineages of Anguillicoloides novaezelandiae in New Zealand and Tasmania. Phylogenetic analysis of the Spirurina grouped the Anguillicolidae together with members of the Gnathostomatidae and Seuratidae. CONCLUSIONS: The Anguillicolidae is part of a complex radiation of parasitic nematodes of vertebrates with wide host diversity (chondrichthyes, teleosts, squamates and mammals), most closely related to other marine vertebrate parasites that also have complex life cycles. Molecular analyses do not support the recent division of Anguillicolidae into two genera. The described species may hide cryptic taxa, identified here by DNA taxonomy, and this DNA barcoding approach may assist in tracking species invasions. The propensity for host switching, and thus the potential for invasive behaviour, is found in A. crassus, A. novaezelandiae and A. papernai, and thus may be common to the group.

The genome sequence of the high brown fritillary, Fabriciana adippe (Dennis & Schiffermüller, 1775).

We present a genome assembly from an individual female Fabriciana adippe (the high brown fritillary; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 485 megabases in span. Most of the assembly (99.98%) is scaffolded into 29 chromosomal pseudomolecules with the Z sex chromosome assembled. The complete mitochondrial genome was also assembled and is 15.1 kilobases in length. Gene annotation of this assembly in Ensembl identified 13,536 protein coding genes.

The genome sequence of the Arran brown, Erebia ligea (Linnaeus, 1758).

We present a genome assembly from an individual male Erebia ligea (Arran brown; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 506 megabases in span. The majority (99.92%) of the assembly is scaffolded into 29 chromosomal pseudomolecules, with the Z sex chromosome assembled. The complete mitochondrial genome was also assembled and is 15.2 kilobases in length.

The genome sequence of the silver-studded blue, Plebejus argus (Linnaeus, 1758).

We present a genome assembly from an individual male Plebejus argus (silver-studded blue; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequence is 382 megabases in span. The entire assembly (100%) is scaffolded into 23 chromosomal pseudomolecules with the Z sex chromosome assembled. The complete mitochondrial genome was also assembled and is 27.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,693 protein coding genes.

The genome sequence of the white admiral, Limenitis camilla (Linnaeus, 1764).

We present a genome assembly from an individual female Limenitis camilla (the white admiral; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 435 megabases in span. Most of the assembly (99.97%) is scaffolded into 31 chromosomal pseudomolecules, corresponding to 29 autosomes plus the W and Z sex chromosomes. The complete mitochondrial genome was also assembled and is 15.2 kilobases in length. Gene annotation of this assembly on Ensembl identified 12,489 protein coding genes.

The genome sequence of the lesser marbled fritillary, Brenthis ino, and evidence for a segregating neo-Z chromosome.

The lesser marbled fritillary, Brenthis ino (Rottemburg, 1775), is a species of Palearctic butterfly. Male Brenthis ino individuals have been reported to have between 12 and 14 pairs of chromosomes, a much-reduced chromosome number than is typical in butterflies. Here, we present a chromosome-level genome assembly for Brenthis ino, as well as gene and transposable element annotations. The assembly is 411.8 Mb in length with a contig N50 of 9.6 Mb and a scaffold N50 of 29.5 Mb. We also show evidence that the male individual from which we generated HiC data was heterozygous for a neo-Z chromosome, consistent with inheriting 14 chromosomes from one parent and 13 from the other. This genome assembly will be a valuable resource for studying chromosome evolution in Lepidoptera, as well as for comparative and population genomics more generally.

The genome sequence of the scarce swallowtail, Iphiclides podalirius.

The scarce swallowtail, Iphiclides podalirius (Linnaeus, 1758), is a species of butterfly in the family Papilionidae. Here, we present a chromosome-level genome assembly for Iphiclides podalirius as well as gene and transposable element annotations. We investigate how the density of genomic features differs between the 30 Iphiclides podalirius chromosomes. We find that shorter chromosomes have higher heterozygosity at four-fold-degenerate sites and a greater density of transposable elements. While the first result is an expected consequence of differences in recombination rate, the second suggests a counter-intuitive relationship between recombination and transposable element evolution. This high-quality genome assembly, the first for any species in the tribe Leptocircini, will be a valuable resource for population genomics in the genus Iphiclides and comparative genomics more generally.