Luana S. Maroja

Butterfly genome reveals promiscuous exchange of mimicry adaptations among species

The evolutionary importance of hybridization and introgression has long been debated. Hybrids are usually rare and unfit, but even infrequent hybridization can aid adaptation by transferring beneficial traits between species. Here we use genomic tools to investigate introgression in Heliconius, a rapidly radiating genus of neotropical butterflies widely used in studies of ecology, behaviour, mimicry and speciation. We sequenced the genome of Heliconius melpomene and compared it with other taxa to investigate chromosomal evolution in Lepidoptera and gene flow among multiple Heliconius species and races. Among 12,669 predicted genes, biologically important expansions of families of chemosensory and Hox genes are particularly noteworthy. Chromosomal organization has remained broadly conserved since the Cretaceous period, when butterflies split from the Bombyx (silkmoth) lineage. Using genomic resequencing, we show hybrid exchange of genes between three co-mimics, Heliconius melpomene, Heliconius timareta and Heliconius elevatus, especially at two genomic regions that control mimicry pattern. We infer that closely related Heliconius species exchange protective colour-pattern genes promiscuously, implying that hybridization has an important role in adaptive radiation.

Genomic hotspots for adaptation: the population genetics of Müllerian mimicry in the Heliconius melpomene clade

Wing patterning in Heliconius butterflies is a longstanding example of both Müllerian mimicry and phenotypic radiation under strong natural selection. The loci controlling such patterns are “hotspots” for adaptive evolution with great allelic diversity across different species in the genus. We characterise nucleotide variation, genotype-by-phenotype associations, linkage disequilibrium, and candidate gene expression at two loci and across multiple hybrid zones in Heliconius melpomene and relatives. Alleles at HmB control the presence or absence of the red forewing band, while alleles at HmYb control the yellow hindwing bar. Across HmYb two regions, separated by ∼100 kb, show significant genotype-by-phenotype associations that are replicated across independent hybrid zones. In contrast, at HmB a single peak of association indicates the likely position of functional sites at three genes, encoding a kinesin, a G-protein coupled receptor, and an mRNA splicing factor. At both HmYb and HmB there is evidence for enhanced linkage disequilibrium (LD) between associated sites separated by up to 14 kb, suggesting that multiple sites are under selection. However, there was no evidence for reduced variation or deviations from neutrality that might indicate a recent selective sweep, consistent with these alleles being relatively old. Of the three genes showing an association with the HmB locus, the kinesin shows differences in wing disc expression between races that are replicated in the co-mimic, Heliconius erato, providing striking evidence for parallel changes in gene expression between Müllerian co-mimics. Wing patterning loci in Heliconius melpomene therefore show a haplotype structure maintained by selection, but no evidence for a recent selective sweep. The complex genetic pattern contrasts with the simple genetic basis of many adaptive traits studied previously, but may provide a better model for most adaptation in natural populations that has arisen over millions rather than tens of years.

Searching for candidate speciation genes using a proteomic approach: seminal proteins in field crickets

In many animals, male seminal proteins influence gamete interactions and fertilization ability and are probably involved in barriers to gene flow between diverging lineages. Here we use a proteomic approach to identify seminal proteins that are transferred to females during copulation and that may be involved in fertilization barriers between two hybridizing field crickets (Gryllus firmus and Gryllus pennsylvanicus). Analyses of patterns of divergence suggest that much of the field cricket genome has remained undifferentiated following the evolution of reproductive isolation. By contrast, seminal protein genes are highly differentiated. Tests of selection reveal that positive selection is likely to be responsible for patterns of differentiation. Together, our observations suggest that some of the loci encoding seminal proteins may indeed play a role in fertilization barriers in field crickets.

Phylogeography of spruce beetles (Dendroctonus rufipennis Kirby) (Curculionidae: Scolytinae) in North America

