Significance
. For millions of years retroviruses, such as HIV in humans, have attacked vertebrates. Occasionally retroviruses infiltrate germ cells, incorporate themselves into the host’s genome, and transmit vertically to the host’s offspring as endogenous retroviruses (ERVs). Consequently, ERVs make up large portions of vertebrate genomes and represent a record of past host–retrovirus interactions. We developed pan-vertebrate ERV analyses to provide an overview of host–retrovirus interactions, generating insights into retroviral evolution, diversity, host-switching, and the factors influencing retroviral transmission. Astoundingly, we found over 36,000 ERV lineages across our sample of vertebrate diversity. The results provide knowledge about host–retrovirus coevolution, suggesting an unprecedented ability of retroviruses to switch between distantly related vertebrates and implying existence of additional, yet unidentified retroviruses.

Significance
. Retroviruses, such as HIV, are important pathogens of vertebrates, including humans. They are capable of crossing species barriers to infect new hosts, but knowledge about the evolutionary history of retroviruses is limited. However, genomic traces of past retrovirus activities known as “endogenous retroviruses” can be screened from sequenced genomes and analyzed to improve understanding of retrovirus evolution. Here we use a unique approach to address the evolution of one group of retroviruses in a screen of 60 diverse vertebrate host genomes. We find evidence of rampant host-switching across mammalian orders by members of this group throughout their evolutionary history. We also find evidence that the spread of infective retroviruses from this group may be facilitated by rats.

ABSTRACTBackgroundTransposable elements (TEs) are found in nearly all eukaryotic genomes and are implicated in a range of evolutionary processes. Despite considerable research attention on TEs, their annotation and characterisation remain challenging, particularly for non-specialists. Current methods of automated TE annotation are subject to several issues that can reduce their overall quality: (i) fragmented and overlapping TE annotations may lead to erroneous estimates of TE count and coverage; (ii) repeat models may represent small proportions of their total length, where 5’ and 3’ regions are poorly captured; (iii) resultant libraries may contain redundancy, with the same TE family represented more than once. Existing pipelines can also be challenging to install, run, and extract data from. To address these issues, we present Earl Grey: a fully automated transposable element annotation pipeline designed for the user-friendly curation and annotation of TEs in eukaryotic genome assemblies.ResultsUsing a simulated genome, three model genome assemblies, and three non-model genome assemblies, Earl Grey outperforms current widely used TE annotation methodologies in ameliorating the issues mentioned above by producing longer TE consensus sequences in non-redundant TE libraries, which are then used to produce less fragmented TE annotations without the presence of overlaps. Earl Grey scores highly in benchmarking for TE annotation (MCC: 0.99) and classification (97% correctly classified) in comparison to existing software.ConclusionsEarl Grey provides a comprehensive and fully automated TE annotation toolkit that provides researchers with paper-ready summary figures and outputs in standard formats compatible with other bioinformatics tools. Earl Grey has a modular format, with great scope for the inclusion of additional modules focussed on further quality control aspects and tailored analyses in future releases.

Abstract
. Background
Transposable elements (TEs) are found in nearly all eukaryotic genomes and are implicated in a range of evolutionary processes. Despite considerable research attention on TEs, their annotation and characterisation remain challenging, particularly for non-specialists. Current methods of automated TE annotation are subject to several issues that can reduce their overall quality: (i) fragmented and overlapping TE annotations may lead to erroneous estimates of TE count and coverage; (ii) repeat models may represent small proportions of their total length, where 5’ and 3’ regions are poorly captured; (iii) resultant libraries may contain redundancy, with the same TE family represented more than once. Existing pipelines can also be challenging to install, run, and extract data from. To address these issues, we present Earl Grey: a fully automated transposable element annotation pipeline designed for the user-friendly curation and annotation of TEs in eukaryotic genome assemblies.
Results
Using a simulated genome, three model genome assemblies, and three non-model genome assemblies, Earl Grey outperforms current widely used TE annotation methodologies in ameliorating the issues mentioned above by producing longer TE consensus sequences in non-redundant TE libraries, which are then used to produce less fragmented TE annotations without the presence of overlaps. Earl Grey scores highly in benchmarking for TE annotation (MCC: 0.99) and classification (97% correctly classified) in comparison to existing software.
Conclusions
Earl Grey provides a comprehensive and fully automated TE annotation toolkit that provides researchers with paper-ready summary figures and outputs in standard formats compatible with other bioinformatics tools. Earl Grey has a modular format, with great scope for the inclusion of additional modules focussed on further quality control aspects and tailored analyses in future releases.

AbstractBackgroundLepidoptera (butterflies and moths) are an important model system in ecology and evolution. A high-quality chromosomal genome assembly is available for the monarch butterfly (Danaus plexippus), but it lacks an in-depth transposable element (TE) annotation, presenting an opportunity to explore monarch TE dynamics and the impact of TEs on shaping the monarch genome.ResultsWe find 6.21% of the monarch genome is comprised of TEs, a reduction of 6.85% compared to the original TE annotation performed on the draft genome assembly. Monarch TE content is low compared to two closely related species with available genomes, Danaus chrysippus (33.97% TE) and Danaus melanippus (11.87% TE). The biggest TE contributions to genome size in the monarch are LINEs and Penelope-like elements, and three newly identified families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE), collectively contribute 34.92% of total TE content. We find evidence of recent TE activity, with two novel Tc1 families rapidly expanding over recent timescales (tc1-1_dPle, tc1- 2_dPle). LINE fragments show signatures of genomic deletions indicating a high rate of TE turnover. We investigate associations between TEs and wing colouration and immune genes and identify a three-fold increase in TE content around immune genes compared to other host genes.ConclusionsWe provide a detailed TE annotation and analysis for the monarch genome, revealing a considerably smaller TE contribution to genome content compared to two closely related Danaus species with available genome assemblies. We identify highly successful novel DNA TE families rapidly expanding over recent timescales, and ongoing signatures of both TE expansion and removal highlight the dynamic nature of repeat content in the monarch genome. Our findings also suggest that insect immune genes are promising candidates for future interrogation of TE-mediated host adaptation.

Background
Lepidoptera (butterflies and moths) are an important model system in ecology and evolution. A high-quality chromosomal genome assembly is available for the monarch butterfly (Danaus plexippus), but it lacks an in-depth transposable element (TE) annotation, presenting an opportunity to explore monarch TE dynamics and the impact of TEs on shaping the monarch genome.

Results
We find 6.21% of the monarch genome is comprised of TEs, a reduction of 6.85% compared to the original TE annotation performed on the draft genome assembly. Monarch TE content is low compared to two closely related species with available genomes, Danaus chrysippus (33.97% TE) and Danaus melanippus (11.87% TE). The biggest TE contributions to genome size in the monarch are LINEs and Penelope-like elements, and three newly identified families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE), collectively contribute 34.92% of total TE content. We find evidence of recent TE activity, with two novel Tc1 families rapidly expanding over recent timescales (tc1-1_dPle, tc1- 2_dPle). LINE fragments show signatures of genomic deletions indicating a high rate of TE turnover. We investigate associations between TEs and wing colouration and immune genes and identify a three-fold increase in TE content around immune genes compared to other host genes.

Conclusions
We provide a detailed TE annotation and analysis for the monarch genome, revealing a considerably smaller TE contribution to genome content compared to two closely related Danaus species with available genome assemblies. We identify highly successful novel DNA TE families rapidly expanding over recent timescales, and ongoing signatures of both TE expansion and removal highlight the dynamic nature of repeat content in the monarch genome. Our findings also suggest that insect immune genes are promising candidates for future interrogation of TE-mediated host adaptation.

