Conditional knockdown of transformer in sheep blow fly suggests a role in repression of dosage compensation and potential for population suppression.

The transformer (tra) gene is essential for female development in many insect species, including the Australian sheep blow fly, Lucilia cuprina. Sex-specific tra RNA splicing is controlled by Sex lethal (Sxl) in Drosophila melanogaster but is auto-regulated in L. cuprina. Sxl also represses X chromosome dosage compensation in female D. melanogaster. We have developed conditional Lctra RNAi knockdown strains using the tet-off system. Four strains did not produce females on diet without tetracycline and could potentially be used for genetic control of L. cuprina. In one strain, which showed both maternal and zygotic tTA expression, most XX transformed males died at the pupal stage. RNAseq and qRT-PCR analyses of mid-stage pupae showed increased expression of X-linked genes in XX individuals. These results suggest that Lctra promotes somatic sexual differentiation and inhibits X chromosome dosage compensation in female L. cuprina. However, XX flies homozygous for a loss-of-function Lctra knockin mutation were fully transformed and showed high pupal eclosion. Two of five X-linked genes examined showed a significant increase in mRNA levels in XX males. The stronger phenotype in the RNAi knockdown strain could indicate that maternal Lctra expression may be essential for initiation of dosage compensation suppression in female embryos.

Coordination between terminal variation of the viral genome and insect microRNAs regulates rice stripe virus replication in insect vectors.

Maintenance of a balance between the levels of viral replication and selective pressure from the immune systems of insect vectors is one of the prerequisites for efficient transmission of insect-borne propagative phytoviruses. The mechanism regulating the adaptation of RNA viruses to insect vectors by genomic variation remains unknown. Our previous study demonstrated an extension of the 3'-untranslated terminal region (UTR) of two genomic segments of rice stripe virus (RSV). In the present study, a reverse genetic system for RSV in human cells and an insect vector, the small brown planthopper Laodelphax striatellus, was used to demonstrate that the 3'-terminal extensions suppressed viral replication in vector insects by inhibiting promoter activity due to structural interference with the panhandle structure formed by viral 3'- and 5'-UTRs. The extension sequence in the viral RNA1 segment was targeted by an endogenous insect microRNA, miR-263a, which decreased the inhibitory effect of the extension sequence on viral promoter activity. Surprisingly, the expression of miR-263a was negatively regulated by RSV infection. This elaborate coordination between terminal variation of the viral genome and endogenous insect microRNAs controls RSV replication in planthopper, thus reflecting a distinct strategy of adaptation of phytoviruses to insect vectors.

Intraspecific variation in immune gene expression and heritable symbiont density.

Host genetic variation plays an important role in the structure and function of heritable microbial communities. Recent studies have shown that insects use immune mechanisms to regulate heritable symbionts. Here we test the hypothesis that variation in symbiont density among hosts is linked to intraspecific differences in the immune response to harboring symbionts. We show that pea aphids (Acyrthosiphon pisum) harboring the bacterial endosymbiont Regiella insecticola (but not all other species of symbionts) downregulate expression of key immune genes. We then functionally link immune expression with symbiont density using RNAi. The pea aphid species complex is comprised of multiple reproductively-isolated host plant-adapted populations. These 'biotypes' have distinct patterns of symbiont infections: for example, aphids from the Trifolium biotype are strongly associated with Regiella. Using RNAseq, we compare patterns of gene expression in response to Regiella in aphid genotypes from multiple biotypes, and we show that Trifolium aphids experience no downregulation of immune gene expression while hosting Regiella and harbor symbionts at lower densities. Using F1 hybrids between two biotypes, we find that symbiont density and immune gene expression are both intermediate in hybrids. We propose that in this system, Regiella symbionts are suppressing aphid immune mechanisms to increase their density, but that some hosts have adapted to prevent immune suppression in order to control symbiont numbers. This work therefore suggests that antagonistic coevolution can play a role in host-microbe interactions even when symbionts are transmitted vertically and provide a clear benefit to their hosts. The specific immune mechanisms that we find are downregulated in the presence of Regiella have been previously shown to combat pathogens in aphids, and thus this work also highlights the immune system's complex dual role in interacting with both beneficial and harmful microbes.

Many, but not all, lineage-specific genes can be explained by homology detection failure.

