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.

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.

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.

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.

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.

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.

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.

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.

The mosquito taxa that make up the Culex pipiens complex are important vectors of the agents of several human diseases such as West Nile and St. Louis encephalitides, and lymphatic filariasis. They are also important vectors of avian malaria, which impacts livestock and wildlife. The development of effective strategies for the control of these mosquitoes requires knowledge of their origins, distribution, dispersal patterns, and the extent to which discreet taxonomic entities within the complex interbreed. To achieve these objectives, it is necessary to compare patterns of genetic diversity across many mosquito samples, which can be cost-prohibitive. To address this limitation, we developed a targeted, gene-based assay that allowed us to cost-effectively genotype a large number of genetic variants from a representative global sampling of individual Cx. pipiens complex mosquitoes. We show that this assay is a powerful tool for examining genetic structure and hybridization among populations. We also explore its utility for surveying alleles previously shown to be associated with insecticide resistance. Future use of this enrichment assay and the bioinformatics methods described here will allow researchers to study evolutionary patterns across the Cx. pipiens complex as well as monitor the presence of genetic variation that could affect control efforts.

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.

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.

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.

Paleogenomics illuminates the evolutionary history of the extinct Holocene "horned" crocodile of Madagascar, Voay robustus

Ancient DNA is transforming our ability to reconstruct historical patterns and mechanisms shaping modern diversity and distributions. In particular, molecular data from extinct Holocene island faunas have revealed surprising biogeographic scenarios. Here, we recovered partial mitochondrial (mt) genomes for 1300-1400 year old specimens (n = 2) of the extinct "horned" crocodile, Voay robustus, collected from Holocene deposits in southwestern Madagascar. Phylogenetic analyses of partial mt genomes and tip-dated timetrees based on molecular, fossil, and stratigraphic data favor a sister group relationship between Voay and Crocodylus (true crocodiles). These well supported trees conflict with recent morphological systematic work that has consistently placed Voay within Osteolaeminae (dwarf crocodiles and kin) and provide evidence for likely homoplasy in crocodylian cranial anatomy and snout shape. The close relationship between Voay and Crocodylus lends additional context for understanding the biogeographic origins of these genera and refines competing hypotheses for the recent extinction of Voay from Madagascar.

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

The use of Cannabis is gaining greater social acceptance for its beneficial medicinal and recreational uses. With this acceptance has come new opportunities for crop management, selective breeding, and the potential for targeted genetic manipulation. However, as an agricultural product Cannabis lags far behind other domesticated plants in knowledge of the genes and genetic variation that influence plant traits of interest such as growth form and chemical composition. Despite this lack of information, there are substantial publicly available resources that document phenotypic traits believed to be associated with particular Cannabis varieties. Such databases could be a valuable resource for developing a greater understanding of genes underlying phenotypic variation if combined with appropriate genetic information. To test this potential, we collated phenotypic data from information available through multiple online databases. We then produced a Cannabis SNP database from 845 strains to examine genome wide associations in conjunction with our assembled phenotypic traits. Our goal was not to locate Cannabis-specific genetic variation that correlates with phenotypic variation as such, but rather to examine the potential utility of these databases more broadly for future, explicit genome wide association studies (GWAS), either in stand-alone analyses or to complement other types of data. For this reason, we examined a very broad array of phenotypic traits. In total, we performed 201 distinct association tests using web-derived phenotype data appended to 290 uniquely named Cannabis strains. Our results indicated that chemical phenotypes, such as tetrahydrocannabinol (THC) and cannabidiol (CBD) content, may have sufficiently high-quality information available through web-based sources to allow for genetic association inferences. In many cases, variation in chemical traits correlated with genetic variation in or near biologically reasonable candidate genes, including several not previously implicated in Cannabis chemical variation. As with chemical phenotypes, we found that publicly available data on growth traits such as height, area of growth, and floral yield may be precise enough for use in future association studies. In contrast, phenotypic information for subjective traits such as taste, physiological affect, neurological affect, and medicinal use appeared less reliable. These results are consistent with the high degree of subjectivity for such trait data found on internet databases, and suggest that future work on these important but less easily quantifiable characteristics of Cannabis may require dedicated, controlled phenotyping.

Genomic Analysis of the Only Blind Cichlid Reveals Extensive Inactivation in Eye and Pigment Formation Genes