AbstractThe Myriapoda including millipedes and centipedes is of major importance in terrestrial ecology and nutrient recycling. Here, we sequenced and assembled two chromosomal-scale genomes of millipedes Helicorthomorpha holstii (182 Mb, N50 18.11 Mb mainly on 8 pseudomolecules) and Trigoniulus corallinus (449 Mb, N50 26.78 Mb mainly on 15 pseudomolecules). Unique defense systems, genomic features, and patterns of gene regulation in millipedes, not observed in other arthropods, are revealed. Millipedes possesses a unique ozadene defensive gland unlike the venomous forcipules in centipedes. Sets of genes associated with anti-microbial activity are identified with proteomics, suggesting that the ozadene gland is not primarily an antipredator adaptation (at least in T. corallinus). Macro-synteny analyses revealed highly conserved genomic blocks between centipede and the two millipedes. Tight Hox and the first loose ecdysozoan ParaHox homeobox clusters are identified, and a myriapod-specific genomic rearrangement including Hox3 is also observed. The Argonaute proteins for loading small RNAs are duplicated in both millipedes, but unlike insects, an argonaute duplicate has become a pseudogene. Evidence of post-transcriptional modification in small RNAs, including species-specific microRNA arm switching that provide differential gene regulation is also obtained. Millipede genomes reveal a series of unique genomic adaptations and microRNA regulation mechanisms have occurred in this major lineage of arthropod diversity. Collectively, the two millipede genomes shed new light on this fascinating but poorly understood branch of life, with a highly unusual body plan and novel adaptations to their environment.

Retroviruses are a virus family of considerable medical and veterinary importance. Additionally, it is now clear that endogenous retroviruses (ERVs) comprise significant portions of vertebrate genomes. Until recently, very little was known about the deep evolutionary origins of retroviruses. However, advances in genomics and informatics have opened the way for great strides in understanding. Recent research employing a wide variety of bioinformatic approaches has demonstrated that retroviruses evolved during the early Palaeozoic Era, between 460-550 million years ago, providing the oldest inferred date estimate for any virus group. This finding presents an important framework to investigate the evolutionary transitions that led to the emergence of the retroviruses, offering potential insights into the infectious origins of a major group of pathogenic viruses.

Chromosome-scale genome assemblies based on ultralong-read sequencing technologies are able to illuminate previously intractable aspects of genome biology such as fine-scale centromere structure and large-scale variation in genome features such as heterochromatin, GC content, recombination rate, and gene content. We present here a new chromosome-scale genome of the Mongolian gerbil (Meriones unguiculatus), which includes the complete sequence of all centromeres. Gerbils are thus the one of the first vertebrates to have their centromeres completely sequenced. Gerbil centromeres are composed of four different repeats of length 6, 37, 127, or 1,747�
bp, which occur in simple alternating arrays and span 1-6�
Mb. Gerbil genomes have both an extensive set of GC-rich genes and chromosomes strikingly enriched for constitutive heterochromatin. We sought to determine if there was a link between these two phenomena and found that the two heterochromatic chromosomes of the Mongolian gerbil have distinct underpinnings: Chromosome 5 has a large block of intraarm heterochromatin as the result of a massive expansion of centromeric repeats, while chromosome 13 is comprised of extremely large (>150�
kb) repeated sequences. In addition to characterizing centromeres, our results demonstrate the importance of including karyotypic features such as chromosome number and the locations of centromeres in the interpretation of genome sequence data and highlight novel patterns involved in the evolution of chromosomes.

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.

BACKGROUND: Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored. RESULTS: Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies. CONCLUSIONS: Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.

Transposable elements (TEs) are mobile genetic sequences present within host genomes. TEs can contribute to the evolution of host traits, since transposition is mutagenic and TEs often contain host regulatory and protein coding sequences. We review cases where TEs influence animal colouration, reporting major patterns and outstanding questions. TE-induced colouration phenotypes typically arise via introduction of novel regulatory sequences and splice sites, affecting pigment cell development or pigment synthesis. We discuss if particular TE types may be more frequently involved in the evolution of colour variation in animals, given that examples involving long terminal repeat (LTR) elements appear to dominate. Currently, examples of TE-induced colouration phenotypes in animals mainly concern model and domesticated insect and mammal species. However, several influential recent examples, coupled with increases in genome sequencing, suggest cases reported from wild species will increase considerably.

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

Accumulation of plastics in the marine environment has widespread detrimental consequences for ecosystems and wildlife. Marine plastics are rapidly colonised by a wide diversity of bacteria, including human pathogens, posing potential risks to health. Here, we investigate the effect of polymer type, residence time and estuarine location on bacterial colonisation of common household plastics, including pathogenic bacteria. We submerged five main household plastic types: low-density PE (LDPE), high-density PE (HDPE), polypropylene (PP), polyvinyl chloride (PVC) and polyethylene terephthalate (PET) at an estuarine site in Cornwall (U.K.) and tracked bacterial colonisation dynamics. Using both culture-dependent and culture-independent approaches, we found that bacteria rapidly colonised plastics irrespective of polymer type, reaching culturable densities of up to 1000 cells cm3 after 7 weeks. Community composition of the biofilms changed over time, but not among polymer types. The presence of pathogenic bacteria, quantified using the insect model Galleria mellonella, increased dramatically over a five-week period, with Galleria mortality increasing from 4% in week one to 65% in week five. No consistent differences in virulence were observed between polymer types. Pathogens isolated from plastic biofilms using Galleria enrichment included Serratia and Enterococcus species and they harboured a wide range of antimicrobial resistance genes. Our findings show that plastics in coastal waters are rapidly colonised by a wide diversity of bacteria independent of polymer type. Further, our results show that marine plastic biofilms become increasingly associated with virulent bacteria over time.

The moth Heortia vitessoides Moore (Lepidoptera: Crambidae) is a major pest of ecologically, commercially and culturally important agarwood-producing trees in the genus Aquilaria. In particular, H. vitessoides is one of the most destructive defoliating pests of the incense tree Aquilaria sinesis, which produces a valuable fragrant wood used as incense and in traditional Chinese medicine [33]. Nevertheless, a genomic resource for H. vitessoides is lacking. Here, we present a chromosomal-level assembly for H. vitessoides, consisting of a 517 megabase (Mb) genome assembly with high physical contiguity (scaffold N50 of 18.2 Mb) and high completeness (97.9% complete BUSCO score). To aid gene annotation, 8 messenger RNA transcriptomes from different developmental stages were generated, and a total of 16,421 gene models were predicted. Expansion of gene families involved in xenobiotic metabolism and development were detected, including duplications of cytosolic sulfotransferase (SULT) genes shared among lepidopterans. In addition, small RNA sequencing of 5 developmental stages of H. vitessoides facilitated the identification of 85 lepidopteran conserved microRNAs, 94 lineage-specific microRNAs, as well as several microRNA clusters. A large proportion of the H. vitessoides genome consists of repeats, with a 29.12% total genomic contribution from transposable elements, of which long interspersed nuclear elements (LINEs) are the dominant component (17.41%). A sharp decrease in the genome-wide percentage of LINEs with lower levels of genetic distance to family consensus sequences suggests that LINE activity has peaked in H. vitessoides. In contrast, opposing patterns suggest a substantial recent increase in DNA and LTR element activity. Together with annotations of essential sesquiterpenoid hormonal pathways, neuropeptides, microRNAs and transposable elements, the high-quality genomic and transcriptomic resources we provide for the economically important moth H. vitessoides provide a platform for the development of genomic approaches to pest management, and contribute to addressing fundamental research questions in Lepidoptera.

