Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus.

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Optix and cortex/ivory/mir-193 again: the repeated use of two mimicry hotspot loci.

The extent to which evolution is repeatable has been a debated topic among evolutionary biologists. Although rewinding the tape of life perhaps would not lead to the same outcome every time, repeated evolution of analogous genes for similar functions has been extensively reported. Wing phenotypes of butterflies and moths have provided a wealth of examples of gene re-use, with certain 'hotspot loci' controlling wing patterns across diverse taxa. Here, we present an example of convergent evolution in the molecular genetic basis of Batesian wing mimicry in two Hypolimnas butterfly species. We show that mimicry is controlled by variation near cortex/ivory/mir-193, a known butterfly hotspot locus. By dissecting the genetic architecture of mimicry in Hypolimnas misippus and Hypolimnas bolina, we present evidence that distinct non-coding regions control the development of white pattern elements in the forewing and hindwing of the two species, suggesting independent evolution, and that no structural variation is found at the locus. Finally, we also show that orange coloration in H. bolina is associated with optix, a well-known patterning gene. Overall, our study once again implicates variation near the hotspot loci cortex/ivory/mir-193 and optix in butterfly wing mimicry and thereby highlights the repeatability of adaptive evolution.

A novel terpene synthase controls differences in anti-aphrodisiac pheromone production between closely related Heliconius butterflies.

Plants and insects often use the same compounds for chemical communication, but not much is known about the genetics of convergent evolution of chemical signals. The terpene (E)-ß-ocimene is a common component of floral scent and is also used by the butterfly Heliconius melpomene as an anti-aphrodisiac pheromone. While the biosynthesis of terpenes has been described in plants and microorganisms, few terpene synthases (TPSs) have been identified in insects. Here, we study the recent divergence of 2 species, H. melpomene and Heliconius cydno, which differ in the presence of (E)-ß-ocimene; combining linkage mapping, gene expression, and functional analyses, we identify 2 novel TPSs. Furthermore, we demonstrate that one, HmelOS, is able to synthesise (E)-ß-ocimene in vitro. We find no evidence for TPS activity in HcydOS (HmelOS ortholog of H. cydno), suggesting that the loss of (E)-ß-ocimene in this species is the result of coding, not regulatory, differences. The TPS enzymes we discovered are unrelated to previously described plant and insect TPSs, demonstrating that chemical convergence has independent evolutionary origins.

Whole-chromosome hitchhiking driven by a male-killing endosymbiont.

Neo-sex chromosomes are found in many taxa, but the forces driving their emergence and spread are poorly understood. The female-specific neo-W chromosome of the African monarch (or queen) butterfly Danaus chrysippus presents an intriguing case study because it is restricted to a single 'contact zone' population, involves a putative colour patterning supergene, and co-occurs with infection by the male-killing endosymbiont Spiroplasma. We investigated the origin and evolution of this system using whole genome sequencing. We first identify the 'BC supergene', a broad region of suppressed recombination across nearly half a chromosome, which links two colour patterning loci. Association analysis suggests that the genes yellow and arrow in this region control the forewing colour pattern differences between D. chrysippus subspecies. We then show that the same chromosome has recently formed a neo-W that has spread through the contact zone within approximately 2,200 years. We also assembled the genome of the male-killing Spiroplasma, and find that it shows perfect genealogical congruence with the neo-W, suggesting that the neo-W has hitchhiked to high frequency as the male-killer has spread through the population. The complete absence of female crossing-over in the Lepidoptera causes whole-chromosome hitchhiking of a single neo-W haplotype, carrying a single allele of the BC supergene and dragging multiple non-synonymous mutations to high frequency. This has created a population of infected females that all carry the same recessive colour patterning allele, making the phenotypes of each successive generation highly dependent on uninfected male immigrants. Our findings show how hitchhiking can occur between the physically unlinked genomes of host and endosymbiont, with dramatic consequences.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths.

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

End Users Want Alternative Intervention Delivery Models: Usability and Acceptability of the REMOTE-CR Exercise-Based Cardiac Telerehabilitation Program.