Trait loss represents an intriguing evolutionary problem, particularly when it occurs across independent lineages. Fishes in light-poor environments often evolve “troglomorphic” traits, including reduction or loss of both pigment and eyes. Here, we investigate the genomic basis of trait loss in a blind and depigmented African cichlid, Lamprologus lethops, and explore evolutionary forces (selection and drift) that may have contributed to these losses. This species, the only known blind cichlid, is endemic to the lower Congo River. Available evidence suggests that it inhabits deep, low-light habitats. Using genome sequencing, we show that genes related to eye formation and pigmentation, as well as other traits associated with troglomorphism, accumulated inactivating mutations rapidly after speciation. A number of the genes affected in L. lethops are also implicated in troglomorphic phenotypes in Mexican cavefish (Astyanax mexicanus) and other species. Analysis of heterozygosity patterns across the genome indicates that L. lethops underwent a significant population bottleneck roughly 1 Ma, after which effective population sizes remained low. Branch-length tests on a subset of genes with inactivating mutations show little evidence of directional selection; however, low overall heterozygosity may reduce statistical power to detect such signals. Overall, genome-wide patterns suggest that accelerated genetic drift from a severe bottleneck, perhaps aided by directional selection for the loss of physiologically expensive traits, caused inactivating mutations to fix rapidly in this species.

Global evaluation of taxonomic relationships and admixture within the Culex pipiens complex of mosquitoes

BACKGROUND

Within the Culex pipiens mosquito complex, there are six contemporarily recognized taxa: Cx. quinquefasciatus, Cx. pipiens f. pipiens, Cx. pipiens f. molestus, Cx. pipiens pallens, Cx. australicus and Cx. globocoxitus. Many phylogenetic aspects within this complex have eluded resolution, such as the relationship of the two Australian endemic taxa to the other four members, as well as the evolutionary origins and taxonomic status of Cx. pipiens pallens and Cx. pipiens f. molestus. Ultimately, insights into lineage relationships within the complex will facilitate a better understanding of differential disease transmission by these mosquitoes. To this end, we have combined publicly available data with our own sequencing efforts to examine these questions.

RESULTS

We found that the two Australian endemic complex members, Cx. australicus and Cx. globocoxitus, comprise a monophyletic group, are genetically distinct, and are most closely related to the cosmopolitan Cx. quinquefasciatus. Our results also show that Cx. pipiens pallens is genetically distinct, but may have arisen from past hybridization. Lastly, we observed complicated patterns of genetic differentiation within and between Cx. pipiens f. pipiens and Cx. pipiens f. molestus.

CONCLUSIONS

Two Australian endemic Culex taxa, Cx. australicus and Cx. globocoxitus, belong within the Cx. pipiens complex, but have a relatively older evolutionary origin. They likely diverged from Cx. quinquefasciatus after its colonization of Australia. The taxon Cx. pipiens pallens is a distinct evolutionary entity that likely arose from past hybridization between Cx. quinquefasciatus and Cx. pipiens f. pipiens/Cx. pipiens f. molestus. Our results do not suggest it derives from ongoing hybridization. Finally, genetic differentiation within the Cx. pipiens f. pipiens and Cx. pipiens f. molestus samples suggests that they collectively form two separate geographic clades, one in North America and one in Europe and the Mediterranean. This may indicate that the Cx. pipiens f. molestus form has two distinct origins, arising from Cx. pipiens f. pipiens in each region. However, ongoing genetic exchange within and between these taxa have obscured their evolutionary histories, and could also explain the absence of monophyly among our samples. Overall, this work suggests many avenues that warrant further investigation.

Genetic characterization of Anaplasma phagocytophilum strains from goats (Capra aegagrus hircus) and water buffalo (Bubalus bubalis) by 16S rRNA gene, ankA gene and multilocus sequence typing

Anaplasma phagocytophilum is a Gram-negative obligate intracellular bacterium that replicates in neutrophil granulocytes. It is transmitted by ticks and causes tick-borne fever in domestic ruminants such as sheep, cattle and goats. However, in contrast to sheep and cattle little is known about the clinical course of infection in goats. We report here on three cases of symptomatic infection with A. phagocytophilum in two goats (Capra aegagrus hircus) and one water buffalo (Bubalus bubalis). The animals showed symptoms and laboratory findings similar to sheep and cattle. To our knowledge, this is the first report on the symptomatic infection of water buffalos with A. phagocytophilum. The infecting strains were genetically characterized by 16S rRNA gene, ankA gene and multilocus sequence typing (MLST). Four other strains from asymptomatically infected goats were also included. The ankA sequences from five goats were part of the formerly described ankA gene clusters I and IV that are known to contain A. phagocytophilum strains from sheep and cattle. However, the sequences from one goat and from the water buffalo belonged to ankA gene cluster II that was formerly described to be restricted to roe deer. A similar observation was made for MLST as three goats clustered with sequences from sheep and cattle, whereas three other goats and the water buffalo were found to be part of the roe deer cluster. However, since most of the strains from sheep and cattle were distinct from the roe deer strains, roe deer might not represent major reservoir hosts for tick-borne fever in domestic ruminants. When differing parts of the 16S rRNA gene were used for typing the results were conflicting. This shows that the use of a standardized typing method such as MLST is highly desirable to generate easily comparable results.

