The genomic and cellular basis of biosynthetic innovation in rove beetles.

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

Patterns of cytonuclear discordance and divergence between subspecies of the scarlet macaw (Ara macao) in Central America.

The scarlet macaw, Ara macao, is a neotropical parrot that contains two described subspecies with broadly discrete geographical distributions. One subspecies, A. m. macao, is found from South America north into southwestern Costa Rica, while the second subspecies, A. m. cyanoptera, is found from eastern Costa Rica north into central Mexico. Our previous research using mitochondrial data to examine phylogeographical divergence across the collective range of these two subspecies concluded that they represent distinct evolutionary entities, with minimal contemporary hybridization between them. Here we further examine phylogenetic relationships and patterns of genetic variation between these two subspecies using a dataset of genetic markers derived from their nuclear genomes. Our analyses show clear nuclear divergence between A. m. macao and A. m. cyanoptera in Central America. Collectively however, samples from this region appear genetically more similar to one another than they do to the examined South American (Brazilian) A. m. macao sample. This observation contradicts our previous assessments based on mitochondrial DNA analyses that A. m. macao in Central and South America represent a single phylogeographical group that is evolutionarily distinct from Central American A. m. cyanoptera. Nonetheless, in agreement with our previous findings, ongoing genetic exchange between the two subspecies appears limited. Rather, our analyses indicate that incomplete lineage sorting is the best supported explanation for cytonuclear discordance within these parrots. High-altitude regions in Central America may act as a reproductive barrier, limiting contemporary hybridization between A. m. macao and A. m. cyanoptera. The phylogeographic complexities of scarlet macaw taxa in this region highlight the need for additional evolutionary examinations of these populations.

Demographic Expansions and the Emergence of Host Specialization in Genetically Distinct Ecotypes of the Tick-Transmitted Bacterium Anaplasma phagocytophilum.

In Europe, genetically distinct ecotypes of the tick-vectored bacterium Anaplasma phagocytophilum circulate among mammals in three discrete enzootic cycles. To date, potential ecological factors that contributed to the emergence of these divergent ecotypes have been poorly studied. Here, we show that the ecotype that predominantly infects roe deer (Capreolus capreolus) is evolutionarily derived. Its divergence from a host generalist ancestor occurred after the last glacial maximum as mammal populations, including roe deer, recolonized the European mainland from southern refugia. We also provide evidence that this host specialist ecotype's effective population size (Ne) has tracked changes in the population of its roe deer host. Specifically, both host and bacterium have undergone substantial increases in Ne over the past 1,500 years. In contrast, we show that while it appears to have undergone a major population expansion starting ~3,500 years ago, in the past 500 years, the contemporary host generalist ecotype has experienced a substantial reduction in genetic diversity levels, possibly as a result of reduced opportunities for transmission between competent hosts. IMPORTANCE The findings of this study reveal specific events important for the evolution of host specialization in a naturally occurring, obligately intracellular bacterial pathogen. Specifically, they show that host range shifts and the emergence of host specialization may occur during periods of population growth in a generalist ancestor. Our results also demonstrate the close correlation between demographic patterns in host and pathogen for a specialist system. These findings have important relevance for understanding the evolution of host range diversity. They may inform future work on host range dynamics, and they provide insights for understanding the emergence of pathogens that have human and veterinary health implications.

Persistence and expansion of cryptic endangered red wolf genomic ancestry along the American Gulf coast.

Admixture and introgression play a critical role in adaptation and genetic rescue that has only recently gained a deeper appreciation. Here, we explored the geographical and genomic landscape of cryptic ancestry of the endangered red wolf that persists within the genome of a ubiquitous sister taxon, the coyote, all while the red wolf has been extinct in the wild since the early 1980s. We assessed admixture across 120,621 single nucleotiode polymorphism (SNP) loci genotyped in 293 canid genomes. We found support for increased red wolf ancestry along a west-to-east gradient across the southern United States associated with historical admixture in the past 100 years. Southwestern Louisiana and southeastern Texas, the geographical zone where the last red wolves were known prior to extinction in the wild, contained the highest and oldest levels of red wolf ancestry. Further, given the paucity of inferences based on chromosome types, we compared patterns of ancestry on the X chromosome and autosomes. We additionally aimed to explore the relationship between admixture timing and recombination rate variation to investigate gene flow events. We found that X-linked regions of low recombination rates were depleted of introgression, relative to the autosomes, consistent with the large X effect and enrichment with loci involved in maintaining reproductive isolation. Recombination rate was positively correlated with red wolf ancestry across coyote genomes, consistent with theoretical predictions. The geographical and genomic extent of cryptic red wolf ancestry can provide novel genomic resources for recovery plans targeting the conservation of the endangered red wolf.

