Comparison of Methods for Species-Tree Inference in the Sawfly Genus Neodiprion (Hymenoptera: Diprionidae)

Conifer-feeding sawflies in the genus Neodiprion provide an excellent opportunity to investigate the origin and maintenance of barriers to reproduction, but obtaining a phylogenetic estimate for comparative studies of Neodiprion speciation has proved difficult. Specifically, nonmonophyly within and discordance between individual gene trees, both of which are common in groups that diverged recently and/or rapidly, make it impossible to infer a species tree using methods that are designed to estimate gene trees. Therefore, in this study, we estimate relationships between members of the lecontei species group using four approaches that are intended to estimate species, not gene, trees: (1) minimize deep coalescences (MDC), (2) shallowest divergences (SD), (3) Bayesian estimation of species trees (BEST), and (4) a novel approach that combines concatenation with monophyly constraints (CMC). Multiple populations are sampled for most species and all four methods incorporate this intraspecific variation into estimates of interspecific relationships. We investigate the sensitivity of each method to taxonomic sampling, and, for the BEST method, we assess the impact of prior choice on species-tree inference. We also compare species-tree estimates to one another and to a morphologically based hypothesis to identify clades that are supported by multiple analyses and lines of evidence. We find that both taxonomic sampling and method choice impact species-tree estimates and that, for these data, the BEST method is strongly influenced by and branch-length priors. We also find that the CMC method is the least sensitive to taxonomic sampling. Finally, although interspecific genetic variation is low due to the recent divergence of the lecontei group, our results to date suggest that incomplete lineage sorting and interspecific gene flow are the main factors complicating species-tree inference in Neodiprion. Based on these analyses, we propose a phylogenetic hypothesis for the lecontei group. Finally, our results suggest that, even for very challenging groups like Neodiprion, an underlying species-tree signal can be extracted from multi-locus data as long as intraspecific variation is adequately sampled and methods that focus on the estimation of species trees are used. [Bayesian estimation of species trees (BEST); concatenation with monophyly constraints (CMC); gene-tree discordance; hybridization; introgression; lineage sorting; minimize deep coalescences (MDC); shallowest divergences (SD).]Organismic and Evolutionary Biolog

Oviposition Traits Generate Extrinsic Postzygotic Isolation Between Two Pine Sawfly Species

Background: Although empirical data indicate that ecological speciation is prevalent in nature, the relative importance of different forms of reproductive isolation and the traits generating reproductive isolation remain unclear. To address these questions, we examined a pair of ecologically divergent pine-sawfly species: while Neodiprion pinetum specializes on a thin-needled pine (Pinus strobus), N. lecontei utilizes thicker-needled pines. We hypothesized that extrinsic postzygotic isolation is generated by oviposition traits. To test this hypothesis, we assayed ovipositor morphology, oviposition behavior, and host-dependent oviposition success in both species and in F1 and backcross females. Results: Compared to N. lecontei, N. pinetum females preferred P. strobus more strongly, had smaller ovipositors, and laid fewer eggs per needle. Additionally, we observed host- and trait-dependent reductions in oviposition success in F1 and backcross females. Hybrid females that had pinetum-like host preference (P. strobus) and lecontei-like oviposition traits (morphology and egg pattern) fared especially poorly. Conclusions: Together, these data indicate that maladaptive combinations of oviposition traits in hybrids contribute to extrinsic postzygotic isolation between N. lecontei and N. pinetum, suggesting that oviposition traits may be an important driver of divergence in phytophagous insects

The population genetics of crypsis in vertebrates: recent insights from mice, hares, and lizards

