A microsatellite-based linkage map of salt tolerant tilapia (Oreochromis mossambicus x Oreochromis spp.) and mapping of sex-determining loci

Regulation of omega-3 fatty acids production by different genes in freshwater fish species: a review

The present study aims to compare the gene expression of three different fish species (common carp, tilapia, and trout) with varying levels of fatty acids (FA). Based on transcriptome analysis and RNA sequencing, various genes and their associated metabolic pathways are identified. Pathways are categorized based on the genes they encode. Genes that were differentially expressed and their promoter's methylation patterns were revealed by RNA-seq analysis in common carp. Furthermore, fatty acid-enriched pathways, such as ARA4 and adipocytokine signaling, were also identified. Many genes and pathways may influence tilapia's growth and omega-3 content. Using the mTOR pathway, trout with differential expression were discovered to be involved in producing omega-3 fatty acids. This study revealed major pathways in fish species to produce omega-3 fatty acids.

Genomics-driven prophylactic measures to increase streptococcosis resistance in tilapia

Streptococcosis caused by Streptococcus agalactiae and S. iniae is a significant problem that affects the success of tilapia aquaculture industries worldwide. In this critical review, we summarize the applicable practical strategies which may effectively enhance the world tilapia aquaculture development. Recently, the effect of vaccination and selective breeding programmes has been recognized as valuable tools to control the target disease and other consequent negative impacts caused by chemical and drug application. Advances in sequencing and molecular technologies are vital helpful factors with which to develop robust vaccines and increase the selective breeding programme's precision against streptococcosis. The genomic selection for streptococcosis-resistant tilapia strains and crucial genomic application for genomics' contribution to the development of novel Streptococcus vaccine, comparative genomics approach identifying vaccine candidates by reverse vaccinology, and next-generation vaccine design were described. Information from our review is encouraging for practical implementation of the development of vaccination and genomic selection in tilapia for streptococcosis resistance, which may be vital factors to sustain the world tilapia aquaculture industry effectively.

Microsatellites as Molecular Markers with Applications in Exploitation and Conservation of Aquatic Animal Populations

A large number of species and taxa has been studied for genetic polymorphism. Microsatellites have been known as hypervariable neutral molecular markers with the highest resolution power in comparison with any other markers. However, the discovery of a new type of molecular marker—single nucleotide polymorphism (SNP) has put the existing applications of microsatellites to the test. To ensure good resolution power in studies of populations and individuals, a number of microsatellite loci from 14 to 20 was often used, which corresponds to about 200 independent alleles. Recently, these numbers have tended to be increased by the application of genomic sequencing of expressed sequence tags (ESTs), and the choice of the most informative loci for genotyping depends on the aims of research. Examples of successful applications of microsatellite molecular markers in aquaculture, fisheries, and conservation genetics in comparison with SNPs have been summarized in this review. Microsatellites can be considered superior markers in such topics as kinship and parentage analysis in cultured and natural populations, the assessment of gynogenesis, androgenesis and ploidization. Microsatellites can be coupled with SNPs for mapping QTL. Microsatellites will continue to be used in research on genetic diversity in cultured stocks, and also in natural populations as an economically advantageous genotyping technique.

Master-Key Regulators of Sex Determination in Fish and Other Vertebrates-A Review

In vertebrates, mainly single genes with an allele ratio of 1:1 trigger sex-determination (SD), leading to initial equal sex-ratios. Such genes are designated master-key regulators (MKRs) and are frequently associated with DNA structural variations, such as copy-number variation and null-alleles. Most MKR knowledge comes from fish, especially cichlids, which serve as a genetic model for SD. We list 14 MKRs, of which dmrt1 has been identified in taxonomically distant species such as birds and fish. The identification of MKRs with known involvement in SD, such as amh and fshr, indicates that a common network drives SD. We illustrate a network that affects estrogen/androgen equilibrium, suggesting that structural variation may exert over-expression of the gene and thus form an MKR. However, the reason why certain factors constitute MKRs, whereas others do not is unclear. The limited number of conserved MKRs suggests that their heterologous sequences could be used as targets in future searches for MKRs of additional species. Sex-specific mortality, sex reversal, the role of temperature in SD, and multigenic SD are examined, claiming that these phenomena are often consequences of artificial hybridization. We discuss the essentiality of taxonomic authentication of species to validate purebred origin before MKR searches.

Absence of Figla-like Gene Is Concordant with Femaleness in Cichlids Harboring the LG1 Sex-Determination System

Oreochromis niloticus has been used as a reference genome for studies of tilapia sex determination (SD) revealing segregating genetic loci on linkage groups (LGs) 1, 3, and 23. The master key regulator genes (MKR) underlying the SD regions on LGs 3 and 23 have been already found. To identify the MKR in fish that segregate for the LG1 XX/XY SD-system, we applied short variant discovery within the sequence reads of the genomic libraries of the Amherst hybrid stock, Coptodon zillii and Sarotherodon galilaeus, which were aligned to a 3-Mbp-region of the O. aureus genome. We obtained 66,372 variants of which six were concordant with the XX/XY model of SD and were conserved across these species, disclosing the male specific figla-like gene. We further validated this observation in O. mossambicus and in the Chitralada hybrid stock. Genome alignment of the 1252-bp transcript showed that the figla-like gene’s size was 2664 bp, and that its three exons were capable of encoding 99 amino acids including a 45-amino-acid basic helix–loop–helix domain that is typical of the ovary development regulator—factor-in-the-germline-alpha (FIGLA). In Amherst gonads, the figla-like gene was exclusively expressed in testes. Thus, the figla-like genomic presence determines male fate by interrupting the female developmental program. This indicates that the figla-like gene is the long-sought SD MKR on LG1.

Sex determination mechanisms and sex control approaches in aquaculture animals

Aquaculture is one of the most efficient modes of animal protein production and plays an important role in global food security. Aquaculture animals exhibit extraordinarily diverse sexual phenotypes and underlying mechanisms, providing an ideal system to perform sex determination research, one of the important areas in life science. Moreover, sex is also one of the most valuable traits because sexual dimorphism in growth, size, and other economic characteristics commonly exist in aquaculture animals. Here, we synthesize current knowledge of sex determination mechanisms, sex chromosome evolution, reproduction strategies, and sexual dimorphism, and also review several approaches for sex control in aquaculture animals, including artificial gynogenesis, application of sex-specific or sex chromosome-linked markers, artificial sex reversal, as well as gene editing. We anticipate that better understanding of sex determination mechanisms and innovation of sex control approaches will facilitate sustainable development of aquaculture.

Transcriptome characterization and expression profile of Coix lacryma-jobi L. in response to drought

Coix lacryma-jobi L. is a very important economic crop widely cultivated in Southeast Asia. Drought affects more than four million square kilometers every year, and is a significant factor limiting agricultural productivity. However, relatively little is known about how Coix lacryma-jobi L. responds to drought treatments. To obtain a detailed and comprehensive understanding of the mechanisms regulating the transcriptional responses of Coix lacryma-jobi L. to drought treatment, we employed high throughput short-read sequencing of cDNA prepared from polyadenylated RNA to explore global gene expression after a seven-day drought treatment. We generated a de novo assembled transcriptome comprising 65,480 unique sequences. Differential expression analysis based on RSEM-estimated transcript abundances identified 5,315 differentially expressed genes (DEGs) when comparing samples from plants following drought-treatment and from the appropriate controls. Among these, the transcripts for 3,460 genes were increased in abundance, whereas 1,855 were decreased. Real-time quantitative PCR for 5 transcripts confirmed the changes identified by RNA-Seq. The results provide a transcriptional overview of the changes in Coix lacryma-jobi L. in response to drought, and will be very useful for studying the function of associated genes and selection of molecular marker of Coix lacryma-jobi L in the future.

Origin of a Giant Sex Chromosome

Chromosome size and morphology vary within and among species, but little is known about the proximate or ultimate causes of these differences. Cichlid fish species in the tribe Oreochromini share an unusual giant chromosome that is ∼3 times longer than the other chromosomes. This giant chromosome functions as a sex chromosome in some of these species. We test two hypotheses of how this giant sex chromosome may have evolved. The first hypothesis proposes that it evolved by accumulating repetitive elements as recombination was reduced around a dominant sex determination locus, as suggested by canonical models of sex chromosome evolution. An alternative hypothesis is that the giant sex chromosome originated via the fusion of an autosome with a highly repetitive B chromosome, one of which carried a sex determination locus. We test these hypotheses using comparative analysis of chromosome-scale cichlid and teleost genomes. We find that the giant sex chromosome consists of three distinct regions based on patterns of recombination, gene and transposable element content, and synteny to the ancestral autosome. The WZ sex determination locus encompasses the last ∼105 Mb of the 134-Mb giant chromosome. The last 47 Mb of the giant chromosome shares no obvious homology to any ancestral chromosome. Comparisons across 69 teleost genomes reveal that the giant sex chromosome contains unparalleled amounts of endogenous retroviral elements, immunoglobulin genes, and long noncoding RNAs. The results favor the B chromosome fusion hypothesis for the origin of the giant chromosome.

Pomc Plays an Important Role in Sexual Size Dimorphism in Tilapia

Sexual dimorphism is common across the animal kingdom. Knowledge of the mechanisms of sexual size dimorphism is limited although it is important in biology and aquaculture. Tilapia is the common name for ~ 100 species of cichlid fish. Some are important aquaculture species and males outgrow females. To gain novel insights into the mechanisms underlying sexual size dimorphism, we analyzed the differences of brain transcriptomes between males and females in Mozambique tilapia and studied the function of the pro-opiomelanocortin (Pomc) gene in tilapia and zebrafish. The transcriptome analysis identified 123, 55, and 2706 sex-biased genes at 5, 30, and 90 dph (days post-hatch), respectively, indicating sexual dimorphism of gene expressions in the brain. The expression of Pomc in the tilapia brain was a female-biased at 30, 90, and 120 dph. An analysis of the DNA sequence located upstream of the tilapia Pomc transcriptional start site identified two estrogenic response elements. In vitro luciferase assay of the two elements revealed that β-estradiol significantly enhanced the expression of luciferase activity, suggesting that the expression of Pomc is mediated by estrogen. We knocked out Pomc in zebrafish using Crispr/Cas-9. The Pomc-knockout zebrafish showed faster growth and higher sensitivity to feeding as compared to the wild-type fish. Taken together, our results indicate that Pomc contributes to sexual size dimorphism and suggest that the high estrogen level in females promotes the expression of Pomc and suppresses feeding in female tilapias, which leads to the slower growth of female tilapias.

A duplication of the Anti-Müllerian hormone gene is associated with genetic sex determination of different Oreochromis niloticus strains

Sex determination (SD) mechanisms are ancient and conserved, yet much diversity is exhibited in primary sex-determining signals that trigger male or female development. In O. niloticus, SD is associated with a male-specific locus on linkage group (LG) 23 which harbors the Y-linked Anti-Müllerian hormone (amh) gene, and a truncated duplication, denoted amhΔy. We have evaluated the possible role of identified indels and SNPs in the amh gene on SD, based on conservation in different O. niloticus strains. A fluorescent assay for the detection of a 5 bp insertion in amhΔy exon VI, efficiently discriminated between XX, XY, and YY genotypes. Concordance rate between amhΔy and sex varied in six Oreochromis strains, from 100% (Ghana) through 90% (Swansea) to 85% (Thai-Chitralada). The association of amhΔy with sex was found to be conserved in all tested O. niloticus strains, and thus supports its key role in SD. However, the previously identified missense SNP (C/T) in amh exon II was found only in the Swansea strain, thus excluding its candidacy for the causal variation of SD across all strains. Effects of markers on LGs 1, 3, and 23 (amhΔy) fully explained sex distribution in one Thai-Chitralada family (R2 = 1.0), whereas in another family only the major effect of LG23 (amhΔy) was significant (R2 = 0.37). Thus, amhΔy on LG23 is associated with genetic SD, either as a single causal gene in different O. niloticus strains, or in combination with segregating genes on LGs 1 and 3 in the Thai-Chitralada hybrid strain.

Vitellogenins and choriogenins are biomarkers for monitoring Oryzias curvinotus juveniles exposed to 17 β - estradiol

The effect of estrogens on Oryzias curvinotus juveniles were investigated by sequencing the transcriptome of O. curvinotus juveniles exposed to 17 β - estradiol for 24 h. A total of 69,071,524 and 71,210,528 raw reads were obtained for the control group (NC) and 17 β - estradiol exposure group (E2), respectively. After de novo assembly, total 133,210 unigenes were identified, and 85,837 unigenes (64.44% of 133,210) were annotated. Analysis of the transcriptome showed that exposure to 2 μg/L 17 β - estradiol led to the up-regulation of 19 genes and down-regulation of 18 genes. The eef1b and rps4x was most suitable as controls for quantitative real-time PCR (qPCR) using Reffinder. Different expression genes enrichment analysis found that exposed to 2 μg/L 17 β - estradiol affected various physiological processes, including spliceosome, phototransduction, amino sugar and nuclear sugar metabolism, hypotaurine metabolism, and renin-angiotensin system, etc. Exposing O. curvinotus juveniles to increasing concentrations of 17 β - estradiol (2 ng/L, 20 ng/L, 200 ng/L and 2 μg/L) led to significant up-regulation of vitellogenins (vtgs) and choriogenins (chgs) mRNA expression. The present study is the first high-throughput transcriptome sequencing of O. curvinotus juveniles, which will be useful for future functional analysis of genes related to environmental estrogen exposed, and development of biomarkers.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

Sex determination in the GIFT strain of tilapia is controlled by a locus in linkage group 23

Background

Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The ‘genetically improved farmed tilapia’ (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins.

Results

A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284).

Conclusions

Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.

Identifying a Long QTL Cluster Across chrLG18 Associated with Salt Tolerance in Tilapia Using GWAS and QTL-seq

Understanding the genetic mechanism of osmoregulation is important for the improvement of salt tolerance in tilapia. In our previous study, we have identified a major quantitative trait locus (QTL) region located at 23.0 Mb of chrLG18 in a Nile tilapia line by QTL-seq. However, the conservation of these QTLs in other tilapia populations or species is not clear. In this study, we successfully investigated the QTLs associated with salt tolerance in a mass cross population from the GIFT line of Nile tilapia (Oreochromis niloticus) using a ddRAD-seq-based genome-wide association study (GWAS) and in a full-sib family from the Malaysia red tilapia strain (Oreochromis spp) using QTL-seq. Our study confirmed the major QTL interval that is located at nearly 23.0 Mb of chrLG18 in Nile tilapia and revealed a long QTL cluster across chrLG18 controlling for the salt-tolerant trait in both red tilapia and Nile tilapia. This is the first GWAS analysis on salt tolerance in tilapia. Our finding provides important insights into the genetic architecture of salinity tolerance in tilapia and supplies a basis for fine mapping QTLs, marker-assisted selection, and further detailed functional analysis of the underlying genes for salt tolerance in tilapia.