Activity of Genes with Functions in Human Williams-Beuren Syndrome Is Impacted by Mobile Element Insertions in the Gray Wolf Genome

In canines, transposon dynamics have been associated with a hyper-social behavioral syndrome, although the functional mechanism has yet to be described. We investigate the epigenetic and transcriptional consequences of these behavior-associated mobile element insertions (MEIs) in dogs and Yellowstone gray wolves. We posit that the transposons themselves may not be the causative feature; rather, their transcriptional regulation may exert the functional impact. We survey four outlier transposons associated with hyper-sociability, with the expectation that they are targeted for epigenetic silencing. We predict hyper-methylation of MEIs, suggestive that the epigenetic silencing of and not the MEIs themselves may be driving dysregulation of nearby genes. We found that transposon-derived sequences are significantly hyper-methylated, regardless of their copy number or species. Further, we have assessed transcriptome sequence data and found evidence that MEIs impact the expression levels of six genes (WBSCR17, LIMK1, GTF2I, WBSCR27, BAZ1B, and BCL7B), all of which have known roles in human Williams-Beuren syndrome due to changes in copy number, typically hemizygosity. Although further evidence is needed, our results suggest that a few insertions alter local expression at multiple genes, likely through a cis-regulatory mechanism that excludes proximal methylation.

The role of latitudinal, genetic and temperature variation in the induction of diapause of Papilio glaucus (Lepidoptera: Papilionidae)

A key adaptation in insects for dealing with variable environmental conditions is the ability to diapause. The tiger swallowtail butterflies, Papilio glaucus and P. canadensis are ideal species to explore the genetic causes and population genetic consequences of diapause because divergence in this trait is believed to be a salient factor in maintaining a hybrid zone between these species. Yet little is known about the factors that influence diapause induction in this system. Here we explored how spatial (latitudinal), environmental (temperature) and genetic (hybridization) factors affect diapause induction in this system. Specifically, a series of growth chamber experiments using wild caught individuals from across the eastern United States were performed to: (1) evaluate how critical photoperiod varies with latitude, (2) isolate the stage in which induction occurs, (3) test whether changes in temperature affected rates of diapause induction, and (4) explore how the incidence of diapause is affected in hybrid offspring. We find that induction occurs in the larval stage, is not sensitive to a relatively broad range of temperatures, appears to have a complex genetic basis (i.e., is not simply a dominant trait following a Mendelian inheritance pattern) and that the critical photoperiod increases by 0.4 h with each increasing degree in latitude. This work deepens our understanding of how spatial, environmental and genetic variation influences a key seasonal adaptation (diapause induction) in a well-developed ecological model system and will make possible future studies that explore how climatic variation affects the population dynamics and genetics of this system.

Phylogenetic incongruence and the evolutionary origins of cardenolide-resistant forms of Na(+) ,K(+) -ATPase in Danaus butterflies

Many distantly related insect species are specialized feeders of cardenolide-containing host plants such as milkweed (Asclepias spp.). Previous studies have revealed frequent, parallel substitution of a functionally important amino acid substitution (N122H) in the alpha subunit of Na(+) ,K(+) -ATPase in a number of these species. This substitution facilitates the ability of these insects to feed on their toxic hosts and sequester cardenolides for their own use in defense. Among milkweed butterflies of the genus Danaus, the previously established phylogeny for this group suggests that N122H arose independently and fixed in two distinct lineages. We reevaluate this conclusion by examining Danaus phylogenetic relationships using >400 orthologous gene sequences assembled from transcriptome data. Our results indicate that the three Danaus species known to harbor the N122H substitution are more closely related than previously thought, consistent with a single, common origin for N122H. However, we also find evidence of both incomplete lineage sorting and post-speciation genetic exchange among these butterfly species, raising the possibility of collateral evolution of cardenolide-insensitivity in this species group.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0335-z) contains supplementary material, which is available to authorized users.

