Ovary transcriptome profiling via artificial intelligence reveals a transcriptomic fingerprint predicting egg quality in striped bass, Morone saxatilis

Combinatorial metabolomic and transcriptomic analysis of muscle growth in hybrid striped bass (female white bass Morone chrysops x male striped bass M. saxatilis)

BACKGROUND

Understanding growth regulatory pathways is important in aquaculture, fisheries, and vertebrate physiology generally. Machine learning pattern recognition and sensitivity analysis were employed to examine metabolomic small molecule profiles and transcriptomic gene expression data generated from liver and white skeletal muscle of hybrid striped bass (white bass Morone chrysops x striped bass M. saxatilis) representative of the top and bottom 10 % by body size of a production cohort.

RESULTS

Larger fish (good-growth) had significantly greater weight, total length, hepatosomatic index, and specific growth rate compared to smaller fish (poor-growth) and also had significantly more muscle fibers of smaller diameter (≤ 20 µm diameter), indicating active hyperplasia. Differences in metabolomic pathways included enhanced energetics (glycolysis, citric acid cycle) and amino acid metabolism in good-growth fish, and enhanced stress, muscle inflammation (cortisol, eicosanoids) and dysfunctional liver cholesterol metabolism in poor-growth fish. The majority of gene transcripts identified as differentially expressed between groups were down-regulated in good-growth fish. Several molecules associated with important growth-regulatory pathways were up-regulated in muscle of fish that grew poorly: growth factors including agt and agtr2 (angiotensins), nicotinic acid (which stimulates growth hormone production), gadd45b, rgl1, zfp36, cebpb, and hmgb1; insulin-like growth factor signaling (igfbp1 and igf1); cytokine signaling (socs3, cxcr4); cell signaling (rgs13, rundc3a), and differentiation (rhou, mmp17, cd22, msi1); mitochondrial uncoupling proteins (ucp3, ucp2); and regulators of lipid metabolism (apoa1, ldlr). Growth factors pttg1, egfr, myc, notch1, and sirt1 were notably up-regulated in muscle of good-growing fish.

CONCLUSION

A combinatorial pathway analysis using metabolomic and transcriptomic data collectively suggested promotion of cell signaling, proliferation, and differentiation in muscle of good-growth fish, whereas muscle inflammation and apoptosis was observed in poor-growth fish, along with elevated cortisol (an anti-inflammatory hormone), perhaps related to muscle wasting, hypertrophy, and inferior growth. These findings provide important biomarkers and mechanisms by which growth is regulated in fishes and other vertebrates as well.

Seasonal modulation of the testis transcriptome reveals insights into hibernation and reproductive adaptation in Onychostoma macrolepis

The Onychostoma macrolepis have a unique survival strategy, overwintering in caves and returning to the river for reproduction in summer. The current knowledge on the developmental status of its testes during winter and summer is still undiscovered. We performed RNA-seq analysis on O. macrolepis testes between January and June, using the published genome (NCBI, ASM1243209v1). Through KEGG and GO enrichment analysis, we were able to identify 2111 differentially expressed genes (DEGs) and demonstrate their functions in signaling networks associated with the development of organism. At the genomic level, we found that during the overwintering phase, genes associated with cell proliferation (ccnb1, spag5, hdac7) were downregulated while genes linked to testicular fat metabolism (slc27a2, scd, pltp) were upregulated. This indicates suppression of both mitosis and meiosis, thereby inhibiting energy expenditure through genetic regulation of testicular degeneration. Furthermore, in January, we observed the regulation of autophagy and apoptosis (becn1, casp13), which may have the function of protecting reproductive organs and ensuring their maturity for the breeding season. The results provide a basis for the development of specialized feed formulations to regulate the expression of specific genes, or editing of genes during the fish egg stage, to ensure that the testes of O. macrolepis can mature more efficiently after overwintering, thereby enhancing reproductive performance.

Transcriptome analysis of Schizothorax oconnori (Cypriniformes: Cyprinidae) oocytes: The role of K+ in promoting yolk globule fusion and regulating oocyte maturation

Schizothorax oconnori (S. oconnori) is an economically important fish in Tibet. Oocyte maturation is a physiological process that is of great significance to reproduction and seed production in S. oconnori, yet little is currently known regarding the molecular mechanisms of oocyte development in this species. To identify candidate genes involved in reproduction of female fish, a combination of PacBio and Illumina HiSeq technologies was employed to provide deep coverage of the oocyte transcriptome. Transcriptome analysis revealed several candidate genes that are potentially involved in the regulation of oocyte maturation in S. oconnori, including GIRK1, CHRM3, NPY2R, GABRA3, GnRH3, mGluR1α, GPER1, GDF9, HSP90, and ESR2. Genes that are significantly expressed during oocyte maturation mainly contribute to the GPCR signaling pathway and the estrogen signaling pathway. Neurotransmitter (Ach, NPY, and GABA) and peptide hormone (GnRH3) binding to G protein-coupled receptors (GPCRs) frees G-protein βγ subunits to interact with the G protein-gated inward rectifier K+ channel 1 (GIRK1). This process helps release K+ from granulosa cells to maturing oocytes, allowing yolk globule fusion. This mechanism may play an important role in oocyte maturation in S. oconnori. In conclusion, this study provides a valuable basis for deciphering the reproductive system in S. oconnori during the oocyte maturation process.

Multiple Tissues Transcriptome of Zig-Zag Eel (Mastacembelus armatus) with Different Growth Rates

In order to explore the main regulatory genes and related pathways of growth traits, transcriptome sequencing was first performed on the brain, liver, and muscle tissues of 3-month-old M. armatus with different growth rates. By comparative transcriptome analysis of fast-growing and slow-growing groups of M. armatus, a total of 2887 DEGs were screened, of which 59 up-regulated genes and 105 down-regulated genes were detected in the brain, 146 up-regulated genes and 202 down-regulated genes were detected in the liver, and 529 up-regulated genes and 1846 down-regulated genes were detected in muscle, including insulin-like growth factor binding protein 1a (IGFBP1A), insulin-like growth factor binding protein 1b (IGFBP1B), myosin, light chain 1 (MYL1), and myoglobin (MB). Through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, we identified a total of 288 significantly enriched GO entries and 68 significantly enriched KEGG pathways related to growth, such as skeletal muscle tissue development, insulin-like growth factor binding, and the mitotic cell cycle. These key genes and signaling pathways may play a key role in regulating the growth of M. armatus. Digging into the regulatory mechanisms of these key genes will provide a theoretical basis for further exploration of the molecular mechanisms related to the growth and development of M. armatus, and help to breed new varieties of M. armatus with rapid growth.

Machine learning reveals common transcriptomic signatures across rat brain and placenta following developmental organophosphate ester exposure