Constructing a High-Density Genetic Linkage Map for Large Yellow Croaker (Larimichthys crocea) and Mapping Resistance Trait Against Ciliate Parasite Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is the most economically important marine cage-farming fish in China in the past decade. However, the sustainable development of large yellow croaker aquaculture has been severely hampered by several diseases, of which, the white spot disease caused by ciliate protozoan parasite Cryptocaryon irritans ranks the most damaging disease in large yellow croaker cage farms. To better understand the genetic basis of parasite infection and disease resistance to C. irritans, it is vital to map the traits and localize the underlying candidate genes in L. crocea genome. Here, we constructed a high-density genetic linkage map using double-digest restriction-site associated DNA (ddRAD)-based high-throughput SNP genotyping data of a F1 mapping family, which had been challenged with C. irritans for resistant trait measure. A total of 5261 SNPs was grouped and oriented into 24 linkage groups (LGs), representing 24 chromosomes of L. crocea. The total genetic map length was 1885.67 cM with an average inter-locus distance of 0.36 cM. Quantitative trait loci (QTL) mapping identified seven significant QTLs in four LGs linked to C. irritans disease resistance. Candidate genes underlying disease resistance were identified from the reference genome, including ifnar1, ifngr2, ikbke, and CD112. Comparative genomic analysis between large yellow croaker and the four closely related species revealed high evolutionary conservation of chromosomes, though inter-chromosomal rearrangements do exist. Especially, the croaker genome structure was closer to the medaka genome than stickleback, indicating that the croaker genome might retain the teleost ancestral genome structure. The high-density genetic linkage map provides an important tool and resource for fine mapping, comparative genome analysis, and molecular selective breeding of large yellow croaker.

Chromosome-scale assemblies reveal the structural evolution of African cichlid genomes

BACKGROUND

African cichlid fishes are well known for their rapid radiations and are a model system for studying evolutionary processes. Here we compare multiple, high-quality, chromosome-scale genome assemblies to elucidate the genetic mechanisms underlying cichlid diversification and study how genome structure evolves in rapidly radiating lineages.

RESULTS

We re-anchored our recent assembly of the Nile tilapia (Oreochromis niloticus) genome using a new high-density genetic map. We also developed a new de novo genome assembly of the Lake Malawi cichlid, Metriaclima zebra, using high-coverage Pacific Biosciences sequencing, and anchored contigs to linkage groups (LGs) using 4 different genetic maps. These new anchored assemblies allow the first chromosome-scale comparisons of African cichlid genomes. Large intra-chromosomal structural differences (∼2-28 megabase pairs) among species are common, while inter-chromosomal differences are rare (<10 megabase pairs total). Placement of the centromeres within the chromosome-scale assemblies identifies large structural differences that explain many of the karyotype differences among species. Structural differences are also associated with unique patterns of recombination on sex chromosomes. Structural differences on LG9, LG11, and LG20 are associated with reduced recombination, indicative of inversions between the rock- and sand-dwelling clades of Lake Malawi cichlids. M. zebra has a larger number of recent transposable element insertions compared with O. niloticus, suggesting that several transposable element families have a higher rate of insertion in the haplochromine cichlid lineage.

CONCLUSION

This study identifies novel structural variation among East African cichlid genomes and provides a new set of genomic resources to support research on the mechanisms driving cichlid adaptation and speciation.

Construction of a High-Density Linkage Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Channel Catfish (Ictalurus punctatus)

A high-density genetic linkage map is of particular importance in the fine mapping for important economic traits and whole genome assembly in aquaculture species. The channel catfish (Ictalurus punctatus), a species native to North America, is one of the most important commercial freshwater fish in the world. Outside of the United States, China has become the major producer and consumer of channel catfish after experiencing rapid development in the past three decades. In this study, based on restriction site associated DNA sequencing (RAD-seq), a high-density genetic linkage map of channel catfish was constructed by using single nucleotide polymorphisms (SNPs) in a F₁ family composed of 156 offspring and their two parental individuals. A total of 4,768 SNPs were assigned to 29 linkage groups (LGs), and the length of the linkage map reached 2,480.25 centiMorgans (cM) with an average distance of 0.55 cM between loci. Based on this genetic linkage map, 223 genomic scaffolds were anchored to the 29 LGs of channel catfish, and a total length of 704.66 Mb was assembled. Quantitative trait locus (QTL) mapping and genome-wide association analysis identified 10 QTLs of sex-related and six QTLs of growth-related traits at LG17 and LG28, respectively. Candidate genes associated with sex dimorphism, including spata2, spata5, sf3, zbtb38, and fox, were identified within QTL intervals on the LG17. A sex-linked marker with simple sequence repeats (SSR) in zbtb38 gene of the LG17 was validated for practical verification of sex in the channel catfish. Thus, the LG17 was considered as a sex-related LG. Potential growth-related genes were also identified, including important regulators such as megf9, npffr1, and gas1. In a word, we constructed the high-density genetic linkage map and developed the sex-linked marker in channel catfish, which are important genetic resources for future marker-assisted selection (MAS) of this economically important teleost.

Transcriptome Analysis Identified Genes for Growth and Omega-3/-6 Ratio in Saline Tilapia

Growth and omega-3/-6 ratio are important traits in aquaculture. The mechanisms underlying quick growth and high omega-3/-6 ratio in fish are not fully understood. The consumption of the meat of tilapia suffers a bad reputation due to its low omega-3/-6 ratio. To facilitate the improvement of these traits and to understand more about the mechanisms underlying quick growth and high omega-3/-6 ratio, we conducted transcriptome analysis in the muscle and liver of fast- and slow-growing hybrid saline tilapia generated by crossing Mozambique tilapia and red tilapia. A transcriptome with an average length of 963 bp was generated by using 486.65 million clean 100 bp paired-end reads. A total of 42,699 annotated unique sequences with an average length of 3.4 kb were obtained. Differentially expressed genes (DEGs) in the muscle and liver were identified between fast- and slow-growing tilapia. Pathway analysis classified these genes into many pathways. Ten genes, including foxK1, sparc, smad3, usp38, crot, fadps, sqlea, cyp7b1, impa1, and gss, from the DEGs were located within QTL for growth and omega-3, which were previously detected content in tilapia, suggesting that these ten genes could be important candidate genes for growth and omega-3 fatty acid content. Analysis of SNPs in introns 1 and 2 of foxK1 revealed that the SNPs were significantly associated with growth and omega-3/-6 ratio. This study lays the groundwork for further investigation of the molecular mechanisms underlying the phenotypic variation of these two traits and provides SNPs for selecting these traits at fingerling stage.

Marker-assisted selection of YY supermales from a genetically improved farmed tilapia-derived strain

Genetically improved farmed tilapia (GIFT) and GIFT-derived strains account for the majority of farmed tilapia worldwide. As male tilapias grow much faster than females, they are often considered more desirable in the aquacultural industry. Sex reversal of females to males using the male sex hormone 17-α-methyltestosterone (MT) is generally used to induce phenotypic males during large-scale production of all male fingerlings. However, the widespread use of large quantities of sex reversal hormone in hatcheries may pose a health risk to workers and ecological threats to surrounding environments. Breeding procedures to produce genetically all-male tilapia with limited or no use of sex hormones are therefore urgently needed. In this study, by applying marker-assisted selection (MAS) for the selection of YY supermales from a GIFT-derived strain, we identified 24 XY pseudofemale and 431 YY supermale tilapias. Further performance evaluation on the progenies of the YY supermales resulted in male rates of 94.1%, 99.5% and 99.6%, respectively, in three populations, and a daily increase in body weight of 1.4 g at 3 months (n=997). Our study established a highly effective MAS procedure in the selection of YY supermales from a GIFT-derived strain. Furthermore, the development of MAS-selected YY supermales will help reduce the utilization of hormones for controlling sex in the tilapia aquaculture.