Tree‐feeding insects that are widespread in north temperate regions are excellent models for studying how past glaciations have impacted differentiation and speciation. We used mitochondrial DNA (mtDNA) sequences and allele frequencies at nine microsatellite loci to examine genetic population structure across the current range of the spruce beetle (Dendroctonus rufipennis), an economically important insect in North America. Two major haplotype groups occur across northern North America, from Newfoundland to Alaska, on white spruce (Picea glauca), and a third distinctive haplotype group occurs throughout the Rocky Mountains on Engelmann spruce (Picea engelmannii). The two mtDNA lineages found in northern populations are 3–4% divergent from each other and from the lineages found in the Rocky Mountains. Analyses of microsatellite data also suggest the existence of major population groupings associated with different geographical regions. In the Pacific Northwest, concordant contact zones for genetically distinct populations of spruce beetles and their principal hosts appear to reflect recent secondary contact. Although we could detect no evidence of historical mtDNA gene flow between allopatric population groups, patterns of variation in the Pacific Northwest suggest recent hybridization and introgression. Together with the pollen record for spruce, they also suggest that beetles have spread from at least three glacial refugia. A minimum estimate of divergence time between the Rocky Mountain and northern populations was 1.7 Myr (million years), presumably reflecting the combined effects of isolation during multiple glacial cycles.

GENEALOGICAL DISCORDANCE AND PATTERNS OF INTROGRESSION AND SELECTION ACROSS A CRICKET HYBRID ZONE

In recently diverged species, ancestral polymorphism and introgression can cause incongruence between gene and species trees. In the face of hybridization, few genomic regions may exhibit reciprocal monophyly, and these regions, usually evolving rapidly under selection, may be important for the maintenance of species boundaries. In animals with internal fertilization, genes encoding seminal protein are candidate barrier genes. Recently diverged hybridizing species such as the field crickets Gryllus firmus and G. pennsylvanicus, offer excellent opportunities to investigate the origins of barriers to gene exchange. These recently diverged species form a well-characterized hybrid zone, and share ancestral polymorphisms across the genome. We analyzed DNA sequence divergence for seminal protein loci, housekeeping loci, and mtDNA, using a combination of analytical approaches and extensive sampling across both species and the hybrid zone. We report discordant genealogical patterns and differential introgression rates across the genome. The most dramatic outliers, showing near-zero introgression and more structured species trees, are also the only two seminal protein loci under selection. These are candidate barrier genes with possible reproductive functions. We also use genealogical data to examine the demographic history of the field crickets and the current structure of the hybrid zone.

The gene cortex controls mimicry and crypsis in butterflies and moths

The wing patterns of butterflies and moths (Lepidoptera) are diverse and striking examples of evolutionary diversification by natural selection1,2. Lepidopteran wing colour patterns are a key innovation, consisting of arrays of coloured scales. We still lack a general understanding of how these patterns are controlled and if there is any commonality across the 160,000 moth and 17,000 butterfly species. Here, we identify a gene, cortex, through fine-scale mapping using population genomics and gene expression analyses, which regulates pattern switches in multiple species across the mimetic radiation in Heliconius butterflies. cortex belongs to a fast evolving subfamily of the otherwise highly conserved fizzy family of cell cycle regulators3, suggesting that it most likely regulates pigmentation patterning through regulation of scale cell development. In parallel with findings in the peppered moth (Biston betularia)4, our results suggest that this mechanism is common within Lepidoptera and that cortex has become a major target for natural selection acting on colour and pattern variation in this group of insects.

Major Improvements to the Heliconius melpomene Genome Assembly Used to Confirm 10 Chromosome Fusion Events in 6 Million Years of Butterfly Evolution

The Heliconius butterflies are a widely studied adaptive radiation of 46 species spread across Central and South America, several of which are known to hybridize in the wild. Here, we present a substantially improved assembly of the Heliconius melpomene genome, developed using novel methods that should be applicable to improving other genome assemblies produced using short read sequencing. First, we whole-genome-sequenced a pedigree to produce a linkage map incorporating 99% of the genome. Second, we incorporated haplotype scaffolds extensively to produce a more complete haploid version of the draft genome. Third, we incorporated ∼20x coverage of Pacific Biosciences sequencing, and scaffolded the haploid genome using an assembly of this long-read sequence. These improvements result in a genome of 795 scaffolds, 275 Mb in length, with an N50 length of 2.1 Mb, an N50 number of 34, and with 99% of the genome placed, and 84% anchored on chromosomes. We use the new genome assembly to confirm that the Heliconius genome underwent 10 chromosome fusions since the split with its sister genus Eueides, over a period of about 6 million yr.