Genes for which homologs can be detected only in a limited group of evolutionarily related species, called "lineage-specific genes," are pervasive: Essentially every lineage has them, and they often comprise a sizable fraction of the group's total genes. Lineage-specific genes are often interpreted as "novel" genes, representing genetic novelty born anew within that lineage. Here, we develop a simple method to test an alternative null hypothesis: that lineage-specific genes do have homologs outside of the lineage that, even while evolving at a constant rate in a novelty-free manner, have merely become undetectable by search algorithms used to infer homology. We show that this null hypothesis is sufficient to explain the lack of detected homologs of a large number of lineage-specific genes in fungi and insects. However, we also find that a minority of lineage-specific genes in both clades are not well explained by this novelty-free model. The method provides a simple way of identifying which lineage-specific genes call for special explanations beyond homology detection failure, highlighting them as interesting candidates for further study.

FlyBase: updates to the Drosophila melanogaster knowledge base.

FlyBase (flybase.org) is an essential online database for researchers using Drosophila melanogaster as a model organism, facilitating access to a diverse array of information that includes genetic, molecular, genomic and reagent resources. Here, we describe the introduction of several new features at FlyBase, including Pathway Reports, paralog information, disease models based on orthology, customizable tables within reports and overview displays ('ribbons') of expression and disease data. We also describe a variety of recent important updates, including incorporation of a developmental proteome, upgrades to the GAL4 search tab, additional Experimental Tool Reports, migration to JBrowse for genome browsing and improvements to batch queries/downloads and the Fast-Track Your Paper tool.

Maternal effect killing by a supergene controlling ant social organization.

Supergenes underlie striking polymorphisms in nature, yet the evolutionary mechanisms by which they arise and persist remain enigmatic. These clusters of linked loci can spread in populations because they captured coadapted alleles or by selfishly distorting the laws of Mendelian inheritance. Here, we show that the supergene haplotype associated with multiple-queen colonies in Alpine silver ants is a maternal effect killer. All eggs from heterozygous queens failed to hatch when they did not inherit this haplotype. Hence, the haplotype specific to multiple-queen colonies is a selfish genetic element that enhances its own transmission by causing developmental arrest of progeny that do not carry it. At the population level, such transmission ratio distortion favors the spread of multiple-queen colonies, to the detriment of the alternative haplotype associated with single-queen colonies. Hence, selfish gene drive by one haplotype will impact the evolutionary dynamics of alternative forms of colony social organization. This killer hidden in a social supergene shows that large nonrecombining genomic regions are prone to cause multifarious effects across levels of biological organization.

The gene cortex controls mimicry and crypsis in butterflies and moths.

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

The industrial melanism mutation in British peppered moths is a transposable element.

Discovering the mutational events that fuel adaptation to environmental change remains an important challenge for evolutionary biology. The classroom example of a visible evolutionary response is industrial melanism in the peppered moth (Biston betularia): the replacement, during the Industrial Revolution, of the common pale typica form by a previously unknown black (carbonaria) form, driven by the interaction between bird predation and coal pollution. The carbonaria locus has been coarsely localized to a 200-kilobase region, but the specific identity and nature of the sequence difference controlling the carbonaria-typica polymorphism, and the gene it influences, are unknown. Here we show that the mutation event giving rise to industrial melanism in Britain was the insertion of a large, tandemly repeated, transposable element into the first intron of the gene cortex. Statistical inference based on the distribution of recombined carbonaria haplotypes indicates that this transposition event occurred around 1819, consistent with the historical record. We have begun to dissect the mode of action of the carbonaria transposable element by showing that it increases the abundance of a cortex transcript, the protein product of which plays an important role in cell-cycle regulation, during early wing disc development. Our findings fill a substantial knowledge gap in the iconic example of microevolutionary change, adding a further layer of insight into the mechanism of adaptation in response to natural selection. The discovery that the mutation itself is a transposable element will stimulate further debate about the importance of 'jumping genes' as a source of major phenotypic novelty.

Transcriptomic basis and evolution of the ant nurse-larval social interactome.

