Contrasting genomic and phenotypic outcomes of hybridization between pairs of mimetic butterfly taxa across a suture zone

Hybrid zones, whereby divergent lineages come into contact and eventually hybridize, can provide insights on the mechanisms involved in population differentiation and reproductive isolation, and ultimately speciation. Suture zones offer the opportunity to compare these processes across multiple species. In this paper we use reduced‐complexity genomic data to compare the genetic and phenotypic structure and hybridization patterns of two mimetic butterfly species, Ithomia salapia and Oleria onega (Nymphalidae: Ithomiini), each consisting of a pair of lineages differentiated for their wing colour pattern and that come into contact in the Andean foothills of Peru. Despite similarities in their life history, we highlight major differences, both at the genomic and phenotypic level, between the two species. These differences include the presence of hybrids, variations in wing phenotype, and genomic patterns of introgression and differentiation. In I. salapia, the two lineages appear to hybridize only rarely, whereas in O. onega the hybrids are not only more common, but also genetically and phenotypically more variable. We also detected loci statistically associated with wing colour pattern variation, but in both species these loci were not over‐represented among the candidate barrier loci, suggesting that traits other than wing colour pattern may be important for reproductive isolation. Our results contrast with the genomic patterns observed between hybridizing lineages in the mimetic Heliconius butterflies, and call for a broader investigation into the genomics of speciation in Ithomiini ‐ the largest radiation of mimetic butterflies.     

Geographical patterns of evolution in Neotropical Lepidoptera: differentiation of the species of Melinaea and Mechanitis (Nymphalidae,Ithomiinae)*

Abstract Biosystematic analysis incorporating abundant new field data from many parts of the Neotropics has led to an ordered revision of the mimetic ithomiine genera Melinaea and Mechanitis. The various polytypic species of these general probably served as prime movers for the differentiation of other mimetic butterflies in Quaternary forest refuges. The revisions are presented in the form of supplements to the works of Richard M. Fox on these genera, with analyses based on his divisions. Seven species (or monophyletic species-groups) and sixty-three well-differentiated geographic subspecies (six of these described here for the first time) are recognized in Melinaea. Specimens are illustrated which demonstrate intergradation between refuge-derived subspecies. Five species and fifty-two differentiated subspecies are recognized in Mechanitis, whose members are more abundant and gregarious, more plastic, and apparently more vagile than those of Melinaea, resulting in fewer clear-cut mimetic associations, more extensive blurring of differentiation patterns, and apparently fewer incipient biological species in this genus than in Melinaea.     

Two sisters in the same dress: <Emphasis Type="Italic">Heliconius</Emphasis> cryptic species

Background  

Sister species divergence and reproductive isolation commonly results from ecological adaptation. In mimetic Heliconius butterflies, shifts in colour pattern contribute to pre- and post-mating reproductive isolation and are commonly correlated with speciation. Closely related mimetic species are therefore not expected, as they should lack several important sources of reproductive isolation.     

Chromosome-level de novo genome assembly of two conifer-parasitic wasps,Megastigmus duclouxiana and Megastigmus sabinae,reveals genomic imprints of adaptation to hosts

Conifers make up about one third of global forests but are threatened by seed parasitoid wasp species. Many of these wasps belong to the genus Megastigmus, yet little is known about their genomic background. In this study, we provide chromosome-level genome assemblies for two oligophagous conifer parasitoid species of Megastigmus, which represent the first two chromosome-level genomes of the genus. The assembled genomes of Megastigmus duclouxiana and M. sabinae are 878.48 Mb (scaffold N50 of 215.60 Mb) and 812.98 Mb (scaffold N50 of 139.16 Mb), respectively, which are larger than the genome size of most hymenopterans due to the expansion of transposable elements. Expanded gene families highlight the difference in sensory-related genes between the two species, reflecting the difference in their hosts. We further found that these two species have fewer family members but more single-gene duplications than polyphagous congeners in the gene families of ATP-binding cassette transporter (ABC), cytochrome P450 (P450) and olfactory receptors (OR). These findings shed light on the pattern of adaptation to a narrow spectrum of hosts in oligophagous parasitoids. Our findings suggest potential drivers underlying genome evolution and parasitism adaptation, and provide valuable resources for understanding the ecology, genetics and evolution of Megastigmus, as well as for the research and biological control of global conifer forest pests.     

