Quality Of Genome Assemblies Research Articles

The Impact of Oxford Nanopore Technologies Based Methodologies on the Genome Sequencing and Assembly of Romanian Strains of Drosophila suzukii.

Drosophila suzukii is a worldwide invasive species with serious economic impacts. Herein, we are presenting the first project of sequencing and assembling the whole genomes of two lines of D. suzukii derived from Romanian local populations using exclusively Oxford Nanopore Technologies data. We implemented both MinION and Flongle flow-cells and tested the impact of various basecalling models and assembly strategies on the quality of the sought-after representative genome assemblies. We demonstrate that the sup-basecalling model significantly improved the read quality and that adding a relatively small collection of reads had a significant positive impact over the assembly quality. The novel dScaff bioinformatics prototype tool allowed us to perform sequence-level quality tests, as well as to represent assembly selections and display both the contig redundancy and the repeats-enriched genomic sub-sequences. Moreover, we used dScaff to propose a minimal assembly variant corresponding to one of our lines, GB-ls-coga4, which assured a basic linear coverage of the genome and exhibited quality parameters comparable with those particular to the current reference genome assembly. The study presents the first sequencing and assembly of a D. suzukii line in Romania and argues the efficiency of long-read sequencing strategies.

Chromosome-level genome assembly and annotation of a sea toad (Chaunax sp.)

The sea toad genus Chaunax is a group of small benthic fishes that predominantly inhabiting the deep seas of the Atlantic, Indian, and Pacific Oceans. Although they have the potential to make excellent systems for studies of evolutionary adaptation to deep-sea environments, genomic research on Chaunax has been hindered by a scarcity of high-quality genomic resources. We present a chromosome-scale genome assembly of a Chaunax specimen generated using PacBio long-read sequencing and high-throughput chromosome conformation capture technology. The size of the assembled genome was 706.94 Mb, with a contig N50 of 15.24 Mb and scaffold N50 of 29.42 Mb. Approximately 96.11% of assembled sequences were anchored and oriented onto 24 pseudo-chromosomes. The genome contained 213.47 Mb repetitive sequences, 25,280 protein-coding genes, and 5,090 non-coding RNAs. The high ratio of complete BUSCO genes (97.20%) indicates high quality of genome assembly. The chromosomal-level reference genome of Chaunax sp. provides a preliminary molecular basis for understanding deep-sea adaptation and phenotypic evolution as well as an important reference for whole-genome sequencing of related species.

A new genome sequence resource for five invasive fruit flies of agricultural concern: Ceratitis capitata, C. quilicii, C. rosa, Zeugodacus cucurbitae and Bactrocera zonata (Diptera, Tephritidae)

Here, we present novel high quality genome assemblies for five invasive tephritid species of agricultural concern: Ceratitis capitata, C. quilicii, C. rosa, Zeugodacus cucurbitae and Bactrocera zonata (read depths between 65 and 78x). Three assemblies (C. capitata, C. quilicii and Z. cucurbitae) were scaffolded with chromosome conformation data and annotated using RNAseq reads. For some species this is the first reference genome available (B. zonata, C. quilicii and C. rosa), for others we have published improved annotated genomes (C. capitata and Z. cucurbitae). Together, the new references provide an important resource to advance research on genetic techniques for population control, develop rapid species identification methods, and explore eco-evolutionary studies.

A reference genome for Trichogramma kaykai: A tiny desert-dwelling parasitoid wasp with competing sex-ratio distorters.