Convergent, modular expression of ebony and tan in the mimetic wing patterns of Heliconius butterflies

The evolution of pigmentation in vertebrates and flies has involved repeated divergence at a small number of genes related to melanin synthesis. Here, we study insect melanin synthesis genes in Heliconius butterflies, a group characterised by its diversity of wing patterns consisting of black (melanin), and yellow and red (ommochrome) pigmented scales. Consistent with their respective biochemical roles in Drosophila melanogaster, ebony is upregulated in non-melanic wing regions destined to be pigmented red whilst tan is upregulated in melanic regions. Wing regions destined to be pigmented yellow, however, are downregulated for both genes. This pattern is conserved across multiple divergent and convergent phenotypes within the Heliconii, suggesting a conserved mechanism for the development of black, red and yellow pattern elements across the genus. Linkage mapping of five melanin biosynthesis genes showed that, in contrast to other organisms, these genes do not control pattern polymorphism. Thus, the pigmentation genes themselves are not the locus of evolutionary change but lie downstream of a wing pattern regulatory factor. The results suggest a modular system in which particular combinations of genes are switched on whenever red, yellow or black pattern elements are favoured by natural selection for diverse and mimetic wing patterns.

Small mammal populations of an agroecosystem in the Atlantic Forest domain, southeastern Brazil.

This study reports 2 years of the population dynamics and reproduction of a small mammal community using the removal method. The study was conducted in a rural area of the Atlantic Forest, in Sumidouro, Rio de Janeiro State, Brazil. The population sizes, age structure and reproduction were studied for the four most common species in the study area. The overall diversity was 1.67 and ranged between 0.8 to 1.67. The species richness was 13 considering the whole study. The most abundant species were the rodents Nectomys squamipes (n = 133), Akodon cursor (n = 74), Oligoryzomys nigripes (n = 25) and the marsupials Didelphis aurita (n = 58) and Philander frenatus (n = 50). Seven other rodents were captured once: Necromys lasiurus, Akodon montensis, Sooretamys angouya, Oecomys catherine, Oxymycterus judex, Euryzygomatomys spinosus and Trinomys iheringi. There were higher peaks for diversity and species richness during the winter (dry) months, probably due to higher food availability. The marsupials had a seasonal reproduction with highest population sizes at the end of the rainy seasons. Nectomys squamipes reproduced mostly during rainy periods. Akodon cursor reproduced predominantly in the winter with the highest population peaks occurring during this season. The analysis of the population dynamics of the rodent species indicated that no species behaved as an agricultural pest, probably due to the heterogeneous landscape of high rotativity of vegetable cultivation. Rodent populations were more susceptible to the removal procedure than marsupial ones.

Genes with Restricted Introgression in a Field Cricket (Gryllus firmus/Gryllus pennsylvanicus) Hybrid Zone Are Concentrated on the X Chromosome and a Single Autosome

Characterizing the extent of genomic differentiation between recently diverged lineages provides an important context for understanding the early stages of speciation. When such lineages form discrete hybrid zones, patterns of differential introgression allow direct estimates of which genome regions are likely involved in speciation and local adaptation. Here we use a backcross experimental design to construct a genetic linkage map for the field crickets Gryllus firmus and Gryllus pennsylvanicus, which interact in a well-characterized hybrid zone in eastern North America. We demonstrate that loci with major allele frequency differences between allopatric populations are not randomly distributed across the genome. Instead, most are either X-linked or map to a few small autosomal regions. Furthermore, the subset of those highly differentiated markers that exhibit restricted introgression across the cricket hybrid zone are also concentrated on the X chromosome (39 of 50 loci) and in a single 7-cM region of one autosome. Although the accumulation on the sex chromosome of genes responsible for postzygotic barriers is a well-known phenomenon, less attention has been given to the genomic distribution of genes responsible for prezygotic barriers. We discuss the implications of our results for speciation, both in the context of the role of sex chromosomes and also with respect to the likely causes of heterogeneous genomic divergence. Although we do not yet have direct evidence for the accumulation of ecological, behavioral, or fertilization prezygotic barrier genes on the X chromosome, faster-X evolution could make these barriers more likely to be X-linked.