Next Article in Journal
Reproductive and Developmental Biology of Acroclisoides sinicus, a Hyperparasitoid of Scelionid Parasitoids
Previous Article in Journal
Changes Occurring on the Activity of Salivary Alpha-Amylase Proteoforms in Two Naturalistic Situations Using a Spectrophotometric Assay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Markers of Genetic Variation in Blue Gourami (Trichogaster trichopterus) as a Model for Labyrinth Fish

by
Gad Degani
1,2,
Isana Veksler-Lublinsky
3 and
Ari Meerson
1,2,*
1
MIGAL–Galilee Research Institute, P.O.B. 831, Kiryat Shmona 1101602, Israel
2
Faculty of Sciences, Tel-Hai Academic College, Upper Galilee 1220800, Israel
3
Department of Software and Information Systems Engineering, Ben Gurion University of the Negev, P.O.B. 653, Beer-Sheva 8410501, Israel
*
Author to whom correspondence should be addressed.
Biology 2021, 10(3), 228; https://doi.org/10.3390/biology10030228
Submission received: 9 February 2021 / Revised: 5 March 2021 / Accepted: 12 March 2021 / Published: 16 March 2021
(This article belongs to the Section Zoology)

Abstract

:

Simple Summary

This review is a summary of recent studies of genes, many of them involved in growth and reproduction, which can be used for distinguishing between species of the Anabantoidei suborder of fish, focusing on the Blue Gourami as a model species. This is important in both basic science and aquaculture applications.

Abstract

Markers of genetic variation between species are important for both applied and basic research. Here, various genes of the blue gourami (Trichogaster trichopterus, suborder Anabantoidei, a model labyrinth fish), many of them involved in growth and reproduction, are reviewed as markers of genetic variation. The genes encoding the following hormones are described: kisspeptins 1 and 2, gonadotropin-releasing hormones 1, 2, and 3, growth hormone, somatolactin, prolactin, follicle- stimulating hormone and luteinizing hormone, as well as mitochondrial genes encoding cytochrome b and 12S rRNA. Genetic markers in blue gourami, representing the suborder Anabantoidei, differ from those in other bony fishes. The sequence of the mitochondrial cytochrome c oxidase subunit 1 (COI) gene of blue gourami is often used to study the Anabantoidei suborder. Among the genes involved in controlling growth and reproduction, the most suitable genetic markers for distinguishing between species of the Anabantoidei have functions in the hypothalamic–pituitary–somatotropic axis: pituitary adenylate cyclase-activating polypeptide and growth hormone, and the 12S rRNA gene.

1. Introduction

Genetic variability between the organisms refers to sequence differences between their genomes, part of which is reflected in the sequence of protein-coding genes. This variation in DNA sequence can be used as a marker to distinguish between organisms, including fishes, at all systemic levels [1,2,3,4,5,6,7,8]. Compared to the much-studied genetic markers in some fish species of economic value, little research has been done with the group of labyrinth fishes. Blue gourami (Trichogaster trichopterus) is a labyrinth fish of the suborder Anabantoidei, order Perciformes. It is a small tropical freshwater fish characterized by the presence of a chamber, or labyrinth, above the gills for the uptake of oxygen from the air for breathing. The labyrinth enables adaptation to life in water with low dissolved oxygen by partial air breathing. Anabantoid fishes are geographically distributed in central Africa, India, and southern Asia [9] (Figure 1). In their natural habitat, they adapt to an unpredictable environment in which water-dissolved oxygen concentration varies throughout the year and can reach very low concentration [10]. The 16 known genera of anabantoid fishes contain about 80 species (FishBase, Nelson et al. Fishes of the World [10]). However, the systematic characters of the labyrinth fishes have not been agreed upon and many synonyms are used. According to Vierke [10], taxonomists classify the labyrinth fishes into four families: Anabantidae (genera Sandelia, Ctenopoma, Anabans), Belontiidae (genera Trichopsis, Trichogaster, Sphaerichthys, Pseudosphromenus, Parosphromenus, Malpulutta, Helostoma, Ctenops, Colisa, Betta, Belontia), Osphronemidae (genus Osphronemus), and Helostomatidae (genus Helostoma).
Labyrinth fishes have two different stages in their life cycle: (i) before labyrinth organ development, from eggs to juveniles, when oxygen is absorbed over the entire surface by diffusion; and (ii) after labyrinth development, when the organ becomes important for breathing [10]. The adaptation to the development of eggs and fry in water with low oxygen concentration involves laying eggs in a bubble nest which supplies oxygen to the eggs and larvae [10]. In natural habitats, when there is a low density of mature males, they become territorial by building a bubble nest and protecting it from other males (Figure 2) [12,13].
After courting and fertilization, the female swims under the bubble nest and spawns eggs into it. The male guards the nest with the eggs. If an egg falls out, the male returns it to the nest. The male also protects the young fish, immediately after hatching (Figure 3).
Ornamental fish populations in their native habitats, mainly in tropical areas, have declined due to overfishing for their sale on the tropical fish market. Fish of the suborder Anabantoidei are important in the ornamental fish industry and have long been produced in aquaculture [10]. The introduction of T. trichopterus in Florida is considered to have failed, according to US Fish and Wildlife Service, June 2019 [11].

2. Sequencing Analysis of the 12S rRNA and Cytochrome b Gene Variations in Blue Gourami

Information on the molecular variation of species belonging to the Anabantoidei in the order Perciformes (perch-like-fishes) is very limited [2,15,16]. Polymorphisms in several enzymes can be utilized as a genetic marker for these Anabantoidei species belonging to the Perciformes [5,17].
Thus, two Anabantoidei species—Trichogaster trichopterus and Trichogaster leerii—are similar in morphology and in their geographical distribution. Colisa lalia does not exhibit an overlapping geographical distribution, but shows a high similarity to Betta betta, which is distributed in the small area covered by the genus Trichogaster. The gene loci of species belonging to the genus Colisa showed a higher degree of similarity to those of the Trichogaster species than to those of the Betta species [5]. Population structures of Trichogaster pectoralis collected from five locations in Thailand [17] were also studied by isozyme analysis. The highest genetic identity coefficient was found between Samutprakan and Pitsanulok populations in Thailand; the lowest between Pitsanulok and Pattanee populations.
Several sets of degenerate oligonucleotides were used by Degani [2] to clone DNA fragments encoding portions of the cytochrome b and 12S rRNA genes. These genes were used to examine the genetic variation between species of the Anabantoidei by mitochondrial gene-sequencing analysis. Results demonstrated a similarity between the gene sequences of the various Belontiidae species, leading to the finding that these genes could serve as molecular markers for the systematic study of Belontiidae reproduction. The cytochrome b sequences of Anabantoidei fish examined in that study are shown in Figure 4.
Based on this gene, the most similar species were Trichogaster trichopterus (gold) and Trichogaster trichopterus (blue) (100%). Trichogaster leerii was less similar to these (86.0%), and even lower similarity was found between the species T. trichopterus and Trichogaster labiosus (85.6%). The least similarity was observed between Betta betta and the genera Colisa (50.2%) and Trichogaster (60.1%). The similarity value for the cytochrome b gene of Macropodus opercularis was between that of Betta and Colisa [2]. The phylogenetic tree for cytochrome b is shown in Figure 5 [2].
Additionally, sequences of the rRNA 12S gene from T. trichopterus (gold), T. trichopterus (blue), T. leerii, C. lalia, and B. betta are shown in Figure 6, and the nucleotide-similarity phylogenetic tree is presented in Figure 7 [2]. The phylogenetic tree results for the 12S rRNA gene were very similar to those for cytochrome b. 12S rRNA gene sequence similarity between Trichogaster species was high (91.4%–100%), and there was less similarity between this genus and C. lalia (88.4%). In addition, there was less similarity between T. trichopterus and C. lalia than between T. trichopterus and B. betta (82.6% and 84.0%, respectively) [2].
The cytochrome b and 12S rRNA gene sequences of the Anabantoid fishes were compared to those of other fishes and presented as a phylogenetic trees of nucleotide similarity (Figure 8 and Figure 9). In this comparison, some fishes that do not belong to the Anabantoidei seemed to have high similarity to the Anabantoid fishes.

