Population genomics of wild and laboratory zebrafish (Danio rerio)

Comparative study of the growth, stress status and reproductive capabilities of four wild-type zebrafish (Danio rerio) lines

BACKGROUND

Zebrafish are widely used in various research fields and to fulfil the diverse research needs, numerous zebrafish lines are available, each with a unique domestication background, potentially resulting in intraspecies differences in specific biological functions. Few studies have compared multiple zebrafish lines under identical conditions to investigate both inter- and intra-line variability related to different functions. However, such variability could pose a challenge for the reproducibility of results in studies utilising zebrafish, particularly when the line used is not clearly specified. This study assessed growth, stress status (cortisol, serotonin) and reproductive capabilities (maturity, fecundity, fertilisation rate, sperm quality) of four commonly used wild-type zebrafish lines (AB, SJD, TU, WIK) using standardized protocols.

RESULTS

The stress markers levels were found to be similar across the lines, indicating that the endocrine stress status is robust to diverse domestication histories. Variations were observed in the growth and reproductive parameters. The lines exhibited differences in the timing of puberty (86 dpf for AB and SJD lines vs. 107 dpf for the WIK line) despite achieving similar sizes, suggesting that there are line-specific variations in the induction of maturation. Additionally, the AB line demonstrated higher sperm quality than did the other lines and higher fecundity and fertilization rates than did the SJD line. The AB line also exhibiting a smaller adult size but a heavier brain relative to its body weight.

CONCLUSION

These findings emphasize the importance of line selection for zebrafish research, indicating that researchers should consider line-specific traits to ensure the biological relevance and reproducibility of the results.

Individual differences in the boldness of female zebrafish are associated with alterations in serotonin function

One of the most prevalent axes of behavioral variation in both humans and animals is risk taking, where individuals that are more willing to take risk are characterized as bold while those that are more reserved are regarded as shy. Brain monoamines (i.e. serotonin, dopamine and noradrenaline) have been found to play a role in a variety of behaviors related to risk taking. Using zebrafish, we investigated whether there was a relationship between monoamine function and boldness behavior during exploration of a novel tank. We found a correlation between serotonin metabolism (5-HIAA:5-HT ratio) and boldness during the initial exposure to the tank in female animals. The DOPAC:DA ratio correlated with boldness behavior on the third day in male fish. There was no relationship between boldness and noradrenaline. To probe differences in serotonergic function in bold and shy fish, we administered a selective serotonin reuptake inhibitor, escitalopram, and assessed exploratory behavior. We found that escitalopram had opposing effects on thigmotaxis in bold and shy female animals: the drug caused bold fish to spend more time near the center of the tank and shy fish spent more time near the periphery. Taken together, our findings indicate that variation in serotonergic function has sex-specific contributions to individual differences in risk-taking behavior.

Copy number variation and population-specific immune genes in the model vertebrate zebrafish

Copy number variation in large gene families is well characterized for plant resistance genes, but similar studies are rare in animals. The zebrafish (Danio rerio) has hundreds of NLR immune genes, making this species ideal for studying this phenomenon. By sequencing 93 zebrafish from multiple wild and laboratory populations, we identified a total of 1513 NLRs, many more than the previously known 400. Approximately half of those are present in all wild populations, but only 4% were found in 80% or more of the individual fish. Wild fish have up to two times as many NLRs per individual and up to four times as many NLRs per population than laboratory strains. In contrast to the massive variability of gene copies, nucleotide diversity in zebrafish NLR genes is very low: around half of the copies are monomorphic and the remaining ones have very few polymorphisms, likely a signature of purifying selection.

Key HPI axis receptors facilitate light adaptive behavior in larval zebrafish

The vertebrate stress response (SR) is mediated by the hypothalamic-pituitary-adrenal (HPA) axis and contributes to generating context appropriate physiological and behavioral changes. Although the HPA axis plays vital roles both in stressful and basal conditions, research has focused on the response under stress. To understand broader roles of the HPA axis in a changing environment, we characterized an adaptive behavior of larval zebrafish during ambient illumination changes. Genetic abrogation of glucocorticoid receptor (nr3c1) decreased basal locomotor activity in light and darkness. Some key HPI axis receptors (mc2r [ACTH receptor], nr3c1), but not nr3c2 (mineralocorticoid receptor), were required to adapt to light more efficiently but became dispensable when longer illumination was provided. Such light adaptation was more efficient in dimmer light. Our findings show that the HPI axis contributes to the SR, facilitating the phasic response and maintaining an adapted basal state, and that certain adaptations occur without HPI axis activity.

River zebrafish combine behavioral plasticity and generalized morphology with specialized sensory and metabolic physiology to survive in a challenging environment

Phenotypes that allow animals to detect, weather, and predict changes efficiently are essential for survival in fluctuating environments. Some phenotypes may remain specialized to suit an environment perfectly, while others become more plastic or generalized, shifting flexibly to match current context or adopting a form that can utilize a wide range of contexts. Here, we tested the differences in behavior, morphology, sensory and metabolic physiology between wild zebrafish (Danio rerio) in highly variable fast-flowing rivers and still-water sites. We found that river zebrafish moved at higher velocities than did still-water fish, had lower oxygen demands, and responded less vigorously to small changes in flow rate, as we might expect for fish that are well-suited to high-flow environments. River zebrafish also had less streamlined bodies and were more behaviorally plastic than were still-water zebrafish, both features that may make them better-suited to a transitional lifestyle. Our results suggest that zebrafish use distinct sensory mechanisms and metabolic physiology to reduce energetic costs of living in fast-flowing water while relying on morphology and behavior to create flexible solutions to a challenging habitat. Insights on animals' reliance on traits with different outcomes provide a framework to better understand their survival in future environmental fluctuations.

The canonical HPA axis facilitates and maintains light adaptive behavior

The vertebrate stress response (SR) is mediated by the hypothalamic-pituitary-adrenal (HPA) axis and contributes to generating context appropriate physiological and behavioral changes. Although the HPA axis plays vital roles both in stressful and basal conditions, research has focused on the response under stress. To understand broader roles of the HPA axis in a changing environment, we characterized an adaptive behavior of larval zebrafish during ambient illumination changes. The glucocorticoid receptor (nr3c1) was necessary to maintain basal locomotor activity in light and darkness. The HPA axis was required to adapt to light more efficiently but became dispensable when longer illumination was provided. Light adaptation was more efficient in dimmer light and did not require the mineralocorticoid receptor (nr3c2). Our findings show that the HPA axis contributes to the SR at various stages, facilitating the phasic response and maintaining an adapted basal state, and that certain adaptations occur without HPA axis activity.

