RNA sequencing to study gene expression and SNP variations associated with growth in zebrafish fed a plant protein-based diet

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.

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.

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.

On the Compatibility of Fish Meal Replacements in Aquafeeds for Rainbow Trout. A Combined Metabolomic, Proteomic and Histological Study

The sustainable development of modern aquaculture must rely on a significant reduction of the fish meal (FM) used in aquafeed formulations. However, FM substitution with alternative ingredients in diets for carnivorous fish species often showed reduced nutrient absorption, significantly perturbed metabolisms, and histological changes at both hepatic and intestinal levels. In the present study, rainbow trout (Oncorhynchus mykiss) were fed three different experimental aquafeeds. A control diet with higher FM content (27.3%) than two test formulations in which FM was substituted with two more sustainable and promising alternatives: insect meal (Hermetia illucens larvae = 10.1%, FM = 11.6%) and poultry by-products meal (PBM = 14.8%; FM = 11.7%). Combined metabolomics and proteomics analyses of fish liver, together with histological examination of liver and intestine demonstrated that a well-balanced formulation of nutrients in the three diets allowed high metabolic compatibility of either substitution, paving the way for a deeper understanding of the impact of novel raw materials for the fish feed industry. Results show that the main metabolic pathways of nutrient absorption and catabolism were essentially unaltered by alternative feed ingredients, and also histological alterations were negligible. It is demonstrated that the substitution of FM with sustainable alternatives does not have a negative impact on fish metabolism, as long as the nutritional requirements of rainbow trout are fulfilled.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Identification and Expression Characterization of the Smad3 Gene and SNPs Associated with Growth Traits in the Hard Clam (Meretrix meretrix)

It has been demonstrated that the sekelsky mothers against decapentaplegic homolog 3 (Smad3) plays an important role in the growth and development of vertebrates. However, little is known about the association between the Smad3 gene and the growth traits of mollusks. In this study, Smad3 from the hard clam Meretrix meretrix (Mm-Smad3) was cloned, characterized, and screened for growth-related single nucleotide polymorphisms (SNPs) in its exons. The full-length cDNA of Mm-Smad3 was 1938 bp, encoding a protein with 428 amino acid residues. The protein sequence included an MH1 (27–135 aa) and MH2 domain (233–404 aa). Promoter analysis showed that the promoter sequence of Mm-Smad3 was 2548 bp, and the binding sites of Pit-1a, Antp, Hb, and other transcription factors are related to the growth and development of hard clams. The phylogenetic tree was divided into two major clusters, including mollusks and vertebrate. The expression level of Mm-Smad3 was predominantly detected in the mantle and foot, while extremely less expression was observed in the digestive gland. The low expression level of Mm-Smad3 was detected at the stages of unfertilized mature eggs, fertilized eggs, four-cell embryos, blastula, gastrulae, trochophore, and D-shaped larvae, whereas an opposite trend was observed regarding the highest expression at the umbo larvae stage (p < 0.05). In the mantle repair experiment, the time-course expression profiles showed that compared to the expression level at 0 h, Mm-Smad3 significantly decreased at 6 h (p < 0.05) but increased at 12 and 48 h. Further, the association analysis identified 11 SNPs in the exons of Mm-Smad3, of which three loci (c.597 C > T, c.660 C > T, c.792 A > T) were significantly related to the growth traits of clam (p < 0.05). Overall, our findings indicated that Mm-Smad3 is a growth-related gene and the detected SNP sites provide growth-related markers for molecular marker-assisted breeding of this species.

Comparative transcriptome analyses of the liver between Xenocypris microlepis and Xenocypris davidi under low copper exposure

Copper is one of the most ubiquitous environmental pollutants worldwide. Previous studies have focused on the toxicology of high copper exposure, while there has been comparatively less research on the biological effects of low copper exposure. Low concentrations of copper often exist in freshwater ecosystems, and its impact on the fish is unclear. Both Xenocypris microlepis and Xenocypris davidi are bottom-feeding fishes widely distributed in freshwater ecosystems of China, and they are more likely to be contaminated by low concentrations of copper. Low copper exposure may have effects on molecular regulation at the level of gene expression in the two Xenocypris species. To investigate gene expression differences involved in the response to low copper concentrations between X. microlepis and X. davidi, we established the responses to low copper exposure of 0.01 mg/L for 14 days at the transcriptional level, and RNA-Seq was used to perform a comparative transcriptomic analysis of the liver. A total of 74,135 and 60,894 unigenes from X. microlepis and X. davidi were assembled by transcriptome profiling, respectively. Among these, 84 genes of X. microlepis and 165 genes of X. davidi were identified as differentially expressed genes (DEGs). There were 60 and 135 up-regulated, 24 and 30 down-regulated genes in the two species, respectively. Comparative transcriptome analyses identified five differentially co-expressed genes (DCGs) related to low copper exposure from the DEGs of the two Xenocypris species. The five DCGs were related to the fishes' growth, antioxidant system, immune system and heavy metal tolerance. The results could help us to understand the molecular mechanisms of the response to low copper exposure, and the data should provide a valuable transcriptomic resource for the genus Xenocypris.