3. Hypothalamus–Pituitary-Gonad (HPG) Axis Genes as Molecular Markers in Blue Gourami and Other Anabantoid Species

The genes involved in the HPG are important in the breeding of fishes and therefore their sequences have been studied. In teleosts, as in other vertebrates, kisspeptin (Kiss) has recently received considerable attention as a potential key player in the indirect neuroendocrine control of reproduction [18]. Kiss is a member of the RFamide peptide family. Originally identified as a metastasis suppressor in mammals, the Kiss1 gene produces several Kiss peptides in mammals. Kiss54 and its endogenous variants, Kiss14, Kiss13, and Kiss10, are generated by proteolytic cleavage of the Kiss precursor derived from the Kiss1 gene. They share a common core 10-amino acid (aa) sequence (Kiss10) at their C-terminal end, which allows them to bind to their cognate G-protein-coupled receptor (GPR54) or Kiss receptor (Kiss1r) [19]. Thus, Kiss1 controls the HPG axis, acts on the caudal hypothalamus and seems to affect receptors of gonadotropin-releasing hormone (GnRH) [20]. It controls the release of the pituitary gonadotropins follicle-stimulating hormone (FSH) and luteinizing hormone (LH), which in turn control gametogenesis [21] through steroids [21,22,23,24]. Studies on Kiss in teleosts have shown variation in their involvement in reproduction; however, more detailed studies are required due to the relatively large size of this class and the inherently high variation in hormones involved in reproduction. The brain of zebrafish (Danio rerio), one of the most intensively studied model fish, has two Kiss genes, Kiss1 and Kiss2, and two Kiss receptors (GPR54): Kiss1r and Kiss2r [19,25]; this is similar to other fishes such as lamprey (Petromyzon marinus) [25], medaka (Oryzias latipes) [26], goldfish (Carassius auratus) [27] and striped bass (Morone saxatilis) [28]. We studied the DNA sequence of the brain Kiss2 and the two Kiss receptors (GPR54), Kiss1r and Kiss2r, in blue gourami [29,30]. The obtained partial sequences of Kiss2 were compared with homologous sequences from a number of other fish species (Figure 10).
There was a low degree of similarity for both nucleotide and amino acid sequences of Kiss2 between the blue gourami and other fish species (Figure 10). A higher degree of similarity was found between the Kiss2r of blue gourami and those of the other fish species (Figure 11). The results show that Trichogster trichopterus Kiss2r’s amino acids sequence is very different from those of the species Tetraodon nigroviridis, Xenopus tropicalis, Scomber japonicus, and Gasterosteus aculeatus and closer to the species Danio rerio, Carassius auratus, and Dicentrarchus labrax.
The obtained sequences of Kiss1r in the blue gourami were also compared with homologous sequences in a number of other fish species. There was a low degree of similarity for both nucleotide and amino acid sequences between the blue gourami sequence and those of other fish species (Figure 12).
GnRH plays a central role in the control of vertebrate reproduction by affecting various other hormones, e.g., gonadotropins and the growth hormone (GH) family, which in turn regulate gametogenesis and steroidogenesis [31,32]. cDNA cloning of three GnRH forms (GnRH1, GnRH2, and GnRH3) of blue gourami was described by Levy et al. [31,32,33] (Figure 13 and Figure 14).
As shown in Table 1, the degree of identity between the blue gourami GnRH3 preprohormone amino acid sequence and that of other fishes was 80–82.2% [31]. There was a higher degree of identity at both the nucleotide and amino acid levels between the blue gourami sequence and the respective sequences of Dicentrarchus labrax and Pagrus major (Order Perciformes). In other perciforms (Cynoscion nebulosus, Rachycentron canadum, Sciaenops ocellatus, and Mugil cephalus), a higher degree of amino acid similarity was found, as compared to the more variable nucleotide sequence [31,32,33].

4. DNA Sequences of FSH and LH as Molecular for Genetic Similarity between Blue Gourami and Other Fish Species

FSH and LH are gonadotropins that control gametogenesis in fishes. The sequences of those genes are described in many fishes. FSH and LH cDNAs from the pituitary gland of blue gourami (T. trichopterus), encoding the α and β subunits of these hormones, were cloned [34]. The two cDNAs were sequenced and analyzed. The deduced amino acid sequences of both FSH and LH cDNAs are presented in Figure 15.
Comparison of FSH-β and LH-β are shown in Figure 16. The blue gourami FSH-β was most similar to its striped bass counterpart, with the two polypeptides sharing 73% of their residues. The lowest similarity was found between blue gourami and of goldfish FSH-β, with only 44% similarity (Figure 17A) [34]. A dendrogram, graphically representing these polypeptide comparisons, showed highest similarity between blue gourami and striped bass FSH-β, with 84% identical residues (Figure 17A). For LH-β, the lowest similarity was found with the Coregonus migratorius, a whitefish of the Family Salmonidae. (only 65% identical residues). The dendrogram showing the relationships between the LH-β polypeptides is presented in Figure 17B [34].

5. Hypothalamic–Pituitary–Somatotropic (HPS) Axis Genes as Markers for Genetic Variation between Blue Gourami and Other Fish Species

Pituitary adenylyl cyclase-activating polypeptide (PACAP) and PACAP-related peptide (PRP) are members of the secretin/glucagon/vasoactive intestinal polypeptide family. PACAP was first isolated from the ovine hypothalamus [36], and its sequence has also been determined in representative species of non-mammalian vertebrates. In fishes, PRP (formerly known as GHRH or GHRH-like peptide) was first isolated from the hypothalamus of common carp [37]. In teleosts, the complete sequence of the PRP–PACAP cDNA has been cloned from several species [36]. In all vertebrates, the PACAP preprohormone consists of a signal peptide, a cryptic peptide, a PRP segment and a PACAP segment. Two processing sites, Lys–Arg and Gly–Arg–Arg, between the PRP and PACAP fragments and in the PACAP sequence, respectively, result in three peptides: PRP, PACAP27 and PACAP38. In several teleosts, there are two transcripts of PRP–PACAP as a result of exon skipping. This stems from alternative splicing resulting in the partial excision of exon 4, which encodes part of PRP (residues 1–32), leaving the PACAP-encoding region in the correct reading frame [38]. PACAP can bind specifically to three G-protein-coupled receptor (GPCR) variants (PACAP-Rs) respectively termed: PAC1-R, VPAC1-R, and VPAC2-R [36]. Differential distribution of PACAP-R has been identified in the brain, pituitary, heart, spleen, liver, gut, gills, kidney, skin, blood, and gonads [36]. PRP-R also belongs to the GPCR family, but has only been identified in non-mammalian species, e.g., in the goldfish pituitary [39]. These findings suggest that evolutionary pressure has acted to preserve the bioactive sequence of PACAP, indicating that the peptide must exert important physiological functions. PACAP and PRP are involved in growth promotion and GH control. In teleosts, PACAP is involved in various physiological processes, such as brain development, ventilation and cardiac baroreflex control, digestive physiology, immune response, food intake, and growth promotion [40].
The structures of the PRP–PACAP cDNA are presented in Figure 18A [36], and the full-length cDNA sequence of the PRP–PACAP precursor, compiled from data obtained from 5’ and 3’ RACE, and its deduced amino acid sequence, are shown in Figure 18B [36]. The nucleotide and amino acid sequences of this gene in blue gourami were compared with homologous sequences from a number of other fish species (Figure 19) [36]. Sequence alignment of the PRP amino acid sequence revealed that only the first 32 aa at the N terminus are highly conserved in closely related fishes. There was 80.0–93.3% sequence identity between the blue gourami PRP sequence and that of other teleosts, whereas there was only 31.6–56.6% sequence identity between other tetrapods and fishes (Figure 19A). On the other hand, the amino acid sequence of PACAP is highly conserved in mammals and fishes.