Assembly and RNA‐free annotation of highly heterozygous genomes: The case of the thick‐billed murre (Uria lomvia)

Thanks to a dramatic reduction in sequencing costs followed by a rapid development of bioinformatics tools, genome assembly and annotation have become accessible to many researchers in recent years. Among tetrapods, birds have genomes that display many features that facilitate their assembly and annotation, such as small genome size, low number of repeats and highly conserved genomic structure. However, we found that high genomic heterozygosity could have a great impact on the quality of the genome assembly of the thick‐billed murre (Uria lomvia), an arctic colonial seabird. In this study, we tested the performance of three genome assemblers, ray /sscape , soapdenovo 2 and platanus , in assembling the highly heterozygous genome of the thick‐billed murre. Our results show that platanus , an assembler specifically designed for heterozygous genomes, outperforms the other two approaches and produces a highly contiguous (N50 = 15.8 Mb) and complete genome assembly (93% presence of genes from the Benchmarking Universal Single Copy Ortholog [BUSCO] gene set). Additionally, we annotated the thick‐billed murre genome using a homology‐based approach that takes advantage of the genomic resources available for birds and other taxa. Our study will be useful for those researchers who are approaching assembly and annotation of highly heterozygous genomes, or genomes of species of conservation concern, and/or who have limited financial resources.     

The diversification of Heliconius butterflies: what have we learned in 150 years?

Research into Heliconius butterflies has made a significant contribution to evolutionary biology. Here, we review our understanding of the diversification of these butterflies, covering recent advances and a vast foundation of earlier work. Whereas no single group of organisms can be sufficient for understanding life's diversity, after years of intensive study, research into Heliconius has addressed a wide variety of evolutionary questions. We first discuss evidence for widespread gene flow between Heliconius species and what this reveals about the nature of species. We then address the evolution and diversity of warning patterns, both as the target of selection and with respect to their underlying genetic basis. The identification of major genes involved in mimetic shifts, and homology at these loci between distantly related taxa, has revealed a surprising predictability in the genetic basis of evolution. In the final sections, we consider the evolution of warning patterns, and Heliconius diversity more generally, within a broader context of ecological and sexual selection. We consider how different traits and modes of selection can interact and influence the evolution of reproductive isolation.     

Pollen flow in Psiguria warscewiczii: a comparison of Heliconius butterflies and hummingbirds

Summary Pollen dispersal by male and female Heliconius butterflies and hermit hummingbirds was examined in a natural population of the tropical vine Psiguria warscewiczii. Fluorescent dyes were used to track flowers visited by individual pollinators. Mean distance of pollen dispersal was significantly higher by males than by female Heliconius. Butterflies dispersed pollen greater distances than did hummingbirds. There were no significant differences in dispersal distance among species of Heliconius. The number of flowers or plants to which dye was dispersed was greater for Heliconius than for hummingbirds.The home range behavior and site specificity exhibited by Heliconius, together with the preference by most Heliconius species in the study site for Psiguria led us to examine the relationship between the butterflies' home range and its daily foraging among Psiguria. A Heliconius mark-recapture program was conducted simultaneously with the dye program. This study showed that individual butterflies, for which extensive recaptures were made, visited an area of Psiguria plants on any one day that corresponded to the area in which those butterflies were caught or observed. The area over which a butterfly moved was termed its home range; within home ranges relative movement by butterflies was estimated as the mean distance between sequential captures. This estimate differed significantly among species of Heliconius; males of all species moving greater distances than females.     

A fake Heliconius (Lepidoptera, Nymphalidae)