The tiny parasitoid wasp Trichogramma kaykai inhabits the Mojave Desert of the southwest United States. Populations of this tiny insect variably host up to two different sex-distorting genetic elements: (1) the endosymbiotic bacterium Wolbachia which induces the parthenogenetic reproduction of females, and (2) a B-chromosome, "Paternal Sex Ratio" (PSR), which converts would-be female offspring to PSR-transmitting males. We report here the genome of a Wolbachia-infected Trichogramma kaykai isofemale colony KSX58. Using Oxford Nanopore sequencing we produced a final genome assembly of 203 Mbp with 45x coverage, consisting of 213 contigs with an N50 of 1.9 Mbp. The assembly is quite complete, with 91.41% complete BUSCOs recovered: a very high score for Trichogrammatids that have been previously characterized for having high levels of core gene losses. We also report a complete mitochondrial genome for T. kaykai, and an assembly of the associated Wolbachia, strain wTkk. We identified copies of the parthenogenesis-inducing genes pifA and pifB in a remnant prophage region of the wTkk genome. The Trichogramma kaykai assembly is the highest quality genome assembly for the genus to-date and will serve as a great resource for understanding the evolution of sex and selfish genetic elements.

Physiological and Molecular Mechanisms of Lepidopteran Insects: Genomic Insights and Applications of Genome Editing for Future Research.

Lepidopteran insects are a major threat to global agriculture, causing significant crop losses and economic damage. Traditional pest control methods are becoming less effective due to the rapid evolution of insecticide resistance. This study explores the current status and genomic characteristics of 1315 Lepidopteran records, alongside an overview of relevant research, utilizing advanced functional genomics techniques, including RNA-seq and CRISPR/Cas9 gene-editing technologies to uncover the molecular mechanisms underlying insecticide resistance. Our genomic analysis revealed significant variability in genome size, assembly quality, and chromosome number, which may influence species' biology and resistance mechanisms. We identified key resistance-associated genes and pathways, including detoxification and metabolic pathways, which help these insects evade chemical control. By employing CRISPR/Cas9 gene-editing techniques, we directly manipulated resistance-associated genes to confirm their roles in resistance, demonstrating their potential for targeted interventions in pest management. These findings emphasize the value of integrating genomic data into the development of effective and sustainable pest control strategies, reducing reliance on chemical insecticides and promoting environmentally friendly integrated pest management (IPM) approaches. Our study highlights the critical role of functional genomics in IPM and its potential to provide long-term solutions to the growing challenge of Lepidopteran resistance.

Leveraging the power of long reads for targeted sequencing.

Long-read sequencing technologies have improved the contiguity and, as a result, the quality of genome assemblies by generating reads long enough to span and resolve complex or repetitive regions of the genome. Several groups have shown the power of long reads in detecting thousands of genomic and epigenomic features that were previously missed by short-read sequencing approaches. While these studies demonstrate how long reads can help resolve repetitive and complex regions of the genome, they also highlight the throughput and coverage requirements needed to accurately resolve variant alleles across large populations using these platforms. At the time of this review, whole-genome long-read sequencing is more expensive than short-read sequencing on the highest throughput short-read instruments; thus, achieving sufficient coverage to detect low-frequency variants (such as somatic variation) in heterogenous samples remains challenging. Targeted sequencing, on the other hand, provides the depth necessary to detect these low-frequency variants in heterogeneous populations. Here, we review currently used and recently developed targeted sequencing strategies that leverage existing long-read technologies to increase the resolution with which we can look at nucleic acids in a variety of biological contexts.

Chromosome-level genome assembly of the sacoglossan sea slug Elysia timida (Risso, 1818)

