Expression of vitellogenin receptor gene in the ovary of wild and captive Atlantic bluefin tuna (Thunnus thynnus)

Comparison of ovarian mRNA expression levels in wild and hatchery-produced greater amberjack Seriola dumerili

The greater amberjack Seriola dumerili is a promising candidate for aquaculture production. This study compares the ovary transcriptome of greater amberjack sampled in the wild (WILD) with hatchery-produced breeders reared in aquaculture sea cages in the Mediterranean Sea. Among the seven sampled cultured fish, three were classified as reproductively dysfunctional (DysF group), while four showed no signs of reproductive alteration (NormalF group). The DysF fish showed 1,166 differentially expressed genes (DEGs) compared to WILD females, and 755 DEGs compared to the NormalF. According to gene ontology (GO) analysis, DysF females exhibited enrichment of genes belonging to the biological categories classified as Secreted, ECM-receptor interaction, and Focal adhesion. Protein-protein interaction analysis revealed proteins involved in the biological categories of ECM-receptor interaction, Enzyme-linked receptor protein signaling, Wnt signal transduction pathways, and Ovulation cycle. KEGG pathway analysis showed DEGs involved in 111 pathways, including Neuroactive ligand-receptor interaction, Steroid hormone biosynthesis, Cell cycle, Oocyte meiosis, Necroptosis, Ferroptosis, Apoptosis, Autophagy, Progesterone-mediated oocyte maturation, Endocytosis and Phagosome, as well as Hedgehog, Apelin, PPAR, Notch, and GnRH signalling pathways. Additionally, DysF females exhibited factors encoded by upregulated genes associated with hypogonadism and polycystic ovary syndrome in mammals. This study -which is part of a broader research effort examining the transcriptome of the entire reproductive axis in greater amberjack of both sexes-, enhances our comprehension of the mechanisms underlying the appearance of reproductive dysfunctions when fish are reared under aquaculture conditions.

The role of multiple vitellogenins in early development of fishes

Functions of vitellogenins have been in the limelight of fish reproductive physiology research for decades. The Vtg system of acanthomorph teleosts consists of two complete forms of Vtgs (VtgAa and VtgAb) and an incomplete form, VtgC. Insufficient uptake and processing of Vtgs and their yolk proteins lead to inadequate oocyte hydration ensuing failure in acquisition of egg buoyancy and early developmental deficiencies. This review presents a summary of our studies on utilization of multiple Vtgs in species with different egg buoyancy characteristics, as examples. Studies of moronids revealed limited degradation of all three forms of lipovitellin heavy chain derived from their three respective forms of Vtg, by which they contribute to the free amino acid pool driving oocyte hydration during oocyte maturation. In later studies, CRISPR/Cas9 was employed to invalidate zebrafish type I, type II and type III Vtgs, which are orthologs of acanthamorph VtgAa, VtgAb and VtgC, respectively. Results revealed type I Vtg to have essential developmental and nutritional functions in both late embryos and larvae. Genomic disturbance of type II Vtg led to high mortalities during the first 24 h of embryonic development. Despite being a minor form of Vtg in zebrafish and most other species, type III Vtg was also found to contribute essentially to the developmental potential of zebrafish zygotes and early embryos. Apart from severe effects on progeny survival, these studies also disclosed previously unreported regulatory effects of Vtgs on fecundity and fertility, and on embryo hatching. We recently utilized parallel reactions monitoring based liquid chromatography tandem mass spectrometry to assess the processing and utilization of lipovitellins derived from different forms of Vtg in Atlantic halibut and European plaice. Results showed the Lv heavy chain of VtgAa (LvHAa) to be consumed during oocyte maturation and the Lv light chain of VtgAb (LvLAb) to be utilized specifically during late larval stages, while all remaining YPs (LvLAa, LvHAb, LvHC, and LvLC) were utilized during or after hatching up until first feeding in halibut. In plaice, all YPs except LvHAa, which similarly to halibut supports oocyte maturation, are utilized from late embryo to late larval development up until first feeding. The collective findings from these studies affirm substantial disparity in modes of utilization of different types of Vtgs among fish species with various egg buoyancy characteristics, and they reveal previously unknown regulatory functions of Vtgs in maintenance of reproductive assets such as maternal fecundity and fertility, and in embryonic hatching. Despite the progress that has been made over the past two decades by examining multiple Vtgs and their functions, a higher complexity of these systems with much greater diversity between species in modes of Vtg utilization is now evident. Further research is needed to reveal novel ways each species has evolved to utilize these complex multiple Vtg systems and to discover unifying principles for this evolution in fishes of diverse lineages, habitats and life history characteristics.

