Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes.

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology.

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Evidence for gene flow and trait reversal during radiation of Mexican Goodeid fish.

Understanding the phylogeographic history of a group and identifying the factors contributing to speciation is an important challenge in evolutionary biology. The Goodeinae are a group of live-bearing fishes endemic to Mexico. Here, we develop genomic resources for species within the Goodeinae and use phylogenomic approaches to characterise their evolutionary history. We sequenced, assembled and annotated the genomes of four Goodeinae species, including Ataeniobius toweri, the only matrotrophic live-bearing fish without a trophotaenia in the group. We estimated timings of species divergence and examined the extent and timing of introgression between the species to assess if this may have occurred during an early radiation, or in more recent episodes of secondary contact. We used branch-site models to detect genome-wide positive selection across Goodeinae, and we specifically asked whether this differs in A. toweri, where loss of placental viviparity has recently occurred. We found evidence of gene flow between geographically isolated species, suggesting vicariant speciation was supplemented by limited post-speciation gene flow, and gene flow may explain previous uncertainties about Goodeid phylogeny. Genes under positive selection in the group are likely to be associated with the switch to live-bearing. Overall, our studies suggest that both volcanism-driven vicariance and changes in reproductive mode influenced radiation in the Goodeinae.

Bringing bioinformatics to schools with the 4273pi project.