BackgroundSequencing and annotating genomes of non-model organisms helps to understand genome architecture, the genetic processes underlying species traits, and how these genes have evolved in closely-related taxa, among many other biological processes. However, many metazoan groups, such as the extremely diverse molluscs, are still underrepresented in the number of sequenced and annotated genomes. Although sequencing techniques have recently improved in quality and quantity, molluscs are still neglected due to difficulties in applying standardized protocols for obtaining genomic data.ResultsIn this study, we present the chromosome-level genome assembly and annotation of the sacoglossan sea slug species Elysia timida, known for its ability to store the chloroplasts of its food algae. In particular, by optimizing the long-read and chromosome conformation capture library preparations, the genome assembly was performed using PacBio HiFi and Arima HiC data. The scaffold and contig N50s, at 41.8 Mb and 1.92 Mb, respectively, are approximately 30-fold and fourfold higher compared to other published sacoglossan genome assemblies. Structural annotation resulted in 19,904 protein-coding genes, which are more contiguous and complete compared to publicly available annotations of Sacoglossa with respect to metazoan BUSCOs. We found no evidence for horizontal gene transfer (HGT), i.e. no photosynthetic genes encoded in the sacoglossan nucleus genome. However, we detected genes encoding polyketide synthases in E. timida, indicating that polypropionates are produced. HPLC–MS/MS analysis confirmed the presence of a large number of polypropionates, including known and yet uncharacterised compounds.ConclusionsWe can show that our methodological approach helps to obtain a high-quality genome assembly even for a "difficult-to-sequence" organism, which may facilitate genome sequencing in molluscs. This will enable a better understanding of complex biological processes in molluscs, such as functional kleptoplasty in Sacoglossa, by significantly improving the quality of genome assemblies and annotations.

Comparative analysis of amphibian genomes: An emerging resource for basic and applied research.

Amphibians are the most threatened group of vertebrates and are in dire need of conservation intervention to ensure their continued survival. They exhibit unique features including a high diversity of reproductive strategies, permeable and specialized skin capable of producing toxins and antimicrobial compounds, multiple genetic mechanisms of sex determination and in some lineages, the ability to regenerate limbs and organs. Although genomic approaches would shed light on these unique traits and aid conservation, sequencing and assembly of amphibian genomes has lagged behind other taxa due to their comparatively large genome sizes. Fortunately, the development of long-read sequencing technologies and initiatives has led to a recent burst of new amphibian genome assemblies. Although growing, the field of amphibian genomics suffers from the lack of annotation resources, tools for working with challenging genomes and lack of high-quality assemblies in multiple clades of amphibians. Here, we analyse 51 publicly available amphibian genomes to evaluate their usefulness for functional genomics research. We report considerable variation in genome assembly quality and completeness and report some of the highest transposable element and repeat contents of any vertebrate. Additionally, we detected an association between transposable element content and climatic variables. Our analysis provides evidence of conserved genome synteny despite the long divergence times of this group, but we also highlight inconsistencies in chromosome naming and orientation across genome assemblies. We discuss sequencing gaps in the phylogeny and suggest key targets for future sequencing endeavours. Finally, we propose increased investment in amphibian genomics research to promote their conservation.

Signatures of Mollicutes-related endobacteria in publicly available Mucoromycota genomes.

Mucoromycota fungi and their Mollicutes-related endobacteria (MRE) are an ideal system for studying bacterial-fungal interactions and evolution due to the long-term and intimate nature of their interactions. However, methods for detecting MRE face specific challenges due to the poor representation of MRE in sequencing databases coupled with the high sequence divergence of their genomes, making traditional similarity searches unreliable. This has precluded estimations on the diversity of MRE associated with Mucoromycota. To determine the prevalence of previously undetected MRE in fungal genome sequences, we scanned 389 Mucoromycota genome assemblies available from the National Center for Biotechnology Information for the presence of MRE sequences using publicly available tools to map contigs from fungal assemblies to publicly available MRE genomes. We demonstrate a higher diversity of MRE genomes than previously described in Mucoromycota and a lack of cophylogeny between MRE and the majority of their fungal hosts. This supports the late invasion hypothesis regarding MRE acquisition across most of the examined fungal families. In contrast with other Mucoromycota lineages, MRE from the Gigasporaceae displayed some degree of cophylogeny with their hosts, which may indicate that horizontal transmission is restricted between members of this family or that transmission is strictly vertical. These results underscore the need for a refined process to capture sequencing data from potential fungal endosymbionts to discern their evolution and transmission. Screens of fungal genomes for MRE can help improve the quality of fungal genome assemblies while identifying new MRE lineages to further test hypotheses on their origin and evolution.IMPORTANCEMollicutes-related endobacteria (MRE) are obligate intracellular bacteria found within Mucoromycota fungi. Despite their frequent detection, MRE roles in host functioning are still unknown. Comparative genomic investigations can improve our understanding of the impact of MRE on their fungal hosts by identifying similarities and differences in MRE genome evolution. However, MRE genomes have only been assembled from a small fraction of Mucoromycota hosts. Here, we demonstrate that MRE can be present yet undetected in publicly available Mucoromycota genome assemblies. We use these newfound sequences to assess the broader diversity of MRE and their phylogenetic relationships with respect to their hosts. We demonstrate that publicly available tools can be used to extract novel MRE sequences from assembled fungal genomes leading to insights on MRE evolution. This work contributes to a greater understanding of the fungal microbiome, which is crucial to improving knowledge on the dynamics and impacts of fungi in microbial ecosystems.