Intestinal Transcriptome Analysis Reveals Enrichment of Genes Associated with Immune and Lipid Mechanisms, Favoring Soybean Meal Tolerance in High-Growth Zebrafish (Danio Rerio)

The molecular mechanisms underlying fish tolerance to soybean meal (SBM) remain unclear. Identifying these mechanisms would be beneficial, as this trait favors growth. Two fish replicates from 19 experimental families were fed fishmeal-(100FM) or SBM-based diets supplemented with saponin (50SBM + 2SPN) from juvenile to adult stages. Individuals were selected from families with a genotype-by-environment interaction higher (HG-50SBM + 2SPN, 170 ± 18 mg) or lower (LG-50SBM + 2SPN, 76 ± 10 mg) weight gain on 50SBM + 2SPN for intestinal transcriptomic analysis. A histological evaluation confirmed middle intestinal inflammation in the LG- vs. HG-50SBM + 2SPN group. Enrichment analysis of 665 differentially expressed genes (DEGs) identified pathways associated with immunity and lipid metabolism. Genes linked to intestinal immunity were downregulated in HG fish (mpx, cxcr3.2, cftr, irg1l, itln2, sgk1, nup61l, il22), likely dampening inflammatory responses. Conversely, genes involved in retinol signaling were upregulated (rbp4, stra6, nr2f5), potentially favoring growth by suppressing insulin responses. Genes associated with lipid metabolism were upregulated, including key components of the SREBP (mbtps1, elov5l, elov6l) and cholesterol catabolism (cyp46a1), as well as the downregulation of cyp7a1. These results strongly suggest that transcriptomic changes in lipid metabolism mediate SBM tolerance. Genotypic variations in DEGs may become biomarkers for improving early selection of fish tolerant to SMB or others plant-based diets.

Immunologic Profiling of the Atlantic Salmon Gill by Single Nuclei Transcriptomics

Anadromous salmonids begin life adapted to the freshwater environments of their natal streams before a developmental transition, known as smoltification, transforms them into marine-adapted fish. In the wild, smoltification is a photoperiod-regulated process, involving radical remodeling of gill function to cope with the profound osmotic and immunological challenges of seawater (SW) migration. While prior work has highlighted the role of specialized "mitochondrion-rich" cells (MRCs) and accessory cells (ACs) in delivering this phenotype, recent RNA profiling experiments suggest that remodeling is far more extensive than previously appreciated. Here, we use single-nuclei RNAseq to characterize the extent of cytological changes in the gill of Atlantic salmon during smoltification and SW transfer. We identify 20 distinct cell clusters, including known, but also novel gill cell types. These data allow us to isolate cluster-specific, smoltification-associated changes in gene expression and to describe how the cellular make-up of the gill changes through smoltification. As expected, we noted an increase in the proportion of seawater mitochondrion-rich cells, however, we also identify previously unknown reduction of several immune-related cell types. Overall, our results provide fresh detail of the cellular complexity in the gill and suggest that smoltification triggers unexpected immune reprogramming.

Changes in testicular gene expression following reduced estradiol synthesis: A complex pathway to increased porcine Sertoli cell proliferation