6. GH and Prolactin (PRL) Family Hormones as Genetic Variation Markers for Blue Gourami and Other Anabantoid Fishes

Growth hormone (GH) is a single-chain polypeptide. Together with prolactin, growth hormone and somatolactin, it forms a family of related polypeptide hormones whose sequences seem to have evolved from a common ancestor [31,33,36]. GH has been studied extensively, and cDNA sequences are available for blue gourami (Figure 20) and many other teleosts. The study of GH expression in relation to growth and the reproductive cycle could contribute to an understanding of the interactions between somatotropic and gonadotropic axes at the pituitary level in order to elucidate the effects of GH on fish reproduction [31,33,36]. Other HPS axis genes that control growth and that are involved in cell division have been sequenced in blue gourami but less in other species of Anabantoid fishes [7,31,36,41,42].
In a comparison of the deduced amino acid sequence of bgGH to other fish sequences found in the GenBank database, bgGH was most similar to that of Sparus aurata (86% identical residues) and least similar to Anguilla japonica (43% identity) (Figure 21). The sequence similarity was in accordance with prevailing systematics [41].
Prolactin is a hormone involved in a large number of biological processes, including growth, reproduction and osmoregulation in fish. The complete cDNA of the blue gourami PRL (bgPRL) was cloned by RACE PCR [43]. The deduced amino acid sequence of bgPRL was compared with that of homologous subunits from a number of other fish species (Table 2). The highest degree of homology was between the gourami and the Perciforms Perca flavescens, Dicentrarchus labrax, and Sparus aurata, followed by members of the order Salmoniformes (66%), Siluriformes, and Cypriniformes (Table 2) [43]. The lowest level of homology was observed with the Anguilliformes (Table 2) [43].
Somatolactin (SL), a specific pituitary hormone belonging to the PRL superfamily, is involved in environmental adaptation, osmoregulation, reproduction, and fatty acid metabolism. The cDNA of SL from blue gourami was cloned and subjected to DNA sequence analysis [44,45]. The partial cDNA sequence of SL was compared to those of other fish species (Figure 22).

7. Mitochondrial Cytochrome c Oxidase Subunit 1 (COI) Gene as a Variation Marker for Blue Gourami and Other Anabantoid Fishes

COI is another genetic marker widely used to identify the kinship of fish species (Figure 24 and Figure 25).

8. Discussion

This review presents sequence analysis of genes involved in various growth and reproductive processes that may serve as markers of genetic variation between labyrinth fishes of the suborder Anabantoidei and other fish species and groups (Summary, Table 3). The various molecular markers in blue gourami are mitochondrial or nuclear DNA sequences, but not microsatellites or single-nucleotide polymorphisms. In the present review, most of the genetic markers were associated with growth (PACAP, GH and PRL) and reproduction (GnRH1, GnRH2, GnRH3, LH, and FSH). Only two mitochondrial genes (cytochrome b and 12S) were studied. Isozyme markers also have been studied in labyrinth fishes [5]. Thus, one can compare genomic markers with mitochondrial ones. The similarity of cytochrome b sequence between labyrinth and other fish species was between 86% and 46%, and of the 12S rRNA gene, between 91% and 40%. Sequence similarities between the blue gourami and other fish species for the hormone-encoding genes seemed to be lower than for the mitochondrial genes (Table 3). Among the genes involved in controlling growth and reproduction, those with the highest sequence similarity between species were those from the HPS axis. The findings based on DNA sequence comparison presented in this review are in agreement with many other studies [1,2,3,4,5,6,7,8,48,49]. All of the tested genes in blue gourami had high similarity to their counterparts in other fishes in the order Perciformes, to which the blue gourami belongs [10], and some of them can be useful as genetic markers in other classes of fish.
For example, the cytochrome b sequence was used to study genetic variation of the genus Capoeta and specifically the species C. damascina [49,50]. C. damascina is one of the most common freshwater fish species found in a wide range of isolated water bodies throughout the Levant, Mesopotamia, Turkey, and Iran. Prior to these studies, C. damascina was not a well-defined species. Likewise, the COI gene is widely used in assessing genetic variation in both geographical distribution and aquaculture contexts [51,52]. This gene is suitable for the separation of different species belonging to the genus Trichogaster, but not of the order Anabantiformes (Figure 23 and Figure 24).

Author Contributions