OBJECTIVE: Evaluate user experiences of an exercise-based cardiac telerehabilitation intervention (REMOTE-CR) that provided near universal access to real-time remote coaching and behavioral support from exercise specialists. DESIGN: Secondary analysis (12-week follow-up) of a parallel group, single blind, randomized controlled noninferiority trial (ACTRN12614000843651). SETTING: Community-based cardiac rehabilitation. PARTICIPANTS: Adults (N=162) with coronary heart disease who were eligible for outpatient cardiac rehabilitation. Eighty-two of 162 trial participants were randomized to receive REMOTE-CR; 67 completed usability and acceptability assessment at 12-week follow-up. INTERVENTION: REMOTE-CR comprised 12 weeks of individualized exercise prescription, real-time physiological monitoring, coaching, and behavioral support, delivered via a bespoke telerehabilitation platform. OUTCOMES: Ease of use, satisfaction with the technology platform and intervention content, and demand for real-world implementation as an alternative to traditional center-based programs were assessed at 12-week follow-up. RESULTS: Components of usability and acceptability were positively evaluated by most participants (44-66 of 67, 66%-99%). Fifty-eight of 67 (87%) would choose REMOTE-CR if it was available as a usual care service, primarily because it provides convenient and flexible access to real-time individualized support from exercise specialists. Technology challenges were rare and had little effect on user experiences or demand for REMOTE-CR. CONCLUSIONS: REMOTE-CR can extend the reach and impact of existing cardiac rehabilitation services by overcoming traditional participation barriers while preserving expert oversight. Adoption of emerging technologies should be accelerated to support dynamic, engaging, individualized intervention delivery models, but optimizing overall cardiac rehabilitation participation rates will require multiple delivery models that are tailored to satisfy diverse participant preferences.

Identification and functional analyses of sex determination genes in the sexually dimorphic stag beetle Cyclommatus metallifer.

BACKGROUND: Genes in the sex determination pathway are important regulators of sexually dimorphic animal traits, including the elaborate and exaggerated male ornaments and weapons of sexual selection. In this study, we identified and functionally analyzed members of the sex determination gene family in the golden metallic stag beetle Cyclommatus metallifer, which exhibits extreme differences in mandible size between males and females. RESULTS: We constructed a C. metallifer transcriptomic database from larval and prepupal developmental stages and tissues of both males and females. Using Roche 454 pyrosequencing, we generated a de novo assembled database from a total of 1,223,516 raw reads, which resulted in 14,565 isotigs (putative transcript isoforms) contained in 10,794 isogroups (putative identified genes). We queried this database for C. metallifer conserved sex determination genes and identified 14 candidate sex determination pathway genes. We then characterized the roles of several of these genes in development of extreme sexual dimorphic traits in this species. We performed molecular expression analyses with RT-PCR and functional analyses using RNAi on three C. metallifer candidate genes--Sex-lethal (CmSxl), transformer-2 (Cmtra2), and intersex (Cmix). No differences in expression pattern were found between the sexes for any of these three genes. In the RNAi gene-knockdown experiments, we found that only the Cmix had any effect on sexually dimorphic morphology, and these mimicked the effects of Cmdsx knockdown in females. Knockdown of CmSxl had no measurable effects on stag beetle phenotype, while knockdown of Cmtra2 resulted in complete lethality at the prepupal period. These results indicate that the roles of CmSxl and Cmtra2 in the sex determination cascade are likely to have diverged in stag beetles when compared to Drosophila. Our results also suggest that Cmix has a conserved role in this pathway. In addition to those three genes, we also performed a more complete functional analysis of the C. metallifer dsx gene (Cmdsx) to identify the isoforms that regulate dimorphism more fully using exon-specific RNAi. We identified a total of 16 alternative splice variants of the Cmdsx gene that code for up to 14 separate exons. Despite the variation in RNA splice products of the Cmdsx gene, only four protein isoforms are predicted. The results of our exon-specific RNAi indicated that the essential CmDsx isoform for postembryonic male differentiation is CmDsxB, whereas postembryonic female specific differentiation is mainly regulated by CmDsxD. CONCLUSIONS: Taken together, our results highlight the importance of studying the function of highly conserved sex determination pathways in numerous insect species, especially those with dramatic and exaggerated sexual dimorphism, because conservation in protein structure does not always translate into conservation in downstream function.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates.

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

A general mechanism for conditional expression of exaggerated sexually-selected traits.

Sexually-selected exaggerated traits tend to be unusually reliable signals of individual condition, as their expression tends to be more sensitive to nutritional history and physiological circumstance than that of other phenotypes. As such, these traits are the foundation for many models of sexual selection and animal communication, such as "handicap" and "good genes" models. Exactly how expression of these traits is linked to the bearer's condition has been a central yet unresolved question, in part because the underlying physiological mechanisms regulating their development have remained largely unknown. Recent discoveries across animals as diverse as deer, beetles, and flies now implicate the widely conserved insulin-like signaling pathway, as a common physiological mechanism regulating condition-sensitive structures with extreme growth. This raises the exciting possibility that one highly conserved pathway may underlie the evolution of trait exaggeration in a multitude of sexually-selected signal traits across the animal kingdom.