Spatio-temporal patterns of ovarian development and VgR gene silencing reduced fecundity in parthenogenetic Artemia

The halophilic zooplankton brine shrimp Artemia has been used as an experimental animal in multidisciplinary studies. However, the reproductive patterns and its regulatory mechanisms in Artemia remain unclear. In this study, the ovarian development process of parthenogenetic Artemia (A. parthenogenetica) was divided into five stages, and oogenesis or egg formation was identified in six phases. The oogenesis mode was assumed to be polytrophic. We also traced the dynamic translocation of candidate germline stem cells (cGSCs) using EdU labelling and elucidated several key cytological events in oogenesis through haematoxylin and eosin staining and fluorescence imaging. Distinguished from the ovary structure of insects and crustaceans, Artemia germarium originated from ovariole buds and are located at the base of the ovarioles. RNA-seq based on five stages of ovarian development identified 2657 upregulated genes related to reproduction by pair-to-pair comparison. Gbb, Dpp, piwi, vasa, nanos, VgA and VgR genes associated with cGSCs recognition and reproductive development were screened and verified using qPCR. Silencing of the VgR gene in A. parthenogenetica (Ap-VgR) at ovarian development Stage II led to a low level of gene expression (less than 10%) within 5 days, which resulted in variations in oogenesis-related gene expression and significantly inhibited vitellogenesis, impeded oocyte maturation, and eventually decreased the number of offspring. In conclusion, we have illustrated the patterns of ovarian development, outlined the key spatio-temporal features of oogenesis and identified the negative impacts of VgR gene knockdown on oogenesis using A. parthenogenetica as an experimental animal. The findings of this study also lay a foundation for the further study of reproductive biology of invertebrates.

Vitellogenin Ab structure of the amazonian Arapaima gigas

The vitellogenin is composed by polypeptides that are precursors of egg yolk proteins that provides embryo and larvae nutrition. The mRNA encoding for vitellogenin Ab (Vtg-Ab; 4,536 bp long and 1,512 amino acids) were obtained by RNA-Seq library sequencing of pirarucu gonads. The Vtg-Ab sequences had high homology with Vtgs of other three teleosts species of the order Osteoglossiformes. The transcript of ovarian Vtg was identified based on structural criteria, and so we classify the Vtg of pirarucu as Vtg-Ab due to the truncated or shortened phosvitin (N-terminal end) and phosvitinless domain (C-terminal end). The Vtg-Ab of pirarucu present two major deletions with 133 amino acids in the Lipovitellin I domain and 89 amino acids in the truncated or shortened Phosvitin domain, both located in the N-terminal end region. The three-dimensional (3-D) structure Vtg-Ab protein shows the presence of a typical 4α-helices bundle protein that runs in anti-parallel. In general, the characterization of Vtg-Ab may be the useful elucidation of the hormonal regulation of vitellogenesis and improve the production of pirarucu for broodstock management in aquaculture and preparation of Vtg antibody production (species-specific) for sex identification.

A Comparison of Reproductive Performances in Young and Old Females: A Case Study on the Atlantic Bluefin Tuna in the Mediterranean Sea

Simple Summary