Toxicogenomics is a critical area of inquiry for hazard identification and to identify both mechanisms of action and potential markers of exposure to toxic compounds. However, data generated by these experiments are highly dimensional and present challenges to standard statistical approaches, requiring strict correction for multiple comparisons. This stringency often fails to detect meaningful changes to low expression genes and/or eliminate genes with small but consistent changes particularly in tissues where slight changes in expression can have important functional differences, such as brain. Machine learning offers an alternative analytical approach for "omics" data that effectively sidesteps the challenges of analyzing highly dimensional data. Using 3 rat RNA transcriptome sets, we utilized an ensemble machine learning approach to predict developmental exposure to a mixture of organophosphate esters (OPEs) in brain (newborn cortex and day 10 hippocampus) and late gestation placenta of male and female rats, and identified genes that informed predictor performance. OPE exposure had sex specific effects on hippocampal transcriptome, and significantly impacted genes associated with mitochondrial transcriptional regulation and cation transport in females, including voltage-gated potassium and calcium channels and subunits. To establish if this holds for other tissues, RNAseq data from cortex and placenta, both previously published and analyzed via a more traditional pipeline, were reanalyzed with the ensemble machine learning methodology. Significant enrichment for pathways of oxidative phosphorylation and electron transport chain was found, suggesting a transcriptomic signature of OPE exposure impacting mitochondrial metabolism across tissue types and developmental epoch. Here we show how machine learning can complement more traditional analytical approaches to identify vulnerable "signature" pathways disrupted by chemical exposures and biomarkers of exposure.

Quality of fish eggs and production of androgenetic and gynogenetic doubled haploids (DHs)

Induced development of haploid embryos (H) with only paternal (androgenesis) or maternal (gynogenesis) chromosomes requires irradiation of eggs before fertilization or activation of eggs with irradiated spermatozoa, respectively. To provide doubled haploids (DHs), androgenetic and gynogenetic haploid zygotes need to be subjected to the thermal or high hydrostatic pressure (HHP) shock to suppress the first mitotic cleavage and to double paternal or maternal haploid set of chromosomes. Androgenesis and mitotic gynogenesis (mito-gynogenesis) result in the generation of fully homozygous individuals in a single generation. DHs have been utilized in selective breeding programs, in studies concerning the phenotypic consequences of recessive alleles and to evaluate the impact of sex chromosomes on the early ontogeny. Moreover, the use of DHs for the NGS approach radically improves de novo the assembly of the genomes. However, reduced survival of the doubled haploids limits the wide application of androgenotes and gynogenotes. The high mortality of DHs may be only partly explained by the expression of recessive traits. Observed inter-clutch variation in the survival of DHs developing in eggs originating from different females make it necessary to take a closer look at the quality of the eggs used during induced androgenesis and gynogenesis. Moreover, the developmental competence of eggs that are subjected to irradiation before fertilization in order to deactivate maternal chromosomes when undergoing induced androgenesis and exposed to the physical shock after fertilization that leads to the duplication of the zygotes in both mito-gynogenesis and androgenesis may be also altered as irradiation and sublethal values of temperatures and hydrostatic pressure are considered as harmful for the cell organelles and biomolecules. Here, recently provided results concerning the morphological, biochemical, genomic, and transcriptomic characteristics of fish eggs showing high and low competence for androgenesis and mito-gynogenesis are reviewed.

Combining supervised machine learning with statistics reveals differential gene expression patterns related to energy metabolism in the jejuna of chickens divergently selected for antibody response to sheep red blood cells

Since the 1970s, 2 lines of White Leghorn chickens, HAS and LAS, have been continuously divergently selected for 5-day postinjection antibody titer to injection with sheep red blood cells (SRBC). Antibody response is a complex genetic trait and characterizing differences in gene expression could facilitate better understanding of physiological changes due to selection and antigen exposure. At 41 d of age, randomly selected HAS and LAS chickens, which had been coraised from hatch, were either injected with SRBC (HASI and LASI) or kept as the noninjected cohort (HASN and LASN). Five days later, all were euthanized, and samples collected from the jejunum for RNA isolation and sequencing. Resulting gene expression data were analyzed combining traditional statistics with machine learning to obtain signature gene lists for functional analysis. Differences in ATP production and cellular processes were observed in the jejunum between lines and following SRBC injection. HASN vs. LASN exhibited upregulation of ATP production, immune cell motility, and inflammation. LASI exhibits upregulation of ATP production and protein synthesis vs. LASN, reflective of what was observed in HASN vs. LASN. In contrast, no corresponding upregulation of ATP production was observed in HASI vs. HASN, and most other cellular processes appear inhibited. Without exposure to SRBC, gene expression in the jejunum indicates HAS generates more ATP than LAS, suggesting HAS maintains a "primed" system; and gene expression of HASI vs. HASN further suggests this basal ATP production is sufficient for robust antibody responses. Conversely, LASI vs. LASN jejunal gene expression implies a physiological need for increased ATP production with only minimal correlating antibody production. The results of this experiment provide insight into energetic resource needs and allocations in the jejunum in response to genetic selection and antigen exposure in HAS and LAS which may help explain phenotypic differences observed in antibody response.

Differential bumble bee gene expression associated with pathogen infection and pollen diet

BACKGROUND

Diet and parasitism can have powerful effects on host gene expression. However, how specific dietary components affect host gene expression that could feed back to affect parasitism is relatively unexplored in many wild species. Recently, it was discovered that consumption of sunflower (Helianthus annuus) pollen reduced severity of gut protozoan pathogen Crithidia bombi infection in Bombus impatiens bumble bees. Despite the dramatic and consistent medicinal effect of sunflower pollen, very little is known about the mechanism(s) underlying this effect. However, sunflower pollen extract increases rather than suppresses C. bombi growth in vitro, suggesting that sunflower pollen reduces C. bombi infection indirectly via changes in the host. Here, we analyzed whole transcriptomes of B. impatiens workers to characterize the physiological response to sunflower pollen consumption and C. bombi infection to isolate the mechanisms underlying the medicinal effect. B. impatiens workers were inoculated with either C. bombi cells (infected) or a sham control (un-infected) and fed either sunflower or wildflower pollen ad libitum. Whole abdominal gene expression profiles were then sequenced with Illumina NextSeq 500 technology.

RESULTS

Among infected bees, sunflower pollen upregulated immune transcripts, including the anti-microbial peptide hymenoptaecin, Toll receptors and serine proteases. In both infected and un-infected bees, sunflower pollen upregulated putative detoxification transcripts and transcripts associated with the repair and maintenance of gut epithelial cells. Among wildflower-fed bees, infected bees downregulated immune transcripts associated with phagocytosis and the phenoloxidase cascade.

CONCLUSIONS

Taken together, these results indicate dissimilar immune responses between sunflower- and wildflower-fed bumble bees infected with C. bombi, a response to physical damage to gut epithelial cells caused by sunflower pollen, and a strong detoxification response to sunflower pollen consumption. Identifying host responses that drive the medicinal effect of sunflower pollen in infected bumble bees may broaden our understanding of plant-pollinator interactions and provide opportunities for effective management of bee pathogens.

Comparative transcriptome provides insights into differentially expressed genes between testis and ovary of Onychostoma macrolepis in reproduction period

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered fishery species inhabiting the river of Qinling Mountains and some flowing freshwaters in China. The declining population of O. macrolepis caused by asynchrony of male and female development prompted us to focus on genetic regulation of its reproduction. In this study, high-throughput RNA-sequencing technology was applied to assemble and annotate the transcriptome of O. macrolepis testis and ovary. The results showed that a number of 338089335 (ovary:163216500, testis:174872835) raw sequences were obtained. After non-redundant analysis, a number of 207826065 (ovary:102334008, testis:105492057) high quality reads were obtained and predicted as unigenes, in which 201,038,682 unigenes were annotated with multiple databases. Taking the ovarian transcriptome as a control, comparative transcriptome analysis showed that 9918 differentially expressed genes (DEGs) up-regulated in the testis and 13,095 DEGs down-regulated. Many DEGs were involved with sex-related GO terms and KEGG pathways, such as oocyte maturation, gonadal development, steroid biosynthesis pathways, MAPK signaling pathway and Wnt signaling pathway. Finally, the expression patterns of 19 unigenes were validated by using quantitative real-time polymerase chain reaction (qRT-PCR). This study illustrates a potential molecular mechanism on the unsynchronized male and female development of the O. macrolepis during the reproduction period in June and provides a theoretical basis for future artificial reproduction.