Effective practices for thermal tolerance polygon experiments using mottled catfish Corydoras paleatus

Temperature is an important environmental factor that affects how organisms allocate metabolic resources to physiological processes. Laboratory experiments that determine absolute thermal limits for representative species are important for understanding how fishes are affected by climate change. Critical Thermal Methodology (CTM) and Chronic Lethal Methodology (CLM) experiments were utilized to construct a complete thermal tolerance polygon for the South American fish species, Mottled catfish (Corydoras paleatus). Mottled catfish showed Chronic Lethal Maxima (CLMax) of 34.9 ± 0.52 °C and Chronic Lethal Minima (CLMin) of 3.8 ± 0.08 °C. Fish were chronically acclimated (∼2 weeks) to 6 temperatures ranging from 7.2 ± 0.05 °C →32.2 ± 0.16 °C (7 °C, 12 °C, 17 °C, 22 °C, 27 °C, and 32 °C), and CTM used to estimate upper and lower acute temperature tolerance. Linear regressions of Critical Thermal Maxima (CTMax) and Minima (CTMin) data with each acclimation temperature were used along with CLMax and CLMin to create a complete thermal tolerance polygon. The highest CTMax was 38.4 ± 0.60 °C for fish acclimated to 32.2 ± 0.16 °C, and the lowest CTMin was 3.36 ± 1.84 °C for fish acclimated to 7.2 ± 0.05 °C. Mottled catfish have a polygon measuring 785.7°C2, and the slope of the linear regressions showed the species gained 0.55 °C and 0.32 °C of upper and lower tolerance per degree of acclimation temperature, respectively. We compared slopes of CTMax or CTMin regression lines to each other using a set of comparisons between 3, 4, 5, or 6 acclimation temperatures. Our data demonstrated that 3 acclimation temperatures were as sufficient as 4 → 6 to pair with estimates of chronic upper and lower thermal limits for accurately determining a complete thermal tolerance polygon. Construction of this species' complete thermal tolerance polygon provides a template for other researchers. The following is sufficient to generate a complete thermal tolerance polygon: Three chronic acclimation temperatures that are spread somewhat evenly across a species' thermal range, include an estimation of CLMax and CLMin, and are followed by CTMax and CTMin measurements.

Correlations begin at home: drivers of co-occurrence patterns in personality and cognitive ability in wild populations of zebrafish

Aquatic habitats are extremely dynamic, with constantly changing ecological factors, which has now been exacerbated due to human-induced rapid environmental change. In such variable environments, it becomes essential to understand how personality and cognition in organisms affect the adaptability of individuals to different habitat conditions. To test this, we studied how personality-related traits as well as cognitive ability differ between populations of wild-caught zebrafish (Danio rerio) from habitats that differed in various environmental factors. We measured emergence into a novel environment as an indicator of boldness, and performance in a spatial task inferred from feeding latencies in a maze over repeated trials to assess learning and memory, as an indicator of cognitive ability. We found that personality affects cognition and although bolder fish are better learners, they show poorer retention of memory across populations. Although personality and cognitive ability varied between habitats, the patterns of their correlations remained similar within each population. However, the individual traits (such as sex and size) that were drivers of personality and cognition differed between the habitats, suggesting that not only do behavioral traits vary between populations, but also the factors that are important in determining them. Personality and cognitive ability and the correlations between these traits determine how well an organism performs in its habitat, as well as how likely it is to find new habitats and adapt to them. Studying these across wild zebrafish populations helps predict performance efficiencies among individuals and also explains how fish adapt to extremely dynamic environments that can lead to variation in behavioral traits and correlations between them. This study not only sheds light on the drivers of interindividual variation and co-occurrence patterns of personality and cognition, but also individual and population factors that might have an effect on them.

The phylogenetic position of zebrafish (Danio rerio) from south african pet shops

BACKGROUND

Zebrafish (Danio rerio), a small freshwater fish that originates from India, Bangladesh, Nepal, Bhutan and northern Myanmar, have been widely used as a model organism for studies of developmental biology and genetics. The current study aimed to determine the origin of South African pet shop stock that are currently being used to establish a laboratory population founded from diverse sources available locally.

METHODS AND RESULTS

Zebrafish DNA was extracted from 65 specimens housed at the University of the Free State (UFS) Department of Genetics. For phylogenetic analysis, cytb sequences were generated from all samples. A further 178 sequences were downloaded from the GenBank database, including sequences of an outgroup species (Danio kyathit). Five microsatellite markers were used to further assess the genetic diversity of the UFS zebrafish specimens. A maximum likelihood analysis was performed for the cytb data. Results of the phylogenetic analyses divided the sequences into three major genetic groups, which was congruent with a previous study on laboratory zebrafish provenance. The SA pet shop fish grouped with the lines from the northern and north-eastern regions of India. High levels of microsatellite genetic diversity were observed for the pet shop sourced population, correlating to what has previously been observed in zebrafish.

CONCLUSION

These results can be used to guide the future development of laboratory lines suited to the needs at the UFS.

Behavioral genetics of alcohol's effects in three zebrafish (Danio rerio) populations

Alcohol abuse is one of the most dangerous and serious problems for patients and society. Interpopulation studies are important in understanding how genetic background contributes to the effects of alcohol. In this study, we applied a chronic alcohol exposure protocol in three zebrafish populations (Danio rerio; both sexes; AB, TU, and outbred fish - OB). We analyzed the behavioral responses and mRNA expression involved in neurotransmitter metabolism - th1, tph1, ache, ada1, gaba1, gad1b, and bdnf. Locomotion patterns were similar between populations (increased speed after acute alcohol and unaltered locomotion after chronic and withdrawal treatments). All populations exhibited increased expression of genes associated with locomotion (th1, gad1b, and gaba1) after acute alcohol exposure. Anxiety-like responses increased in AB and TU fish during withdrawal and decreased in AB fish after acute alcohol exposure. Genes related to anxiety-like behavior (tph1 and ada1) were overexpressed in AB and TU fish after acute and withdrawal treatments, while OB fish exhibited unaltered responses. Bdnf levels decreased during withdrawal in AB and OB fish, while TU showed upregulated levels in both chronic and withdrawal treatments. Our results suggest that zebrafish populations respond differently to alcohol exposure, which may contribute to understanding the mechanisms underlying alcohol use and dependence. Moreover, we found that a more diverse genetic background (OB) was related to higher variability in behavioral and mRNA expression, demonstrating that inbred populations (AB and TU) may be useful tools in identifying alcohol use and abuse mechanisms.

Genomic basis of fishing-associated selection varies with population density