Conceptualization, G.D.; methodology, G.D. and I.V.-L.; software, I.V.-L.; formal analysis, G.D. and I.V.-L.; investigation, G.D.; writing—original draft preparation, G.D.; writing—review and editing, A.M.; visualization, G.D. and I.V.-L.; supervision, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the students who contributed directly or indirectly to this paper: Gal Levy, Karen Jackson, Svetlana Yom Din, Tzchori Itai, David Kushnirov, Avshalom Hurvitz, Oren Pearlson, and Tali Goldberg; colleagues at MIGAL and Tel-Hai Academic College: Dani Bercovich, Doron Goldberg; past and current collaborators on the study of fish endocrinology: Berta Levavi-Sivan, Yoav Gothilf, and Yonathan Zohar.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maqsood, H.M.; Ahmad, S.M. Advances in molecular markers and their applications in aquaculture and fisheries. Genet. Aquat. Organ. 2017, 1, 27–41. [Google Scholar]
  2. Degani, G. Mitochondrial DNA sequence analysis in Anabantoidei fish. Adv. Biol. Chem. 2013, 3, 347–355. [Google Scholar] [CrossRef] [Green Version]
  3. Okumus, F.; Ciftci, Y. Fish population genetics and molecular markers II molecular markers and their applications in fisheries and aquaculture. Turk. J. Fishe. 2020, 51–79. [Google Scholar]
  4. Degani, G.; Alon, A.; Hajouj, A.; Meerson, A. Vitellogenesis in the blue gourami (Trichogaster trichopterus) is accompanied by changes in the brain transcriptome. Fishes 2019, 4, 54. [Google Scholar]
  5. Degani, G.; Veit, M. Electrophoretic variations of isozyme systems in the muscle and liver of Anabantidae fish. Isr. J. Aquacul. 1990, 42, 67–76. [Google Scholar]
  6. Chauhan, T.; Rajiv, K. Molecular markers and their applications in fisheries and aquaculture Advan. Biosc. Biotechn. 2010, 1, 281–291. [Google Scholar]
  7. Degani, G.; Jackson, K.; Yom-Din, S.; Goldberg, D. cDNA cloning and mRNA expression of growth hormone in Belontiidae (Anabantoidei suborder) fish. Isr. J. Aquacul. 2006, 58, 124–136. [Google Scholar]
  8. Degani, G.; Hurvitz, A.; Eliraz, Y.; Meerson, A. Sex-related gonadal gene expression differences in the Russian sturgeon (Acipenser gueldenstaedtii) grown in stable aquaculture conditions. Anim. Reprod. Sci. 2019, 199, 75–85. [Google Scholar] [CrossRef]
  9. Forselius, S. Studies of anabantid fishes. Parts I, II, III. Zoologiska Bidrag Fran Uppsala 1975, 32, 597. [Google Scholar]
  10. Vierke, J. Betta Gouramis and Other Anabantoid labyrinth Fishes of the World; T.F.H. Publication: Neptune, NJ, USA, 1988; pp. 1–102. [Google Scholar]
  11. Froese, R.; Pauly, D. Trichopodus trichopterus (Pallas, 1770). FishBase 2019. Trichogaster labiosus distribution U.S. Fish and Wildlife Service. Available online: http://www.fishbase.org/summary/Trichopodus-trichopterus (accessed on 29 February 2021).
  12. Degani, G. The effect of temperature, light, fish size and container size on breeding of Trichogaster Trichopterus Isr. J. Aquacul. 1989, 41, 67–73. [Google Scholar]
  13. Degani, G.; Ziv, Z. Male blue gourami (Trichogaster trichopterus) nest-building behavior is affected.by other males and females. Open J. Anim. Sci. 2016, 6, 195–201. [Google Scholar] [CrossRef] [Green Version]
  14. Degani, G. The effect of sexual behavior on oocyte development and steroid changes in Trichogaster Trichopterus. Copeia 1993, 4, 1091–1096. [Google Scholar] [CrossRef]
  15. Ruber, L.; Britz, R.; Zarddya, R. Molecular phylogenetics and evolutionary diversification of labyrinth fishes (Perciformes: Anabantoidei). Syst. Biol. 2006, 55, 374–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Degani, G. Blue Gourami (Trichogaster trichopterus) Model for Labyrinth Fish; Laser Pages Publishing: Jerusalem, Israel, 2001; pp. 1–134. [Google Scholar]
  17. Pongthana, N. Aquaculture genetics research in Thailand. In Fish Genetics Research in Member Countries and Institutions of the International Network on Genetics in Aquaculture; Gupta, M.V., Acosta, B.O., Eds.; The World Fish Center, Penang: Bayan Lepas, Malaysia, 2001; pp. 77–89. [Google Scholar]
  18. Oakley, A.E.; Clifton, D.K.; Steiner, R.A. Kisspeptin signaling in the brain. Endoc. Rev. Esp. Herp. 2009, 30, 713–743. [Google Scholar] [CrossRef] [PubMed]
  19. Servili, A.; Le Page, Y.; Leprince, J.; Caraty, A.; Escobar, S.; Parhar, I.S.; Seong, J.Y.; Vaudry, H.; Olivier Kah, O. Organization of two Independent kisspeptin systems derived from evolutionary-ancient kiss genesin the brain of zebrafish. Endocrinology 2011, 152, 1527–1540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Shahjahan, M.; Motohashi, E.; Doi, H.; Hironori Ando, H. Elevation of Kiss2 and its receptor gene expression in the brain and pituitary of grass puffer during the spawning season. Gen. Comp. Endocrinol. 2010, 169, 48–57. [Google Scholar] [CrossRef] [PubMed]
  21. Yaron, Z.; Levavi-Sivan, B. Endocrine regulation of fish reproduction In Encyclopedia of Fish Physiology: From Genome to Environment; Farrell, A.P., Ed.; Elsevier Inc.: San Diego, CA, USA, 2011; Volume 2, pp. 1500–1508. [Google Scholar]
  22. Degani, G.; Boker, R. Sensitivity to maturation inducing steroids and gonadotropin in the oocytes of Blue gourami Trichogaster trichopterus, (Anabantidae, Pallas, 1770) Gen. Comp. Endinocrinol. 1992, 85, 430–439. [Google Scholar] [CrossRef]
  23. Degani, G.; Boker, R. Vitellogenesis level and the induction of maturation in the ovary of the blue gourami Trichogaster trichopterus (Anabantidae, Pallas 1770). J. Exp. Zool. 1992, 263, 330–337. [Google Scholar] [CrossRef]
  24. Degani, G. The effects of human chorionic gonadotropin on steroid changes in Trichogaster trichopterus. Comp. Biochem. Physiol. 1990, 96A, 525–528. [Google Scholar] [CrossRef]
  25. Lee, Y.R.; Tsunekawa, K.; Moon, M.J.; Um, H.N.; Hwang, J.I.; Osugi, T.; Otaki, N.; Sunakawa, Y.; Kim, K.; Vaudry, H.; et al. Molecular evolution of multiple forms of kisspeptins and GPR5 4 receptors in vertebrates. Endocrinology 2009, 150, 2837–2846. [Google Scholar] [CrossRef] [Green Version]
  26. Kitahashi, T.; Ogawa, S.; Parhar, I.S. Cloning and expression of kiss2 in the zebrafish and medaka. Endocrinology 2009, 150, 821–831. [Google Scholar] [CrossRef]
  27. Huang, T.; Huang, X.; LuD, W.; Zhu, P.; Shi, Y.; Cheng, C.H.; Liu, X.; Lin, H. Structural and functional multiplicity of the kisspeptin/GPR54 system in goldfish (Carassius auratus). 8. J. Endocrinol. 2009, 201, 407–441. [Google Scholar]
  28. Espigares, F.; Rocha, A.; Gómez, A.; Carrill, M.; Zanuy, S. Photoperiod modulates the reproductive axis of European sea bass through regulation of kiss1 and gnrh2 neuronal expression. Gen. Com. Endo. 2017, 240, 35–45. [Google Scholar] [CrossRef] [PubMed]
  29. Degani, G.; Alon, A.; Stoler, A.; Bercocvich, D. Evidence of a reproduction-related function for brine Kisspeptin 2 and its receptors in Anabantidae fish (Trichogaster trichopterus) Nternat. J. Zool. Invest. 2017, 2, 106–122. [Google Scholar]
  30. Sower, S.A.; Freamat, M.; Kavanaugh, S.I. The origins of the vertebrate hypothalamic-pituitary-gonadal (HPG) and hypothalamicpituitary-thyroid (HPT) endocrine systems: New insights from lampreys. Gen. Comp. Endocrinol. 2009, 161, 20–29. [Google Scholar] [CrossRef] [PubMed]
  31. Levy, G.; Degani, G. Involvement of GnRH, PACAP and PRP in the reproduction of blue gourami females (Trichogaster trichopterus). J. Mol. Neurosci. 2012, 48, 9730–9738. [Google Scholar] [CrossRef] [PubMed]
  32. Levy, G.; Gothilf, Y.; Degani, G. Brain gonadotropin releasing hormone3 expression variation during oogenesis and sexual behavior and its effect on pituitary hormonal expression in the blue gourami. Com. Biochem. Physio. Part A 2009, 154, 241–248. [Google Scholar] [CrossRef] [PubMed]
  33. Levy, G.; Degani, G. The role of brain peptides in the reproduction of blue gourami males (Trichogaster trichopterus). J. Exp. Zool. 2013, 9999A, 1–10. [Google Scholar] [CrossRef] [PubMed]
  34. Jackson, K.; Goldberg, G.; Ofir, R.; Abraham, M.; Degani, G. Blue gourami (Trichogaster trichopterus) gonadotropic subunits (I & II) cDNA sequences and expression during oogenesis J. Mol. Endocrinol. 1999, 23, 177–187. [Google Scholar]
  35. Yom Din, S.; Hurvitz, A.; Goldberg, D.; Jackson, K.; Levavi-Sivan, B.; Degani, G. Cloning of Russian sturgeon (Acipenser gueldenstaedtii) growth hormone and insulin-like growth factor 1 and their expression in male and female fish during the first period of growth. J. Endocrinol. Investig. 2008, 31, 201–210. [Google Scholar] [CrossRef]
  36. Levy, G.; Goldberg, D.; Jackson, K.; Degani, G. Association between pituitary adenylate cyclase activating polypeptide and reproduction in the blue gourami. Gen. Comp. Endo. 2010, 166, 83–93. [Google Scholar] [CrossRef]
  37. Vaughan, J.M.; Rivier, J.; Spiess, J.; Peng, C.; Chang, J.P.; Peter, R.E.; Vale, W. Isolation and characterization of hypothalamic growth-hormone releasing factor from common carp Cyprinus carpio. Neuroendocrinology 1992, 56, 539–549. [Google Scholar] [CrossRef] [PubMed]
  38. Small, B.C.; Nonneman, D. Sequence and expression of a cDNA encoding both pituitary adenylate cyclase activating polypeptide and growth hormone-releasing hormone-like peptide in channel catfish (Ictalurus punctatus). Gen. Comp. Endocrinol. 2001, 122, 354–363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Castro, A.; Becerra, M.; Manso, M.J.; Tello, J.; Sherwood, N.M.; Anadon, R. Distribution of growth hormone-releasing hormone-like peptide: Immunoreactivity in the central nervous system of the adult zebrafish (Danio rerio). J. Comp. Neurol. 2009, 513, 685–701. [Google Scholar] [CrossRef] [PubMed]
  40. Lancien, F.; Mimassi, N.; Conlon, J.M.; Le Mevel, J.C. Central pituitary adenylate cyclase-activating polypeptide ( PACAP) and vasoactive intestinal peptide (VIP) decrease the baroreflex sensitivity in trout. Gen. Comp. Endocrinol. 2011, 171, 245–251. [Google Scholar] [CrossRef] [PubMed]
  41. Goldberg, D.; Jackson, K.; Yom-Din, S.; Degani, G. Growth hormone of Trichogaster trichopterus: cDNA cloning, sequencing and analysis of mRNA expression during oogenesis. J. Aqua. Trop. 2004, 19, 215–229. [Google Scholar]
  42. Levy, G.; David, D.; Degani, G. Effect of environmental temperature on growth and reproduction-related hormone gene expression in the female blue gourami (Trichogaster trichopterus). Com. Biochem. Physiol. Part A 2011, 160, 381–389. [Google Scholar] [CrossRef]
  43. Degani, G.; Yom-Din, S.; Goldberg, D.; Jackson, K. cDNA cloning of blue gourami (Trichogaster trichopterus) prolactin and its expression during the gonadal cycles of males and females. J. Endocrinol. Investig. 2010, 33, 7–12. [Google Scholar] [CrossRef] [PubMed]
  44. Degani, G. Somatolactin transcription during oogenesis in female blue gourami (Trichogaster trichopterus). Advan. Biolog. Chemis. 2015, 5, 279–285. [Google Scholar] [CrossRef] [Green Version]
  45. Degani, G. Brain control reproduction by the endocrine system of female blue gourami (Trichogaster trichopterus). Biology 2020, 9, 109. [Google Scholar] [CrossRef]
  46. Davey, J.W.; Hohenlohe, P.A.; Etter, P.D.; Boone, J.Q.; Catchen, J.M.; Blaxter, M.L. Genome-wide genetic marker discovery and genotyping using next generation sequencing. Nat. Rev. Genet. 2011, 12, 499–510. [Google Scholar] [CrossRef] [PubMed]
  47. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K.T. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  48. Degani, G.; Yom-Din, S.; Hurvitz, A. Transcription of insulin-like growth factor receptor in russian sturgeon (Acipenser Gueldenstaedtii) ovary during oogenesis. Univers. J. Agric. Res. 2017, 5, 119–124. [Google Scholar] [CrossRef] [Green Version]
  49. Degani, G. DNA Variation of Capoeta damascina (Valenciennes, 1842) in three rivers in northern Israel. J. Biophys. Chem. 2014, 5, 107–117. [Google Scholar] [CrossRef]
  50. Levin, B.A.; Freyhof, J.; Lajbner, Z.; Perea, S.; Abdoli, A.; Gaffaroglu, M.; Ozulug, M.; Rubenyan, H.R.; Salnikov, V.B. and Doadrio, I. Phylogenetic relationships of the algae scraping cyprinid genus Capoeta (Teleostei Cyprinidae). Mol. Phylogenetics Evol. 2012, 62, 542–549. [Google Scholar] [CrossRef]
  51. Syaifudin, M.; Jubaedah, D.; Yonarta, D.; Hastuti, Z. DNA barcoding of snakeskin gourami Trichogaster pectoralis and blue gourami Trichogaster trichopterus based on cythocrome c oxidase subunit I (COI) gene. Earth. Enviro. Sci. 2019, 348, 012031. [Google Scholar] [CrossRef]
  52. Ratnasingham, S.; Hebert, P.D. Bold: The Barcode of Life Data System (http://www.barcodinglife.org). Mol. Ecol. Notes. 2007, 7, 355–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Known global distribution of Trichogaster trichopterus. Locations are in Australia, Papua- New Guinea, Indonesia, Malaysia, the Philippines, Taiwan, Vietnam, Cambodia, Laos, Thailand, Myanmar, India, Namibia, Colombia, Brazil, and the United States. Map from the GBIF Secretariat [11].
Figure 1. Known global distribution of Trichogaster trichopterus. Locations are in Australia, Papua- New Guinea, Indonesia, Malaysia, the Philippines, Taiwan, Vietnam, Cambodia, Laos, Thailand, Myanmar, India, Namibia, Colombia, Brazil, and the United States. Map from the GBIF Secretariat [11].
Biology 10 00228 g001
Figure 2. A scheme of two habitats, with high (A) and low (B) densities of males. At high density, the male does not build a nest. At low density, the male builds a nest and sexual behavior is initiated [10,13,14].
Figure 2. A scheme of two habitats, with high (A) and low (B) densities of males. At high density, the male does not build a nest. At low density, the male builds a nest and sexual behavior is initiated [10,13,14].
Biology 10 00228 g002
Figure 3. Sexual behavior of male blue gourami during the reproductive cycle. (1) The male builds a nest. (2) and (3) The male courts the female under the nest. (4) Fertilization takes place and the fertilized eggs float up and stick to the bubble nest. (5) The male guards the eggs in the nest. (6) The male further protects the fry immediately after hatching while they swim in the nest area [10,13].
Figure 3. Sexual behavior of male blue gourami during the reproductive cycle. (1) The male builds a nest. (2) and (3) The male courts the female under the nest. (4) Fertilization takes place and the fertilized eggs float up and stick to the bubble nest. (5) The male guards the eggs in the nest. (6) The male further protects the fry immediately after hatching while they swim in the nest area [10,13].
Biology 10 00228 g003
Figure 4. Comparison of cytochrome b nucleotide sequences (150 bp) in T. trichopterus (blue and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), T. labiosus (Colisa labiosa), B. betta (fighting fish), and M. opercularis [2].
Figure 4. Comparison of cytochrome b nucleotide sequences (150 bp) in T. trichopterus (blue and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), T. labiosus (Colisa labiosa), B. betta (fighting fish), and M. opercularis [2].
Biology 10 00228 g004
Figure 5. Composite phylogenetic tree constructed from the analysis of cytochrome b sequences of six species: T. trichopterus (blue gourami and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), T. labiosus (C. labiosa), B. betta (fighting fish), and M. opercularis. The phylogenetic tree was constructed by Clustal W and analysis alignment methods in the MegAlign program (DNASTAR) [2]. Branch lengths represent evolutionary distances. Percent DNA sequence identities are also shown [2].
Figure 5. Composite phylogenetic tree constructed from the analysis of cytochrome b sequences of six species: T. trichopterus (blue gourami and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), T. labiosus (C. labiosa), B. betta (fighting fish), and M. opercularis. The phylogenetic tree was constructed by Clustal W and analysis alignment methods in the MegAlign program (DNASTAR) [2]. Branch lengths represent evolutionary distances. Percent DNA sequence identities are also shown [2].
Biology 10 00228 g005
Figure 6. Comparison of the 12S rRNA gene nucleotide sequences (150 bp) in Trichogaster trichopterus (blue gourami and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), and B. betta (fighting fish) [2].
Figure 6. Comparison of the 12S rRNA gene nucleotide sequences (150 bp) in Trichogaster trichopterus (blue gourami and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), and B. betta (fighting fish) [2].
Biology 10 00228 g006
Figure 7. Phylogenetic tree constructed from analysis of 12S rRNA gene sequences of T. trichopterus (blue gourami and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), and B. betta (fighting fish). The phylogenetic tree was constructed by ClustalW and analysis alignment methods in the MegAlign program (DNASTAR). Branch lengths represent evolutionary distances. Percent DNA sequence identities are also shown [2].
Figure 7. Phylogenetic tree constructed from analysis of 12S rRNA gene sequences of T. trichopterus (blue gourami and gold gourami), T. leerii (pearl gourami), C. lalia (dwarf gourami), and B. betta (fighting fish). The phylogenetic tree was constructed by ClustalW and analysis alignment methods in the MegAlign program (DNASTAR). Branch lengths represent evolutionary distances. Percent DNA sequence identities are also shown [2].
Biology 10 00228 g007
Figure 8. Phylogenetic tree of the cytochrome b fragment based on nucleotide sequences of various fish species. The length of each pair of branches represents the distance between sequence pairs, while the units at the bottom of the tree indicate the number of substitution events. The phylogenetic tree was constructed by Clustal W and analysis alignment in the MegAlign program (DNASTAR). Percent DNA sequence identities are also shown [2].
Figure 8. Phylogenetic tree of the cytochrome b fragment based on nucleotide sequences of various fish species. The length of each pair of branches represents the distance between sequence pairs, while the units at the bottom of the tree indicate the number of substitution events. The phylogenetic tree was constructed by Clustal W and analysis alignment in the MegAlign program (DNASTAR). Percent DNA sequence identities are also shown [2].
Biology 10 00228 g008
Figure 9. Phylogenetic tree of the 12S rRNA gene fragment based on nucleotide sequences. The length of each pair of branches represents the distance between sequence pairs, while the units at the bottom of the tree indicate the number of substitution events. The phylogenetic tree was constructed by Clustal W and analysis alignment in the MegAlign program (DNASTAR). Percent DNA sequence identities are also shown [2].
Figure 9. Phylogenetic tree of the 12S rRNA gene fragment based on nucleotide sequences. The length of each pair of branches represents the distance between sequence pairs, while the units at the bottom of the tree indicate the number of substitution events. The phylogenetic tree was constructed by Clustal W and analysis alignment in the MegAlign program (DNASTAR). Percent DNA sequence identities are also shown [2].
Biology 10 00228 g009
Figure 10. (A): Phylogenetic tree showing the relationships between fish Kiss2 cDNA sequences. The tree was generated using MAFFT alignment in DNASTAR MegAlign PRO. (B): Nucleotide and deduced amino-acid sequences [29].
Figure 10. (A): Phylogenetic tree showing the relationships between fish Kiss2 cDNA sequences. The tree was generated using MAFFT alignment in DNASTAR MegAlign PRO. (B): Nucleotide and deduced amino-acid sequences [29].
Biology 10 00228 g010
Figure 11. Phylogenetic tree showing the relationship between vertebrate Kiss2r amino acid sequences. The tree was generated by maximum Clustal W alignment in the MegAlign program (DNASTAR). Sequence alignment was conducted using DNA Star MegAlign Pro Clustal Omega. The Trichogaster trichopterus sequence was obtained from [29]. All other sequences were obtained from NCBI GenBank.
Figure 11. Phylogenetic tree showing the relationship between vertebrate Kiss2r amino acid sequences. The tree was generated by maximum Clustal W alignment in the MegAlign program (DNASTAR). Sequence alignment was conducted using DNA Star MegAlign Pro Clustal Omega. The Trichogaster trichopterus sequence was obtained from [29]. All other sequences were obtained from NCBI GenBank.
Biology 10 00228 g011
Figure 12. Phylogenetic tree and DNA sequences showing the relationships between vertebrate Kiss1r cDNA sequences. The tree was generated using MAFFT alignment in DNASTAR MegAlign PRO [30]. T. trichopterus sequence was obtained from [29]. All other sequences were obtained from NCBI GenBank according to the accession numbers provided in [30]. The units of branch length represent nucleotide substitutions per site.
Figure 12. Phylogenetic tree and DNA sequences showing the relationships between vertebrate Kiss1r cDNA sequences. The tree was generated using MAFFT alignment in DNASTAR MegAlign PRO [30]. T. trichopterus sequence was obtained from [29]. All other sequences were obtained from NCBI GenBank according to the accession numbers provided in [30]. The units of branch length represent nucleotide substitutions per site.
Biology 10 00228 g012
Figure 13. The schematic structure of the prepro GnRH of blue gourami (A); cDNA sequences of GnRH1 (B) and GnRH2 (C). The signal peptide is underlined in black; the GnRH sequence is double-underlined, the cleavage-site position is underlined with a bold black line; the GAP sequence is underlined with a dashed line; and the 3’ poly(a) signal is boxed [31].
Figure 13. The schematic structure of the prepro GnRH of blue gourami (A); cDNA sequences of GnRH1 (B) and GnRH2 (C). The signal peptide is underlined in black; the GnRH sequence is double-underlined, the cleavage-site position is underlined with a bold black line; the GAP sequence is underlined with a dashed line; and the 3’ poly(a) signal is boxed [31].
Biology 10 00228 g013
Figure 14. Schematic structure of the Prepro GnRH3 (A) and the cDNA and deduced amino acid sequence of GnRH3 (B) of the blue gourami. The signal peptide is underlined in black; the GnRH3 sequence is double-underlined and the cleavage-site position is underlined with a bold black line; the GAP sequence is underlined with a dashed black line and the 3′ poly(a) signal is boxed [31].
Figure 14. Schematic structure of the Prepro GnRH3 (A) and the cDNA and deduced amino acid sequence of GnRH3 (B) of the blue gourami. The signal peptide is underlined in black; the GnRH3 sequence is double-underlined and the cleavage-site position is underlined with a bold black line; the GAP sequence is underlined with a dashed black line and the 3′ poly(a) signal is boxed [31].
Biology 10 00228 g014
Figure 15. Nucleotide sequences of blue gourami FSH-β (A) and LH-β (B). The amino acid sequence. of each hormone appears in three-letter code. The GenBank accession numbers of these sequences are AF157630 and AF157631, respectively [34].
Figure 15. Nucleotide sequences of blue gourami FSH-β (A) and LH-β (B). The amino acid sequence. of each hormone appears in three-letter code. The GenBank accession numbers of these sequences are AF157630 and AF157631, respectively [34].
Biology 10 00228 g015
Figure 16. A comparison of the amino acid sequences of FSH-β (gth1) and LH-β (gth2). Identical amino acids are. boxed in black; residues with similar physicochemical properties are boxed in gray [34].
Figure 16. A comparison of the amino acid sequences of FSH-β (gth1) and LH-β (gth2). Identical amino acids are. boxed in black; residues with similar physicochemical properties are boxed in gray [34].
Biology 10 00228 g016
Figure 17. Phylogenetic tree showing the relationships between fish FSH-β (A) and LH-β (B) sequences. The dendrograms were created by the UPGMA method from similarity matrices that were produced from the corresponding sequence alignments. The scale bar is a measure of the estimated number of amino acid substitutions per 100 residues in pairwise comparisons [34,35].
Figure 17. Phylogenetic tree showing the relationships between fish FSH-β (A) and LH-β (B) sequences. The dendrograms were created by the UPGMA method from similarity matrices that were produced from the corresponding sequence alignments. The scale bar is a measure of the estimated number of amino acid substitutions per 100 residues in pairwise comparisons [34,35].
Biology 10 00228 g017
Figure 18. Schematic structures and sequences of PRP–PACAP. (A). Schematic structures of the prepro PRP–PACAP cDNA of blue gourami. (B). cDNA sequence of PRP–PACAP. The signal peptide is underlined with a dashed black line; the cryptic peptide is underlined with a dashed double line; the PRP sequence is underlined in black; the cleavage-site positions are boxed and the PACAP sequence is double-underlined. The gray box designates the deleted region in the PRP–PACAP short form. The triangle denotes the 3’ poly(A) signal ATTAAA [31,33,36].
Figure 18. Schematic structures and sequences of PRP–PACAP. (A). Schematic structures of the prepro PRP–PACAP cDNA of blue gourami. (B). cDNA sequence of PRP–PACAP. The signal peptide is underlined with a dashed black line; the cryptic peptide is underlined with a dashed double line; the PRP sequence is underlined in black; the cleavage-site positions are boxed and the PACAP sequence is double-underlined. The gray box designates the deleted region in the PRP–PACAP short form. The triangle denotes the 3’ poly(A) signal ATTAAA [31,33,36].
Biology 10 00228 g018
Figure 19. PRP (A) and PACAP (B) amino acid sequence comparison between blue gourami and other vertebrates. The sequence alignment was conducted using ClustalW and the MegAlign program (DNASTAR) on PRP and PACAP mature peptide sequences. Identical amino acids are represented by dots, and differing residues are shown in the single letter code for amino acids. In contrast to PRP, PACAP is highly conserved [36].
Figure 19. PRP (A) and PACAP (B) amino acid sequence comparison between blue gourami and other vertebrates. The sequence alignment was conducted using ClustalW and the MegAlign program (DNASTAR) on PRP and PACAP mature peptide sequences. Identical amino acids are represented by dots, and differing residues are shown in the single letter code for amino acids. In contrast to PRP, PACAP is highly conserved [36].
Biology 10 00228 g019
Figure 20. Nucleotide sequence of blue gourami growth hormone (GH) cDNA together with the hormone’s putative amino acid sequence (in three letter code) [42].
Figure 20. Nucleotide sequence of blue gourami growth hormone (GH) cDNA together with the hormone’s putative amino acid sequence (in three letter code) [42].
Biology 10 00228 g020
Figure 21. Phylogenetic trees based on the comparison of the deduced amino acid sequence of bgGH to similar Belontiidae family (Anabantoidei suborder) fishes (A) and various fish species (B) GH found in the GenBank database [41]. Tree was constructed using Neighbor-joining method using Kimura 2 Distance Model, with pairwise deletion and MUSCLE alignment in the MEGA-X software [41].
Figure 21. Phylogenetic trees based on the comparison of the deduced amino acid sequence of bgGH to similar Belontiidae family (Anabantoidei suborder) fishes (A) and various fish species (B) GH found in the GenBank database [41]. Tree was constructed using Neighbor-joining method using Kimura 2 Distance Model, with pairwise deletion and MUSCLE alignment in the MEGA-X software [41].
Biology 10 00228 g021
Figure 22. Comparison of cDNA sequences of somatolactin (SL) from blue gourami and various other fish species [44,45]. The phylogenetic tree showing the relationship between SL amino acid sequences compared to those of different teleosts is presented in Figure 23 [45].
Figure 22. Comparison of cDNA sequences of somatolactin (SL) from blue gourami and various other fish species [44,45]. The phylogenetic tree showing the relationship between SL amino acid sequences compared to those of different teleosts is presented in Figure 23 [45].
Biology 10 00228 g022
Figure 23. Phylogenetic tree showing the relationship between SL amino acid sequences from blue gourami and other fish species. The tree was generated by maximum Clustal W alignment and the MegAlign program (DNASTAR). All sequences were obtained from NCBI GenBank [44,45]. The highest degree of homology for SL was between T. trichopterus and the Perciform Perca flavescens and Cyclopterus lumpus. The lowest level of homology was observed with Anguilla anguilla.
Figure 23. Phylogenetic tree showing the relationship between SL amino acid sequences from blue gourami and other fish species. The tree was generated by maximum Clustal W alignment and the MegAlign program (DNASTAR). All sequences were obtained from NCBI GenBank [44,45]. The highest degree of homology for SL was between T. trichopterus and the Perciform Perca flavescens and Cyclopterus lumpus. The lowest level of homology was observed with Anguilla anguilla.
Biology 10 00228 g023
Figure 24. Phylogenetic tree of the COI-5p gene records of the genus Trichogaster. Twelve nucleotide sequences from 8 species were retrieved from BoldSystem [46], such that every species is represented by sequences from unique bins. Tree was constructed using Neighbour-joining method using Kimura 2 Distance Model, with pairwise deletion and MUSCLE alignment in the MEGA-X software [47]. A sequence from the genus Ctenops was used as an outgroup. Leaves are marked with genbank_accession|Bold_bin_id|species. Sequences that belong to the species that appear more than once in the tree are underlined with the same color. Bins that appear more than once are underlined with a dotted line.
Figure 24. Phylogenetic tree of the COI-5p gene records of the genus Trichogaster. Twelve nucleotide sequences from 8 species were retrieved from BoldSystem [46], such that every species is represented by sequences from unique bins. Tree was constructed using Neighbour-joining method using Kimura 2 Distance Model, with pairwise deletion and MUSCLE alignment in the MEGA-X software [47]. A sequence from the genus Ctenops was used as an outgroup. Leaves are marked with genbank_accession|Bold_bin_id|species. Sequences that belong to the species that appear more than once in the tree are underlined with the same color. Bins that appear more than once are underlined with a dotted line.
Biology 10 00228 g024
Figure 25. Phylogenetic tree of the COI-5p gene records of the order Anabantiformes. Nineteen nucleotide sequences from different species were retrieved from BoldSystem [46], such that every genus is represented by one sequence. Tree was constructed using Neighbor-joining method using Kimura 2 Distance Model, with pairwise deletion and MUSCLE alignment in the MEGA-X software [47]. A sequence from the species Mastacembelus brachyrhinus that belongs to the order Synbranchiformes was used as an outgroup. Leaves are marked with genbank_accession |Bold_bin_id|species|genus|subfamily|family. Sequences that belong to the same taxonomic family are marked with the same color.
Figure 25. Phylogenetic tree of the COI-5p gene records of the order Anabantiformes. Nineteen nucleotide sequences from different species were retrieved from BoldSystem [46], such that every genus is represented by one sequence. Tree was constructed using Neighbor-joining method using Kimura 2 Distance Model, with pairwise deletion and MUSCLE alignment in the MEGA-X software [47]. A sequence from the species Mastacembelus brachyrhinus that belongs to the order Synbranchiformes was used as an outgroup. Leaves are marked with genbank_accession |Bold_bin_id|species|genus|subfamily|family. Sequences that belong to the same taxonomic family are marked with the same color.
Biology 10 00228 g025
Table 1. Nucleotide and amino acid similarities (%) between cDNA sequences of GnRH1 of blue gourami and other teleosts [31]. Notes: a Pejerrey-type GnRH sequence, b Whitefish-type GnRH sequence.
Table 1. Nucleotide and amino acid similarities (%) between cDNA sequences of GnRH1 of blue gourami and other teleosts [31]. Notes: a Pejerrey-type GnRH sequence, b Whitefish-type GnRH sequence.
OrderSpeciesAccession No.% Amino Acid Identity% Nucleotide Identity
Perciformes (Scombridae)Thunnus thynnusEU23950070.374.3
Perciformes (Serranidae)Epinephelus bruneusFJ38004764.872.3
Perciformes (Cichlidae)Oreochromis niloticusAB1016656768.2
Perciformes (Moronidae)Morone chrysopsDQ00067261.562.6
PleuronectiformesParalichthys olivaceusDQ07469350.551.4
Verasper variegatusHM13160062.670.5
AtherinomorphaOdontesthes bonariensisaAY7446893322.6
Fundulus heteroclitusAB30226565.960.8
Coregonus clupeaformisb 47.354.2
MugilomorphaMugil cephalusAY37345062.659
Table 2. Degree of homology between blue gourami (bp) PRL and PRL of other fish classes.
Table 2. Degree of homology between blue gourami (bp) PRL and PRL of other fish classes.
SpeciesClass/OrderAccession no.bgPRL (%)
Perca flavescensPerciformesAY33249179
Dientrarchus labraxX7872379
Spaurus aurataAF06054177
Paralichthys olivaceusAF047616 75
Onchorhynchus mykissSalmoniformesM2473866
Coregonus autummalisZ2311466
S. salarX8478766
Heteropneustes fossilisSiluriformesAF37265362
I. punctatusAF26799062
Hypoththalmichtys molitrixCypriniformesX6105262
Danio rerioAY13514961
Cyprinus carpioX1254161
A. japonicaAnguiliformesAY15800959
A. anguillaX6914959
Table 3. Degree of nucleotide sequence similarity (%) of blue gourami (Trichogaster trichopterus) to different species [2,32,35,36,41,43].
Table 3. Degree of nucleotide sequence similarity (%) of blue gourami (Trichogaster trichopterus) to different species [2,32,35,36,41,43].
Cytochrome bMitochondrial
RNA 12S Gene
Growth
Hormone
ProlactinPACAPGnRH1GnRH2GnRH3
Trichopterus trichopterus (gold) 100%Trichogaster trichopterus (gold) 100%Lates calcarifer
84%
Perca flavescens
79%
Oreochromis mossambicus
94%
Thunnus thynnus
74.3%
Epinephelus bruneus
78.5%
Dicentrar chus labrax
77.4%
Colisa lalia
86.6%
Trichogaster leeri
91.4%
Seriola dumerili
84%
Dicentrarchus labrax
79%
Gadus morhua
97.4%
Epinephelu sruneus
72.3%
Verasper variegatus
74.9%
Pagrus major
70.2%
Trichogaster leerii
86.0%
Colisa lalia
88.4%
Sparus aurata
83%
Sparus aurata
77%
Takifugu rubripes
97.4%
Verasper variegatus
70.5%
Paralichthys olivaceus
73.4%
Cynoscion nebulosus
58.1%
Trichogaster labiosus
85.6%
Betta betta
82.6%
Acanthopagrus butcheri
82%
Paralichthys olivaceus
77%
Haplocho misburtoni
97.4%
Oreochromis niloticus
68.2%
Thunnus thynnus
72.8%
Rachycentron canadum
57.7%
Macropodus opercularis
81.6%
Colisa chuna
41.4%
Oreochromis niloticus
79%
Onchorhynechus mykiss
66%
Epinephelu soioides
97.4%
Morone saxatilis
62.6%
Morone saxatilis
72.2%
Micropog onias undulates
57.5%
Cyprinus carpio
77.9%
Trichogaster pectoralis
40.5%
Morone saxatilis
79%
Coregonus autumnalis
66%
Fundulus heteroclitus
60.8%
Odontesthes bonariensis
65.2%
Sciaenops ocellatus
57.5%
Betta betta
58.9%
Cternopoma setherici
41.1%
Caranx delicatissimus
74%
Salmo salar
66%
Mugil cephalus
59%
Mugil cephalus
65%
Salmo salar
76.9%
Macropodus opercularis
41.1%
Oncorhynchus tshawytscha
69%
Heteropneustes fossilis
62%
Coregonus clupeaformis
54.2%
Fundulus heteroclitus
61.8%
Trichogaster fasciatus
46.2%
Osphronemus goramy
40.5%
Salmon salar
68%
Ictalurus punctatus
62%
Odontesthes bonariensis
22.6%
Coregonus clupeaformis
59.6%
Pangasius pangasius
64%
Hypophthalmi chthys molitrix
62%
Pangasinodon gigas
64%
Danio rerio
61%
Cyprinus carpio
64%
Cyprinus carpio
61%
Carassius auratus
63%
Anguilla japonica
59%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Degani, G.; Veksler-Lublinsky, I.; Meerson, A. Markers of Genetic Variation in Blue Gourami (Trichogaster trichopterus) as a Model for Labyrinth Fish. Biology 2021, 10, 228. https://doi.org/10.3390/biology10030228

AMA Style

Degani G, Veksler-Lublinsky I, Meerson A. Markers of Genetic Variation in Blue Gourami (Trichogaster trichopterus) as a Model for Labyrinth Fish. Biology. 2021; 10(3):228. https://doi.org/10.3390/biology10030228

Chicago/Turabian Style

Degani, Gad, Isana Veksler-Lublinsky, and Ari Meerson. 2021. "Markers of Genetic Variation in Blue Gourami (Trichogaster trichopterus) as a Model for Labyrinth Fish" Biology 10, no. 3: 228. https://doi.org/10.3390/biology10030228

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop