Supergenes are not necessary to explain the maintenance of complex alternative phenotypes

Evolutionary biology aims to explain the diversity seen in nature. Evolutionary theory provides frameworks to understand how simple polymorphisms or continuous variation are maintained, but phenotypes inherited as discrete suites of quantitative traits are difficult to fit into this framework. Supergenes have been proposed as a solution to this problem-if causal genes are co-located, they can be inherited as if a single gene, thus bridging the gap between simple polymorphisms and continuous traits. We develop models to ask: how are critical supergenes for maintaining phenotypic diversity? In our simplest model, without explicit genetic architectures, three alternative reproductive morphs are maintained in many of the parameter combinations we evaluated. For these same parameter values, models with demographic stochasticity, recombination and mutation (but without explicit genetic architecture) maintained only two of these three morphs, with stochasticity determining which morphs persisted. With explicit genetic architectures, regardless of whether causal loci were co-located in a supergene or distributed randomly, this stochasticity in which morphs are maintained was reduced. Even when phenotypic variation was lost, genetic diversity was maintained. Altogether, categorical traits with polygenic bases exhibited similar evolutionary dynamics to those determined by supergenes. Our work suggests that supergenes are not the only answer to the puzzle of how discrete polygenic phenotypic variation is maintained.

Transcriptomics of differences in thermal plasticity associated with selection for an exaggerated male sexual trait

The information about the magnitude of differences in thermal plasticity both between and within populations, as well as identification of the underlying molecular mechanisms are key to understanding the evolution of thermal plasticity. In particular, genes underlying variation in the physiological response to temperature can provide raw material for selection acting on plastic traits. Using RNAseq, we investigate the transcriptional response to temperature in males and females from bulb mite populations selected for the increased frequency of one of two discrete male morphs (fighter- and scrambler-selected populations) that differ in relative fitness depending on temperature. We show that different mechanisms underlie the divergence in thermal response between fighter- and scrambler-selected populations at decreased vs. increased temperature. Temperature decrease to 18 °C was associated with higher transcriptomic plasticity of males with more elaborate armaments, as indicated by a significant selection-by-temperature interaction effect on the expression of 40 genes, 38 of which were upregulated in fighter-selected populations in response to temperature decrease. In response to 28 °C, no selection-by-temperature interaction in gene expression was detected. Hence, differences in phenotypic response to temperature increase likely depended on genes associated with their distinct morph-specific thermal tolerance. Selection of males also drove gene expression patterns in females. These patterns could be associated with temperature-dependent fitness differences between females from fighter- vs. scrambler-selected populations reported in previous studies. Our study shows that selection for divergent male sexually selected morphologies and behaviors has a potential to drive divergence in metabolic pathways underlying plastic response to temperature in both sexes.

Sex-biased gene expression at single-cell resolution: cause and consequence of sexual dimorphism

Gene expression differences between males and females are thought to be key for the evolution of sexual dimorphism, and sex-biased genes are often used to study the molecular footprint of sex-specific selection. However, gene expression is often measured from complex aggregations of diverse cell types, making it difficult to distinguish between sex differences in expression that are due to regulatory rewiring within similar cell types and those that are simply a consequence of developmental differences in cell-type abundance. To determine the role of regulatory versus developmental differences underlying sex-biased gene expression, we use single-cell transcriptomic data from multiple somatic and reproductive tissues of male and female guppies, a species that exhibits extensive phenotypic sexual dimorphism. Our analysis of gene expression at single-cell resolution demonstrates that nonisometric scaling between the cell populations within each tissue and heterogeneity in cell-type abundance between the sexes can influence inferred patterns of sex-biased gene expression by increasing both the false-positive and false-negative rates. Moreover, we show that, at the bulk level, the subset of sex-biased genes that are the product of sex differences in cell-type abundance can significantly confound patterns of coding-sequence evolution. Taken together, our results offer a unique insight into the effects of allometry and cellular heterogeneity on perceived patterns of sex-biased gene expression and highlight the power of single-cell RNA-sequencing in distinguishing between sex-biased genes that are the result of regulatory change and those that stem from sex differences in cell-type abundance, and hence are a consequence rather than a cause of sexual dimorphism.

Different waves of postglacial recolonisation and genomic structure of bank vole populations in NE Poland

Previous studies indicated that in some species phylogeographic patterns obtained in the analysis of nuclear and mitochondrial DNA (mtDNA) markers can be different. Such mitonuclear discordance can have important evolutionary and ecological consequences. In the present study, we aimed to check whether there was any discordance between mtDNA and nuclear DNA in the bank vole population in the contact zone of its two mtDNA lineages. We analysed the population genetic structure of bank voles using genome-wide genetic data (SNPs) and diversity of sequenced heart transcriptomes obtained from selected individuals from three populations inhabiting areas outside the contact zone. The SNP genetic structure of the populations confirmed the presence of at least two genetic clusters, and such division was concordant with the patterns obtained in the analysis of other genetic markers and functional genes. However, genome-wide SNP analyses revealed the more detailed structure of the studied population, consistent with more than two bank vole recolonisation waves, as recognised previously in the study area. We did not find any significant differences between individuals representing two separate mtDNA lineages of the species in functional genes coding for protein-forming complexes, which are involved in the process of cell respiration in mitochondria. We concluded that the contemporary genetic structure of the populations and the width of the contact zone were shaped by climatic and environmental factors rather than by genetic barriers. The studied populations were likely isolated in separate Last Glacial Maximum refugia for insufficient amount of time to develop significant genetic differentiation.

Immunoecology of Species with Alternative Reproductive Tactics and Strategies

Alternative reproductive tactics and strategies (ARTS) refer to polymorphic reproductive behaviours in which in addition to the usual two sexes, there are one or more alternative morphs, usually male, that have evolved the ability to circumvent direct intra-sexual competition. Each morph has its own morphological, ecological, developmental, behavioural, life-history, and physiological profile that shifts the balance between reproduction and self-maintenance, one aspect being immunity. Immunoecological work on species with ARTS, which is the topic of this review, is particularly interesting because the alternative morphs make it possible to separate the effects of sex per se from other factors that in other species are inextricably linked with sex. We first summarize the evolution, development, and maintenance of ARTS. We then review immunoecological hypotheses relevant to species with ARTS, dividing them into physiological, life-history, and ecological hypotheses. In context of these hypotheses, we critically review in detail all immunoecological studies we could find on species with ARTS. Several interesting patterns emerge. Oddly, there is a paucity of studies on insects, despite the many benefits that arise from working with insects: larger sample sizes, simple immune systems, and countless forms of alternative reproductive strategies and tactics. Of all the hypotheses considered, the immunocompetence handicap hypothesis has generated the greatest amount of work, but not necessarily the greatest level of understanding. Unfortunately, it is often used as a general guiding principle rather than a source of explicitly articulated predictions. Other hypotheses are usually considered a posteriori, but perhaps they should take centre stage. Whereas blanket concepts such as “immunocompetence” and “androgens” might be useful to develop a rationale, predictions need to be far more explicitly articulated. Integration so far has been a one-way street, with ecologists delving deeper into physiology, sometimes at the cost of ignoring their organisms’ evolutionary history and ecology. One possible useful framework is to divide ecological and evolutionary factors affecting immunity into those that stimulate the immune system, and those that depress it. Finally, the contributions of genomics to ecology are being increasingly recognized and sometimes applied to species with ARTS, but we must ensure that evolutionary and ecological hypotheses drive the effort, as there is no grandeur in the strict reductionist view of life.

Intralocus conflicts associated with a supergene

Chromosomal inversions frequently underlie major phenotypic variation maintained by divergent selection within and between sexes. Here we examine whether and how intralocus conflicts contribute to balancing selection stabilizing an autosomal inversion polymorphism in the ruff Calidris pugnax. In this lekking shorebird, three male mating morphs (Independents, Satellites and Faeders) are controlled by an inversion-based supergene. We show that in a captive population, Faeder females, who are smaller and whose inversion haplotype has not undergone recombination, have lower average reproductive success in terms of laying rate, egg size, and offspring survival than Independent females, who lack the inversion. Satellite females, who carry a recombined inversion haplotype and have intermediate body size, more closely resemble Independent than Faeder females in reproductive performance. We inferred that the lower reproductive output of Faeder females is most likely balanced by higher than average reproductive success of individual Faeder males. These findings suggest that intralocus conflicts may play a major role in the evolution and maintenance of supergene variants.

Comparative Transcriptome Analysis Reveals the Mechanism Related to Fluazinam Stress of Panonychus citri (Acarina: Tetranychidae)

Simple Summary

The citrus red mite, Panonychus citri, is an important pest that causes serious citrus production losses in China. The insecticide fluazinam has a good control effect on the pest mites; however, its mechanism of action on mites remains unclear. In this study, we analyzed the transcriptomic sequencing and differential expression genes in P. citri treated with fluazinam, and identified some of the genes potential involved in detoxification metabolism related with the fluazinam exposure. Evaluating the efficacy of fluazinam, and analyzing the transcriptome data of P. citri under fluazinam stress, potentially provide a new agent for prevention and control of P. citri, and also preliminary research results for exploring the mechanism of action of fluazinam on P. citri. Given the up-regulated expression levels of genes for Mn-superoxide dismutase and catalase, we speculate that they play an important role in fluazinam-stress action on P. citri.

Abstract

The use of a large number of chemical acaricides to control these pest mites has led to an increasing problem of pesticide resistance, which has always been the difficulty in integrated pest management (IPM). Fluazinam has a good control effect on Panonychus citri, the serious pest on citrus; however, we only know the mechanism of action of fluazinam as a fungicide and its mechanism of action on mites remains unclear. Through analysis using Illumina high-throughput transcriptomic sequencing and differential expression genes in P. citri treated with fluazinam, 59 cytochrome P450 genes, 23 glutathione s-transferase genes, five carboxylate esterase genes, 11 superoxide dismutase genes and 15 catalase genes were identified. The Gene Ontology enrichment and the enrichment of KEGG results showed that the treatment were enrichment for redox enzyme pathways. Evaluating the efficacy of fluazinam, and analyzing the transcriptome data of P. citri under fluazinam stress, potentially provide a new agent for prevention and control of P. citri, and also preliminary research results for exploring the mechanism of action of fluazinam on P. citri. Given the up-regulated expression levels of genes for Mn-superoxide dismutase and catalase, we speculate that they play an important role in fluazinam-stress action on P. citri.

Release from intralocus sexual conflict? Evolved loss of a male sexual trait demasculinizes female gene expression

The loss of sexual ornaments is observed across taxa, and pleiotropic effects of such losses provide an opportunity to gain insight into underlying dynamics of sex-biased gene expression and intralocus sexual conflict (IASC). We investigated this in a Hawaiian field cricket, Teleogryllus oceanicus, in which an X-linked genotype ( flatwing) feminizes males' wings and eliminates their ability to produce sexually selected songs. We profiled adult gene expression across somatic and reproductive tissues of both sexes. Despite the feminizing effect of flatwing on male wings, we found no evidence of feminized gene expression in males. Instead, female transcriptomes were more strongly affected by flatwing than males', and exhibited demasculinized gene expression. These findings are consistent with a relaxation of IASC constraining female gene expression through loss of a male sexual ornament. In a follow-up experiment, we found reduced testes mass in flatwing males, whereas female carriers showed no reduction in egg production. By contrast, female carriers exhibited greater measures of body condition. Our results suggest sex-limited phenotypic expression offers only partial resolution to IASC, owing to pleiotropic effects of the loci involved. Benefits conferred by release from intralocus conflict could help explain widespread loss of sexual ornaments across taxa.

Female Mimicry by Sneaker Males Has a Transcriptomic Signature in Both the Brain and the Gonad in a Sex-Changing Fish

Phenotypic plasticity represents an elegant adaptive response of individuals to a change in their environment. Bluehead wrasses (Thalassoma bifasciatum) exhibit astonishing sexual plasticity, including female-to-male sex change and discrete male morphs that differ strikingly in behavior, morphology, and gonadal investment. Using RNA-seq transcriptome profiling, we examined the genes and physiological pathways underlying flexible behavioral and gonadal differences among female, dominant (bourgeois) male, and female-mimic (sneaker) male blueheads. For the first time in any organism, we find that female mimicry by sneaker males has a transcriptional signature in both the brain and the gonad. Sneaker males shared striking similarity in neural gene expression with females, supporting the idea that males with alternative reproductive phenotypes have "female-like brains." Sneaker males also overexpressed neuroplasticity genes, suggesting that their opportunistic reproductive strategy requires a heightened capacity for neuroplasticity. Bourgeois males overexpressed genes associated with socio-sexual behaviors (e.g., isotocin), but also neuroprotective genes and biomarkers of oxidative stress and aging, indicating a hitherto unexplored cost to these males of attaining the reproductively privileged position at the top of the social hierarchy. Our novel comparison of testicular transcriptomes in a fish with male sexual polymorphism associates greater gonadal investment by sneaker males with overexpression of genes involved in cell proliferation and sperm quality control. We propose that morphological female-mimicry by sneaker male teleosts entails pervasive downregulation of androgenesis genes, consistent with low androgen production in males lacking well-developed secondary sexual characters.

The role of genetic diversity in the evolution and maintenance of environmentally-cued, male alternative reproductive tactics

BACKGROUND

Alternative reproductive tactics (ARTs) are taxonomically pervasive strategies adopted by individuals to maximize reproductive success within populations. Even for conditionally-dependent traits, consensus postulates most ARTs involve both genetic and environmental interactions (GEIs), but to date, quantifying genetic variation underlying the threshold disposing an individual to switch phenotypes in response to an environmental cue has been a difficult undertaking. Our study aims to investigate the origins and maintenance of ARTs within environmentally disparate populations of the microscopic bulb mite, Rhizoglyphus robini, that express 'fighter' and 'scrambler' male morphs mediated by a complex combination of environmental and genetic factors.

RESULTS

Using never-before-published individual genetic profiling, we found all individuals across populations are highly inbred with the exception of scrambler males in stressed environments. In fact within the poor environment, scrambler males and females showed no significant difference in genetic differentiation (Fst) compared to all other comparisons, and although fighters were highly divergent from the rest of the population in both poor or rich environments (e.g., Fst, STRUCTURE), fighters demonstrated approximately three times less genetic divergence from the population in poor environments. AMOVA analyses further corroborated significant genetic differentiation across subpopulations, between morphs and sexes, and among subpopulations within each environment.

CONCLUSION

Our study provides new insights into the origin of ARTs in the bulb mite, highlighting the importance of GEIs: genetic correlations, epistatic interactions, and sex-specific inbreeding depression across environmental stressors. Asymmetric reproductive output, coupled with the purging of highly inbred individuals during environmental oscillations, also facilitates genetic variation within populations, despite evidence for strong directional selection. This cryptic genetic variation also conceivably facilitates stable population persistence even in the face of spatially or temporally unstable environmental challenges. Ultimately, understanding the genetic context that maintains thresholds, even for conditionally-dependent ARTs, will enhance our understanding of within population variation and our ability to predict responses to selection.

Genomics of end-Pleistocene population replacement in a small mammal

Current species distributions at high latitudes are the product of expansion from glacial refugia into previously uninhabitable areas at the end of the last glaciation. The traditional view of postglacial colonization is that southern populations expanded their ranges into unoccupied northern territories. Recent findings on mitochondrial DNA (mtDNA) of British small mammals have challenged this simple colonization scenario by demonstrating a more complex genetic turnover in Britain during the Pleistocene-Holocene transition where one mtDNA clade of each species was replaced by another mtDNA clade of the same species. Here, we provide evidence from one of those small mammals, the bank vole (Clethrionomys glareolus), that the replacement was genome-wide. Using more than 10 000 autosomal SNPs we found that similar to mtDNA, bank vole genomes in Britain form two (north and south) clusters which admix. Therefore, the genome of the original postglacial colonists (the northern cluster) was probably replaced by another wave of migration from a different continental European population (the southern cluster), and we gained support for this by modelling with approximate Bayesian computation. This finding emphasizes the importance of analysis of genome-wide diversity within species under changing climate in creating opportunities for sophisticated testing of population history scenarios.

The ART of brain expression

When trying to find a mate, one might think about using a dating app. Imagine that someone else has installed the same app and tries to access the same potential mates that you have chosen, and that this someone uses false facts about himself/herself to increase the chance of dating someone before you or anyone else with honest status information does. Sounds familiar? It actually is, and in no way is such comportment restricted to human courtship behaviour. Alternative reproductive tactics (ARTs) are widespread in the animal kingdom. In this issue of Molecular Ecology, Cardoso, Gonçalves, Goesmann, Canário, and Oliveira () investigate plastic ARTs of the peacock blenny (Salaria pavo), in which males occur in three morphs: nestholders, sneakers and transitionals between the two former. They apply transcriptome sequencing to answer the question how brain gene expression contributes to sex role-specific behaviour and to intersex phenotypes.