Quantitative proteome profiling reveals molecular hallmarks of egg quality in Atlantic halibut: impairments of transcription and protein folding impede protein and energy homeostasis during early development

Background

Tandem mass tag spectrometry (TMT labeling-LC-MS/MS) was utilized to examine the global proteomes of Atlantic halibut eggs at the 1-cell-stage post fertilization. Comparisons were made between eggs judged to be of good quality (GQ) versus poor quality (BQ) as evidenced by their subsequent rates of survival for 12 days. Altered abundance of selected proteins in BQ eggs was confirmed by parallel reaction monitoring spectrometry (PRM-LC-MS/MS). Correspondence of protein levels to expression of related gene transcripts was examined via qPCR. Potential mitochondrial differences between GQ and BQ eggs were assessed by transmission electron microscopy (TEM) and measurements of mitochondrial DNA (mtDNA) levels.

Results

A total of 115 proteins were found to be differentially abundant between GQ and BQ eggs. Frequency distributions of these proteins indicated higher protein folding activity in GQ eggs compared to higher transcription and protein degradation activities in BQ eggs. BQ eggs were also significantly enriched with proteins related to mitochondrial structure and biogenesis. Quantitative differences in abundance of several proteins with parallel differences in their transcript levels were confirmed in egg samples obtained over three consecutive reproductive seasons. The observed disparities in global proteome profiles suggest impairment of protein and energy homeostasis related to unfolded protein response and mitochondrial stress in BQ eggs. TEM revealed BQ eggs to contain significantly higher numbers of mitochondria, but differences in corresponding genomic mtDNA (mt-nd5 and mt-atp6) levels were not significant. Mitochondria from BQ eggs were significantly smaller with a more irregular shape and a higher number of cristae than those from GQ eggs.

Conclusion

The results of this study indicate that BQ Atlantic halibut eggs are impaired at both transcription and translation levels leading to endoplasmic reticulum and mitochondrial disorders. Observation of these irregularities over three consecutive reproductive seasons in BQ eggs from females of diverse background, age and reproductive experience indicates that they are a hallmark of poor egg quality. Additional research is needed to discover when in oogenesis and under what circumstances these defects may arise. The prevalence of this suite of markers in BQ eggs of diverse vertebrate species also begs investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08859-0.

Neurodevelopment vs. the immune system: Complementary contributions of maternally-inherited gene transcripts and proteins to successful embryonic development in fish

The aim of this study was to investigate the respective contribution of maternally-inherited mRNAs and proteins to egg molecular cargo and to its developmental competence in fish using pikeperch as a model. Our study provides novel insights into the understanding of type-specific roles of maternally-inherited molecules in fish. Here we show, for the first time, that transcripts and proteins have distinct, yet complementary, functions in the egg of teleost fish. Maternally-inherited mRNAs would shape embryo neurodevelopment, while maternally-inherited proteins would rather be responsible for protecting the embryo against pathogens. Additionally, we observed that processes directly preceding ovulation may considerably affect the reproductive success by modifying expression level of genes crucial for proper embryonic development, being novel fish egg quality markers (e.g., smarca4 or h3f3a). These results are of major importance for understanding the influence of external factors on reproductive fitness in both captive and wild-type fish species.

Comparisons among rainbow trout, Oncorhynchus mykiss, populations of maternal transcript profile associated with egg viability

Background

Transcription is arrested in the late stage oocyte and therefore the maternal transcriptome stored in the oocyte provides nearly all the mRNA required for oocyte maturation, fertilization, and early cleavage of the embryo. The transcriptome of the unfertilized egg, therefore, has potential to provide markers for predictors of egg quality and diagnosing problems with embryo production encountered by fish hatcheries. Although levels of specific transcripts have been shown to associate with measures of egg quality, these differentially expressed genes (DEGs) have not been consistent among studies. The present study compares differences in select transcripts among unfertilized rainbow trout eggs of different quality based on eyeing rate, among 2 year classes of the same line (A1, A2) and a population from a different hatchery (B). The study compared 65 transcripts previously reported to be differentially expressed with egg quality in rainbow trout.

Results

There were 32 transcripts identified as DEGs among the three groups by regression analysis. Group A1 had the most DEGs, 26; A2 had 15, 14 of which were shared with A1; and B had 12, 7 of which overlapped with A1 or A2. Six transcripts were found in all three groups, dcaf11, impa2, mrpl39_like, senp7, tfip11 and uchl1.

Conclusions

Our results confirmed maternal transcripts found to be differentially expressed between low- and high-quality eggs in one population of rainbow trout can often be found to overlap with DEGs in other populations. The transcripts differentially expressed with egg quality remain consistent among year classes of the same line. Greater similarity in dysregulated transcripts within year classes of the same line than among lines suggests patterns of transcriptome dysregulation may provide insight into causes of decreased viability within a hatchery population. Although many DEGs were identified, for each of the genes there is considerable variability in transcript abundance among eggs of similar quality and low correlations between transcript abundance and eyeing rate, making it highly improbable to predict the quality of a single batch of eggs based on transcript abundance of just a few genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07773-1.

Vision-based egg quality prediction in Pacific bluefin tuna (Thunnus orientalis) by deep neural network

Closed-cycle aquaculture using hatchery produced seed stocks is vital to the sustainability of endangered species such as Pacific bluefin tuna (Thunnus orientalis) because this aquaculture system does not depend on aquaculture seeds collected from the wild. High egg quality promotes efficient aquaculture production by improving hatch rates and subsequent growth and survival of hatched larvae. In this study, we investigate the possibility of a simple, low-cost, and accurate egg quality prediction system based only on photographic images using deep neural networks. We photographed individual eggs immediately after spawning and assessed their qualities, i.e., whether they hatched normally and how many days larvae survived without feeding. The proposed system predicted normally hatching eggs with higher accuracy than human experts. It was also successful in predicting which eggs would produce longer-surviving larvae. We also analyzed the image aspects that contributed to the prediction to discover important egg features. Our results suggest the applicability of deep learning techniques to efficient egg quality prediction, and analysis of early developmental stages of development.

The Molecular Machinery of Gametogenesis in Geodia Demosponges (Porifera): Evolutionary Origins of a Conserved Toolkit across Animals

