RNA-Seq identifies SNP markers for growth traits in rainbow trout

Identification of Single-Nucleotide Polymorphisms in Differentially Expressed Genes Favoring Soybean Meal Tolerance in Higher-Growth Zebrafish (Danio rerio)

Genetic variability within the same fish species could confer soybean meal (SBM) tolerance in some individuals, thus favoring growth. This study investigates the single-nucleotide polymorphisms (SNPs) in differentially expressed genes (DEGs) favoring SBM tolerance in higher-growth zebrafish (Danio rerio). In a previous work, nineteen families of zebrafish were fed a fish meal diet (100FM control diet) or SBM-based diets supplemented with saponin (50SBM + 2SPN-experimental diet), from juvenile to adult stages. Individuals were selected from families with a genotype-by-environment interaction higher (170 ± 18 mg) or lower (76 ± 10 mg) weight gain on 50SBM + 2SPN in relation to 100FM. Intestinal transcriptomic analysis using RNA-seq revealed six hundred and sixty-five differentially expressed genes in higher-growth fish fed 50SBM + 2SPN diet. In this work, using these results, 47 SNPs in DEGs were selected. These SNPs were genotyped by Sequenom in 340 zebrafish that were fed with a 50SBM + 2SPN diet or with 100FM diet. Marker-trait analysis revealed 4 SNPs associated with growth in 3 immunity-related genes (aif1l, arid3c, and cst14b.2) in response to the 50SBM + 2SPN diet (p-value < 0.05). Two SNPs belonging to aif1l y arid3c produce a positive (+19 mg) and negative (-26 mg) effect on fish growth, respectively. These SNPs can be used as markers to improve the early selection of tolerant fish to SBM diet or other plant-based diets. These genes can be used as biomarkers to identify SNPs in commercial fish, thus contributing to the aquaculture sustainability.

Combination of Transcriptomics and Metabolomics Analyses Provides Insights into the Mechanisms of Growth Differences in Spotted Seabass (Lateolabrax maculatus) Fed a Low-Phosphorus Diet

To analyze the potential mechanisms of growth differences in spotted seabass (Lateolabrax maculatus) fed a low-phosphorus diet, a total of 150 spotted seabass with an initial body weight of 4.49 ± 0.01 g were used (50 fish per tank) and fed a low-phosphorus diet for eight weeks. At the end of the experiment, five of the heaviest and five of the lightest fish were selected from each tank as fast-growing spotted seabass (FG) and slow-growing spotted seabass (SG), respectively, and their livers were analyzed by metabolomics and transcriptomics. The hepatic antioxidant capacity of the FG fed a low-phosphorus diet was significantly higher than that of the SG. A total of 431 differentially expressed genes (DEGs) were determined in the two groups, and most of the DEGs were involved in metabolism-related pathways such as steroid biosynthesis, glycolysis/gluconeogenesis, and protein digestion and absorption. Substance transport-related regulators and transporters were predominantly up-regulated. Furthermore, a large number of metabolites in the liver of FG were significantly up-regulated, especially amino acids, decanoyl-L-carnitine and dehydroepiandrosterone. The integration analysis of differential metabolites and genes further revealed that the interaction between protein digestion and absorption, as well as phenylalanine metabolism pathways were significantly increased in the liver of FG compared to those of the SG. In general, FG fed a low-phosphorus diet induced an enhancement in hepatic immune response, substance transport, and amino acid metabolism. This study provides new information on genetic mechanisms and regulatory pathways underlying differential growth rate and provides a basis for the foundation of efficient utilization of low-phosphorus diets and selective breeding programs for spotted seabass.

Development and application of a 1K functional liquid chip for lactation performance in Bactrian camels

Introduction

The advancement of high-throughput, high-quality, flexible, and cost-effective genotyping platforms is crucial for the progress of dairy breeding in Bactrian camels. This study focuses on developing and evaluating a 1K functional liquid single nucleotide polymorphism (SNP) array specifically designed for milk performance in Bactrian camels.

Methods

We utilized RNA sequencing data from 125 lactating camels to identify and select 1,002 loci associated with milk production traits for inclusion in the SNP array. The array's performance was then assessed using 24 randomly selected camels. Additionally, the array was employed to genotype 398 individuals, which allowed for population validation to assess the polymorphism of SNP sites.

Results

The SNP array demonstrated high overall SNP call rates (> 99%) and a remarkable 100% consistency in genotyping. Population validation results indicate that camels from six breeding areas in Northwest China share a similar genetic background regarding lactation functionality.

Discussion

This study highlights the potential of the SNP array to accelerate the breeding process of lactating Bactrian camels and provides a robust technical foundation for improving lactation performance.

Identification of pigmentation genes in skin, muscle and tail of a Thai-flag variety of Siamese fighting fish Betta splendens

Pigmentation genes expressed in skin, body muscle and tail of Thai-flag compared with Blue, White and Red varieties of Siamese fighting fish Betta splendens were identified. In total, 22,919 new unigenes were found. Pearson correlation and PCA analysis revealed that expression profiles of genes in muscle, skin and tail across solid color variety were similar. In contrast, those in skin and red tail part of Thai-flag were closely related but they showed different expression profiles with the white tail part. Moreover, 21,347-64,965 SNPs were identified in exonic regions of identified genes. In total, 28,899 genes were differentially expressed between paired comparisons of libraries where 13,907 genes (48.12 %) were upregulated and 14,992 genes (51.88 %) were downregulated. DEGs between paired libraries were 106-5775 genes relative to the compared libraries (56-2982 and 50-2782 for upregulated and downregulated DEGs). Interestingly, 432 pigmentation genes of B. splendens were found. Of these, 297 DEGs showed differential expression between varieties. Many DEGs in melanogenesis (Bsmcr1r, Bsmcr5r, and Bsslc2a15b), tyrosine metabolism (Bstyr, Bstyrp1b and Bsdct), stripe repressor (BsAsip1 and BsAsip2b), pteridine (Bsgch2) and carotenoid (BsBco2) biosynthesis were downregulated in the Thai-flag compared with solid color varieties. Expression of Bsbco1l, Bsfrem2b, Bskcnj13, Bszic2a and Bspah in skin, muscle and tail of Thai-flag, Blue, Red and White varieties was analyzed by qRT-PCR and revealed differential expression between fish varieties and showed anatomical tissue-preferred expression patterns in the same fish variety. The information could be applied to assist genetic-based development of new B. splendens varieties in the future.

FishSNP: a high quality cross-species SNP database of fishes