A high-density genetic linkage map and QTL mapping for growth and sex of yellow drum (Nibea albiflora)

A high-density genetic linkage map is essential for the studies of comparative genomics and gene mapping, and can facilitate assembly of reference genome. Herein, we constructed a high-density genetic linkage map with 8,094 SNPs selected from 113 sequenced fish of a F1 family. Ultimately, the consensus map spanned 3818.24 cM and covered nearly the whole genome (99.4%) with a resolution of 0.47 cM. 1,457 scaffolds spanning 435.15 Mb were anchored onto 24 linkage groups, accounting for 80.7% of the draft genome assembly of the yellow drum. Comparative genomic analyses with medaka and zebrafish genomes showed superb chromosome-scale synteny between yellow drum and medaka. QTL mapping and association analysis congruously revealed 22 QTLs for growth-related traits and 13 QTLs for sex dimorphism. Some important candidate genes such as PLA2G4A, BRINP3 and P2RY1 were identified from these growth-related QTL regions. A gene family including DMRT1, DMRT2 and DMRT3 was identified from these sex-related QTL regions on the linkage group LG9. We demonstrate that this linkage map can facilitate the ongoing marker-assisted selection and genomic and genetic studies for yellow drum.

Development and Validation of 58K SNP-Array and High-Density Linkage Map in Nile Tilapia (O. niloticus)

Despite being the second most important aquaculture species in the world accounting for 7.4% of global production in 2015, tilapia aquaculture has lacked genomic tools like SNP-arrays and high-density linkage maps to improve selection accuracy and accelerate genetic progress. In this paper, we describe the development of a genotyping array containing more than 58,000 SNPs for Nile tilapia (Oreochromis niloticus). SNPs were identified from whole genome resequencing of 32 individuals from the commercial population of the Genomar strain, and were selected for the SNP-array based on polymorphic information content and physical distribution across the genome using the Orenil1.1 genome assembly as reference sequence. SNP-performance was evaluated by genotyping 4991 individuals, including 689 offspring belonging to 41 full-sib families, which revealed high-quality genotype data for 43,588 SNPs. A preliminary genetic linkage map was constructed using Lepmap2 which in turn was integrated with information from the O_niloticus_UMD1 genome assembly to produce an integrated physical and genetic linkage map comprising 40,186 SNPs distributed across 22 linkage groups (LGs). Around one-third of the LGs showed a different recombination rate between sexes, with the female being greater than the male map by a factor of 1.2 (1632.9 to 1359.6 cM, respectively), with most LGs displaying a sigmoid recombination profile. Finally, the sex-determining locus was mapped to position 40.53 cM on LG23, in the vicinity of the anti-Müllerian hormone (amh) gene. These new resources has the potential to greatly influence and improve the genetic gain when applying genomic selection and surpass the difficulties of efficient selection for invasively measured traits in Nile tilapia.

Genetic Mapping of Head Size Related Traits in Common Carp (Cyprinus carpio)

Head size is important economic trait for many aquaculture fish which is directly linked to their carcass yield. The genetic basis of head size trait remains unclear in many widely cultured fish species. Common carp (Cyprinus carpio) is one of the most widely studied fish due to its importance on both economic and environmental aspects. In this study, we performed genome-wide association study using 433 Yellow River carp individuals from multiple families to identify loci and genes potentially associated with head size related traits including head length (HL), head length/body length ratio (HBR), eye diameter (ED), and eye cross (EC). QTL mapping was utilized to filter the effects of population stratification and improve power for the candidates identification in the largest surveyed family with a published genetic linkage map. Twelve SNPs showed significant for head size traits in GWAS and 18 QTLs were identified in QTL mapping. Our study combining both GWAS and QTL mapping could compensate the deficiency from each other and advance our understanding of head size traits in common carp. To acquire a better understanding of the correlation between head size and body growth, we also performed comparisons between QTLs of head size traits and growth-related traits. Candidate genes underlying head size traits were identified surrounding the significant SNPs, including parvalbumin, srpk2, fsrp5, igf1, igf3, grb10, igf1r, notch2, sfrp2. Many of these genes have been identified with potential functions on bone formation and growth. Igf1 was a putative gene associated with both head size and body growth in Yellow River carp. The teleost-specific igf3 was a candidate head size related gene, related to both HL and HBR. Our study also indicated the importance of Igf signaling pathway for both growth and head size determination in common carp, which could be potentially used in future selective breeding in common carp as well as other species.

Quantitative trait loci on LGs 9 and 14 affect the reproductive interaction between two Oreochromis species, O. niloticus and O. aureus

Effective farming of tilapia requires all-male culture, characterized by uniformity and high growth rate. Males of O. aureus (Oa) and females of O. niloticus (On) produce all-male offspring, but there is a behavioral reproductive barrier between the two species that prevents mass production. In crosses between Oa and On broodstocks, few hybrid females are attracted to the Oa male nests (denoted responders), and if they harbor the On alleles for the sex determination (SD) sites on linkage groups (LGs) 1, 3, and 23, all-male progeny are produced. Yet, without controlling for the alleles underlying SD, the parental stocks gradually lose their capability for all-male production. Hypothesizing that marker-assisted selection for female responders would allow production of sustainable broodstocks, we applied genotyping-by-sequencing to generate 4983 informative SNPs from 13 responding and 28 non-responding females from two full-sib families. Accounting for multiple comparisons in a genome-wide association study, seven SNPs met a false discovery rate of 0.061. Lowest nominal probabilities were on LGs 9 and 14, for which microsatellite DNA markers were designed within the candidate genes PTGDSL and CASRL, respectively. By increasing the sample size to 22 responders and 47 non-responders and by genotyping additional established microsatellites, we confirmed the association of these LGs with female responsiveness. The combined effects of microsatellites GM171 and CARSL-LOC100690618 on LGs 9 and 14 explained 37% of the phenotypic variance of reproductive interaction (p < 0.0001). Based on these findings, we propose a strategy for mass production of all-male tilapia hybrids through selection for genomic loci affecting SD and female responsiveness.

Losing cichlid fish biodiversity: genetic and morphological homogenization of tilapia following colonization by introduced species

Among the many negative impacts of invasive species, hybridization with indigenous species has increasingly become recognized as a major issue. However, relatively few studies have characterized the phenotypic outcomes of hybridization following biological invasions. Here we investigate the genetic and morphological consequences of stocking invasive tilapia species in two water bodies in central Tanzania. We sampled individuals from the Mindu Reservoir on the Ruvu river system, and at Kidatu on the Great Ruaha-Rufiji river system. We screened individuals at 16 microsatellite loci, and quantified morphology using geometric morphometrics and linear measurements. In both the Mindu and Kidatu systems, we identified evidence of hybridization between indigenous Wami tilapia (Oreochromis urolepis) and the introduced Nile tilapia (Oreochromis niloticus) or blue-spotted tilapia (Oreochromis leucostictus). At both sites, purebred individuals could largely be separated using geometric morphometric variables, with hybrids occupying a broad morphospace among the parental species. Our data demonstrate that the gene pools and phenotypic identity of the indigenous O. urolepis have been severely impacted by the stocking of the invasive species. Given the lack of evidence for clear commercial benefits from stocking invasive tilapia species in waters already populated by indigenous congenerics, we suggest further spread of introduced species should be undertaken with considerable caution.

Transcriptome-seq provides insights into sex-preference pattern of gene expression between testis and ovary of the crucifix crab (Charybdis feriatus)

