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.

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.

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.

ExprAlign--the identification of ESTs in non-model species by alignment of cDNA microarray expression profiles.

BACKGROUND: Sequence identification of ESTs from non-model species offers distinct challenges particularly when these species have duplicated genomes and when they are phylogenetically distant from sequenced model organisms. For the common carp, an environmental model of aquacultural interest, large numbers of ESTs remained unidentified using BLAST sequence alignment. We have used the expression profiles from large-scale microarray experiments to suggest gene identities. RESULTS: Expression profiles from approximation 700 cDNA microarrays describing responses of 7 major tissues to multiple environmental stressors were used to define a co-expression landscape. This was based on the Pearsons correlation coefficient relating each gene with all other genes, from which a network description provided clusters of highly correlated genes as 'mountains'. We show that these contain genes with known identities and genes with unknown identities, and that the correlation constitutes evidence of identity in the latter. This procedure has suggested identities to 522 of 2701 unknown carp ESTs sequences. We also discriminate several common carp genes and gene isoforms that were not discriminated by BLAST sequence alignment alone. Precision in identification was substantially improved by use of data from multiple tissues and treatments. CONCLUSION: The detailed analysis of co-expression landscapes is a sensitive technique for suggesting an identity for the large number of BLAST unidentified cDNAs generated in EST projects. It is capable of detecting even subtle changes in expression profiles, and thereby of distinguishing genes with a common BLAST identity into different identities. It benefits from the use of multiple treatments or contrasts, and from the large-scale microarray data.

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.

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.

Ancient and modern duplication events and the evolution of stearoyl-CoA desaturases in teleost fishes.

Stearoyl-CoA desaturases (SCDs) are key enzymes of fatty acid biosynthesis whose regulation underpins responses to dietary, thermal, and hormonal treatment. Although two isoforms are known to exist in the common carp and human and four in mouse, there is no coherent view on how this gene family evolved to generate functionally diverse members. Here we identify numerous new SCD homologs in teleost fishes, using sequence data from expressed sequence tag (EST) and cDNA collections and genomic model species. Phylogenetic analyses of the deduced coding sequences produced only partially resolved molecular trees. The multiple SCD isoforms were, however, consistent with having arisen by an ancient gene duplication event in teleost fishes together with a more recent duplication in the tetraploid carp and possibly also salmonid lineages. Critical support for this interpretation comes from comparison across all vertebrate groups of the gene order in the genomic environments of the SCD isoforms. Using syntenically aligned chromosomal fragments from large-insert clones of common carp and grass carp together with those from genomically sequenced model species, we show that the ancient and modern SCD duplication events in the carp lineage were each associated with large chromosomal segment duplications, both possibly linked to whole genome duplications. By contrast, the four mouse isoforms likely arose by tandem duplications. Each duplication in the carp lineage gave rise to differentially expressed SCD isoforms, either induced by cold or diet as previously shown for the recent duplicated carp isoforms or tissue specific as demonstrated here for the ancient duplicate zebrafish isoforms.

Effects of fluoxetine on the reproductive axis of female goldfish (Carassius auratus).

We investigated the effects of fluoxetine, a selective serotonin reuptake inhibitor, on neuroendocrine function and the reproductive axis in female goldfish. Fish were given intraperitoneal injections of fluoxetine twice a week for 14 days, resulting in five injections of 5 microg fluoxetine/g body wt. We measured the monoamine neurotransmitters serotonin, dopamine, and norepinephrine in addition to their metabolites with HPLC. Homovanillic acid, a metabolite in the dopaminergic pathway, increased significantly in the hypothalamus. Plasma estradiol levels were measured by radioimmunoassay and were significantly reduced approximately threefold after fluoxetine treatment. We found that fluoxetine also significantly reduced the expression of estrogen receptor (ER)beta1 mRNA by 4-fold in both the hypothalamus and the telencephalon and ERalpha mRNA by 1.7-fold in the telencephalon. Fluoxetine had no effect on the expression of ERbeta2 mRNA in the hypothalamus or telencephalon. Microarray analysis identified isotocin, a neuropeptide that stimulates reproductive behavior in fish, as a candidate gene affected by fluoxetine treatment. Real-time RT-PCR verified that isotocin mRNA was downregulated approximately sixfold in the hypothalamus and fivefold in the telencephalon. Intraperitoneal injection of isotocin (1 microg/g) increased plasma estradiol, providing a potential link between changes in isotocin gene expression and decreased circulating estrogen in fluoxetine-injected fish. Our results reveal targets of serotonergic modulation in the neuroendocrine brain and indicate that fluoxetine has the potential to affect sex hormones and modulate genes involved in reproductive function and behavior in the brain of female goldfish. We discuss these findings in the context of endocrine disruption because fluoxetine has been detected in the environment.

Novel candidate genes identified in the brain during nociception in common carp (Cyprinus carpio) and rainbow trout (Oncorhynchus mykiss).