By combining well-established population genetic theory with high-throughput sequencing data from natural populations, major strides have recently been made in understanding how, why, and when vertebrate populations evolve crypsis. Here, we focus on background matching, a particular facet of crypsis that involves the ability of an organism to conceal itself through matching its color to the surrounding environment. While interesting in and of itself, the study of this phenotype has also provided fruitful population genetic insights into the interplay of strong positive selection with other evolutionary processes. Specifically, and predicated upon the findings of previous candidate gene association studies, a primary focus of this recent literature involves the realization that the inference of selection from DNA sequence data first requires a robust model of population demography in order to identify genomic regions which do not conform to neutral expectations. Moreover, these demographic estimates provide crucial information about the origin and timing of the onset of selective pressures associated with, for example, the colonization of a novel environment. Furthermore, such inference has revealed crypsis to be a particularly useful phenotype for investigating the interplay of migration and selection-with examples of gene flow constraining rates of adaptation, or alternatively providing the genetic variants that may ultimately sweep through the population. Here, we evaluate the underlying evidence, review the strengths and weaknesses of the many population genetic methodologies used in these studies, and discuss how these insights have aided our general understanding of the evolutionary process

Variation in DNA Methylation Is Not Consistently Reflected by Sociality in Hymenoptera

Changes in gene regulation that underlie phenotypic evolution can be encoded directly in the DNA sequence or mediated by chromatin modifications such as DNA methylation. It has been hypothesized that the evolution of eusocial division of labor is associated with enhanced gene regulatory potential, which may include expansions in DNA methylation in the genomes of Hymenoptera (bees, ants, wasps, and sawflies). Recently, this hypothesis garnered support from analyses of a commonly used metric to estimate DNA methylation in silico, CpG content. Here, we test this hypothesis using direct, nucleotide-level measures of DNA methylation across nine species of Hymenoptera. In doing so, we generated new DNA methylomes for three species of interest, including one solitary and one facultatively eusocial halictid bee and a sawfly. We demonstrate that the strength of correlation between CpG content and DNA methylation varies widely among hymenopteran taxa, highlighting shortcomings in the utility of CpG content as a proxy for DNA methylation in comparative studies of taxa with sparse DNA methylomes. We observed strikingly high levels of DNA methylation in the sawfly relative to other investigated hymenopterans. Analyses of molecular evolution suggest the relatively distinct sawfly DNA methylome may be associated with positive selection on functional DNMT3 domains. Sawflies are an outgroup to all ants, bees, and wasps, and no sawfly species are eusocial. We find no evidence that either global expansions or variation within individual ortholog groups in DNA methylation are consistently associated with the evolution of social behavior

Gregariousness Does Not Vary with Geography, Developmental Stage, or Group Relatedness in Feeding Redheaded Pine Sawfly Larvae

Aggregations are widespread across the animal kingdom, yet the underlying proximate and ultimate causes are still largely unknown. An ideal system to investigate this simple, social behavior is the pine sawfly genus Neodiprion, which is experimentally tractable and exhibits interspecific variation in larval gregariousness. To assess intraspecific variation in this trait, we characterized aggregative tendency within a single widespread species, the redheaded pine sawfly (N. lecontei). To do so, we developed a quantitative assay in which we measured interindividual distances over a 90-min video. This assay revealed minimal behavioral differences: (1) between early-feeding and late-feeding larval instars, (2) among larvae derived from different latitudes, and (3) between groups composed of kin and those composed of nonkin. Together, these results suggest that, during the larval feeding period, the benefits individuals derive from aggregating outweigh the costs and that this cost-to-benefit ratio does not vary dramatically across space (geography) or ontogeny (developmental stage). In contrast to the feeding larvae, our assay revealed a striking reduction in gregariousness following the final larval molt in N. lecontei. We also found some intriguing interspecific variation: While N. lecontei and N. maurus feeding larvae exhibit significant aggregative tendencies, feeding N. compar larvae do not aggregate at all. These results set the stage for future work investigating the proximate and ultimate mechanisms underlying developmental and interspecific variation in larval gregariousness across Neodiprion

Eda haplotypes in three-spined stickleback are associated with variation in immune gene expression