Sexual dimorphism and heightened conditional expression in a sexually selected weapon in the Asian rhinoceros beetle

Among the most dramatic examples of sexual selection are the weapons used in battles between rival males over access to females. As with ornaments of female choice, the most "exaggerated" sexually selected weapons vary from male to male more widely than other body parts (hypervariability), and their growth tends to be more sensitive to nutritional state or physiological condition compared with growth of other body parts ("heightened" conditional expression). Here, we use RNAseq analysis to build on recent work exploring these mechanisms in the exaggerated weapons of beetles, by examining patterns of differential gene expression in exaggerated (head and thorax horns) and non-exaggerated (wings, genitalia) traits in the Asian rhinoceros beetle, Trypoxylus dichotomus. Our results suggest that sexually dimorphic expression of weaponry involves large-scale changes in gene expression, relative to other traits, while nutrition-driven changes in gene expression in these same weapons are less pronounced. However, although fewer genes overall were differentially expressed in high- vs. low-nutrition individuals, the number of differentially expressed genes varied predictably according to a trait's degree of condition dependence (head horn > thorax horn > wings > genitalia). Finally, we observed a high degree of similarity in direction of effects (vectors) for subsets of differentially expressed genes across both sexually dimorphic and nutritionally responsive growth. Our results are consistent with a common set of mechanisms governing sexual size dimorphism and condition dependence.

Costs of weaponry: Unarmed males sire more offspring than armed males in a male-dimorphic mite

Morphological structures used as weapons in male–male competition are not only costly to develop but are also probably costly to maintain during adulthood. Therefore, having weapons could reduce the energy available for other fitness‐enhancing actions, such as post‐copulatory investment. We tested the hypothesis that armed males make lower post‐copulatory investments than unarmed males, and that this difference will be most pronounced under food‐limited conditions. We performed two experiments using the male‐dimorphic bulb mite Rhizoglyphus robini, in which males are either armed “fighters” or unarmed “scramblers.” Firstly, we tested whether fighters and scramblers differed in their reproductive output after being starved or fed for 1 or 2 weeks. Secondly, we measured the reproductive output of scramblers and fighters (starved or fed) after one, two or three consecutive matings. Scramblers sired more offspring than fighters after 1 week, but scramblers and fighters only sired a few offspring after 2 weeks. Scramblers also sired more offspring than fighters at the first mating, and males rarely sired offspring after consecutive matings. Contrary to our hypothesis, the fecundity of starved and fed males did not differ. The higher reproductive output of scramblers suggests that, regardless of nutritional state, scramblers make larger post‐copulatory investments than fighters. Alternatively, (cryptic) female choice generally favours scramblers. Why the morphs differed in their reproductive output is unclear. Neither morph performed well relatively late in life or after multiple matings. It remains to be investigated to what extent the apparent scrambler advantage contributes to the maintenance and evolution of male morph expression.

Life-history consequences of bidirectional selection for male morph in a male-dimorphic bulb mite

Intralocus sexual conflict (IASC) arises when males and females have different trait optima. Some males pursue different alternative reproductive tactics (ARTs) with different trait optima, resulting in different strengths of IASC. Consequently, for instance daughter fitness is differentially affected by her sire's morph. We tested if-and which-other life-history traits correlatively change in bidirectional, artificial selection experiments for ARTs. We used the male-dimorphic bulb mite Rhizoglyphus robini, the males of which are high-fitness 'fighters' or low-fitness 'scramblers'. Twice in each of the five generations of selection, we assessed clutch composition (number of mites of the various life stages present) and size (total number of offspring). Furthermore, we tracked offspring from egg to adulthood in the first and final generation to detect differences between selection lines in the size and duration of stages, and in maturation time. We found that selection for male morph increased the frequency of that morph. Furthermore, compared to fighter lines, scrambler lines produced more females, which laid larger eggs (in the final generations), and maintained a higher egg-laying rate for longer. Otherwise, our results showed no consistent differences between the selection lines in clutch size and composition, life stage size or duration, or maturation time. Though we found few correlated life-history trait changes in response to selection on male morph, the differences in egg laying rate and egg size suggest that IASC between fighters is costlier to females than IASC with scramblers. We hypothesize that these differences in reproductive traits allow fighter-offspring to perform better in small, declining populations but scrambler-offspring to perform better in large, growing populations.

Masculinization of gene expression is associated with male quality in Drosophila melanogaster

The signature of sexual selection has been revealed through the study of differences in patterns of genome-wide gene expression, both between the sexes and between alternative reproductive morphs within a single sex. What remains unclear, however, is whether differences in gene expression patterns between individuals of a given sex consistently map to variation in individual quality. Such a pattern, particularly if found in males, would provide unambiguous evidence that the phenotypic response to sexual selection is shaped through sex-specific alterations to the transcriptome. To redress this knowledge gap, we explored whether patterns of sex-biased gene expression are associated with variation in male reproductive quality in Drosophila melanogaster. We measured two male reproductive phenotypes, and their association with sex-biased gene expression, across a selection of inbred lines from the Drosophila Genetic Reference Panel. Genotypes with higher expression of male-biased genes produced males exhibiting shorter latencies to copulation, and higher capacity to inseminate females. Conversely, female-biased genes tended to show negative associations with these male reproductive traits across genotypes. We uncovered similar patterns, by reanalyzing a published dataset from a second D. melanogaster population. Our results reveal the footprint of sexual selection in masculinising the male transcriptome.

