@article{6710,
  abstract     = {Sexual dimorphism in morphology, physiology or life history traits is common in dioecious plants at reproductive maturity, but it is typically inconspicuous or absent in juveniles. Although plants of different sexes probably begin to diverge in gene expression both before their reproduction commences and before dimorphism becomes readily apparent, to our knowledge transcriptome-wide differential gene expression has yet to be demonstrated for any angiosperm species.},
  author       = {Cossard, Guillaume and Toups, Melissa A and Pannell, John },
  issn         = {1095-8290},
  journal      = {Annals of botany},
  number       = {7},
  pages        = {1119--1131},
  publisher    = {Oxford University Press},
  title        = {{Sexual dimorphism and rapid turnover in gene expression in pre-reproductive seedlings of a dioecious herb}},
  doi          = {10.1093/aob/mcy183},
  volume       = {123},
  year         = {2019},
}

@article{6755,
  abstract     = {Differentiated sex chromosomes are accompanied by a difference in gene dose between X/Z-specific and autosomal genes. At the transcriptomic level, these sex-linked genes can lead to expression imbalance, or gene dosage can be compensated by epigenetic mechanisms and results into expression level equalization. Schistosoma mansoni has been previously described as a ZW species (i.e., female heterogamety, in opposition to XY male heterogametic species) with a partial dosage compensation, but underlying mechanisms are still unexplored. Here, we combine transcriptomic (RNA-Seq) and epigenetic data (ChIP-Seq against H3K4me3, H3K27me3,andH4K20me1histonemarks) in free larval cercariae and intravertebrate parasitic stages. For the first time, we describe differences in dosage compensation status in ZW females, depending on the parasitic status: free cercariae display global dosage compensation, whereas intravertebrate stages show a partial dosage compensation. We also highlight regional differences of gene expression along the Z chromosome in cercariae, but not in the intravertebrate stages. Finally, we feature a consistent permissive chromatin landscape of the Z chromosome in both sexes and stages. We argue that dosage compensation in schistosomes is characterized by chromatin remodeling mechanisms in the Z-specific region.},
  author       = {Picard, Marion A L and Vicoso, Beatriz and Roquis, David and Bulla, Ingo and Augusto, Ronaldo C. and Arancibia, Nathalie and Grunau, Christoph and Boissier, Jérôme and Cosseau, Céline},
  issn         = {1759-6653},
  journal      = {Genome biology and evolution},
  number       = {7},
  pages        = {1909--1922},
  publisher    = {Oxford Academic Press},
  title        = {{Dosage compensation throughout the Schistosoma mansoni lifecycle: Specific chromatin landscape of the Z chromosome}},
  doi          = {10.1093/gbe/evz133},
  volume       = {11},
  year         = {2019},
}

@article{6821,
  abstract     = {To determine the visual sensitivities of an organism of interest, quantitative reverse transcription–polymerase chain reaction (qRT–PCR) is often used to quantify expression of the light‐sensitive opsins in the retina. While qRT–PCR is an affordable, high‐throughput method for measuring expression, it comes with inherent normalization issues that affect the interpretation of results, especially as opsin expression can vary greatly based on developmental stage, light environment or diurnal cycles. We tested for diurnal cycles of opsin expression over a period of 24 hr at 1‐hr increments and examined how normalization affects a data set with fluctuating expression levels using qRT–PCR and transcriptome data from the retinae of the cichlid Pelmatolapia mariae. We compared five methods of normalizing opsin expression relative to (a) the average of three stably expressed housekeeping genes (Ube2z, EF1‐α and β‐actin), (b) total RNA concentration, (c) GNAT2, (the cone‐specific subunit of transducin), (d) total opsin expression and (e) only opsins expressed in the same cone type. Normalizing by proportion of cone type produced the least variation and would be best for removing time‐of‐day variation. In contrast, normalizing by housekeeping genes produced the highest daily variation in expression and demonstrated that the peak of cone opsin expression was in the late afternoon. A weighted correlation network analysis showed that the expression of different cone opsins follows a very similar daily cycle. With the knowledge of how these normalization methods affect opsin expression data, we make recommendations for designing sampling approaches and quantification methods based upon the scientific question being examined.},
  author       = {Yourick, Miranda R. and Sandkam, Benjamin A. and Gammerdinger, William J and Escobar-Camacho, Daniel and Nandamuri, Sri Pratima and Clark, Frances E. and Joyce, Brendan and Conte, Matthew A. and Kocher, Thomas D. and Carleton, Karen L.},
  issn         = {1755-0998},
  journal      = {Molecular Ecology Resources},
  number       = {6},
  pages        = {1447--1460},
  publisher    = {Wiley},
  title        = {{Diurnal variation in opsin expression and common housekeeping genes necessitates comprehensive normalization methods for quantitative real-time PCR analyses}},
  doi          = {10.1111/1755-0998.13062},
  volume       = {19},
  year         = {2019},
}

@article{6831,
  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 information is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life‐cycle dynamics.
* Here, we investigated 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.},
  issn         = {1469-8137},
  journal      = {New Phytologist},
  number       = {3},
  pages        = {1108--1120},
  publisher    = {Wiley},
  title        = {{Variation in sexual dimorphism in a wind-pollinated plant: The influence of geographical context and life-cycle dynamics}},
  doi          = {10.1111/nph.16050},
  volume       = {224},
  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{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{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{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},
}

@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},
}

@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},
}