Heat-induced masculinization in domesticated zebrafish is family-specific and yields a set of different gonadal transcriptomes.

Understanding environmental influences on sex ratios is important for the study of the evolution of sex-determining mechanisms and for evaluating the effects of global warming and chemical pollution. Fishes exhibit sexual plasticity, but the underlying mechanisms of environmental effects on their reproduction are unclear even in the well-established teleost research model, the zebrafish. Here we established the conditions to study the effects of elevated temperature on zebrafish sex. We showed that sex ratio response to elevated temperature is family-specific and typically leads to masculinization (female-to-male sex reversal), resulting in neomales. These results uncovered genotype-by-environment interactions that support a polygenic sex determination system in domesticated (laboratory) zebrafish. We found that some heat-treated fish had gene expression profiles similar to untreated controls of the same sex, indicating that they were resistant to thermal effects. Further, most neomales had gonadal transcriptomes similar to that of regular males. Strikingly, we discovered heat-treated females that displayed a normal ovarian phenotype but with a "male-like" gonadal transcriptome. Such major transcriptomic reprogramming with preserved organ structure has never been reported. Juveniles were also found to have a male-like transcriptome shortly after exposure to heat. These findings were validated by analyzing the expression of genes and signaling pathways associated with sex differentiation. Our results revealed a lasting thermal effect on zebrafish gonads, suggesting new avenues for detection of functional consequences of elevated temperature in natural fish populations in a global warming scenario.

Correction: Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

[This corrects the article DOI: 10.1371/journal.pgen.1005954.].

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding.

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Small-scale transcriptomics reveals differences among gonadal stages in Asian seabass (Lates calcarifer).

BACKGROUND: The Asian seabass (Lates calcarifer) is a protandrous hermaphrodite that typically matures as a male at approximately 2-4 years of age and then changes sex in subsequent years. Although several sexual maturation stages have been described histologically for both testis and ovary, the underlying gene expression profiles remain lacking. The development of a gene expression platform is therefore necessary to improve our understanding of the gonad development of this cultured teleost species. METHODS: Thirty Asian seabass gonads were collected from farms in Singapore, examined histologically and staged according to their sex and gonadal maturation status. Partial coding sequences of 24 sex-related genes were cloned using degenerate primers and were sequenced. Additional 13 cDNA sequences were obtained through next-generation sequencing. A real-time qPCR was then performed using the microfluidic-based Fluidigm 48.48 Dynamic arrays. RESULTS: We obtained 17 ovaries and 13 testes at various stages of sexual maturation. Of the 37 genes that were tested, 32 (86%) showed sexually dimorphic expression. These genes included sex-related genes, sox9, wt1, amh, nr5a2, dmrt1 and nr0b1, which showed testis-enhanced expression similar to other vertebrate species. Known male- and female-enhanced germ cells markers, which were established from studies in other species, similarly showed testis- and ovary-enhanced expression, respectively, in the Asian seabass. Three pro-Wnt signaling genes were also upregulated in the ovary, consistent with existing studies that suggested the role of Wnt signaling in ovarian differentiation in teleosts and mammals. The expression patterns of genes involved in steroidogenesis, retinoic acid metabolism, apoptosis and NF-κB signaling were also described. We were able to classify gonads according to sex and gonadal maturation stages, based on their small-scale transcriptomic profiles, and to uncover a wide variation in expression profiles among individuals of the same sex. CONCLUSIONS: The analysis of a selected set of genes related to reproduction and in sufficient number of individuals using a qPCR array can elucidate new insights into the molecular mechanisms involved in Asian seabass gonad development. Given the conservation of gene expression patterns found in this study, these insights may also help us draw parallels with other teleosts.

Infection of Tilapia tilapinevirus in Mozambique Tilapia (Oreochromis mossambicus), a Globally Vulnerable Fish Species.

