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.

Acute and chronic stress prevents responses to pain in zebrafish: evidence for stress-induced analgesia.

The state of an animal prior to the application of a noxious stimulus can have a profound effect on their nociceptive threshold and subsequent behaviour. In mammals, the presence of acute stress preceding a painful event can have an analgesic effect whereas the presence of chronic stress can result in hyperalgesia. While considerable research has been conducted on the ability of stress to modulate mammalian responses to pain, relatively little is known about fish. This is of particular concern given that zebrafish (Danio rerio) are an extensively used model organism subject to a wide array of invasive procedures where the level of stress prior to experimentation could pose a major confounding factor. This study, therefore, investigated the impact of both acute and chronic stress on the behaviour of zebrafish subjected to a potentially painful laboratory procedure, the fin clip. In stress-free individuals, those subjected to the fin clip spent more time in the bottom of the tank, had reduced swimming speeds and less complex swimming trajectories; however, these behavioural changes were absent in fin-clipped fish that were first subject to either chronic or acute stress, suggesting the possibility of stress-induced analgesia (SIA). To test this, the opioid antagonist naloxone was administered to fish prior to the application of both the stress and fin-clip procedure. After naloxone, acutely stressed fin-clipped zebrafish exhibited the same behaviours as stress-free fin-clipped fish. This indicates the presence of SIA and the importance of opioid signalling in this mechanism. As stress reduced nociceptive responses in zebrafish, this demonstrates the potential for an endogenous analgesic system akin to the mammalian system. Future studies should delineate the neurobiological basis of stress-induced analgesia in fish.

Automated monitoring of behaviour in zebrafish after invasive procedures.

Fish are used in a variety of experimental contexts often in high numbers. To maintain their welfare and ensure valid results during invasive procedures it is vital that we can detect subtle changes in behaviour that may allow us to intervene to provide pain-relief. Therefore, an automated method, the Fish Behaviour Index (FBI), was devised and used for testing the impact of laboratory procedures and efficacy of analgesic drugs in the model species, the zebrafish. Cameras with tracking software were used to visually track and quantify female zebrafish behaviour in real time after a number of laboratory procedures including fin clipping, PIT tagging, and nociceptor excitation via injection of acetic acid subcutaneously. The FBI was derived from activity and distance swum measured before and after these procedures compared with control and sham groups. Further, the efficacy of a range of drugs with analgesic properties to identify efficacy of these agents was explored. Lidocaine (5 mg/L), flunixin (8 mg/L) and morphine (48 mg/L) prevented the associated reduction in activity and distance swum after fin clipping. From an ethical perspective, the FBI represents a significant refinement in the use of zebrafish and could be adopted across a wide range of biological disciplines.

Host selectively contributes to shaping intestinal microbiota of carnivorous and omnivorous fish.

Fish production is increasingly important to global food security. A major factor in maintaining health, productivity and welfare of farmed fish is the establishment and promotion of a stable and beneficial intestinal microbiota. Understanding the effects of factors such as host and environment on gut microbial community structure is essential for developing strategies for stimulating the establishment of a health-promoting gut-microbiota. We compared intestinal microbiota of common carp and rainbow trout, two fish with different dietary habits, sourced from various farm locations. There were distinct differences in the gut microbiota of carp and trout intestine. The microbiota of carp was dominated by Fusobacteriia and Gammaproteobacteria, while the trout microbiota consisted predominantly of Mollicutes and Betaproteobacteria. The majority of bacterial sequences clustered into a relatively low number of operational taxonomic units (OTUs) revealing a comparatively simple microbiota, with Cetobacterium, Aeromonas and Mycoplasma being highly abundant. Within each species, fish from different facilities were found to have markedly similar predominant bacterial populations despite distinctly different rearing environments, demonstrating intra-species uniformity and significant influence of host selectivity. This study demonstrates that in fish the host species imparts substantial impact in shaping the community structure of the intestinal microbiota.

Mating system manipulation and the evolution of sex-biased gene expression in Drosophila.

Sex differences in dioecious animals are pervasive and result from gene expression differences. Elevated sexual selection has been predicted to increase the number and expression of male-biased genes, and experimentally imposing monogamy on Drosophila melanogaster has led to a relative feminisation of the transcriptome. Here, we test this hypothesis further by subjecting another polyandrous species, D. pseudoobscura, to 150 generations of experimental monogamy or elevated polyandry. We find that sex-biased genes do change in expression but, contrary to predictions, there is usually masculinisation of the transcriptome under monogamy, although this depends on tissue and sex. We also identify and describe gene expression changes following courtship experience. Courtship often influences gene expression, including patterns in sex-biased gene expression. Our results confirm that mating system manipulation disproportionately influences sex-biased gene expression but show that the direction of change is dynamic and unpredictable.