Sex-biased expression between guppies varying in the presence of ornamental coloration

Sex-biased gene expression provides a means to achieve sexual dimorphism across a genome largely shared by both sexes. Trinidadian guppies are ideal to examine questions of sex-bias as they exhibit sexual dimorphism in ornamental coloration with male only expression. Here we use RNA-sequencing to quantify whole transcriptome gene expression differences, with a focus on differential expression of color genes between the sexes. We determine whether males express genes positively correlated with coloration at higher levels than females. We find that all the differentially expressed color genes were more highly expressed by males. Males also expressed all known black melanin synthesis genes at higher levels than females, regardless of whether the gene was significantly differentially expressed in the analysis. These differences correlated with the visual color differences between sexes at the stage sampled, as all males had ornamental black coloration apparent. We propose that sexual dimorphism in ornamental coloration is caused by male-biased expression of color genes.

Relative costs and benefits of alternative reproductive phenotypes at different temperatures - genotype-by-environment interactions in a sexually selected trait

Background

The maintenance of considerable genetic variation in sexually selected traits (SSTs) is puzzling given directional selection expected to act on these traits. A possible explanation is the existence of a genotype-by-environment (GxE) interaction for fitness, by which elaborate SSTs are favored in some environments but selected against in others. In the current study, we look for such interactions for fitness-related traits in the bulb mite, a male-dimorphic species with discontinuous expression of a heritable SST in the form of enlarged legs that are used as weapons.

Results

We show that evolution at 18 °C resulted in populations with a higher prevalence of this SST compared to evolution at 24 °C. We further demonstrate that temperature modified male reproductive success in a way that was consistent with these changes. There was a genotype-by-environment interaction for reproductive success – at 18 °C the relative reproductive success of armored males competing with unarmored ones was higher than at the moderate temperature of 24 °C. However, male morph did not have interactive effects with temperature with respect to other life history traits (development time and longevity).

Conclusions

A male genotype that is associated with the expression of a SST interacted with temperature in determining male reproductive success. This interaction caused an elaborate SST to evolve in different directions (more or less prevalent) depending on the thermal environment. The implication of this finding is that seasonal temperature fluctuations have the potential to maintain male polymorphism within populations. Furthermore, spatial heterogeneity in thermal conditions may cause differences among populations in SST selection. This could potentially cause selection against male immigrants from populations in different environments and thus strengthen barriers to gene flow.

Male-limited secondary sexual trait interacts with environment in determining female fitness

Selection for secondary sexual trait (SST) elaboration may increase intralocus sexual conflict over the optimal values of traits expressed from shared genomes. This conflict can reduce female fitness, and the resulting gender load can be exacerbated by environmental stress, with consequences for a population's ability to adapt to novel environments. However, how the evolution of SSTs interacts with environment in determining female fitness is not well understood. Here, we investigated this question using replicate lines of bulb mites selected for increased or decreased prevalence of a male SST-thickened legs used as weapons. The fitness of females from these lines was measured at a temperature to which the mites were adapted (24°C), as well as at two novel temperatures: 18°C and 28°C. We found the prevalence of the SST interacted with temperature in determining female fecundity. At 28°C, females from populations with high SST prevalence were less fecund than females from populations in which the SST was rare, but the reverse was true at 18°C. Thus, a novel environment does not universally depress female fitness more in populations with a high degree of sexually selected dimorphism. We discuss possible consequences of the interaction we detected for adaptation to novel environments.

De novo RNA-seq and functional annotation of Ornithonyssus bacoti

Ornithonyssus bacoti (Hirst) (Acari: Macronyssidae) is a vector and reservoir of pathogens causing serious infectious diseases, such as epidemic hemorrhagic fever, endemic typhus, tularemia, and leptospirosis. Its genome and transcriptome data are lacking in public databases. In this study, total RNA was extracted from live O. bacoti to conduct RNA-seq, functional annotation, coding domain sequence (CDS) prediction and simple sequence repeats (SSRs) detection. The results showed that 65.8 million clean reads were generated and assembled into 72,185 unigenes, of which 49.4% were annotated by seven functional databases. 23,121 unigenes were annotated and assigned to 457 species by non-redundant protein sequence database. The BLAST top-two hit species were Metaseiulus occidentalis and Ixodes scapularis. The procedure detected 12,426 SSRs, of which tri- and di-nucleotides were the most abundant types and the representative motifs were AAT/ATT and AC/GT. 26,936 CDS were predicted with a mean length of 711 bp. 87 unigenes of 30 functional genes, which are usually involved in stress responses, drug resistance, movement, metabolism and allergy, were further identified by bioinformatics methods. The unigenes putatively encoding cytochrome P450 proteins were further analyzed phylogenetically. In conclusion, this study completed the RNA-seq and functional annotation of O. bacoti successfully, which provides reliable molecular data for its future studies of gene function and molecular markers.

Recent advances in understanding the mechanisms of sexually dimorphic plasticity: insights from beetle weapons and future directions

Many traits that are sexually dimorphic, appearing either differently or uniquely in one sex, are also sensitive to an organism's condition. This phenomenon seems to have evolved to limit genetic conflict between traits that are under different selective pressures in each sex. Recent work has shed light on the molecular and developmental mechanisms that govern this condition sensitive growth, and this work has now expanded to encompass both sexual dimorphism as well as conditionally plastic growth, as it seems the two phenomena are linked on a molecular level. In all cases studied the gene doublesex, a conserved regulator of sex differentiation, controls both sexual dimorphism as well as the condition-dependent plastic responses common to these traits. However, the advent of next-generation -omics technologies has allowed researchers to decipher the common and diverged mechanisms of sexually dimorphic plasticity and expand investigations beyond the foundation laid by studies utilizing beetle weapons.