Fine-scale spatial genetic structure of common and declining bumble bees across an agricultural landscape.

Land-use changes have threatened populations of many insect pollinators, including bumble bees. Patterns of dispersal and gene flow are key determinants of species' ability to respond to land-use change, but have been little investigated at a fine scale (<10 km) in bumble bees. Using microsatellite markers, we determined the fine-scale spatial genetic structure of populations of four common Bombus species (B. terrestris, B. lapidarius, B. pascuorum and B. hortorum) and one declining species (B. ruderatus) in an agricultural landscape in Southern England, UK. The study landscape contained sown flower patches representing agri-environment options for pollinators. We found that, as expected, the B. ruderatus population was characterized by relatively low heterozygosity, number of alleles and colony density. Across all species, inbreeding was absent or present but weak (FIS  = 0.01-0.02). Using queen genotypes reconstructed from worker sibships and colony locations estimated from the positions of workers within these sibships, we found that significant isolation by distance was absent in B. lapidarius, B. hortorum and B. ruderatus. In B. terrestris and B. pascuorum, it was present but weak; for example, in these two species, expected relatedness of queens founding colonies 1 m apart was 0.02. These results show that bumble bee populations exhibit low levels of spatial genetic structure at fine spatial scales, most likely because of ongoing gene flow via widespread queen dispersal. In addition, the results demonstrate the potential for agri-environment scheme conservation measures to facilitate fine-scale gene flow by creating a more even distribution of suitable habitats across landscapes.

The remote exercise monitoring trial for exercise-based cardiac rehabilitation (REMOTE-CR): a randomised controlled trial protocol.

BACKGROUND: Exercise is an essential component of contemporary cardiac rehabilitation programs for the secondary prevention of coronary heart disease. Despite the benefits associated with regular exercise, adherence with supervised exercise-based cardiac rehabilitation remains low. Increasingly powerful mobile technologies, such as smartphones and wireless physiological sensors, may extend the capability of exercise-based cardiac rehabilitation by enabling real-time exercise monitoring for those with coronary heart disease. This study compares the effectiveness of technology-assisted, home-based, remote monitored exercise-based cardiac rehabilitation (REMOTE) to standard supervised exercise-based cardiac rehabilitation in New Zealand adults with a diagnosis of coronary heart disease. METHODS/DESIGN: A two-arm, parallel, non-inferiority, randomised controlled trial will be conducted at two sites in New Zealand. One hundred and sixty two participants will be randomised at a 1:1 ratio to receive a 12-week program of technology-assisted, home-based, remote monitored exercise-based cardiac rehabilitation (intervention), or an 8-12 program of standard supervised exercise-based cardiac rehabilitation (control).The primary outcome is post-treatment maximal oxygen uptake (V̇O2max). Secondary outcomes include cardiovascular risk factors (blood lipid and glucose concentrations, blood pressure, anthropometry), self-efficacy, intentions and motivation to be active, objectively measured physical activity, self-reported leisure time exercise and health-related quality of life. Cost information will also be collected to compare the two modes of delivery. All outcomes are assessed at baseline, post-treatment, and 6 months, except for V̇O2max, blood lipid and glucose concentrations, which are assessed at baseline and post-treatment only. DISCUSSION: This novel study will compare the effectiveness of technology-supported exercise-based cardiac rehabilitation to a traditional supervised approach. If the REMOTE program proves to be as effective as traditional cardiac rehabilitation, it has potential to augment current practice by increasing access for those who cannot utilise existing services. TRIAL REGISTRATION: Australian New Zealand Clinical Trials RegistryStudy ID number: ACTRN12614000843651. Registered 7 August 2014.

Cloning and characterization of an mRNA encoding an insulin receptor from the horned scarab beetle Onthophagus nigriventris (Coleoptera: Scarabaeidae).

The insulin signaling pathway is the primary signaling pathway coupling growth with nutritional condition in all animals. Sensitivity to circulating levels of insulin has been shown to regulate the growth of specific traits in a dose-dependent manner in response to environmental conditions in a diversity of insect species. Alternative phenotypes in insects manifest in a variety of morphologies such as the sexually dimorphic and male dimorphic horned beetles. Large males of the sexually dimorphic dung beetle Onthophagus nigriventris develop a thoracic horn up to twice the length of the body whereas small males and females never develop this horn. The regulation of this dimorphism is known to be nutrition dependent for males. We focused on the insulin signaling pathway as a potential regulator of this dimorphism. We sequenced a full-length gene transcript encoding the O. nigriventris insulin receptor (OnInR), which is the receptor for circulating insulin and insulin-like peptides in animals. We show that the predicted OnInR protein is similar in overall amino acid identity to other insulin receptors (InRs) and is most closely related phylogenetically to insect InRs. Expression of the OnInR transcript was found during development of imaginal tissues in both males and females. However, expression of OnInR in the region where a horn would grow of small males and female was significantly higher than in the horn tissues of large males at the end of growth. This variation in OnInR expression between sexes and morphs indicates a role for the InR in polymorphic horn development.

