@misc{6060,
  author       = {Vicoso, Beatriz},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Supplementary data for "Sex-biased gene expression and dosage compensation on the Artemia franciscana Z-chromosome" (Huylman, Toups et al., 2019). }},
  doi          = {10.15479/AT:ISTA:6060},
  year         = {2019},
}

@article{6089,
  abstract     = {Pleiotropy is the well-established idea that a single mutation affects multiple phenotypes. If a mutation has opposite effects on fitness when expressed in different contexts, then genetic conflict arises. Pleiotropic conflict is expected to reduce the efficacy of selection by limiting the fixation of beneficial mutations through adaptation, and the removal of deleterious mutations through purifying selection. Although this has been widely discussed, in particular in the context of a putative “gender load,” it has yet to be systematically quantified. In this work, we empirically estimate to which extent different pleiotropic regimes impede the efficacy of selection in Drosophila melanogaster. We use whole-genome polymorphism data from a single African population and divergence data from D. simulans to estimate the fraction of adaptive fixations (α), the rate of adaptation (ωA), and the direction of selection (DoS). After controlling for confounding covariates, we find that the different pleiotropic regimes have a relatively small, but significant, effect on selection efficacy. Specifically, our results suggest that pleiotropic sexual antagonism may restrict the efficacy of selection, but that this conflict can be resolved by limiting the expression of genes to the sex where they are beneficial. Intermediate levels of pleiotropy across tissues and life stages can also lead to maladaptation in D. melanogaster, due to inefficient purifying selection combined with low frequency of mutations that confer a selective advantage. Thus, our study highlights the need to consider the efficacy of selection in the context of antagonistic pleiotropy, and of genetic conflict in general.},
  author       = {Fraisse, Christelle and Puixeu Sala, Gemma and Vicoso, Beatriz},
  issn         = {1537-1719},
  journal      = {Molecular biology and evolution},
  number       = {3},
  pages        = {500--515},
  publisher    = {Oxford University Press},
  title        = {{Pleiotropy modulates the efficacy of selection in drosophila melanogaster}},
  doi          = {10.1093/molbev/msy246},
  volume       = {36},
  year         = {2019},
}

@article{6418,
  abstract     = {Males and females of Artemia franciscana, a crustacean commonly used in the aquarium trade, are highly dimorphic. Sex is determined by a pair of ZW chromosomes, but the nature and extent of differentiation of these chromosomes is unknown. Here, we characterize the Z chromosome by detecting genomic regions that show lower genomic coverage in female than in male samples, and regions that harbor an excess of female-specific SNPs. We detect many Z-specific genes, which no longer have homologs on the W, but also Z-linked genes that appear to have diverged very recently from their existing W-linked homolog. We assess patterns of male and female expression in two tissues with extensive morphological dimorphism, gonads, and heads. In agreement with their morphology, sex-biased expression is common in both tissues. Interestingly, the Z chromosome is not enriched for sex-biased genes, and seems to in fact have a mechanism of dosage compensation that leads to equal expression in males and in females. Both of these patterns are contrary to most ZW systems studied so far, making A. franciscana an excellent model for investigating the interplay between the evolution of sexual dimorphism and dosage compensation, as well as Z chromosome evolution in general.},
  author       = {Huylmans, Ann K and Toups, Melissa A and Macon, Ariana and Gammerdinger, William J and Vicoso, Beatriz},
  issn         = {1759-6653},
  journal      = {Genome biology and evolution},
  number       = {4},
  pages        = {1033--1044},
  publisher    = {Oxford University Press},
  title        = {{Sex-biased gene expression and dosage compensation on the Artemia franciscana Z-chromosome}},
  doi          = {10.1093/gbe/evz053},
  volume       = {11},
  year         = {2019},
}

