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.

Genetic Coupling of Mate Recognition Systems in the Genomic Era.

The concept of "genetic coupling" in mate recognition systems arose in the 1960s as a potential mechanism to maintain coordination between signals and receivers during evolutionary divergence. At its most basic it proposed that the same genes might influence trait and preference, and therefore mutations could result in coordinated changes in both traits. Since then, the concept has expanded in scope and is often used to include linkage or genetic correlation between recognition system components. Here we review evidence for genetic coupling, concentrating on proposed examples of a common genetic basis for signals and preferences. Mapping studies have identified several examples of tight genetic linkage between genomic regions influencing signals and preferences, or assortative mating. Whether this extends as far as demonstrating pleiotropy remains a more open question. Some studies, notably of Drosophila, have identified genes in the sex determination pathway and in pheromonal communication where single loci can influence both signals and preferences. This may be based on isoform divergence, in which sex- and tissue-specific effects are facilitated by alternative spicing, or on regulatory divergence. Hence it is not clear that such examples provide compelling evidence of pleiotropy in the sense that "magic mutations" could maintain trait coordination. Rather, coevolution may be facilitated by regulatory divergence but require different mutations or coevolution across isoforms. Reconsidering the logic of genetic coupling, it may be that pleiotropy could actually be less effective than linkage if distinct but associated variants allow molecular coevolution to occur more readily than potentially "unbalanced" mutations in single genes. Genetic manipulation or studies of mutation order effects during divergence are challenging but perhaps the only way to disentangle the role of pleiotropy versus close linkage in coordinated trait divergence.

How Important Is Sexual Isolation to Speciation?

A central role for sexual isolation in the formation of new species and establishment of species boundaries has been noticed since Darwin and is frequently emphasized in the modern literature on speciation. However, an objective evaluation of when and how sexual isolation plays a role in speciation has been carried out in few taxa. We discuss three approaches for assessing the importance of sexual isolation relative to other reproductive barriers, including the relative evolutionary rate of sexual trait differentiation, the relative strength of sexual isolation in sympatry, and the role of sexual isolation in the long-term persistence of diverging forms. First, we evaluate evidence as to whether sexual isolation evolves faster than other reproductive barriers during the early stages of divergence. Second, we discuss available evidence as to whether sexual isolation is as strong or stronger than other barriers between closely related sympatric species. Finally, we consider the effect of sexual isolation on long-term species persistence, relative to other reproductive barriers. We highlight challenges to our knowledge of and opportunities to improve upon our understanding of sexual isolation from different phases of the speciation process.

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.

Selection on the Fly: Short-Term Adaptation to an Altered Sexual Selection Regime in Drosophila pseudoobscura.

Experimental evolution studies are powerful approaches to examine the evolutionary history of lab populations. Such studies have shed light on how selection changes phenotypes and genotypes. Most of these studies have not examined the time course of adaptation under sexual selection manipulation, by resequencing the populations' genomes at multiple time points. Here, we analyze allele frequency trajectories in Drosophila pseudoobscura where we altered their sexual selection regime for 200 generations and sequenced pooled populations at 5 time points. The intensity of sexual selection was either relaxed in monogamous populations (M) or elevated in polyandrous lines (E). We present a comprehensive study of how selection alters population genetics parameters at the chromosome and gene level. We investigate differences in the effective population size-Ne-between the treatments, and perform a genome-wide scan to identify signatures of selection from the time-series data. We found genomic signatures of adaptation to both regimes in D. pseudoobscura. There are more significant variants in E lines as expected from stronger sexual selection. However, we found that the response on the X chromosome was substantial in both treatments, more pronounced in E and restricted to the more recently sex-linked chromosome arm XR in M. In the first generations of experimental evolution, we estimate Ne to be lower on the X in E lines, which might indicate a swift adaptive response at the onset of selection. Additionally, the third chromosome was affected by elevated polyandry whereby its distal end harbors a region showing a strong signal of adaptive evolution especially in E lines.

Divergence and introgression among the virilis group of Drosophila.

Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and more recent gene flow between closely related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.

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.

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.

Divergence and correlated evolution of male wing spot and courtship display between Drosophila nepalensis and D. trilutea.

Male-specific wing spots are usually associated with wing displays in the courtship behavior of Drosophila and may play important roles in sexual selection. Two closely related species, D. nepalensis and D. trilutea, differ in wing spots and scissoring behavior. Here, we compare male morphological characters, pigmentation intensity of male wing spots, wing-scissoring behavior, courtship songs, and reproductive isolation between 2 species. F1 fertile females and sterile males result from the cross between females of D. nepalensis and males of D. trilutea. The pigmentation of wing spots is significantly weaker in D. trilutea than in D. nepalensis and the F1 hybrid. Males scissor both wings in front of the female during courtship, with a posture spreading wings more widely, and at a faster frequency in D. nepalensis than in D. trilutea and the F1s. Males of D. trilutea vibrate wings to produce 2 types (A and B) of pulse songs, whereas D. nepalensis and the F1s sing only type B songs. The incidence of wing vibration and scissoring during courtship suggests that wing vibration is essential but scissoring is a facultative courtship element for successful mating in both species. The association between the darker wing spots with more elaborate scissoring might be the consequence of correlated evolution of these traits in D. nepalensis; however, D. trilutea retains wing scissoring during courtship despite having weaker pigmentation of wing spots. The genetic architecture of 2 traits differs in the F1s, consistent with maternal or sex-linked effects for spots but nonadditive effects for scissoring.

DrosoPhyla: Resources for Drosophilid Phylogeny and Systematics.

The vinegar fly Drosophila melanogaster is a pivotal model for invertebrate development, genetics, physiology, neuroscience, and disease. The whole family Drosophilidae, which contains over 4,400 species, offers a plethora of cases for comparative and evolutionary studies. Despite a long history of phylogenetic inference, many relationships remain unresolved among the genera, subgenera, and species groups in the Drosophilidae. To clarify these relationships, we first developed a set of new genomic markers and assembled a multilocus data set of 17 genes from 704 species of Drosophilidae. We then inferred a species tree with highly supported groups for this family. Additionally, we were able to determine the phylogenetic position of some previously unplaced species. These results establish a new framework for investigating the evolution of traits in fruit flies, as well as valuable resources for systematics.

The discovery, distribution, and diversity of DNA viruses associated with Drosophila melanogaster in Europe.

Drosophila melanogaster is an important model for antiviral immunity in arthropods, but very few DNA viruses have been described from the family Drosophilidae. This deficiency limits our opportunity to use natural host-pathogen combinations in experimental studies, and may bias our understanding of the Drosophila virome. Here, we report fourteen DNA viruses detected in a metagenomic analysis of 6668 pool-sequenced Drosophila, sampled from forty-seven European locations between 2014 and 2016. These include three new nudiviruses, a new and divergent entomopoxvirus, a virus related to Leptopilina boulardi filamentous virus, and a virus related to Musca domestica salivary gland hypertrophy virus. We also find an endogenous genomic copy of galbut virus, a double-stranded RNA partitivirus, segregating at very low frequency. Remarkably, we find that Drosophila Vesanto virus, a small DNA virus previously described as a bidnavirus, may be composed of up to twelve segments and thus represent a new lineage of segmented DNA viruses. Two of the DNA viruses, Drosophila Kallithea nudivirus and Drosophila Vesanto virus are relatively common, found in 2 per cent or more of wild flies. The others are rare, with many likely to be represented by a single infected fly. We find that virus prevalence in Europe reflects the prevalence seen in publicly available datasets, with Drosophila Kallithea nudivirus and Drosophila Vesanto virus the only ones commonly detectable in public data from wild-caught flies and large population cages, and the other viruses being rare or absent. These analyses suggest that DNA viruses are at lower prevalence than RNA viruses in D.melanogaster, and may be less likely to persist in laboratory cultures. Our findings go some way to redressing an earlier bias toward RNA virus studies in Drosophila, and lay the foundation needed to harness the power of Drosophila as a model system for the study of DNA viruses.

Experimental evolution supports signatures of sexual selection in genomic divergence.