Tilapia tilapinevirus, or tilapia lake virus (TiLV), is a highly contagious virus found in tilapia and its hybrid species that has been reported worldwide, including in Asia, the Americas, and Africa. In this study, we experimentally challenged Mozambique tilapia (Oreochromis mossambicus) with a virulent TiLV strain, VETKU-TV01, at both low (1 × 103 TCID50/mL) and high (1 × 105 TCID50/mL) concentration. After the challenge, the Mozambique tilapia showed pale skin with some hemorrhage and erosion, lethargy, abdominal swelling, congestion around the eye, and exophthalmos; there was a cumulative mortality rate at 48.89% and 77.78% in the groups that received the low and high concentration, respectively. Quantitative PCR and in situ hybridization confirmed the presence of TiLV in the internal organs of moribund fish. Notably, severe histopathological changes, including glycogen depletion, syncytial hepatic cells containing multiple nuclei and intracytoplasmic inclusion bodies, and infiltration of melanomacrophage into the spleen, were frequently found in the Mozambique tilapia challenged with high TiLV concentration. Comparatively, the infectivity and pathology of the TiLV infection in Mozambique tilapia and red hybrid tilapia (Oreochromis spp.) were found to be similar. Our results confirmed the susceptibility of Mozambique tilapia, which has recently been determined to be a vulnerable species, to TiLV infection, expanding knowledge that the virus can cause disease in this fish species.

Masculinization of Zebrafish Through Partial Depletion of Primordial Germ Cells by Injecting Diluted Morpholino Oligonucleotides into Embryos.

The regulation of reproduction in zebrafish, the prime model of fish research, is not fully understood. An efficient tool to gain a better understanding of this complicated process is utilization of severely sex-biased families or groups. Here, we describe a method for partial depletion of primordial germ cells (PGCs) that leads to eventual masculinization of zebrafish. The technique is based on injecting early embryos with diluted morpholino oligonucleotides that temporarily interfere with the production of Dead end (Dnd), an RNA-binding protein essential for PGC survival. In addition, we also propose the use of eviscerated trunk, as a suitable alternative for examining gonadal expression in juvenile zebrafish.

Heat Shock Factor 5 Is Essential for Spermatogenesis in Zebrafish.

Heat shock factors (Hsfs) are transcription factors that regulate responses to heat shock and other environmental stimuli. Four heat shock factors (Hsf1-4) have been characterized from vertebrates to date. In addition to stress response, they also play important roles in development and gametogenesis. Here, we study the fifth member of heat shock factor family, Hsf5, using zebrafish as a model organism. Mutant hsf5-/- males, generated by CRISPR/Cas9 technique, were infertile with drastically reduced sperm count, increased sperm head size, and abnormal tail architecture, whereas females remained fertile. We show that Hsf5 is required for progression through meiotic prophase 1 during spermatogenesis as suggested by the accumulation of cells in the leptotene and zygotene-pachytene stages and increased apoptosis in post-meiotic cells. hsf5-/- mutants show gonadal misregulation of a substantial number of genes with roles in cell cycle, apoptosis, protein modifications, and signal transduction, indicating an important role of Hsf5 in early stages of spermatogenesis.

The complete mitochondrial genome of a basal teleost, the Asian arowana (Scleropages formosus, Osteoglossidae).

BACKGROUND: Mitochondrial DNA-derived sequences have become popular markers for evolutionary studies, as their comparison may yield significant insights into the evolution of both the organisms and their genomes. From the more than 24,000 teleost species, only 254 complete mtDNA sequences are available (GenBank status on 06 Sep 2006). In this paper, we report the complete mitochondrial genome sequence of Asian arowana, a basal bonytongue fish species, which belongs to the order of Osteoglossiformes. RESULTS: The complete mitochondrial genomic sequence (mtDNA) of Asian arowana (Scleropages formosus) was determined by using shotgun sequencing method. The length of Asian arowana mtDNA is ca. 16,650 bp (its variation is due to polymorphic repeats in the control region), containing 13 protein-coding genes, 22 tRNA and 2 rRNA genes. Twelve of the thirteen protein coding genes were found to be encoded by the heavy strand in the order typically observed for vertebrate mitochondrial genomes, whereas only nad6 was located on the light strand. An interesting feature of Asian arowana mitogenome is that two different repeat arrays were identified in the control region: a 37 bp tandem repeat at the 5' end and an AT-type dinucleotide microsatellite at the 3' end. Both repeats show polymorphism among the six individuals tested; moreover the former one is present in the mitochondrial genomes of several other teleost groups. The TACAT motif described earlier only from mammals and lungfish was found in the tandem repeat of several osteoglossid and eel species. Phylogenetic analysis of fish species representing Actinopterygii and Sarcopterygii taxa has shown that the Asian arowana is located near the baseline of the teleost tree, confirming its status among the ancestral teleost lineages. CONCLUSION: The mitogenome of Asian arowana is very similar to the typical vertebrate mitochondrial genome in terms of gene arrangements, codon usage and base composition. However its control region contains two different types of repeat units at both ends, an interesting feature that to our knowledge has never been reported before for other vertebrate mitochondrial control regions. Phylogenetic analysis using the complete mtDNA sequence of Asian arowana confirmed that it belongs to an ancestral teleost lineage.

