De Novo Whole-Genome Assembly of the Swede Midge (Contarinia nasturtii), a Specialist of Brassicaceae, Using Linked-Read Sequencing

The chromosome-level genome assembly of Aphidoletes aphidimyza Rondani (Diptera: Cecidomyiidae)

Aphidoletes aphidimyza is widely recognized as an effective predator of aphids in agricultural systems. However, there is limited understanding of its predation mechanisms. In this study, we generated a high-quality chromosome level of the A. aphidimyza genome by combining PacBio, Illumina, and Hi-C data. The genome has a size of 192.08 Mb, with a scaffold N50 size of 46.85 Mb, and 99.08% (190.35 Mb) of the assembly is located on four chromosomes. The BUSCO analysis of our assembly indicates a completeness of 97.8% (n = 1,367), including 1,307 (95.6%) single-copy BUSCOs and 30 (2.2%) duplicated BUSCOs. Additionally, we annotated a total of 13,073 protein-coding genes, 18.43% (35.40 Mb) repetitive elements, and 376 non-coding RNAs. Our study is the first time to report the chromosome-scale genome for the species of A. aphidimyza. It provides a valuable genomic resource for the molecular study of A. aphidimyza.

Exploring the contribution of the salivary gland and midgut to digestion in the swede midge (Contarinia nasturtii) through a genomics-guided approach

The larvae of Contarinia nasturtii (Kieffer) (Diptera: Cecidomyiidae), the swede midge, targets the meristem of brassica crops where they induce the formation of galls and disrupt seed and vegetable production. Previously, we examined the salivary gland transcriptome of newly-hatched first instar larvae as they penetrated the host and initiated gall formation. Here we examine the salivary gland and midgut transcriptome of third instar larvae and provide evidence for cooperative nutrient acquisition beginning with secretion of enzymes and feeding facilitators followed by gastrointestinal digestion. Sucrose, presumably obtained from the phloem, appeared to be a major nutrient source as several α-glucosidases (sucrases, maltases) and β-fructofuranosidases (invertases) were identified. Genes encoding β-fructofuranosidases/invertases were among the most highly expressed in both tissues and represented two distinct gene families that may have originated via horizontal gene transfer from bacteria. The importance of the phloem as a nutrient source is underscored by the expression of genes encoding regucalcin and ARMET (arginine-rich mutated in early stages of tumor) which interfere with calcium signalling and prevent sieve tube occlusion. Lipids, proteins, and starch appear to serve as a secondary nutrient sources. Genes encoding enzymes involved in the detoxification of glucosinolates (myrosinases, arylsulfatases, and glutathione-S-transferases) were expressed indicative of Brassicaceae host specialization. The midgut expressed simple peritrophins and mucins typical of those found in Type II peritrophic matrices, the first such description for a gall midge.

Characterization of the swede midge, Contarinia nasturtii, first instar larval salivary gland transcriptome

Proteins in saliva of gall-forming insect larvae govern insect-host plant interactions. Contarinia nasturtii, the swede midge, is a pest of brassicaceous vegetables (cabbage, cauliflower, broccoli) and canola. We examined the salivary gland (SG) transcriptome of first instar larvae reared on Brassica napus and catalogued genes encoding secreted proteins that may contribute to the initial stages of larval establishment, the synthesis of plant growth hormones, extra-oral digestion and evasion of host defenses. A significant portion of the secreted proteins with unknown functions were unique to C. nasturtii and were often members of larger gene families organized in genomic clusters with conservation patterns suggesting that they are undergoing selection.

The genome of the soybean gall midge (Resseliella maxima)

The cecidomyiid fly, soybean gall midge, Resseliella maxima Gagné, is a recently discovered insect that feeds on soybean plants in the Midwestern United States. R. maxima larvae feed on soybean stems that may induce plant death and can cause considerable yield losses, making it an important agricultural pest. From three pools of 50 adults each, we used long-read nanopore sequencing to assemble a R. maxima reference genome. The final genome assembly is 206 Mb with 64.88× coverage, consisting of 1,009 contigs with an N50 size of 714 kb. The assembly is high quality with a Benchmarking Universal Single-Copy Ortholog (BUSCO) score of 87.8%. Genome-wide GC level is 31.60%, and DNA methylation was measured at 1.07%. The R. maxima genome is comprised of 21.73% repetitive DNA, which is in line with other cecidomyiids. Protein prediction annotated 14,798 coding genes with 89.9% protein BUSCO score. Mitogenome analysis indicated that R. maxima assembly is a single circular contig of 15,301 bp and shares highest identity to the mitogenome of the Asian rice gall midge, Orseolia oryzae Wood-Mason. The R. maxima genome has one of the highest completeness levels for a cecidomyiid and will provide a resource for research focused on the biology, genetics, and evolution of cecidomyiids, as well as plant-insect interactions in this important agricultural pest.