The Atlantic bluefin tuna Thunnus thynnus is a species characterized by complex trans-oceanic migrations linked to size, which rely on the delicate trade-off between somatic growth and reproduction before and during the migratory movements to reach spawning grounds. Therefore, understanding the processes that drive reproduction and elucidating its age-related regulation is essential in the context of sustainable fishery management. In this study, carried out in the Mediterranean Sea, older bluefin tuna females were found to have greater reproductive performances than younger females according to a molecular biology approach (i.e., gene expression), a result that likely mirrors a better physical condition, different habitat usage or migratory behaviour. This result highlights the importance of preserving large females for their major reproductive contribution at a stock level. Furthermore, the gonad-specific mir-202, which belongs to a class of non-coding RNA, called miRNA, that regulate the post-transcription of protein-coding genes, was identified as a potential candidate to play a role in egg quality and quantity (i.e., fecundity) during ovarian maturation through age- or stage-dependent reproductive processes. Overall, the present study contributes to improve the sustainability of the Atlantic bluefin tuna fishery in the Mediterranean Sea.

Abstract

In the Mediterranean Sea, a demographic substructure of the Atlantic bluefin tuna Thunnus thynnus has emerged over the last decade, with old and young individuals exhibiting different horizontal movements and spatial–temporal patterns of gonad maturation. In the present study, histology and molecular reproductive markers were integrated with the gonad-specific mir-202 gene expression and ovarian localization to provide a comprehensive picture of the reproductive performances in young and old females and investigate the role played by the mir-202 during gonadal maturation. During the reproductive period, old females (>100 kg; 194.6 ± 33.9 cm straight fork length; 11.3 ± 2.7 years old) were found to have greater reproductive performances than younger females (<80 kg; 139.3 ± 18.8 cm straight fork length; 8.4 ± 1.1 years old) according to gene expression results, suggesting a prolonged spawning season, earlier arrival on spawning grounds and/or better condition in older females. The mir-202-5p showed no global changes; it was abundantly expressed in granulosa cells and faintly present in the ooplasm. On the other hand, the mir-202-3p expression profile reflected levels of oocyte maturation molecular markers (star, lhr) and both histological and molecular (casp3) levels of follicular atresia. Overall, old females exhibited greater reproductive performances than younger females, likely reflecting different reproductive dynamics linked to the physical condition, habitat usage and migratory behaviour. These results highlight the importance of preserving large and old females in the context of fishery management. Finally, the mir-202 appears to be a good candidate to regulate the reproductive output of this species in an autocrine/paracrine manner through either stage- or age-dependent processes.

Atresia of ovarian follicles in fishes, and implications and uses in aquaculture and fisheries

Atresia of ovarian follicles, that is the degenerative process of germ cells and their associated somatic cells, is a complex process involving apoptosis, autophagy and heterophagy. Follicular atresia is a normal component of fish oogenesis and it is observed throughout the ovarian cycle, although it is more frequent in regressing ovaries during the postspawning period. An increased occurrence of follicular atresia above physiological rates reduces fish fecundity and even causes reproductive failure in both wild and captive-reared fish stocks, and hence, this phenomenon has a wide range of implications in applied sciences such as fisheries and aquaculture. The present article reviews the available literature on both basic and applied traits of oocyte loss by atresia, including its morpho-physiological aspects and factors that cause a supraphysiological increase of follicular atresia. Finally, the review presents the use of early follicular atresia identification in the selection process of induced spawning in aquaculture and the implications of follicular atresia in fisheries management.

Identification and structural characterization of the factors involved in vitellogenesis and its regulation in the African Osteoglossiforme of aquacultural interest Heterotis niloticus (Cuvier, 1829)