Nutrigenomic and Nutritional Analyses Reveal the Effects of Pelleted Feeds on Asian Seabass (Lates calcarifer).

As nutrition-related expenses constitute the majority of the costs for aquaculture farms, it is essential for them to use feeds that provide an ideal combination of nutrients for the species of choice. In this study, the relative effect of consuming three different pelleted feeds (B, C and D) in comparison to frozen baitfish (A; control) were compared on juvenile Asian seabass (77.3 ± 22.4g) that were selected for increased growth rate over two generations. Our objectives were: 1) to evaluate the effects of different pelleted feeds based on overall physiological changes and nutritional quality of fillets; 2) improve our understanding of the underlying mechanisms with transcriptomic analysis; 3) if possible, identify the feed type that supports the growth of these fishes without substantially reducing the nutritional quality of fillet. The growth performance, fatty acid composition of fillet, hepatic histology and transcriptome of the fishes (Groups A-D) were analyzed. The majority of fatty acids of the fillets, except γ-linolenic acid (GLA, C18:3n6), correlated significantly with the respective diets. Asian seabass fed Feed C showed highest specific growth rate (SGR) and feed conversion efficiency (FCE) with closest histology and transcriptomic profile to control, but their fillet contained the highest n6/n3 ratio. When the liver-based transcriptomes were analyzed, a complex set of differentially expressed genes were detected between groups fed pelleted feeds and controls as well as among the pellet-fed groups themselves. Significant enrichment of genes with growth-related function tallied with the morphological data measured. When compared with control (Group A), 'Biosynthesis of unsaturated fatty acids' and 'Steroid biosynthesis' pathways were significantly enriched in pellet-fed groups. Reduced goblet cell numbers were observed in the gut of pellet-fed fish compared to controls and fads6 was found to be a suitable candidate gene to separate wild-caught Asian seabass, from pellet-fed ones. These results provide insights for researchers on the various effects of feeds on the biochemistry and global gene expression of the fish and potentially for seabass farms to make more informed feed choices.

Zebrafish sex: a complicated affair.

In this review, we provide a detailed overview of studies on the elusive sex determination (SD) and gonad differentiation mechanisms of zebrafish (Danio rerio). We show that the data obtained from most studies are compatible with polygenic sex determination (PSD), where the decision is made by the allelic combinations of several loci. These loci are typically dispersed throughout the genome, but in some teleost species a few of them might be located on a preferential pair of (sex) chromosomes. The PSD system has a much higher level of variation of SD genotypes both at the level of gametes and the sexual genotype of individuals, than that of the chromosomal sex determination systems. The early sexual development of zebrafish males is a complicated process, as they first develop a 'juvenile ovary', that later undergoes a transformation to give way to a testis. To date, three major developmental pathways were shown to be involved with gonad differentiation through the modulation of programmed cell death. In our opinion, there are more pathways participating in the regulation of zebrafish gonad differentiation/transformation. Introduction of additional powerful large-scale genomic approaches into the analysis of zebrafish reproduction will result in further deepening of our knowledge as well as identification of additional pathways and genes associated with these processes in the near future.

Differential transcriptomic response in the spleen and head kidney following vaccination and infection of Asian seabass with Streptococcus iniae.