Porcine Sertoli cell number including number present at puberty is increased if testicular estradiol synthesis is reduced during the neonatal interval. Evaluating the changes in gene expression during the crucial interval of suppressed estradiol that leads to the increased Sertoli cell population will increase our understanding of Sertoli cell biology but this evaluation first required a more precise determination of the critical interval for treatment and timing of a detectable response. Previously, reduced testicular estrogens from 1 week of age were accompanied by increased Sertoli cell number at 6.5 weeks of age but the age at which Sertoli cell numbers were initially increased was unknown, one of the current objectives. Additional experiments were designed to further delineate the essential timing of treatment for the Sertoli cell response. Finally, changes in gene expression induced by the reduced estradiol synthesis were evaluated to elucidate molecular mechanisms. Experimental design typically consisted of one member of littermate pairs of boars treated with the aromatase inhibitor, letrozole, beginning at 1 week of age and the remaining member treated with canola oil vehicle. Weekly treatments continued through 5 weeks of age or tissue collection, whichever came first. Increases in Sertoli cell numbers were not detectable prior to 6.5 weeks of age and persistent treatment through 5 weeks of age was required to induce the increase in Sertoli cell numbers. This increase resulted from prolonging the first interval of Sertoli cell proliferation in the treated animals. Few genes exhibited dramatically altered transcription and similarities in pathway analysis or principal modified genes were quite limited in 2, 3, and 5-week-old boars. The critical timing and prolonged treatment required and the sequential changes in gene expression suggest a complex mechanism is involved in this model of increased proliferation of Sertoli cells.

Plant-Based Diets Induce Transcriptomic Changes in Muscle of Zebrafish and Atlantic Salmon

With the expansion of the aquaculture industry in the last two decades, there has been a large increase in the use of plant ingredients in aquafeeds, which has created new challenges in fish growth, health and welfare. Fish muscle growth is an important trait that is strongly affected by diet, but our knowledge on the effect of plant protein-based diets on global gene expression in muscle is still scant. The present study evaluated nutrigenomic effects of the inclusion of proteins from pea, soy and wheat into aquafeeds, compared to a control diet with fishmeal as the main protein source using the zebrafish model by RNA-seq; these results were extended to an important aquaculture species by analyzing selected differentially expressed genes identified in the zebrafish model on on-growing Atlantic salmon fed with equivalent plant protein-based diets. Expression of selected Atlantic salmon paralogues of the zebrafish homologs was analyzed using paralogue-specific qPCR assays. Global gene expression changes in muscle of zebrafish fed with plant-based diets were moderate, with the highest changes observed in the soy diet-fed fish, and no change for the pea diet-fed fish compared to the control diet. Among the differentially expressed genes were mylpfb, hsp90aa1.1, col2a1a, and odc1, which are important in regulating muscle growth, maintaining muscle structure and function, and muscle tissue homeostasis. Furthermore, those genes and their paralogues were differentially expressed in Atlantic salmon fed with the equivalent percentage of soy or wheat protein containing diets. Some of these genes were similarly regulated in both species while others showed species-specific regulation. The present study expands our understanding on the molecular effects of plant ingredients in fish muscle. Ultimately, the knowledge gained would be of importance for the improved formulation of sustainable plant-based diets for the aquaculture industry.

Current basis and future directions of zebrafish nutrigenomics

This review investigates the current state of nutrigenomics in the zebrafish animal models. The zebrafish animal model has been used extensively in the study of disease onset and progression and associated molecular changes. In this review, we provide a synopsis of nutrigenomics using the zebrafish animal model. Obesity and dyslipidemia studies describe the genomics of dietary-induced obesity in relation to high-fat/high-calorie diets. Inflammation and cardiovascular studies describe dietary effects on the expression of acute inflammatory markers and resulting chronic inflammatory issues including atherosclerosis. We also evaluated the genomic response to bioactive dietary compounds associated with metabolic disorders. Carbohydrate metabolism and β-cell function studies describe the impacts of high-carbohydrate dietary challenges on nutritional programming. We also report tumorigenesis in relation to dietary carcinogen exposure studies that can result in permanent genomic changes. Vitamin and mineral deficiency studies demonstrate transgenerational genomic impacts of micronutrients in the diet and temporal expression changes. Circadian rhythm studies describe the relation between metabolism and natural temporal cycles of gene expression that impacts health. Bone formation studies describe the role of dietary composition that influences bone reabsorption regulation. Finally, this review provides future directions in the use of the zebrafish model for nutrigenomic and nutrigenetic research.

Serologic, but Not Genetic, Markers Are Associated With Impaired Anthropometrics at Diagnosis of Pediatric Crohn's Disease

OBJECTIVES

Children with Crohn's disease may present with malnutrition and linear growth impairment, which can be secondary to insufficient caloric intake, chronic inflammation, malabsorption, and suppression of growth-promoting hormones. We evaluated clinical, serologic, and genetic data to determine risk factors for impaired anthropometrics in Crohn's disease at diagnosis.

METHODS