@article{6467,
  abstract     = {Fitness interactions between mutations can influence a population’s evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA (small nucleolar RNA). Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations.},
  author       = {Fraisse, Christelle and Welch, John J.},
  issn         = {1744957X},
  journal      = {Biology Letters},
  number       = {4},
  publisher    = {Royal Society of London},
  title        = {{The distribution of epistasis on simple fitness landscapes}},
  doi          = {10.1098/rsbl.2018.0881},
  volume       = {15},
  year         = {2019},
}

@misc{9803,
  abstract     = {Understanding the mechanisms causing phenotypic differences between females and males has long fascinated evolutionary biologists. An extensive literature exists on animal sexual dimorphism but less is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life-cycle dynamics. Here, we investigate patterns of genetically-based sexual dimorphism in vegetative and reproductive traits of a wind-pollinated dioecious plant, Rumex hastatulus, across three life-cycle stages using open-pollinated families from 30 populations spanning the geographic range and chromosomal variation (XY and XY1Y2) of the species. The direction and degree of sexual dimorphism was highly variable among populations and life-cycle stages. Sex-specific differences in reproductive function explained a significant amount of temporal change in sexual dimorphism. For several traits, geographical variation in sexual dimorphism was associated with bioclimatic parameters, likely due to the differential responses of the sexes to climate. We found no systematic differences in sexual dimorphism between chromosome races. Sex-specific trait differences in dioecious plants largely result from a balance between sexual and natural selection on resource allocation. Our results indicate that abiotic factors associated with geographical context also play a role in modifying sexual dimorphism during the plant life cycle.},
  author       = {Puixeu Sala, Gemma and Pickup, Melinda and Field, David and Barrett, Spencer C.H.},
  publisher    = {Dryad},
  title        = {{Data from: Variation in sexual dimorphism in a wind-pollinated plant: the influence of geographical context and life-cycle dynamics}},
  doi          = {10.5061/dryad.n1701c9},
  year         = {2019},
}

@article{190,
  abstract     = {The German cockroach, Blattella germanica, is a worldwide pest that infests buildings, including homes, restaurants, and hospitals, often living in unsanitary conditions. As a disease vector and producer of allergens, this species has major health and economic impacts on humans. Factors contributing to the success of the German cockroach include its resistance to a broad range of insecticides, immunity to many pathogens, and its ability, as an extreme generalist omnivore, to survive on most food sources. The recently published genome shows that B. germanica has an exceptionally high number of protein coding genes. In this study, we investigate the functions of the 93 significantly expanded gene families with the aim to better understand the success of B. germanica as a major pest despite such inhospitable conditions. We find major expansions in gene families with functions related to the detoxification of insecticides and allelochemicals, defense against pathogens, digestion, sensory perception, and gene regulation. These expansions might have allowed B. germanica to develop multiple resistance mechanisms to insecticides and pathogens, and enabled a broad, flexible diet, thus explaining its success in unsanitary conditions and under recurrent chemical control. The findings and resources presented here provide insights for better understanding molecular mechanisms that will facilitate more effective cockroach control.},
  author       = {Harrison, Mark and Arning, Nicolas and Kremer, Lucas and Ylla, Guillem and Belles, Xavier and Bornberg Bauer, Erich and Huylmans, Ann K and Jongepier, Evelien and Puilachs, Maria and Richards, Stephen and Schal, Coby},
  journal      = {Journal of Experimental Zoology Part B: Molecular and Developmental Evolution},
  pages        = {254--264},
  publisher    = {Wiley},
  title        = {{Expansions of key protein families in the German cockroach highlight the molecular basis of its remarkable success as a global indoor pest}},
  doi          = {10.1002/jez.b.22824},
  volume       = {330},
  year         = {2018},
}

@article{199,
  abstract     = {Sex-biased genes are central to the study of sexual selection, sexual antagonism, and sex chromosome evolution. We describe a comprehensive de novo assembled transcriptome in the common frog Rana temporaria based on five developmental stages and three adult tissues from both sexes, obtained from a population with karyotypically homomorphic but genetically differentiated sex chromosomes. This allows the study of sex-biased gene expression throughout development, and its effect on the rate of gene evolution while accounting for pleiotropic expression, which is known to negatively correlate with the evolutionary rate. Overall, sex-biased genes had little overlap among developmental stages and adult tissues. Late developmental stages and gonad tissues had the highest numbers of stage-or tissue-specific genes. We find that pleiotropic gene expression is a better predictor than sex bias for the evolutionary rate of genes, though it often interacts with sex bias. Although genetically differentiated, the sex chromosomes were not enriched in sex-biased genes, possibly due to a very recent arrest of XY recombination. These results extend our understanding of the developmental dynamics, tissue specificity, and genomic localization of sex-biased genes.},
  author       = {Ma, Wen and Veltsos, Paris and Toups, Melissa A and Rodrigues, Nicolas and Sermier, Roberto and Jeffries, Daniel and Perrin, Nicolas},
  journal      = {Genes},
  number       = {6},
  publisher    = {MDPI AG},
  title        = {{Tissue specificity and dynamics of sex biased gene expression in a common frog population with differentiated, yet homomorphic, sex chromosomes}},
  doi          = {10.3390/genes9060294},
  volume       = {9},
  year         = {2018},
}

@article{542,
  abstract     = {The t-haplotype, a mouse meiotic driver found on chromosome 17, has been a model for autosomal segregation distortion for close to a century, but several questions remain regarding its biology and evolutionary history. A recently published set of population genomics resources for wild mice includes several individuals heterozygous for the t-haplotype, which we use to characterize this selfish element at the genomic and transcriptomic level. Our results show that large sections of the t-haplotype have been replaced by standard homologous sequences, possibly due to occasional events of recombination, and that this complicates the inference of its history. As expected for a long genomic segment of very low recombination, the t-haplotype carries an excess of fixed nonsynonymous mutations compared to the standard chromosome. This excess is stronger for regions that have not undergone recent recombination, suggesting that occasional gene flow between the t and the standard chromosome may provide a mechanism to regenerate coding sequences that have accumulated deleterious mutations. Finally, we find that t-complex genes with altered expression largely overlap with deleted or amplified regions, and that carrying a t-haplotype alters the testis expression of genes outside of the t-complex, providing new leads into the pathways involved in the biology of this segregation distorter.},
  author       = {Kelemen, Réka K and Vicoso, Beatriz},
  journal      = {Genetics},
  number       = {1},
  pages        = {365 -- 375},
  publisher    = {Genetics Society of America},
  title        = {{Complex history and differentiation patterns of the t-haplotype, a mouse meiotic driver}},
  doi          = {10.1534/genetics.117.300513},
  volume       = {208},
  year         = {2018},
}

@misc{5586,
  abstract     = {Input files and scripts from "Evolution of gene dosage on the Z-chromosome of schistosome parasites" by Picard M.A.L., et al (2018).},
  author       = {Vicoso, Beatriz},
  keywords     = {schistosoma, Z-chromosome, gene expression},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Input files and scripts from "Evolution of gene dosage on the Z-chromosome of schistosome parasites" by Picard M.A.L., et al (2018)}},
  doi          = {10.15479/AT:ISTA:109},
  year         = {2018},
}

@misc{5757,
  abstract     = {File S1. Variant Calling Format file of the ingroup: 197 haploid sequences of D. melanogaster from Zambia (Africa) aligned to the D. melanogaster 5.57 reference genome.

File S2. Variant Calling Format file of the outgroup: 1 haploid sequence of D. simulans aligned to the D. melanogaster 5.57 reference genome.

File S3. Annotations of each transcript in coding regions with SNPeff: Ps (# of synonymous polymorphic sites); Pn (# of non-synonymous polymorphic sites); Ds (# of synonymous divergent sites); Dn (# of non-synonymous divergent sites); DoS; ⍺ MK . All variants were included.

File S4. Annotations of each transcript in non-coding regions with SNPeff: Ps (# of synonymous polymorphic sites); Pu (# of UTR polymorphic sites); Ds (# of synonymous divergent sites); Du (# of UTR divergent sites); DoS; ⍺ MK . All variants were included.

File S5. Annotations of each transcript in coding regions with SNPGenie: Ps (# of synonymous polymorphic sites); πs (synonymous diversity); Ss_p (total # of synonymous sites in the polymorphism data); Pn (# of non-synonymous polymorphic sites); πn (non-synonymous diversity); Sn_p (total # of non-synonymous sites in the polymorphism data); Ds (# of synonymous divergent sites); ks (synonymous evolutionary rate); Ss_d (total # of synonymous sites in the divergence data); Dn (# of non-synonymous divergent sites); kn (non-synonymous evolutionary rate); Sn_d (total # of non-
synonymous sites in the divergence data); DoS; ⍺ MK . All variants were included.

File S6. Gene expression values (RPKM summed over all transcripts) for each sample. Values were quantile-normalized across all samples.

File S7. Final dataset with all covariates, ⍺ MK , ωA MK and DoS for coding sites, excluding variants below 5% frequency.

File S8. Final dataset with all covariates, ⍺ MK , ωA MK and DoS for non-coding sites, excluding variants below 5%
frequency.

File S9. Final dataset with all covariates, ⍺ EWK , ωA EWK and deleterious SFS for coding sites obtained with the Eyre-Walker and Keightley method on binned data and using all variants.},
  author       = {Fraisse, Christelle},
  keywords     = {(mal)adaptation, pleiotropy, selective constraint, evo-devo, gene expression, Drosophila melanogaster},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Supplementary Files for "Pleiotropy modulates the efficacy of selection in Drosophila melanogaster"}},
  doi          = {10.15479/at:ista:/5757},
  year         = {2018},
}

@article{5989,
  abstract     = {Schistosomes are the causative agents of schistosomiasis, a neglected tropical disease affecting over 230 million people worldwide.Additionally to their major impact on human health, they are also models of choice in evolutionary biology. These parasitic flatwormsare unique among the common hermaphroditic trematodes as they have separate sexes. This so-called “evolutionary scandal”displays a female heterogametic genetic sex-determination system (ZZ males and ZW females), as well as a pronounced adult sexualdimorphism. These phenotypic differences are determined by a shared set of genes in both sexes, potentially leading to intralocussexual conflicts. To resolve these conflicts in sexually selected traits, molecular mechanisms such as sex-biased gene expression couldoccur, but parent-of-origin gene expression also provides an alternative. In this work we investigated the latter mechanism, that is,genes expressed preferentially from either the maternal or the paternal allele, inSchistosoma mansonispecies. To this end, tran-scriptomes from male and female hybrid adults obtained by strain crosses were sequenced. Strain-specific single nucleotide poly-morphism (SNP) markers allowed us to discriminate the parental origin, while reciprocal crosses helped to differentiate parentalexpression from strain-specific expression. We identified genes containing SNPs expressed in a parent-of-origin manner consistentwith paternal and maternal imprints. Although the majority of the SNPs was identified in mitochondrial and Z-specific loci, theremaining SNPs found in male and female transcriptomes were situated in genes that have the potential to explain sexual differencesin schistosome parasites. Furthermore, we identified and validated four new Z-specific scaffolds.},
  author       = {Kincaid-Smith, Julien and Picard, Marion A L and Cosseau, Céline and Boissier, Jérôme and Severac, Dany and Grunau, Christoph and Toulza, Eve},
  issn         = {1759-6653},
  journal      = {Genome Biology and Evolution},
  number       = {3},
  pages        = {840--856},
  publisher    = {Oxford University Press},
  title        = {{Parent-of-Origin-Dependent Gene Expression in Male and Female Schistosome Parasites}},
  doi          = {10.1093/gbe/evy037},
  volume       = {10},
  year         = {2018},
}

@article{63,
  abstract     = {African cichlids display a remarkable assortment of jaw morphologies, pigmentation patterns, and mating behaviors. In addition to this previously documented diversity, recent studies have documented a rich diversity of sex chromosomes within these fishes. Here we review the known sex-determination network within vertebrates, and the extraordinary number of sex chromosomes systems segregating in African cichlids. We also propose a model for understanding the unusual number of sex chromosome systems within this clade.},
  author       = {Gammerdinger, William J and Kocher, Thomas},
  journal      = {Genes},
  number       = {10},
  publisher    = {MDPI AG},
  title        = {{Unusual diversity of sex chromosomes in African cichlid fishes}},
  doi          = {10.3390/genes9100480},
  volume       = {9},
  year         = {2018},
}

@article{131,
  abstract     = {XY systems usually show chromosome-wide compensation of X-linked genes, while in many ZW systems, compensation is restricted to a minority of dosage-sensitive genes. Why such differences arose is still unclear. Here, we combine comparative genomics, transcriptomics and proteomics to obtain a complete overview of the evolution of gene dosage on the Z-chromosome of Schistosoma parasites. We compare the Z-chromosome gene content of African (Schistosoma mansoni and S. haematobium) and Asian (S. japonicum) schistosomes and describe lineage-specific evolutionary strata. We use these to assess gene expression evolution following sex-linkage. The resulting patterns suggest a reduction in expression of Z-linked genes in females, combined with upregulation of the Z in both sexes, in line with the first step of Ohno’s classic model of dosage compensation evolution. Quantitative proteomics suggest that post-transcriptional mechanisms do not play a major role in balancing the expression of Z-linked genes. },
  author       = {Picard, Marion A and Cosseau, Celine and Ferré, Sabrina and Quack, Thomas and Grevelding, Christoph and Couté, Yohann and Vicoso, Beatriz},
  journal      = {eLife},
  publisher    = {eLife Sciences Publications},
  title        = {{Evolution of gene dosage on the Z-chromosome of schistosome parasites}},
  doi          = {10.7554/eLife.35684},
  volume       = {7},
  year         = {2018},
}

@article{139,
  abstract     = {Genome-scale diversity data are increasingly available in a variety of biological systems, and can be used to reconstruct the past evolutionary history of species divergence. However, extracting the full demographic information from these data is not trivial, and requires inferential methods that account for the diversity of coalescent histories throughout the genome. Here, we evaluate the potential and limitations of one such approach. We reexamine a well-known system of mussel sister species, using the joint site frequency spectrum (jSFS) of synonymousmutations computed either fromexome capture or RNA-seq, in an Approximate Bayesian Computation (ABC) framework. We first assess the best sampling strategy (number of: individuals, loci, and bins in the jSFS), and show that model selection is robust to variation in the number of individuals and loci. In contrast, different binning choices when summarizing the jSFS, strongly affect the results: including classes of low and high frequency shared polymorphisms can more effectively reveal recent migration events. We then take advantage of the flexibility of ABC to compare more realistic models of speciation, including variation in migration rates through time (i.e., periodic connectivity) and across genes (i.e., genome-wide heterogeneity in migration rates). We show that these models were consistently selected as the most probable, suggesting that mussels have experienced a complex history of gene flow during divergence and that the species boundary is semi-permeable. Our work provides a comprehensive evaluation of ABC demographic inference in mussels based on the coding jSFS, and supplies guidelines for employing different sequencing techniques and sampling strategies. We emphasize, perhaps surprisingly, that inferences are less limited by the volume of data, than by the way in which they are analyzed.},
  author       = {Fraisse, Christelle and Roux, Camille and Gagnaire, Pierre and Romiguier, Jonathan and Faivre, Nicolas and Welch, John and Bierne, Nicolas},
  journal      = {PeerJ},
  number       = {7},
  publisher    = {PeerJ},
  title        = {{The divergence history of European blue mussel species reconstructed from Approximate Bayesian Computation: The effects of sequencing techniques and sampling strategies}},
  doi          = {10.7717/peerj.5198},
  volume       = {2018},
  year         = {2018},
}