Vaccination is an important strategy in the protection of aquaculture species from major diseases. However, we still do not have a good understanding of the mechanisms underlying vaccine-induced disease resistance. This is further complicated by the presence of several lymphoid organs that play different roles when mounting an immune response. In this study, we attempt to elucidate some of these mechanisms using a microarray-based approach. Asian seabass (Lates calcarifer) were vaccinated against Streptococcus iniae and the transcriptomic changes within the spleen and head kidney at one and seven days post-vaccination were profiled. We subsequently challenged the seabass at three weeks post-vaccination with live S. iniae and similarly profiled the transcriptomes of the two organs after the challenge. We found that vaccination induced an early, but transient transcriptomic change in the spleens and a delayed response in the head kidneys, which became more similar to one another compared to un-vaccinated ones. When challenged with the pathogen, the spleen, but not the head kidneys, responded transcriptomically at 25-29 hours post-challenge. A unique set of genes, in particular those involved in the activation of NF-κB signaling, was up-regulated in the vaccinated spleens upon pathogen challenge but not in the un-vaccinated spleens. A semi-quantitative PCR detection of S. iniae using metagenomic DNA extracted from the water containing the seabass also revealed that vaccination resulted in reduction of pathogen shedding. This result indicated that vaccination not only led to a successful immune defense against the infection, but also reduced the chances for horizontal transmission of the pathogen. In conclusion, we have provided a transcriptomic analysis of how the teleost spleen and head kidneys responded to vaccination and subsequent infection. The different responses from the two organs are suggestive of their unique roles in establishing a vaccine-induced disease resistance.

Polygenic sex determination system in zebrafish.

BACKGROUND: Despite the popularity of zebrafish as a research model, its sex determination (SD) mechanism is still unknown. Most cytogenetic studies failed to find dimorphic sex chromosomes and no primary sex determining switch has been identified even though the assembly of zebrafish genome sequence is near to completion and a high resolution genetic map is available. Recent publications suggest that environmental factors within the natural range have minimal impact on sex ratios of zebrafish populations. The primary aim of this study is to find out more about how sex is determined in zebrafish. METHODOLOGY/PRINCIPAL FINDINGS: Using classical breeding experiments, we found that sex ratios across families were wide ranging (4.8% to 97.3% males). On the other hand, repeated single pair crossings produced broods of very similar sex ratios, indicating that parental genotypes have a role in the sex ratio of the offspring. Variation among family sex ratios was reduced after selection for breeding pairs with predominantly male or female offspring, another indication that zebrafish sex is regulated genetically. Further examinations by a PCR-based "blind assay" and array comparative genomic hybridization both failed to find universal sex-linked differences between the male and female genomes. Together with the ability to increase the sex bias of lines by selective breeding, these data suggest that zebrafish is unlikely to utilize a chromosomal sex determination (CSD) system. CONCLUSIONS/SIGNIFICANCE: Taken together, our study suggests that zebrafish sex is genetically determined with limited, secondary influences from the environment. As we have not found any sign for CSD in the species, we propose that the zebrafish has a polygenic sex determination system.

FluoMEP: a new genotyping method combining the advantages of randomly amplified polymorphic DNA and amplified fragment length polymorphism.

PCR-based identification of differences between two unknown genomes often requires complex manipulation of the templates prior to amplification and/or gel electrophoretic separation of a large number of samples with manual methods. Here, we describe a new genotyping method, called fluorescent motif enhanced polymorphism (fluoMEP). The fluoMEP method is based on random amplified polymorphic DNA (RAPD) assay, but combines the advantages of the large collection of unlabelled 10mer primers (ca. 5000) from commercial sources and the power of the automated CE devices used for the detection of amplified fragment length polymorphism (AFLP) patterns. The link between these two components is provided by a fluorescently labeled "common primer" that is used in a two-primer PCR together with an unlabeled RAPD primer. By using the same "common primer" and a series of RAPD primers, DNA templates can be screened quickly and effectively for polymorphisms. Our manuscript describes the optimization of the method and its characterization on different templates. We demonstrate by using several different approaches that the addition of the "common primer" to the PCR changes the profile of amplified fragments, allowing for screening various parts of the genome with the same set of unlabeled primers. We also present an in silico analysis of the genomic localization of fragments amplified by a RAPD primer with two different "common primers" and alone.