Recent studies have demonstrated that teleost fish possess nociceptors that detect potentially painful stimuli and that the physiological properties of these fibres are markedly similar to those found in mammals. This finding led to suggestions of possible pain perception in fish, contrary to the view that the sensory response in these animals is limited to the spinal cord and hindbrain and as such is reflexive. Therefore, the aim of this study was to determine if the brain is active at the molecular level by using a microarray analysis of gene expression in the forebrain, midbrain and hindbrain of two fish species. A comparison between the two species at different time points showed that many genes were differentially regulated in response to a noxious stimulus compared with controls. A number of genes that are involved in mammalian nociception, such as brain-derived neurotrophic factor (BDNF) and the cannabinoid CB1 receptor were regulated in the fish brain after a nociceptive event. Novel candidates that showed significant regulation in both species were also identified. In particular, the Van Gogh-like 2 gene, was regulated in both carp and trout and should be pursued to establish its precise role in nociception.

Beyond the lipid hypothesis: mechanisms underlying phenotypic plasticity in inducible cold tolerance.

The physiological adjustment of organisms in response to temperature variation is a crucial part of coping with environmental stress. An important component of the cold response is the increase in membrane lipid unsaturation, and this has been linked to an enhanced resistance to the debilitating or lethal effects of cold. Underpinning the lipid response is the upregulation of fatty acid desaturases (des), particularly those introducing double bonds at the 9-10 position of saturated fatty acids. For plants and microbes there is good genetic evidence that regulation of des genes, and the consequent changes in lipid saturation, are causally linked to generation of a cold-tolerant phenotype. In animals, however, supporting evidence is almost entirely limited to correlations of saturation with cold conditions. We describe our recent attempts to provide a direct test of this relationship by genetic manipulation of the nematode Caenorhabditis elegans. We show that this species displays a strong cold tolerant phenotype induced by prior conditioning to cold, and that this is directly linked to upregulated des activity. However, whilst genetic disruption of des activity and lipid unsaturation significantly reduced cold tolerance, animals retained a substantial component of their stress tolerant phenotype produced by cold conditioning. This indicates that mechanisms other than lipid unsaturation play an important role in cold adaptation.

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.

Encoded microcarriers for high-throughput multiplexed detection.

Since the decoding of the human genome, the quest to obtain more and more molecular information from smaller and smaller samples is intensifying. Today the burden of this challenge is being borne by planar arrays, but the quality of the data provided by this approach is limited by variations in performance between different arrays. Suspension arrays of encoded microspheres provide higher quality data, but the amount of molecular information that can be acquired with them is limited by the number of codes that can be distinguished in the same sample. New methods of preparing encoded particles promise to alleviate this problem, but in the face of a growing number of new technologies it is sometimes difficult to decide which, if any, will succeed. Herein we appraise these new forms of encoded particle critically, and ask if they can deliver the necessary multiplexing power and whether they will perform well in multiplexed assays.

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.

Seasonally hibernating phenotype assessed through transcript screening.

Hibernation is a seasonally entrained and profound phenotypic transition to conserve energy in winter. It involves significant biochemical reprogramming, although our understanding of the underpinning molecular events is fragmentary and selective. We have conducted a large-scale gene expression screen of the golden-mantled ground squirrel, Spermophilus lateralis, to identify transcriptional responses associated specifically with the summer-winter transition and the torpid-arousal transition in winter. We used 112 cDNA microarrays comprising 12,288 probes that cover at least 5,109 genes. In liver, the profiles of torpid and active states in the winter were almost identical, although we identified 102 cDNAs that were differentially expressed between winter and summer, 90% of which were downregulated in the winter states. By contrast, in cardiac tissue, 59 and 115 cDNAs were elevated in interbout arousal and torpor, respectively, relative to the summer active condition, but only 7 were common to both winter states, and during arousal none was downregulated. In brain, 78 cDNAs were found to change in winter, 44 of which were upregulated. Thus transcriptional changes associated with hibernation are qualitatively modest and, since these changes are generally less than twofold, also quantitatively modest. Unbiased Gene Ontology profiling of the transcripts suggests a winter switch to beta-oxidation of lipids in liver and heart, a reduction in metabolism of toxic compounds and the urea cycle in liver, and downregulated electron transport in the brain. We identified just one strongly winter-induced transcript common to all tissues, namely an RNA-binding protein, RBM3. This analysis clearly differentiates responses of the principal tissues, identifies a large number of new genes undergoing regulation, and broadens our understanding of affected cellular processes that, in part, account for the winter-adaptive hibernating phenotype.

The use of transcriptomics to address questions in behaviour: production of a suppression subtractive hybridisation library from dominance hierarchies of rainbow trout.

Microarrays, or gene chips, are transforming the way that gene expression is measured by allowing us to determine the expression of thousands of genes from a sample. This gives immense power to examine gene expression on a global scale within individual animals and between animals. The scope for analysing complex animal functions at the molecular level is within our grasp. Relatively few studies have examined complex behaviours and correlated them with gene expression in the central nervous system. Here, we review the use of microarray technology in the dissection of behaviour and focus specifically on dominance status. A cDNA library using suppression subtraction hybridisation on rainbow trout Oncorhynchus mykiss of differing status has been produced to enrich the cDNA library for genes that are differentially expressed between individuals of different dominance status. A preliminary analysis demonstrated that there were 1,165 genes that differed between fish of different dominance status. Therefore, there is the potential of correlating gene expression profile with rank position within dominance hierarchies, thus identifying targets for candidate gene approaches.

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.