Contribution of trans regulatory eQTL to cryptic genetic variation in C. elegans.

BACKGROUND: Cryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population. RESULTS: We measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40-57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL. CONCLUSIONS: These results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.

The Role of Omics in the Application of Adverse Outcome Pathways for Chemical Risk Assessment.

In conjunction with the second International Environmental Omics Symposium (iEOS) conference, held at the University of Liverpool (United Kingdom) in September 2014, a workshop was held to bring together experts in toxicology and regulatory science from academia, government and industry. The purpose of the workshop was to review the specific roles that high-content omics datasets (eg, transcriptomics, metabolomics, lipidomics, and proteomics) can hold within the adverse outcome pathway (AOP) framework for supporting ecological and human health risk assessments. In light of the growing number of examples of the application of omics data in the context of ecological risk assessment, we considered how omics datasets might continue to support the AOP framework. In particular, the role of omics in identifying potential AOP molecular initiating events and providing supportive evidence of key events at different levels of biological organization and across taxonomic groups was discussed. Areas with potential for short and medium-term breakthroughs were also discussed, such as providing mechanistic evidence to support chemical read-across, providing weight of evidence information for mode of action assignment, understanding biological networks, and developing robust extrapolations of species-sensitivity. Key challenges that need to be addressed were considered, including the need for a cohesive approach towards experimental design, the lack of a mutually agreed framework to quantitatively link genes and pathways to key events, and the need for better interpretation of chemically induced changes at the molecular level. This article was developed to provide an overview of ecological risk assessment process and a perspective on how high content molecular-level datasets can support the future of assessment procedures through the AOP framework.

Transcriptomic signatures differentiate survival from fatal outcomes in humans infected with Ebola virus.

BACKGROUND: In 2014, Western Africa experienced an unanticipated explosion of Ebola virus infections. What distinguishes fatal from non-fatal outcomes remains largely unknown, yet is key to optimising personalised treatment strategies. We used transcriptome data for peripheral blood taken from infected and convalescent recovering patients to identify early stage host factors that are associated with acute illness and those that differentiate patient survival from fatality. RESULTS: The data demonstrate that individuals who succumbed to the disease show stronger upregulation of interferon signalling and acute phase responses compared to survivors during the acute phase of infection. Particularly notable is the strong upregulation of albumin and fibrinogen genes, which suggest significant liver pathology. Cell subtype prediction using messenger RNA expression patterns indicated that NK-cell populations increase in patients who survive infection. By selecting genes whose expression properties discriminated between fatal cases and survivors, we identify a small panel of responding genes that act as strong predictors of patient outcome, independent of viral load. CONCLUSIONS: Transcriptomic analysis of the host response to pathogen infection using blood samples taken during an outbreak situation can provide multiple levels of information on both disease state and mechanisms of pathogenesis. Host biomarkers were identified that provide high predictive value under conditions where other predictors, such as viral load, are poor prognostic indicators. The data suggested that rapid analysis of the host response to infection in an outbreak situation can provide valuable information to guide an understanding of disease outcome and mechanisms of disease.

Transcriptome sequencing of human breast cancer reveals aberrant intronic transcription in amplicons and dysregulation of alternative splicing with major therapeutic implications.

Advances in genomic and transcriptome sequencing are revealing the massive scale of previously unrecognised alterations occurring during neoplastic transformation. Breast cancers are genetically and phenotypically heterogeneous. Each of the three major subtypes [ERBB2 amplified, estrogen receptor (ESR)-positive and triple-negative] poses diagnostic and therapeutic challenges. Here we show, using high-resolution next-generation transcriptome sequencing, that in all three breast cancer subtypes, but not matched controls, there was significant overexpression of transcripts from intronic and untranslated regions in addition to exons from specific genes, particularly amplified oncogenes and hormone receptors. For key genes ERBB2 and ESR1, we demonstrate that overexpression is linked to the production of highly modified and truncated splice variants in tumours, but not controls, correlated with tumour subtype. Translation of these tumour-specific splice variants generates truncated proteins with altered subcellular locations and functions, modifying the phenotype, affecting tumour biology, and targeted antitumour therapies. In contrast, tumour suppressors TP53, BRCA1/2 and NF1 did not show intronic overexpression or truncated splice variants in cancers. These findings emphasize the detection of intronic as well as exonic changes in the transcriptional landscapes of cancers have profound therapeutic implications.

