Genome-wide discovery of gene-related SNPs in Barramundi Lates calcarifer

Analysis of genetic diversity among three Triplophysa tenuis populations by RAD-seq

To investigate the genetic diversity of Triplophysa tenuis in the Shule River Basin of Gansu province, three populations were sequenced via RAD-seq technology. Twenty-nine microsatellite (SSR) markers with polymorphisms were finally screened to access the genetic diversity among the populations, of which 15 had high polymorphisms. The quantity of the alleles detected in the three populations of T. tenuis varied from 2 to 24. The locus with the most alleles was SSRC1, which had 24 alleles. Among the 29 SSRs, the range of effective allele number, observed heterozygosity, expected heterozygosity, and polymorphic information content were 1.246-16.615, 0.222-1, 0.198-0.940, and 0.178-0.937, respectively. Most of the identified loci were in the Hardy-Weinberg equilibrium. Analysis of the population structure revealed that the Yumen and Changma populations shared the same origin, while the Qiaowan population was different from them. The developed SSR markers discovered in this study will contribute to the conservation research on T. tenuis and the conservation of the fishery resources of the Shule River, providing scientific guidance for the development and utilization of T. tenuis resources and environmental protection.

Genomic Basis of Striking Fin Shapes and Colors in the Fighting Fish

Resolving the genomic basis underlying phenotypic variations is a question of great importance in evolutionary biology. However, understanding how genotypes determine the phenotypes is still challenging. Centuries of artificial selective breeding for beauty and aggression resulted in a plethora of colors, long-fin varieties, and hyper-aggressive behavior in the air-breathing Siamese fighting fish (Betta splendens), supplying an excellent system for studying the genomic basis of phenotypic variations. Combining whole-genome sequencing, quantitative trait loci mapping, genome-wide association studies, and genome editing, we investigated the genomic basis of huge morphological variation in fins and striking differences in coloration in the fighting fish. Results revealed that the double tail, elephant ear, albino, and fin spot mutants each were determined by single major-effect loci. The elephant ear phenotype was likely related to differential expression of a potassium ion channel gene, kcnh8. The albinotic phenotype was likely linked to a cis-regulatory element acting on the mitfa gene and the double-tail mutant was suggested to be caused by a deletion in a zic1/zic4 coenhancer. Our data highlight that major loci and cis-regulatory elements play important roles in bringing about phenotypic innovations and establish Bettas as new powerful model to study the genomic basis of evolved changes.

Scale Drop Disease Virus (SDDV) and Lates calcarifer Herpes Virus (LCHV) Coinfection Downregulate Immune-Relevant Pathways and Cause Splenic and Kidney Necrosis in Barramundi Under Commercial Farming Conditions

Marine farming of barramundi (Lates calcarifer) in Southeast Asia is currently severely affected by viral diseases. To better understand the biological implications and gene expression response of barramundi in commercial farming conditions during a disease outbreak, the presence of pathogens, comparative RNAseq, and histopathology targeting multiple organs of clinically “sick” and “healthy” juveniles were investigated. Coinfection of scale drop disease virus (SDDV) and L. calcarifer herpes virus (LCHV) were detected in all sampled fish, with higher SDDV viral loads in sick than in healthy fish. Histopathology showed that livers in sick fish often had moderate to severe abnormal fat accumulation (hepatic lipidosis), whereas the predominant pathology in the kidneys shows moderate to severe inflammation and glomerular necrosis. The spleen was the most severely affected organ, with sick fish presenting severe multifocal and coalescing necrosis. Principal component analysis (PC1 and PC2) explained 70.3% of the observed variance and strongly associated the above histopathological findings with SDDV loads and with the sick phenotypes, supporting a primary diagnosis of the fish being impacted by scale drop disease (SDD). Extracted RNA from kidney and spleen of the sick fish were also severely degraded likely due to severe inflammation and tissue necrosis, indicating failure of these organs in advanced stages of SDD. RNAseq of sick vs. healthy barramundi identified 2,810 and 556 differentially expressed genes (DEGs) in the liver and muscle, respectively. Eleven significantly enriched pathways (e.g., phagosome, cytokine-cytokine-receptor interaction, ECM-receptor interaction, neuroactive ligand-receptor interaction, calcium signaling, MAPK, CAMs, etc.) and gene families (e.g., tool-like receptor, TNF, lectin, complement, interleukin, chemokine, MHC, B and T cells, CD molecules, etc.) relevant to homeostasis and innate and adaptive immunity were mostly downregulated in sick fish. These DEGs and pathways, also previously identified in L. calcarifer as general immune responses to other pathogens and environmental stressors, suggest a failure of the clinically sick fish to cope and overcome the systemic inflammatory responses and tissue degeneration caused by SDD.