The progress of aquaculture heavily depends on the efficient utilization of diverse genetic resources to enhance production efficiency and maximize profitability. Single nucleotide polymorphisms (SNPs) have been widely used in the study of aquaculture genomics, genetics, and breeding research since they are the most prevalent molecular markers on the genome. Currently, a large number of SNP markers from cultured fish species are scattered in individual studies, making querying complicated and data reuse problematic. We compiled relevant SNP data from literature and public databases to create a fish SNP database, FishSNP ( http://bioinfo.ihb.ac.cn/fishsnp ), and also used a unified analysis pipeline to process raw data that the author of the literature did not perform SNP calling on to obtain SNPs with high reliability. This database presently contains 45,690,243 (45 million) nonredundant SNP data for 13 fish species, with 30,288,958 (30 million) of those being high-quality SNPs. The main function of FishSNP is to search, browse, annotate and download SNPs, which provide researchers various and comprehensive associated information.

Comparative transcriptomic analyses of brain-liver-muscle in channel catfish (Ictalurus punctatus) with differential growth rate

Growth is an important economically trait for aquatic animals. The popularity of farmed channel catfish (Ictalurus punctatus) in China has recently surged, prompting a need for research into the genetic mechanisms that drive growth and development to expedite the selection of fast-growing variants. In this study, the brain, liver and muscle transcriptomes of channel catfish between fast-growing and slow-growing groups were analyzed using RNA-Seq. Totally, 63, 110 and 86 differentially expressed genes (DEGs) were from brain, liver and muscle tissues. DEGs are primarily involved in growth, development, metabolism and immunity, which are related to the growth regulation of channel catfish, such as growth hormone receptor b (ghrb), fibroblast growth factor receptor 4 (fgfr4), bone morphogenetic protein 1a (bmp1a), insulin-like growth factor 2a (igf2a), collagen, type I, alpha 1a (col1a1a), acyl-CoA synthetase long chain family member 2 (acsl2) and caveolin 1 (cav1). This study advances our knowledge of the genetic mechanisms accounting for differences in growth rate and offers crucial gene resources for future growth-related molecular breeding programs in channel catfish.

Differential responses of hepatopancreas transcriptome between fast and slow growth in giant freshwater prawns (Macrobrachium rosenbergii) fed a plant-based diet

Efficient utilisation of plant-based diets in the giant freshwater prawn, Marcrobrachium rosenbergii, varies according to individual, suggesting that it might be associated with differences in physiological and metabolic responses. Therefore, we aimed to investigate the individual differences in the growth response of shrimp fed to a soybean-based diet (SBM). Two hundred shrimp were fed SBM for 90 days, and specific growth rate (SGR) was determined individually. Fast- and slow-growing shrimp (F-shrimp vs. S-shrimp), with the highest and lowest 5% SGRs, respectively, were sampled to determine haemolymph chemistry and carcass composition. The hepatopancreas of these shrimps were used for transcriptome analysis through RNA sequencing (RNA-Seq). The results showed no significant differences in haemolymph chemistry parameters. In terms of carcass proximate composition, F-shrimp exhibited higher protein composition than did S-shrimp, suggesting that F-shrimp have higher protein anabolism. Using RNA-seq and real-time reverse transcription polymerase chain reaction (qRT-PCR), the expression levels of several genes encoding physiologic and metabolic enzymes were found to be upregulated in F-shrimp compared to in S-shrimp, suggesting that these enzymes/proteins mediated the efficient use of SBM-based diets for growth promotion in shrimp. Various DEGs associated with the immune system were observed, indicating a difference in immune processes between F- and S-shrimp. The expression of several housekeeping genes was found to be upregulated in S-shrimp. Collectively, the upregulated expression of several enzymes associated with physiological and/or metabolic processes and increased protein anabolism may be attributed to the efficient use of SBM for maximal growth in shrimp.

Single nucleotide polymorphism SNP19140160 A > C is a potential breeding locus for fast-growth largemouth bass (Micropterus salmoides)

BACKGROUND

Largemouth bass (Micropterus salmoides) has significant economic value as a high-yielding fish species in China's freshwater aquaculture industry. Determining the major genes related to growth traits and identifying molecular markers associated with these traits serve as the foundation for breeding strategies involving gene pyramiding. In this study, we screened restriction-site associated DNA sequencing (RAD-seq) data to identify single nucleotide polymorphism (SNP) loci potentially associated with extreme growth differences between fast-growth and slow-growth groups in the F1 generation of a largemouth bass population.

RESULTS

We subsequently identified associations between these loci and specific candidate genes related to four key growth traits (body weight, body length, body height, and body thickness) based on SNP genotyping. In total, 4,196,486 high-quality SNPs were distributed across 23 chromosomes. Using a population-specific genotype frequency threshold of 0.7, we identified 30 potential SNPs associated with growth traits. Among the 30 SNPs, SNP19140160, SNP9639603, SNP9639605, and SNP23355498 showed significant associations; three of them (SNP9639603, SNP9639605, and SNP23355498) were significantly associated with one trait, body length, in the F1 generation, and one (SNP19140160) was significantly linked with four traits (body weight, height, length, and thickness) in the F1 generation. The markers SNP19140160 and SNP23355498 were located near two growth candidate genes, fam174b and ppip5k1b, respectively, and these candidate genes were closely linked with growth, development, and feeding. The average body weight of the group with four dominant genotypes at these SNP loci in the F1 generation population (703.86 g) was 19.63% higher than that of the group without dominant genotypes at these loci (588.36 g).

CONCLUSIONS

Thus, these four markers could be used to construct a population with dominant genotypes at loci related to fast growth. These findings demonstrate how markers can be used to identify genes related to fast growth, and will be useful for molecular marker-assisted selection in the breeding of high-quality largemouth bass.

Exploring the influence of DNA methylation and single nucleotide polymorphisms of the Myostatin gene on growth traits in the hybrid grouper (Epinephelus fuscoguttatus (female) × Epinephelus polyphekadion (male))

Investigations into the correlation between growth characteristics and DNA methylation levels, along with genetic variations, can provide fundamental insights to enhance growth performance in groupers. The Myostatin (mstn) gene plays a vital role in regulating skeletal muscle development and growth. This study scrutinized the DNA methylation levels of the mstn gene across hybrid groupers (E. fuscoguttatus (♀) × E. polyphekadion (♂)) and their parental species, to evaluate its impact on growth attributes in grouper fish. The nucleotide sequence of the mstn gene was directly sequenced in the hybrid grouper, exhibiting different growth performance to identify the single nucleotide polymorphisms (SNPs) of the mstn gene and explore their correlation with growth characteristics. The findings revealed no significant differences in global DNA methylation levels within muscle tissue among the hybrid grouper and parents. However, significant differences in DNA methylation sites were discovered between the hybrid grouper and E. polyphekadion at sites 824 and 1521 (located at exon 2 and intron 2, respectively), and between E. fuscoguttatus and E. polyphekadion at site 1521. These variations could potentially influence the mRNA expression of the mstn gene. The study also identified that SNP g.1003 T > C in exon 2 of the mstn gene was significantly associated with various growth traits including body weight, total length, body length, head length, caudal peduncle height, and body height (p < 0.01). Specimens with the TT genotype at site 1003 demonstrated superior growth performance compared to those with the TC genotype. Furthermore, microstructural analyses of muscle tissue showed that the average area and diameter of muscle fibers in TT genotype individuals were significantly greater than those in TC genotype individuals. Therefore, this research provides robust evidence linking the DNA methylation level and polymorphisms of the mstn gene with growth traits, which could be beneficial for grouper breeding programs.

Identification of differentially expressed genes and SNPs linked to harvest body weight of genetically improved rohu carp, Labeo rohita

In most of the aquaculture selection programs, harvest body weight has been a preferred performance trait for improvement. Molecular interplay of genes linked to higher body weight is not elucidated in major carp species. The genetically improved rohu carp with 18% average genetic gain per generation with respect to harvest body weight is a promising candidate for studying genes' underlying performance traits. In the present study, muscle transcriptome sequencing of two groups of individuals, with significant difference in breeding value, belonging to the tenth generation of rohu carp was performed using the Illumina HiSeq 2000 platform. A total of 178 million paired-end raw reads were generated to give rise to 173 million reads after quality control and trimming. The genome-guided transcriptome assembly and differential gene expression produced 11,86,119 transcripts and 451 upregulated and 181 downregulated differentially expressed genes (DEGs) between high-breeding value and low-breeding value (HB & LB) groups, respectively. Similarly, 39,158 high-quality coding SNPs were identified with the Ts/Tv ratio of 1.23. Out of a total of 17 qPCR-validated transcripts, eight were associated with cellular growth and proliferation and harbored 13 SNPs. The gene expression pattern was observed to be positively correlated with RNA-seq data for genes such as myogenic factor 6, titin isoform X11, IGF-1 like, acetyl-CoA, and thyroid receptor hormone beta. A total of 26 miRNA target interactions were also identified to be associated with significant DETs (p-value < 0.05). Genes such as Myo6, IGF-1-like, and acetyl-CoA linked to higher harvest body weight may serve as candidate genes in marker-assisted breeding and SNP array construction for genome-wide association studies and genomic selection.

Transcriptome comparison for identification of pigmentation-related genes in different color varieties of Siamese fighting fish Betta splendens

Transcriptome comparison was performed to identify genes expressed in skin, muscle and tails of mono-color (Red, Blue, Black, White and Yellow), bi-color (Cambodian) and multi-color (Marble) varieties of Siamese fighting fish Betta splendens. In total, 163,140 unigenes covering 26.348 Gb were found. Of these, 93,899 (57.55 %) unigenes significantly matched at least one database. In total, 5039 differentially expressed genes (DEGs) were found where 2415 genes (47.93 %) showed higher expression and 2624 genes (52.07 %) showed lower expression for all pairwise comparisons. DEGs between paired color varieties were 133-443. Of these, 38-220 genes were more highly expressed while 37-280 genes were more lowly expressed relative to the compared varieties. A total of 897 sequences (148 genes) significantly matched pigmentation-related genes of Danio rerio (E-value < 1e-06). Of these, 19 DEGs were identified. Examples are tyrosinase-related protein 1a (BsTyrp1a), epidermal growth factor receptor (BsEgfr) and neurofibronin 1a (BsNf1a). Moreover, 711,123 SNPs were identified and 1365 of these were located in pigmentation-related genes. Interestingly, an A > C₄₇₄ SNP in the gene BsTrpm7 and an indel (position 3571) in the BsItgb1a gene were found only in Cambodian. A C > T₂₅₂₀ SNP in BsFzd4 and 10 of 11 SNPs in BsTyrp1a were found only in Black. Different expression levels (P < 0.05) were found for tyrosinase (BsTyr), BsTyrp1a, BsNf1a and BsEgf1 among skin, body muscle and tails of the same variety and among the same tissues of different varieties (Red, Green, Blue, Black, Cambodian and Multi-colors, N = 5 each).

Development of Disease-Resistance-Associated Microsatellite DNA Markers for Selective Breeding of Tilapia (Oreochromis spp.) Farmed in Taiwan

There are numerous means to improve the tilapia aquaculture industry, and one is to develop disease resistance through selective breeding using molecular markers. In this study, 11 disease-resistance-associated microsatellite markers including 3 markers linked to hamp2, 4 linked to hamp1, 1 linked to pgrn2, 2 linked to pgrn1, and 1 linked to piscidin 4 (TP4) genes were established for tilapia strains farmed in Taiwan after challenge with Streptococcus inae. The correlation analysis of genotypes and survival revealed a total of 55 genotypes related to survival by the chi-square and Z-test. Although fewer markers were found in B and N2 strains compared with A strain, they performed well in terms of disease resistance. It suggested that this may be due to the low potency of some genotypes and the combinatorial arrangement between them. Therefore, a predictive model was built by the genotypes of the parental generation and the mortality rate of different combinations was calculated. The results show the same trend of predicted mortality in the offspring of three new disease-resistant strains as in the challenge experiment. The present findings is a nonkilling method without requiring the selection by challenge with bacteria or viruses and might increase the possibility of utilization of selective breeding using SSR markers in farms.

Transcriptome analysis of growth variation in early juvenile stage sandfish Holothuria scabra

The sandfish Holothuria scabra is a high-value tropical sea cucumber species representing a major mariculture prospect across the Indo-Pacific. Advancements in culture technology, rearing, and processing present options for augmenting capture production, stock restoration, and sustainable livelihood activities from hatchery-produced sandfish. Further improvements in mariculture production may be gained from the application of genomic technologies to improve performance traits such as growth. In this study, we performed de novo transcriptome assembly and characterization of fast- and slow-growing juvenile H. scabra from three Philippine populations. Analyses revealed 66 unigenes that were consistently differentially regulated in fast-growing sandfish and found to be associated with immune response and metabolism. Further, we identified microsatellite and single nucleotide polymorphism markers potentially associated with fast growth. These findings provide insight on potential genomic determinants underlying growth regulation in early juvenile sandfish which will be useful for further functional studies.

De novo assembly, transcriptome characterization and marker discovery in Indian major carp, Labeo rohita through pyrosequencing

Labeo rohita, one of the Indian major carps, is the most popular culture species in Indian subcontinent due to its consumer preference and delicacy. A selective breeding program for harvest body weight has resulted in an average genetic gain of 17% per generation. Transcriptome resource for this species is scanty. Here, we have characterized the liver and muscle transcriptomes of rohu using Roche 454 GS-FLX next generation sequencing platform. In total, 1.2 million reads were generated, de novo assembly and clustering resulted in 4171 transcripts. Out of these, 4171 had significant blast hit against NCBI nr database, and 2130 transcripts were successfully annotated. In total, 289 SSRs were identified with an identification rate of 5.8%, and dinucleotide repeat motifs were observed to be the most abundant SSRs. Further, 2231 putative SNPs were identified with high confidence. Validation of eight putative SNPs using Sanger sequencing resulted in 100% true SNPs. Significant allelic imbalance of M1, M4 and M5 loci between growth selected and control individual were observed. Furthermore, 13 transcription factors were identified in the present study belonging to six different transcription factor families. The present study demonstrated the utility of RNAseq to develop genomics resources in non-model fish species, and the marker resources developed would support the genetic improvement program of this species.

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.

Detection of selection signatures in the genome of a farmed population of anadromous rainbow trout (Oncorhynchus mykiss)

Domestication processes and artificial selection are likely to leave signatures that can be detected at a molecular level in farmed rainbow trout (Oncorhynchus mykiss). These signatures of selection are genomic regions that contain functional genetic variants conferring a higher fitness to their bearers. We genotyped 749 rainbow trout from a commercial population using a rainbow trout Axiom 57 K SNP array panel and identified putative genomic regions under selection using the pcadapt, Composite Likelihood Ratio (CLR) and Integrated Haplotype Score (iHS) methods. After applying quality-control pipelines and statistical analyses, we detected 12, 96 and 16 SNPs putatively under selection, associated with 96, 781 and 115 candidate genes, respectively. Several of these candidate genes were associated with growth, early development, reproduction, behavior and immune system traits. In addition, some of the SNPs were found in interesting regions located in autosomal inversions on Omy05 and Omy20. These findings could represent a genome-wide map of selection signatures in farmed rainbow trout and could be important in explaining domestication and selection for genetic traits of commercial interest.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

Construction of a chromosome-level genome assembly for genome-wide identification of growth-related quantitative trait loci in Sinocyclocheilus grahami (Cypriniformes, Cyprinidae)

The Dianchi golden-line barbel, Sinocyclocheilus grahami (Regan, 1904), is one of the “Four Famous Fishes” of Yunnan Province, China. Given its economic value, this species has been artificially bred successfully since 2007, with a nationally selected breed (“S. grahami, Bayou No. 1”) certified in 2018. For the future utilization of this species, its growth rate, disease resistance, and wild adaptability need to be improved, which could be achieved with the help of molecular marker-assisted selection (MAS). In the current study, we constructed the first chromosome-level genome of S. grahami, assembled 48 pseudo-chromosomes, and obtained a genome assembly of 1.49 Gb. We also performed QTL-seq analysis of S. grahami using the highest and lowest bulks (i.e., largest and smallest size) in both a sibling and random population. We screened two quantitative trait loci (QTLs) (Chr3, 14.9–39.1 Mb and Chr17, 4.1–27.4 Mb) as the major growth-related locations. Several candidate genes (e.g., map2k5, stat1, phf21a, sox6, and smad6) were also identified, with functions related to growth, such as cell differentiation, neuronal development, skeletal muscle development, chondrogenesis, and immunity. These results built a solid foundation for in-depth MAS studies on the growth traits of S. grahami.

Growth genes are implicated in the evolutionary divergence of sympatric piscivorous and insectivorous rainbow trout (Oncorhynchus mykiss)

Background

Identifying ecologically significant phenotypic traits and the genomic mechanisms that underly them are crucial steps in understanding traits associated with population divergence. We used genome-wide data to identify genomic regions associated with key traits that distinguish two ecomorphs of rainbow trout (Oncorhynchus mykiss)—insectivores and piscivores—that coexist for the non-breeding portion of the year in Kootenay Lake, southeastern British Columbia. “Gerrards” are large-bodied, rapidly growing piscivores with high metabolic rates that spawn north of Kootenay Lake in the Lardeau River, in contrast to the insectivorous populations that are on average smaller in body size, with lower growth and metabolic rates, mainly forage on aquatic insects, and spawn in tributaries immediately surrounding Kootenay Lake. We used pool-seq data representing ~ 60% of the genome and 80 fish per population to assess the level of genomic divergence between ecomorphs and to identify and interrogate loci that may play functional or selective roles in their divergence.

Results

Genomic divergence was high between sympatric insectivores and piscivores (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST = 0.188), and in fact higher than between insectivorous populations from Kootenay Lake and the Blackwater River (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST = 0.159) that are > 500 km apart. A window-based \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST analysis did not reveal “islands” of genomic differentiation; however, the window with highest \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_{\text{ST}}$$\end{document}FST estimate did include a gene associated with insulin secretion. Although we explored the use of the “Local score” approach to identify genomic outlier regions, this method was ultimately not used because simulations revealed a high false discovery rate (~ 20%). Gene ontology (GO) analysis identified several growth processes as enriched in genes occurring in the ~ 200 most divergent genomic windows, indicating many loci of small effect involved in growth and growth-related metabolic processes are associated with the divergence of these ecomorphs.

Conclusion

Our results reveal a high degree of genomic differentiation between piscivorous and insectivorous populations and indicate that the large body piscivorous phenotype is likely not due to one or a few loci of large effect. Rather, the piscivore phenotype may be controlled by several loci of small effect, thus highlighting the power of whole-genome resequencing in identifying genomic regions underlying population-level phenotypic divergences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01795-9.

Identifying Candidate Genes Involved in the Regulation of Early Growth Using Full-Length Transcriptome and RNA-Seq Analyses of Frontal and Parietal Bones and Vertebral Bones in Bighead Carp (Hypophthalmichthys nobilis)

Growth, one of the most important traits monitored in domestic animals, is essentially associated with bone development. To date, no large-scale transcriptome studies investigating bone development in bighead carp have been reported. In this study, we applied Isoform-sequencing technology to uncover the entire transcriptomic landscape of the bighead carp (Hypophthalmichthys nobilis) in early growth stage, and obtained 63,873 non-redundant transcripts, 20,907 long non-coding RNAs, and 1,579 transcription factors. A total of 381 alternative splicing events were seen in the frontal and parietal bones with another 784 events simultaneously observed in the vertebral bones. Coupling this to RNA sequencing (RNA-seq) data, we identified 27 differentially expressed unigenes (DEGs) in the frontal and parietal bones and 45 DEGs in the vertebral bones in the fast-growing group of fish, when compared to the slow-growing group of fish. Finally, 15 key pathways and 20 key DEGs were identified and found to be involved in regulation of early growth such as energy metabolism, immune function, and cytoskeleton function and important cellular pathways such as the arginine and proline metabolic pathway (p4ha1), FoxO signaling pathway (sgk1), cell adhesion molecules (b2m, ptprc, and mhcII), and peroxisome proliferator-activated receptor signaling pathway (scd). We established a novel full-length transcriptome resource and combined it with RNA-seq to elucidate the mechanism of genetic regulation of differential growth in bighead carp. The key DEGs identified in this study could fuel further studies investigating associations between growth and bone development and serve as a source of potential candidate genes for marker-assisted breeding programs.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish

Grass carp is one of the most important freshwater aquaculture species in China. However, the mechanisms underlying the growth of muscle tissue in the fish are unclear. High-throughput RNA-Seq was used to analyze the transcriptome of grass carp muscle tissue between fast- and slow-growing fish family groups. Twenty-four individuals each from 4 fast-growing families and 4 slow-growing families were used to reduce background noise. 71 up-regulated and 35 down-regulated genes were identified in the differentially expressed genes (DEGs). GO and KEGG enrichment analyses revealed the DEGs were involved in the GH/IGF axis, calcium metabolism, protein and glycogen synthesis, oxygen transport, cytoskeletal and myofibrillar components. IGFBP1 was up-regulated in big fish while GHR2 was down-regulated. Glutamic pyruvate transaminase 2, an indicator of liver tissue damage, was down-regulated in big grass carp, which indicates that the fish was better adapted to an artificially formulated diet. GAPDH, the rate-limiting enzyme in glycolytic flux was highly expressed in fast-growing grass carp, reflecting enhanced carbohydrate metabolism. Higher expression of ALAS2 and myoglobin 1 in big grass carp, related to oxygen transport might promote aerobic exercise along with food intake and muscle growth. Genes for cytoskeletal and myofibrillar components such as tropomyosin, meromyosin, and troponin I were also up-regulated in big grass carp. These results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for grass carp and contribute to our understanding of the molecular mechanisms and regulative pathways regulating growth in fish.

Genome-wide identification of loci associated with growth in rainbow trout

Background

Growth is a major economic production trait in aquaculture. Improvements in growth performance will reduce time and cost for fish to reach market size. However, genes underlying growth have not been fully explored in rainbow trout.

Results

A previously developed 50 K gene-transcribed SNP chip, containing ~ 21 K SNPs showing allelic imbalances potentially associated with important aquaculture production traits including body weight, muscle yield, was used for genotyping a total of 789 fish with available phenotypic data for bodyweight gain. Genotyped fish were obtained from two consecutive generations produced in the NCCCWA growth-selection breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform a genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with bodyweight gain. Using genomic sliding windows of 50 adjacent SNPs, 247 SNPs associated with bodyweight gain were identified. SNP-harboring genes were involved in cell growth, cell proliferation, cell cycle, lipid metabolism, proteolytic activities, chromatin modification, and developmental processes. Chromosome 14 harbored the highest number of SNPs (n = 50). An SNP window explaining the highest additive genetic variance for bodyweight gain (~ 6.4%) included a nonsynonymous SNP in a gene encoding inositol polyphosphate 5-phosphatase OCRL-1. Additionally, based on a single-marker GWA analysis, 33 SNPs were identified in association with bodyweight gain. The highest SNP explaining variation in bodyweight gain was identified in a gene coding for thrombospondin-1 (THBS1) (R² = 0.09).

Conclusion

The majority of SNP-harboring genes, including OCRL-1 and THBS1, were involved in developmental processes. Our results suggest that development-related genes are important determinants for growth and could be prioritized and used for genomic selection in breeding programs.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.

Single Nucleotide Polymorphism Analysis Using KASP Assay Reveals Genetic Variability in Rotylenchulus reniformis

This study employed single nucleotide polymorphisms (SNPs) to determine the genetic variability present in 26 isolates of Rotylenchulus reniformis from Louisiana, Mississippi, Arkansas, South Carolina, Georgia, Hawaii, and Alabama. Genomic DNA from reniform nematode was extracted and increased quantitatively using the process of whole genome amplification. More than 162 putative SNPs were identified, 31 of which were tested using a KASP kompetitive allele-specific PCR genotyping assay. Of the SNPs tested, 13, 17, and 19 SNPs revealed genetic variability within reniform nematode isolates from Louisiana, Mississippi, and Arkansas, respectively. Seven SNPs elucidated genetic differences among isolates of reniform nematode from Louisiana, Mississippi, and Arkansas. Eight SNPs determined genetic variability among individual isolates from South Carolina, Georgia, Hawaii, and Alabama. This study is the first to report genetic variability in geographic isolates of reniform nematode employing a SNP assay. This study also demonstrated that SNP markers can be used to evaluate isolates of R. reniformis and could be useful to assess their genetic diversity, origin, and distribution. Such information would be extremely useful in resistance breeding programs.

Comparative Transcriptome Analysis Reveals Molecular Basis Underlying Fast Growth of the Selectively Bred Pacific Oyster, Crassostrea gigas

Fast growth is one of the most desired traits for all food animals, which affects the profitability of animal production. The Pacific oyster, Crassostrea gigas, is an important aquaculture shellfish around the world with the largest annual production. Growth of the Pacific oyster has been greatly improved by artificial selection breeding, but molecular mechanisms underlying growth remains poorly understood, which limited the molecular integrative breeding of fast growth with other superior traits. In this study, comparative transcriptome analyses between the fast-growing selectively bred Pacific oyster and unselected wild Pacific oysters were conducted by RNA-Seq. A total of 1,303 protein-coding genes differentially expressed between fast-growing oysters and wild controls were identified, of which 888 genes were expressed at higher levels in the fast-growing oysters. Functional analysis of the differentially expressed genes (DEGs) indicated that genes involved in microtubule motor activity and biosynthesis of nucleotides and proteins are potentially important for growth in the oyster. Positive selection analysis of genes at the transcriptome level showed that a significant number of ribosomal protein genes had undergone positive selection during the artificial selection breeding process. These results also indicated the importance of protein biosynthesis and metabolism for the growth of oysters. The alternative splicing (AS) of genes was also compared between the two groups of oysters. A total of 3,230 differential alternative splicing events (DAS) were identified, involved in 1,818 genes. These DAS genes were associated with specific functional pathways related to growth, such as “long-term potentiation,” “salivary secretion,” and “phosphatidylinositol signaling system.” The findings of this study will be valuable resources for future investigation to unravel molecular mechanisms underlying growth regulation in the oyster and other marine invertebrates and to provide solid support for breeding application to integrate fast growth with other superior traits in the Pacific oyster.

Identification and genetic analysis of alternative splicing of long non-coding RNAs in tomato initial flowering stage

Alternative splicing (AS) is a key modulator of development in many eukaryotic organisms. Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs that play essential regulatory roles in various developmental processes and stress responses. However, the functions of AS lncRNAs during the initial flowering of tomato are largely unknown. This study was designed to investigate the AS pattern of lncRNAs in tomato flower, leaf, and root tissues at the initial flowering stage. Using RNA-Seq, we found that 72.55% of lncRNAs underwent AS in these tissues, yielding a total of 16,995 AS events. Among them, the main type of AS event is alternative first exon (AFE), followed by retained intron (RI). We performed candidate target genes analysis on tissue-specific AS lncRNA, and the results indicated that the candidate target genes of these lncRNAs may be involved in the regulation of circadian rhythm, plant immunity, cellulose synthesis and phosphate-containing compound metabolic process. Moreover, a total of 73,085 putative SNPs and 15,679 InDels were detected, and the potential relationship between the AS of lncRNAs and interesting SNP and InDel loci, as well as their numbers, revealed their effects on tomato genetic diversity and genomic stability. Our data provide new insights into the complexity of the transcriptome and the regulation of AS.

Comparative analysis of growth performance and liver transcriptome response of juvenile Ancherythroculter nigrocauda fed diets with different protein levels

This study aimed at investigating the effects of dietary protein levels on the growth and liver transcriptome in juvenile Ancherythroculter nigrocauda. Six semi-purified diets were formulated containing 25 (control), 30, 35, 40, 45, and 50% protein. Each diet was fed to three groups of 35 fish (mean initial weight: 5.86 ± 0.10 g) for 56 days. The rate of weight gain and specific growth rate increased with dietary protein levels from 25% to 40%, but remained unchanged when fed with 45 or 50% dietary protein. The feed conversion ratio was significantly influenced by the dietary protein levels, being the lowest in fish fed 40% protein. Illumina RNA-seq analysis was performed to investigate liver gene expression changes under different dietary protein treatments. A total of 367.78 million clean reads were obtained from the six libraries. Compared with 25% protein treatment library, there were 734, 1946, 1755, 2726, and 1523 upregulated genes, and 407, 1882, 1865, 2216 and 1624 downregulated genes in the 30, 35, 40, 45, and 50% protein treatment libraries, respectively. Trend analysis of these differentially expressed genes (DEGs) identified six statistically significant trends. A series of DEGs that related to protein metabolism, growth and development, lipid metabolism and immune and stress response were identified. Moreover, gene ontology enrichment analysis of the DEGs demonstrated that cellular process, single-organism process, metabolic process and biological regulation were the most highly overrepresented biological processes. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that protein processing in endoplasmic reticulum, PPAR signaling pathway, complement and coagulation cascades, and cytochrome P450 (CYP450s) were significantly enriched in the dietary protein treatment groups. Furthermore, qPCR results showed excellent agreement on those of RNA-seq for both up- and down-regulated genes (including fasn, accα, SCD, CPT-I, igf1, ST, AST, trdmt1, hsp70, cyp450, MHC-II, C4, tgfβ, ube4b, apoE and abcb7). Thus, our results provide the baseline information for the feed formulation and nutritional research for A. nigrocauda.

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

Cloning, SNP detection, and growth correlation analysis of the 5' flanking regions of two myosin heavy chain-7 genes in Mandarin fish (Siniperca chuatsi)

Myosin heavy chains (MYHs) play important roles in muscle growth and contraction. In fish, MYHs contribute to hyperplasia and hypertrophy of muscle fibers, which can continue into adult life and thus result in indeterminate growth in some species. We previously identified two MYH genes, MYH-7a and MYH-7b, that are differentially expressed in Mandarin fish (Siniperca chuatsi) and appear to function in early growth. However, the regulatory role of their 5' flanking regions is unknown. To examine the effects of single nucleotide polymorphisms (SNPs) in these regions, we used genome walking to amplify their flanking sequences and analyzed the regulatory elements and binding sites. A single SNP locus was found in the flanking sequence of each gene. These SNP loci are located in the conserved glucocorticoid receptor binding region (MYH-7a: G-614A; Allele frequency: G:A = 94.9:5.1; GG (89.76) and AG (10.24) genotypes) and the LIM homeobox domain transcription factor binding sequence (MYH-7b: C-1933A; Allele frequency: C:A = 54.8:45.2; AA (20.82), AC (48.81), and CC (30.37) genotypes). At the G-614A loci, the GG genotype exhibited more superior growth traits (total length, body length, body height, etc.) than the AG genotype, with the exception of caudal peduncle length. Alternatively, at the C-1933A loci, the AC and AA genotypes showed significant differences in all growth traits, except for head length, with AC exhibiting superior traits. The AA and CC genotypes showed significant differences in caudal peduncle length and height, while no differences were observed between the AC and CC genotypes. Thus, these SNPs in the 5' flanking regions of MYH-7a and MYH-7b are correlated with superior growth and can be used for selecting Mandarin fish during breeding.

Genome-Wide Association Analysis With a 50K Transcribed Gene SNP-Chip Identifies QTL Affecting Muscle Yield in Rainbow Trout

Detection of coding/functional SNPs that change the biological function of a gene may lead to identification of putative causative alleles within QTL regions and discovery of genetic markers with large effects on phenotypes. This study has two-fold objectives, first to develop, and validate a 50K transcribed gene SNP-chip using RNA-Seq data. To achieve this objective, two bioinformatics pipelines, GATK and SAMtools, were used to identify ~21K transcribed SNPs with allelic imbalances associated with important aquaculture production traits including body weight, muscle yield, muscle fat content, shear force, and whiteness in addition to resistance/susceptibility to bacterial cold-water disease (BCWD). SNPs ere identified from pooled RNA-Seq data collected from ~620 fish, representing 98 families from growth- and 54 families from BCWD-selected lines with divergent phenotypes. In addition, ~29K transcribed SNPs without allelic-imbalances were strategically added to build a 50K Affymetrix SNP-chip. SNPs selected included two SNPs per gene from 14K genes and ~5K non-synonymous SNPs. The SNP-chip was used to genotype 1728 fish. The average SNP calling-rate for samples passing quality control (QC; 1,641 fish) was ≥ 98.5%. The second objective of this study was to test the feasibility of using the new SNP-chip in GWA (Genome-wide association) analysis to identify QTL explaining muscle yield variance. GWA study on 878 fish (representing 197 families from 2 consecutive generations) with muscle yield phenotypes and genotyped for 35K polymorphic markers (passing QC) identified several QTL regions explaining together up to 28.40% of the additive genetic variance for muscle yield in this rainbow trout population. The most significant QTLs were on chromosomes 14 and 16 with 12.71 and 10.49% of the genetic variance, respectively. Many of the annotated genes in the QTL regions were previously reported as important regulators of muscle development and cell signaling. No major QTLs were identified in a previous GWA study using a 57K genomic SNP chip on the same fish population. These results indicate improved detection power of the transcribed gene SNP-chip in the target trait and population, allowing identification of large-effect QTLs for important traits in rainbow trout.

Population genetic analysis of aquaculture salmonid populations in China using a 57K rainbow trout SNP array

Various salmonid species are cultivated in cold water aquaculture. However, due to limited genomic data resources, specific high-throughput genotyping tools are not available to many of the salmonid species. In this study, a 57K single nucleotide polymorphism (SNP) array for rainbow trout (Oncorhynchus mykiss) was utilized to detect polymorphisms in seven salmonid species, including Hucho taimen, Oncorhynchus masou, Salvelinus fontinalis, Brachymystax lenok, Salvelinus leucomaenis, O. kisutch, and O. mykiss. The number of polymorphic markers per population ranged from 3,844 (O. kisutch) to 53,734 (O. mykiss), indicating that the rainbow trout SNP array was applicable as a universal genotyping tool for other salmonid species. Among the six other salmonid populations from four genera, 28,882 SNPs were shared, whereas 525 SNPs were polymorphic in all four genera. The genetic diversity and population relationships of the seven salmonid species were studied by principal component analysis (PCA). The phylogenetic relationships among populations were analyzed using the maximum likelihood method, which indicated that the shared SNP markers provide reliable genomic information for population genetic analyses in common aquaculture salmonid fishes. Furthermore, this obtained genomic information may be applicable for population genetic evaluation, marker-assisted breeding, and propagative parent selection in fry production.

Integrated analysis of lncRNA and mRNA expression in rainbow trout families showing variation in muscle growth and fillet quality traits

Muscle yield and quality traits are important for the aquaculture industry and consumers. Genetic selection for these traits is difficult because they are polygenic and result from multifactorial interactions. To study the genetic architecture of these traits, phenotypic characterization of whole body weight (WBW), muscle yield, fat content, shear force and whiteness were measured in ~500 fish representing 98 families from a growth-selected line. RNA-Seq was used to sequence the muscle transcriptome of different families exhibiting divergent phenotypes for each trait. We have identified 240 and 1,280 differentially expressed (DE) protein-coding genes and long noncoding RNAs (lncRNAs), respectively, in fish families exhibiting contrasting phenotypes. Expression of many DE lncRNAs (n = 229) was positively correlated with overlapping, neighboring or distantly located protein-coding genes (n = 1,030), resulting in 3,392 interactions. Three DE antisense lncRNAs were co-expressed with sense genes known to impact muscle quality traits. Forty-four DE lncRNAs had potential sponge functions to miRNAs that affect muscle quality traits. This study (1) defines muscle quality associated protein-coding and noncoding genes and (2) provides insight into non-coding RNAs involvement in regulating growth and fillet quality traits in rainbow trout.

Population resequencing reveals candidate genes associated with salinity adaptation of the Pacific oyster Crassostrea gigas

The Pacific oyster Crassostrea gigas is an important cultivated shellfish. As a euryhaline species, it has evolved adaptive mechanisms responding to the complex and changeable intertidal environment that it inhabits. To investigate the genetic basis of this salinity adaptation mechanism, we conducted a genome-wide association study using phenotypically differentiated populations (hyposalinity and hypersalinity adaptation populations, and control population), and confirmed our results using an independent population, high-resolution melting, and mRNA expression analysis. For the hyposalinity adaptation, we determined 24 genes, including Cg_CLCN7 (chloride channel protein 7) and Cg_AP1 (apoptosis 1 inhibitor), involved in the ion/water channel and transporter mechanisms, free amino acid and reactive oxygen species metabolism, immune responses, and chemical defence. Three SNPs located on these two genes were significantly differentiated between groups, as was Cg_CLCN7. For the hypersalinity adaptation, the biological process for positive regulating the developmental process was enriched. Enriched gene functions were focused on transcriptional regulation, signal transduction, and cell growth and differentiation, including calmodulin (Cg_CaM) and ficolin-2 (Cg_FCN2). These genes and polymorphisms possibly play an important role in oyster hyposalinity and hypersalinity adaptation. They not only further our understanding of salinity adaptation mechanisms but also provide markers for highly adaptable oyster strains suitable for breeding.

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.

Transcriptome profiling of claw muscle of the mud crab (Scylla paramamosain) at different fattening stages

In crustaceans, muscle growth and development is complicated, and to date substantial knowledge gaps exist. In this study, the claw muscle, hepatopancreas and nervous tissue of the mud crab (Scylla paramamosain) were collected at three fattening stages for sequence by the Illumina sequencing. A total of 127.87 Gb clean data with no less than 3.94 Gb generated for each sample and the cycleQ30 percentages were more than 86.13% for all samples. De Bruijn assembly of these clean data produced 94,853 unigenes, thereinto, 50,059 unigenes were found in claw muscle. A total of 121 differentially expressed genes (DEGs) were revealed in claw muscle from the three fattening stages with a Padj value < 0.01, including 63 genes with annotation. Functional annotation and enrichment analysis showed that the DEGs clusters represented the predominant gene catalog with roles in biochemical processes (glycolysis, phosphorylation and regulation of transcription), molecular function (ATP binding, 6-phosphofructokinase activity, and sequence-specific DNA binding) and cellular component (6-phosphofructokinase complex, plasma membrane, and integral component of membrane). qRT-PCR was employed to further validate certain DEGs. Single nucleotide polymorphism (SNP) analysis obtained 159,322, 125,963 and 166,279 potential SNPs from the muscle transcriptome at stage B, stage C and stage D, respectively. In addition, there were sixteen neuropeptide transcripts being predicted in the claw muscle. The present study provides a comprehensive transcriptome of claw muscle of S. paramamosain during fattening, providing a basis for screening the functional genes that may affect muscle growth of S. paramamosain.

Differential gene expression and SNP association between fast- and slow-growing turbot (Scophthalmus maximus)

Growth is among the most important traits for animal breeding. Understanding the mechanisms underlying growth differences between individuals can contribute to improving growth rates through more efficient breeding schemes. Here, we report a transcriptomic study in muscle and brain of fast- and slow-growing turbot (Scophthalmus maximus), a relevant flatfish in European and Asian aquaculture. Gene expression and allelic association between the two groups were explored. Up-regulation of the anaerobic glycolytic pathway in the muscle of fast-growing fish was observed, indicating a higher metabolic rate of white muscle. Brain expression differences were smaller and not associated with major growth-related genes, but with regulation of feeding-related sensory pathways. Further, SNP variants showing frequency differences between fast- and slow-growing fish pointed to genomic regions likely involved in growth regulation, and three of them were individually validated through SNP typing. Although different mechanisms appear to explain growth differences among families, general mechanisms seem also to be involved, and thus, results provide a set of useful candidate genes and markers to be evaluated for more efficient growth breeding programs and to perform comparative genomic studies of growth in fish and vertebrates.

Transcriptome analysis of the pectoral muscles of local chickens and commercial broilers using Ribo-Zero ribonucleic acid sequencing

BACKGROUND

The molecular mechanisms underlying meat quality and muscle growth are not clear. The meat quality and growth rates of local chickens and commercial broilers are very different. The Ribo-Zero RNA-Seq technology is an effective means of analyzing transcript groups to clarify molecular mechanisms. The aim of this study was to provide a reference for studies of the differences in the meat quality and growth of different breeds of chickens.

RESULTS

Ribo-Zero RNA-Seq technology was used to analyze the pectoral muscle transcriptomes of Gushi chickens and AA broilers. Compared with AA broilers, 1649 genes with annotated information were significantly differentially expressed (736 upregulated and 913 downregulated) in Gushi chickens with Q≤0.05 (Q is the P-value corrected by multiple assumptions test) at a fold change ≥2 or ≤0.5. In addition, 2540 novel significantly differentially expressed (SDE) genes (1405 upregulated and 1135 downregulated) were discovered. The results showed that the main signal transduction pathways that differed between Gushi chickens and AA broilers were related to amino acid metabolism. Amino acids are important for protein synthesis, and they regulate key metabolic pathways to improve the growth, development and reproduction of organisms.

CONCLUSION

This study showed that differentially expressed genes in the pectoral tissues of Gushi chickens and AA broilers were related to fat metabolism, which affects meat. Additionally, a large number of novel genes were found that may be involved in fat metabolism and thus may affect the formation of meat, which requires further study. The results of this study provide a reference for further studies of the molecular mechanisms of meat formation.

RNA-Seq reveals MicroRNA expression signature and genetic polymorphism associated with growth and muscle quality traits in rainbow trout

The role of microRNA expression and genetic variation in microRNA-binding sites of target genes on growth and muscle quality traits is poorly characterized. We used RNA-Seq approach to investigate their importance on 5 growth and muscle quality traits: whole body weight (WBW), muscle yield, muscle crude-fat content, muscle shear force and whiteness. Phenotypic data were collected from 471 fish, representing 98 families (~5 fish/family) from a growth-selected line. Muscle microRNAs and mRNAs were sequenced from 22 families showing divergent phenotypes. Ninety microRNAs showed differential expression between families with divergent phenotypes, and their expression was strongly associated with variation in phenotypes. A total of 204 single nucleotide polymorphisms (SNPs) present in 3′ UTR of target genes either destroyed or created novel illegitimate microRNA target sites; of them, 78 SNPs explained significant variation in the aforementioned 5 muscle traits. Majority of the phenotype-associated SNPs were present in microRNA-binding sites of genes involved in energy metabolism and muscle structure. These findings suggest that variation in microRNA expression and/or sequence variation in microRNA binding sites in target genes play an important role in mediating differences in fish growth and muscle quality phenotypes.

Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout

Background

Coding/functional SNPs change the biological function of a gene and, therefore, could serve as “large-effect” genetic markers. In this study, we used two bioinformatics pipelines, GATK and SAMtools, for discovering coding/functional SNPs with allelic-imbalances associated with total body weight, muscle yield, muscle fat content, shear force, and whiteness. Phenotypic data were collected for approximately 500 fish, representing 98 families (5 fish/family), from a growth-selected line, and the muscle transcriptome was sequenced from 22 families with divergent phenotypes (4 low- versus 4 high-ranked families per trait).

Results

GATK detected 59,112 putative SNPs; of these SNPs, 4798 showed allelic imbalances (>2.0 as an amplification and <0.5 as loss of heterozygosity). SAMtools detected 87,066 putative SNPs; and of them, 4962 had allelic imbalances between the low- and high-ranked families. Only 1829 SNPs with allelic imbalances were common between the two datasets, indicating significant differences in algorithms. The two datasets contained 7930 non-redundant SNPs of which 4439 mapped to 1498 protein-coding genes (with 6.4% non-synonymous SNPs) and 684 mapped to 295 lncRNAs. Validation of a subset of 92 SNPs revealed 1) 86.7–93.8% success rate in calling polymorphic SNPs and 2) 95.4% consistent matching between DNA and cDNA genotypes indicating a high rate of identifying SNPs with allelic imbalances. In addition, 4.64% SNPs revealed random monoallelic expression. Genome distribution of the SNPs with allelic imbalances exhibited high density for all five traits in several chromosomes, especially chromosome 9, 20 and 28. Most of the SNP-harboring genes were assigned to important growth-related metabolic pathways.