@misc{9841,
  abstract     = {Around 150 million years ago, eusocial termites evolved from within the cockroaches, 50 million years before eusocial Hymenoptera, such as bees and ants, appeared. Here, we report the 2-Gb genome of the German cockroach, Blattella germanica, and the 1.3-Gb genome of the drywood termite Cryptotermes secundus. We show evolutionary signatures of termite eusociality by comparing the genomes and transcriptomes of three termites and the cockroach against the background of 16 other eusocial and non-eusocial insects. Dramatic adaptive changes in genes underlying the production and perception of pheromones confirm the importance of chemical communication in the termites. These are accompanied by major changes in gene regulation and the molecular evolution of caste determination. Many of these results parallel molecular mechanisms of eusocial evolution in Hymenoptera. However, the specific solutions are remarkably different, thus revealing a striking case of convergence in one of the major evolutionary transitions in biological complexity.},
  author       = {Harrison, Mark C. and Jongepier, Evelien and Robertson, Hugh M. and Arning, Nicolas and Bitard-Feildel, Tristan and Chao, Hsu and Childers, Christopher P. and Dinh, Huyen and Doddapaneni, Harshavardhan and Dugan, Shannon and Gowin, Johannes and Greiner, Carolin and Han, Yi and Hu, Haofu and Hughes, Daniel S. T. and Huylmans, Ann K and Kemena, Carsten and Kremer, Lukas P. M. and Lee, Sandra L. and Lopez-Ezquerra, Alberto and Mallet, Ludovic and Monroy-Kuhn, Jose M. and Moser, Annabell and Murali, Shwetha C. and Muzny, Donna M. and Otani, Saria and Piulachs, Maria-Dolors and Poelchau, Monica and Qu, Jiaxin and Schaub, Florentine and Wada-Katsumata, Ayako and Worley, Kim C. and Xie, Qiaolin and Ylla, Guillem and Poulsen, Michael and Gibbs, Richard A. and Schal, Coby and Richards, Stephen and Belles, Xavier and Korb, Judith and Bornberg-Bauer, Erich},
  publisher    = {Dryad},
  title        = {{Data from: Hemimetabolous genomes reveal molecular basis of termite eusociality}},
  doi          = {10.5061/dryad.51d4r},
  year         = {2018},
}

@article{9915,
  abstract     = {The evolution of assortative mating is a key part of the speciation process. Stronger assortment, or greater divergence in mating traits, between species pairs with overlapping ranges is commonly observed, but possible causes of this pattern of reproductive character displacement are difficult to distinguish. We use a multidisciplinary approach to provide a rare example where it is possible to distinguish among hypotheses concerning the evolution of reproductive character displacement. We build on an earlier comparative analysis that illustrated a strong pattern of greater divergence in penis form between pairs of sister species with overlapping ranges than between allopatric sister-species pairs, in a large clade of marine gastropods (Littorinidae). We investigate both assortative mating and divergence in male genitalia in one of the sister-species pairs, discriminating among three contrasting processes each of which can generate a pattern of reproductive character displacement: reinforcement, reproductive interference and the Templeton effect. We demonstrate reproductive character displacement in assortative mating, but not in genital form between this pair of sister species and use demographic models to distinguish among the different processes. Our results support a model with no gene flow since secondary contact and thus favor reproductive interference as the cause of reproductive character displacement for mate choice, rather than reinforcement. High gene flow within species argues against the Templeton effect. Secondary contact appears to have had little impact on genital divergence.},
  author       = {Hollander, Johan and Montaño-Rendón, Mauricio and Bianco, Giuseppe and Yang, Xi and Westram, Anja M and Duvaux, Ludovic and Reid, David G. and Butlin, Roger K.},
  issn         = {2056-3744},
  journal      = {Evolution Letters},
  number       = {6},
  pages        = {557--566},
  publisher    = {Wiley},
  title        = {{Are assortative mating and genital divergence driven by reinforcement?}},
  doi          = {10.1002/evl3.85},
  volume       = {2},
  year         = {2018},
}

@article{9917,
  abstract     = {Adaptive divergence and speciation may happen despite opposition by gene flow. Identifying the genomic basis underlying divergence with gene flow is a major task in evolutionary genomics. Most approaches (e.g., outlier scans) focus on genomic regions of high differentiation. However, not all genomic architectures potentially underlying divergence are expected to show extreme differentiation. Here, we develop an approach that combines hybrid zone analysis (i.e., focuses on spatial patterns of allele frequency change) with system-specific simulations to identify loci inconsistent with neutral evolution. We apply this to a genome-wide SNP set from an ideally suited study organism, the intertidal snail Littorina saxatilis, which shows primary divergence between ecotypes associated with different shore habitats. We detect many SNPs with clinal patterns, most of which are consistent with neutrality. Among non-neutral SNPs, most are located within three large putative inversions differentiating ecotypes. Many non-neutral SNPs show relatively low levels of differentiation. We discuss potential reasons for this pattern, including loose linkage to selected variants, polygenic adaptation and a component of balancing selection within populations (which may be expected for inversions). Our work is in line with theory predicting a role for inversions in divergence, and emphasizes that genomic regions contributing to divergence may not always be accessible with methods purely based on allele frequency differences. These conclusions call for approaches that take spatial patterns of allele frequency change into account in other systems.},
  author       = {Westram, Anja M and Rafajlović, Marina and Chaube, Pragya and Faria, Rui and Larsson, Tomas and Panova, Marina and Ravinet, Mark and Blomberg, Anders and Mehlig, Bernhard and Johannesson, Kerstin and Butlin, Roger},
  issn         = {2056-3744},
  journal      = {Evolution Letters},
  number       = {4},
  pages        = {297--309},
  publisher    = {Wiley},
  title        = {{Clines on the seashore: The genomic architecture underlying rapid divergence in the face of gene flow}},
  doi          = {10.1002/evl3.74},
  volume       = {2},
  year         = {2018},
}

@misc{9929,
  abstract     = {The evolution of assortative mating is a key part of the speciation process. Stronger assortment, or greater divergence in mating traits, between species pairs with overlapping ranges is commonly observed, but possible causes of this pattern of reproductive character displacement are difficult to distinguish. We use a multidisciplinary approach to provide a rare example where it is possible to distinguish among hypotheses concerning the evolution of reproductive character displacement. We build on an earlier comparative analysis that illustrated a strong pattern of greater divergence in penis form between pairs of sister species with overlapping ranges than between allopatric sister-species pairs, in a large clade of marine gastropods (Littorinidae). We investigate both assortative mating and divergence in male genitalia in one of the sister-species pairs, discriminating among three contrasting processes each of which can generate a pattern of reproductive character displacement: reinforcement, reproductive interference and the Templeton effect. We demonstrate reproductive character displacement in assortative mating, but not in genital form between this pair of sister species and use demographic models to distinguish among the different processes. Our results support a model with no gene flow since secondary contact and thus favour reproductive interference as the cause of reproductive character displacement for mate choice, rather than reinforcement. High gene flow within species argues against the Templeton effect. Secondary contact appears to have had little impact on genital divergence.},
  author       = {Hollander, Johan and Montaño-Rendón, Mauricio and Bianco, Giuseppe and Yang, Xi and Westram, Anja M and Duvaux, Ludovic and Reid, David G. and Butlin, Roger K.},
  publisher    = {Dryad},
  title        = {{Data from: Are assortative mating and genital divergence driven by reinforcement?}},
  doi          = {10.5061/dryad.51sd2p5},
  year         = {2018},
}