Genomic signatures of spatially divergent selection at clownfish range margins.

Understanding how evolutionary forces interact to drive patterns of selection and distribute genetic variation across a species' range is of great interest in ecology and evolution, especially in an era of global change. While theory predicts how and when populations at range margins are likely to undergo local adaptation, empirical evidence testing these models remains sparse. Here, we address this knowledge gap by investigating the relationship between selection, gene flow and genetic drift in the yellowtail clownfish, Amphiprion clarkii, from the core to the northern periphery of the species range. Analyses reveal low genetic diversity at the range edge, gene flow from the core to the edge and genomic signatures of local adaptation at 56 single nucleotide polymorphisms in 25 candidate genes, most of which are significantly correlated with minimum annual sea surface temperature. Several of these candidate genes play a role in functions that are upregulated during cold stress, including protein turnover, metabolism and translation. Our results illustrate how spatially divergent selection spanning the range core to the periphery can occur despite the potential for strong genetic drift at the range edge and moderate gene flow from the core populations.

The urban-adapted underground mosquito Culex pipiens form molestus maintains exogenously influenceable circadian rhythms.

Genes known to affect circadian rhythms (i.e. 'clock genes') also influence the photoperiodic induction of overwintering reproductive diapause in the northern house mosquito, Culex pipiens f. pipiens. This suggests that molecular changes in one or more clock genes could contribute to the inability to diapause in a second form of this mosquito, Culex pipiens f. molestus. Temperate populations of Cx. pipiens f. molestus inhabit underground locations generally devoid of predictable photoperiods. For this reason, there could be limited fitness consequences if the hypothesized molecular changes to its clock genes also eliminated this mosquito's ability to regulate circadian rhythms in response to photoperiod variation. Here, we demonstrate that in contrast to this prediction, underground derived Cx. pipiens f. molestus retain exogenously influenceable circadian rhythms. Nonetheless, our genetic analyses indicate that the gene Helicase domino (dom) has a nine-nucleotide, in-frame deletion specific to Cx. pipiens f. molestus. Previous work has shown that splice variants in this gene differentially influence circadian behavior in Drosophila melanogaster. We also find derived, non-synonymous single nucleotide polymorphisms (SNPs) in eight genes that may also affect circadian rhythms and/or diapause induction in Cx. pipiens f. molestus. Finally, four putative circadian genes were found to have no quantifiable expression during any examined life stage, suggesting potential regulatory effects. Collectively, our findings indicate that the distinct, but molecularly interconnected life-history traits of diapause induction and circadian rhythms are decoupled in Cx. pipiens f. molestus and suggest this taxon may be a valuable tool for exploring exogenously influenced phenotypes in mosquitoes more broadly.

Updating the Bibliography of Interbreeding among Canis in North America.

This bibliography provides a collection of references that documents the evolution of studies evidencing interbreeding among Canis species in North America. Over the past several decades, advances in biology and genomic technology greatly improved our ability to detect and characterize species interbreeding, which has significance for understanding species in a changing landscape as well as for endangered species management. This bibliography includes a discussion within each category of interbreeding, the timeline of developing evidence, and includes a review of past research conducted on experimental crosses. Research conducted in the early 20th century is rich with detailed records and photographs of hybrid offspring development and behavior. With the progression of molecular methods, studies can estimate historical demographic parameters and detect chromosomal patterns of ancestry. As these methods continue to increase in accessibility, the field will gain a deeper and richer understanding of the evolutionary history of North American Canis.

Prosthetic Valve Endocarditis with Bartonella washoensis in a Human European Patient and Its Detection in Red Squirrels (Sciurus vulgaris).

Members of the genus Bartonella are fastidious Gram-negative facultative intracellular bacteria that are typically transmitted by arthropod vectors. Several Bartonella spp. have been found to cause culture-negative endocarditis in humans. Here, we report the case of a 75-year-old German woman with prosthetic valve endocarditis due to Bartonella washoensis The infecting agent was characterized by sequencing of six housekeeping genes (16S rRNA, ftsZ, gltA, groEL, ribC, and rpoB), applying a multilocus sequence typing (MLST) approach. The 5,097 bp of the concatenated housekeeping gene sequence from the patient were 99.0% identical to a sequence from a B. washoensis strain isolated from a red squirrel (Sciurus vulgarisorientis) from China. A total of 39% (24/62) of red squirrel (S. vulgaris) samples from the Netherlands were positive for the B. washoensisgltA gene variant detected in the patient. This suggests that the red squirrel is the reservoir host for human infection in Europe.