Conclusion

These results demonstrate utility of RNA-Seq in assessing phenotype-associated allelic imbalances in pooled RNA-Seq samples. The SNPs identified in this study were included in a new SNP-Chip design (available from Affymetrix) for genomic and genetic analyses in rainbow trout.

Electronic supplementary material

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

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.)

Background

Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish.

Results

In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions.

Conclusions

These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.

Electronic supplementary material

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

Blood transcriptomes and de novo identification of candidate loci for mating success in lekking great snipe (Gallinago media)

We assembled the great snipe blood transcriptome using data from fourteen lekking males, in order to de novo identify candidate genes related to sexual selection, and determined the expression profiles in relation to mating success. The three most highly transcribed genes were encoding different haemoglobin subunits. All tended to be overexpressed in males with high mating success. We also called single nucleotide polymorphisms (SNPs) from the transcriptome data and found considerable genetic variation for many genes expressed during lekking. Among these, we identified 14 polymorphic candidate SNPs that had a significant genotypic association with mating success (number of females mated with) and/or mating status (mated or not). Four of the candidate SNPs were found in HBAA (encoding the haemoglobin α-chain). Heterozygotes for one of these and one SNP in the gene PABPC1 appeared to enjoy higher mating success compared to males homozygous for either of the alleles. In a larger data set of individuals, we genotyped 38 of the identified SNPs but found low support for consistent selection as only one of the zygosities of previously identified candidate SNPs and none of their genotypes were associated with mating status. However, candidate SNPs generally showed lower levels of spatial genetic structure compared to noncandidate markers. We also scored the prevalence of avian malaria in a subsample of birds. Males infected with avian malaria parasites had lower mating success in the year of sampling than noninfected males. Parasite infection and its interaction with specific genes may thus affect performance on the lek.

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.