@misc{9930,
  abstract     = {Adaptive divergence and speciation may happen despite opposition by gene flow. Identifying the genomic basis underlying divergence with gene flow is a major task in evolutionary genomics. Most approaches (e.g. outlier scans) focus on genomic regions of high differentiation. However, not all genomic architectures potentially underlying divergence are expected to show extreme differentiation. Here, we develop an approach that combines hybrid zone analysis (i.e. focuses on spatial patterns of allele frequency change) with system-specific simulations to identify loci inconsistent with neutral evolution. We apply this to a genome-wide SNP set from an ideally-suited study organism, the intertidal snail Littorina saxatilis, which shows primary divergence between ecotypes associated with different shore habitats. We detect many SNPs with clinal patterns, most of which are consistent with neutrality. Among non-neutral SNPs, most are located within three large putative inversions differentiating ecotypes. Many non-neutral SNPs show relatively low levels of differentiation. We discuss potential reasons for this pattern, including loose linkage to selected variants, polygenic adaptation and a component of balancing selection within populations (which may be expected for inversions). Our work is in line with theory predicting a role for inversions in divergence, and emphasises that genomic regions contributing to divergence may not always be accessible with methods purely based on allele frequency differences. These conclusions call for approaches that take spatial patterns of allele frequency change into account in other systems.},
  author       = {Westram, Anja M and Rafajlović, Marina and Chaube, Pragya and Faria, Rui and Larsson, Tomas and Panova, Marina and Ravinet, Mark and Blomberg, Anders and Mehlig, Bernhard and Johannesson, Kerstin and Butlin, Roger},
  publisher    = {Dryad},
  title        = {{Data from: Clines on the seashore: the genomic architecture underlying rapid divergence in the face of gene flow}},
  doi          = {10.5061/dryad.bp25b65},
  year         = {2018},
}

@article{448,
  abstract     = {Around 150 million years ago, eusocial termites evolved from within the cockroaches, 50 million years before eusocial Hymenoptera, such as bees and ants, appeared. Here, we report the 2-Gb genome of the German cockroach, Blattella germanica, and the 1.3-Gb genome of the drywood termite Cryptotermes secundus. We show evolutionary signatures of termite eusociality by comparing the genomes and transcriptomes of three termites and the cockroach against the background of 16 other eusocial and non-eusocial insects. Dramatic adaptive changes in genes underlying the production and perception of pheromones confirm the importance of chemical communication in the termites. These are accompanied by major changes in gene regulation and the molecular evolution of caste determination. Many of these results parallel molecular mechanisms of eusocial evolution in Hymenoptera. However, the specific solutions are remarkably different, thus revealing a striking case of convergence in one of the major evolutionary transitions in biological complexity.},
  author       = {Harrison, Mark and Jongepier, Evelien and Robertson, Hugh and Arning, Nicolas and Bitard Feildel, Tristan and Chao, Hsu and Childers, Christopher and Dinh, Huyen and Doddapaneni, Harshavardhan and Dugan, Shannon and Gowin, Johannes and Greiner, Carolin and Han, Yi and Hu, Haofu and Hughes, Daniel and Huylmans, Ann K and Kemena, Karsten and Kremer, Lukas and Lee, Sandra and López Ezquerra, Alberto and Mallet, Ludovic and Monroy Kuhn, Jose and Moser, Annabell and Murali, Shwetha and Muzny, Donna and Otani, Saria and Piulachs, Maria and Poelchau, Monica and Qu, Jiaxin and Schaub, Florentine and Wada Katsumata, Ayako and Worley, Kim and Xie, Qiaolin and Ylla, Guillem and Poulsen, Michael and Gibbs, Richard and Schal, Coby and Richards, Stephen and Belles, Xavier and Korb, Judith and Bornberg Bauer, Erich},
  journal      = {Nature Ecology and Evolution},
  number       = {3},
  pages        = {557--566},
  publisher    = {Springer Nature},
  title        = {{Hemimetabolous genomes reveal molecular basis of termite eusociality}},
  doi          = {10.1038/s41559-017-0459-1},
  volume       = {2},
  year         = {2018},
}