Fisheries-associated selection is recognized as one of the strongest potential human drivers of contemporary evolution in natural populations. The results of this study show that while simulated commercial fishing techniques consistently remove fish with traits associated with growth, metabolism, and social behavior, the specific genes under fishing selection differ depending on the density of the targeted population. This finding suggests that different fish populations of varying sizes will respond differently to fishing selection at the genetic level. Furthermore, as a population is fished over time, the genes under selection may change as the population diminishes. This could have repercussions on population resilience. This study highlights the importance of selection but also environmental and density effects on harvested fish populations.

Fisheries induce one of the strongest anthropogenic selective pressures on natural populations, but the genetic effects of fishing remain unclear. Crucially, we lack knowledge of how capture-associated selection and its interaction with reductions in population density caused by fishing can potentially shift which genes are under selection. Using experimental fish reared at two densities and repeatedly harvested by simulated trawling, we show consistent phenotypic selection on growth, metabolism, and social behavior regardless of density. However, the specific genes under selection—mainly related to brain function and neurogenesis—varied with the population density. This interaction between direct fishing selection and density could fundamentally alter the genomic responses to harvest. The evolutionary consequences of fishing are therefore likely context dependent, possibly varying as exploited populations decline. These results highlight the need to consider environmental factors when predicting effects of human-induced selection and evolution.

Density influences the heritability and genetic correlations of fish behaviour under trawling-associated selection

Fishing-associated selection is one of the most important human-induced evolutionary pressures for natural populations. However, it is unclear whether fishing leads to heritable phenotypic changes in the targeted populations, as the heritability and genetic correlations of traits potentially under selection have received little attention. In addition, phenotypic changes could arise from fishing-associated environmental effects, such as reductions in population density. Using fish reared at baseline and reduced group density and repeatedly harvested by simulated trawling, we show that trawling can induce direct selection on fish social behaviour. As sociability has significant heritability and is also genetically correlated with activity and exploration, trawling has the potential to induce both direct selection and indirect selection on a variety of fish behaviours, potentially leading to evolution over time. However, while trawling selection was consistent between density conditions, the heritability and genetic correlations of behaviours changed according to the population density. Fishing-associated environmental effects can thus modify the evolutionary potential of fish behaviour, revealing the need to use a more integrative approach to address the evolutionary consequences of fishing.

Microinjection quality control in zebrafish model for genetic manipulations

Microinjection technique is one of the essential methodologies that are used widely in zebrafish research. Microinjection is utilized to perform genetic manipulations within the developing zebrafish model. Further, this technique is used to study a wide range of genetic diseases and gene of interest role in early developmental processes. Thus, quality control for microinjection is an essential factor to ensure experimental reproducibility and consistency. In this technical note, in vitro transcribed synthetic mRNA encoding green fluorescence protein (eGFP), and red fluorescent protein (m-cherry) as well as fluorescein and rhodamine fluorescent dyes were injected into a single-cell zebrafish embryo for volume quality control. Given the importance of having quality control system and methodology to yield similar genetic manipulation within the zebrafish embryo:•

We aimed to establish the unified delivery of injected material into zebrafish one cell stage embryo.

•

We aimed to establish consistency of the injected volume into mineral oil droplets that will serve as a quality control parameter to conforms a quality control practice to ensure the reproducibility of the microinjection technique.

•

The calibration of microinjection droplet size resulted in the visualization of fluorescent protein and dyes in the zebrafish embryo with precise volumes of delivered materials under the control of needle opening, injection pressure and time.

Wild Zebrafish Sentinels: Biological Monitoring of Site Differences Using Behavior and Morphology

Environmental change poses a devastating risk to human and environmental health. Rapid assessment of water conditions is necessary for monitoring, evaluating, and addressing this global health danger. Sentinels or biological monitors can be deployed in the field using minimal resources to detect water quality changes in real time, quickly and cheaply. Zebrafish (Danio rerio) are ideal sentinels for detecting environmental changes due to their biomedical tool kit, widespread geographic distribution, and well-characterized phenotypic responses to environmental disturbances. Here, we demonstrate the utility of zebrafish sentinels by characterizing phenotypic differences in wild zebrafish between two field sites in India. Site 1 was a rural environment with flowing water, low-hypoxic conditions, minimal human-made debris, and high iron and lead concentrations. Site 2 was an urban environment with still water, hypoxic conditions, plastic pollution, and high arsenic, iron, and chromium concentrations. We found that zebrafish from Site 2 were smaller, more cohesive, and less active than Site 1 fish. We also found sexually dimorphic body shapes within the Site 2, but not the Site 1, population. Advancing zebrafish sentinel research and development will enable rapid detection, evaluation, and response to emerging global health threats.

Thermal tolerance and routine oxygen consumption of convict cichlid, Archocentrus nigrofasciatus, acclimated to constant temperatures (20 °C and 30 °C) and a daily temperature cycle (20 °C → 30 °C)

Organismal temperature tolerance and metabolic responses are correlated to recent thermal history, but responses to thermal variability are less frequently assessed. There is great interest in whether organisms that experience greater thermal variability can gain metabolic or tolerance advantages through phenotypic plasticity. We compared thermal tolerance and routine aerobic metabolism of Convict cichlid acclimated for 2 weeks to constant 20 °C, constant 30 °C, or a daily cycle of 20 → 30 °C (1.7 °C/h). Acute routine mass-specific oxygen consumption ([Formula: see text]O₂) and critical thermal maxima/minima (CTMax/CTMin) were compared between groups, with cycle-acclimated fish sampled from the daily minimum (20 °C, 0900 h) and maximum (30 °C, 1600 h). Cycle-acclimated fish demonstrated statistically similar CTMax at the daily minimum and maximum (39.0 °C, 38.6 °C) but distinct CTMin values, with CTMin 2.4 °C higher for fish sampled from the daily 30 °C maximum (14.8 °C) compared to the daily 20 °C minimum (12.4 °C). Measured acutely at 30 °C, [Formula: see text]O₂ decreased with increasing acclimation temperature; 20 °C acclimated fish had an 85% higher average [Formula: see text]O₂ than 30 °C acclimated fish. Similarly, acute [Formula: see text]O₂ at 20 °C was 139% higher in 20 °C acclimated fish compared to 30 °C acclimated fish. Chronic [Formula: see text]O₂ was measured in separate fish continually across the 20 → 30 °C daily cycle for all 3 acclimation groups. Chronic [Formula: see text]O₂ responses were very similar between groups between average individual hourly values, as temperatures increased or decreased (1.7 °C/h). Acute [Formula: see text]O₂ and thermal tolerance responses highlight "classic" trends, but dynamic, chronic trials suggest acclimation history has little effect on the relative change in oxygen consumption during a thermal cycle. Our results strongly suggest that the minimum and maximum temperatures experienced more strongly influence fish physiology, rather than the thermal cycle itself. This research highlights the importance of collecting data in both cycling and static (constant) thermal conditions, and further research should seek to understand whether ectotherm metabolism does respond uniquely to fluctuating temperatures.

The Laboratory Domestication of Zebrafish: From Diverse Populations to Inbred Substrains

We know from human genetic studies that practically all aspects of biology are strongly influenced by the genetic background, as reflected in the advent of “personalized medicine.” Yet, with few exceptions, this is not taken into account when using laboratory populations as animal model systems for research in these fields. Laboratory strains of zebrafish (Danio rerio) are widely used for research in vertebrate developmental biology, behavior, and physiology, for modeling diseases, and for testing pharmaceutic compounds in vivo. However, all of these strains are derived from artificial bottleneck events and therefore are likely to represent only a fraction of the genetic diversity present within the species. Here, we use restriction site-associated DNA sequencing to genetically characterize wild populations of zebrafish from India, Nepal, and Bangladesh, and to compare them to previously published data on four common laboratory strains. We measured nucleotide diversity, heterozygosity, and allele frequency spectra, and find that wild zebrafish are much more diverse than laboratory strains. Further, in wild zebrafish, there is a clear signal of GC-biased gene conversion that is missing in laboratory strains. We also find that zebrafish populations in Nepal and Bangladesh are most distinct from all other strains studied, making them an attractive subject for future studies of zebrafish population genetics and molecular ecology. Finally, isolates of the same strains kept in different laboratories show a pattern of ongoing differentiation into genetically distinct substrains. Together, our findings broaden the basis for future genetic, physiological, pharmaceutic, and evolutionary studies in Danio rerio.

Low potential for evolutionary rescue from climate change in a tropical fish

Climate change is increasing global temperatures and intensifying the frequency and severity of extreme heat waves. How organisms will cope with these changes depends on their inherent thermal tolerance, acclimation capacity, and ability for evolutionary adaptation. Yet, the potential for adaptation of upper thermal tolerance in vertebrates is largely unknown. We artificially selected offspring from wild-caught zebrafish (Danio rerio) to increase (Up-selected) or decrease (Down-selected) upper thermal tolerance over six generations. Selection to increase upper thermal tolerance was also performed on warm-acclimated fish to test whether plasticity in the form of inducible warm tolerance also evolved. Upper thermal tolerance responded to selection in the predicted directions. However, compared to the control lines, the response was stronger in the Down-selected than in the Up-selected lines in which evolution toward higher upper thermal tolerance was slow (0.04 ± 0.008 °C per generation). Furthermore, the scope for plasticity resulting from warm acclimation decreased in the Up-selected lines. These results suggest the existence of a hard limit in upper thermal tolerance. Considering the rate at which global temperatures are increasing, the observed rates of adaptation and the possible hard limit in upper thermal tolerance suggest a low potential for evolutionary rescue in tropical fish living at the edge of their thermal limits.

A fish is not a mouse: understanding differences in background genetics is critical for reproducibility

Poorly controlled background genetics in animal models contributes to the lack of reproducibility that is increasingly recognized in biomedical research. The laboratory zebrafish, Danio rerio, has been an important model organism for decades in many research areas, yet inbred strains and traditionally managed outbred stocks are not available for this species. Sometimes incorrectly referred to as 'inbred strains' or 'strains', zebrafish wild-type lines possess background genetics that are often not well characterized, and breeding practices for these lines have not been consistent over time or among institutions. In this Perspective, we trace key milestones in the history of one of the most widely used genetic backgrounds, the AB line, to illustrate the dynamic complexity within an example background that is largely invisible when reading the scientific literature. Failure to adequately control for genetic background compromises the validity of experimental outcomes. We therefore propose that authors provide as much specific detail about the origin and genetic makeup of zebrafish lines as is reasonable and possible, and that the terms used to describe background genetics be applied in a way that is consistent with other fish and mammalian model organisms. We strongly encourage the adoption of genetic monitoring for the characterization of existing zebrafish lines, to help detect genetic contamination in breeding colonies and to verify the level of genetic heterogeneity in breeding colonies over time. Careful attention to background genetics will improve transparency and reproducibility, therefore improving the utility of the zebrafish as a model organism.

Shedding new light on early sex determination in zebrafish

In contrast to established zebrafish gene annotations, the question of sex determination has still not been conclusively clarified for developing zebrafish, Danio rerio, larvae, 28 dpf or earlier. Recent studies indicate polygenic sex determination (PSD), with the genes being distributed throughout the genome. Early genetic markers of sex in zebrafish help unravel co-founding sex-related differences to apply to human health and environmental toxicity studies. A qPCR-based method was developed for six genes: cytochrome P450, family 17, subfamily A, polypeptide 1 (cyp17a1); cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a); cytochrome P450, family 19, subfamily A, polypeptides 1b (cyp19a1b); vitellogenin 1 (vtg1); nuclear receptor subfamily 0, group B, member 1 (nr0b1), sry (sex-determining region Y)-box 9b (sox9b) and actin, beta 1 (actb1), the reference gene. Sry-box 9a (Sox9a), insulin-like growth factor 3 (igf3) and double sex and mab-3 related transcription factor 1 (dmrt1), which are also known to be associated with sex determination, were used in gene expression tests. Additionally, Next-Generation-Sequencing (NGS) sequenced the genome of two adult female and male and two juveniles. PCR analysis of adult zebrafish revealed sex-specific expression of cyp17a1, cyp19a1a, vtg1, igf3 and dmrt1, the first four strongly expressed in female zebrafish and the last one highly expressed in male conspecifics. From NGS, nine female and four male-fated genes were selected as novel for assessing zebrafish sex, 28 dpf. Differences in transcriptomes allowed allocation of sex-specific genes also expressed in juvenile zebrafish.

Development and Validation of an Open Access SNP Array for Nile Tilapia (Oreochromis niloticus)

Tilapia are among the most important farmed fish species worldwide, and are fundamental for the food security of many developing countries. Several genetically improved Nile tilapia (Oreochromis niloticus) strains exist, such as the iconic Genetically Improved Farmed Tilapia (GIFT), and breeding programs typically follow classical pedigree-based selection. The use of genome-wide single-nucleotide polymorphism (SNP) data can enable an understanding of the genetic architecture of economically important traits and the acceleration of genetic gain via genomic selection. Due to the global importance and diversity of Nile tilapia, an open access SNP array would be beneficial for aquaculture research and production. In the current study, a ∼65K SNP array was designed based on SNPs discovered from whole-genome sequence data from a GIFT breeding nucleus population and the overlap with SNP datasets from wild fish populations and several other farmed Nile tilapia strains. The SNP array was applied to clearly distinguish between different tilapia populations across Asia and Africa, with at least ∼30,000 SNPs segregating in each of the diverse population samples tested. It is anticipated that this SNP array will be an enabling tool for population genetics and tilapia breeding research, facilitating consistency and comparison of results across studies.