All animals are capable of undergoing gametogenesis. The ability of forming haploid cells from diploid cells through meiosis and recombination appeared early in eukaryotes, whereas further gamete differentiation is mostly a metazoan signature. Morphologically, the gametogenic process presents many similarities across animal taxa, but little is known about its conservation at the molecular level. Porifera are the earliest divergent animals and therefore are an ideal phylum to understand evolution of the gametogenic toolkits. Although sponge gametogenesis is well known at the histological level, the molecular toolkits for gamete production are largely unknown. Our goal was to identify the genes and their expression levels which regulate oogenesis and spermatogenesis in five gonochoristic and oviparous species of the genus Geodia, using both RNAseq and proteomic analyses. In the early stages of both female and male gametogenesis, genes involved in germ cell fate and cell-renewal were upregulated. Then, molecular signals involved in retinoic acid pathway could trigger the meiotic processes. During later stages of oogenesis, female sponges expressed genes involved in cell growth, vitellogenesis, and extracellular matrix reassembly, which are conserved elements of oocyte maturation in Metazoa. Likewise, in spermatogenesis, genes regulating the whole meiotic cycle, chromatin compaction, and flagellum axoneme formation, that are common across Metazoa were overexpressed in the sponges. Finally, molecular signals possibly related to sperm capacitation were identified during late stages of spermatogenesis for the first time in Porifera. In conclusion, the activated molecular toolkit during gametogenesis in sponges was remarkably similar to that deployed during gametogenesis in vertebrates.

Distinctive Cellular and Metabolic Reprogramming in Porcine Lung Mononuclear Phagocytes Infected With Type 1 PRRSV Strains

Porcine reproductive and respiratory syndrome (PRRS) has an extensive impact on pig production. The causative virus (PRRSV) is divided into two species, PRRSV-1 (European origin) and PRRSV-2 (North American origin). Within PRRSV-1, PRRSV-1.3 strains, such as Lena, are more pathogenic than PRRSV-1.1 strains, such as Flanders 13 (FL13). To date, the molecular interactions of PRRSV with primary lung mononuclear phagocyte (MNP) subtypes, including conventional dendritic cells types 1 (cDC1) and 2 (cDC2), monocyte-derived DCs (moDC), and pulmonary intravascular macrophages (PIM), have not been thoroughly investigated. Here, we analyze the transcriptome profiles of in vivo FL13-infected parenchymal MNP subpopulations and of in vitro FL13- and Lena-infected parenchymal MNP. The cell-specific expression profiles of in vivo sorted cells correlated with their murine counterparts (AM, cDC1, cDC2, moDC) with the exception of PIM. Both in vivo and in vitro, FL13 infection altered the expression of a low number of host genes, and in vitro infection with Lena confirmed the higher ability of this strain to modulate host response. Machine learning (ML) and gene set enrichment analysis (GSEA) unraveled additional relevant genes and pathways modulated by FL13 infection that were not identified by conventional analyses. GSEA increased the cellular pathways enriched in the FL13 data set, but ML allowed a more complete comprehension of functional profiles during FL13 in vitro infection. Data indicates that cellular reprogramming differs upon Lena and FL13 infection and that the latter might keep antiviral and inflammatory macrophage/DC functions silent. Although the slow replication kinetics of FL13 likely contribute to differences in cellular gene expression, the data suggest distinct mechanisms of interaction of the two viruses with the innate immune system during early infection.

Fatty-acid Profiles and Fingerprints of Seven Types of Fish Roes as Determined by Chemometric Methods

The fatty acids in seven species of fish roes were determined by GC-MS in combination with principal component and cluster analyses in order to derive their fatty-acid profiles and fingerprints. Twenty-three common chromatography peaks were identified in the fatty-acid fingerprints of the seven fish roes. A total of 19 typical fatty acids were identified in the fish roes studied. The fatty acid contents of the roes were significantly different, with saturated-fatty-acid contents in the seven roes ranging from 26.69% to 41.81%, and the unsaturated-fatty-acid contents ranging from 57.65% to 72.21%, the total EPA and DHA content (37.20%) is high in E. cypselurus roe, especially. The seven roe species were clearly distinguished according to fatty-acid composition and content by principal component analysis (PCA) and divided into two groups by cluster analysis (CA). PCA of the fatty acid data yielded three significant PCs , which together account for 94% of the total variance; with PC1 contributing 54% of the total.

Biallelic variant in cyclin B3 is associated with failure of maternal meiosis II and recurrent digynic triploidy

Background

Triploidy is one of the most common chromosome abnormalities affecting human gestation and accounts for an important fraction of first-trimester miscarriages. Triploidy has been demonstrated in a few cases of recurrent pregnancy loss (RPL) but its molecular mechanisms are unknown. This study aims to identify the genetic cause of RPL associated with fetus triploidy.

Methods

We investigated genomic imprinting, genotyped sequence-tagged site (STS) markers and performed exome sequencing in a family including two sisters with RPL. Moreover, we evaluated oocyte maturation in vivo and in vitro and effect of the candidate protein variant in silico.

Results

While features of hydatidiform mole were excluded, the presence of triploidy of maternal origin was demonstrated in the fetuses. Oocyte maturation was deficient and all the maternally inherited pericentromeric STS alleles were homozygous in the fetuses. A deleterious missense variant (p.V1251D) of the cyclin B3 gene (CCNB3) affecting a residue conserved in placental mammals and located in a region that can interact with the cyclin-dependent kinase 1 or cyclin-dependent kinase 2 cosegregated in homozygosity with RPL.

Conclusion

Here, we report a family in which a damaging variant in cyclin B3 is associated with the failure of oocyte meiosis II and recurrent fetus triploidy, implicating a rationale for CCNB3 testing in RPL.

Differential impacts of carp and salmon pituitary extracts on induced oogenesis, egg quality, molecular ontogeny and embryonic developmental competence in European eel

Low egg quality and embryonic survival are critical challenges in aquaculture, where assisted reproduction procedures and other factors may impact egg quality. This includes European eel (Anguilla anguilla), where pituitary extract from carp (CPE) or salmon (SPE) is applied to override a dopaminergic inhibition of the neuroendocrine system, preventing gonadotropin secretion and gonadal development. The present study used either CPE or SPE to induce vitellogenesis in female European eel and compared impacts on egg quality and offspring developmental competence with emphasis on the maternal-to-zygotic transition (MZT). Females treated with SPE produced significantly higher proportions of floating eggs with fewer cleavage abnormalities and higher embryonic survival. These findings related successful embryogenesis to higher abundance of mRNA transcripts of genes involved in cell adhesion, activation of MZT, and immune response (dcbld1, epcam, oct4, igm) throughout embryonic development. The abundance of mRNA transcripts of cldnd, foxr1, cea, ccna1, ccnb1, ccnb2, zar1, oct4, and npm2 was relatively stable during the first eight hours, followed by a drop during MZT and low levels thereafter, indicating transfer and subsequent clearance of maternal mRNA. mRNA abundance of zar1, epcam, and dicer1 was associated with cleavage abnormalities, while mRNA abundance of zar1, sox2, foxr1, cldnd, phb2, neurod4, and neurog1 (before MZT) was associated with subsequent embryonic survival. In a second pattern, low initial mRNA abundance with an increase during MZT and higher levels persisting thereafter indicating the activation of zygotic transcription. mRNA abundance of ccna1, npm2, oct4, neurod4, and neurog1 during later embryonic development was associated with hatch success. A deviating pattern was observed for dcbld1, which mRNA levels followed the maternal-effect gene pattern but only for embryos from SPE treated females. Together, the differences in offspring production and performance reported in this study show that PE composition impacts egg quality and embryogenesis and in particular, the transition from initial maternal transcripts to zygotic transcription.

Domestication is associated with differential expression of pikeperch egg proteins involved in metabolism, immune response and protein folding