Over the last few decades, the nature of life sciences research has changed enormously, generating a need for a workforce with a variety of computational skills such as those required to store, manage, and analyse the large biological datasets produced by next-generation sequencing. Those with such expertise are increasingly in demand for employment in both research and industry. Despite this, bioinformatics education has failed to keep pace with advances in research. At secondary school level, computing is often taught in isolation from other sciences, and its importance in biological research is not fully realised, leaving pupils unprepared for the computational component of Higher Education and, subsequently, research in the life sciences. The 4273pi Bioinformatics at School project (https://4273pi.org) aims to address this issue by designing and delivering curriculum-linked, hands-on bioinformatics workshops for secondary school biology pupils, with an emphasis on equitable access. So far, we have reached over 180 schools across Scotland through visits or teacher events, and our open education resources are used internationally. Here, we describe our project, our aims and motivations, and the practical lessons we have learned from implementing a successful bioinformatics education project over the last 5 years.

Experimental sexual selection reveals rapid evolutionary divergence in sex-specific transcriptomes and their interactions following mating.

Post copulatory interactions between the sexes in internally fertilizing species elicits both sexual conflict and sexual selection. Macroevolutionary and comparative studies have linked these processes to rapid transcriptomic evolution in sex-specific tissues and substantial transcriptomic post mating responses in females, patterns of which are altered when mating between reproductively isolated species. Here, we tested multiple predictions arising from sexual selection and conflict theory about the evolution of sex-specific and tissue-specific gene expression and the post mating response at the microevolutionary level. Following over 150 generations of experimental evolution under either reduced (enforced monogamy) or elevated (polyandry) sexual selection in Drosophila pseudoobscura, we found a substantial effect of sexual selection treatment on transcriptomic divergence in virgin male and female reproductive tissues (testes, male accessory glands, the female reproductive tract and ovaries). Sexual selection treatment also had a dominant effect on the post mating response, particularly in the female reproductive tract - the main arena for sexual conflict - compared to ovaries. This effect was asymmetric with monandry females typically showing more post mating responses than polyandry females, with enriched gene functions varying across treatments. The evolutionary history of the male partner had a larger effect on the post mating response of monandry females, but females from both sexual selection treatments showed unique patterns of gene expression and gene function when mating with males from the alternate treatment. Our microevolutionary results mostly confirm comparative macroevolutionary predictions on the role of sexual selection on transcriptomic divergence and altered gene regulation arising from divergent coevolutionary trajectories between sexual selection treatments.

Purifying Selection in Corvids Is Less Efficient on Islands.

Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size Ne. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species.

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.

A genome-wide investigation of adaptive signatures in protein-coding genes related to tool behaviour in New Caledonian and Hawaiian crows.

Very few animals habitually manufacture and use tools. It has been suggested that advanced tool behaviour co-evolves with a suite of behavioural, morphological and life history traits. In fact, there are indications for such an adaptive complex in tool-using crows (genus Corvus species). Here, we sequenced the genomes of two habitually tool-using and ten non-tool-using crow species to search for genomic signatures associated with a tool-using lifestyle. Using comparative genomic and population genetic approaches, we screened for signals of selection in protein-coding genes in the tool-using New Caledonian and Hawaiian crows. While we detected signals of recent selection in New Caledonian crows near genes associated with bill morphology, our data indicate that genetic changes in these two lineages are surprisingly subtle, with little evidence at present for convergence. We explore the biological explanations for these findings, such as the relative roles of gene regulation and protein-coding changes, as well as the possibility that statistical power to detect selection in recently diverged lineages may have been insufficient. Our study contributes to a growing body of literature aiming to decipher the genetic basis of recently evolved complex behaviour.

Sex-specific responses to cold in a very cold-tolerant, northern Drosophila species.

Organisms can plastically alter resource allocation in response to changing environmental factors. For example, in harsh conditions, organisms are expected to shift investment from reproduction toward survival; however, the factors and mechanisms that govern the magnitude of such shifts are relatively poorly studied. Here we compared the impact of cold on males and females of the highly cold-tolerant species Drosophila montana at the phenotypic and transcriptomic levels. Although both sexes showed similar changes in cold tolerance and gene expression in response to cold treatment, indicating that the majority of changes are concordant between the sexes, we identified a clear reduction in sexually dimorphic gene expression, suggesting that preparing for the colder season involves reducing investment in sex-specific traits. This reduction was larger in males than females, as expected if male sexual traits are more condition-dependent than female traits, as predicted by theory. Gene expression changes were primarily associated with shifts in metabolic profile, which likely play a role in increasing cold tolerance. Finally, we found that the expression of immune genes was reduced following cold treatment, suggesting that reduced investment in costly immune function may be important in helping flies survive colder periods.

Sexual selection and population divergence III: Interspecific and intraspecific variation in mating signals.

A major challenge for studying the role of sexual selection in divergence and speciation is understanding the relative influence of different sexually selected signals on those processes in both intra- and interspecific contexts. Different signals may be more or less susceptible to co-option for species identification depending on the balance of sexual and ecological selection acting upon them. To examine this, we tested three predictions to explain geographic variation in long- versus short-range sexual signals across a 3,500 + km transect of two related Australian field cricket species (Teleogryllus spp.): (a) selection for species recognition, (b) environmental adaptation and (c) stochastic divergence. We measured male calling song and male and female cuticular hydrocarbons (CHCs) in offspring derived from wild populations, reared under common garden conditions. Song clearly differentiated the species, and no hybrids were observed suggesting that hybridization is rare or absent. Spatial variation in song was not predicted by geography, genetics or climatic factors in either species. In contrast, CHC divergence was strongly associated with an environmental gradient supporting the idea that the climatic environment selects more directly upon these chemical signals. In light of recently advocated models of diversification via ecological selection on secondary sexual traits, the different environmental associations we found for song and CHCs suggest that the impact of ecological selection on population divergence, and how that influences speciation, might be different for acoustic versus chemical signals.

Behavioural mechanisms of sexual isolation involving multiple modalities and their inheritance.

Speciation research dissects the genetics and evolution of reproductive barriers between parental species. Hybrids are the "gatekeepers" of gene flow, so it is also important to understand the behavioural mechanisms and genetics of any potential isolation from their parental species. We tested the role of multiple behavioural barriers in reproductive isolation among closely related field crickets and their hybrids (Teleogryllus oceanicus and Teleogryllus commodus). These species hybridize in the laboratory, but the behaviour of hybrids is unusual and there is little evidence for gene flow in the wild. We found that heterospecific pairs exhibited reduced rates of courtship behaviour due to discrimination by both sexes, and that this behavioural isolation was symmetrical. However, hybrids were not sexually selected against and exhibited high rates of courtship behaviour even though hybrid females are sterile. Using reciprocal hybrid crosses, we characterized patterns of interspecific divergence and inheritance in key sexual traits that might underlie the mating patterns we found: calling song, courtship song and cuticular hydrocarbons (CHCs). Song traits exhibited both sex linkage and transgressive segregation, whereas CHCs exhibited only the latter. Calculations of the strength of isolation exerted by these sexual traits suggest that close-range signals are as important as long-distance signals in contributing to interspecific sexual isolation. The surprisingly weak mating barriers observed between hybrids and parental species highlight the need to examine reproductive isolating mechanisms and their genetic bases across different potential stages of introgressive hybridization.

Evolution and diversity of the courtship repertoire in the Drosophila montium species group (Diptera: Drosophilidae).

Changes in elements of courtship behaviour can influence sexual isolation between species. Large-scale analyses of changes, including loss and gain of courtship elements, across a relatively complete phylogenetic group are rare but needed to understand the significance of such changes, for example whether the gain and loss of courtship elements are essentially arbitrary or equally reversible. In most species of Drosophila, courtship, including singing, mainly occurs before mounting as premounting courtship. The Drosophila montium species group is unusual because loss of premounting courtship and gain of post-mounting one has been detected in this group. Here, we provide an extensive analysis on the courtship repertoire and songs of 42 species in this group. Synchronously captured video and audio recordings were analysed to describe courtship patterns and male courtship songs, and changes were analysed in a phylogenetic context. Ancestral state reconstruction suggests that a gain of post-mounting courtship singing at the ancestor of this species group has been accompanied by a concurrent decrease in the incidence of premounting courtship singing and has led to subsequent further decrease and eventually complete loss of premounting courtship song in several lineages. Alongside this evolutionary trend towards post-mounting courtship, sine song and a special type of "high pulse repetition song" have become more widely used for courtship during species diversification in the montium group. It is likely that the elaboration of post-mounting courtship behaviours is associated with changes in the relative importance of pre- and post-mounting components of mating systems, such as sperm competition or cryptic female choice.

Selection for reproduction under short photoperiods changes diapause-associated traits and induces widespread genomic divergence.

The incidence of reproductive diapause is a critical aspect of life history in overwintering insects from temperate regions. Much has been learned about the timing, physiology and genetics of diapause in a range of insects, but how the multiple changes involved in this and other photoperiodically regulated traits are inter-related is not well understood. We performed quasinatural selection on reproduction under short photoperiods in a northern fly species, Drosophila montana, to trace the effects of photoperiodic selection on traits regulated by the photoperiodic timer and/or by a circadian clock system. Selection changed several traits associated with reproductive diapause, including the critical day length for diapause (CDL), the frequency of diapausing females under photoperiods that deviate from daily 24 h cycles and cold tolerance, towards the phenotypes typical of lower latitudes. However, selection had no effect on the period of free-running locomotor activity rhythm regulated by the circadian clock in fly brain. At a genomic level, selection induced extensive divergence from the control line in 16 gene clusters involved in signal transduction, membrane properties, immunologlobulins and development. These changes resembled those detected between latitudinally divergent D. montana populations in the wild and involved SNP divergence associated with several genes linked with diapause induction. Overall, our study shows that photoperiodic selection for reproduction under short photoperiods affects diapause-associated traits without disrupting the central clock network generating circadian rhythms in fly locomotor activity.

Increased socially mediated plasticity in gene expression accompanies rapid adaptive evolution.

Recent theory predicts that increased phenotypic plasticity can facilitate adaptation as traits respond to selection. When genetic adaptation alters the social environment, socially mediated plasticity could cause co-evolutionary feedback dynamics that increase adaptive potential. We tested this by asking whether neural gene expression in a recently arisen, adaptive morph of the field cricket Teleogryllus oceanicus is more responsive to the social environment than the ancestral morph. Silent males (flatwings) rapidly spread in a Hawaiian population subject to acoustically orienting parasitoids, changing the population's acoustic environment. Experimental altering crickets' acoustic environments during rearing revealed broad, plastic changes in gene expression. However, flatwing genotypes showed increased socially mediated plasticity, whereas normal-wing genotypes exhibited negligible expression plasticity. Increased plasticity in flatwing crickets suggests a coevolutionary process coupling socially flexible gene expression with the abrupt spread of flatwing. Our results support predictions that phenotypic plasticity should rapidly evolve to be more pronounced during early phases of adaptation.

Sexual selection predicts species richness across the animal kingdom.

Our improving knowledge of the animal tree of life consistently demonstrates that some taxa diversify more rapidly than others, but what contributes to this variation remains poorly understood. An influential hypothesis proposes that selection arising from competition for mating partners plays a key role in promoting speciation. However, empirical evidence showing a link between proxies of this sexual selection and species richness is equivocal. Here, we collected standardized metrics of sexual selection for a broad range of animal taxa, and found that taxonomic families characterized by stronger sexual selection on males show relatively higher species richness. Thus, our data support the hypothesis that sexual selection elevates species richness. This could occur either by promoting speciation and/or by protecting species against extinction.

Opposing patterns of intraspecific and interspecific differentiation in sex chromosomes and autosomes.

Linking intraspecific and interspecific divergence is an important challenge in speciation research. X chromosomes are expected to evolve faster than autosomes and disproportionately contribute to reproductive barriers, and comparing genetic variation on X and autosomal markers within and between species can elucidate evolutionary processes that shape genome variation. We performed RADseq on a 16 population transect of two closely related Australian cricket species, Teleogryllus commodus and T. oceanicus, covering allopatry and sympatry. This classic study system for sexual selection provides a rare exception to Haldane's rule, as hybrid females are sterile. We found no evidence of recent introgression, despite the fact that the species coexist in overlapping habitats in the wild and interbreed in the laboratory. Putative X-linked loci showed greater differentiation between species compared with autosomal loci. However, population differentiation within species was unexpectedly lower on X-linked markers than autosomal markers, and relative X-to-autosomal genetic diversity was inflated above neutral expectations. Populations of both species showed genomic signatures of recent population expansions, but these were not strong enough to account for the inflated X/A diversity. Instead, most of the excess polymorphism on the X could better be explained by sex-biased processes that increase the relative effective population size of the X, such as interspecific variation in the strength of sexual selection among males. Taken together, the opposing patterns of diversity and differentiation at X versus autosomal loci implicate a greater role for sex-linked genes in maintaining species boundaries in this system.

Insect capa neuropeptides impact desiccation and cold tolerance.

The success of insects is linked to their impressive tolerance to environmental stress, but little is known about how such responses are mediated by the neuroendocrine system. Here we show that the capability (capa) neuropeptide gene is a desiccation- and cold stress-responsive gene in diverse dipteran species. Using targeted in vivo gene silencing, physiological manipulations, stress-tolerance assays, and rationally designed neuropeptide analogs, we demonstrate that the Drosophila melanogaster capa neuropeptide gene and its encoded peptides alter desiccation and cold tolerance. Knockdown of the capa gene increases desiccation tolerance but lengthens chill coma recovery time, and injection of capa peptide analogs can reverse both phenotypes. Immunohistochemical staining suggests that capa accumulates in the capa-expressing Va neurons during desiccation and nonlethal cold stress but is not released until recovery from each stress. Our results also suggest that regulation of cellular ion and water homeostasis mediated by capa peptide signaling in the insect Malpighian (renal) tubules is a key physiological mechanism during recovery from desiccation and cold stress. This work augments our understanding of how stress tolerance is mediated by neuroendocrine signaling and illustrates the use of rationally designed peptide analogs as agents for disrupting protective stress tolerance.

Genome-wide DNA methylation patterns in wild samples of two morphotypes of threespine stickleback (Gasterosteus aculeatus).

Epigenetic marks such as DNA methylation play important biological roles in gene expression regulation and cellular differentiation during development. To examine whether DNA methylation patterns are potentially associated with naturally occurring phenotypic differences, we examined genome-wide DNA methylation within Gasterosteus aculeatus, using reduced representation bisulfite sequencing. First, we identified highly methylated regions of the stickleback genome, finding such regions to be located predominantly within genes, and associated with genes functioning in metabolism and biosynthetic processes, cell adhesion, signaling pathways, and blood vessel development. Next, we identified putative differentially methylated regions (DMRs) of the genome between complete and low lateral plate morphs of G. aculeatus. We detected 77 DMRs that were mainly located in intergenic regions. Annotations of genes associated with these DMRs revealed potential functions in a number of known divergent adaptive phenotypes between G. aculeatus ecotypes, including cardiovascular development, growth, and neuromuscular development.

A test of genetic models for the evolutionary maintenance of same-sex sexual behaviour.

The evolutionary maintenance of same-sex sexual behaviour (SSB) has received increasing attention because it is perceived to be an evolutionary paradox. The genetic basis of SSB is almost wholly unknown in non-human animals, though this is key to understanding its persistence. Recent theoretical work has yielded broadly applicable predictions centred on two genetic models for SSB: overdominance and sexual antagonism. Using Drosophila melanogaster, we assayed natural genetic variation for male SSB and empirically tested predictions about the mode of inheritance and fitness consequences of alleles influencing its expression. We screened 50 inbred lines derived from a wild population for male-male courtship and copulation behaviour, and examined crosses between the lines for evidence of overdominance and antagonistic fecundity selection. Consistent variation among lines revealed heritable genetic variation for SSB, but the nature of the genetic variation was complex. Phenotypic and fitness variation was consistent with expectations under overdominance, although predictions of the sexual antagonism model were also supported. We found an unexpected and strong paternal effect on the expression of SSB, suggesting possible Y-linkage of the trait. Our results inform evolutionary genetic mechanisms that might maintain low but persistently observed levels of male SSB in D. melanogaster, but highlight a need for broader taxonomic representation in studies of its evolutionary causes.