Exploring Chihuahuan Desert diversification in the gray-banded kingsnake, Lampropeltis alterna (Serpentes: Colubridae).

Within many biomes, the cause of phylogeographic structure remains unknown even across regions throughout North America, including within the biodiverse Chihuahuan Desert. For example, little is known about population structure or the timing of diversification of Chihuahuan endemics. This is due largely to the lack of population genomic studies within this region. We generated ultra-conserved element data for the gray-banded kingsnake (Lampropeltis alterna) to investigate lineage divergence and historical demography across the Chihuahuan Desert. We found three unique lineages corresponding to the Trans-Pecos and Mapimian biogeographic regions of the Chihuahuan Desert, and a distinct population in the Sierra Madre Occidental. Using several mutation rates to calibrate the timing of divergence among these lineages, we show that lineage divergence likely occurred during the Pleistocene, which indicates that careful consideration needs to be used when applying mutation rates to ultra-conserved elements. We suggest that biogeographic provinces within the Chihuahuan Desert may have served as allopatric refugia during climatic fluctuations of the Quaternary. This work serves as an important template for further testing biogeographic hypotheses within the region.

Sexual Dimorphism and Retinal Mosaic Diversification following the Evolution of a Violet Receptor in Butterflies.

Numerous animal lineages have expanded and diversified the opsin-based photoreceptors in their eyes underlying color vision behavior. However, the selective pressures giving rise to new photoreceptors and their spectral tuning remain mostly obscure. Previously, we identified a violet receptor (UV2) that is the result of a UV opsin gene duplication specific to Heliconius butterflies. At the same time the violet receptor evolved, Heliconius evolved UV-yellow coloration on their wings, due to the pigment 3-hydroxykynurenine (3-OHK) and the nanostructure architecture of the scale cells. In order to better understand the selective pressures giving rise to the violet receptor, we characterized opsin expression patterns using immunostaining (14 species) and RNA-Seq (18 species), and reconstructed evolutionary histories of visual traits in five major lineages within Heliconius and one species from the genus Eueides. Opsin expression patterns are hyperdiverse within Heliconius. We identified six unique retinal mosaics and three distinct forms of sexual dimorphism based on ommatidial types within the genus Heliconius. Additionally, phylogenetic analysis revealed independent losses of opsin expression, pseudogenization events, and relaxation of selection on UVRh2 in one lineage. Despite this diversity, the newly evolved violet receptor is retained across most species and sexes surveyed. Discriminability modeling of behaviorally preferred 3-OHK yellow wing coloration suggests that the violet receptor may facilitate Heliconius color vision in the context of conspecific recognition. Our observations give insights into the selective pressures underlying the origins of new visual receptors.

Varying influences of selection and demography in host-adapted populations of the tick-transmitted bacterium, Anaplasma phagocytophilum.

BACKGROUND: The host range of a pathogenic bacterial strain likely influences its effective population size, which in turn affects the efficacy of selection. Transmission between competent hosts may occur more frequently for host generalists than for specialists. This could allow higher bacterial population densities to persist within an ecological community and increase the efficacy of selection in these populations. Conversely, specialist strains may be better adapted to their hosts and consequently achieve greater within-host population densities, with corresponding increases in selection efficacy. To assess these different hypotheses, we examined the effective population sizes of three strains of the bacterium Anaplasma phagocytophilum and categorized the varying roles of selection and demography on patterns of genetic diversity and divergence in these populations. A. phagocytophilum is a tick-transmitted, obligately intracellular pathogen. Strains of A. phagocytophilum display varying degrees of host specialization, making this a good species for exploring questions regarding host range, effective population size and selection efficacy. RESULTS: We found that a roe deer specialist harbored the most genetic diversity of the three A. phagocytophilum strains and correspondingly had the largest effective population size. Another strain that is ecologically specialized on rodents and insectivores had the smallest effective population size. However, these mammalian hosts are distantly related evolutionarily. The third strain, a host generalist, was intermediate in its effective population size between the other two strains. Evolutionary constraint on non-synonymous sites was pervasive in all three strains, although some slightly deleterious mutations may also be segregating in these populations. We additionally found evidence of genome-wide selective sweeps in the generalist strain, whereas signals of repeated bottlenecks were detected in the strain with the smallest effective population size. CONCLUSIONS: A. phagocytophilum is a diverse bacterial species that differs among distinct strains in its effective population size, as well as how genetic diversity and divergence have been influenced by selection and demographic changes. In this species, host specialization may facilitate increased population growth and allow more opportunities for selection to act. These results provide insights into how host range has influenced evolutionary patterns of strain divergence in an emerging zoonotic bacterium.