Machine-learning-based detection of adaptive divergence of the stream mayfly Ephemera strigata populations

Adaptive divergence is a key mechanism shaping the genetic variation of natural populations. A central question linking ecology with evolutionary biology is how spatial environmental heterogeneity can lead to adaptive divergence among local populations within a species. In this study, using a genome scan approach to detect candidate loci under selection, we examined adaptive divergence of the stream mayfly Ephemera strigata in the Natori River Basin in northeastern Japan. We applied a new machine-learning method (i.e., random forest) besides traditional distance-based redundancy analysis (dbRDA) to examine relationships between environmental factors and adaptive divergence at non-neutral loci. Spatial autocorrelation analysis based on neutral loci was employed to examine the dispersal ability of this species. We conclude the following: (a) E. strigata show altitudinal adaptive divergence among the populations in the Natori River Basin; (b) random forest showed higher resolution for detecting adaptive divergence than traditional statistical analysis; and (c) separating all markers into neutral and non-neutral loci could provide full insight into parameters such as genetic diversity, local adaptation, and dispersal ability.

Collective Behavior in Wild Zebrafish

Anthropogenic change is expected to alter environments at alarming rates. To predict the impact of modified environments on social behavior, we must study the relationship between environmental features and collective behavior in a genetically tractable model, zebrafish (Danio rerio). Here, we conducted a field study to examine the relationship between salient environmental features and collective behavior in four populations of zebrafish. We found zebrafish in flowing water formed volatile groups, whereas those in still water had more consistent membership and leadership. Groups in fast-flowing water were large (up to 2000 fish) and tightly knit with short nearest neighbor distances, whereas group sizes were smaller (11 fish/group) with more space between individual fish in still and slow-flowing water. These observations point to a possible profound role of water flow in influencing collective behavior in wild zebrafish.

Genome-wide association analyses based on whole-genome sequencing of Protosalanx hyalocranius provide insights into sex determination of Salangid fishes

Identification of sex determination system and sex-determining genes have important implications in conservation, ecology and evolution. However, much remains to be discovered about the evolution of different sexual determination systems in teleost fishes, of which the mechanisms of sex determination are remarkably variable. In the present study, the whole genomes of 20 males and 20 females of a Salangid fish, Protosalanx hyalocranius, were sequenced and genome wide association analyses were conducted to uncover its sex determination system and putative sex-determining genes. A total of 150 SNPs were significantly associated with sex, which showed high differentiation between sexes (F_ST ranged from 0.245 to 0.556). Of the 150 sex-associated SNPs, 76 SNPs displayed sex specificity with even coverage of depth and were female heterogametic, which suggested a ZZ/ZW sex determination system. Interestingly, one scaffold containing sex-specific SNPs displayed synteny to the sex chromosome of medaka. Annotations of sex-associated loci suggested that both transcriptional regulators (e.g., FOX genes) and secreted hormones and their receptors might be involved in the sex determination/differentiation of P. hyalocranius. More strikingly, we found a nonsense mutation in one copy of GALNT homology gene of all females, which suggested that "Z dosage" effect might play a vital role in the processes of sex determination/differentiation. These sex-specific loci could be a valuable resource for further research on sex determination of Salangid fishes and the results could contribute to the understanding of sex determination mechanisms and the evolution of sex chromosome in teleost fishes.

A Comparison of Isogenic Homozygous Clone and Wildtype Zebrafish (Danio rerio): Survival and Developmental Responses to Low pH Conditions

The value of bioassays as analytical methods for assessing the potency of particular stressors on live animal models depends on the precision of their results, which are greatly influenced by the choice of test subjects. The genetic makeup of experimental subjects varies, and, as such, so will their responses to the test environment. Genetic diversity of test populations may contribute to statistical variability; therefore, the use of genetically similar subjects may enhance the utility of bioassays. This study addresses the efficacy of using isogenic homozygous zebrafish (Danio rerio) as subjects for bioassays. Stress responses (acidic conditions) were compared during early development for gynogenetically produced isogenic homozygous line of zebrafish (C32) and wildtype (WT) zebrafish. Experiments evaluated early life stage milestones after exposure to low pH in water of a different electrolyte composition. Because the isogenic homozygous clonal (IHC) fish possessed far less genetic variability than the WT fish tested, it was predicted that the IHC fish would exhibit less variability in their response to stress. Although we found no significant differences in the variability between the responses of the IHC and WT fish, pH and water hardness level had a differential effect on the two groups. Simple strain differences may be the probable cause of the response differences to environmental stress. Factors that may affect stress response, such as heterogeneity, co-adapted gene complexes, and domestication, are discussed. Our findings and review of recent zebrafish literature stress the need for researchers to carefully consider breeding histories and trait characteristics for each potential test subject to maximize the sensitivity of the assay.

Zebrafish: Housing and husbandry recommendations

This article provides recommendations for the care of laboratory zebrafish (Danio rerio) as part of the further implementation of Annex A to the European Convention on the protection of vertebrate animals used for experimental and other scientific purposes, EU Commission Recommendation 2007/526/EC and the fulfilment of Article 33 of EU Directive 2010/63, both concerning the housing and care of experimental animals. The recommendations provide guidance on best practices and ranges of husbandry parameters within which zebrafish welfare, as well as reproducibility of experimental procedures, are assured.

Husbandry procedures found today in zebrafish facilities are numerous. While the vast majority of these practices are perfectly acceptable in terms of zebrafish physiology and welfare, the reproducibility of experimental results could be improved by further standardisation of husbandry procedures and exchange of husbandry information between laboratories. Standardisation protocols providing ranges of husbandry parameters are likely to be more successful and appropriate than the implementation of a set of fixed guidance values neglecting the empirically successful daily routines of many facilities and will better reflect the wide range of environmental parameters that characterise the natural habitats occupied by zebrafish.

A joint working group on zebrafish housing and husbandry recommendations, with members of the European Society for Fish Models in Biology and Medicine (EUFishBioMed) and of the Federation of European Laboratory Animal Science Associations (FELASA) has been given a mandate to provide guidelines based on a FELASA list of parameters, ‘Terms of Reference’.

Building a DNA barcode library for the freshwater fishes of Bangladesh