Infection from pathogens is consequential in host survival and fitness. Genetic differences in immune components can provide benefits in defending a host from pathogens. The Toll-like receptors (TLRs) are a family on innate immune receptors which detect pathogen components and have been found to show signs of positive selection across taxa. Genetic variation within TLRs has also been shown to relate to susceptibility to infection. In the first chapter I examine signatures of positive selection across the order Carnivora. The chapter aims to look for trends within and between TLRs using tests for selection, location of positively selected codons (PSCs), and TLR cellular location. I also investigate strength of selection over each species’ evolutionary history to find relationships with ecological and life history characteristics. The majority of TLRs show robust evidence of positive selection, with PSCs located within areas of biological significance for pathogen interaction. Stronger signs of positive selection were found in extra-cellular receptors than intra-cellular. Significant relationships between strength of selection and ecological characteristics and life history traits were found in multiples TLRs. Most notably, these were with the reproductive traits of gestation length, interbirth interval, and litter size, which potentially illustrate theories of life history trade off in investment between reproduction and immune function. The second chapter of this thesis focuses on genetic and behavioural defences which are crucial in host-pathogen relationships but not often combined into one study. This study looks at population-level TLR polymorphism in the bank vole (Myodes glareolus) and boldness and exploration personality scores. The haplotypes of the voles formed three clusters (A, B and C) but no significant relationships were found between the clusters and personality. However, there was a trend between B cluster haplotypes and increased explorations scores. These prospective relationships have implications in the maintenance of personality and genetic variation within populations and their roles in the spread of disease. The study highlights the importance of considering multiple aspects of host immune defence and their interactions rather than individualising the components. Overall, interaction with pathogens plays a large role in patterns in selection seen across taxa and within species. This thesis shows how the inclusion and consideration of ecological and behavioural factors into immunogenetics can provide a more representative idea of the complex interactions that are at work in host immune defences.

BACKGROUND: Schistosomiasis, or bilharzia, is a parasitic disease caused by trematode flatworms of the genus Schistosoma. Infection by Schistosoma mansoni in humans results when cercariae emerge into water from freshwater snails in the genus Biomphalaria and seek out and penetrate human skin. The snail Biomphalaria straminea is native to South America and is now also present in Central America and China, and represents a potential vector host for spreading schistosomiasis. To date, genomic information for the genus is restricted to the neotropical species Biomphalaria glabrata. This limits understanding of the biology and management of other schistosomiasis vectors, such as B. straminea. FINDINGS: Using a combination of Illumina short-read, 10X Genomics linked-read, and Hi-C sequencing data, our 1.005 Gb B. straminea genome assembly is of high contiguity, with a scaffold N50 of 25.3 Mb. Transcriptomes from adults were also obtained. Developmental homeobox genes, hormonal genes, and stress-response genes were identified, and repeat content was annotated (40.68% of genomic content). Comparisons with other mollusc genomes (including Gastropoda, Bivalvia, and Cephalopoda) revealed syntenic conservation, patterns of homeobox gene linkage indicative of evolutionary changes to gene clusters, expansion of heat shock protein genes, and the presence of sesquiterpenoid and cholesterol metabolic pathway genes in Gastropoda. In addition, hormone treatment together with RT-qPCR assay reveal a sesquiterpenoid hormone responsive system in B. straminea, illustrating that this renowned insect hormonal system is also present in the lophotrochozoan lineage. CONCLUSION: This study provides the first genome assembly for the snail B. straminea and offers an unprecedented opportunity to address a variety of phenomena related to snail vectors of schistosomiasis, as well as evolutionary and genomics questions related to molluscs more widely.

Insects are capable of extraordinary feats of long-distance movement that have profound impacts on the function of terrestrial ecosystems. The ability to undertake these movements arose multiple times through the evolution of a suite of traits that make up the migratory syndrome, however the underlying genetic pathways involved remain poorly understood. Migratory hoverflies (Diptera: Syrphidae) are an emerging model group for studies of migration. They undertake seasonal movements in huge numbers across large parts of the globe and are important pollinators, biological control agents and decomposers. Here, we assembled a high-quality draft genome of the marmalade hoverfly (Episyrphus balteatus). We leveraged this genomic resource to undertake a genome-wide transcriptomic comparison of actively migrating Episyrphus, captured from a high mountain pass as they flew south to overwinter, with the transcriptomes of summer forms which were non-migratory. We identified 1543 genes with very strong evidence for differential expression. Interrogation of this gene set reveals a remarkable range of roles in metabolism, muscle structure and function, hormonal regulation, immunity, stress resistance, flight and feeding behaviour, longevity, reproductive diapause and sensory perception. These features of the migrant phenotype have arisen by the integration and modification of pathways such as insulin signalling for diapause and longevity, JAK/SAT for immunity, and those leading to octopamine production and fuelling to boost flight capabilities. Our results provide a powerful genomic resource for future research, and paint a comprehensive picture of global expression changes in an actively migrating insect, identifying key genomic components involved in this important life-history strategy.