Comparative genomics has contributed to the growing evidence that sexual selection is an important component of evolutionary divergence and speciation. Divergence by sexual selection is implicated in faster rates of divergence of the X chromosome and of genes thought to underlie sexually selected traits, including genes that are sex biased in expression. However, accurately inferring the relative importance of complex and interacting forms of natural selection, demography, and neutral processes that occurred in the evolutionary past is challenging. Experimental evolution provides an opportunity to apply controlled treatments for multiple generations and examine the consequent genomic divergence. Here, we altered sexual selection intensity, elevating sexual selection in polyandrous lines and eliminating it in monogamous lines, and examined patterns of allele frequency divergence in the genome of Drosophila pseudoobscura after more than 160 generations of experimental evolution. Divergence is not uniform across the genome but concentrated in "islands," many of which contain candidate genes implicated in mating behaviors and other sexually selected phenotypes. These are more often seen on the X chromosome, which also shows greater divergence in F ST than neutral expectations. There are characteristic signatures of selection seen in these regions, with lower diversity on the X chromosome than the autosomes, and differences in diversity on the autosomes between selection regimes. Reduced Tajima's D within some of the divergent regions may imply that selective sweeps have occurred, despite considerable recombination. These changes are associated with both differential gene expression between the lines and sex-biased gene expression within the lines. Our results are very similar to those thought to implicate sexual selection in divergence between species and natural populations, and hence provide experimental support for the likely role of sexual selection in driving such types of genetic divergence, but also illustrate how variable outcomes can be for different genomic regions.

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.

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.

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.

Antibiotics for Uncomplicated Diverticulitis.

Learn more about the benefits and harms of antibiotics for uncomplicated diverticulitis.

Within-population sperm competition intensity does not predict asymmetry in conpopulation sperm precedence.

Postcopulatory sexual selection can generate evolutionary arms races between the sexes resulting in the rapid coevolution of reproductive phenotypes. As traits affecting fertilization success diverge between populations, postmating prezygotic (PMPZ) barriers to gene flow may evolve. Conspecific sperm precedence is a form of PMPZ isolation thought to evolve early during speciation yet has mostly been studied between species. Here, we show conpopulation sperm precedence (CpSP) between Drosophila montana populations. Using Pool-seq genomic data we estimate divergence times and ask whether PMPZ isolation evolved in the face of gene flow. We find models incorporating gene flow fit the data best indicating populations experienced considerable gene flow during divergence. We find CpSP is asymmetric and mirrors asymmetry in non-competitive PMPZ isolation, suggesting these phenomena have a shared mechanism. However, we show asymmetry is unrelated to the strength of postcopulatory sexual selection acting within populations. We tested whether overlapping foreign and coevolved ejaculates within the female reproductive tract altered fertilization success but found no effect. Our results show that neither time since divergence nor sperm competitiveness predicts the strength of PMPZ isolation. We suggest that instead cryptic female choice or mutation-order divergence may drive divergence of postcopulatory phenotypes resulting in PMPZ isolation. This article is part of the theme issue 'Fifty years of sperm competition'.

Dinner and a show: The role of male copulatory courtship song and female blood-feeding in the reproductive success of Lutzomyia longipalpis from Lapinha, Brazil.

Lutzomyia longipalpis is the main vector of visceral Leishmaniasis in the Americas and is composed of a species complex. Males of this sand-fly produce acoustic signals during copulation and different patterns are observed among Brazilian populations. Such acoustic signals are commonly involved in species recognition. However, since the song is only produced during copulation it is not clear how it affects mating success or contributes to sexual isolation. Another aspect that may affect reproductive success is the presence of food. Since hematophagy is such an important aspect of L. longipalpis biology, we wanted to test if blood-feeding can influence the reproductive behaviour of this insect. We performed crossing experiments removing males' wings (silencing them) and playing back either the homo-specific or the hetero-specific song to either unfed or blood-fed females. Our results showed that both songs and blood-feeding affect insemination success, but not the frequency of copulation. In trials where females were not blood-fed song clearly affected insemination; males with wings, and males with homo-specific song playback had a higher insemination success than wingless males (no song) and trials with hetero-specific song. Blood-feeding females prior to the trials increased insemination in all groups including the control group which suggests that mating happens simultaneously with, or immediately after, the blood meal. Blood-fed females also seemed to discriminate less against the wrong song or the lack of song (wingless) one day after feeding, however trials with the correct song still had higher insemination rates. Altogether, our results show that both the male copulatory courtship songs and female blood-feeding are important for reproductive success and as such are important components of the sexual behaviour of L. longipalpis.