Domestication is a condition in which the breeding, care and feeding of animals are, at least in part, controlled by humans. Information regarding the changes in the protein composition of eggs in response to domestication is very limited. Such data are prerequisite for improvements in the reproduction of domesticated fish. The aim of this study was to examine the impact of domestication on the proteome of pikeperch eggs using two-dimensional differential in-gel electrophoresis. We analysed high-quality eggs from domesticated and wild pikeperch fish to reveal proteins that were presumably only related to the domestication process and not to the quality of eggs. Here, we show that domestication has a profound impact on the protein profile of pikeperch eggs. We identified 66 differentially abundant protein spots, including 27 spots that were more abundant in wild-caught pikeperch eggs and 39 spots that were enriched in eggs collected from domesticated females. Eggs originating from wild-caught females showed higher expression levels of proteins involved in folding, apoptotic process, purine metabolism and immune response, whereas eggs of domesticated females showed higher expression levels of proteins that participated mainly in metabolism. The changes in metabolic proteins in eggs from domesticated females can reflect the adaptation of pikeperch to commercial diets, which have profoundly distinct compositions compared with natural diets. The decrease in the abundance of proteins related to immune response in eggs from the domesticated population suggests that domestication may lead to disturbances in defence mechanisms. In turn, the lower abundance of heat shock proteins in eggs of domesticated fish may indicate their adaptation to stable farming conditions and reduced environmental stressors or their better tolerance of stress from breeding. The proteins identified in this study can increase our knowledge concerning the mechanism of the pikeperch domestication process.

Heat stress-induced alterations in the expression of genes associated with gonadal integrity of the teleost Puntius sophore

Temperature plays an important role on reproductive physiology of vertebrates including mammals, fish, and birds. It has varying effects on fish reproduction depending on the species; higher temperatures favor the spring-spawning species, while lower temperatures stimulate reproduction in autumn spawners. To evaluate the impact of high temperature on the reproductive physiology of minnow Puntius sophore, we carried out expression analysis of selected genes associated with gamete quality (hsp60, hsp70, hsp90, hsf1, vtg), pleuripotency (sox2, oct4, nanog), and sex determination (dmrt1) in gonads (ovary and testis) of P. sophore, heat stressed for different time periods (36 °C/7 days or 60 days) using real-time quantitative polymerase chain reaction (RT-qPCR). Expression of most of the hsp, vtg, and pleuripotency marker genes sox-2, oct-4, and nanog genes was downregulated in both ovary and testis of heat-stressed fish. The expression of dmrt-1 was upregulated in testis but downregulated in ovary of the heat-stressed fish which could be a male favoring effect of high temperature in P. sophore. This study suggests that the reproductive physiology and health of the nutrient dense P. sophore would be negatively affected by high temperature stress.

Transcriptomic Analysis Revealed the Regulatory Mechanisms of Oocyte Maturation and Hydration in Orange-Spotted Grouper (Epinephelus coioides)

Oocyte maturation and hydration are regulated by a complex interplay of various hormones and local factors. We have investigated the morphological changes of follicles and serum steroid levels during the HCG (human choionic gonadotophin)-induced oocyte maturation in the orange-spotted grouper. For the first time, a large-scale transcriptomic analysis of follicles during the maturation has been conducted in a fish species which produce pelagic oocytes. Eight cDNA libraries of follicle samples, from full-grown immature follicles to mature follicles, were constructed. A total of 402,530,284 high-quality clean reads were obtained after filtering, 79.66% of which perfectly mapped to the orange-spotted grouper genome. Real-time PCR results of 12 representative genes related to oocyte maturation and hydration verified the reliability of the RNA-seq data. A large number of genes related to oocyte maturation and hydration were identified in the transcriptome dataset. And the transcriptomic analysis revealed the dynamic changes of the steroid synthesis pathway and the pathway of hydration during oocyte maturation. The present study will facilitate future study on the oocyte maturation and hydration in the orange-spotted grouper and other marine pelagic egg spawner.

Transcriptome analysis of egg viability in rainbow trout, Oncorhynchus mykiss

Background

Maternal transcripts are accumulated in the oocyte during oogenesis to provide for protein synthesis from oocyte maturation through early embryonic development, when nuclear transcription is silenced. The maternal mRNAs have short poly(A) tails after undergoing post-transcriptional processing necessary for stabilizing them for storage. The transcripts undergo cytoplasmic polyadenylation when they are to be translated. Transcriptome analyses comparing total mRNA and elongated poly(A) mRNA content among eggs of different quality can provide insight into molecular mechanisms affecting egg developmental competence in rainbow trout. The present study used RNA-seq to compare transcriptomes of unfertilized eggs of rainbow trout females yielding different eyeing rates, following rRNA removal and poly(A) retention for construction of the libraries.

Results

The percentage of embryos to reach the 32-cell stage at 24 h post fertilization was significantly correlated to family eyeing rate, indicating that inviable embryos were developmentally compromised before zygotic genome activation. RNA sequencing identified 2 differentially expressed transcripts (DETs) from total mRNA sequencing comparing females with low-quality (< 5% eyeing), medium-quality (30–50% eyeing), and high-quality (> 80% eyeing) eggs. In contrast, RNA sequencing from poly(A) captured transcripts identified 945 DETs between low- and high-quality eggs, 1012 between low- and medium-quality eggs, and only 2 between medium- and high-quality eggs. The transcripts of mitochondrial genes were enriched with polyadenylated transcript sequencing and they were significantly reduced in low-quality eggs. Similarly, mitochondrial DNA was reduced in low-quality eggs compared with medium- and high-quality eggs. The functional gene analysis classified the 945 DETs between low- and high-quality eggs into 31 functional modules, many of which were related to ribosomal and mitochondrial functions. Other modules involved transcription, translation, cell division, apoptosis, and immune responses.

Conclusions

Our results indicate that differences in egg quality may be derived from differences in maternal nuclear transcript activation and cytoplasmic polyadenylation before ovulation, as opposed to accumulation and storage of maternal nuclear transcripts during oogenesis. Transcriptome comparisons suggest low-quality eggs suffered from impaired oxidative phosphorylation and translation. The DETs identified in this study provide insight into developmental competence in rainbow trout eggs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5690-5) contains supplementary material, which is available to authorized users.

Transcriptome Profile Analysis on Ovarian Tissues of Autotetraploid Fish and Diploid Red Crucian Carp

Polyploidization can significantly alter the size of animal gametes. Autotetraploid fish (RRRR, 4nRR = 200) (4nRR) possessing four sets of chromosomes were derived from whole-genome duplication in red crucian carp (RR, 2n = 100) (RCC). The diploid eggs of the 4nRR fish were significantly larger than the eggs of RCC. To explore the differences between the ovaries of these two ploidies of fishes at the molecular level, we compared the ovary transcriptome profiles of 4nRR fish and RCC and identified differentially expressed genes (DEGs). A total of 19,015 unigenes were differentially expressed between 4nRR fish and RCC, including 12,591 upregulated and 6,424 downregulated unigenes in 4nRR fish. Functional analyses revealed that eight genes (CDKL1, AHCY, ARHGEF3, TGFβ, WNT11, CYP27A, GDF7, and CKB) were involved in the regulation of cell proliferation, cell division, gene transcription, ovary development and energy metabolism, suggesting that these eight genes were related to egg size in 4nRR fish and RCC. We validated the expression levels of these eight DEGs in 4nRR fish and RCC using quantitative PCR. The study results provided insights into the regulatory mechanisms underlying the differences in crucian carp egg sizes.