Chlomito: a novel tool for precise elimination of organelle genome contamination from nuclear genome assembly.

Accurate reference genomes are fundamental to understanding biological evolution, biodiversity, hereditary phenomena and diseases. However, many assembled nuclear chromosomes are often contaminated by organelle genomes, which will mislead bioinformatic analysis, and genomic and transcriptomic data interpretation. To address this issue, we developed a tool named Chlomito, aiming at precise identification and elimination of organelle genome contamination from nuclear genome assembly. Compared to conventional approaches, Chlomito utilized new metrics, alignment length coverage ratio (ALCR) and sequencing depth ratio (SDR), thereby effectively distinguishing true organelle genome sequences from those transferred into nuclear genomes via horizontal gene transfer (HGT). The accuracy of Chlomito was tested using sequencing data from Plum, Mango and Arabidopsis. The results confirmed that Chlomito can accurately detect contigs originating from the organelle genomes, and the identified contigs covered most regions of the organelle reference genomes, demonstrating efficiency and precision of Chlomito. Considering user convenience, we further packaged this method into a Docker image, simplified the data processing workflow. Overall, Chlomito provides an efficient, accurate and convenient method for identifying and removing contigs derived from organelle genomes in genomic assembly data, contributing to the improvement of genome assembly quality.

Genome sequencing and assembly of Indian major carp, Cirrhinus mrigala (Hamilton, 1822)

The freshwater aquaculture species, Cirrhinus mrigala, commonly known as mrigal, holds significance in the carp-dominated aquaculture system, globally. Despite constituting 1.08% of the total freshwater aquaculture production, mrigal is the third most important Indian major carp. However, its genome and associated information is not available. This study aims to address this gap by generation high quality genome assembly using PacBio long reads, Illumina short reads and Hi-C scaffolding. The characterization of assembled highly contiguous genome, 1.057 Gb in size, revealed 39,091 genes with functional annotations. The orthology analysis based on direct orthologs and single copy ortholgs places C. mrigala in a distinct position within the Otophysi clade. Additionally, the study delves into Hox gene clusters, identifying 38 Hox genes distributed in seven clusters. The present genomic information offers potential applications for sustainable aquaculture management, including selection programs for economic traits.

AssemblyQC: a Nextflow pipeline for reproducible reporting of assembly quality.

Genome assembly projects have grown exponentially due to breakthroughs in sequencing technologies and assembly algorithms. Evaluating the quality of genome assemblies is critical to ensure the reliability of downstream analysis and interpretation. To fulfil this task, we have developed the AssemblyQC pipeline that performs file-format validation, contaminant checking, contiguity measurement, gene- and repeat-space completeness quantification, telomere inspection, taxonomic assignment, synteny alignment, scaffold examination through Hi-C contact-map visualization, and assessments of completeness, consensus quality and phasing through k-mer analysis. It produces a comprehensive HTML report with method descriptions, tables, and visualizations. The pipeline uses Nextflow for workflow orchestration and adheres to the best-practice established by the nf-core community. This pipeline offers a reproducible, scalable, and portable method to assess the quality of genome assemblies-the code is available online at GitHub: https://github.com/Plant-Food-Research-Open/assemblyqc.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes

BackgroundThe honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses.ResultsPacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants.ConclusionsOur new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

The Chromosome-Scale Genome of Chitala ornata Illuminates the Evolution of Early Teleosts.

Teleosts are the most prolific vertebrates, occupying the vast majority of aquatic environments, and their pectoral fins have undergone remarkable physiological transformations throughout their evolution. Studying early teleost fishes, such as those belonging to the Osteoglossiformes order, could offer crucial insights into the adaptive evolution of pectoral fins within this group. In this study, we have assembled a chromosomal-level genome for the Clown featherback (Chitala ornata), achieving the highest quality genome assembly for Osteoglossiformes to date, with a contig N50 of 32.78 Mb and a scaffold N50 of 40.73 Mb. By combining phylogenetic analysis, we determined that the Clown featherback diverged approximately 202 to 203 million years ago (Ma), aligning with continental separation events. Our analysis revealed the intriguing discovery that a unique deletion of regulatory elements is adjacent to the Gli3 gene, specifically in teleosts. This deletion might be tied to the specialized adaptation of their pectoral fins. Furthermore, our findings indicate that specific contractions and expansions of transposable elements (TEs) in teleosts, including the Clown featherback, could be connected to their adaptive evolution. In essence, this study not only provides a high-quality genomic resource for Osteoglossiformes but also sheds light on the evolutionary trajectory of early teleosts.

Variability in Chromosome 1 of Select Moroccan Pyrenophora teres f. teres Isolates Overcomes a Highly Effective Barley Chromosome 6H Source of Resistance.

Barley net form net blotch (NFNB) is a destructive foliar disease caused by Pyrenophora teres f. teres. Barley line CIho5791, which harbors the broadly effective chromosome 6H resistance gene Rpt5, displays dominant resistance to P. teres f. teres. To genetically characterize P. teres f. teres avirulence/virulence on the barley line CIho5791, we generated a P. teres f. teres mapping population using a cross between the Moroccan CIho5791-virulent isolate MorSM40-3 and the avirulent reference isolate 0-1. Full genome sequences were generated for 103 progenies. Saturated chromosome-level genetic maps were generated, and quantitative trait locus (QTL) mapping identified two major QTL associated with P. teres f. teres avirulence/virulence on CIho5791. The most significant QTL mapped to chromosome (Ch) 1, where the virulent allele was contributed by MorSM40-3. A second QTL mapped to Ch8; however, this virulent allele was contributed by the avirulent parent 0-1. The Ch1 and Ch8 loci accounted for 27 and 15% of the disease variation, respectively, and the avirulent allele at the Ch1 locus was epistatic over the virulent allele at the Ch8 locus. As a validation, we used a natural P. teres f. teres population in a genome-wide association study that identified the same Ch1 and Ch8 loci. We then generated a new reference quality genome assembly of parental isolate MorSM40-3 with annotation supported by deep transcriptome sequencing of infection time points. The annotation identified candidate genes predicted to encode small, secreted proteins, one or more of which are likely responsible for overcoming the CIho5791 resistance. [Formula: see text] The author(s) have dedicated the work to the public domain under the Creative Commons CC0 "No Rights Reserved" license by waiving all of his or her rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law, 2024.

Assembly collapsing versus heterozygosity oversizing: detection of homokaryotic and heterokaryotic Laccaria trichodermophora strains by hybrid genome assembly.