The African bonytongue (Heterotis niloticus) is an excellent candidate for fish farming because it has outstanding biological characteristics and zootechnical performances. However, the absence of sexual dimorphism does not favor its reproduction in captivity or the understanding of its reproductive behavior. Moreover, no molecular data related to its reproduction is yet available. This study therefore focuses on the structural identification of the different molecular actors of vitellogenesis expressed in the pituitary gland, the liver and the ovary of H. niloticus. A transcriptomic approach based on de novo RNA sequencing of the pituitary gland, ovary and liver of females in vitellogenesis led to the creation of three transcriptomes. In silico analysis of these transcriptomes identified the sequences of pituitary hormones such as prolactin (PRL), luteinizing hormone (LH) and follicle-stimulating hormone (FSH) and their ovarian receptors (PRLR, FSHR, LHR). In the liver and ovary, estrogen receptors (ER) beta and gamma, liver vitellogenins (VtgB and VtgC) and their ovarian receptors (VLDLR) were identified. Finally, the partial transcript of an ovarian Vtg weakly expressed compared to hepatic Vtg was identified based on structural criteria. Moreover, a proteomic approach carried out from mucus revealed the presence of one Vtg exclusively in females in vitellogenesis. In this teleost fish that does not exhibit sexual dimorphism, mucus Vtg could be used as a sexing biomarker based on a non-invasive technique compatible with the implementation of experimental protocols in vivo.

Increased melanomacrophage centres in the liver of reproductively dysfunctional female greater amberjack Seriola dumerili (Risso, 1810)

The greater amberjack Seriola dumerili is a new aquaculture fish that may display reproductive dysfunctions. During extensive follicular atresia, which is a common reproductive dysfunction in females during vitellogenesis, part of the reabsorbed yolk returns to the liver to be metabolized and recycled. Melanomacrophage centres (MMCs) are aggregates of macrophage-like cells that play a role in the destruction, detoxification and recycling of endogenous and exogenous materials, and have been associated with systemic stress. Wild and captive-reared greater amberjack were sampled in the Mediterranean Sea during two different phases of the reproductive cycle. The liver of reproductively dysfunctional captive-reared females sampled during the spawning season showed a high density of both MMCs and apoptotic cells. A weak liver anti-cytochrome P450 monooxygenase 1A immunoreactivity was observed, suggesting that the examined fish were not exposed to environmental pollutants. We propose that the observed increase in MMCs and apoptosis in captive-reared fish was related to the hepatic overload associated to the metabolism of yolk proteins reabsorbed during extensive follicular atresia. Since follicular atresia is a frequent physiological and pathological event in teleosts, we suggest that the reproductive state should be always assessed when MMCs are used as markers of exposure to stress or pollutants.

Vitellogenin receptor and fatty acid profiles of individual lipid classes of oocytes from wild and captive-reared greater amberjack (Seriola dumerili) during the reproductive cycle

The greater amberjack Seriola dumerili (Risso, 1810) is a large migratory pelagic fish occurring in tropical and temperate waters with a great potential for the world aquaculture industry. Previous studies showed that wild-caught female greater amberjack reared in sea cages and handled during the reproductive season, underwent extensive ovarian atresia. This atresia, however, was not related to an insufficient liver transcription or oocyte uptake of vitellogenin (Vtg). In the present study, the structure of two greater amberjack vitellogenin receptors, namely Vtgr (Lr8-) and Lrp13, was characterized. Moreover, vtgr and lrp13 gene expression and the fatty acid profiles of specific phospholipids and neutral lipids were compared in the ovaries of wild and captive-reared greater amberjack during different phases of the reproductive cycle (i.e. early gametogenesis, advanced gametogenesis and spawning). Ovarian vtgr and lrp13 transcription was more active during early gametogenesis, suggesting that vitellogenin receptor transcripts were synthesized by previtellogenic oocytes and remained in the cellular mRNA pool until oocytes resumed meiosis and entered into secondary growth (i.e. vitellogenesis). Rearing of wild-caught greater amberjack in captivity together with handling during the reproductive season was associated with a reduced vtgr and lrp13 transcription and with a diminished capacity of oocytes in the early phase of gametogenesis (primary oocyte growth) to enter into vitellogenesis. During early gametogenesis, remarkable differences in the fatty acid composition were observed between wild and captive-reared individuals: all phospholipids of captive fish displayed dramatic increases of saturates (16:0 and 18:0) and decreases of arachidonic acid (ARA) and docosahexaenoic acid (DHA). The present study confirms the susceptibility of greater amberjack reproductive function to handling stress and suggests that the consequent extensive atresia of vitellogenic follicles originated during the primary oocytes growth when the capacity of oocytes to synthesize vitellogenin receptors was reduced. The study also suggests that this reduced capacity was associated with an altered oocyte phospholipid fatty acid composition during early gametogenesis.