Temporal variation in brain transcriptome is associated with the expression of female mimicry as a sequential male alternative reproductive tactic in fish

Distinct patterns of gene expression often underlie intra- and intersexual differences, and the study of this set of coregulated genes is essential to understand the emergence of complex behavioural phenotypes. Here, we describe the development of a de novo transcriptome and brain gene expression profiles of wild-caught peacock blenny, Salaria pavo, an intertidal fish with sex-role reversal in courtship behaviour (i.e., females are the courting sex) and sequential alternative reproductive tactics in males (i.e., larger and older nest-holder males and smaller and younger sneaker males occur). Sneakers mimic both female's courtship behaviour and nuptial coloration to get access to nests and sneak fertilizations, and later in life transition into nest-holder males. Thus, this species offers the unique opportunity to study how the regulation of gene expression can contribute to intersex phenotypes and to the sequential expression of male and female behavioural phenotypes by the same individual. We found that at the whole brain level, expression of the sneaker tactic was paralleled by broader and divergent gene expression when compared to either females or nest-holder males, which were more similar between themselves. When looking at sex-biased transcripts, sneaker males are intersex rather than being either nest-holder or female-like, and their transcriptome is simultaneously demasculinized for nest-holder-biased transcripts and feminized for female-biased transcripts. These results indicate that evolutionary changes in reproductive plasticity can be achieved through regulation of gene expression, and in particular by varying the magnitude of expression of sex-biased genes, throughout the lifetime of the same individual.

cDNA library construction of two human Demodexspecies

The research of Demodex, a type of pathogen causing various dermatoses in animals and human beings, is lacking at RNA level. This study aims at extracting RNA and constructing cDNA library for Demodex. First, P. cuniculiand D. farinaewere mixed to establish homogenization method for RNA extraction. Second, D. folliculorumand D. breviswere collected and preserved in Trizol, which were mixed with D. farinaerespectively to extract RNA. Finally, cDNA library was constructed and its quality was assessed. The results indicated that for D. folliculorum& D. farinae, the recombination rate of cDNA library was 90.67% and the library titer was 7.50 × 104 pfu/ml. 17 of the 59 positive clones were predicted to be of D. folliculorum; For D. brevis& D. farinae, the recombination rate was 90.96% and the library titer was 7.85 x104 pfu/ml. 40 of the 59 positive clones were predicted to be of D. brevis. Further detection by specific primers demonstrated that mtDNA cox1, cox3and ATP6 detected from cDNA libraries had 96.52%-99.73% identities with the corresponding sequences in GenBank. In conclusion, the cDNA libraries constructed for Demodexmixed with D. farinaewere successful and could satisfy the requirements for functional genes detection.

The transcriptional architecture of phenotypic dimorphism

The profound differences in gene expression between the sexes are increasingly used to study the molecular basis of sexual dimorphism, sexual selection and sexual conflict. Studies of transcriptional architecture, based on comparisons of gene expression, have also been implemented for a wide variety of other intra-specific polymorphisms. These efforts are based on key assumptions regarding the relationship between transcriptional architecture, phenotypic variation and the target of selection. Some of these assumptions are better supported by available evidence than others. In all cases, the evidence is largely circumstantial, leaving considerable gaps in our understanding of the relationship between transcriptional and phenotypic dimorphism.

Selective Landscapes in newt Immune Genes Inferred from Patterns of Nucleotide Variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome.

Inferring regulatory change from gene expression: the confounding effects of tissue scaling

Comparative studies of gene expression are often designed with the aim of identifying regulatory changes associated with phenotypic variation. In recent years, large-scale transcriptome sequencing methods have increasingly been applied to nonmodel organisms to ask important ecological or evolutionary questions. Although experimental design varies, many of these studies have been based on RNA libraries obtained from heterogeneous tissue samples, for example homogenized whole bodies. Comparisons between groups of samples that vary in tissue composition can introduce sufficient variation in RNA abundance to produce patterns of differential expression that are mistakenly interpreted as evidence of regulatory differences. Here, we present a simple model that demonstrates this effect. The model describes the relationship between transcript abundance and tissue composition in a two-tissue system, and how this relationship varies under different scaling relationships. Using a range of biologically realistic variables, including real biological examples, to parameterize the model we highlight the potentially severe influence of tissue scaling on relative transcript abundance. We use these results to identify key aspects of experimental design and analysis that can help to limit the influence of tissue scaling on the inference of regulatory difference from comparative studies of gene expression.

Transcriptomics of Intralocus Sexual Conflict: Gene Expression Patterns in Females Change in Response to Selection on a Male Secondary Sexual Trait in the Bulb Mite

Intralocus sexual conflict (IASC) prevents males and females from reaching their disparate phenotypic optima and is widespread, but little is known about its genetic underpinnings. In Rhizoglyphus robini, a mite species with alternative male morphs, elevated sexual dimorphism of the armored fighter males (compared to more feminized scramblers males) was previously reported to be associated with increased IASC. Because IASC persists if gene expression patterns are correlated between sexes, we compared gene expression patterns of males and females from the replicate lines selected for increased proportion of fighter or scrambler males (F- and S-lines, respectively). Specifically, we tested the prediction that selection for fighter morph caused correlated changes in gene expression patterns in females. We identified 532 differentially expressed genes (FDR < 0.05) between the F-line and S-line males. Consistent with the prediction, expression levels of these genes also differed between females from respective lines. Thus, significant proportion of genes differentially expressed between sexually selected male phenotypes showed correlated expression levels in females, likely contributing to elevated IASC in F-lines reported in a previous study.