Sex chromosome mosaicism and hybrid speciation among tiger swallowtail butterflies.

Hybrid speciation, or the formation of a daughter species due to interbreeding between two parental species, is a potentially important means of diversification, because it generates new forms from existing variation. However, factors responsible for the origin and maintenance of hybrid species are largely unknown. Here we show that the North American butterfly Papilio appalachiensis is a hybrid species, with genomic admixture from Papilio glaucus and Papilio canadensis. Papilio appalachiensis has a mosaic phenotype, which is hypothesized to be the result of combining sex-linked traits from P. glaucus and P. canadensis. We show that P. appalachiensis' Z-linked genes associated with a cooler thermal habitat were inherited from P. canadensis, whereas its W-linked mimicry and mitochondrial DNA were inherited from P. glaucus. Furthermore, genome-wide AFLP markers showed nearly equal contributions from each parental species in the origin of P. appalachiensis, indicating that it formed from a burst of hybridization between the parental species, with little subsequent backcrossing. However, analyses of genetic differentiation, clustering, and polymorphism based on molecular data also showed that P. appalachiensis is genetically distinct from both parental species. Population genetic simulations revealed P. appalachiensis to be much younger than the parental species, with unidirectional gene flow from P. glaucus and P. canadensis into P. appalachiensis. Finally, phylogenetic analyses, combined with ancestral state reconstruction, showed that the two traits that define P. appalachiensis' mosaic phenotype, obligatory pupal diapause and mimicry, evolved uniquely in P. canadensis and P. glaucus, respectively, and were then recombined through hybridization to form P. appalachiensis. These results suggest that natural selection and sex-linked traits may have played an important role in the origin and maintenance of P. appalachiensis as a hybrid species. In particular, ecological barriers associated with a steep thermal cline appear to maintain the distinct, mosaic genome of P. appalachiensis despite contact and occasional hybridization with both parental species.

Genomic analyses indicate the North American Ap-ha variant of the tick-vectored bacterium Anaplasma phagocytophilum was introduced from Europe.

BACKGROUND: Anaplasma phagocytophilum is a tick-vectored, obligately intracellular bacterium that infects a diversity of vertebrate hosts. In North America, the Ap-ha variant of A. phagocytophilum can cause dangerous infections in humans, whereas symptomatic human infections in Europe are rare. Conversely, the European host-generalist ecotype of A. phagocytophilum frequently causes illness in domestic ruminants while no comparable infections have been recorded from North America. Despite these differences in pathogenicity, the Ap-ha variant is closely aligned phylogenetically with the European host-generalist ecotype. Furthermore, North American populations of A. phagocytophilum are less genetically diverse than those in Europe. Taken together, these observations suggest that the North American Ap-ha variant may represent an introduced population of this bacterium. METHODS: Data from publicly available whole genomes of A. phagocytophilum were used to compare phylogeographic patterns and the extent of genetic divergence between the North American Ap-ha variant and the European host-generalist ecotype. RESULTS: The results confirm that North American Ap-ha samples are phylogenetically nested within the diversity of the European host-generalist ecotype, and that Ap-ha likely radiated within the last 100 years. As expected, the Ap-ha variant also exhibited relatively low genetic diversity levels compared to the European host-generalist ecotype. Finally, North American Ap-ha harbored significantly more derived alleles than the European host-generalist A. phagocytophilum population. CONCLUSIONS: Collectively, these results support the hypothesis that the Ap-ha variant was recently introduced to North America from Europe and underwent a strong genetic bottleneck during this process (i.e. a 'founder event'). Adaptation to novel vectors may have also played a role in shaping genetic diversity and divergence patterns in these pathogenic bacteria. These findings have implications for future studies aimed at understanding evolutionary patterns and pathogenicity variation within A. phagocytophilum.

Granulocytic anaplasmosis in cats from central Europe and molecular characterization of feline Anaplasma phagocytophilum strains by ankA gene, groEL gene and multilocus sequence typing.