We sequenced the standard DNA barcode gene fragment in 694 newly collected specimens, representing 243 species level Operational Barcode Units (OBUs) of freshwater fishes from Bangladesh. We produced coi sequences for 149 out of the 237 species already recorded from Bangladesh. Another 83 species sequenced were not previously recorded for the country, and include about 30 undescribed or potentially undescribed species. Several of the taxa that we could not sample represent erroneous records for the country, or sporadic occurrences. Species identifications were classified at confidence levels 1(best) to 3 (worst). We propose the new term Operational Barcode Unit (OBU) to simplify references to would-be DNA barcode sequences and sequence clusters. We found one case where there were two mitochondrial lineages present in the same species, several cases of cryptic species, one case of introgression, one species yielding a pseudogene to standard barcoding primers, and several cases of taxonomic uncertainty and need for taxonomic revision. Large scale national level DNA barcode prospecting in high diversity regions may suffer from lack of taxonomic expertise that cripples the result. Consequently, DNA barcoding should be performed in the context of taxonomic revision, and have a defined, competent end-user.

Are model organisms representative for climate change research? Testing thermal tolerance in wild and laboratory zebrafish populations

Model organisms can be useful for studying climate change impacts, but it is unclear whether domestication to laboratory conditions has altered their thermal tolerance and therefore how representative of wild populations they are. Zebrafish in the wild live in fluctuating thermal environments that potentially reach harmful temperatures. In the laboratory, zebrafish have gone through four decades of domestication and adaptation to stable optimal temperatures with few thermal extremes. If maintaining thermal tolerance is costly or if genetic traits promoting laboratory fitness at optimal temperature differ from genetic traits for high thermal tolerance, the thermal tolerance of laboratory zebrafish could be hypothesized to be lower than that of wild zebrafish. Furthermore, very little is known about the thermal environment of wild zebrafish and how close to their thermal limits they live. Here, we compared the acute upper thermal tolerance (critical thermal maxima; CT_max) of wild zebrafish measured on-site in West Bengal, India, to zebrafish at three laboratory acclimation/domestication levels: wild-caught, F₁ generation wild-caught and domesticated laboratory AB-WT line. We found that in the wild, CT_max increased with increasing site temperature. Yet at the warmest site, zebrafish lived very close to their thermal limit, suggesting that they may currently encounter lethal temperatures. In the laboratory, acclimation temperature appeared to have a stronger effect on CT_max than it did in the wild. The fish in the wild also had a 0.85-1.01°C lower CT_max compared to all laboratory populations. This difference between laboratory-held and wild populations shows that environmental conditions can affect zebrafish's thermal tolerance. However, there was no difference in CT_max between the laboratory-held populations regardless of the domestication duration. This suggests that thermal tolerance is maintained during domestication and highlights that experiments using domesticated laboratory-reared model species can be appropriate for addressing certain questions on thermal tolerance and global warming impacts.

Using zebrafish as a model to study the role of epigenetics in hearing loss

INTRODUCTION

The rapid progress of bioinformatics and high-throughput screening techniques in recent years has led to the identification of many candidate genes and small-molecule drugs that have the potential to make significant contributions to our understanding of the developmental and pathological processes of hearing, but it remains unclear how these genes and regulatory factors are coordinated. Increasing evidence suggests that epigenetic mechanisms are essential for establishing gene expression profiles and likely play an important role in the development of inner ear and in the pathology of hearing-associated diseases. Zebrafish are a valuable and tractable in vivo model organism for monitoring changes in the epigenome and for identifying new epigenetic processes and drug molecules that can influence vertebrate development. Areas covered: In this review, the authors focus on zebrafish as a model to summarize recent findings concerning the roles of epigenetics in the development, regeneration, and protection of hair cells. Expert opinion: Using the zebrafish model in combination with high-throughput screening and genome-editing technologies to investigate the function of epigenetics in hearing is crucial to help us better understand the molecular and genetic mechanisms of auditory development and function. It will also contribute to the development of new strategies to restore hearing loss.

Breaking RAD: an evaluation of the utility of restriction site-associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site-associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome-scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq-based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD-tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best-case modelling scenario, we found that RADseq using six- or eight-base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD-tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.

Body length rather than routine metabolic rate and body condition correlates with activity and risk-taking in juvenile zebrafish Danio rerio

In this study, the following hypotheses were explored using zebrafish Danio rerio: (1) individuals from the same cohort differ consistently in activity and risk-taking and (2) variation in activity and risk-taking is linked to individual differences in metabolic rate, body length and body condition. To examine these hypotheses, juvenile D. rerio were tested for routine metabolic rate and subsequently exposed to an open field test. Strong evidence was found for consistent among-individual differences in activity and risk-taking, which were overall negatively correlated with body length, i.e. larger D. rerio were found to be less active in a potentially dangerous open field and a similar trend was found with respect to a more direct measure of their risk-taking tendency. In contrast, routine metabolic rate and body condition were uncorrelated with both activity and risk-taking of juvenile D. rerio. These findings suggest that body length is associated with risk-related behaviours in juvenile D. rerio for which larger, rather than smaller, individuals may have a higher risk of predation, while the role for routine metabolic rate is relatively limited or non-existent, at least under the conditions of the present study.

Aquaculture, husbandry, and shipping at the Zebrafish International Resource Center

This chapter provides an overview of the Zebrafish International Resource Center (ZIRC) aquaculture systems, husbandry, and live fish shipping methods. The ZIRC has a recirculating water system with mechanical and biological filtration in its main fish facility, and a flow-through water system for its quarantine room. I describe basic husbandry methods for breeding, rearing larvae, feeding and shipping of fish. Whereas the procedures presented here are merely an overview and only the most essential methods are included, they offer starting points to set up, develop, refine, or troubleshoot methods in other fish facilities.

High major histocompatibility complex class I polymorphism despite bottlenecks in wild and domesticated populations of the zebra finch (Taeniopygia guttata)

BACKGROUND

Two subspecies of zebra finch, Taeniopygia guttata castanotis and T. g. guttata are native to Australia and the Lesser Sunda Islands, respectively. The Australian subspecies has been domesticated and is now an important model system for research. Both the Lesser Sundan subspecies and domesticated Australian zebra finches have undergone population bottlenecks in their history, and previous analyses using neutral markers have reported reduced neutral genetic diversity in these populations. Here we characterize patterns of variation in the third exon of the highly variable major histocompatibility complex (MHC) class I α chain. As a benchmark for neutral divergence, we also report the first mitochondrial NADH dehydrogenase 2 (ND2) sequences in this important model system.

RESULTS