Genomic Response to Selection for Predatory Behavior in a Mammalian Model of Adaptive Radiation

If genetic architectures of various quantitative traits are similar, as studies on model organisms suggest, comparable selection pressures should produce similar molecular patterns for various traits. To test this prediction, we used a laboratory model of vertebrate adaptive radiation to investigate the genetic basis of the response to selection for predatory behavior and compare it with evolution of aerobic capacity reported in an earlier work. After 13 generations of selection, the proportion of bank voles (Myodes [=Clethrionomys] glareolus) showing predatory behavior was five times higher in selected lines than in controls. We analyzed the hippocampus and liver transcriptomes and found repeatable changes in allele frequencies and gene expression. Genes with the largest differences between predatory and control lines are associated with hunger, aggression, biological rhythms, and functioning of the nervous system. Evolution of predatory behavior could be meaningfully compared with evolution of high aerobic capacity, because the experiments and analyses were performed in the same methodological framework. The number of genes that changed expression was much smaller in predatory lines, and allele frequencies changed repeatably in predatory but not in aerobic lines. This suggests that more variants of smaller effects underlie variation in aerobic performance, whereas fewer variants of larger effects underlie variation in predatory behavior. Our results thus contradict the view that comparable selection pressures for different quantitative traits produce similar molecular patterns. Therefore, to gain knowledge about molecular-level response to selection for complex traits, we need to investigate not only multiple replicate populations but also multiple quantitative traits.

Genomic heterogeneity of historical gene flow between two species of newts inferred from transcriptome data

The role of gene flow in species formation is a major unresolved issue in speciation biology. Progress in this area requires information on the long-term patterns of gene flow between diverging species. Here, we used thousands of single-nucleotide polymorphisms derived from transcriptome resequencing and a method modeling the joint frequency spectrum of these polymorphisms to reconstruct patterns of historical gene flow between two Lissotriton newts: L. vulgaris (Lv) and L. montandoni (Lm). We tested several models of divergence including complete isolation and various scenarios of historical gene flow. The model of secondary contact received the highest support. According to this model, the species split from their common ancestor ca. 5.5 million years (MY) ago, evolved in isolation for ca. 2 MY, and have been exchanging genes for the last 3.5 MY Demographic changes have been inferred in both species, with the current effective population size of ca. 0.7 million in Lv and 0.2 million in Lm. The postdivergence gene flow resulted in two-directional introgression which affected the genomes of both species, but was more pronounced from Lv to Lm. Interestingly, we found evidence for genomic heterogeneity of interspecific gene flow. This study demonstrates the complexity of long-term gene flow between distinct but incompletely reproductively isolated taxa which divergence was initiated millions of years ago.

Male Sexual Development and the Androgenic Gland: Novel Insights through the de novo Assembled Transcriptome of the Eastern Spiny Lobster, Sagmariasus verreauxi

The Eastern spiny lobster, Sagmariasus verreauxi, is commercially important in fisheries, with growing aquaculture potential, driving an interest to better understand male sexual differentiation. Amongst the Decapoda, the androgenic gland (AG) and the insulin-like androgenic gland hormone (IAG) have a well-defined function in male sexual differentiation. However, IAG is not a sex determinant and therefore must be considered as part of a broader, integrated pathway. This work uses a transcriptomic, multi-tissue approach to provide an integrated description of male-biased expression as mediated through the AG. Transcriptomic analyses demonstrate that IAG expression is stage- and eyestalk-regulated (low in immature, high in mature and 6-times higher in hypertrophied glands), with IAG being the predominant AG-specific factor. The low expression of this key factor in immature males suggests the involvement of other tissues in male sexual differentiation. Across tissues, the gonad (87.8%) and antennal gland (73.5%) show the highest male-biased differential expression of transcripts and also express 4 sex-determination regulators, known as Dmrts, with broader expression of Sv-Sxl and Sv-TRA-2. In order to better understand male sexual differentiation, tissues other than the AG must also be considered. This research highlights the gonad and antennal gland as showing significant male-biased expression patterns, including the Sv-Dmrts.

A Paradox of Genetic Variance in Epigamic Traits: Beyond "Good Genes" View of Sexual Selection

Maintenance of genetic variance in secondary sexual traits, including bizarre ornaments and elaborated courtship displays, is a central problem of sexual selection theory. Despite theoretical arguments predicting that strong sexual selection leads to a depletion of additive genetic variance, traits associated with mating success show relatively high heritability. Here we argue that because of trade-offs associated with the production of costly epigamic traits, sexual selection is likely to lead to an increase, rather than a depletion, of genetic variance in those traits. Such trade-offs can also be expected to contribute to the maintenance of genetic variation in ecologically relevant traits with important implications for evolutionary processes, e.g. adaptation to novel environments or ecological speciation. However, if trade-offs are an important source of genetic variation in sexual traits, the magnitude of genetic variation may have little relevance for the possible genetic benefits of mate choice.

The locus of sexual selection: moving sexual selection studies into the post-genomics era

Sexual selection drives fundamental evolutionary processes such as trait elaboration and speciation. Despite this importance, there are surprisingly few examples of genes unequivocally responsible for variation in sexually selected phenotypes. This lack of information inhibits our ability to predict phenotypic change due to universal behaviours, such as fighting over mates and mate choice. Here, we discuss reasons for this apparent gap and provide recommendations for how it can be overcome by adopting contemporary genomic methods, exploiting underutilized taxa that may be ideal for detecting the effects of sexual selection and adopting appropriate experimental paradigms. Identifying genes that determine variation in sexually selected traits has the potential to improve theoretical models and reveal whether the genetic changes underlying phenotypic novelty utilize common or unique molecular mechanisms. Such a genomic approach to sexual selection will help answer questions in the evolution of sexually selected phenotypes that were first asked by Darwin and can furthermore serve as a model for the application of genomics in all areas of evolutionary biology.