Development is often strongly regulated by interactions among close relatives, but the underlying molecular mechanisms are largely unknown. In eusocial insects, interactions between caregiving worker nurses and larvae regulate larval development and resultant adult phenotypes. Here, we begin to characterize the social interactome regulating ant larval development by collecting and sequencing the transcriptomes of interacting nurses and larvae across time. We find that the majority of nurse and larval transcriptomes exhibit parallel expression dynamics across larval development. We leverage this widespread nurse-larva gene co-expression to infer putative social gene regulatory networks acting between nurses and larvae. Genes with the strongest inferred social effects tend to be peripheral elements of within-tissue regulatory networks and are often known to encode secreted proteins. This includes interesting candidates such as the nurse-expressed giant-lens, which may influence larval epidermal growth factor signaling, a pathway known to influence various aspects of insect development. Finally, we find that genes with the strongest signatures of social regulation tend to experience relaxed selective constraint and are evolutionarily young. Overall, our study provides a first glimpse into the molecular and evolutionary features of the social mechanisms that regulate all aspects of social life.

Functional genetic validation of key genes conferring insecticide resistance in the major African malaria vector, Anopheles gambiae.

Resistance in Anopheles gambiae to members of all 4 major classes (pyrethroids, carbamates, organochlorines, and organophosphates) of public health insecticides limits effective control of malaria transmission in Africa. Increase in expression of detoxifying enzymes has been associated with insecticide resistance, but their direct functional validation in An. gambiae is still lacking. Here, we perform transgenic analysis using the GAL4/UAS system to examine insecticide resistance phenotypes conferred by increased expression of the 3 genes-Cyp6m2, Cyp6p3, and Gste2-most often found up-regulated in resistant An. gambiae We report evidence in An. gambiae that organophosphate and organochlorine resistance is conferred by overexpression of GSTE2 in a broad tissue profile. Pyrethroid and carbamate resistance is bestowed by similar Cyp6p3 overexpression, and Cyp6m2 confers only pyrethroid resistance when overexpressed in the same tissues. Conversely, such Cyp6m2 overexpression increases susceptibility to the organophosphate malathion, presumably due to conversion to the more toxic metabolite, malaoxon. No resistant phenotypes are conferred when either Cyp6 gene overexpression is restricted to the midgut or oenocytes, indicating that neither tissue is involved in insecticide resistance mediated by the candidate P450s examined. Validation of genes conferring resistance provides markers to guide control strategies, and the observed negative cross-resistance due to Cyp6m2 gives credence to proposed dual-insecticide strategies to overcome pyrethroid resistance. These transgenic An. gambiae-resistant lines are being used to test the "resistance-breaking" efficacy of active compounds early in their development.

HapSolo: an optimization approach for removing secondary haplotigs during diploid genome assembly and scaffolding.

BACKGROUND: Despite marked recent improvements in long-read sequencing technology, the assembly of diploid genomes remains a difficult task. A major obstacle is distinguishing between alternative contigs that represent highly heterozygous regions. If primary and secondary contigs are not properly identified, the primary assembly will overrepresent both the size and complexity of the genome, which complicates downstream analysis such as scaffolding. RESULTS: Here we illustrate a new method, which we call HapSolo, that identifies secondary contigs and defines a primary assembly based on multiple pairwise contig alignment metrics. HapSolo evaluates candidate primary assemblies using BUSCO scores and then distinguishes among candidate assemblies using a cost function. The cost function can be defined by the user but by default considers the number of missing, duplicated and single BUSCO genes within the assembly. HapSolo performs hill climbing to minimize cost over thousands of candidate assemblies. We illustrate the performance of HapSolo on genome data from three species: the Chardonnay grape (Vitis vinifera), with a genome of 490 Mb, a mosquito (Anopheles funestus; 200 Mb) and the Thorny Skate (Amblyraja radiata; 2650 Mb). CONCLUSIONS: HapSolo rapidly identified candidate assemblies that yield improvements in assembly metrics, including decreased genome size and improved N50 scores. Contig N50 scores improved by 35%, 9% and 9% for Chardonnay, mosquito and the thorny skate, respectively, relative to unreduced primary assemblies. The benefits of HapSolo were amplified by down-stream analyses, which we illustrated by scaffolding with Hi-C data. We found, for example, that prior to the application of HapSolo, only 52% of the Chardonnay genome was captured in the largest 19 scaffolds, corresponding to the number of chromosomes. After the application of HapSolo, this value increased to ~ 84%. The improvements for the mosquito's largest three scaffolds, representing the number of chromosomes, were from 61 to 86%, and the improvement was even more pronounced for thorny skate. We compared the scaffolding results to assemblies that were based on PurgeDups for identifying secondary contigs, with generally superior results for HapSolo.

Space colonization by branching trachea explains the morphospace of a simple respiratory organ.