Haplotypes underlying local adaptation and speciation are predicted to have numerous phenotypic effects, but few genes involved have been identified, with much work to date concentrating on visible, morphological, phenotypes. The link between genes controlling these adaptive morphological phenotypes and the immune system has seldom been investigated, even though changes in the immune system could have profound adaptive consequences. The Eda gene in three-spined stickleback is one of the best studied major adaptation genes; it directly controls bony plate architecture and has been associated with additional aspects of adaptation to freshwater. Here, we exposed F2 hybrids, used to separate Eda genotype from genetic background, to contrasting conditions in semi-natural enclosures. We demonstrate an association between the Eda haplotype block and the expression pattern of key immune system genes. Furthermore, low plated fish grew less and experienced higher burdens of a common ectoparasite with fitness consequences. Little is currently known about the role of the immune system in facilitating adaptation to novel environments, but this study provides an indication of its potential importance

Challenges in Estimating Insecticide Selection Pressures from Mosquito Field Data

Insecticide resistance has the potential to compromise the enormous effort put into the control of dengue and malaria vector populations. It is therefore important to quantify the amount of selection acting on resistance alleles, their contributions to fitness in heterozygotes (dominance) and their initial frequencies, as a means to predict the rate of spread of resistance in natural populations. We investigate practical problems of obtaining such estimates, with particular emphasis on Mexican populations of the dengue vector Aedes aegypti. Selection and dominance coefficients can be estimated by fitting genetic models to field data using maximum likelihood (ML) methodology. This methodology, although widely used, makes many assumptions so we investigated how well such models perform when data are sparse or when spatial and temporal heterogeneity occur. As expected, ML methodologies reliably estimated selection and dominance coefficients under idealised conditions but it was difficult to recover the true values when datasets were sparse during the time that resistance alleles increased in frequency, or when spatial and temporal heterogeneity occurred. We analysed published data on pyrethroid resistance in Mexico that consists of the frequency of a Ile1,016 mutation. The estimates for selection coefficient and initial allele frequency on the field dataset were in the expected range, dominance coefficient points to incomplete dominance as observed in the laboratory, although these estimates are accompanied by strong caveats about possible impact of spatial and temporal heterogeneity in selection

GB virus-C – a virus without a disease: We cannot give it chronic fatigue syndrome

BACKGROUND: Chronic fatigue syndrome (CFS) is an illness in search of an infectious etiology. GB virus-C (GBV-C) virus is a flavivirus with cell tropism and host defense induction qualities compatible with a role in producing the syndrome. The GBV-C genome is detectable in 4% of the population and 12% of the population is seropositive. The present study evaluated the association between infection with GBV and CFS. METHODS: We used a commercial EIA to detect antibodies against the GBV-C E2 protein and a quantitative real-time RT-PCR assay to detect active GBV-C infection. Sera were from a case control study of CFS in Atlanta, Georgia. The Fisher's exact two-tailed test was used for statistical analysis. RESULTS: Two of 12 CFS patients and one of 21 controls were seropositive for prior GBV-C infection and one control had viral RNA detected, indicating active infection. The results are not statistically different. CONCLUSION: We found no evidence that active or past infection with GBV is associated with CFS

Convergent Evolution of Escape from Hepaciviral Antagonism in Primates

Escape from antagonism by hepatitis C and related viruses has repeatedly evolved in antiviral factor MAVS via convergent evolution, revealing an ancient history of previous viral encounters in primates

Identification of GBV-D, a Novel GB-like Flavivirus from Old World Frugivorous Bats (Pteropus giganteus) in Bangladesh

Bats are reservoirs for a wide range of zoonotic agents including lyssa-, henipah-, SARS-like corona-, Marburg-, Ebola-, and astroviruses. In an effort to survey for the presence of other infectious agents, known and unknown, we screened sera from 16 Pteropus giganteus bats from Faridpur, Bangladesh, using high-throughput pyrosequencing. Sequence analyses indicated the presence of a previously undescribed virus that has approximately 50% identity at the amino acid level to GB virus A and C (GBV-A and -C). Viral nucleic acid was present in 5 of 98 sera (5%) from a single colony of free-ranging bats. Infection was not associated with evidence of hepatitis or hepatic dysfunction. Phylogenetic analysis indicates that this first GBV-like flavivirus reported in bats constitutes a distinct species within the Flaviviridae family and is ancestral to the GBV-A and -C virus clades