Four ‘doubtful species’ are listed in catalogues of the genus Heliconius. All are non–existent: (1) H. Carolina (Herbst) is a forgery made from trimmed fragments of at least three butterflies, the two identifiable ones both beingH. erato; (2) H.cinereo–fuscus (Goeze) is a faded specimen of H.ricini; (3) H.arcuatus (Goeze) is so inadequately described that its identity cannot be determined, although it is apparently some species of heliconiine or ithomiine;(4)H.euclea (Godart) is the valid senior name of an ithomiid in the genus Hypothyris.     

Variable Selection and the Coexistence of Multiple mimetic forms of the Butterfly Heliconius numata

Polymorphism in aposematic animals and coexistence of multiple mimicry rings within a habitat are not predicted by classical Müllerian mimicry. The butterfly Heliconius numata Cramer (Lepidoptera: Nymphalidae; Heliconiinae) is both polymorphic and aposematic. The polymorphism is due to variation at a single locus (or `supergene') which determines colour patterns involved in Müllerian mimicry. We sampled 11 sites in a small area (approx. 60×30km) of North-eastern Peru for H. numata and its co-mimics in the genus Melinaea and Athyrtis (Ithomiinae), and examined the role of temporal and spatial heterogeneity in the maintenance of polymorphism. Colour-patterns of Melinaea communities, which constitute the likely `mimetic environment' for H. numata, are differentiated on a more local scale than morphs of H. numata, but the latter do show a strong and significant response to local selection for colour-pattern. In contrast, analysis of enzyme polymorphism in H. numata across the region revealed no spatial structure, which is consistent with a high mobility of this species. Differences in spatial variability in the two taxa may have caused H. numata to become polymorphic, while temporal variability, not significant in this study, probably has a lesser effect. The mimetic polymorphism is therefore explained by means of multiple selection-migration clines at a single locus, a similar process to that which explains narrow hybrid zones between geographic races of other Heliconius butterflies. This revised version was published online in July 2006 with corrections to the Cover Date.     

Ithomiini Butterflies (Lepidoptera: Nymphalidae) of Antioquia, Colombia

Colombia is one of the most biodiverse countries on the planet. However, economic and scientific investment in completing inventories of its biodiversity has been relatively poor in comparison with other Neotropical countries. Butterflies are the best studied group of invertebrates, with the highest proportion of known to expected species. More than 3,200 species of butterflies have been recorded in Colombia, although the study of the still many unexplored areas will presumably increase this number. This work provides a list of Ithomiini butterflies collected in the department of Antioquia and estimates the total number of species present, based on revision of entomological collections, records in the literature and field work performed between 2003 and 2011. The list includes 99 species and 32 genera, representing 27% of all Ithomiini species. We report 50 species of Ithomiini not formerly listed from Antioquia, and found the highest diversity of ithomiine species to be at middle elevations (900–1,800 m). The mean value of the Chao2 estimator for number of species in Antioquia is 115 species, which is close to a predicted total of 109 based on known distributions of other Ithomiini not yet recorded from the department. Nine species are potentially of particular conservation importance because of their restricted distributions, and we present range maps for each species. We also highlight areas in Antioquia with a lack of biodiversity knowledge to be targeted in future studies. This paper contributes to mapping the distribution of the Lepidoptera of Antioquia department in particular and of Colombia in general.     

The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft‐ and hard‐seeded cultivars

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single‐molecule sequencing and high‐throughput chromosome conformation capture techniques to produce a high‐quality and long‐range contiguity chromosome‐scale genome assembly of the soft‐seeded pomegranate cultivar ‘Tunisia’. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein‐coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft‐ and hard‐seeded pomegranate varieties. A set of selective loci containing SUC8‐like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard‐ and soft‐seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft‐ and hard‐seeded pomegranate varieties.     

Insights into identifying resistance genes for cold and disease stresses through chromosome-level reference genome analyses of Poncirus polyandra

Poncirus polyandra, a plant species with extremely small populations in China, has become extinct in the wild. This study aimed to identify functional genes that improve tolerance to abiotic and biotic stresses. Here, we present a high-quality chromosome-scale reference genome of P. polyandra. The reference genome is 315.78 Mb in size, with an N50 scaffold size of 32.07 Mb, and contains nine chromosomes with 20,815 protein-coding genes, covering 97.82% of the estimated gene space. We identified 17 rapidly evolving nucleotide-binding-site (NBS) genes, three C-repeat-binding factors (CBF) genes, 19 citrus greening disease (Huanglongbing, HLB) tolerance genes, 11 citrus tristeza virus (CTV) genes, and one citrus nematode resistance gene. A divergence time of 1.96 million years ago was estimated between P. polyandra and P. trifoliata. This is the first genome-scale assembly and annotation of P. polyandra, which will be useful for genetic, genomic, and molecular research and provide guidance for the development of conservation strategies.     

nGASP – the nematode genome annotation assessment project

Background  

While the C. elegans genome is extensively annotated, relatively little information is available for other Caenorhabditis species. The nematode genome annotation assessment project (nGASP) was launched to objectively assess the accuracy of protein-coding gene prediction software in C. elegans, and to apply this knowledge to the annotation of the genomes of four additional Caenorhabditis species and other nematodes. Seventeen groups worldwide participated in nGASP, and submitted 47 prediction sets across 10 Mb of the C. elegans genome. Predictions were compared to reference gene sets consisting of confirmed or manually curated gene models from WormBase.     

A set of microsatellite markers for Heliconius melpomene and closely related species

The butterflies in the genus Heliconius offer an exceptional opportunity for the study of the ecology and genetics of an adaptive radiation due to their extensive intra‐ and interspecific variation in wing colour patterns and mimetic associations. Here, we characterize 22 polymorphic microsatellite loci in Heliconius melpomene that have been shown to be useful for linkage mapping and population studies in this and other species. Levels of variation were high, although heterozygosity deficiencies were found in most loci, probably due to null alleles. The loci showed broad amplification success on six other species across the genus.     

Discrimination of flying mimetic,passion-vine butterflies Heliconius

Wing-beat frequency and the degree of asymmetry in wing motion were more similar among mimics than among sister species of passion-vine butterflies in the genus Heliconius. Asymmetry in wing motion is not attributed to lift production, and serves as the first clear example of a mimetic behavioural signal for a flying organism. Because the similarities in wing motion are too subtle for humans to observe with the naked eye, they serve as a previously unexplored mimetic signal.     

Comparative genomics of Chlamydomonas

Despite its role as a reference organism in the plant sciences, the green alga Chlamydomonas reinhardtii entirely lacks genomic resources from closely related species. We present highly contiguous and well-annotated genome assemblies for three unicellular C. reinhardtii relatives: Chlamydomonas incerta, Chlamydomonas schloesseri, and the more distantly related Edaphochlamys debaryana. The three Chlamydomonas genomes are highly syntenous with similar gene contents, although the 129.2 Mb C. incerta and 130.2 Mb C. schloesseri assemblies are more repeat-rich than the 111.1 Mb C. reinhardtii genome. We identify the major centromeric repeat in C. reinhardtii as a LINE transposable element homologous to Zepp (the centromeric repeat in Coccomyxa subellipsoidea) and infer that centromere locations and structure are likely conserved in C. incerta and C. schloesseri. We report extensive rearrangements, but limited gene turnover, between the minus mating type loci of these Chlamydomonas species. We produce an eight-species core-Reinhardtinia whole-genome alignment, which we use to identify several hundred false positive and missing genes in the C. reinhardtii annotation and >260,000 evolutionarily conserved elements in the C. reinhardtii genome. In summary, these resources will enable comparative genomics analyses for C. reinhardtii, significantly extending the analytical toolkit for this emerging model system.

High-quality genome assemblies and annotations for three of the closest relatives of Chlamydomonas reinhardtii enable comparative genomics analyses.     

Genome of an iconic Australian bird: High‐quality assembly and linkage map of the superb fairy‐wren (Malurus cyaneus)

The superb fairy‐wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high‐quality reference. To rectify that, we present the first high‐quality genome assembly of a Meliphagides representative: the superb fairy‐wren. We combined Illumina shotgun and mate‐pair sequences, PacBio long‐reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07‐Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high‐quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree‐based bird genetic maps, we find that the fairy‐wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high‐quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy‐wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.