Transposable elements (TEs) are DNA sequences with the capability to move within a genome. Whilst often detrimental, there are also multiple examples demonstrating the potential of TEs for the evolution of their hosts, including in the evolution of insecticide resistance, mammalian pregnancy, and adaptive immunity. Despite these examples, open questions remain regarding the extent to which TEs provide a general mechanism for host evolution, and the importance of TEs as a source of genetic variation compared to other contributions.
To address broadscale questions on TE-host dynamics, comparative genomic studies can be used to assess the contributions of TEs to host genomes, whilst also uncovering the interactions which can often lead to large differences in TE abundance and diversity among species, even within a single genus. Such studies are timely due to the increasing availability of high-quality genome assemblies for non-model organisms. However, processing large numbers of genome assemblies to understand TE-host dynamics requires effective automated TE annotation methodologies, as manual annotation is unviable for studies considering more than a handful of genome assemblies. To address this, in Chapter Two, a fully-automated TE curation and annotation pipeline named ‘Earl Grey’ is developed, which aims to address some of the core issues regarding TE annotation, namely poorly-defined TE boundaries and redundant consensus sequences in TE consensus libraries, and fragmented and overlapping TE annotations. Earl Grey aims to outperform other widely-used TE annotation tools and produces outputs in standard formats for compatibility with downstream analyses, along with summary figures. Earl Grey is capable of analysing large numbers of genomes without any intervention, and is an effective tool for large-scale comparative analyses incorporating large numbers of genome assemblies. The aim is for Earl Grey to continue to incorporate new modules using feedback from users and advances in TE annotation methodologies, with the aim of becoming a community-led TE annotation and curation tool.
in Chapter Three, a recent high-quality chromosomal assembly of the monarch butterfly (Danaus plexippus) is used for an in-depth exploration into the impacts of TEs in shaping the host genome by examining TE expansions, host interactions, and removal, whilst the availability of two other Danaus genomes provide a comparative context. The TE content of the monarch was found to be much lower than the content annotated in the original draft genome assembly (This study: 6.21%, Draft Genome: 13.06%), and also to be much lower than in other Danaus genomes (D. melanippus: 11.87%, D. chrysippus: 33.97%). The reduction in annotated TE content compared to the original draft genome is attributed to an improved understanding of TE structure resulting in the exclusion of erroneous annotations previously annotated as TE, along with the conservative annotation approach taken in this study. The differences among Danaus species were not due to variation in DNA deletion rates, but are hypothesised to be due to expansions in lineage-specific TEs. Three newly-identified TE families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE) contribute over one third of the total TE content in the monarch. Historical bursts of LINE activity are evident in the monarch genome. However, just two novel Tc1 families (tc1-1_dPle and tc1-2_dPle) have rapidly expanded over the last ~500,000 years. TE content was found to be unevenly distributed between different genomic compartments, with a strong negative correlation between gene density and TE density. Six gene hotspots containing putatively important host functions were significantly depleted of TEs, potentially reflecting the deleterious consequences of TE insertions and the selection against TEs detrimental to host fitness. There is evidence of LINE and Penelope-like element removal via both dissociation of transcription machinery and genomic deletions, and the presence of swathes of small fragments suggests rapid turnover of Non-LTR TEs. This contrasts with patterns observed in mammals, where lower rates of TE turnover result in the accumulation of more ancient TEs. Together, the findings presented in Chapter Three demonstrate the ongoing dynamic nature of the interactions between TEs and their host genomes, with ongoing TE activity having the potential to considerably alter host genomes over evolutionary timescales, which can drive significant differences even among very closely related species.
Aphids are destructive crop pests, and several species have evolved multiple resistance to insecticides. Meanwhile, TEs have been implicated in the evolution of insecticide resistance. Consequently, in Chapter Four, 21 available aphid genomes are analysed to explore the extent to which TEs act as a general source of genomic novelty for host evolution. TEs were found to be significantly enriched at xenobiotic resistance loci (XRL) when compared to other genes, and showed enrichment levels similar to housekeeping genes. However, unlike at housekeeping genes, the maintenance of TEs around XRL is unlikely to arise through constitutive expression and associated open chromatin, and is more likely to arise via selection for insertions in these regions, as XRL are not expressed in germline cells. Further, TEs are also enriched around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides in three aphids of agricultural importance (A. pisum, M. cerasi, and M. persicae). Together, these findings suggest that TEs are selected around XRL in aphids for beneficial purposes.
in Chapter Five, 88 high quality genome assemblies are analysed to uncover the physiological and ecological determinants leading to variation in the TE content across butterflies. TE content, as a proportion of total genome size, is highly variable across butterflies, ranging between 6.21% in Danaus plexippus to 67.55% in Satyrium esculi. A strong phylogenetic signal was found in both TE abundance and diversity, indicating that TE abundance and diversity in extant butterflies are good indicators of the evolutionary past. Three life history traits were found to strongly correlate with TE abundance (forewing size, voltinism, and species distribution), whilst two were found to strongly correlate with TE diversity (voltinism and species distribution). All strongly correlated life history traits impacted TE abundance and diversity in the negative direction. For voltinism, this confirmed the hypothesis, where multivoltinism is predicted to result in more efficient purging of TE insertions from host genomes. However, this directionality contrasts hypotheses for forewing size and species distribution, where larger butterflies and those with larger distributions were expected to have higher TE abundances and diversities through more ecological opportunities raising the potential for invasion of novel TEs into host genomes through processes such as horizontal transfer. Overall, these findings highlight the significant impacts that ecological and physiological characteristics of species can have on TE abundance and diversity.
This thesis aims to provide methodologies for the automated annotation and curation of TEs, and apply these to further our understanding of TE-host dynamics. The development of an improved automated TE annotation tool highlights the continued need for enhanced methodologies to advance the field of TE biology. Community-led efforts have the potential to provide a significant benefit on this front through the combination of expertise to produce a consistent TE annotation tool of benefit both within the field of TE biology and more widely. Limited understanding of the processes leading to large differences in TE content among closely related species, maintenance of TEs around genes associated with processes under strong selection, and the ecological drivers influencing TE diversity and abundance, highlight the need for further large-scale comparative studies.

Animals display a fascinating diversity of body plans. Correspondingly, genomic analyses have revealed dynamic evolution of gene gains and losses among animal lineages. Here we sequence six new myriapod genomes (three millipedes, three centipedes) at key phylogenetic positions within this major but understudied arthropod lineage. We combine these with existing genomic resources to conduct a comparative analysis across all available myriapod genomes. We find that millipedes generally have considerably smaller genomes than centipedes, with the repeatome being a major contributor to genome size, driven by independent large gains of transposons in three centipede species. In contrast to millipedes, centipedes gained a large number of gene families after the subphyla diverged, with gains contributing to sensory and locomotory adaptations that facilitated their ecological shift to predation. We identify distinct horizontal gene transfer (HGT) events from bacteria to millipedes and centipedes, with no identifiable HGTs shared among all myriapods. Loss of juvenile hormone O-methyltransferase, a key enzyme in catalysing sesquiterpenoid hormone production in arthropods, was also revealed in all millipede lineages. Our findings suggest that the rapid evolution of distinct genomic pathways in centipede and millipede lineages following their divergence from the myriapod ancestor, was shaped by differing ecological pressures.

We present a genome assembly from an individual male Pyrgus malvae (the grizzled skipper; Arthropoda; Insecta; Lepidoptera; Hesperiidae). The genome sequence is 725 megabases in span. The majority (99.97%) of the assembly is scaffolded into 31 chromosomal pseudomolecules, with the Z sex chromosome assembled.

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 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.

Transposable elements are known by many names, including 'transposons', 'interspersed repeats', 'selfish genetic elements', 'jumping genes', and 'parasitic DNA', but here we will refer to them simply as transposable elements. Many biologists will have heard of transposable elements and their ability to transpose (change position) within the genome. But fewer may be aware of their varied influences on host biology, including contributions to the evolution of diverse host traits such as internal gestation, memory, colouration, and adaptive immunity. Transposable elements are a near ubiquitous feature of eukaryotic genomes, and they often comprise a substantial proportion of total genomic content. Consequently, transposable element genes are considered among the most abundant coding sequences in nature. Recent advances in genome sequencing have ushered in a golden age for transposable-element research, providing opportunities to greatly improve our understanding of the effects of transposable elements on host evolution and disease. However, our ability to detect and analyse transposable elements still faces significant challenges, impairing efforts to decipher their evolution, characterise their diversity, and elucidate their myriad host influences. Below, we summarise key aspects of transposable element biology in eukaryotes and discuss major outstanding research questions.

The aphid Myzus persicae is a destructive agricultural pest that displays an exceptional ability to develop resistance to both natural and synthetic insecticides. To investigate the evolution of resistance in this species we generated a chromosome-scale genome assembly and living panel of >110 fully sequenced globally sampled clonal lines. Our analyses reveal a remarkable diversity of resistance mutations segregating in global populations of M. persicae. We show that the emergence and spread of these mechanisms is influenced by host-plant associations, uncovering the widespread co-option of a host-plant adaptation that also offers resistance against synthetic insecticides. We identify both the repeated evolution of independent resistance mutations at the same locus, and multiple instances of the evolution of novel resistance mechanisms against key insecticides. Our findings provide fundamental insights into the genomic responses of global insect populations to strong selective forces, and hold practical relevance for the control of pests and parasites.