BACKGROUND: Anaplasma phagocytophilum is a Gram-negative obligate intracellular bacterium that replicates in neutrophil granulocytes. It is transmitted by ticks of the Ixodes ricinus complex and causes febrile illness called granulocytic anaplasmosis primarily in humans, horses, dogs, sheep, cattle and goats. In comparison, clinically apparent disease has been described rarely in cats especially compared to dogs and horses. It is currently unknown whether cats are less susceptible to A. phagocytophilum or whether granulocytic anaplasmosis might be underdiagnosed in cats. METHODS: To address this question, we examined clinical signs and laboratory findings in seven A. phagocytophilum infected cats from Germany and Switzerland. We then genetically characterized feline A. phagocytophilum strains and compared them to those from other hosts showing clinically apparent disease. For this purpose, ankA-based, groEL-based and multilocus sequence typing (MLST) were applied. Furthermore, the concordance between these typing methods was assessed. RESULTS: Fever, lethargy and anorexia were the most common clinical signs in cats suffering from granulocytic anaplasmosis. The most frequent laboratory finding was thrombocytopenia. All three typing methods consistently indicated that the A. phagocytophilum strains found infecting cats are the same as those that cause disease in humans, dogs and horses. In general, the three typing methods applied exhibited high concordance. CONCLUSIONS: The genetic characterization of the feline A. phagocytophilum strains indicates that strain divergence is not the explanation for the fact that granulocytic anaplasmosis is much less frequently diagnosed in cats than in dogs and horses. Otherwise, it may be possible that cats are less susceptible to the same strains than dogs and horse are. However, due to the unspecific clinical signs, it should be considered that granulocytic anaplasmosis may be under-diagnosed in cats.

Predatory fireflies and their toxic firefly prey have evolved distinct toxin resistance strategies.

Toxic cardiotonic steroids (CTS) act as a defense mechanism in many firefly species (Lampyridae) by inhibiting a crucial enzyme called Na+,K+-ATPase (NKA). While most fireflies produce these toxins internally, species of the genus Photuris acquire them from a surprising source: predation on other fireflies. The contrasting physiology of toxin exposure and sequestration between Photuris and other firefly genera suggests that distinct strategies may be required to prevent self-intoxication. Our study demonstrates that both Photuris and their firefly prey have evolved highly-resistant NKAs. Using an evolutionary analysis of the specific target of CTS (ATPα) in fireflies, and gene-editing in Drosophila, we find that the initial steps towards resistance were shared among Photuris and other firefly lineages. However, the Photuris lineage subsequently underwent multiple rounds of gene duplication and neofunctionalization, resulting in the development of ATPα paralogs that are differentially expressed and exhibit increasing resistance to CTS. In contrast, other firefly species have maintained a single copy. Our results implicate gene duplication as a facilitator in the transition of Photuris to its distinct ecological role as predator of toxic firefly prey.

Predatory fireflies and their toxic firefly prey have evolved distinct toxin resistance strategies.

Toxic cardiotonic steroids (CTSs) act as a defense mechanism in many firefly species (Lampyridae) by inhibiting a crucial enzyme called Na+,K+-ATPase (NKA). Although most fireflies produce these toxins internally, species of the genus Photuris acquire them from a surprising source: predation on other fireflies. The contrasting physiology of toxin exposure and sequestration between Photuris and other firefly genera suggests that distinct strategies may be required to prevent self-intoxication. Our study demonstrates that both Photuris and their firefly prey have evolved highly resistant NKAs. Using an evolutionary analysis of the specific target of CTS (ATPα) in fireflies and gene editing in Drosophila, we find that the initial steps toward resistance were shared among Photuris and other firefly lineages. However, the Photuris lineage subsequently underwent multiple rounds of gene duplication and neofunctionalization, resulting in the development of ATPα paralogs that are differentially expressed and exhibit increasing resistance to CTS. By contrast, other firefly species have maintained a single copy. Our results implicate gene duplication as a facilitator in the transition of Photuris to its distinct ecological role as a predator of toxic firefly prey.

The genomic and cellular basis of biosynthetic innovation in rove beetles.