Branching morphogenesis helps increase the efficiency of gas and liquid transport in many animal organs. Studies in several model organisms have highlighted the molecular and cellular complexity behind branching morphogenesis. To understand this complexity, computational models have been developed with the goal of identifying the "major rules" that globally explain the branching patterns. These models also guide further experimental exploration of the biological processes that execute and maintain these rules. In this paper we introduce the tracheal gills of mayfly (Ephemeroptera) larvae as a model system to study the generation of branched respiratory patterns. First, we describe the gills of the mayfly Cloeon dipterum, and quantitatively characterize the geometry of its branching trachea. We next extend this characterization to those of related species to generate the morphospace of branching patterns. Then, we show how an algorithm based on the "space colonization" concept (SCA) can generate this branching morphospace via growth towards a hypothetical attractor molecule (M). SCA differs from other branch-generating algorithms in that the geometry generated depends to a great extent on its perception of the "external" space available for branching, uses few rules and, importantly, can be easily translated into a realistic "biological patterning algorithm". We identified a gene in the C. dipterum genome (Cd-bnl) that is orthologous to the fibroblast growth factor branchless (bnl), which stimulates growth and branching of embryonic trachea in Drosophila. In C. dipterum, this gene is expressed in the gill margins and areas of finer tracheolar branching from thicker trachea. Thus, Cd-bnl may perform the function of M in our model. Finally, we discuss this general mechanism in the context of other branching pattern-generating algorithms.

even-skipped acts as a pair-rule gene in germ band stages of Tribolium development.

The pair-rule gene even-skipped (eve) is essential for insect segmentation, yet its function varies among insect clades. While loss of eve results in typical pair-rule phenotypes in Drosophila, knock-down of eve orthologs shows segmental, gap-like, or asegmental phenotypes in non-Drosophila insects. In Tribolium, knock-down of the eve ortholog (Tc-eve) resulted in a graded phenotypic series ranging from strong to weak, the most informative of which was intermediate phenotypes. The strong knock-down embryos displayed asegmental phenotypes and severely disorganized germ bands which have prevented determination of Tc-eve function in later stages. In order to understand the segmentation function of Tc-eve during later germ band elongation stages, we analyzed intermediate Tc-eveRNAi embryos in which germ band elongation was less affected. Most intermediate Tc-eveRNAi germ bands displayed segmentation defects with a double segmental periodicity in the abdomen. In these intermediate embryos, Tc-engrailed (Tc-en) stripes were ectopically expanded into large bands with a double segmental periodicity, while the remaining Tc-en stripes between the expanded Tc-en stripes were absent or barely formed. The expanded Tc-en stripes seemed to be activated by primary Tc-eve stripes and Tc-paired, both of which failed to resolve into secondary segmental stripes. The absence of Tc-en stripes appeared to be a consequence of the absence of the secondary stripes of Tc-runt that were required for the activation of Tc-en stripes. These results suggest that Tc-eve functions as a pair-rule gene at least in the germ band stages of Tribolium development.

Wnt gene regulation and function during maxillary palp development in Drosophila melanogaster.

Wnt genes encode secreted ligands that play many important roles in the development of metazoans. There are thirteen known Wnt gene subfamilies and seven of these are represented in Drosophila melanogaster. While wingless (wg) is the best understood and most widely studied Wnt gene in Drosophila, the functions of many of the other Drosophila Wnt genes are less well understood. For example, relatively little is known about Wnt6, which is an ancient paralog of wg and they form a conserved Wnt cluster together with Wnt9 (Dwnt4) and Wnt10. Wg and Wnt6 encode similar proteins and exhibit overlapping expression in several tissues during development. Both wg and Wnt6 were previously shown to regulate the development of maxillary palps, important olfactory organs in flies, but it remained unclear how these two ligands may combine to carry out specific functions and how this is regulated. Here, we have further analysed Wnt6 function in the context of maxillary palp development. Surprisingly, we found that Wnt6 does not appear to be necessary for development of maxillary palps. While a deletion of the 5' region of Wnt6 results in very small maxillary palps, we show that this effect is more likely to be a consequence of removing cis-regulatory elements that may regulate wg expression in this tissue rather than through the loss of Wnt6 function. Although, we cannot completely exclude the possibility that Wnt6 may subtly regulate maxillary palp development in combination with wg, our analysis of Wnt6 loss of function mutants suggests this ligand plays a more general role in regulating growth during development. Taken together our results provide new insights into maxillary palp formation and Wnt6 functions in Drosophila, and further evidence for a complex cis-regulatory landscape in the Wnt9-wg-Wnt6-Wnt10 cluster, which may help explain its evolutionary conservation.