Transcriptome assembly and identification of genes and SNPs associated with growth traits in largemouth bass (Micropterus salmoides)

Growth is one of the most crucial economic traits of all aquaculture species, but the molecular mechanisms involved in growth of largemouth bass (Micropterus salmoides) are poorly understood. The objective of this study was to screen growth-related genes of M. salmoides by RNA sequencing and identify growth-related single-nucleotide polymorphism (SNP) markers through a growth association study. The muscle transcriptomes of fast- and slow-growing largemouth bass were obtained using the RNA-Seq technique. A total of 54,058,178 and 54,742,444 qualified Illumina read pairs were obtained for the fast-growing and slow-growing groups, respectively, giving rise to 4,865,236,020 and 4,926,819,960 total clean bases, respectively. Gene expression profiling showed that 3,530 unigenes were differentially expressed between the fast-growing and slow-growing phenotypes (false discovery rate ≤0.001, the absolute value of log₂ (fold change) ≥1), including 1,441 up-regulated and 2,889 down-regulated unigenes in the fast-growing largemouth bass. Analysis of these genes revealed that several signalling pathways, including the growth hormone-insulin-like growth factor 1 axis and signalling pathway, the glycolysis pathway, and the myostatin/transforming growth factor beta signalling pathway, as well as heat shock protein, cytoskeleton, and myofibril component genes might be associated with muscle growth. From these genes, 10 genes with putative SNPs were selected, and 17 SNPs were genotyped successfully. Marker-trait analysis in 340 individuals of Youlu No. 1 largemouth bass revealed three SNPs associated with growth in key genes (phosphoenolpyruvate carboxykinase 1, FOXO3b, and heat shock protein beta-1). This research provides information about key genes and SNPs related to growth, providing new clues to understanding the molecular basis of largemouth bass growth.