The crucifix crab, Charybdis feriatus, which mainly inhabits Indo-Pacific region, is regarded as one of the most high-potential species for domestication and incorporation into the aquaculture sector. However, the regulatory mechanisms of sex determination and differentiation of this species remain unclear. To identify candidate genes involved in sex determination and differentiation, high throughput sequencing of transcriptome from the testis and ovary of C. feriatus was performed by the Illumina platform. After removing adaptor primers, low-quality sequences and very short (<50 nt) reads, we obtained 80.9 million and 66.2 million clean reads from testis and ovary, respectively. A total of 86,433 unigenes were assembled, and ~43% (37,500 unigenes) were successfully annotated to the NR, NT, Swiss-Prot, KEGG, COG, GO databases. By comparing the testis and ovary libraries, we obtained 27,636 differentially expressed genes. Some candidate genes involved in the sex determination and differentiation of C. feriatus were identified, such as vasa, pgds, vgr, hsp90, dsx-f, fem-1, and gpr. In addition, 88,608 simple sequence repeats were obtained, and 61,929 and 77,473 single nucleotide polymorphisms from testis and ovary were detected, respectively. The transcriptome profiling was validated by quantitative real-time PCR in 30 selected genes, which showed a good consistency. The present study is the first high-throughput transcriptome sequencing of C. feriatus. These findings will be useful for future functional analysis of sex-associated genes and molecular marker-assisted selections in C. feriatus.

Mapping QTL for Omega-3 Content in Hybrid Saline Tilapia

Tilapia is one of most important foodfish species. The low omega-3 to omega-6 fatty acid ratio in freshwater tilapia meat is disadvantageous for human health. Increasing omega-3 content is an important task in breeding to increase the nutritional value of tilapia. However, conventional breeding to increase omega-3 content is difficult and slow. To accelerate the increase of omega-3 through marker-assisted selection (MAS), we conducted QTL mapping for fatty acid contents and profiles in a F₂ family of saline tilapia generated by crossing red tilapia and Mozambique tilapia. The total omega-3 content in F₂ hybrid tilapia was 2.5 ± 1.0 mg/g, higher than that (2.00 mg/g) in freshwater tilapia. Genotyping by sequencing (GBS) technology was used to discover and genotype SNP markers, and microsatellites were also genotyped. We constructed a linkage map with 784 markers (151 microsatellites and 633 SNPs). The linkage map was 2076.7 cM long and consisted of 22 linkage groups. Significant and suggestive QTL for total lipid content were mapped on six linkage groups (LG3, -4, -6, -8, -13, and -15) and explained 5.8-8.3% of the phenotypic variance. QTL for omega-3 fatty acids were located on four LGs (LG11, -18, -19, and -20) and explained 5.0 to 7.5% of the phenotypic variance. Our data suggest that the total lipid and omega-3 fatty acid content were determined by multiple genes in tilapia. The markers flanking the QTL for omega-3 fatty acids can be used in MAS to accelerate the genetic improvements of these traits in salt-tolerant tilapia.

Identifying a Major QTL Associated with Salinity Tolerance in Nile Tilapia Using QTL-Seq

Selection of new lines with high salinity tolerance allows for economically feasible production of tilapias in brackish water areas. Mapping QTLs and identifying the markers linked to salinity-tolerant traits are the first steps in the improvement of the tolerance in tilapia through marker-assisted selection techniques. By using QTL-seq strategy and linkage-based analysis, two significant QTL intervals (chrLG4 and chrLG18) on salinity-tolerant traits were firstly identified in the Nile tilapia. Fine mapping with microsatellite and SNP markers suggested a major QTL region that located at 23.0 Mb of chrLG18 and explained 79% of phenotypic variation with a LOD value of 95. Expression analysis indicated that at least 10 genes (e.g., LACTB2, KINH, NCOA2, DIP2C, LARP4B, PEX5R, and KCNJ9) near or within the QTL interval were significantly differentially expressed in intestines, brains, or gills under 10, 15, or 20 ppt challenges. Our findings suggest that QTL-seq can be effectively utilized in QTL mapping of salinity-tolerant traits in fish. The identified major QTL is a promising locus to improve our knowledge on the genetic mechanism of salinity tolerance in tilapia.

Genome-Wide QTL Analysis Identified Significant Associations Between Hypoxia Tolerance and Mutations in the GPR132 and ABCG4 Genes in Nile Tilapia

Exposure to hypoxia induces both acute and chronic stress responses, which plays an important role in health of cultured organisms including growth, reproduction, immunity, and other energy demanding activities. Application of advanced genomic technologies allows rapid identification of hypoxia trait-associated genes and precise selection of superior brood stocks with high tolerance in tilapia. By applying QTL-seq and double-digest restriction-site associated DNA sequencing (ddRAD-seq) techniques, we identified four genome-wide significant quantitative trait loci (QTLs) for hypoxia tolerance and many suggestive QTLs in Nile tilapia. These QTLs explained 6.6-14.7% of the phenotypic variance. Further analysis revealed that single nucleotide polymorphisms (SNPs) in exons of both GPR132 and ABCG4 genes located in genome-wide QTL intervals were significantly associated with hypoxia-tolerant traits. Expression analysis of both genes suggested that they were strong candidate genes involved into hypoxia tolerance in tilapia. Our findings suggest that both QTL-seq and ddRAD-seq techniques can be effectively utilized in QTL mapping of hypoxia traits in fish. Our data supply a basis for further marker-assisted selection of super lines with a high level of tolerance against low oxygen stress in the tilapia.

A High-Density Genetic Linkage Map and QTL Fine Mapping for Body Weight in Crucian Carp (Carassius auratus) Using 2b-RAD Sequencing

A high-resolution genetic linkage map is essential for a wide range of genetics and genomics studies such as comparative genomics analysis and QTL fine mapping. Crucian carp (Carassius auratus) is widely distributed in Eurasia, and is an important aquaculture fish worldwide. In this study, a high-density genetic linkage map was constructed for crucian carp using 2b-RAD technology. The consensus map contains 8487 SNP markers, assigning to 50 linkage groups (LGs) and spanning 3762.88 cM, with an average marker interval of 0.44 cM and genome coverage of 98.8%. The female map had 4410 SNPs, and spanned 3500.42 cM (0.79 cM/marker), while the male map had 4625 SNPs and spanned 3346.33 cM (0.72 cM/marker). The average recombination ratio of female to male was 2.13:1, and significant male-biased recombination suppressions were observed in LG47 and LG49. Comparative genomics analysis revealed a clear 2:1 syntenic relationship between crucian carp LGs and chromosomes of zebrafish and grass carp, and a 1:1 correspondence, but extensive chromosomal rearrangement, between crucian carp and common carp, providing evidence that crucian carp has experienced a fourth round of whole genome duplication (4R-WGD). Eight chromosome-wide QTL for body weight at 2 months after hatch were detected on five LGs, explaining 10.1-13.2% of the phenotypic variations. Potential candidate growth-related genes, such as an EGF-like domain and TGF-β, were identified within the QTL intervals. This high-density genetic map and QTL analysis supplies a basis for genome evolutionary studies in cyprinid fishes, genome assembly, and QTL fine mapping for complex traits in crucian carp.

Genome-wide mapping and characterization of microsatellites in the swamp eel genome