A de novo transcriptome assembly approach elucidates the dynamics of ovarian maturation in the swordfish (Xiphias gladius)

The Mediterranean swordfish (Xiphias gladius) has been recently classified as overfished and in 2016, the International Commission for the Conservation of the Atlantic Tunas (ICCAT) established a multi-annual management plan to recover this stock. To successfully achieve this goal, knowledge about swordfish biology is needed. To date, few studies on swordfish have been performed and none of them has provided useful insights into the reproductive biology at molecular level. Here we set to characterise the molecular dynamics underlying ovarian maturation by employing a de novo transcriptome assembly approach. Differential gene expression analysis in mature and immature ovaries identified a number of differentially expressed genes associated with biological processes driving ovarian maturation. Focusing on ovarian steroidogenesis and vitellogenin uptake, we depict the molecular dynamics characterizing these processes while a phylogenetic analysis let us identify a candidate vitellogenin receptor. This is the first swordfish transcriptome assembly and these findings provide in-depth understanding of molecular processes describing ovarian maturation. Moreover, the establishment of a publicly available database containing information on the swordfish transcriptome aims to boost research on this species with the long-term of developing more comprehensive and successful stock management plans.

Molecular and functional characterization of the vitellogenin receptor in oriental river prawn, Macrobrachium nipponense

A complementary DNA (cDNA) that encodes the vitellogenin receptor (VgR) in the oriental river prawn, Macrobrachium nipponense, was cloned using expressed sequence tag analysis and a rapid amplification of cDNA ends approach. The coding region consists of 5920 base pairs (bp) that encode a 1902 amino acid protein, with a predicted molecular mass of 209 kDa. The coding region is flanked by a 45 bp 5'-untranslated region (UTR) and a 166 bp 3'-UTR. The deduced amino acid sequence of the M. nipponense VgR cDNA had typically conserved domains, such as an extracellular, lipoprotein-binding domain, epidermal growth factor-like and O-glycosylation domains, a transmembrane domain and a short C-terminal, cytosolic tail. Quantitative real-time PCR (qPCR) indicated that Mn-VgR is highly expressed in the female ovary. Expression analysis by qPCR demonstrated the larval and ovarian developmental stage-specific expression pattern. As the ovaries developed, the expression level of Mn-VgR gradually increased during the reproductive cycle (stage I), to reach a peak in stage III. Levels then dropped as a new development cycle was entered after reproduction molting. Eyestalk ablation led to a significant increase in the expression of Mn-VgR during the ovarian development stages (P<0.05), when compared with the eyestalk-intact group. The investigation revealed that eyestalk ablation initially affected Mn-VgR expression and then influenced vitellogenesis. In adult females, VgR RNA interference (RNAi) dramatically delayed the maturation of the ovary, in accordance with the gonad somatic index. In addition, Mn-VgR RNAi led to vitellin depletion in the oocytes and the accumulation of vitellin in the hepatopancreas.

The Essential Role of Vitellogenin Receptor in Ovary Development and Vitellogenin Uptake in Bactrocera dorsalis (Hendel)

The vitellogenin receptor (VgR) functions as an essential component in uptaking and transporting vitellogenin (Vg) in female adults, which is involved in ovary development and oviposition. This study aimed to clarify the molecular characteristics and function of VgR in the oriental fruit fly Bactrocera dorsalis (Hendel). Here, we identified the full-length of BdVgR (GenBank Accession No. JX469118), encoding a 1925 residue (aa) protein with a 214.72 kDa molecular mass and several typical motifs of low-density lipoprotein receptor superfamily (LDLR). Phylogenic analysis suggested that BdVgR was evolutionary conserved with other Dipteran VgRs. The expression of BdVgR was exclusively detected in the ovaries rather than head, thorax or other tissues. The developmental expression patterns showed that the signal of BdVgR was detectable in very beginning of adult stage, and positively correlated with the growth rate of ovaries and the expression levels of its ligands. In addition, we also demonstrated that the expression level of BdVgR, and ovary development were significantly suppressed after being injected with BdVgR-targeted dsRNA. Together, all of these results indicated that BdVgR was critical for yolk protein absorption and ovary maturation in B. dorsalis, playing a vital role in female reproduction.