Despite natural and human-mediated population bottlenecks, we find that high MHC class I polymorphism persists across all populations. As expected, we find higher levels of nucleotide diversity in the MHC locus relative to neutral loci, and strong evidence of positive selection acting on important residues forming the peptide-binding region (PBR). Clear population differentiation of MHC allele frequencies is also evident, and this may be due to adaptation to new habitats and associated pathogens and/or genetic drift. Whereas the MHC Class I locus shows broad haplotype sharing across populations, ND2 is the first locus surveyed to date to show reciprocal monophyly of the two subspecies.

CONCLUSIONS

Despite genetic bottlenecks and genetic drift, all surveyed zebra finch populations have maintained high MHC Class I diversity. The diversity at the MHC Class I locus in the Lesser Sundan subspecies contrasts sharply with the lack of diversity in previously examined neutral loci, and may thus be a result of selection acting to maintain polymorphism. Given uncertainty in historical population demography, however, it is difficult to rule out neutral processes in maintaining the observed diversity. The surveyed populations also differ in MHC Class I allele frequencies, and future studies are needed to assess whether these changes result in functional immune differences.

Parallelism and Epistasis in Skeletal Evolution Identified through Use of Phylogenomic Mapping Strategies

The identification of genetic mechanisms underlying evolutionary change is critical to our understanding of natural diversity, but is presently limited by the lack of genetic and genomic resources for most species. Here, we present a new comparative genomic approach that can be applied to a broad taxonomic sampling of nonmodel species to investigate the genetic basis of evolutionary change. Using our analysis pipeline, we show that duplication and divergence of fgfr1a is correlated with the reduction of scales within fishes of the genus Phoxinellus. As a parallel genetic mechanism is observed in scale-reduction within independent lineages of cypriniforms, our finding exposes significant developmental constraint guiding morphological evolution. In addition, we identified fixed variation in fgf20a within Phoxinellus and demonstrated that combinatorial loss-of-function of fgfr1a and fgf20a within zebrafish phenocopies the evolved scalation pattern. Together, these findings reveal epistatic interactions between fgfr1a and fgf20a as a developmental mechanism regulating skeletal variation among fishes.

The mind, the lab, and the field: Three kinds of populations in scientific practice

Scientists use models to understand the natural world, and it is important not to conflate model and nature. As an illustration, we distinguish three different kinds of populations in studies of ecology and evolution: theoretical, laboratory, and natural populations, exemplified by the work of R. A. Fisher, Thomas Park, and David Lack, respectively. Biologists are rightly concerned with all three types of populations. We examine the interplay between these different kinds of populations, and their pertinent models, in three examples: the notion of "effective" population size, the work of Thomas Park on Tribolium populations, and model-based clustering algorithms such as Structure. Finally, we discuss ways to move safely between three distinct population types while avoiding confusing models and reality.

A population genomic scan in Chorthippus grasshoppers unveils previously unknown phenotypic divergence

Understanding the genetics of speciation and the processes that drive it is a central goal of evolutionary biology. Grasshoppers of the Chorthippus species group differ strongly in calling song (and corresponding female preferences) but are exceedingly similar in other characteristics such as morphology. Here, we performed a population genomic scan on three Chorthippus species (Chorthippus biguttulus, C. mollis and C. brunneus) to gain insight into the genes and processes involved in divergence and speciation in this group. Using an RNA-seq approach, we examined functional variation between the species by calling SNPs for each of the three species pairs and using FST -based approaches to identify outliers. We found approximately 1% of SNPs in each comparison to be outliers. Between 37% and 40% of these outliers were nonsynonymous SNPs (as opposed to a global level of 17%) indicating that we recovered loci under selection. Among the outliers were several genes that may be involved in song production and hearing as well as genes involved in other traits such as food preferences and metabolism. Differences in food preferences between species were confirmed with a behavioural experiment. This indicates that multiple phenotypic differences implicating multiple evolutionary processes (sexual selection and natural selection) are present between the species.

Characterizing the neurotranscriptomic states in alternative stress coping styles

BACKGROUND

Animals experience stress in many contexts and often successfully cope. Individuals exhibiting the proactive versus reactive stress coping styles display qualitatively different behavioral and neuroendocrine responses to stressors. The predisposition to exhibiting a particular coping style is due to genetic and environmental factors. In this study we explore the neurotranscriptomic and gene network biases that are associated with differences between zebrafish (Danio rerio) lines selected for proactive and reactive coping styles and reared in a common garden environment.

RESULTS

Using RNA-sequencing we quantified the basal transcriptomes from the brains of wild-derived zebrafish lines selectively bred to exhibit the proactive or reactive stress coping style. We identified 1953 genes that differed in baseline gene expression levels. Weighted gene coexpression network analyses identified one gene module associated with line differences. Together with our previous pharmacological experiment, we identified a core set of 62 genes associated with line differences. Gene ontology analyses reveal that many of these core genes are implicated in neurometabolism (e.g. organic acid biosynthetic and fatty acid metabolic processes).

CONCLUSIONS

Our results show that proactive and reactive stress coping individuals display distinct basal neurotranscriptomic states. Differences in baseline expression of select genes or regulation of specific gene modules are linked to the magnitude of the behavioral response and the display of a coping style, respectively. Our results expand the molecular mechanisms of stress coping from one focused on the neurotransmitter systems to a more complex system that involves an organism's capability to handle neurometabolic loads and allows for comparisons with other animal taxa to uncover potential conserved mechanisms.

The evolutionary legacy of size-selective harvesting extends from genes to populations

Size-selective harvesting is assumed to alter life histories of exploited fish populations, thereby negatively affecting population productivity, recovery, and yield. However, demonstrating that fisheries-induced phenotypic changes in the wild are at least partly genetically determined has proved notoriously difficult. Moreover, the population-level consequences of fisheries-induced evolution are still being controversially discussed. Using an experimental approach, we found that five generations of size-selective harvesting altered the life histories and behavior, but not the metabolic rate, of wild-origin zebrafish (Danio rerio). Fish adapted to high positively size selective fishing pressure invested more in reproduction, reached a smaller adult body size, and were less explorative and bold. Phenotypic changes seemed subtle but were accompanied by genetic changes in functional loci. Thus, our results provided unambiguous evidence for rapid, harvest-induced phenotypic and evolutionary change when harvesting is intensive and size selective. According to a life-history model, the observed life-history changes elevated population growth rate in harvested conditions, but slowed population recovery under a simulated moratorium. Hence, the evolutionary legacy of size-selective harvesting includes populations that are productive under exploited conditions, but selectively disadvantaged to cope with natural selection pressures that often favor large body size.

SNPfisher: tools for probing genetic variation in laboratory-reared zebrafish