Genome assembly and annotation using short-paired reads is challenging for eukaryotic organisms due to their large size, variable ploidy and large number of repetitive elements. However, the use of single-molecule long reads improves assembly quality (completeness and contiguity), but haplotype duplications still pose assembly challenges. To address the effect of read length on genome assembly quality, gene prediction and annotation, we compared genome assemblers and sequencing technologies with four strains of the ectomycorrhizal fungus Laccaria trichodermophora. By analysing the predicted repertoire of carbohydrate enzymes, we investigated the effects of assembly quality on functional inferences. Libraries were generated using three different sequencing platforms (Illumina Next-Seq, Mi-Seq and PacBio Sequel), and genomes were assembled using single and hybrid assemblies/libraries. Long reads or hybrid assemby resolved the collapsing of repeated regions, but the nuclear heterozygous versions remained unresolved. In dikaryotic fungi, each cell includes two nuclei and each nucleus has differences not only in allelic gene version but also in gene composition and synteny. These heterokaryotic cells produce fragmentation and size overestimation of the genome assembly of each nucleus. Hybrid assembly revealed a wider functional diversity of genomes. Here, several predicted oxidizing activities on glycosyl residues of oligosaccharides and several chitooligosaccharide acetylase activities would have passed unnoticed in short-read assemblies. Also, the size and fragmentation of the genome assembly, in combination with heterozygosity analysis, allowed us to distinguish homokaryotic and heterokaryotic strains isolated from L. trichodermophora fruit bodies.

High quality assemblies of four indigenous chicken genomes and related functional data resources

Many lines of evidence indicate that red jungle fowl (RJF) is the primary ancestor of domestic chickens. Although multiple versions of RJF (galgal2-galgal5 and GRCg6a) and commercial chickens (GRCg7b/w and Huxu) genomes have been assembled since 2004, no high-quality indigenous chicken genomes have been assembled, hampering the understanding of chicken domestication and evolution. To fill the gap, we sequenced the genomes of four indigenous chickens with distinct morphological traits in southwest China, using a combination of short, long and Hi-C reads. We assembled each genome (~1.0 Gb) into 42 chromosomes with chromosome N50 90.5–90.9 Mb, amongst the highest quality of chicken genome assemblies. To provide resources for gene annotation and functional analysis, we also sequenced transcriptomes of 10 tissues for each of the four chickens. Moreover, we corrected many mis-assemblies and assembled missing micro-chromosomes 29 and 34–39 for GRCg6a. Our assemblies, sequencing data and the correction of GRCg6a can be valuable resources for studying chicken domestication and evolution.

Development of Genome-Wide Intron Length Polymorphism (ILP) Markers in Tea Plant (Camellia sinensis) and Related Applications for Genetics Research.

The market value of tea is largely dependent on the tea species and cultivar. Therefore, it is important to develop efficient molecular markers covering the entire tea genome that can be used for the identification of tea varieties, marker-assisted breeding, and mapping important quantitative trait loci for beneficial traits. In this study, genome-wide molecular markers based on intron length polymorphism (ILP) were developed for tea trees. A total of 479, 1393, and 1342 tea ILP markers were identified using the PCR method in silico from the 'Shuchazao' scaffold genome, the chromosome-level genome of 'Longjing 43', and the ancient tea DASZ chromosome-level genome, respectively. A total of 230 tea ILP markers were used to amplify six tea tree species. Among these, 213 pairs of primers successfully characterize products in all six species, with 112 primer pairs exhibiting polymorphism. The polymorphism rate of primer pairs increased with the improvement in reference genome assembly quality level. The cross-species transferability analysis of 35 primer pairs of tea ILP markers showed an average amplification rate of 85.17% through 11 species in 6 families, with high transferability in Camellia reticulata and tobacco. We also used 40 pairs of tea ILP primers to evaluate the genetic diversity and population structure of C. tetracocca with 176 plants from Puan County, Guizhou Province, China. These genome-wide markers will be a valuable resource for genetic diversity analysis, marker-assisted breeding, and variety identification in tea, providing important information for the tea industry.

High quality genome assembly of the brown hare (Lepus europaeus) with chromosome-level scaffolding

High quality genome assembly of the brown hare (Lepus europaeus) with chromosome-level scaffolding

The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes.

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.