Remarkably Divergent Regions Punctuate the Genome Assembly of the Caenorhabditis elegans Hawaiian Strain CB4856.

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion-deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.

Life without oxygen: gene regulatory responses of the crucian carp (Carassius carassius) heart subjected to chronic anoxia.

Crucian carp are unusual among vertebrates in surviving extended periods in the complete absence of molecular oxygen. During this time cardiac output is maintained though these mechanisms are not well understood. Using a high-density cDNA microarray, we have defined the genome-wide gene expression responses of cardiac tissue after exposing the fish at two temperatures (8 and 13 °C) to one and seven days of anoxia, followed by seven days after restoration to normoxia. At 8 °C, using a false discovery rate of 5%, neither anoxia nor re-oxygenation elicited appreciable changes in gene expression. By contrast, at 13 °C, 777 unique genes responded strongly. Up-regulated genes included those involved in protein turnover, the pentose phosphate pathway and cell morphogenesis while down-regulated gene categories included RNA splicing and transcription. Most genes were affected between one and seven days of anoxia, indicating gene regulation over the medium term but with few early response genes. Re-oxygenation for 7 days was sufficient to completely reverse these responses. Glycolysis displayed more complex responses with anoxia up-regulated transcripts for the key regulatory enzymes, hexokinase and phosphofructokinase, but with down-regulation of most of the non-regulatory genes. This complex pattern of responses in genomic transcription patterns indicates divergent cardiac responses to anoxia, with the transcriptionally driven reprogramming of cardiac function seen at 13 °C being largely completed at 8 °C.

Single-step selection of bivalent aptamers validated by comparison with SELEX using high-throughput sequencing.

The identification of nucleic acid aptamers would be advanced if they could be obtained after fewer rounds of selection and amplification. In this paper the identification of bivalent aptamers for thrombin by SELEX and single-step selection are compared using next generation sequencing and motif finding informatics. Results show that similar aptamers are identified by both methods. This is significant because it shows that next generation sequencing and motif finding informatics have the potential to simplify the selection of aptamers by avoiding multiple rounds of enzymatic transcription and amplification.

A rapid and massive gene expression shift marking adolescent transition in C. elegans.

Organismal development is the most dynamic period of the life cycle, yet we have only a rough understanding of the dynamics of gene expression during adolescent transition. Here we show that adolescence in Caenorhabditis elegans is characterized by a spectacular expression shift of conserved and highly polymorphic genes. Using a high resolution time series we found that in adolescent worms over 10,000 genes changed their expression. These genes were clustered according to their expression patterns. One cluster involved in chromatin remodelling showed a brief up-regulation around 50 h post-hatch. At the same time a spectacular shift in expression was observed. Sequence comparisons for this cluster across many genotypes revealed diversifying selection. Strongly up-regulated genes showed signs of purifying selection in non-coding regions, indicating that adolescence-active genes are constrained on their regulatory properties. Our findings improve our understanding of adolescent transition and help to eliminate experimental artefacts due to incorrect developmental timing.

Molecular basis of chill resistance adaptations in poikilothermic animals.

Chill and freeze represent very different components of low temperature stress. Whilst the principal mechanisms of tissue damage and of acquired protection from freeze-induced effects are reasonably well established, those for chill damage and protection are not. Non-freeze cold exposure (i.e. chill) can lead to serious disruption to normal life processes, including disruption to energy metabolism, loss of membrane perm-selectivity and collapse of ion gradients, as well as loss of neuromuscular coordination. If the primary lesions are not relieved then the progressive functional debilitation can lead to death. Thus, identifying the underpinning molecular lesions can point to the means of building resistance to subsequent chill exposures. Researchers have focused on four specific lesions: (i) failure of neuromuscular coordination, (ii) perturbation of bio-membrane structure and adaptations due to altered lipid composition, (iii) protein unfolding, which might be mitigated by the induced expression of compatible osmolytes acting as 'chemical chaperones', (iv) or the induced expression of protein chaperones along with the suppression of general protein synthesis. Progress in all these potential mechanisms has been ongoing but not substantial, due in part to an over-reliance on straightforward correlative approaches. Also, few studies have intervened by adoption of single gene ablation, which provides much more direct and compelling evidence for the role of specific genes, and thus processes, in adaptive phenotypes. Another difficulty is the existence of multiple mechanisms, which often act together, thus resulting in compensatory responses to gene manipulations, which may potentially mask disruptive effects on the chill tolerance phenotype. Consequently, there is little direct evidence of the underpinning regulatory mechanisms leading to induced resistance to chill injury. Here, we review recent advances mainly in lower vertebrates and in arthropods, but increasingly in genetic model species from a broader range of taxa.