The HIF1αn gene and its association with hypoxia tolerance in the Asian seabass

Hypoxia is one of the major challenges in aquaculture industry. Breeding of fish tolerant to hypoxia is an important task in genetic improvement of aquaculture species. The Asian seabass, Lates calcarifer, is an important foodfish species. We identified and characterized the hypoxia-inducible factor inhibitor (HIF1αn) gene in the Asian seabass. The full-length cDNA sequence of the HIF1αn was 3425 bp, with an ORF of 1065 bp, encoding 354 amino acids. The genomic sequence of the gene was 8667 bp in length, and contained eight exons and seven introns. Phylogenetic analysis of the gene in fish and tetrapods revealed that the HIF1αn in the Asian seabass was closely related to that of tilapia (Oreochromis niloticus). The HIF1αn was highly up-regulated in the gill, spleen and heart after 3.5-hours hypoxia treatment. We identified three SNPs in the third and fourth introns of the HIF1αn gene. The SNP (i.e. SNP 9332241 (C/T)) in the fourth intron was significantly (P < 0.01) associated with hypoxia tolerance. This SNP might be useful in selecting Asian seabass for improved tolerance to hypoxia. Our data also provide useful information for further detailed analysis of the function of the HIF1αn gene in hypoxia tolerance.

Constructing High-Density Genetic Maps and Developing Sexing Markers in Northern Snakehead (Channa argus)

High-density genetic maps are essential for mapping QTL, improving genome assembly, comparative genomics, and studying sex chromosome evolution. The northern snakehead (Channa argus) is an economically important foodfish species with significant sexual dimorphism, where the males grow much faster and bigger than the females. However, to date, the sex determination pattern is still not clear, limiting identification of sex chromosomes, even sex determination genes and development of monosex populations that are valuable for both sex evolution of vertebrates and aquaculture practices. Here, a sex-averaged map and two sex-specific genetic maps were constructed with 2974, 2323, and 2338 SNPs, respectively. Little difference was observed in the pattern of sex-specific recombination between female- and male-specific genetic maps. Genome scan identified a major locus for sex determination at LG16. Females and males are, respectively, homogametic and heterogametic, suggesting an XY sex determination system for this species. By resequencing genomes, InDels in the sex-associated QTL region were discovered and used for developing sex-specific PCR assays for fast sexing of snakehead. These high-density genetic maps provide useful resources for future genomic studies in snakehead and its related species. The PCR assays for sexing are of importance in developing all male populations for aquaculture.

Does Circadian Rhythm Disorder Always Harm Organisms? Positive Effects of Circadian Rhythm Disorder in Lates calcarifer

BACKGROUND AND OBJECTIVE

Circadian rhythm disorder was proved that may make a result of a wide variety of adverse health consequences, however, does circadian rhythm disorder always harm organisms in their whole life process? The aim of the present study was to evaluate the influence of circadian rhythm disorder on barramundi Lates calcarifer at different life stages.

MATERIALS AND METHODS

The growth and immune responses of larval and juvenile barramundi Lates calcarifer under various circadian rhythm disorder situation (13, 16, 19 and 24 h light time) were investigated in this experiment.

RESULTS

The results showed that total length of larval and juvenile significantly increased comparing to the control group when the light time was extended to 16 h, but the length of fish in 24 h light time group showed significantly lower than the control. The cortisol level of larvae in 16 and 19 h light time groups significantly differed from other 2 groups. To juvenile, circadian rhythm disorder didn't make any influence on both T-SOD and ACP activity, the same result also appeared in cortisol concentration. The activity of three immune enzymes of larval from 16 and 19 h light time groups were all stimulated to increase by the homologous circadian rhythm disorder compared to the control.

CONCLUSION

In brief, lightly circadian rhythm disorder (prolonged 3 h light times) can benefit in fish growth. Although it brought some immune response to larvae which may cause adverse effects, all these negative effects would be eliminated when they grow up to juvenile.