We described genome-wide screening and characterization of microsatellites in the swamp eel genome. A total of 99,293 microsatellite loci were identified in the genome with an overall density of 179 microsatellites per megabase of genomic sequences. The dinucleotide microsatellites were the most abundant type representing 71% of the total microsatellite loci and the AC-rich motifs were the most recurrent in all repeat types. Microsatellite frequency decreased as numbers of repeat units increased, which was more obvious in long than short microsatellite motifs. Most of microsatellites were located in non-coding regions, whereas only approximately 1% of the microsatellites were detected in coding regions. Trinucleotide repeats were most abundant microsatellites in the coding regions, which represented amino acid repeats in proteins. There was a chromosome-biased distribution of microsatellites in non-coding regions, with the highest density of 203.95/Mb on chromosome 8 and the least on chromosome 7 (164.06/Mb). The most abundant dinucleotides (AC)n was mainly located on chromosome 8. Notably, genomic mapping showed that there was a chromosome-biased association of genomic distributions between microsatellites and transposon elements. Thus, the novel dataset of microsatellites in swamp eel provides a valuable resource for further studies on QTL-based selection breeding, genetic resource conservation and evolutionary genetics.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Mapping QTL for Sex and Growth Traits in Salt-Tolerant Tilapia (Oreochromis spp. X O. mossambicus)

In aquaculture, growth and sex are economically important traits. To accelerate genetic improvement in increasing growth in salt-tolerant tilapia, we conducted QTL mapping for growth traits and sex with an F2 family, including 522 offspring and two parents. We used 144 polymorphic microsatellites evenly covering the genome of tilapia to genotype the family. QTL analyses were carried out using interval mapping for all individuals, males and females in the family, respectively. Using all individuals, three suggestive QTL for body weight, body length and body thickness respectively were detected in LG20, LG22 and LG12 and explained 2.4% to 3.1% of phenotypic variance (PV). When considering only males, five QTL for body weight were detected on five LGs, and explained 4.1 to 6.3% of PV. Using only females from the F2 family, three QTL for body weight were detected on LG1, LG6 and LG8, and explained 7.9-14.3% of PV. The QTL for body weight in males and females were located in different LGs, suggesting that in salt-tolerant tilapia, different set of genes 'switches' control the growth in males and females. QTL for sex were mapped on LG1 and LG22, indicating multigene sex determination in the salt-tolerant tilapia. This study provides new insights on the locations and effects of QTL for growth traits and sex, and sets the foundation for fine mapping for future marker-assisted selection for growth and sex in salt-tolerant tilapia aquaculture.

Genome-wide methylation analysis identified sexually dimorphic methylated regions in hybrid tilapia

Sexual dimorphism is an interesting biological phenomenon. Previous studies showed that DNA methylation might play a role in sexual dimorphism. However, the overall picture of the genome-wide methylation landscape in sexually dimorphic species remains unclear. We analyzed the DNA methylation landscape and transcriptome in hybrid tilapia (Oreochromis spp.) using whole genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq). We found 4,757 sexually dimorphic differentially methylated regions (DMRs), with significant clusters of DMRs located on chromosomal regions associated with sex determination. CpG methylation in promoter regions was negatively correlated with the gene expression level. MAPK/ERK pathway was upregulated in male tilapia. We also inferred active cis-regulatory regions (ACRs) in skeletal muscle tissues from WGBS datasets, revealing sexually dimorphic cis-regulatory regions. These results suggest that DNA methylation contribute to sex-specific phenotypes and serve as resources for further investigation to analyze the functions of these regions and their contributions towards sexual dimorphisms.

An ultra-high density linkage map and QTL mapping for sex and growth-related traits of common carp (Cyprinus carpio)

High density genetic linkage maps are essential for QTL fine mapping, comparative genomics and high quality genome sequence assembly. In this study, we constructed a high-density and high-resolution genetic linkage map with 28,194 SNP markers on 14,146 distinct loci for common carp based on high-throughput genotyping with the carp 250 K single nucleotide polymorphism (SNP) array in a mapping family. The genetic length of the consensus map was 10,595.94 cM with an average locus interval of 0.75 cM and an average marker interval of 0.38 cM. Comparative genomic analysis revealed high level of conserved syntenies between common carp and the closely related model species zebrafish and medaka. The genome scaffolds were anchored to the high-density linkage map, spanning 1,357 Mb of common carp reference genome. QTL mapping and association analysis identified 22 QTLs for growth-related traits and 7 QTLs for sex dimorphism. Candidate genes underlying growth-related traits were identified, including important regulators such as KISS2, IGF1, SMTLB, NPFFR1 and CPE. Candidate genes associated with sex dimorphism were also identified including 3KSR and DMRT2b. The high-density and high-resolution genetic linkage map provides an important tool for QTL fine mapping and positional cloning of economically important traits, and improving common carp genome assembly.

Charactering the ZFAND3 gene mapped in the sex-determining locus in hybrid tilapia (Oreochromis spp.)

Zinc finger AN1-type domain 3 (ZFAND3) is essential for spermatogenesis in mice. However, its function in teleosts remains unclear. In this study, we characterized the ZFAND3 gene (termed as OsZFAND3) in an important food fish, tilapia. The OsZFAND3 cDNA sequence is 1,050 bp in length, containing an ORF of 615 bp, which encodes a putative peptide of 204 amino acid residues. Quantitative real-time PCR revealed that the OsZFAND3 transcripts were exclusively expressed in the testis and ovary. In situ hybridization showed that the high expression of OsZFAND3 transcripts was predominantly localized in the spermatocyte and spermatid. These results suggest that OsZFAND3 is involved in male germ cell maturation. Three single nucleotide polymorphisms (SNPs) were detected in the introns of OsZFAND3. The OsZFAND3 gene was mapped in the sex-determining locus on linkage group 1 (LG1). The three SNPs in the OsZFAND3 gene were strictly associated with sex phenotype, suggesting that the OsZFAND3 gene is tightly linked to the sex-determining locus. Our study provides new insights into the functions of the OsZFAND3 gene in tilapia and a foundation for further detailed analysis of the OsZFAND3 gene in sex determination and differentiation.

Preliminary genetic linkage map of Indian major carp, Labeo rohita (Hamilton 1822) based on microsatellite markers

Linkage map with wide marker coverage is an essential resource for genetic improvement study for any species. Sex-averaged genetic linkage map of Labeo rohita, popularly known as 'rohu', widely cultured in the Indian subcontinent, was developed by placing 68 microsatellite markers generated by a simplified method. The parents and their F1 progeny (92 individuals) were used as segregating populations. The genetic linkage map spans a sex-averaged total length of 1462.2 cM, in 25 linkage groups. The genome length of rohu was estimated to be 3087.9 cM. This genetic linkage map may facilitate systematic searches of the genome to identify genes associated with commercially important characters and marker-assisted selection programmes of this species.

Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish

Background

Teleost fishes exhibit remarkably diverse and plastic sexual developmental patterns. One of the most astonishing is the rapid socially controlled female-to-male (protogynous) sex change observed in bluehead wrasses (Thalassoma bifasciatum). Such functional sex change is widespread in marine fishes, including species of commercial importance, yet its underlying molecular basis remains poorly explored.

Methods

RNA sequencing was performed to characterize the transcriptomic profiles and identify genes exhibiting sex-biased expression in the brain (forebrain and midbrain) and gonads of bluehead wrasses. Functional annotation and enrichment analysis were carried out for the sex-biased genes in the gonad to detect global differences in gene products and genetic pathways between males and females.

Results