Transcriptome Analysis and Differential Gene Expression on the Testis of Orange Mud Crab, Scylla olivacea, during Sexual Maturation

Adequate genetic information is essential for sustainable crustacean fisheries and aquaculture management. The commercially important orange mud crab, Scylla olivacea, is prevalent in Southeast Asia region and is highly sought after. Although it is a suitable aquaculture candidate, full domestication of this species is hampered by the lack of knowledge about the sexual maturation process and the molecular mechanisms behind it, especially in males. To date, data on its whole genome is yet to be reported for S. olivacea. The available transcriptome data published previously on this species focus primarily on females and the role of central nervous system in reproductive development. De novo transcriptome sequencing for the testes of S. olivacea from immature, maturing and mature stages were performed. A total of approximately 144 million high-quality reads were generated and de novo assembled into 160,569 transcripts with a total length of 142.2 Mb. Approximately 15–23% of the total assembled transcripts were annotated when compared to public protein sequence databases (i.e. UniProt database, Interpro database, Pfam database and Drosophila melanogaster protein database), and GO-categorised with GO Ontology terms. A total of 156,181 high-quality Single-Nucleotide Polymorphisms (SNPs) were mined from the transcriptome data of present study. Transcriptome comparison among the testes of different maturation stages revealed one gene (beta crystallin like gene) with the most significant differential expression—up-regulated in immature stage and down-regulated in maturing and mature stages. This was further validated by qRT-PCR. In conclusion, a comprehensive transcriptome of the testis of orange mud crabs from different maturation stages were obtained. This report provides an invaluable resource for enhancing our understanding of this species’ genome structure and biology, as expressed and controlled by their gonads.

Searching for an Accurate Marker-Based Prediction of an Individual Quantitative Trait in Molecular Plant Breeding

Molecular plant breeding with the aid of molecular markers has played an important role in modern plant breeding over the last two decades. Many marker-based predictions for quantitative traits have been made to enhance parental selection, but the trait prediction accuracy remains generally low, even with the aid of dense, genome-wide SNP markers. To search for more accurate trait-specific prediction with informative SNP markers, we conducted a literature review on the prediction issues in molecular plant breeding and on the applicability of an RNA-Seq technique for developing function-associated specific trait (FAST) SNP markers. To understand whether and how FAST SNP markers could enhance trait prediction, we also performed a theoretical reasoning on the effectiveness of these markers in a trait-specific prediction, and verified the reasoning through computer simulation. To the end, the search yielded an alternative to regular genomic selection with FAST SNP markers that could be explored to achieve more accurate trait-specific prediction. Continuous search for better alternatives is encouraged to enhance marker-based predictions for an individual quantitative trait in molecular plant breeding.