Developing genome-wide SNPs and constructing an ultrahigh-density linkage map in oil palm

Oil palm (Elaeis guineensis Jacq.) is the leading oil-producing crops and the most important edible oil resource worldwide. DNA markers and genetic linkage maps are essential resources for marker-assisted selection to accelerate genetic improvement. We conducted RAD-seq on an Illumina NextSeq500 to discover genome-wide SNPs, and used the SNPs to construct a linkage map for an oil palm (Tenera) population derived from a cross between a Deli Dura and an AVROS Pisifera. The RAD-seq produced 1,076 million single-end reads across the breeding population containing 155 trees. Mining this dataset detected 510,251 loci. After filtering out loci with low accuracy and more than 20% missing data, 11,394 SNPs were retained. Using these SNPs, in combination with 188 anchor SNPs and 123 microsatellites, we constructed a linkage map containing 10,023 markers covering 16 chromosomes. The map length is 2,938.2 cM with an average marker space of 0.29 cM. The large number of SNPs will supply ample choices of DNA markers in analysing the genetic diversity, population structure and evolution of oil palm. This high-density linkage map will contribute to mapping quantitative trait loci (QTL) for important traits, thus accelerating oil palm genetic improvement.

Genome-Wide Association Study Identifies Loci Associated with Resistance to Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), is one major threat to mariculture. Identifying loci and understanding the mechanisms associated with resistance to VNN are important in selective breeding programs. We performed a genome-wide association study (GWAS) using genotyping-by-sequencing (GBS) to study the genomic architecture of resistance to NNV infection in Asian seabass. We genotyped 986 individuals from 43 families produced by 15 founders with 44498 bi-allelic genetic variants using GBS. The GWAS identified three genome-wide significant loci on chromosomes 16, 19, and 20, respectively, and six suggestive loci on chromosomes 1, 8, 14, 15, 21, and 24, respectively, associated with resistance to NNV infection measured as binary and quantitative traits. Using the 500 most significant markers in combination with a training population of 800 samples could reach a genomic prediction accuracy of 0.7. Candidate genes significantly associated with resistance to NNV, including lysine-specific demethylase 2A, beta-defensin 1, and cystatin-B, which play important roles in immune responses against virus infection, were identified. Almost all the candidate genes were differentially expressed in different tissues against NNV infection. The significant genetic variants can be used in genomic selection and help understand the mechanism of resistance to VNN. Future studies should use populations of large effective size and whole genome resequencing to identify more useful genetic variants.

Viral encephalopathy and retinopathy in aquaculture: a review

Viral encephalopathy and retinopathy (VER), otherwise known as viral nervous necrosis (VNN), is a major devastating threat for aquatic animals. Betanodaviruses have been isolated in at least 70 aquatic animal species in marine and in freshwater environments throughout the world, with the notable exception of South America. In this review, the main features of betanodavirus, including its diversity, its distribution and its transmission modes in fish, are firstly presented. Then, the existing diagnosis and detection methods, as well as the different control procedures of this disease, are reviewed. Finally, the potential of selective breeding, including both conventional and genomic selection, as an opportunity to obtain resistant commercial populations, is examined.

Construction of high-resolution recombination maps in Asian seabass

Background

A high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing.

Results

A high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males.

Conclusions

These high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3462-z) contains supplementary material, which is available to authorized users.

Construction of a high-density linkage map and fine mapping of QTL for growth in Asian seabass

A high-density genetic map is essential for comparative genomic studies and fine mapping of QTL, and can also facilitate genome sequence assembly. Here, a high density genetic map of Asian seabass was constructed with 3321 SNPs generated by sequencing 144 individuals in a F₂ family. The length of the map was 1577.67 cM with an average marker interval of 0.52 cM. A high level of genomic synteny among Asian seabass, European seabass, Nile tilapia and stickleback was detected. Using this map, one genome-wide significant and five suggestive QTL for growth traits were detected in six linkage groups (i.e. LG4, LG5, LG11, LG13, LG14 and LG15). These QTL explained 10.5–16.0% of phenotypic variance. A candidate gene, ACOX1 within the significant QTL on LG5 was identified. The gene was differentially expressed between fast- and slow-growing Asian seabass. The high-density SNP-based map provides an important tool for fine mapping QTL in molecular breeding and comparative genome analysis.