The doublesex (dsx) gene functions as a molecular switch at the base of the insect sex determination cascade and triggers male or female somatic sexual differentiation in Drosophila. Having been reported from only seven current insect orders, the exact phylogenetic distribution of dsx within the largest Arthropod sub-phylum, the Hexapoda, is unknown. To understand the evolution of this integral gene relative to other arthropods, we tested for the presence of dsx within public EST and genome sequencing projects representative of all 32 hexapod orders. We find the dsx gene to be ubiquitous, with putative orthologs recovered from 30 orders. Additionally, we recovered both alternatively spliced and putative paralogous dsx transcripts from several orders of hexapods, including basal lineages, indicating the likely presence of these characteristics in the hexapod common ancestor. Of note, other arthropods such as chelicerates and crustaceans express two dsx genes, both of which are shown to lack alternative splicing. Furthermore, we discovered a large degree of length heterogeneity in the common region of dsx coding sequences within and among orders, possibly resulting from lineage-specific selective pressures inherent to each taxon. Our work serves as a valuable resource for understanding the evolution of sex determination in insects.

The doublesex gene controls somatic sexual differentiation of many metazoan species, including the malaria mosquito Anopheles gambiae and the dengue and yellow fever vector Aedes aegypti (Diptera: Culicidae). As in other studied dipteran dsx homologs, the gene maintains functionality via evolutionarily conserved protein domains and sex-specific alternative splicing. The upstream factors that regulate splicing of dsx and the manner in which they do so however remain variable even among closely related organisms. As the induction of sex ratio biases is a central mode of action in many emerging molecular insecticides, it is imperative to elucidate as much of the sex determination pathway as possible in the mosquito disease vectors. Here we report the full-length gene sequence of the doublesex gene in Culex quinquefasciatus (Cxqdsx) and its male and female-specific isoforms. Cxqdsx maintains characteristics possibly derived in the Culicinae and present in the Aedes aegypti dsx gene (Aeadsx) such as gain of exon 3b and the presence of Rbp1 cis-regulatory binding sites, and also retains presumably ancestral attributes present in Anopheles gambiae such as maintenance of a singular female-specific exon 5. Unlike in Aedes aegypti, we find no evidence for intron gain in the female transcript(s), yet recover a second female isoform generated via selection of an alternate splice donor. Utilizing next-gen sequence (NGS) data, we complete the Aeadsx gene model and identify a putative core promoter region in both Aeadsx and Cxqdsx. Also utilizing NGS data, we construct a full-length gene sequence for the dsx homolog of the northern house mosquito Culex pipiens form pipiens (Cxpipdsx). Analysis of peptide evolutionary rates between Cxqdsx and Cxpipdsx (both members of the Culex pipiens complex) shows the male-specific portion of the transcript to have evolved rapidly with respect to female-specific and common regions. As in other studied insects, doublesex maintains sex-specific splicing and conserved doublesex/mab-3 domains in the mosquitoes Culex quinquefasciatus and Cx. pipiens. The cis-regulated splicing of Cxqdsx does not appear to follow either currently described mosquito model (for An. gambiae and Ae. aegypti); each of the three mosquito genera exhibit evidence of unique cis-regulatory mechanisms. The male-specific dsx terminus exhibits rapid peptide evolutionary rates, even among closely related sibling species.

Culex pipiens, an invasive mosquito and vector of West Nile virus in the US, has two morphologically indistinguishable forms that differ dramatically in behavior and physiology. Cx. pipiens form pipiens is primarily a bird-feeding temperate mosquito, while the sub-tropical Cx. pipiens form molestus thrives in sewers and feeds on mammals. Because the feral form can diapause during the cold winters but the domestic form cannot, the two Cx. pipiens forms are allopatric in northern Europe and, although viable, hybrids are rare. Cx. pipiens form molestus has spread across all inhabited continents and hybrids of the two forms are common in the US. Here we elucidate the genes and gene families with the greatest divergence rates between these phenotypically diverged mosquito populations, and discuss them in light of their potential biological and ecological effects. After generating and assembling novel transcriptome data for each population, we performed pairwise tests for nonsynonymous divergence (Ka) of homologous coding sequences and examined gene ontology terms that were statistically over-represented in those sequences with the greatest divergence rates. We identified genes involved in digestion (serine endopeptidases), innate immunity (fibrinogens and α-macroglobulins), hemostasis (D7 salivary proteins), olfaction (odorant binding proteins) and chitin binding (peritrophic matrix proteins). By examining molecular divergence between closely related yet phenotypically divergent forms of the same species, our results provide insights into the identity of rapidly-evolving genes between incipient species. Additionally, we found that families of signal transducers, ATP synthases and transcription regulators remained identical at the amino acid level, thus constituting conserved components of the Cx. pipiens proteome. We provide a reference with which to gauge the divergence reported in this analysis by performing a comparison of transcriptome sequences from conspecific (yet allopatric) populations of another member of the Cx. pipiens complex, Cx. quinquefasciatus.

The documentation of endogenous viral elements (EVEs; virus-derived genetic material integrated into the genome of a nonviral host) has offered insights into how arthropods respond to viral infection via RNA interference pathways. Small non-coding RNAs derived from EVE loci serve to direct RNAi pathways in limiting replication and infection from cognate viruses, thus benefiting the host’s fitness and, potentially, vectorial capacity. Here we use informatic approaches to analyze nine available genome sequences of hard ticks (Acari: Ixodidae; Rhipicephalus sanguineus, R. microplus, R. annulatus, Ixodes ricinus, I. persulcatus, I. scapularis, Hyalomma asiaticum, Haemaphysalis longicornis, and Dermacentor silvarum) to identify endogenous viral elements and to illustrate the shared ancestry of all elements identified. Our results highlight a broad diversity of viral taxa as having given rise to 1234 identified EVEs in ticks, with Mononegavirales (specifically Rhabdoviridae) well-represented in this subset of hard ticks. Further investigation revealed extensive adintovirus integrations in several Ixodes species, the prevalence of Bunyavirales EVEs (notably not observed in mosquitoes), and the presence of several elements similar to known emerging human and veterinary pathogens. These results will inform subsequent work on current and past associations with tick species with regard to the viruses from which their “viral fossils” are derived and may serve as a reference for quality control of various tick-omics data that may suffer from misidentification of EVEs as viral genetic material.

Lone star ticks (Amblyomma americanum L.) are expanding within the northeast United States, a region historically focused on Ixodes scapularis-transmitted diseases. In Monmouth County, NJ, the shift has been dramatic, and lone star ticks now vastly outnumber blacklegged ticks. As a result, there is an enhanced need to focus on the potential health risks of A. americanum-transmitted pathogens, such as the emerging Heartland (HRTV) and Bourbon (BRBV) viruses. We screened 1,205 nymphal lone star ticks for HRTV and BRBV using RT-qPCR assays and detected BRBV in 3 ticks collected in Monmouth County, NJ, in 2021. Additionally, we sequenced a complete BRBV genome from a single infected specimen, finding 99.4% identity with human pathogenic isolates from the eastern-central United States. Our results have important public health implications for a region only recently becoming aware of public health risks posed by lone star ticks. Of note, we report successful detection of viral RNA in samples that were stored and intended for DNA preservation, for example, kept in ethanol at room temperature, which may reduce barriers for public health agencies seeking to expand their tick testing to include viruses.