High level of novelty under the hood of convergent evolution.

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

Repeated genetic adaptation to altitude in two tropical butterflies.

Repeated evolution can provide insight into the mechanisms that facilitate adaptation to novel or changing environments. Here we study adaptation to altitude in two tropical butterflies, Heliconius erato and H. melpomene, which have repeatedly and independently adapted to montane habitats on either side of the Andes. We sequenced 518 whole genomes from altitudinal transects and found many regions differentiated between highland (~ 1200 m) and lowland (~ 200 m) populations. We show repeated genetic differentiation across replicate populations within species, including allopatric comparisons. In contrast, there is little molecular parallelism between the two species. By sampling five close relatives, we find that a large proportion of divergent regions identified within species have arisen from standing variation and putative adaptive introgression from high-altitude specialist species. Taken together our study supports a role for both standing genetic variation and gene flow from independently adapted species in promoting parallel local adaptation to the environment.

A high-quality, chromosome-level genome assembly of the Black Soldier Fly (Hermetia illucens L.).

Hermetia illucens L. (Diptera: Stratiomyidae), the Black Soldier Fly (BSF) is an increasingly important species for bioconversion of organic material into animal feed. We generated a high-quality chromosome-scale genome assembly of the BSF using Pacific Bioscience, 10X Genomics linked read and high-throughput chromosome conformation capture sequencing technology. Scaffolding the final assembly with Hi-C data produced a highly contiguous 1.01 Gb genome with 99.75% of scaffolds assembled into pseudochromosomes representing seven chromosomes with 16.01 Mb contig and 180.46 Mb scaffold N50 values. The highly complete genome obtained a Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness of 98.6%. We masked 67.32% of the genome as repetitive sequences and annotated a total of 16,478 protein-coding genes using the BRAKER2 pipeline. We analyzed an established lab population to investigate the genomic variation and architecture of the BSF revealing six autosomes and an X chromosome. Additionally, we estimated the inbreeding coefficient (1.9%) of the lab population by assessing runs of homozygosity. This provided evidence for inbreeding events including long runs of homozygosity on chromosome 5. The release of this novel chromosome-scale BSF genome assembly will provide an improved resource for further genomic studies, functional characterization of genes of interest and genetic modification of this economically important species.

Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies.

The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.

Cortex cis-regulatory switches establish scale colour identity and pattern diversity in Heliconius.

In Heliconius butterflies, wing colour pattern diversity and scale types are controlled by a few genes of large effect that regulate colour pattern switches between morphs and species across a large mimetic radiation. One of these genes, cortex, has been repeatedly associated with colour pattern evolution in butterflies. Here we carried out CRISPR knockouts in multiple Heliconius species and show that cortex is a major determinant of scale cell identity. Chromatin accessibility profiling and introgression scans identified cis-regulatory regions associated with discrete phenotypic switches. CRISPR perturbation of these regions in black hindwing genotypes recreated a yellow bar, revealing their spatially limited activity. In the H. melpomene/timareta lineage, the candidate CRE from yellow-barred phenotype morphs is interrupted by a transposable element, suggesting that cis-regulatory structural variation underlies these mimetic adaptations. Our work shows that cortex functionally controls scale colour fate and that its cis-regulatory regions control a phenotypic switch in a modular and pattern-specific fashion.