Mutation (G275E) of nAChR subunit Foα6 associated with spinetoram resistance in Australian western flower thrips, Frankliniella occidentalis (Pergande).

Western flower thrips, Frankliniella occidentalis is an economically important agricultural pest. It causes damage by feeding and oviposition or indirectly by plant virus transmission. Australian F. occidentalis are resistant to many insecticides including spinosad and the related chemical spinetoram. Spinetoram resistance to F. occidentalis has been recently reported in three different Australian States, however, mechanisms conferring that resistance have not been investigated. To identify the mechanisms underlying resistance to spinetoram in F. occidentalis, we sequenced the genomic region of nicotinic acetylcholine receptor Foα6 in number of spinosad and spinetoram resistant field-populations. We found that a single nucleotide substitution (G to A) in exon 9 of the α6 subunit was present in resistant strains (G275E) and absent from susceptible. By examining field populations we consider the G275E mutation is the major cause of resistance to spinetoram in Australian F. occidentalis. We developed a real-time PCR diagnostic assay to quickly identify resistant alleles in field-populations. The method was used to test spinetoram resistant F. occidentalis collected from Australian cotton during the 2018-2019. Results show thrips tested carried the G275E mutation and the resistance allele was unusually widely distributed. The wide distribution of G275E mutation was not expected because spinetoram is not extensively used in Australian cotton. We speculate resistance may relate to extensive chemical use in crops nearby such as horticulture where thrips are often targeted for control. Our molecular diagnostic assay can provide timely and precise resistance frequency information that can support sustainable chemical use including spinetoram based IPM.

Genome-wide identification of zero nucleotide recursive splicing in Drosophila.

Recursive splicing is a process in which large introns are removed in multiple steps by re-splicing at ratchet points--5' splice sites recreated after splicing. Recursive splicing was first identified in the Drosophila Ultrabithorax (Ubx) gene and only three additional Drosophila genes have since been experimentally shown to undergo recursive splicing. Here we identify 197 zero nucleotide exon ratchet points in 130 introns of 115 Drosophila genes from total RNA sequencing data generated from developmental time points, dissected tissues and cultured cells. The sequential nature of recursive splicing was confirmed by identification of lariat introns generated by splicing to and from the ratchet points. We also show that recursive splicing is a constitutive process, that depletion of U2AF inhibits recursive splicing, and that the sequence and function of ratchet points are evolutionarily conserved in Drosophila. Finally, we identify four recursively spliced human genes, one of which is also recursively spliced in Drosophila. Together, these results indicate that recursive splicing is commonly used in Drosophila, occurs in humans, and provides insight into the mechanisms by which some large introns are removed.

Consequences of resistance evolution in a Cas9-based sex conversion-suppression gene drive for insect pest management.

The use of a site-specific homing-based gene drive for insect pest control has long been discussed, but the easy design of such systems has become possible only with the recent establishment of CRISPR/Cas9 technology. In this respect, novel targets for insect pest management are provided by new discoveries regarding sex determination. Here, we present a model for a suppression gene drive designed to cause an all-male population collapse in an agricultural pest insect. To evaluate the molecular details of such a sex conversion-based suppression gene drive experimentally, we implemented this strategy in Drosophila melanogaster to serve as a safe model organism. We generated a Cas9-based homing gene-drive element targeting the transformer gene and showed its high efficiency for sex conversion from females to males. However, nonhomologous end joining increased the rate of mutagenesis at the target site, which resulted in the emergence of drive-resistant alleles and therefore curbed the gene drive. This confirms previous studies that simple homing CRISPR/Cas9 gene-drive designs will be ineffective. Nevertheless, by performing population dynamics simulations using the parameters we obtained in D. melanogaster and by adjusting the model for the agricultural pest Ceratitis capitata, we were able to identify adequate modifications that could be successfully applied for the management of wild Mediterranean fruit fly populations using our proposed sex conversion-based suppression gene-drive strategy.

TGFß signaling related genes are involved in hormonal mediation during termite soldier differentiation.