Construction, De-Novo Assembly and Analysis of Transcriptome for Identification of Reproduction-Related Genes and Pathways from Rohu, Labeo rohita (Hamilton)

Rohu is a leading candidate species for freshwater aquaculture in South-East Asia. Unlike common carp the monsoon breeding habit of rohu restricts its seed production beyond season indicating strong genetic control over spawning. Genetic information is limited in this regard. The problem is exacerbated by the lack of genomic-resources. We identified 182 reproduction-related genes previously by Sanger-sequencing which were less to address the issue of seasonal spawning behaviour of this important carp. Therefore, the present work was taken up to generate transcriptome profile by mRNAseq. 16 GB, 72 bp paired end (PE) data was generated from the pooled-RNA of twelve-tissues from pre-spawning rohu using IlluminaGA-II-platform. There were 64.97 million high-quality reads producing 62,283 contigs and 88,612 numbers of transcripts using velvet and oases programs, respectively. Gene ontology annotation identified 940 reproduction-related genes consisting of 184 mainly associated with reproduction, 223 related to hormone-activity and receptor-binding, 178 receptor-activity and 355 embryonic-development related-proteins. The important reproduction-relevant pathways found in KEGG analysis were GnRH-signaling, oocyte-meiosis, steroid-biosynthesis, steroid-hormone biosynthesis, progesterone-mediated oocyte-maturation, retinol-metabolism, neuroactive-ligand-receptor interaction, neurotrophin-signaling and photo-transduction. Twenty nine simple sequence repeat containing sequences were also found out of which 12 repeat loci were polymorphic with mean expected-&-observed heterozygosity of 0.471 and 0.983 respectively. Quantitative RT-PCR analyses of 13-known and 6-unknown transcripts revealed differences in expression level between preparatory and post-spawning phase. These transcriptomic sequences have significantly increased the genetic-&-genomic resources for reproduction-research in Labeo rohita.

Fin spine bone resorption in atlantic bluefin tuna, Thunnus thynnus, and comparison between wild and captive-reared specimens

Bone resorption in the first spine of the first dorsal fin of Atlantic bluefin tuna (ABFT) has long been considered for age estimation studies. In the present paper spine bone resorption was assessed in wild (aged 1 to 13 years) and captive-reared (aged 2 to 11 years) ABFT sampled from the Mediterranean Sea. Total surface (TS), solid surface (SS) and reabsorbed surface (RS) were measured in spine transverse sections in order to obtain proportions of SS and RS. The spine section surface was found to be isometrically correlated to the fish fork length by a power equation. The fraction of solid spine bone progressively decreased according to a logarithmic equation correlating SS/TS to both fish size and age. The values ranged from 57% in the smallest examined individuals to 37% in the largest specimens. This phenomenon was further enhanced in captive-reared ABFT where SS/TS was 22% in the largest measured specimen. The difference between the fraction of SS of wild and captive-reared ABFT was highly significant. In each year class from 1- to 7-year-old wild specimens, the fraction of spine reabsorbed surface was significantly higher in specimens collected from March to May than in those sampled during the rest of the year. In 4-year-old fish the normal SS increase during the summer did not occur, possibly coinciding with their first sexual maturity. According to the correlations between SS/TS and age, the rate of spine bone resorption was significantly higher, even almost double, in captive-reared specimens. This could be attributed to the wider context of systemic dysfunctions occurring in reared ABFT, and may be related to a number of factors, including nutritional deficiencies, alteration of endocrine profile, cortisol-induced stress, and loss of spine functions during locomotion in rearing conditions.