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.

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.

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.

Nitrite regulates hypoxic vasodilation via myoglobin-dependent nitric oxide generation.

BACKGROUND: Hypoxic vasodilation is a physiological response to low oxygen tension that increases blood supply to match metabolic demands. Although this response has been characterized for >100 years, the underlying hypoxic sensing and effector signaling mechanisms remain uncertain. We have shown that deoxygenated myoglobin in the heart can reduce nitrite to nitric oxide (NO·) and thereby contribute to cardiomyocyte NO· signaling during ischemia. On the basis of recent observations that myoglobin is expressed in the vasculature of hypoxia-tolerant fish, we hypothesized that endogenous nitrite may contribute to physiological hypoxic vasodilation via reactions with vascular myoglobin to form NO·. METHODS AND RESULTS: We show in the present study that myoglobin is expressed in vascular smooth muscle and contributes significantly to nitrite-dependent hypoxic vasodilation in vivo and ex vivo. The generation of NO· from nitrite reduction by deoxygenated myoglobin activates canonical soluble guanylate cyclase/cGMP signaling pathways. In vivo and ex vivo vasodilation responses, the reduction of nitrite to NO·, and the subsequent signal transduction mechanisms were all significantly impaired in mice without myoglobin. Hypoxic vasodilation studies in myoglobin and endothelial and inducible NO synthase knockout models suggest that only myoglobin contributes to systemic hypoxic vasodilatory responses in mice. CONCLUSIONS: Endogenous nitrite is a physiological effector of hypoxic vasodilation. Its reduction to NO· via the heme globin myoglobin enhances blood flow and matches O(2) supply to increased metabolic demands under hypoxic conditions.

Functional differentiation of myoglobin isoforms in hypoxia-tolerant carp indicates tissue-specific protective roles.

Because of a recent whole genome duplication, the hypoxia-tolerant common carp and goldfish are the only vertebrates known to possess two myoglobin (Mb) paralogs. One of these, Mb1, occurs in oxidative muscle but also in several other tissues, including capillary endothelial cells, whereas the other, Mb2, is a unique isoform specific to brain neurons. To help understand the functional roles of these diverged isoforms in the tolerance to severe hypoxia in the carp, we have compared their O(2) equilibria, carbon monoxide (CO) and O(2) binding kinetics, thiol S-nitrosation, nitrite reductase activities, and peroxidase activities. Mb1 has O(2) affinity and nitrite reductase activity comparable to most vertebrate muscle Mbs, consistent with established roles for Mbs in O(2) storage/delivery and in maintaining nitric oxide (NO) homeostasis during hypoxia. Both Mb1 and Mb2 can be S-nitrosated to similar extent, but without oxygenation-linked allosteric control. When compared with Mb1, Mb2 displays faster O(2) and CO kinetics, a lower O(2) affinity, and is slower at converting nitrite into NO. Mb2 is therefore unlikely to be primarily involved in either O(2) supply to mitochondria or the generation of NO from nitrite during hypoxia. However, Mb2 proved to be significantly faster at eliminating H(2)O(2,) a major in vivo reactive oxygen species (ROS), suggesting that this diverged Mb isoform may have a specific protective role against H(2)O(2) in the carp brain. This property might be of particular significance during reoxygenation following extended periods of hypoxia, when production of H(2)O(2) and other ROS is highest.

Identification of candidate genes and physiological pathways involved in gonad deformation in whitefish (Coregonus spp.) from Lake Thun, Switzerland.

In 2000, fishermen reported the appearance of deformed reproductive organs in whitefish (Coregonus spp.) from Lake Thun, Switzerland. Despite intensive investigations, the causes of these abnormalities remain unknown. Using gene expression profiling, we sought to identify candidate genes and physiological processes possibly associated with the observed gonadal deformations, in order to gain insights into potential causes. Using in situ-synthesized oligonucleotide arrays, we compared the expression levels at 21,492 unique transcript probes in liver and head kidney tissue of male whitefish with deformed and normally developed gonads, respectively. The fish had been collected on spawning sites of two genetically distinct whitefish forms of Lake Thun. We contrasted the gene expression profiles of 56 individuals, i.e., 14 individuals of each phenotype and of each population. Gene-by-gene analysis revealed weak expression differences between normal and deformed fish, and only one gene, ictacalcin, was found to be up-regulated in head kidney tissue of deformed fish from both whitefish forms, However, this difference could not be confirmed with quantitative real-time qPCR. Enrichment analysis on the level of physiological processes revealed (i) the involvement of immune response genes in both tissues, particularly those linked to complement activation in the liver, (ii) proteolysis in the liver and (iii) GTPase activity and Ras protein signal transduction in the head kidney. In comparison with current literature, this gene expression pattern signals a chronic autoimmune disease in the testes. Based on the recent observations that gonad deformations are induced through feeding of zooplankton from Lake Thun we hypothesize that a xenobiotic accumulated in whitefish via the plankton triggering autoimmunity as the likely cause of gonad deformations. We propose several experimental strategies to verify or reject this hypothesis.

Molecular correlates of social dominance: a novel role for ependymin in aggression.

Theoretical and empirical studies have sought to explain the formation and maintenance of social relationships within groups. The resulting dominance hierarchies have significant fitness and survival consequences dependent upon social status. We hypothesised that each position or rank within a group has a distinctive brain gene expression profile that correlates with behavioural phenotype. Furthermore, transitions in rank position should determine which genes shift in expression concurrent with the new dominance status. We used a custom cDNA microarray to profile brain transcript expression in a model species, the rainbow trout, which forms tractable linear hierarchies. Dominant, subdominant and submissive individuals had distinctive transcript profiles with 110 gene probes identified using conservative statistical analyses. By removing the dominant, we characterised the changes in transcript expression in sub-dominant individuals that became dominant demonstrating that the molecular transition occurred within 48 hours. A strong, novel candidate gene, ependymin, which was highly expressed in both the transcript and protein in subdominants relative to dominants, was tested further. Using antibody injection to inactivate ependymin in pairs of dominant and subdominant zebrafish, the subdominant fish exhibited a substantial increase in aggression in parallel with an enhanced competitive ability. This is the first study to characterise the molecular signatures of dominance status within groups and the first to implicate ependymin in control of aggressive behaviour. It also provides evidence for indirect genetic effect models in which genotype/phenotype of an individual is influenced by conspecific interactions within a group. The variation in the molecular profile of each individual within a group may offer a new explanation of intraspecific variation in gene expression within undefined groups of animals and provides new candidates for empirical study.

An information-rich alternative, chemicals testing strategy using a high definition toxicogenomics and zebrafish (Danio rerio) embryos.

Large-scale toxicogenomic screening approaches offer great promise for generating a bias-free system-wide view of toxicological effects and modes-of-action of chemicals and ecotoxicants. However, early applications of microarray technology have identified relatively small groups of responding genes with which to define new targets for analysis by conventional means. We have trialled a more intensive approach to the design and interpretation of array experiments incorporating a balanced interwoven ANOVA design with higher levels of biological replication, a more thorough analysis of errors and false discovery rates, and an analysis of response patterns using gene network models. Zebrafish embryos were exposed from 1.5 h post-fertilization for 72 h to ecotoxicants representing different classes--2,4-dichlorophenol, 3,4-dichloroaniline, pentachlorophenol, and cadmium chloride--at low concentrations producing a developmental disturbance to 10% of embryos and half of this dose. Extracted whole embryo RNA was then analyzed on microarrays. Analysis revealed responses of 3000-5000 genes, which is 10-1000 times greater than previously reported, with significance at lower levels of fold change. Some gene responses were common to multiple toxicants, and others were restricted to just one or two toxicants. The gene expression profiles for the different toxicants were distinctive, and analysis using network-based models provided a high level of detail of affected processes, some of which were novel. This approach provides a more highly refined view of toxic effects, from which meaningful patterns of response can be discerned and related to functional deficits and from which more reliable indicators of toxicological effect can be predicted.