Whole genome duplication (WGD) has occurred in relatively few sexually reproducing invertebrates. Consequently, the WGD that occurred in the common ancestor of horseshoe crabs ~135 million years ago provides a rare opportunity to decipher the evolutionary consequences of a duplicated invertebrate genome. Here, we present a high-quality genome assembly for the mangrove horseshoe crab Carcinoscorpius rotundicauda (1.7 Gb, N50 = 90.2 Mb, with 89.8% sequences anchored to 16 pseudomolecules, 2n = 32), and a resequenced genome of the tri-spine horseshoe crab Tachypleus tridentatus (1.7 Gb, N50 = 109.7 Mb). Analyses of gene families, microRNAs, and synteny show that horseshoe crabs have undergone three rounds (3R) of WGD. Comparison of C. rotundicauda and T. tridentatus genomes from populations from several geographic locations further elucidates the diverse fates of both coding and noncoding genes. Together, the present study represents a cornerstone for improving our understanding of invertebrate WGD events on the evolutionary fates of genes and microRNAs, at both the individual and population level. We also provide improved genomic resources for horseshoe crabs, of applied value for breeding programs and conservation of this fascinating and unusual invertebrate lineage.

The green peach aphid, Myzus persicae, is a globally distributed highly damaging crop pest. This species has demonstrated an exceptional ability to evolve resistance to both synthetic insecticides used for control, and natural insecticides produced by certain plants as a chemical defense against insect attack. Here we review work characterizing the evolution of resistance in M. persicae to the natural insecticide nicotine and the structurally related class of synthetic neonicotinoid insecticides. We outline how research on this topic has provided insights into long-standing questions of both evolutionary and applied importance. These include questions pertaining to the origins of novel traits, the number and nature of mutational events or 'adaptive steps' underlying the evolution of new phenotypes, and whether host plant adaptations can be co-opted to confer resistance to synthetic insecticides. Finally, research on the molecular mechanisms underlying insecticide resistance in M. persicae has generated several outstanding questions on the genetic architecture of resistance to both natural and synthetic xenobiotics, and we conclude by identifying key knowledge gaps for future research. © 2021 the Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

A fundamental question in evolutionary biology is how microevolutionary processes translate into species diversification. Cophylogeny provides an appropriate framework to address this for symbiotic associations, but historically has been primarily limited to unveiling patterns. We argue that it is essential to integrate advances from ecology and evolutionary biology into cophylogeny, to gain greater mechanistic insights and transform cophylogeny into a platform to advance understanding of interspecific interactions and diversification more widely. We discuss key directions, such as incorporating trait reconstruction and considering multiple scales of network organization, and highlight recent developments for implementation. A new quantitative framework is proposed to allow integration of relevant information, such as quantitative traits and assessment of the contribution of individual mechanisms to cophylogenetic patterns.

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.