Heliconius butterflies have bright patterns on their wings that tell potential predators that they are toxic. As a result, predators learn to avoid eating them. Over time, unrelated species of butterflies have evolved similar patterns to avoid predation through a process known as Müllerian mimicry. Worldwide, there are over 180,000 species of butterflies and moths, most of which have different wing patterns. How do genes create this pattern diversity? And do butterflies use similar genes to create similar wing patterns? One of the genes involved in creating wing patterns is called cortex. This gene has a large region of DNA around it that does not code for proteins, but instead, controls whether cortex is on or off in different parts of the wing. Changes in this non-coding region can act like switches, turning regions of the wing into different colours and creating complex patterns, but it is unclear how these switches have evolved. Butterfly wings get their colour from tiny structures called scales, which each have their own unique set of pigments. In Heliconius butterflies, there are three types of scales: yellow/white scales, black scales, and red/orange/brown scales. Livraghi et al. used a DNA editing technique called CRISPR to find out whether the cortex gene affects scale type. First, Livraghi et al. confirmed that deleting cortex turned black and red scales yellow. Next, they used the same technique to manipulate the non-coding DNA around the cortex gene to see the effect on the wing pattern. This manipulation turned a black-winged butterfly into a butterfly with a yellow wing band, a pattern that occurs naturally in Heliconius butterflies. The next step was to find the mutation responsible for the appearance of yellow wing bands in nature. It turns out that a bit of extra genetic code, derived from so-called 'jumping genes', had inserted itself into the non-coding DNA around the cortex gene, 'flipping' the switch and leading to the appearance of the yellow scales. Genetic information contains the instructions to generate shape and form in most organisms. These instructions evolve over millions of years, creating everything from bacteria to blue whales. Butterfly wings are visual evidence of evolution, but the way their genes create new patterns isn't specific to butterflies. Understanding wing patterns can help researchers to learn how genetic switches control diversity across other species too.

Germline transformation of the stalk-eyed fly, Teleopsis dalmanni.

BACKGROUND: Stalk-eyed flies of the family Diopsidae have proven to be an excellent model organism for studying the evolution of ornamental sexual traits. In diopsid flies the eyes and antennae are borne at the end of lateral head projections called 'eye-stalks'. Eyespan, the distance between the eyes, and the degree of sexual dimorphism in eyespan vary considerably between species and several sexually dimorphic species show sexual selection through female mate preference for males with exaggerated eyespan. Relatively little is known about the molecular genetic basis of intra- or inter-species variation in eyespan, eye-stalk development or growth regulation in diopsids. Molecular approaches including comparative developmental analyses, EST screening and QTL mapping have identified potential candidate loci for eyespan regulation in the model species Teleopsis dalmanni. Functional analyses of these genes to confirm and fully characterise their roles in eye-stalk growth require the development of techniques such as germline transformation to manipulate gene activity in vivo. RESULTS: We used in vivo excision assays to identify transposon vector systems with the activity required to mediate transgenesis in T. dalmanni. Mariner based vectors showed no detectable excision while both Minos and piggyBac were active in stalk-eyed fly embryos. Germline transformation with an overall efficiency of 4% was achieved using a Minos based vector and the 3xP3-EGFP marker construct. Chromosomal insertion of constructs was confirmed by Southern blot analysis. Both autosomal and X-linked inserts were recovered. A homozygous stock, established from one of the X-linked inserts, has maintained stable expression for eight generations. CONCLUSIONS: We have performed stable germline transformation of a stalk-eyed fly, T. dalmanni. This is the first transgenic protocol to be developed in an insect species that exhibits an exaggerated male sexual trait. Transgenesis will enable the development of a range of techniques for analysing gene function in this species and so provide insight into the mechanisms underlying the development of a morphological trait subject to sexual selection. Our X-linked insertion line will permit the sex of live larvae to be determined. This will greatly facilitate the identification of genes which are differentially expressed during eye-stalk development in males and females.

A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning.

BACKGROUND: Diploid genome assembly is typically impeded by heterozygosity because it introduces errors when haplotypes are collapsed into a consensus sequence. Trio binning offers an innovative solution that exploits heterozygosity for assembly. Short, parental reads are used to assign parental origin to long reads from their F1 offspring before assembly, enabling complete haplotype resolution. Trio binning could therefore provide an effective strategy for assembling highly heterozygous genomes, which are traditionally problematic, such as insect genomes. This includes the wood tiger moth (Arctia plantaginis), which is an evolutionary study system for warning colour polymorphism. FINDINGS: We produced a high-quality, haplotype-resolved assembly for Arctia plantaginis through trio binning. We sequenced a same-species family (F1 heterozygosity ∼1.9%) and used parental Illumina reads to bin 99.98% of offspring Pacific Biosciences reads by parental origin, before assembling each haplotype separately and scaffolding with 10X linked reads. Both assemblies are contiguous (mean scaffold N50: 8.2 Mb) and complete (mean BUSCO completeness: 97.3%), with annotations and 31 chromosomes identified through karyotyping. We used the assembly to analyse genome-wide population structure and relationships between 40 wild resequenced individuals from 5 populations across Europe, revealing the Georgian population as the most genetically differentiated with the lowest genetic diversity. CONCLUSIONS: We present the first invertebrate genome to be assembled via trio binning. This assembly is one of the highest quality genomes available for Lepidoptera, supporting trio binning as a potent strategy for assembling heterozygous genomes. Using our assembly, we provide genomic insights into the geographic population structure of A. plantaginis.