Transcriptome-wide expression variation associated with environmental plasticity and mating success in cactophilic Drosophila mojavensis.

Ecological speciation occurs with the adaptation of populations to different environments and concurrent evolution of reproductive isolation. Phenotypic plasticity might influence both ecological adaptation and reproductive traits. We examined environment-specific gene expression and male mating success in cactophilic Drosophila mojavensis using transcriptome sequencing. This species exhibits cactus-dependent mating success across different species of host plants, with genotype-by-environment interactions for numerous traits. We cultured flies from egg to eclosion on two natural cactus hosts and surveyed gene expression in adult males that were either successful or unsuccessful in achieving copulation in courtship trials. We identified gene expression differences that included functions involved with metabolism, most likely related to chemical differences between host cactus species. Several epigenetic-related functions were identified that might play a role in modulating gene expression in adults due to host cactus effects on larvae, and mating success. Cactus-dependent mating success involved expression differences of genes implicated in translation, transcription, and nervous system development. This suggests a role of neurological function genes in the mating success of D. mojavensis males. Together, these results suggest that the influence of environmental variation on mating success via regulation of gene expression might be an important aspect of ecological speciation.

Evolution of mammalian diving capacity traced by myoglobin net surface charge.

Extended breath-hold endurance enables the exploitation of the aquatic niche by numerous mammalian lineages and is accomplished by elevated body oxygen stores and adaptations that promote their economical use. However, little is known regarding the molecular and evolutionary underpinnings of the high muscle myoglobin concentration phenotype of divers. We used ancestral sequence reconstruction to trace the evolution of this oxygen-storing protein across a 130-species mammalian phylogeny and reveal an adaptive molecular signature of elevated myoglobin net surface charge in diving species that is mechanistically linked with maximal myoglobin concentration. This observation provides insights into the tempo and routes to enhanced dive capacity evolution within the ancestors of each major mammalian aquatic lineage and infers amphibious ancestries of echidnas, moles, hyraxes, and elephants, offering a fresh perspective on the evolution of this iconic respiratory pigment.

Transcriptomics and in vivo tests reveal novel mechanisms underlying endocrine disruption in an ecological sentinel, Nucella lapillus.

Anthropogenic endocrine disruptors now contaminate all environments globally, with concomitant deleterious effects across diverse taxa. While most studies on endocrine disruption (ED) have focused on vertebrates, the superimposition of male sexual characteristics in the female dogwhelk, Nucella lapillus (imposex), caused by organotins, provides one of the most clearcut ecological examples of anthropogenically induced ED in aquatic ecosystems. To identify the underpinning mechanisms of imposex for this 'nonmodel' species, we combined Roche 454 pyrosequencing with custom oligoarray fabrication inexpensively to both generate gene models and identify those responding to chronic tributyltin (TBT) treatment. The results supported the involvement of steroid, neuroendocrine peptide hormone dysfunction and retinoid mechanisms, but suggested additionally the involvement of putative peroxisome proliferator-activated receptor (PPAR) pathways. Application of rosiglitazone, a well-known vertebrate PPARγ ligand, to dogwhelks induced imposex in the absence of TBT. Thus, while TBT-induced imposex is linked to the induction of many genes and has a complex phenotype, it is likely also to be driven by PPAR-responsive pathways, hitherto not described in invertebrates. Our findings provide further evidence for a common signalling pathway between invertebrate and vertebrate species that has previously been overlooked in the study of endocrine disruption.

The selection of DNA aptamers for two different epitopes of thrombin was not due to different partitioning methods.

Nearly all aptamers identified so far for any given target molecule have been specific for the same binding site (epitope). The most notable exception to the 1 aptamer per target molecule rule is the pair of DNA aptamers that bind to different epitopes of thrombin. This communication refutes the suggestion that these aptamers exist because different partitioning methods were used when they were selected. The possibility that selection of these aptamers was biased by conflicting secondary structures was also investigated and found not to contribute. The preparation of protein-coated magnetic beads for systematic evolution of ligands by exponential enrichment (SELEX) and the different specificities of the thrombin aptamers for the α and ß forms of thrombin are also reported.