We present a genome assembly from an individual female Melitaea athalia (also known as Mellicta athalia; the heath fritillary; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 610 megabases in span. In total, 99.98% of the assembly is scaffolded into 32 chromosomal pseudomolecules, with the W and Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 12,824 protein coding genes.

We present a genome assembly from an individual male Celastrina argiolus) (the holly blue; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequence is 499 megabases in span. The majority (99.99%) of the assembly is scaffolded into 26 chromosomal pseudomolecules, with the Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 12,199 protein coding genes.

We present genome assemblies from a male and female Pieris napi (the green-veined white; Arthropoda; Insecta; Lepidoptera; Pieridae). The genome sequences of the male and female are 320 and 319 megabases in span, respectively. The majority of the assembly (99.79% of the male assembly, 99.88% of the female) is scaffolded into 24 autosomal pseudomolecules, with the Z sex chromosome assembled for the male and Z and W chromosomes assembled for the female. Gene annotation of the male assembly on Ensembl has identified 13,221 protein coding genes.

The evolution of resistance to insecticides threatens the sustainable control of many of the world's most damaging insect crop pests and disease vectors. To effectively combat resistance, it is important to understand its underlying genetic architecture, including the type and number of genetic variants affecting resistance and their interactions with each other and the environment. While significant progress has been made in characterizing the individual genes or mutations leading to resistance, our understanding of how genetic variants interact to influence its phenotypic expression remains poor. Here, we uncover a mechanism of insecticide resistance resulting from transposon-mediated insertional mutagenesis of a genetically dominant but insecticide-susceptible allele that enables the adaptive potential of a previously unavailable recessive resistance allele to be unlocked. Specifically, we identify clones of the aphid pest Myzus persicae that carry a resistant allele of the essential voltage-gated sodium channel (VGSC) gene with the recessive M918T and L1014F resistance mutations, in combination with an allele lacking these mutations but carrying a Mutator-like element transposon insertion that disrupts the coding sequence of the VGSC. This results in the down-regulation of the dominant susceptible allele and monoallelic expression of the recessive resistant allele, rendering the clones resistant to the insecticide bifenthrin. These findings are a powerful example of how transposable elements can provide a source of evolutionary potential that can be revealed by environmental and genetic perturbation, with applied implications for the control of highly damaging insect pests.

AbstractTrees in the genus Aquilaria (Thymelaeaceae) are known as lign aloes, and are native to the forests of southeast Asia. Lign aloes produce agarwood as an antimicrobial defence. Agarwood has a long history of cultural and medicinal use, and is of considerable commercial value. However, due to habitat destruction and over collection, lign aloes are threatened in the wild. We present a chromosomal‐level assembly for Aquilaria sinensis, a lign aloe endemic to China known as the incense tree, based on Illumina short‐read, 10X Genomics linked‐read, and Hi‐C sequencing data. Our 783.8 Mbp A. sinensis genome assembly is of high physical contiguity, with a scaffold N50 of 87.6 Mbp, and high completeness, with a 95.8% BUSCO score for eudicotyledon genes. We include 17 transcriptomes from various plant tissues, providing a total of 35,965 gene models. We reveal the first complete set of genes involved in sesquiterpenoid production, plant defence, and agarwood production for the genus Aquilaria, including genes involved in the biosynthesis of sesquiterpenoids via the mevalonic acid (MVA), 1‐deoxy‐D‐xylulose‐5‐phosphate (DXP), and methylerythritol phosphate (MEP) pathways. We perform a detailed repeat content analysis, revealing that transposable elements account for ~61% of the genome, with major contributions from gypsy‐like and copia‐like LTR retroelements. We also provide a comparative analysis of repeat content across sequenced species in the order Malvales. Our study reveals the first chromosomal‐level genome assembly for a tree in the genus Aquilaria and provides an unprecedented opportunity to address a variety of applied, genomic and evolutionary questions in the Thymelaeaceae more widely.

Background: Teleost fish play important roles in aquatic ecosystems and aquaculture. Threadfins (Perciformes: Polynemidae) show a range of interesting biology, and are of considerable importance for both wild fisheries and aquaculture. Additionally, the four-finger threadfin Eleutheronema tetradactylum is of conservation relevance since its populations are considered to be in rapid decline and it is classified as endangered. However, no genomic resources are currently available for the threadfin family Polynemidae. Results: We sequenced and assembled the first threadfin fish genome, the four-finger threadfin E. tetradactylum. We provide a genome assembly for E. tetradactylum with high contiguity (scaffold N50 = 56.3 kb) and high BUSCO completeness at 96.5%. The assembled genome size of E. tetradactylum is just 610.5 Mb, making it the second smallest perciform genome assembled to date. Just 9.07–10.91% of the genome sequence was found to consist of repetitive elements (standard RepeatMasker analysis vs custom analysis), making this the lowest repeat content identified to date for any perciform fish. A total of 37,683 protein-coding genes were annotated, and we include analyses of developmental transcription factors, including the Hox, ParaHox, and Sox families. MicroRNA genes were also annotated and compared with other chordate lineages, elucidating the gains and losses of chordate microRNAs. Conclusions: the four-finger threadfin E. tetradactylum genome presented here represents the first available genome sequence for the ecologically, biologically, and commercially important clade of threadfin fish. Our findings provide a useful genomic resource for future research into the interesting biology and evolution of this valuable group of food fish.

The Myriapoda, composed of millipedes and centipedes, is a fascinating but poorly understood branch of life, including species with a highly unusual body plan and a range of unique adaptations to their environment. Here, we sequenced and assembled 2 chromosomal-level genomes of the millipedes Helicorthomorpha holstii (assembly size = 182 Mb; shortest scaffold/contig length needed to cover 50% of the genome [N50] = 18.11 Mb mainly on 8 pseudomolecules) and Trigoniulus corallinus (assembly size = 449 Mb, N50 = 26.78 Mb mainly on 17 pseudomolecules). Unique genomic features, patterns of gene regulation, and defence systems in millipedes, not observed in other arthropods, are revealed. Both repeat content and intron size are major contributors to the observed differences in millipede genome size. Tight Hox and the first loose ecdysozoan ParaHox homeobox clusters are identified, and a myriapod-specific genomic rearrangement including Hox3 is also observed. The Argonaute (AGO) proteins for loading small RNAs are duplicated in both millipedes, but unlike in insects, an AGO duplicate has become a pseudogene. Evidence of post-transcriptional modification in small RNAs-including species-specific microRNA arm switching-providing differential gene regulation is also obtained. Millipedes possesses a unique ozadene defensive gland unlike the venomous forcipules found in centipedes. We identify sets of genes associated with the ozadene that play roles in chemical defence as well as antimicrobial activity. Macro-synteny analyses revealed highly conserved genomic blocks between the 2 millipedes and deuterostomes. Collectively, our analyses of millipede genomes reveal that a series of unique adaptations have occurred in this major lineage of arthropod diversity. The 2 high-quality millipede genomes provided here shed new light on the conserved and lineage-specific features of millipedes and centipedes. These findings demonstrate the importance of the consideration of both centipede and millipede genomes-and in particular the reconstruction of the myriapod ancestral situation-for future research to improve understanding of arthropod evolution, and animal evolutionary genomics more widely.

BACKGROUND: Tc1/mariner transposons are widespread DNA transposable elements (TEs) that have made important contributions to the evolution of host genomic complexity in metazoans. However, the evolution and diversity of the Tc1/mariner superfamily remains poorly understood. Following recent developments in genome sequencing and the availability of a wealth of new genomes, Tc1/mariner TEs have been identified in many new taxa across the eukaryotic tree of life. To date, the majority of studies focussing on Tc1/mariner elements have considered only a single host lineage or just a small number of host lineages. Thus, much remains to be learnt about the evolution of Tc1/mariner TEs by performing analyses that consider elements that originate from across host diversity. RESULTS: We mined the non-redundant database of NCBI using BLASTp searches, with transposase sequences from a diverse set of reference Tc1/mariner elements as queries. A total of 5158 Tc1/mariner elements were retrieved and used to reconstruct evolutionary relationships within the superfamily. The resulting phylogeny is well resolved and includes several new groups of Tc1/mariner elements. In particular, we identify a new family of plant-genome restricted Tc1/mariner elements, which we call PlantMar. We also show that the pogo family is much larger and more diverse than previously appreciated, and we review evidence for a potential revision of its status to become a separate superfamily. CONCLUSIONS: Our study provides an overview of Tc1-mariner phylogeny and summarises the impressive diversity of Tc1-mariner TEs among sequenced eukaryotes. Tc1/mariner TEs are successful in a wide range of eukaryotes, especially unikonts (the taxonomic supergroup containing Amoebozoa, Opisthokonta, Breviatea, and Apusomonadida). In particular, ecdysozoa, and especially arthropods, emerge as important hosts for Tc1/mariner elements (except the PlantMar family). Meanwhile, the pogo family, which is by far the largest Tc1/mariner family, also includes many elements from fungal and chordate genomes. Moreover, there is evidence of the repeated exaptation of pogo elements in vertebrates, including humans, in addition to the well-known example of CENP-B. Collectively, our findings provide a considerable advancement in understanding of Tc1/mariner elements, and more generally they suggest that much work remains to improve understanding of the diversity and evolution of DNA TEs.

Background: Homeobox-containing genes encode crucial transcription factors involved in animal, plant and fungal development, and changes to homeobox genes have been linked to the evolution of novel body plans and morphologies. In animals, some homeobox genes are clustered together in the genome, either as remnants from ancestral genomic arrangements, or due to coordinated gene regulation. Consequently, analyses of homeobox gene organization across animal phylogeny provide important insights into the evolution of genome organization and developmental gene control, and their interaction. However, homeobox gene organization remains to be fully elucidated in several key animal ancestors, including those of molluscs, lophotrochozoans and bilaterians. Results: Here, we present a high-quality chromosome-level genome assembly of the Hong Kong oyster, Magallana hongkongensis (2n = 20), for which 93.2% of the genomic sequences are contained on 10 pseudomolecules (~ 758 Mb, scaffold N50 = 72.3 Mb). Our genome assembly was scaffolded using Hi-C reads, facilitating a larger scaffold size compared to the recently published M. hongkongensis genome of Peng et al. (Mol Ecol Resources, 2020), which was scaffolded using the Crassostrea gigas assembly. A total of 46,963 predicted gene models (45,308 protein coding genes) were incorporated in our genome, and genome completeness estimated by BUSCO was 94.6%. Homeobox gene linkages were analysed in detail relative to available data for other mollusc lineages. Conclusions: the analyses performed in this study and the accompanying genome sequence provide important genetic resources for this economically and culturally valuable oyster species, and offer a platform to improve understanding of animal biology and evolution more generally. Transposable element content is comparable to that found in other mollusc species, contrary to the conclusion of another recent analysis. Also, our chromosome-level assembly allows the inference of ancient gene linkages (synteny) for the homeobox-containing genes, even though a number of the homeobox gene clusters, like the Hox/ParaHox clusters, are undergoing dispersal in molluscs such as this oyster.

A complex series of mutational events protected a mutualistic symbiosis during the shift of an insect to a toxic host plant.

Dictyocaulus nematode worms live as parasites in the lower airways of ungulates and can cause significant disease in both wild and farmed hosts. This study represents the first population genetic analysis of large lungworms in wildlife. Specifically, we quantify genetic variation in Dictyocaulus lungworms from wild deer (red deer, fallow deer and roe deer) in Hungary, based on mitochondrial cytochrome c oxidase subunit 1 (cox1) sequence data, using population genetic and phylogenetic analyses. The studied Dictyocaulus taxa display considerable genetic diversity. At least one cryptic species and a new parasite-host relationship are revealed by our molecular study. Population genetic analyses for Dictyocaulus eckerti revealed high gene flow amongst weakly structured spatial populations that utilise the three host deer species considered here. Our results suggest that D. eckerti is a widespread generalist parasite in ungulates, with a diverse genetic backround and high evolutionary potential. In contrast, evidence of cryptic genetic structure at regional geographic scales was observed for Dictyocaulus capreolus, which infects just one host species, suggesting it is a specialist within the studied area. D. capreolus displayed lower genetic diversity overall, with only moderate gene flow compared to the closely related D. eckerti. We suggest that the differing vagility and dispersal behaviour of hosts are important contributing factors to the population structure of lungworms, and possibly other nematode parasites with single-host life cycles. Our findings are of relevance for the management of lungworms in deer farms and wild deer populations.

1. Insects with complete metamorphosis (holometaboly) are extremely successful, constituting over 60% of all described animal species. Complete metamorphosis confers significant advantages because it enables organisms to optimise life-history components through temporal partitioning, and thereby to exploit multiple ecological niches. Yet holometaboly can also impose costs, and several lineages have evolved life cycle modifications to avoid complete metamorphosis. 2. In this review, we discuss different strategies that have evolved that result in the loss of complete metamorphosis (type I and type II paedomorphosis). In addition, the ecological pressures and developmental modifications that facilitate this avoidance are considered, as well as the importance of life cycle complexity in life-history evolution. 3. Interestingly, only female holometabolous insects have entirely avoided complete metamorphosis, and it is always the ancestrally juvenile morphology that is retained. These findings point to a strong sex-biased trade-off between investment in reproduction and development. While the loss of complete metamorphosis in females has occurred independently on several occasions across holometabolous insects, only a small number of species possessing this ability have been described. 4. Thus, complete metamorphosis, which originated only once in insects, appears to have been almost fully retained. This indicates that significant modifications to the holometabolan metamorphic ground plan are highly constrained, and suggests that the transition to complete metamorphosis is evolutionarily irreversible.

. MicroRNAs (miRNAs) have recently risen to prominence as novel factors responsible for post-transcriptional regulation of gene expression. miRNA genes have been posited as highly conserved in the clades in which they exist. Consequently, miRNAs have been used as rare genome change characters to estimate phylogeny by tracking their gain and loss. However, their short length (21–23 bp) has limited their perceived utility in sequenced-based phylogenetic inference. Here, using reference taxa with established phylogenetic relationships, we demonstrate that miRNA sequences are of high utility in quantitative, rather than in qualitative, phylogenetic analysis. The clear orthology among miRNA genes from different species makes it straightforward to identify and align these sequences from even fragmentary datasets. We also identify significant sequence conservation in the regions directly flanking miRNA genes, and show that this too is of utility in phylogenetic analysis, as well as highlighting conserved regions that will be of interest to other fields. Employing miRNA sequences from 12 sequenced drosophilid genomes, together with a
. Tribolium castaneum
. outgroup, we demonstrate that this approach is robust using Bayesian and maximum-likelihood methods. The utility of these characters is further demonstrated in the rhabditid nematodes and primates. As next-generation sequencing makes it more cost-effective to sequence genomes and small RNA libraries, this methodology provides an alternative data source for phylogenetic analysis. The approach allows rapid resolution of relationships between both closely related and rapidly evolving species, and provides an additional tool for investigation of relationships within the tree of life.
.

Significance
. For millions of years retroviruses, such as HIV in humans, have attacked vertebrates. Occasionally retroviruses infiltrate germ cells, incorporate themselves into the host’s genome, and transmit vertically to the host’s offspring as endogenous retroviruses (ERVs). Consequently, ERVs make up large portions of vertebrate genomes and represent a record of past host–retrovirus interactions. We developed pan-vertebrate ERV analyses to provide an overview of host–retrovirus interactions, generating insights into retroviral evolution, diversity, host-switching, and the factors influencing retroviral transmission. Astoundingly, we found over 36,000 ERV lineages across our sample of vertebrate diversity. The results provide knowledge about host–retrovirus coevolution, suggesting an unprecedented ability of retroviruses to switch between distantly related vertebrates and implying existence of additional, yet unidentified retroviruses.

Significance
. Retroviruses, such as HIV, are important pathogens of vertebrates, including humans. They are capable of crossing species barriers to infect new hosts, but knowledge about the evolutionary history of retroviruses is limited. However, genomic traces of past retrovirus activities known as “endogenous retroviruses” can be screened from sequenced genomes and analyzed to improve understanding of retrovirus evolution. Here we use a unique approach to address the evolution of one group of retroviruses in a screen of 60 diverse vertebrate host genomes. We find evidence of rampant host-switching across mammalian orders by members of this group throughout their evolutionary history. We also find evidence that the spread of infective retroviruses from this group may be facilitated by rats.

AbstractInvestigating complex parasitic life cycles is important for understanding the major fitness components that drive the evolution of host–parasite systems. The rare condition of heterotrophic heteronomy, in which the sexes utilize disparate host taxa, is a poorly understood complex parasitic lifestyle. One of only two known examples occurs in the Myrmecolacidae, an unusual family of the parasitoid order Strepsiptera (Insecta), in which males parasitize ants while females parasitize grasshoppers, crickets, and praying mantids. Here, we reconstruct the evolutionary pattern and timescale of host‐use in a set of morphologically cryptic myrmecolacid taxa currently identified as Caenocholax fenyesi. We find that (i) C. fenyesi contains at least ten cryptic lineages consistent with separate species; (ii) Fossil evidence suggests a very low molecular clock rate and an ancient origin for cryptic lineages; (iii) Diversity among Caenocholax species is partitioned by geography and host association of the female; and (iv) Switches in host usage are uncoupled between the sexes, with changes in female host preference accompanying diversification. This study represents the first phylogeographical analysis of any strepsipteran, and the first molecular examination of host‐use for a heterotrophic heteronomous taxon. Our results have implications for the understanding of evolution, host usage and estimated species richness in parasitic taxa.

BACKGROUND: Biological invasions provide a window on the process of community assembly. In particular, tracking natural enemy recruitment to invading hosts can reveal the relative roles of co-evolution (including local adaptation) and ecological sorting. We use molecular data to examine colonisation of northern Europe by the parasitoid Megastigmus stigmatizans following invasions of its herbivorous oak gallwasp hosts from the Balkans. Local host adaptation predicts that invading gallwasp populations will have been tracked primarily by sympatric Balkan populations of M. stigmatizans (Host Pursuit Hypothesis). Alternatively, ecological sorting allows parasitoid recruitment from geographically distinct populations with no recent experience of the invading hosts (Host Shift Hypothesis). Finally, we test for long-term persistence of parasitoids introduced via human trade of their hosts' galls (Introduction Hypothesis). RESULTS: Polymorphism diagnostic of different southern refugial regions was present in both mitochondrial and nuclear microsatellite markers, allowing us to identify the origins of northern European invaded range M. stigmatizans populations. As with their hosts, some invaded range populations showed genetic variation diagnostic of Balkan sources, supporting the Host Pursuit Hypothesis. In contrast, other invading populations had an Iberian origin, unlike their hosts in northern Europe, supporting the Host Shift Hypothesis. Finally, both British and Italian M. stigmatizans populations show signatures compatible with the Introduction Hypothesis from eastern Mediterranean sources. CONCLUSIONS: These data reveal the continental scale of multi-trophic impacts of anthropogenic disturbance and highlight the fact that herbivores and their natural enemies may face very different constraints on range expansion. The ability of natural enemies to exploit ecologically-similar hosts with which they have had no historical association supports a major role for ecological sorting processes in the recent assembly of these communities. The multitude of origins of invading natural enemy populations in this study emphasises the diversity of mechanisms requiring consideration when predicting consequences of other biological invasions or biological control introductions.

Vitellogenins and other large lipid transfer proteins (LLTP) are well known to play significant roles in the development, metabolism and reproduction of animals. Comparative genomics and phylogenetic analyses of LLTPs using the most comprehensive dataset in metazoans to date are carried out. Our analyses demonstrate that LLTP genes arose significantly earlier, and are more widespread than previously proposed – being present in numerous additional bilaterian and non‐bilaterian lineages. A hypothesis is advanced that the most ancestral animal LLTP gene is Vtg, while loss of domains occurred at the bilaterians stem giving rise to apolipoprotein and microsomal triglyceride transfer proteins genes.

Little is known about the evolutionary history of most complex multi-trophic insect communities. Widespread species from different trophic levels might evolve in parallel, showing similar spatial patterns and either congruent temporal patterns (Contemporary Host-tracking) or later divergence in higher trophic levels (Delayed Host-tracking). Alternatively, host shifts by natural enemies among communities centred on different host resources could disrupt any common community phylogeographic pattern. We examined these alternative models using two Megastigmus parasitoid morphospecies associated with oak cynipid galls sampled throughout their Western Palaearctic distributions. Based on existing host cynipid data, a parallel evolution model predicts that eastern regions of the Western Palaearctic should contain ancestral populations with range expansions across Europe about 1.6 million years ago and deeper species-level divergence at both 8-9 and 4-5 million years ago. Sequence data from mitochondrial cytochrome b and multiple nuclear genes showed similar phylogenetic patterns and revealed cryptic genetic species within both morphospecies, indicating greater diversity in these communities than previously thought. Phylogeographic divergence was apparent in most cryptic species between relatively stable, diverse, putatively ancestral populations in Asia Minor and the Middle East, and genetically depauperate, rapidly expanding populations in Europe, paralleling patterns in host gallwasp species. Mitochondrial and nuclear data also suggested that Europe may have been colonized multiple times from eastern source populations since the late Miocene. Temporal patterns of lineage divergence were congruent within and across trophic levels, supporting the Contemporary Host-tracking Hypothesis for community evolution.

The widespread utilization of molecular markers has revealed that a broad spectrum of taxa contain sets of morphologically cryptic, but genetically distinct lineages (Bickford et al. 2007). The identification of cryptic taxa is important as an accurate appreciation of diversity is crucial for a proper understanding of evolutionary and ecological processes. An example is the study of host specificity in parasitic taxa, where an apparent generalist may be found to contain a complex of several more specific species (Smith et al. 2006). Host specificity is a key life history trait that varies greatly among parasites (Poulin & Keeney 2007). While some can exploit a wide range of hosts, others are confined to just a single species. Access to additional hosts increases the resources available to a parasite. However, physiological or ecological constraints can restrict the extension of host range. Furthermore, there may be a trade‐off between relaxed specificity and performance: generalism can decrease a parasites ability to adapt to each individual host species, and increase exposure to competition from other parasites (Poulin 1998). Despite the central role that host specificity plays in parasite life history, relatively little is known about how host range is determined in natural systems, and data from field studies are required to evaluate among competing ideas. In this issue, an exciting paper by Locke et al. (2010) makes a valuable contribution toward the understanding of host specificity in an important group of trematode flatworms. Using molecular methods, Locke et al. reveal an almost four‐fold increase in the appreciated diversity of their focal group. In combination with a large and elegant sampling design this allows them to accurately assess host specificity for each taxon, and thus draw key insights into the factors that control host range in a dominant parasite group.

AbstractOak gallwasps (Hymenoptera, Cynipidae, Cynipini) are one of seven major animal taxa that commonly reproduce by cyclical parthenogenesis (CP). A major question in research on CP taxa is the frequency with which lineages lose their sexual generations, and diversify as purely asexual radiations. Most oak gallwasp species are only known from an asexual generation, and secondary loss of sex has been conclusively demonstrated in several species, particularly members of the holarctic genus Andricus. This raises the possibility of widespread secondary loss of sex in the Cynipini, and of diversification within purely parthenogenetic lineages. We use two approaches based on analyses of allele frequency data to test for cryptic sexual generations in eight apparently asexual European species distributed through a major western palaearctic lineage of the gallwasp genus Andricus. All species showing adequate levels of polymorphism (7/8) showed signatures of sex compatible with cyclical parthenogenesis. We also use DNA sequence data to test the hypothesis that ignorance of these sexual generations (despite extensive study on this group) results from failure to discriminate among known but morphologically indistinguishable sexual generations. This hypothesis is supported: 35 sequences attributed by leading cynipid taxonomists to a single sexual adult morphospecies, Andricus burgundus, were found to represent the sexual generations of at least six Andricus species. We confirm cryptic sexual generations in a total of 11 Andricus species, suggesting that secondary loss of sex is rare in Andricus.

AbstractHuman dispersal of organisms is an important process modifying natural patterns of biodiversity. Such dispersal generates new patterns of genetic diversity that overlie natural phylogeographical signatures, allowing discrimination between alternative dispersal mechanisms. Here we use allele frequency and DNA sequence data to distinguish between alternative scenarios (unassisted range expansion and long range introduction) for the colonization of northern Europe by an oak‐feeding gallwasp, Andricus kollari. Native to Mediterranean latitudes from Portugal to Iran, this species became established in northern Europe following human introduction of a host plant, the Turkey oak Quercus cerris. Colonization of northern Europe is possible through three alternative routes: (i) unassisted range expansion from natural populations in the Iberian Peninsula; (ii) unassisted range expansion from natural populations in Italy and Hungary; or (iii) descent from populations imported to the UK as trade goods from the eastern Mediterranean in the 1830s. We show that while populations in France were colonized from sources in Italy and Hungary, populations in the UK and neighbouring parts of coastal northern Europe encompass allozyme and sequence variation absent from the known native range. Further, these populations show demographic signatures expected for large stable populations, rather than signatures of rapid population growth from small numbers of founders. The extent and spatial distribution of genetic diversity in the UK suggests that these A. kollari populations are derived from introductions of large numbers of individuals from each of two genetically divergent centres of diversity in the eastern Mediterranean. The strong spatial patterning in genetic diversity observed between different regions of northern Europe, and between sites in the UK, is compatible with leptokurtic models of population establishment.

AbstractInsect parasitoids are important components of many terrestrial ecosystems. However, relatively little is known about the mechanisms responsible for structuring their populations. Here we investigate the ability of Megastigmus stigmatizans, an oak gall wasp parasitoid, to track its host Andricus kollari over two different timescales, and examine its current population structure across a divide in host population structure. The divide represents a transition in gall wasp host‐plant species and offers the opportunity to examine whether the split, which divides gall wasp populations, manifests itself in the next trophic level. Analysis of mitochondrial haplotype data for parasitoid and host reveals: (i) a similar phylogeographic population structure for both, with Iberian populations more derived with respect to more eastern populations. (ii) it is likely that the host colonized the Iberian refuge earlier than the parasitoid, probably by at least one glacial cycle. (iii) Recent range expansion of central European host populations northwards has resulted in pursuit by parasitoids from the same geographic origin. (iv) in addition, Iberian parasitoid populations have crossed a major divide in host population structure to invade northern Europe. Such human‐facilitated escape from natural refugial distributions may have important implications for the composition and structure of northern European gall wasp communities.