The Genome of the Soybean Gall Midge ( Resseliella maxima )

The cecidomyiid fly, soybean gall midge, Resseliella maxima Gagné, is a recently discovered insect that feeds on soybean plants in the Midwest US. Resseliella maxima larvae feed on soybean stems which may induce plant death and can cause considerable yield losses, making it an important agricultural pest. From three pools of 50 adults each, we used long-read nanopore sequencing to assemble a R. maxima reference genome. The final genome assembly is 206 Mb with 64.88X coverage, consisting of 1009 contigs with an N50 size of 714 kb. The assembly is high quality with a BUSCO score of 87.8%. Genome-wide GC level is 31.60% and DNA methylation was measured at 1.07%. The R. maxima genome is comprised of 21.73% repetitive DNA, which is in line with other cecidomyiids. Protein prediction annotated 14,798 coding genes with 89.9% protein BUSCO score. Mitogenome analysis indicated that R. maxima assembly is a single circular contig of 15,301 bp and shares highest identity to the mitogenome of the Asian rice gall midge, Orseolia oryzae (Wood-Mason). The R. maxima genome has one of the highest completeness levels for a cecidomyiid and will provide a resource for research focused on the biology, genetics, and evolution of cecidomyiids, as well as plant-insect interactions in this important agricultural pest.

A chromosome-level genome assembly of the orange wheat blossom midge, Sitodiplosis mosellana Géhin (Diptera: Cecidomyiidae) provides insights into the evolution of a detoxification system

The orange wheat blossom midge Sitodiplosis mosellana Géhin (Diptera: Cecidomyiidae), an economically important pest, has caused serious yield losses in most wheat-growing areas worldwide in the past half-century. A high-quality chromosome-level genome for S. mosellana was assembled using PacBio long read, Illumina short read, and Hi-C sequencing technologies. The final genome assembly was 180.69 Mb, with contig and scaffold N50 sizes of 998.71 kb and 44.56 Mb, respectively. Hi-C scaffolding reliably anchored 4 pseudochromosomes, accounting for 99.67% of the assembled genome. In total, 12,269 protein-coding genes were predicted, of which 91% were functionally annotated. Phylogenetic analysis indicated that S. mosellana and its close relative, the swede midge Contarinia nasturtii, diverged about 32.7 MYA. The S. mosellana genome showed high chromosomal synteny with the genome of Drosophila melanogaster and Anopheles gambiae. The key gene families involved in the detoxification of plant secondary chemistry were analyzed. The high-quality S. mosellana genome data will provide an invaluable resource for research in a broad range of areas, including the biology, ecology, genetics, and evolution of midges, as well as insect–plant interactions and coevolution.

Horizontal Transfer of Microbial Toxin Genes to Gall Midge Genomes

A growing body of evidence has underscored the role of horizontal gene transfer (HGT) in animal evolution. Previously, we discovered the horizontal transfer of the gene encoding the eukaryotic genotoxin cytolethal distending toxin B (cdtB) from the pea aphid Acyrthosiphon pisum secondary endosymbiont (APSE) phages to drosophilid and aphid nuclear genomes. Here, we report cdtB in the nuclear genome of the gall-forming "swede midge" Contarinia nasturtii (Diptera: Cecidomyiidae) via HGT. We searched all available gall midge genome sequences for evidence of APSE-to-insect HGT events and found five toxin genes (aip56, cdtB, lysozyme, rhs, and sltxB) transferred horizontally to cecidomyiid nuclear genomes. Surprisingly, phylogenetic analyses of HGT candidates indicated APSE phages were often not the ancestral donor lineage of the toxin gene to cecidomyiids. We used a phylogenetic signal statistic to test a transfer-by-proximity hypothesis for animal HGT, which suggested that microbe-to-insect HGT was more likely between taxa that share environments than those from different environments. Many of the toxins we found in midge genomes target eukaryotic cells, and catalytic residues important for toxin function are conserved in insect copies. This class of horizontally transferred, eukaryotic cell-targeting genes is potentially important in insect adaptation.