How evolution at the cellular level potentiates change at the macroevolutionary level is a major question in evolutionary biology. With >66,000 described species, rove beetles (Staphylinidae) comprise the largest metazoan family. Their exceptional radiation has been coupled to pervasive biosynthetic innovation whereby numerous lineages bear defensive glands with diverse chemistries. Here, we combine comparative genomic and single-cell transcriptomic data from across the largest rove beetle clade, Aleocharinae. We retrace the functional evolution of two novel secretory cell types that together comprise the tergal gland-a putative catalyst behind Aleocharinae's megadiversity. We identify key genomic contingencies that were critical to the assembly of each cell type and their organ-level partnership in manufacturing the beetle's defensive secretion. This process hinged on evolving a mechanism for regulated production of noxious benzoquinones that appears convergent with plant toxin release systems, and synthesis of an effective benzoquinone solvent that weaponized the total secretion. We show that this cooperative biosynthetic system arose at the Jurassic-Cretaceous boundary, and that following its establishment, both cell types underwent ∼150 million years of stasis, their chemistry and core molecular architecture maintained almost clade-wide as Aleocharinae radiated globally into tens of thousands of lineages. Despite this deep conservation, we show that the two cell types have acted as substrates for the emergence of adaptive, biochemical novelties-most dramatically in symbiotic lineages that have infiltrated social insect colonies and produce host behavior-manipulating secretions. Our findings uncover genomic and cell type evolutionary processes underlying the origin, functional conservation and evolvability of a chemical innovation in beetles.

The evolution of cardenolide-resistant forms of Na⁺,K⁺ -ATPase in Danainae butterflies.

Cardenolides are a class of plant secondary compounds that inhibit the proper functioning of the Na(+) , K(+) -ATPase enzyme in susceptible animals. Nonetheless, many insect species are able to sequester cardenolides for their own defence. These include butterflies in the subfamily Danainae (Family: Nymphalidae) such as the monarch (Danaus plexippus). Previous studies demonstrated that monarchs harbour an asparagine (N) to histidine (H) substitution (N122H) in the α subunit of Na(+) , K(+) -ATPase (ATPα) that reduces this enzyme's sensitivity to cardenolides. More recently, it has been suggested that at ATPα position 111, monarchs may also harbour a leucine (L)/glutamine (Q) polymorphism. This later amino acid could also contribute to cardenolide insensitivity. However, here we find that incorrect annotation of the initially reported DNA sequence for ATPα has led to several erroneous conclusions. Using a population genetic and phylogenetic analysis of monarchs and their close relatives, we show that an ancient Q111L substitution occurred prior to the radiation of all Danainae, followed by a second substitution at the same site to valine (V), which arose before the diversification of the Danaus genus. In contrast, N122H appears to be a recent substitution specific to monarchs. Surprisingly, examination of a broader insect phylogeny reveals that the same progression of amino acid substitutions (Q111L → L111V + N122H) has also occurred in Chyrsochus beetles (Family: Chrysomelidae, Subfamily: Eumolpinae) that feed on cardenolide-containing host plants. The parallel pattern of amino acid substitution in these two distantly related lineages is consistent with an adaptive role for these substitutions in reducing cardenolide sensitivity and suggests that their temporal order may be limited by epistatic interactions.

A gene-based capture assay for surveying patterns of genetic diversity and insecticide resistance in a worldwide group of invasive mosquitoes.

Understanding patterns of diversification, genetic exchange, and pesticide resistance in arthropod disease vectors is necessary for effective population management. With the availability of next-generation sequencing technologies, one of the best approaches for surveying such patterns involves the simultaneous genotyping of many samples for a large number of genetic markers. To this end, the targeting of gene sequences of known function can be a cost-effective strategy. One insect group of substantial health concern are the mosquito taxa that make up the Culex pipiens complex. Members of this complex transmit damaging arboviruses and filariae worms to humans, as well as other pathogens such as avian malaria parasites that are detrimental to birds. Here we describe the development of a targeted, gene-based assay for surveying genetic diversity and population structure in this mosquito complex. To test the utility of this assay, we sequenced samples from several members of the complex, as well as from distinct populations of the relatively under-studied Culex quinquefasciatus. The data generated was then used to examine taxonomic divergence and population clustering between and within these mosquitoes. We also used this data to investigate genetic variants present in our samples that had previously been shown to correlate with insecticide-resistance. Broadly, our gene capture approach successfully enriched the genomic regions of interest, and proved effective for facilitating examinations of taxonomic divergence and geographic clustering within the Cx. pipiens complex. It also allowed us to successfully survey genetic variation associated with insecticide resistance in Culex mosquitoes. This enrichment protocol will be useful for future studies that aim to understand the genetic mechanisms underlying the evolution of these ubiquitous and increasingly damaging disease vectors.

Correction: Can public online databases serve as a source of phenotypic information for Cannabis genetic association studies?

[This corrects the article DOI: 10.1371/journal.pone.0247607.].