Data concatenation, Bayesian concordance and coalescent-based analyses of the species tree for the rapid radiation of Triturus newts

The phylogenetic relationships for rapid species radiations are difficult to disentangle. Here we study one such case, namely the genus Triturus, which is composed of the marbled and crested newts. We analyze data for 38 genetic markers, positioned in 3-prime untranslated regions of protein-coding genes, obtained with 454 sequencing. Our dataset includes twenty Triturus newts and represents all nine species. Bayesian analysis of population structure allocates all individuals to their respective species. The branching patterns obtained by data concatenation, Bayesian concordance analysis and coalescent-based estimations of the species tree differ from one another. The data concatenation based species tree shows high branch support but branching order is considerably affected by allele choice in the case of heterozygotes in the concatenation process. Bayesian concordance analysis expresses the conflict between individual gene trees for part of the Triturus species tree as low concordance factors. The coalescent-based species tree is relatively similar to a previously published species tree based upon morphology and full mtDNA and any conflicting internal branches are not highly supported. Our findings reflect high gene tree discordance due to incomplete lineage sorting (possibly aggravated by hybridization) in combination with low information content of the markers employed (as can be expected for relatively recent species radiations). This case study highlights the complexity of resolving rapid radiations and we acknowledge that to convincingly resolve the Triturus species tree even more genes will have to be consulted.

Single nucleotide polymorphisms reveal genetic structuring of the carpathian newt and provide evidence of interspecific gene flow in the nuclear genome

Genetic variation within species is commonly structured in a hierarchical manner which may result from superimposition of processes acting at different spatial and temporal scales. In organisms of limited dispersal ability, signatures of past subdivision are detectable for a long time. Studies of contemporary genetic structure in such taxa inform about the history of isolation, range changes and local admixture resulting from geographically restricted hybridization with related species. Here we use a set of 139 transcriptome-derived, unlinked nuclear single nucleotide polymorphisms (SNP) to assess the genetic structure of the Carpathian newt (Lissotriton montandoni, Lm) and introgression from its congener, the smooth newt (L. vulgaris, Lv). Two substantially differentiated groups of Lm populations likely originated from separate refugia, both located in the Eastern Carpathians. The colonization of the present range in north-western and south-western directions was accompanied by a modest loss of variation; admixture between the two groups has occurred in the middle of the Eastern Carpathians. Local, apparently recent introgression of Lv alleles into several Lm populations was detected, demonstrating increased power for admixture detection in comparison to a previous study based on a limited number of microsatellite markers. The level of introgression was higher in Lm populations classified as admixed than in syntopic populations. We discuss the possible causes and propose further tests to distinguish between alternatives. Several outlier loci were identified in tests of interspecific differentiation, suggesting genomic heterogeneity of gene flow between species.

Selection for alternative male reproductive tactics alters intralocus sexual conflict

Intralocus sexual conflict (IASC) arises when fitness optima for a shared trait differ between the sexes; such conflict may help maintain genetic variation within populations. Sex-limited expression of sexually antagonistic traits may help resolve the conflict, but the extent of this resolution remains a subject of debate. In species with alternative male reproductive tactics, unresolved conflict should manifest more in a more sexually dimorphic male phenotype. We tested this prediction in the bulb mite (Rhizoglyphus robini), a species in which aggressive fighters coexist with benign scramblers. To do this, we established replicated lines in which we increased the proportion of each of the alternative male morphs using artificial selection. After approximately 40 generations, the proportion of fighters and scramblers stabilized at >0.9 in fighter- and scrambler-selected lines, respectively. We then measured several female fitness components. As predicted by IASC theory, female fecundity and longevity were lower in lines selected for fighters and higher in lines selected for scramblers. This finding indicates that sexually selected phenotypes are associated with an ontogenetic conflict that is not easily resolved. Furthermore, we suggest that IASC may be an important mechanism contributing to the maintenance of genetic variation in the expression of alternative reproductive tactics.

Parallel tagged amplicon sequencing of transcriptome-based genetic markers for Triturus newts with the Ion Torrent next-generation sequencing platform

Next-generation sequencing is a fast and cost-effective way to obtain sequence data for nonmodel organisms for many markers and for many individuals. We describe a protocol through which we obtain orthologous markers for the crested newts (Amphibia: Salamandridae: Triturus), suitable for analysis of interspecific hybridization. We use transcriptome data of a single Triturus species and design 96 primer pairs that amplify c. 180 bp fragments positioned in 3-prime untranslated regions. Next, these markers are tested with uniplex PCR for a set of species spanning the taxonomical width of the genus Triturus. The 52 markers that consistently show a single band of expected length at gel electrophoreses for all tested crested newt species are then amplified in five multiplex PCRs (with a plexity of ten or eleven) for 132 individual newts: a set of 84 representing the seven (candidate) species and a set of 48 from a presumed hybrid population. After pooling multiplexes per individual, unique tags are ligated to link amplicons to individuals. Subsequently, individuals are pooled equimolar and sequenced on the Ion Torrent next-generation sequencing platform. A bioinformatics pipeline identifies the alleles and recodes these to a genotypic format. Next, we test the utility of our markers. baps allocates the 84 crested newt individuals representing (candidate) species to their expected (candidate) species, confirming the markers are suitable for species delineation. newhybrids, a hybrid index and hiest confirm the 48 individuals from the presumed hybrid population to be genetically admixed, illustrating the potential of the markers to identify interspecific hybridization. We expect the set of markers we designed to provide a high resolving power for analysis of hybridization in Triturus.