Analysis of the population structure of Anaplasma phagocytophilum using multilocus sequence typing

Anaplasma phagocytophilum is a Gram-negative obligate intracellular bacterium that replicates in neutrophils. It is transmitted via tick-bite and causes febrile disease in humans and animals. Human granulocytic anaplasmosis is regarded as an emerging infectious disease in North America, Europe and Asia. However, although increasingly detected, it is still rare in Europe. Clinically apparent A. phagocytophilum infections in animals are mainly found in horses, dogs, cats, sheep and cattle. Evidence from cross-infection experiments that A. phagocytophilum isolates of distinct host origin are not uniformly infectious for heterologous hosts has led to several approaches of molecular strain characterization. Unfortunately, the results of these studies are not always easily comparable, because different gene regions and fragment lengths were investigated. Multilocus sequence typing is a widely accepted method for molecular characterization of bacteria. We here provide for the first time a universal typing method that is easily transferable between different laboratories. We validated our approach on an unprecedented large data set of almost 400 A. phagocytophilum strains from humans and animals mostly from Europe. The typability was 74% (284/383). One major clonal complex containing 177 strains was detected. However, 54% (49/90) of the sequence types were not part of a clonal complex indicating that the population structure of A. phagocytophilum is probably semiclonal. All strains from humans, dogs and horses from Europe belonged to the same clonal complex. As canine and equine granulocytic anaplasmosis occurs frequently in Europe, human granulocytic anaplasmosis is likely to be underdiagnosed in Europe. Further, wild boars and hedgehogs may serve as reservoir hosts of the disease in humans and domestic animals in Europe, because their strains belonged to the same clonal complex. In contrast, as they were only distantly related, roe deer, voles and shrews are unlikely to harbor A. phagocytophilum strains infectious for humans, domestic or farm animals.

Parallel molecular evolution in an herbivore community

Numerous insects have independently evolved the ability to feed on plants that produce toxic secondary compounds called cardenolides and can sequester these compounds for use in their defense. We surveyed the protein target for cardenolides, the alpha subunit of the sodium pump, Na(+),K(+)-ATPase (ATPα), in 14 species that feed on cardenolide-producing plants and 15 outgroups spanning three insect orders. Despite the large number of potential targets for modulating cardenolide sensitivity, amino acid substitutions associated with host-plant specialization are highly clustered, with many parallel substitutions. Additionally, we document four independent duplications of ATPα with convergent tissue-specific expression patterns. We find that unique substitutions are disproportionately associated with recent duplications relative to parallel substitutions. Together, these findings support the hypothesis that adaptation tends to take evolutionary paths that minimize negative pleiotropy.

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.

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.

Hybridization leads to host-use divergence in a polyphagous butterfly sibling species pair

Climate warming has lead to increased genetic introgression across a narrow hybrid zone separating the eastern and Canadian tiger swallowtails (Papilio glaucus and Papilio canadensis). This situation has led to the formation of an allochronically separated hybrid population with a delayed emerging phenotype or "late flight". Here, we assess how the recombination of the parental genomes that lead to this phenotype may have facilitated another major ecological shift, host-use divergence. We first contrast the ovipositional profiles of the late flight population to that of the parental species P. glaucus and P. canadensis. Subsequently we contrast the larval survival and growth of the late flight, a P. canadensis and a P. glaucus population, and a population from the northern edge of the hybrid zone on five hosts. Our results indicate that the ovipositional preference of this hybrid swarm is identical to that of the introgressing parental species, P. glaucus. Due to the absence of the preferred hosts of P. glaucus (Liriodendron tulipifera L. and Ptelea trifoliata L.) where the late flight occurs, this ovipositional pattern implies a functional specialization onto a secondary host of both parental species, Fraxinus americana L. In contrast, the larval host-use abilities represent a mixture of P. glaucus and P. canadensis, indicating divergence in larval host-use abilities has not taken place. However, high genetic variability (genetic coefficient of variation) is present for growth on F. americana in the late flight hybrid swarm and tradeoffs for larval performance on the preferred hosts of the parental species are evident; indicating a strong potential for future specialization in larval host-use abilities. This current scenario represents an instance where a shift in a major ecological trait, host use, is likely occurring as a byproduct of a shift in an unrelated trait (delayed emergence) leading to partial reproductive isolation.