Our study evaluated 772 children newly diagnosed with Crohn's disease, inflammatory phenotype, enrolled in the RISK Stratification Project to determine the factors associated with anthropometric impairment. Data were collected on demographics, growth parameters, disease location, serologic and immunologic markers, and disease severity. We performed a genome-wide association study of genetic polymorphisms associated with inflammatory bowel disease. Regression analysis determined associations between anthropometrics and clinical, serologic, and genetic variables.

RESULTS

There were 59 (7%) children with height z score <-2, 126 (14%) with a weight z score <-2, and 156 (17%) with a body mass index z score <-2. Linear growth impairment was associated with hypoalbuminemia (P = 0.0052), elevated granulocyte-macrophage colony stimulating factor autoantibodies (P = 0.0110), and elevated CBir antibodies against flagellin (P = 0.0117). Poor weight gain was associated with female sex (P = 0.0401), hypoalbuminemia (P = 0.0162), and thrombocytosis (P = 0.0081). Malnutrition was associated with hypoalbuminemia (P = 0.0061) and thrombocytosis (P = 0.0011). Children with moderate or severe disease had lower weight (P = 0.02 and P = 1.16×10, respectively) and body mass index z scores (P = 2.7 × 10 and P = 1.01 × 10, respectively) than children with quiescent and mild disease. There was no association between age of diagnosis, Tanner stage, or disease location and having impaired anthropometrics. There was no genome-wide association between the genetic polymorphisms and the serologic variables and anthropometric measurements.

CONCLUSIONS

This is the largest study evaluating growth in treatment-naïve children with Crohn's disease, inflammatory phenotype. It is the first study to use genome-wide sequencing to assess for genetic determinants of growth impairment. Granulocyte-macrophage colony stimulating factor autoantibodies and CBir antibodies are more likely to be elevated in children with growth impairment. Future investigations should evaluate the relationship between genetic polymorphisms, pathologic immune responses, and the biological pathways regulating growth.

Identification of Single Nucleotide Polymorphism in Red Clover (Trifolium pratense L.) Using Targeted Genomic Amplicon Sequencing and RNA-seq

Red clover (Trifolium pratense L.) is a diploid, naturally cross-pollinated, cool-season species. As a perennial forage legume, red clover is mostly cultivated in temperate regions worldwide. Being a non-model crop species, genomic resources for red clover have been underdeveloped. Thus far, genomic analysis used in red clover has mainly relied on simple sequence repeat (SSR) markers. However, SSR markers are sparse in the genome and it is often difficult to unambiguously map them using short reads generated by next generation sequencing technology. Single nucleotide polymorphisms (SNPs) have been successfully applied in genomics assisted breeding in several agriculturally important species. Due to increasing importance of legumes in forage production, there is a clear need to develop SNP based markers for red clover that can be applied in breeding applications. In this study, we first developed an analytical pipeline that can confidently identify SNPs in a set of 72 different red clover genotypes using sequences generated by targeted amplicon sequencing. Then, with the same filtering stringency used in this pipeline, we used sequences from publicly available RNA-seq data to identify confident SNPs in different red clover varieties. Using this strategy, we have identified a total of 69,975 SNPs across red clover varieties. Among these, 28% (19,116) of them are missense mutations. Using Medicago truncatula as the reference, we annotated the regions affected by these missense mutations. We identified 2,909 protein coding regions with missense mutations. Pathway analysis of these coding regions indicated several biological processes impacted by these mutations. Specifically, three domains (homeobox domain, pentatricopeptide repeat containing plant-like, and regulator of Vps4 activity) were identified with five or more missense SNPs. These domain might also be a functional contributor in the molecular mechanisms of self-incompatibility in red clover. Future in-depth sequence diversity analysis of these three genes may yield valuable insights into the molecular mechanism involved in self-incompatibility in red clover.

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.

Integrative functional analyses using rainbow trout selected for tolerance to plant diets reveal nutrigenomic signatures for soy utilization without the concurrence of enteritis