Here we report the first transcriptomic analysis for a protogynous fish. Expression comparison between males and females reveals a large set of genes with sex-biased expression in the gonad, but relatively few such sex-biased genes in the brain. Functional annotation and enrichment analysis suggested that ovaries are mainly enriched for metabolic processes and testes for signal transduction, particularly receptors of neurotransmitters and steroid hormones. When compared to other species, many genes previously implicated in male sex determination and differentiation pathways showed conservation in their gonadal expression patterns in bluehead wrasses. However, some critical female-pathway genes (e.g., rspo1 and wnt4b) exhibited unanticipated expression patterns. In the brain, gene expression patterns suggest that local neurosteroid production and signaling likely contribute to the sex differences observed.

Conclusions

Expression patterns of key sex-related genes suggest that sex-changing fish predominantly use an evolutionarily conserved genetic toolkit, but that subtle variability in the standard sex-determination regulatory network likely contributes to sexual plasticity in these fish. This study not only provides the first molecular data on a system ideally suited to explore the molecular basis of sexual plasticity and tissue re-engineering, but also sheds some light on the evolution of diverse sex determination and differentiation systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-015-0044-8) contains supplementary material, which is available to authorized users.

Signatures of selection in tilapia revealed by whole genome resequencing

Natural selection and selective breeding for genetic improvement have left detectable signatures within the genome of a species. Identification of selection signatures is important in evolutionary biology and for detecting genes that facilitate to accelerate genetic improvement. However, selection signatures, including artificial selection and natural selection, have only been identified at the whole genome level in several genetically improved fish species. Tilapia is one of the most important genetically improved fish species in the world. Using next-generation sequencing, we sequenced the genomes of 47 tilapia individuals. We identified a total of 1.43 million high-quality SNPs and found that the LD block sizes ranged from 10–100 kb in tilapia. We detected over a hundred putative selective sweep regions in each line of tilapia. Most selection signatures were located in non-coding regions of the tilapia genome. The Wnt signaling, gonadotropin-releasing hormone receptor and integrin signaling pathways were under positive selection in all improved tilapia lines. Our study provides a genome-wide map of genetic variation and selection footprints in tilapia, which could be important for genetic studies and accelerating genetic improvement of tilapia.

De novo Transcriptome Analysis of Portunus trituberculatus Ovary and Testis by RNA-Seq: Identification of Genes Involved in Gonadal Development

The swimming crab Portunus trituberculatus is a commercially important crab species in East Asia countries. Gonadal development is a physiological process of great significance to the reproduction as well as commercial seed production for P. trituberculatus. However, little is currently known about the molecular mechanisms governing the developmental processes of gonads in this species. To open avenues of molecular research on P. trituberculatus gonadal development, Illumina paired-end sequencing technology was employed to develop deep-coverage transcriptome sequencing data for its gonads. Illumina sequencing generated 58,429,148 and 70,474,978 high-quality reads from the ovary and testis cDNA library, respectively. All these reads were assembled into 54,960 unigenes with an average sequence length of 879 bp, of which 12,340 unigenes (22.45% of the total) matched sequences in GenBank non-redundant database. Based on our transcriptome analysis as well as published literature, a number of candidate genes potentially involved in the regulation of gonadal development of P. trituberculatus were identified, such as FAOMeT, mPRγ, PGMRC1, PGDS, PGER4, 3β-HSD and 17β-HSDs. Differential expression analysis generated 5,919 differentially expressed genes between ovary and testis, among which many genes related to gametogenesis and several genes previously reported to be critical in differentiation and development of gonads were found, including Foxl2, Wnt4, Fst, Fem-1 and Sox9. Furthermore, 28,534 SSRs and 111,646 high-quality SNPs were identified in this transcriptome dataset. This work represents the first transcriptome analysis of P. trituberculatus gonads using the next generation sequencing technology and provides a valuable dataset for understanding molecular mechanisms controlling development of gonads and facilitating future investigation of reproductive biology in this species. The molecular markers obtained in this study will provide a fundamental basis for population genetics and functional genomics in P. trituberculatus and other closely related species.

A comprehensive transcriptome provides candidate genes for sex determination/differentiation and SSR/SNP markers in yellow catfish

Sex dimorphic growth pattern has significant theory and application implications in fish. Recently, a Y- and X-specific allele marker-assisted sex control technique has been developed for mass production of all-male population in yellow catfish (Pelteobagrus fulvidraco), but the genetic information for sex determination and sex control breeding has remained unclear. Here, we attempted to provide the first insight into a comprehensive transcriptome covering multiple tissues from XX females, XY males, and YY super-males of yellow catfish by using 454 GS-FLX platform, for a better assembly and gene coverage. A total of 1,202,933 high quality reads (about 540 Mbp) were obtained and assembled into 28,297 contigs and 141,951 singletons. BLASTX searches against the NCBI non-redundant protein database (nr) led a total of 52,564 unique sequences including 18,748 contigs and 33,816 singletons to match 25,669 known or predicted unique proteins. All of them with annotated function were categorized by gene ontology (GO) analysis, and 712 were assigned to reproduction and reproductive process. Some potential genes relevant to reproductive system including steroid hormone biosynthesis and GnRH (gonadotropin-releasing hormone) signaling pathway were further identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis; and at least 21 sex determination and differentiation-related genes, such as Dmrt1, Sox9a/b, Cyp19b, WT1, and AMH were identified and characterized. Additionally, a total of 82,794 simple sequence repeats (SSRs), 26,450 single nucleotide polymorphisms (SNPs), and 4,145 insertions and deletions (INDELs) were revealed from the transcriptome data. Therefore, the current transcriptome resources highlight further studies on sex-control breeding in yellow catfish and will benefit future studies on reproduction and sex determination in teleost fish.

A transcriptome derived female-specific marker from the invasive Western mosquitofish (Gambusia affinis)

Sex-specific markers are a prerequisite for understanding reproductive biology, genetic factors involved in sex differences, mechanisms of sex determination, and ultimately the evolution of sex chromosomes. The Western mosquitofish, Gambusia affinis, may be considered a model species for sex-chromosome evolution, as it displays female heterogamety (ZW/ZZ), and is also ecologically interesting as a worldwide invasive species. Here, de novo RNA-sequencing on the gonads of sexually mature G. affinis was used to identify contigs that were highly transcribed in females but not in males (i.e., transcripts with ovary-specific expression). Subsequently, 129 primer pairs spanning 79 contigs were tested by PCR to identify sex-specific transcripts. Of those primer pairs, one female-specific DNA marker was identified, Sanger sequenced and subsequently validated in 115 fish. Sequence analyses revealed a high similarity between the identified sex-specific marker and the 3´ UTR of the aminomethyl transferase (amt) gene of the closely related platyfish (Xiphophorus maculatus). This is the first time that RNA-seq has been used to successfully characterize a sex-specific marker in a fish species in the absence of a genome map. Additionally, the identified sex-specific marker represents one of only a handful of such markers in fishes.

Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish

The traits of cultured fish must continually be genetically improved to supply high-quality animal protein for human consumption. Economically important fish traits are controlled by multiple gene quantitative trait loci (QTL), most of which have minor effects, but a few genes may have major effects useful for molecular breeding. In this review, we chose relevant studies on some of the most intensively cultured fish and concisely summarize progress on identifying and verifying QTLs for such traits as growth, disease and stress resistance and sex in recent decades. The potential applications of these major-effect genes and their associated markers in marker-assisted selection and molecular breeding, as well as future research directions are also discussed. These genetic and genomic analyses will be valuable for elucidating the mechanisms modulating economically important traits and to establish more effective molecular breeding techniques in fish.

A syntenic region conserved from fish to Mammalian x chromosome