Oogenesis and Egg Quality in Finfish: Yolk Formation and Other Factors Influencing Female Fertility

Egg quality in fishes has been a topic of research in aquaculture and fisheries for decades as it represents an important life history trait and is critical for captive propagation and successful recruitment. A major factor influencing egg quality is proper yolk formation, as most fishes are oviparous and the developing offspring are entirely dependent on stored egg yolk for nutritional sustenance. These maternally derived nutrients consist of proteins, carbohydrates, lipids, vitamins, minerals, and ions that are transported from the liver to the ovary by lipoprotein particles including vitellogenins. The yolk composition may be influenced by broodstock diet, husbandry, and other intrinsic and extrinsic conditions. In addition, a number of other maternal factors that may influence egg quality also are stored in eggs, such as gene transcripts, that direct early embryonic development. Dysfunctional regulation of gene or protein expression may lead to poor quality eggs and failure to thrive within hours of fertilization. These gene transcripts may provide important markers as their expression levels may be used to screen broodstock for potential spawning success. In addition to such intrinsic factors, stress may lead to ovarian atresia or reproductive failure and can impact fish behavior, fecundity, and ovulation rate. Finally, postovulatory aging may occur when eggs become overripe and the fish fails to spawn in a timely fashion, leading to low fertility, often encountered during manual strip spawning of fish.

Full genome survey and dynamics of gene expression in the greater amberjack Seriola dumerili

Background

Teleosts of the genus Seriola, commonly known as amberjacks, are of high commercial value in international markets due to their flesh quality and worldwide distribution. The Seriola species of interest to Mediterranean aquaculture is the greater amberjack (Seriola dumerili). This species holds great potential for the aquaculture industry, but in captivity, reproduction has proved to be challenging, and observed growth dysfunction hinders their domestication. Insights into molecular mechanisms may contribute to a better understanding of traits like growth and sex, but investigations to unravel the molecular background of amberjacks have begun only recently.

Findings

Illumina HiSeq sequencing generated a high-coverage greater amberjack genome sequence comprising 45 909 scaffolds. Comparative mapping to the Japanese yellowtail (Seriola quinqueriadiata) and to the model species medaka (Oryzias latipes) allowed the generation of in silico groups. Additional gonad transcriptome sequencing identified sex-biased transcripts, including known sex-determining and differentiation genes. Investigation of the muscle transcriptome of slow-growing individuals showed that transcripts involved in oxygen and gas transport were differentially expressed compared with fast/normal-growing individuals. On the other hand, transcripts involved in muscle functions were found to be enriched in fast/normal-growing individuals.

Conclusion

The present study provides the first insights into the molecular background of male and female amberjacks and of fast- and slow-growing fish. Therefore, valuable molecular resources have been generated in the form of a first draft genome and a reference transcriptome. Sex-biased genes, which may also have roles in sex determination or differentiation, and genes that may be responsible for slow growth are suggested.

Comparison of proteomic profiles in the ovary of Sterlet sturgeon (Acipenser ruthenus) during vitellogenic stages

One of the challenges of sturgeon aquaculture is that sturgeon takes an extended amount of time to reach sexual maturity. The pattern of the protein expression in relation to the late maturity of sturgeon can help to better understand changes in sexual maturity. 17β-estradiol (E2), testosterone (T) and vitellogenin (Vtg) levels were examined at all stages of sexual maturation in Sterlet sturgeon (Acipenser ruthenus). Two-dimensional gel electrophoresis and mass spectrometry analysis were used to show the pattern of the ovarian proteins. The T levels increased from the previtellogenic to the postvitellogenic stages (P < 0.05) and Vtg showed a decremental pattern in pre- and postvitellogenic, and atresia (not significantly). The analysis showed 900 protein spots, 19 of which were successfully identified and had significant differences between the previtellogenic and the vitellogenic groups (P < 0.05). Among the identified proteins, 40% involved in cell defense (heat shock protein, Glutathione peroxidase, natural killer enhancing factor, peroxiredoxin-2), 30% in transcription and translation (constitutive photomorphogenesis 9 and Ybx2), 20% in metabolism and energy production (triose-phosphate isomerase (TPI)) and 10% in transport (glycolipid transfer protein). In the vitellogenic stage, the proteins were related to metabolism and energy production (TPI, ES1, creatin kinase, enolase, nucleoside diphosphate kinase, 50%), cell defense (thioredoxin and dislophid isomerase, 20%) and transport (fatty acid binding protein, 10%). Our findings show changes in protein expression pattern from cell defense to metabolism during egg development.

DNA methylation profiles correlated to striped bass sperm fertility

Background

Striped bass (Morone saxatilis) spermatozoa are used to fertilize in vitro the eggs of white bass (M. chrysops) to produce the preferred hybrid for the striped bass aquaculture industry. Currently, only one source of domestic striped bass juveniles is available to growers that is not obtained from wild-caught parents and is thus devoid of any genetic improvement in phenotypic traits of importance to aquaculture. Sperm epigenetic modification has been predicted to be associated with fertility, which could switch genes on and off without changing the DNA sequence itself. DNA methylation is one of the most common epigenetic modification types and changes in sperm epigenetics can be correlated to sub-fertility or infertility in male striped bass. The objective of this study was to find the differentially methylated regions (DMRs) between high-fertility and sub-fertility male striped bass, which could potentially regulate the fertility performance.

Results

In our present study, we performed DNA methylation analysis of high-fertility and sub-fertility striped bass spermatozoa through MBD-Seq methods. A total of 171 DMRs were discovered in striped bass sperm correlated to fertility. Based on the annotation of these DMRs, we conducted a functional classification analysis and two important groups of genes including the WDR3/UTP12 and GPCR families, were discovered to be related to fertility performance of striped bass. Proteins from the WDR3/UTP12 family are involved in forming the sperm flagella apparatus in vertebrates and GPCRs are involved in hormonal signaling and regulation of tissue development, proliferation and differentiation.

Conclusions

Our results contribute insights into understanding the mechanism of fertility in striped bass, which will provide powerful tools to maximize reproductive efficiencies and to identify those males with superior gametes for this important aquaculture species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4548-6) contains supplementary material, which is available to authorized users.

RNA sequencing of early round goby embryos reveals that maternal experiences can shape the maternal RNA contribution in a wild vertebrate

Background

It has been proposed that non-genetic inheritance could promote species fitness. Non-genetic inheritance could allow offspring to benefit from the experience of their parents, and could advocate pre-adaptation to prevailing and potentially selective conditions. Indeed, adaptive parental effects have been modeled and observed, but the molecular mechanisms behind them are far from understood.

Results

In the present study, we investigated whether maternal RNA can carry information about environmental conditions experienced by the mother in a wild vertebrate. Maternal RNA directs the development of the early embryo in many non-mammalian vertebrates and invertebrates. However, it is not known whether vertebrate maternal RNA integrates information about the parental environment. We sequenced the maternal RNA contribution from a model that we expected to rely on parental effects: the invasive benthic fish species Neogobius melanostomus (Round Goby). We found that maternal RNA expression levels correlated with the water temperature experienced by the mother before oviposition, and identified temperature-responsive gene groups such as core nucleosome components or the microtubule cytoskeleton.

Conclusions