A working knowledge of the proximate factors intrinsic to sterile caste differentiation is necessary to understand the evolution of eusocial insects. Genomic and transcriptomic analyses in social hymenopteran insects have resulted in the hypothesis that sterile castes are generated by the novel function of co-opted or recruited universal gene networks found in solitary ancestors. However, transcriptome analysis during caste differentiation has not been tested in termites, and evolutionary processes associated with acquiring the caste are still unknown. Termites possess the soldier caste, which is regarded as the first acquired permanently sterile caste in the taxon. In this study, we performed a comparative transcriptome analysis in termite heads during 3 molting processes, i.e., worker, presoldier and soldier molts, under natural conditions in an incipient colony of the damp-wood termite Zootermopsis nevadensis. Although similar expression patterns were observed during each molting process, more than 50 genes were shown to be highly expressed before the presoldier (intermediate stage of soldier) molt. We then performed RNA interference (RNAi) of the candidate 13 genes, including transcription factors and uncharacterized protein genes, during presoldier differentiation induced by juvenile hormone (JH) analog treatment. Presoldiers induced after RNAi of two genes related to TGFß (Transforming growth factor beta) signaling were extremely unusual and possessed soldier-like phenotypes. These individuals also displayed aggressive behaviors similar to natural soldiers when confronted with Formica ants as hypothetical enemies. These presoldiers never molted into the next instar, presumably due to the decreased expression levels of the molting hormone (20-hydroxyecdysone; 20E) signaling genes. These results suggest that TGFß signaling was acquired for the novel function of regulating between JH and 20E signaling during soldier differentiation in termites.

Associated patterns of insecticide resistance in field populations of malaria vectors across Africa.

The development of insecticide resistance in African malaria vectors threatens the continued efficacy of important vector control methods that rely on a limited set of insecticides. To understand the operational significance of resistance we require quantitative information about levels of resistance in field populations to the suite of vector control insecticides. Estimation of resistance is complicated by the sparsity of observations in field populations, variation in resistance over time and space at local and regional scales, and cross-resistance between different insecticide types. Using observations of the prevalence of resistance in mosquito species from the Anopheles gambiae complex sampled from 1,183 locations throughout Africa, we applied Bayesian geostatistical models to quantify patterns of covariation in resistance phenotypes across different insecticides. For resistance to the three pyrethroids tested, deltamethrin, permethrin, and λ-cyhalothrin, we found consistent forms of covariation across sub-Saharan Africa and covariation between resistance to these pyrethroids and resistance to DDT. We found no evidence of resistance interactions between carbamate and organophosphate insecticides or between these insecticides and those from other classes. For pyrethroids and DDT we found significant associations between predicted mean resistance and the observed frequency of kdr mutations in the Vgsc gene in field mosquito samples, with DDT showing the strongest association. These results improve our capacity to understand and predict resistance patterns throughout Africa and can guide the development of monitoring strategies.

Analysis of Differentially Expressed Genes of Chrysoperla sinica Related to Flight Capacity by Transcriptome.

The lacewing Chrysoperla sinica (Tjeder) is a common natural enemy of many insect pests in China and is frequently employed for biological control programs. Adults make migratory flights after emergence, which reduces their effectiveness as biological control agents. Previously, we proved that 2-d-old unmated females exhibited significantly stronger flight ability than 3-d-old ones. Meanwhile, 3-d-old unmated adults flew significantly longer distances than mated ones. In this study, Illumina RNA sequencing was performed to characterize differentially expressed genes (DEGs) between virgin and mated adults of different ages in a single female strain of C. sinica. In total, 713,563,726 clean reads were obtained and de novo assembled into 109,165 unigenes with an average length of 847 bp (N50 of 1,754 bp), among which 4,382 (4.01%) unigenes matched known proteins. Based on these annotations, many putative transcripts were related to C. sinica's flight capacity and muscle structure, energy supply, growth, development, environmental adaptability, and metabolism of nutritional components and bioactive components. In addition, the differential expression of transcripts between different ages and mating status were analyzed, and DEGs participating in flight capacity and muscles were detected, including glutathione hydrolase, NAD-specific glutamate dehydrogenase, aminopeptidase, and acidic amino acid decarboxylase. The DEGs with functions associated with flight capacity and muscles exhibited higher transcript levels for younger (2 d--old) virgins. This comprehensive C. sinica transcriptomic data provide a foundation for a better understanding of the molecular mechanisms underlying the flight capacity to meet the physiological demands of flight muscles in C. sinica.