Sex chromosomes bearing the sex-determining gene initiate development along the male or female pathway, no matter which sex is determined by XY male or ZW female heterogamety. Sex chromosomes originate from ancient autosomes but evolved rapidly after the acquisition of sex-determining factors which are highly divergent between species. In the heterogametic male system (XY system), the X chromosome is relatively evolutionary silent and maintains most of its ancestral genes, in contrast to its Y counterpart that has evolved rapidly and degenerated. Sex in a teleost fish, the Nile tilapia (Oreochromis niloticus), is determined genetically via an XY system, in which an unpaired region is present in the largest chromosome pair. We defined the differences in DNA contents present in this chromosome with a two-color comparative genomic hybridization (CGH) and the random amplified polymorphic DNA (RAPD) approach in XY males. We further identified a syntenic segment within this region that is well conserved in several teleosts. Through comparative genome analysis, this syntenic segment was also shown to be present in mammalian X chromosomes, suggesting a common ancestral origin of vertebrate sex chromosomes.

Genetic mapping and QTL analysis of growth traits in the large yellow croaker Larimichthys crocea

Large yellow croaker (Larimichthys crocea) is an important maricultured species in China. A genetic linkage map of the large yellow croaker was constructed using type II microsatellites and expressed sequence tag (EST)-derived microsatellites in two half-sib families (two females and one male). A total of 289 microsatellite markers (contained 93 EST-SSRs) were integrated into 24 linkage groups, which agreed with the haploid chromosome number. The map spanned a length of 1,430.8 cm with an average interval of 5.4 cm, covering 83.9 % of the estimated genome size (1,704.8 cm). A total of seven quantitative trait locis (QTLs) were detected for growth traits on five linkage groups, including two 1 % and five 5 % chromosome-wide significant QTLs, and explained from 2.33 to 5.31 % of the trait variation. The identified QTLs can be applied in marker-assisted selection programs to improve the growth traits.

A genome scan revealed significant associations of growth traits with a major QTL and GHR2 in tilapia

Growth is an important trait in animal breeding. However, the genetic effects underpinning fish growth variability are still poorly understood. QTL mapping and analysis of candidate genes are effective methods to address this issue. We conducted a genome-wide QTL analysis for growth in tilapia. A total of 10, 7 and 8 significant QTLs were identified for body weight, total length and standard length at 140 dph, respectively. The majority of these QTLs were sex-specific. One major QTL for growth traits was identified in the sex-determining locus in LG1, explaining 71.7%, 67.2% and 64.9% of the phenotypic variation (PV) of body weight, total length and standard length, respectively. In addition, a candidate gene GHR2 in a QTL was significantly associated with body weight, explaining 13.1% of PV. Real-time qPCR revealed that different genotypes at the GHR2 locus influenced the IGF-1 expression level. The markers located in the major QTL for growth traits could be used in marker-assisted selection of tilapia. The associations between GHR2 variants and growth traits suggest that the GHR2 gene should be an important gene that explains the difference in growth among tilapia species.

The LBP gene and its association with resistance to Aeromonas hydrophila in tilapia

Resistance to pathogens is important for the sustainability and profitability of food fish production. In immune-related genes, the lipopolysaccharide-binding protein (LBP) gene is an important mediator of the inflammatory reaction. We analyzed the cDNA and genomic structure of the LBP gene in tilapia. The full-length cDNA (1901 bp) of the gene contained a 1416 bp open reading frame, encoding 471 amino acid residues. Its genomic sequence was 5577 bp, comprising 15 exons and 14 introns. Under normal conditions, the gene was constitutively expressed in all examined tissues. The highest expression was detected in intestine and kidney. We examined the responses of the gene to challenges with two bacterial pathogens Streptcoccus agalactiae and Aeromonas hydrophila. The gene was significantly upregulated in kidney and spleen post-infection with S. agalactiae and A. hydrophila, respectively. However, the expression profiles of the gene after the challenge with the two pathogens were different. Furthermore, we identified three SNPs in the gene. There were significant associations (p < 0.05) of two of the three SNPs with the resistance to A. hydrophila, but not with the resistance to S. agalactiae or growth performance. These results suggest that the LBP gene is involved in the acute-phase immunologic response to the bacterial infections, and the responses to the two bacterial pathogens are different. The two SNPs associated with the resistance to A. hydrophila may be useful in the selection of tilapia resistant to A. hydrophila.

Genetic architecture of sex determination in fish: applications to sex ratio control in aquaculture

Controlling the sex ratio is essential in finfish farming. A balanced sex ratio is usually good for broodstock management, since it enables to develop appropriate breeding schemes. However, in some species the production of monosex populations is desirable because the existence of sexual dimorphism, primarily in growth or first time of sexual maturation, but also in color or shape, can render one sex more valuable. The knowledge of the genetic architecture of sex determination (SD) is convenient for controlling sex ratio and for the implementation of breeding programs. Unlike mammals and birds, which show highly conserved master genes that control a conserved genetic network responsible for gonad differentiation (GD), a huge diversity of SD mechanisms has been reported in fish. Despite theory predictions, more than one gene is in many cases involved in fish SD and genetic differences have been observed in the GD network. Environmental factors also play a relevant role and epigenetic mechanisms are becoming increasingly recognized for the establishment and maintenance of the GD pathways. Although major genetic factors are frequently involved in fish SD, these observations strongly suggest that SD in this group resembles a complex trait. Accordingly, the application of quantitative genetics combined with genomic tools is desirable to address its study and in fact, when applied, it has frequently demonstrated a multigene trait interacting with environmental factors in model and cultured fish species. This scenario has notable implications for aquaculture and, depending upon the species, from chromosome manipulation or environmental control techniques up to classical selection or marker assisted selection programs, are being applied. In this review, we selected four relevant species or fish groups to illustrate this diversity and hence the technologies that can be used by the industry for the control of sex ratio: turbot and European sea bass, two reference species of the European aquaculture, and salmonids and tilapia, representing the fish for which there are well established breeding programs.

Allelic variant in the anti-Müllerian hormone gene leads to autosomal and temperature-dependent sex reversal in a selected Nile tilapia line

Owing to the demand for sustainable sex-control protocols in aquaculture, research in tilapia sex determination is gaining momentum. The mutual influence of environmental and genetic factors hampers disentangling the complex sex determination mechanism in Nile tilapia (Oreochromis niloticus). Previous linkage analyses have demonstrated quantitative trait loci for the phenotypic sex on linkage groups 1, 3, and 23. Quantitative trait loci for temperature-dependent sex reversal similarly reside on linkage group 23. The anti-Müllerian hormone gene (amh), located in this genomic region, is important for sexual fate in higher vertebrates, and shows sexually dimorphic expression in Nile tilapia. Therefore this study aimed at detecting allelic variants and marker-sex associations in the amh gene. Sequencing identified six allelic variants. A significant effect on the phenotypic sex for SNP ss831884014 (p<0.0017) was found by stepwise logistic regression. The remaining variants were not significantly associated. Functional annotation of SNP ss831884014 revealed a non-synonymous amino acid substitution in the amh protein. Consequently, a fluorescence resonance energy transfer (FRET) based genotyping assay was developed and validated with a representative sample of fish. A logistic linear model confirmed a highly significant effect of the treatment and genotype on the phenotypic sex, but not for the interaction term (treatment: p<0.0001; genotype: p<0.0025). An additive genetic model proved a linear allele substitution effect of 12% in individuals from controls and groups treated at high temperature, respectively. Moreover, the effect of the genotype on the male proportion was significantly higher in groups treated at high temperature, giving 31% more males on average of the three genotypes. In addition, the groups treated at high temperature showed a positive dominance deviation (+11.4% males). In summary, marker-assisted selection for amh variant ss831884014 seems to be highly beneficial to increase the male proportion in Nile tilapia, especially when applying temperature-induced sex reversal.