Our findings suggest that the maternal RNA contribution may incorporate environmental information. Maternal RNA should therefore be considered a potentially relevant pathway for non-genetic inheritance. Also, the ability of a species to integrate environmental information in the maternal RNA contribution could potentially contribute to species fitness and may also play a role in extraordinary adaptive success stories of invasive species such as the round goby.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1132-2) contains supplementary material, which is available to authorized users.

Leptin Stimulates Cellular Glycolysis Through a STAT3 Dependent Mechanism in Tilapia

We assessed if leptin, a cytokine hormone known to enhance energy expenditure by promoting lipid and carbohydrate catabolism in response to physiologic stress, might directly regulate cellular glycolysis. A transcriptomic analysis of prolactin cells in the tilapia (Oreochromis mossambicus) pituitary rostral pars distalis (RPD) revealed that recombinant leptin (rtLep) differentially regulates 1,995 genes, in vitro. Machine learning algorithms and clustering analyses show leptin influences numerous cellular gene networks including metabolism; protein processing, transport, and metabolism; cell cycle and the hypoxia response. Leptin stimulates transcript abundance of the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (gapdh) in a covariate manner to the hypoxic stress gene network. Orthogonal tests confirm that rtLepA dose-dependently increases gapdh gene expression in the RPD along with transcript abundance of 6-phosphofructo-1-kinase (pfk1), the rate limiting glycolytic enzyme. Functional testing demonstrated that leptin stimulates PFK activity and glycolytic output, while Stattic (a STAT3 blocker) was sufficient to suppress these responses, indicating leptin stimulates glycolysis through a STAT3-dependent mechanism. Leptin also stimulated pfk1 gene expression and lactate production in primary hepatocyte incubations in a similar manner to those shown for the pituitary RPD. This work characterizes a critical metabolic action of leptin to directly stimulate glycolysis across tissue types in a teleost model system, and suggest that leptin may promote energy expenditure, in part, by stimulating glycolysis. These data in a teleost fish, suggest that one of leptin's ancient, highly-conserved functions among vertebrates may be stimulation of glycolysis to facilitate the energetic needs associated with various stressors.

Scrambled eggs: Proteomic portraits and novel biomarkers of egg quality in zebrafish (Danio rerio)

Egg quality is a complex biological trait and a major determinant of reproductive fitness in all animals. This study delivered the first proteomic portraits of egg quality in zebrafish, a leading biomedical model for early development. Egg batches of good and poor quality, evidenced by embryo survival for 24 h, were sampled immediately after spawning and used to create pooled or replicated sample sets whose protein extracts were subjected to different levels of fractionation before liquid chromatography and tandem mass spectrometry. Obtained spectra were searched against a zebrafish proteome database and detected proteins were annotated, categorized and quantified based on normalized spectral counts. Manually curated and automated enrichment analyses revealed poor quality eggs to be deficient of proteins involved in protein synthesis and energy and lipid metabolism, and of some vitellogenin products and lectins, and to have a surfeit of proteins involved in endo-lysosomal activities, autophagy, and apoptosis, and of some oncogene products, lectins and egg envelope proteins. Results of pathway and network analyses suggest that this aberrant proteomic profile results from failure of oocytes giving rise to poor quality eggs to properly transit through final maturation, and implicated Wnt signaling in the etiology of this defect. Quantitative comparisons of abundant proteins in good versus poor quality eggs revealed 17 candidate egg quality markers. Thus, the zebrafish egg proteome is clearly linked to embryo developmental potential, a phenomenon that begs further investigation to elucidate the root causes of poor egg quality, presently a serious and intractable problem in livestock and human reproductive medicine.

Transcriptomic signatures for ovulation in vertebrates

The central roles of luteinizing hormone (LH), progestin and their receptors for initiating ovulation have been well established. However, signaling pathways and downstream targets such as proteases that are essential for the rupture of follicular cells are still unclear. Recently, we found anovulation in nuclear progestin receptor (Pgr) knockout (Pgr-KO) zebrafish, which offers a new model for examining genes and pathways that are important for ovulation and fertility. In this study, we examined expression of all transcripts using RNA-Seq in preovulatory follicular cells collected following the final oocyte maturation, but prior to ovulation, from wild-type (WT) or Pgr-KO fish. Differential expression analysis revealed 3567 genes significantly differentially expressed between WT and Pgr-KO fish (fold change⩾2, p<0.05). Among those, 1543 gene transcripts were significantly more expressed, while 2024 genes were significantly less expressed, in WT than those in Pgr-KO. We then retrieved and compared transcriptional data from online databases and further identified 661 conserved genes in fish, mice, and humans that showed similar levels of high (283 genes) or low (387) expression in animals that were ovulating compared to those with no ovulation. For the first time, ovulatory genes and their involved biological processes and pathways were also visualized using Enrichment Map and Cytoscape. Intriguingly, enrichment analysis indicated that the genes with higher expression were involved in multiple ovulatory pathways and processes such as inflammatory response, angiogenesis, cytokine production, cell migration, chemotaxis, MAPK, focal adhesion, and cytoskeleton reorganization. In contrast, the genes with lower expression were mainly involved in DNA replication, DNA repair, DNA methylation, RNA processing, telomere maintenance, spindle assembling, nuclear acid transport, catabolic processes, and nuclear and cell division. Our results indicate that a large set of genes (>3000) is differentially regulated in the follicular cells in zebrafish prior to ovulation, terminating programs such as growth and proliferation, and beginning processes including the inflammatory response and apoptosis. Further studies are required to establish relationships among these genes and an ovulatory circuit in the zebrafish model.

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.

Tissue localization of piscidin host-defense peptides during striped bass (Morone saxatilis) development

Infectious diseases are a major cause of larval mortality in finfish aquaculture. Understanding ontogeny of the fish immune system and thus developmental timing of protective immune tissues and cells, may help to decrease serious losses of larval fishes when they are particularly vulnerable to infection. One component of the innate immune system of fishes is the host-defense peptides, which include the piscidins. Piscidins are small, amphipathic, α-helical peptides with a broad-spectrum of action against viral, bacterial, fungal, and protozoan pathogens. We describe for the first time the cellular and tissue localization of three different piscidins (1, 3, and 4) during striped bass (Morone saxatilis) larval ontogeny using immunofluorescent histochemistry. From 16 days post hatch to 12 months of age, piscidin staining was observed in cells of the epithelial tissues of gill, digestive tract, and skin, mainly in mast cells. Staining was also seen in presumptive hematopoietic cells in the head kidney. The three piscidins showed variable cellular and tissue staining patterns, possibly relating to differences in tissue susceptibility or pathogen specificity. This furthers our observation that the piscidins are not a monolithic family of antimicrobials, but that different AMPs have different (more specialized) functions. Furthermore, no immunofluorescent staining of piscidins was observed in post-vitellogenic oocytes, embryos, or larvae from hatch to 14 days post hatch, indicating that this critical component of the innate immune system is inactive in pre-hatch and young larval striped bass.

Transcriptomic Profiling of Egg Quality in Sea Bass (Dicentrarchus labrax) Sheds Light on Genes Involved in Ubiquitination and Translation