Single nucleotide polymorphisms (SNPs) are the benchmark molecular markers for modern genomics. Until recently, relatively few SNPs were known in the zebrafish genome. The use of next-generation sequencing for the positional cloning of zebrafish mutations has increased the number of known SNP positions dramatically. Still, the identified SNP variants remain under-utilized, owing to scant annotation of strain specificity and allele frequency. To address these limitations, we surveyed SNP variation in three common laboratory zebrafish strains using whole-genome sequencing. This survey identified an average of 5.04 million SNPs per strain compared with the Zv9 reference genome sequence. By comparing the three strains, 2.7 million variants were found to be strain specific, whereas the remaining variants were shared among all (2.3 million) or some of the strains. We also demonstrate the broad usefulness of our identified variants by validating most in independent populations of the same laboratory strains. We have made all of the identified SNPs accessible through 'SNPfisher', a searchable online database (snpfisher.nichd.nih.gov). The SNPfisher website includes the SNPfisher Variant Reporter tool, which provides the genomic position, alternate allele read frequency, strain specificity, restriction enzyme recognition site changes and flanking primers for all SNPs and Indels in a user-defined gene or region of the zebrafish genome. The SNPfisher site also contains links to display our SNP data in the UCSC genome browser. The SNPfisher tools will facilitate the use of SNP variation in zebrafish research as well as vertebrate genome evolution.

Advancing biology through a deeper understanding of zebrafish ecology and evolution

Over the last two decades, the zebrafish has joined the ranks of premier model organisms for biomedical research, with a full suite of tools and genomic resources. Yet we still know comparatively little about its natural history. Here I review what is known about the natural history of the zebrafish, where significant gaps in our knowledge remain, and how a fuller appreciation of this organism's ecology and behavior, population genetics, and phylogeny can inform a variety of research endeavors.

Applications in the search for genomic selection signatures in fish

Selection signatures are genomic regions harboring DNA sequences functionally involved in the genetic variation of traits subject to selection. Selection signatures have been intensively studied in recent years because of their relevance to evolutionary biology and their potential association with genes that control phenotypes of interest in wild and domestic populations. Selection signature research in fish has been confined to a smaller scale, due in part to the relatively recent domestication of fish species and limited genomic resources such as molecular markers, genetic mapping, DNA sequences, and reference genomes. However, recent genomic technology advances are paving the way for more studies that may contribute to the knowledge of genomic regions underlying phenotypes of biological and productive interest in fish.

Population genomic analyses of early-phase Atlantic Salmon (Salmo salar) domestication/captive breeding

Domestication can have adverse genetic consequences, which may reduce the fitness of individuals once released back into the wild. Many wild Atlantic salmon (Salmo salarL.) populations are threatened by anthropogenic influences, and they are supplemented with captively bred fish. The Atlantic salmon is also widely used in selective breeding programs to increase the mean trait values for desired phenotypic traits. We analyzed a genomewide set of SNPs in three domesticated Atlantic salmon strains and their wild conspecifics to identify loci underlying domestication. The genetic differentiation between domesticated strains and wild populations was low (F_ST < 0.03), and domesticated strains harbored similar levels of genetic diversity compared to their wild conspecifics. Only a few loci showed footprints of selection, and these loci were located in different linkage groups among the different wild population/hatchery strain comparisons. Simulated scenarios indicated that differentiation in quantitative trait loci exceeded that in neutral markers during the early phases of divergence only when the difference in the phenotypic optimum between populations was large. This study indicates that detecting selection using standard approaches in the early phases of domestication might be challenging unless selection is strong and the traits under selection show simple inheritance patterns.

Phylogeny of zebrafish, a "model species," within Danio, a "model genus"

Zebrafish (Danio rerio) is an important model for vertebrate development, genomics, physiology, behavior, toxicology, and disease. Additionally, work on numerous Danio species is elucidating evolutionary mechanisms for morphological development. Yet, the relationships of zebrafish and its closest relatives remain unclear possibly due to incomplete lineage sorting, speciation with gene flow, and interspecies hybridization. To clarify these relationships, we first constructed phylogenomic data sets from 30,801 restriction-associated DNA (RAD)-tag loci (483,026 variable positions) with clear orthology to a single location in the sequenced zebrafish genome. We then inferred a well-supported species tree for Danio and tested for gene flow during the diversification of the genus. An approach independent of the sequenced zebrafish genome verified all inferred relationships. Although identification of the sister taxon to zebrafish has been contentious, multiple RAD-tag data sets and several analytical methods provided strong evidence for Danio aesculapii as the most closely related extant zebrafish relative studied to date. Data also displayed patterns consistent with gene flow during speciation and postspeciation introgression in the lineage leading to zebrafish. The incorporation of biogeographic data with phylogenomic analyses put these relationships in a phylogeographic context and supplied additional support for D. aesculapii as the sister species to D. rerio. The clear resolution of this study establishes a framework for investigating the evolutionary biology of Danio and the heterogeneity of genome evolution in the recent history of a model organism within an emerging model genus for genetics, development, and evolution.

Wild sex in zebrafish: loss of the natural sex determinant in domesticated strains

Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RAD-tag population genomics to identify sex-linked polymorphisms. After verifying this "RAD-sex" method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, "male genotypes" became males and some "female genotypes" also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.

Footprints of directional selection in wild Atlantic salmon populations: evidence for parasite-driven evolution?

Mechanisms of host-parasite co-adaptation have long been of interest in evolutionary biology; however, determining the genetic basis of parasite resistance has been challenging. Current advances in genome technologies provide new opportunities for obtaining a genome-scale view of the action of parasite-driven natural selection in wild populations and thus facilitate the search for specific genomic regions underlying inter-population differences in pathogen response. European populations of Atlantic salmon (Salmo salar L.) exhibit natural variance in susceptibility levels to the ectoparasite Gyrodactylus salaris Malmberg 1957, ranging from resistance to extreme susceptibility, and are therefore a good model for studying the evolution of virulence and resistance. However, distinguishing the molecular signatures of genetic drift and environment-associated selection in small populations such as land-locked Atlantic salmon populations presents a challenge, specifically in the search for pathogen-driven selection. We used a novel genome-scan analysis approach that enabled us to i) identify signals of selection in salmon populations affected by varying levels of genetic drift and ii) separate potentially selected loci into the categories of pathogen (G. salaris)-driven selection and selection acting upon other environmental characteristics. A total of 4631 single nucleotide polymorphisms (SNPs) were screened in Atlantic salmon from 12 different northern European populations. We identified three genomic regions potentially affected by parasite-driven selection, as well as three regions presumably affected by salinity-driven directional selection. Functional annotation of candidate SNPs is consistent with the role of the detected genomic regions in immune defence and, implicitly, in osmoregulation. These results provide new insights into the genetic basis of pathogen susceptibility in Atlantic salmon and will enable future searches for the specific genes involved.