Finding suitable alternative protein sources for diets of carnivorous fish species remains a major concern for sustainable aquaculture. Through genetic selection, we created a strain of rainbow trout that outperforms parental lines in utilizing an all-plant protein diet and does not develop enteritis in the distal intestine, as is typical with salmonids on long-term plant protein-based feeds. By incorporating this strain into functional analyses, we set out to determine which genes are critical to plant protein utilization in the absence of gut inflammation. After a 12-week feeding trial with our selected strain and a control trout strain fed either a fishmeal-based diet or an all-plant protein diet, high-throughput RNA sequencing was completed on both liver and muscle tissues. Differential gene expression analyses, weighted correlation network analyses and further functional characterization were performed. A strain-by-diet design revealed differential expression ranging from a few dozen to over one thousand genes among the various comparisons and tissues. Major gene ontology groups identified between comparisons included those encompassing central, intermediary and foreign molecule metabolism, associated biosynthetic pathways as well as immunity. A systems approach indicated that genes involved in purine metabolism were highly perturbed. Systems analysis among the tissues tested further suggests the interplay between selection for growth, dietary utilization and protein tolerance may also have implications for nonspecific immunity. By combining data from differential gene expression and co-expression networks using selected trout, along with ontology and pathway analyses, a set of 63 candidate genes for plant diet tolerance was found. Risk loci in human inflammatory bowel diseases were also found in our datasets, indicating rainbow trout selected for plant-diet tolerance may have added utility as a potential biomedical model.

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.

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.

Transcriptomic Analysis of Compromise Between Air-Breathing and Nutrient Uptake of Posterior Intestine in Loach (Misgurnus anguillicaudatus), an Air-Breathing Fish

Dojo loach (Misgurnus anguillicaudatus) is an air-breathing fish species by using its posterior intestine to breathe on water surface. So far, the molecular mechanism about accessory air-breathing in fish is seldom addressed. Five cDNA libraries were constructed here for loach posterior intestines form T01 (the initial stage group), T02 (mid-stage of normal group), T03 (end stage of normal group), T04 (mid-stage of air-breathing inhibited group), and T05 (the end stage of air-breathing inhibited group) and subjected to perform RNA-seq to compare their transcriptomic profilings. A total of 92,962 unigenes were assembled, while 37,905 (40.77 %) unigenes were successfully annotated. 2298, 1091, and 3275 differentially expressed genes (fn1, ACE, EGFR, Pxdn, SDF, HIF, VEGF, SLC2A1, SLC5A8 etc.) were observed in T04/T02, T05/T03, and T05/T04, respectively. Expression levels of many genes associated with air-breathing and nutrient uptake varied significantly between normal and intestinal air-breathing inhibited group. Intraepithelial capillaries in posterior intestines of loaches from T05 were broken, while red blood cells were enriched at the surface of intestinal epithelial lining with 241 ± 39 cells per millimeter. There were periodic acid-schiff (PAS)-positive epithelial mucous cells in posterior intestines from both normal and air-breathing inhibited groups. Results obtained here suggested an overlap of air-breathing and nutrient uptake function of posterior intestine in loach. Intestinal air-breathing inhibition in loach would influence the posterior intestine's nutrient uptake ability and endothelial capillary structure stability. This study will contribute to our understanding on the molecular regulatory mechanisms of intestinal air-breathing in loach.

Identification of Immune-Related Genes and Development of SSR/SNP Markers from the Spleen Transcriptome of Schizothorax prenanti

Schizothorax prenanti (S. prenanti) is mainly distributed in the upstream regions of the Yangtze River and its tributaries in China. This species is indigenous and commercially important. However, in recent years, wild populations and aquacultures have faced the serious challenges of germplasm variation loss and an increased susceptibility to a range of pathogens. Currently, the genetics and immune mechanisms of S. prenanti are unknown, partly due to a lack of genome and transcriptome information. Here, we sought to identify genes related to immune functions and to identify molecular markers to study the function of these genes and for trait mapping. To this end, the transcriptome from spleen tissues of S. prenanti was analyzed and sequenced. Using paired-end reads from the Illumina Hiseq2500 platform, 48,517 transcripts were isolated from the spleen transcriptome. These transcripts could be clustered into 37,785 unigenes with an N50 length of 2,539 bp. The majority of the unigenes (35,653, 94.4%) were successfully annotated using non-redundant nucleotide sequence analysis (nt), and the non-redundant protein (nr), Swiss-Prot, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. KEGG pathway assignment identified more than 500 immune-related genes. Furthermore, 7,545 putative simple sequence repeats (SSRs), 857,535 single nucleotide polymorphisms (SNPs), and 53,481 insertion/deletion (InDels) were detected from the transcriptome. This is the first reported high-throughput transcriptome analysis of S. prenanti, and it provides valuable genetic resources for the investigation of immune mechanisms, conservation of germplasm, and molecular marker-assisted breeding of S. prenanti.