Variable and low egg quality is a major limiting factor for the development of efficient aquaculture production. This stems from limited knowledge on the mechanisms underlying egg quality in cultured fish. Molecular analyses, such as transcriptomic studies, are valuable tools to identify the most important processes modulating egg quality. However, very few studies have been devoted to this aspect so far. Within this study, the microarray-based transcriptomic analysis of eggs (of different quality) of sea bass (Dicentrarchus labrax) was performed. An Agilent oligo microarray experiment was performed on labelled mRNA extracted from 16 batches of eggs (each batch obtained from a different female) of sea bass, in which over 24,000 published probe arrays were used. We identified 39 differentially expressed genes exhibiting a differential expression between the groups of low (fertilization rate < 60 %) and high (fertilization rate > 60 %) quality. The mRNA levels of eight genes were further analyzed by quantitative PCR. Seven genes were confirmed by qPCR to be differentially expressed in eggs of low and high quality. This study confirmed the importance of some of the genes already reported to be potential molecular quality indicators (mainly rnf213 and irf7), but we also found new genes (mainly usp5, mem-prot, plec, cenpf), which had not yet been reported to be quality-dependent in fish. These results suggest the importance of genes involved in several important processes, such as protein ubiquitination, translation, DNA repair, and cell structure and architecture; these probably being the mechanisms that contribute to egg developmental competence in sea bass.

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.

Proteomic analysis of early-stage embryos: implications for egg quality in hapuku (Polyprion oxygeneios)

In order to develop biomarkers that may help predict the egg quality of captive hapuku (Polyprion oxygeneios) and provide potential avenues for its manipulation, the present study (1) sequenced the proteome of early-stage embryos using isobaric tag for relative and absolute quantification analysis, and (2) aimed to establish the predictive value of the abundance of identified proteins with regard to egg quality through regression analysis. Egg quality was determined for eight different egg batches by blastomere symmetry scores. In total, 121 proteins were identified and assigned to one of nine major groups according to their function/pathway. A mixed-effects model analysis revealed a decrease in relative protein abundance that correlated with (decreasing) egg quality in one major group (heat-shock proteins). No differences were found in the other protein groups. Linear regression analysis, performed for each identified protein separately, revealed seven proteins that showed a significant decrease in relative abundance with reduced blastomere symmetry: two correlates that have been named in other studies (vitellogenin, heat-shock protein-70) and a further five new candidate proteins (78 kDa glucose-regulated protein, elongation factor-2, GTP-binding nuclear protein Ran, iduronate 2-sulfatase and 6-phosphogluconate dehydrogenase). Notwithstanding issues associated with multiple statistical testing, we conclude that these proteins, and especially iduronate 2-sulfatase and the generic heat-shock protein group, could serve as biomarkers of egg quality in hapuku.

Transcriptomics of mRNA and egg quality in farmed fish: Some recent developments and future directions

Maternal mRNA transcripts deposited in growing oocytes regulate early development and are under intensive investigation as determinants of egg quality. The research has evolved from single gene studies to microarray and now RNA-Seq analyses in which mRNA expression by virtually every gene can be assessed and related to gamete quality. Such studies have mainly focused on genes changing two- to several-fold in expression between biological states, and have identified scores of candidate genes and a few gene networks whose functioning is related to successful development. However, ever-increasing yields of information from high throughput methods for detecting transcript abundance have far outpaced progress in methods for analyzing the massive quantities of gene expression data, and especially for meaningful relation of whole transcriptome profiles to gamete quality. We have developed a new approach to this problem employing artificial neural networks and supervised machine learning with other novel bioinformatics procedures to discover a previously unknown level of ovarian transcriptome function at which minute changes in expression of a few hundred genes is highly predictive of egg quality. In this paper, we briefly review the progress in transcriptomics of fish egg quality and discuss some future directions for this field of study.

Machine learning reveals sex-specific 17β-estradiol-responsive expression patterns in white perch (Morone americana) plasma proteins

With growing abundance and awareness of endocrine disrupting compounds (EDCs) in the environment, there is a need for accurate and reliable detection of EDC exposure. Our objective in the present study was to observe differences within and between the global plasma proteomes of sexually mature male and female white perch (Morone americana) before (Initial Control, IC) and after 17β-estradiol (E2 ) induction. Semiquantitative nanoLC-MS/MS data were analyzed by machine learning support vector machines (SVMs) and by two-way ANOVA. By ANOVA, the expression levels of 44, 77, and 57 proteins varied significantly by gender, treatment, and the interaction of gender and treatment, respectively. SVMs perfectly classified male and female perch IC and E2 -induced plasma samples using the protein expression data. E2 -induced male and female perch plasma proteomes contained significantly higher levels of the yolk precursors vitellogenin Aa and Ab (VtgAa, VtgAb), as well as latrophilin and seven transmembrane domain-containing protein 1 (Eltd1) and kininogen 1 (Kng1). This is the first report that Eltd1 and Kng1 may be E2 -responsive proteins in fishes and therefore may be useful indicators of estrogen induction.

Transcriptome Profile Analysis of Ovarian Tissues from Diploid and Tetraploid Loaches Misgurnus anguillicaudatus

RNA sequencing and short-read assembly was utilized to produce a transcriptome of ovarian tissues from three-year-old diploid and tetraploid loaches (Misgurnus anguillicaudatus). A total of 28,369 unigenes were obtained, comprising 10,546 unigenes with length longer than 1000 bp. More than 73% of the unigenes were annotated through sequence comparison with databases. The RNA-seq data revealed that 2253 genes were differentially expressed between diploid and tetraploid loaches, including 1263 up-regulated and 990 down-regulated genes in tetraploid loach. Some differentially expressed genes, such as vitellogenin (Vtg), gonadotropin releasing hormone receptor type A (GnRHRA), steroidogenic acute regulatory protein (StAR), mitogen-activated protein kinase 14a (MAPK14a), ATP synthase subunit alpha (atp5a), and synaptonemal complex protein 1 (Scp1), were involved in regulation of cell proliferation, division, gene transcription, ovarian development and energy metabolism, suggesting that these genes were related to egg diameter of the loach. Results of transcriptome profiling here were validated using real time quantitative PCR in ten selected genes. The present study provided insights into the transcriptome profile of ovarian tissues from diploid and tetraploid loaches Misgurnus anguillicaudatus, which was made available to the research community for functional genomics, comparative genomics, polyploidy evolution and molecular breeding of this loach and other related species.

Considerations for the use of transcriptomics in identifying the ‘genes that matter’ for environmental adaptation

Transcriptomics has emerged as a powerful approach for exploring physiological responses to the environment. However, like any other experimental approach, transcriptomics has its limitations. Transcriptomics has been criticized as an inappropriate method to identify genes with large impacts on adaptive responses to the environment because: (1) genes with large impacts on fitness are rare; (2) a large change in gene expression does not necessarily equate to a large effect on fitness; and (3) protein activity is most relevant to fitness, and mRNA abundance is an unreliable indicator of protein activity. In this review, these criticisms are re-evaluated in the context of recent systems-level experiments that provide new insight into the relationship between gene expression and fitness during environmental stress. In general, these criticisms remain valid today, and indicate that exclusively using transcriptomics to screen for genes that underlie environmental adaptation will overlook constitutively expressed regulatory genes that play major roles in setting tolerance limits. Standard practices in transcriptomic data analysis pipelines may also be limiting insight by prioritizing highly differentially expressed and conserved genes over those genes that undergo moderate fold-changes and cannot be annotated. While these data certainly do not undermine the continued and widespread use of transcriptomics within environmental physiology, they do highlight the types of research questions for which transcriptomics is best suited and the need for more gene functional analyses. Such information is pertinent at a time when transcriptomics has become increasingly tractable and many researchers may be contemplating integrating transcriptomics into their research programs.