Skip to main content

rDNA and mtDNA analysis for the identification of genetic characters in the hybrid grouper derived from hybridization of Cromileptes altivelis (female) × Epinephelus lanceolatus (male)

Abstract

Background

Hybridization is a useful strategy to produce offspring with more desirable phenotypic characteristics than those of parents. The hybrid grouper derived from the cross of Cromileptes altivelis (♀, 2n = 48) with Epinephelus lanceolatus (♂, 2n = 48) exhibits improved growth compared with its female parent, which makes it valuable to aquaculture. However, the genetic traits of the hybrid grouper are poorly understood.

Results

The observations showed that the hybrid grouper was diploid (2n = 48) and displayed intermediate morphology with the parent's measurable characteristics. The ribosomal DNA (rDNA) and mitochondria DNA (mtDNA) were characterized at molecular and phylogenetic level. High similarity and low genetic distance of 5S rDNA and mtDNA sequences between the hybrid grouper and C. altivelis showed that the hybrid grouper had a closer genetic relationship with female parents. The reconstructed phylogenetic tree based on COI gene and D-loop region of mtDNA recovered that mtDNA was maternally inherited in the hybrid grouper. Additionally, the DNA methylation level of 5S rDNA intergenic spacers (IGS) sequence was tested in here. The results showed that the DNA methylation status of the hybrid grouper was significantly lower than that of C. altivelis.

Conclusion

Results of this study provide important data on the genetic characteristics of the hybrid derived from the cross of C. altivelis and E. lanceolatus, and contribute the knowledge of both evolution and marine fish breeding.

Peer Review reports

Introduction

Cromileptes altivelis (Serranidae: Cromileptes) and Epinephelus lanceolatus (Serranidae: Epinephelus) have been widely cultivated in China and Southeast Asian countries, possessing commercial importance fishing market species [1, 2]. The appearance of C. altivelis is distinct from fish of the subfamily Epinephelinae, because it has a concave slope on the back of head, hump-like bulge on the back and dark gray-brown body with black spots, so it can also be called “humpback grouper” [3]. Juvenile C. altivelis is kept as ornamental fish, while the adult is important to commercial mariculture. E. lanceolatus (also named as giant grouper) is regarded as a high-value species in fish markets due to its fast growing, large size and high nutritional value [4]. It is distinguished by small eye, wide interorbital area, body with irregular whitish blotches and numerous small black spots, and fins slightly yellowish with irregular blackish and whitish marking [5]. Groupers are usually considered to be protogynous hermaphrodites, and their sex type changes from female to male as they mature. Mature males are hard to obtain because of their huge body size and long time to maturity, so it is quite difficult to artificially reproduce them [6]. Hybridization is the most effective and widely used technique in the artificial breeding of grouper, which can produce offspring with more desirable phenotypic characteristics than those of the parent [7]. Many excellent offspring have been obtained using C. altivelis or E. lanceolatus as parent: C. altivelis (♀) × E. tukula (♂) [8], E. moara (♀) × E. lanceolatus (♂) [9], E. coioides (♀) × E. lanceolatus (♂) [10]. The hybrid grouper derived from the crossing of C. altivelis (♀) and E. lanceolatus (♂) exhibits the combined advantages of the parent species, with faster growth and improved taste when compared with wild groupers. However, the genetic traits in the hybrid grouper are poorly understood.

Morphological research relies heavily on statistics and classification methods using measurable, countable, and partly descriptive traits in fish [11]. Where parents differ physically from one another, morphological analysis is a useful tool in revealing the genetic relationship between hybrid progenies and parents [10]. However, it cannot effectively identify the genetic relationship between offspring and parents only by phenotypic characteristics. When combined with morphological analysis, molecular markers can greatly improve the accuracy of hybrid identification. In eukaryotes, the ribosomal DNA multigene family is organized in coding regions and non-transcribed spacer (NTS) regions. Because concerted and birth-and-death evolution mechanisms act simultaneously in rDNA, the coding regions show high conservation, while the NTS regions exhibit different rates of variation within and between taxa. This trait of rDNA sequence structure makes it a good molecular marker for the investigation of rapid evolutionary events [12,13,14]. Intergenic spacers (IGS) can separate these rDNA repeats, and present between two successive genes [15]. IGS transcripts play an important role in the epigenetic control of the rDNA locus [16]. DNA methylation, essentially the methylation of cytosine nucleotides, is the first identified epigenetic mechanism [17] and has been extensively studied. Accumulated studies suggested that DNA methylation was closely correlated with the heterosis in animals. Ou et al. have reported that the association of DNA methylation with the growth heterosis in snakehead fish [18]. Jiang et al. have found that DNA methylation can be involved in the heterosis formation in pig hybrids [19]. In the Pacific oyster Crassostrea gigas, the DNA methylation level is associated with the superior growth of the hybrid crosses [20].

Mitochondria DNA (mtDNA) is also valued for tracking the ancestry of breeds back hundreds of generations [21], because the organization in most fish mtDNA genomes is quite conserved [22]. Many researchers analyzed the mtDNA D-loop region and cytochrome oxidase subunit 1 (COI) to assess phylogenetic relationships and maternal origin of different fish populations [23,24,25]. In this study, we investigated the phenotypic characteristics, chromosomal numbers, rDNA and mtDNA sequences, the DNA methylation level of 5S rDNA IGS sequences between the hybrid grouper and its parents. Based on our data, the genetic traits of the hybrid grouper (C. altivelis (♀) × E. lanceolatus (♂)) were explored.

Materials and methods

Ethics statement

All experiments were conducted in accordance with the guidelines statement of the Administration of Affairs Concerning Animal Experimentation of China. The health of the fish and environmental conditions were monitored daily.

Collection of experimental fish

The fry of E. lanceolatus were derived from Delin Chengxin Aquaculture Co., Ltd, Hainan Province, China. Two-month-old C. altivelis and two-month-old the hybrid groupers were collected from Hainan Chenhai Aquatic Co., Ltd, Hainan Province. All samples were raised up to two years old at Delin Chengxin Aquaculture Co., Ltd. Total genomic DNA was isolated from peripheral blood cells.

Measurement of morphological data

A total of 30 individuals aged 2 years of each from C. altivelis, E. lanceolatus and the hybrid grouper were selected at random. Fish were anaesthetized with 3-aminobenzoic acid ethyl ester methanesulfonate before measuring. The total length, standard length, body height, head length, caudal peduncle length and caudal peduncle height were measured. The average ratios of total length to standard length, head length to standard length, body height to standard length, caudal peduncle length to standard length and caudal peduncle height to standard length were recorded.

Preparation of chromosome spreads

The fin cells of C. altivelis, E. lanceolatus and the hybrid grouper were used for chromosome preparation. The method of fin cell culture and chromosome preparation described by Alvarez et al. was used with some modification [26]. Before collecting cells, concanavalin was added to the cells three times in one-hour intervals. The final concentration was 0.1 μg/mL. The collected cells underwent hypotonic treatment with 0.075 mol/L KCl at 37 °C for 15–20 min and were then fixed in methanol-acetic acid (3:1) for 20 min with two changes. Chromosome preparations were examined under an oil lens. Two-hundred metaphase spreads were examined for each fish sample.

Analysis of the genetic traits of rDNA and mtDNA

Thirty C. altivelis, thirty E. lanceolatus and fifty the hybrid grouper were selected at random. Genomic DNA was extracted from the blood and fin tissue. All PCR primers were listed in Table S1. The PCR amplification program was conducted for 30 cycles with an annealing temperature of 30 s at 57 °C for 5S rDNA, 35 s at 55 °C for COI gene and 30 s at 50 °C for D-loop region. All PCR products underwent a series of applications such as purification, cloning, and sequencing. Multiple sequence alignment analysis was conducted using BioEdit and Clustal W [27, 28]. The RepeatMasker program (http://www.repeatmasker.org/) was used to search for 5S gene repeated elements inside spacers. All sequences were submitted to the NCBI.

Methylation analysis of 5S rDNA IGS sequence

The DNA methylation level in 1021 bp upstream sequence of coding region was tested and analyzed to investigate the DNA methylation patterns of IGS regions in the hybrid grouper and its parents (Supplementary Figure 1). Muscle tissues were applied to extract genomic DNA in C. altivelis (five individuals), E. lanceolatus (seven individuals) and the hybrid grouper (five individuals). MethPrimer software package (http://www.urogene.org/methprimer/index.html) was applied for designing PCR primers (Table S1). PCR amplification was conducted with an annealing temperature of 30 s at 50 ~ 60 °C by KAPA HiFi HotStart Uracil + ReadyMix PCR Kit (Kapa Biosystems, Wilmington, MA, USA). Steps of bisulfite modification and pyrosequencing were showed as follow: firstly, genomic DNA (1 μg) was converted (by the ZYMO EZ DNA Methylation-Gold Kit (Zymo Research, Irvine, CA, USA)) and one-tenth of the elution products were used as templates. Bisulfite sequencing PCR products of IGS sequences were pooled equally, 5'-phosphorylated, 3' -dA-tailed and ligated to barcoded adapter using T4 DNA ligase (NEB). Secondly, Barcoded libraries from all samples were sequenced on Illumina platform [29]. Thirdly, the raw reads needed to remove adapters and filter the low-quality sequences (Trimmomatic-0.36). Finally, clear reads were obtained. The methylation levels of individual cytosines were calculated as the ratio of the total number of methylated CpG cytosines to the number of sequenced clones. A two-tailed Fisher’s Exact Test was used to test the methylated and unmethylated counts for each cytosine between two groups [30].

Phylogenetic analysis and genetic distances

Mitochondrial genome sequences were obtained from the NCBI GenBank. All accession numbers of the sequences were listed in Table S2. Phylogenetic trees were constructed by a maximum-likelihood approach (ML) in the program PhyML 3.0 (http://www.atgc-montpellier.fr/phyml/) and Bayesian inference (BI) using MrBayes 3.1.2 (http://mrbayes.sourceforge.net/). A generalized time reversible (GTR) model with a proportion of invariable sites (I) and a gamma distribution (G) was selected for the concatenated dataset in ML. Nodal robustness was estimated by bootstrap percentage values (BP). The BP values of nodes of the subtree were obtained after 1000 replicates. MrBayes used Markov chain Monte Carlo (MCMC) methods to estimate the posterior distribution of model parameters. Convergence occurred when average standard deviation of split frequencies fell below 0.01. Posterior probability (PP) values were applied to estimate branch support. The phylogenetic trees were shown using the Figtree program v.1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/).

Between-group mean distances were calculated using MEGA 7 with a p-distance model.

Results

Morphological traits

There were obvious phenotypic differences between the hybrid grouper and its parents (C. altivelis and E. lanceolatus) (Fig. 1). Regarding body color, C. altivelis is white with scattered black dots and E. lanceolatus is dark brown with scattered black spots. The hybrid grouper is dark brown without any obvious black spots, which differed visibly from both parents.

Fig. 1
figure 1

The appearances of Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper

The measurable traits in C. altivelis, E. lanceolatus, and the hybrid grouper are presented in Table 1. The values of total length, standard length, head length, body height, and caudal peduncle height for E. lanceolatus were greater than in C. altivelis. However, the caudal peduncle length was greater for C. altivelis than for E. lanceolatus. All measurable traits differed significantly between E. lanceolatus and C. altivelis (p < 0.05). For the hybrid grouper, all measurable traits had values higher than those of C. altivelis and lower than those of E. lanceolatus (excepting for caudal peduncle length). All measurable traits with the exception of head length and body height showed significant differences between the hybrid and C. altivelis (p < 0.05).

Table 1 Comparison of the measurable traits between Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper

To control for differences in growth rate between species, the ratios of the measurable traits were recorded in the hybrid grouper and its parents. The ratios of measurable traits did not differ significantly between the hybrid grouper, C. altivelis, and E. lanceolatus (p > 0.05), which indicates that all three groupers show similar morphological characteristics in terms of total length, head length, standard length, caudal peduncle length, body height, and caudal peduncle height.

Examination of chromosome number

The distribution of chromosome numbers in C. altivelis, E. lanceolatus, and the hybrid grouper is illustrated in Fig. 2. Among C. altivelis, 92% of the chromosomal metaphases had 48 chromosomes. Among E. lanceolatus, 90% of the chromosomal metaphases had 48 chromosomes. Among the hybrid grouper, 89% of the chromosomal metaphases had 48 chromosomes.

Fig. 2
figure 2

Chromosome spreads at metaphase in Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper. a. The 48 chromosomes of Cromileptes altivelis. b. The 48 chromosomes of Epinephelus lanceolatus. c. The 48 chromosomes of the hybrid grouper

Sequence organization and analysis of rDNA and mtDNA

A single band of 5S rDNA was observed in the genomes of C. altivelis (403 bp, accession numbers: OM289959), E. lanceolatus (402 bp, accession numbers: OM289957) and the hybrid grouper. Based on the nucleotide compositions of NTS, 5S sequences in the hybrid grouper were divided into two categories (type I: 402 bp, accession numbers: OM289958; type II: 404 bp, accession numbers: OP244353). All 5S clones in samples had the same coding region sequences (Fig. 3a). The internal control regions (ICRs, i.e., the promoters for transcription)-Box A (positions 48–62), Box B (positions 65–70), and Box C (positions 78–95)-were detected in the 5S coding region. In 5S NTS region, several base substitutions or insertions/deletions were found in the hybrids and their parents. Sequence alignments of the NTS region with BLASTn showed 96.2% similarity between the hybrid grouper (type I) and C. altivelis, and 99.0% similarity between the hybrid grouper (type I) and E. lanceolatus. Type II of the hybrid grouper displayed higher conservation with C. altivelis (99.2%) compared with E. lanceolatus (95.0%) (Table 2). The TATA control element that regulates 5S gene transcription was generally located from -25 to -30 from the 5S, was also identifiable in the hybrid and parents (at -29 in all NTS sequences, where it was modified to TAAA) (Fig. 3b). In addition, long terminal repeat (LTR) retrotransposons were identified by the RepeatMasker program in the NTS regions of C. altivelis (pos. 18–122) and the hybrid grouper (type I: pos. 51–122; type II: pos. 117–223), and the (T)n simple repeat was identified in that of E. lanceolatus (pos. 88–122) (Fig. 4).

Fig. 3
figure 3

Representative sequences of 5S rDNA. a. Complete 5S coding regions from Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper. Internal control regions of the coding region are shaded. b. Comparison of the NTS sequences from Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper, the NTS upstream TATA elements are shaded; asterisks mark variable sites in the NTS

Table 2 The similarity and genetic distance between Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper in 5S rDNA
Fig. 4
figure 4

Schematic representation of the 5S rDNA sequences showing the different types of elements found in the NTS of Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper

The length of COI gene sequences was 1545 bp in C. altivelis, E. lanceolatus, and the hybrid grouper (Supplementary Figure 2). Sequence alignments with BLASTn showed that the COI gene sequences of C. altivelis and E. lanceolatus in this study was identical with the sequences in NCBI database (Cromileptes altivelis: NC_021614.1, Epinephelus lanceolatus: KM386619.1). The similarity between the hybrid grouper and C. altivelis was 99.8% (Table 3). Partial D-loop region sequences (-800 bp) were obtained in C. altivelis, E. lanceolatus, and the hybrid grouper (Supplementary Figure 3). The similarity between the hybrid grouper and C. altivelis was 99.0% (Table 3).

Table 3 The similarity and genetic distance between Cromileptes altivelis, Epinephelus lanceolatus and the hybrid grouper in COI gene and D-loop region

DNA methylation analysis of 5S rDNA IGS sequence

Seven CpG sites were found in the 1021 bp upstream sequence of coding region between the hybrid grouper and its parents (Fig. 5). Only two CpG sites (the position of 107 and 178 in the sequence) exhibited the significant differences between the hybrid grouper and its parents (Fig. 5A, B). In site 107, the DNA methylation level of the hybrid grouper (9.66%) was higher than that of E. lanceolatus (7.59%). There was no significant difference between the hybrid grouper and C. altivelis (Fig. 5C). In site 178, the DNA methylation level of the hybrid grouper (14.60%) was higher than that of E. lanceolatus (6.31%), but lower than that of C. altivelis (27.06%) (Fig. 5D).

Fig. 5
figure 5

DNA methylation of 5S rDNA intergenic spacers (IGS) sequence in Cromileptes altivelis, Epinephelus lanceolatus, the hybrid grouper. a. Schematic representation of the DNA methylation level in CpG sites. The box indicates methylation level of CpG sites with significant difference between the hybrid grouper and parents. b. The partial sequence of IGS. The red lines indicate CpG sites. c. The DNA methylation level in site 107. d. The DNA methylation level in site 178. The symbol * indicates significant differences

Genetic distances and phylogenetic analysis

Genetic distances were estimated between C. altivelis, E. lanceolatus and the hybrid grouper in mtDNA and rDNA (Tables 2 and 3). The hybrids and C. altivelis had the minimal genetic distance (COI gene: 0.002; D-loop region: 0.010). In 5S rDNA, the genetic distance between the hybrids and C. altivelis ranged from 0.007 to 0.050 (average value = 0.0285) and was less than the distance between the hybrids and E. lanceolatus, which ranged from 0.014 to 0.057 (average value = 0.0355).

Phylogenetic analyses generated similar tree topologies in ML (Fig. 6) and BI (Fig. 7). In the phylogenetic tree, the hybrid grouper and C. altivelis were closely clustered together (COI gene: BP = 100%, PP = 1; D-loop region: BP = 100%, PP = 1). E. lanceolatus and other Epinephelus species were clustered together.

Fig. 6
figure 6

Maximum likelihood phylogenetic tree reconstructed based on COI gene (a) and D-loop region (b)

Fig. 7
figure 7

Bayesian phylogenetic tree reconstructed based on COI gene (a) and D-loop region (b)

Discussion

Geometric morphometrics is considered to be one of the most powerful tools in fish shape analysis [31, 32], and can most intuitively express germplasm characteristics. In our results, obvious differences in the body color and blot shape were visible between the hybrid offspring and its parents (Cromileptes altivelis and Epinephelus lanceolatus) at 2 years of age. The hybrid grouper displayed intermediate morphology compared with the parents’ measurable characteristics. In total length, standard length, caudal peduncle length, and caudal peduncle height, there were significant differences between the offspring and C. altivelis. These observations indicated that the offspring's phenotypic traits exhibited the hybridity, and supported the viewpoint that the hybrid offspring exhibits faster growth than its female parent (C. altivelis) [33].

The chromosome number determines the characteristics of the species, and alterations in the chromosome number can reflect changes in genetic material during species evolution. Previous studies have shown that the ploidy level of C. altivelis and E. lanceolatus is diploid (2n = 48, [34, 35]). The ploidy level of the hybrid grouper has not been reported yet. In the present study, the hybrid grouper (2n = 48) was confirmed to be diploid by chromosome count. Such diploid offspring have previously been observed in other crosses between members of Epinephelinae: (Epinephelus fuscoguttatus, ♀ × Epinephelus lanceolatus, ♂), (Epinephelus moara, ♀ × Epinephelus lanceolatus, ♂), and Epinephelus awoara, ♀ × Epinephelus tukula, ♂) [35, 36].

To further analyze genetic characteristics, rDNA and mtDNA sequences were compared between the hybrid grouper and its parents. In rDNA coding region sequences, a low intra-specific and a high rate of homogenization were observed between the hybrid grouper, the female parent (C. altivelis), and the male parent (E. lanceolatus), in which the sequence similarity of 5S was 100%, indicating the presence of concerted evolution in the rDNA of these species. In the concerted evolution model, all the members of a multi-gene family evolve in concert. Variation in the repeat unit extends to all the member genes through gene conversion and unequal crossing over, which leads to homogenization in the units of the multigene family [37, 38].

Different 5S rDNA types have been reported in freshwater fish and several plants [39,40,41,42]. In the bony fish, there are often two 5S gene fragments [43], such as in Carassius auratus (2n = 100) and Cyprinus carpio (2n = 100). In our study, only one 5S gene band was found in the hybrid (2n = 48), C. altivelis (2n = 48), and E. lanceolatus (2n = 48). This result was previously observed in the genome of Megalobrama amblycephala (2n = 48). Researchers [38] have proposed that the genomes of diploid Carassius auratus and Cyprinus carpio with 100 chromosomes likely experienced a polyploidization process and are ancient polyploidy fish. During the process, the variations and reorganizations of the genomes may have resulted in the appearance of new 5S rDNA. There was no evidence to suggest that C. altivelis and E. lanceolatus with 48 chromosomes went through the polyploidization process, and it seems reasonable to observe only one 5S gene type in their genome. In this study, all 5S gene sequences contained essential internal control regions (Box A, Box B, Box C, and TATA control element) for correct gene expression; thus, these 5S genes were likely to correspond to functional genes [42, 44].

In the 5S NTS region, nucleotide variations (including base substitutions or insertions/deletions) discovered in the hybrids, C. altivelis and E. lanceolatus demonstrate the divergence exerted in this region, as predicated by a birth-and-death evolution model [45]. Birth-and-death evolution leads to heterogeneity and pseudogenes in an rDNA multigene family [46, 47]. A mixed-mediated model between concerted evolution and birth-and-death has been described for the fish species Diplodus sargus [48], and Halobatrachus didactylus [49]. LTR retrotransposons were identified in the NTS regions of C. altivelis and the hybrid grouper. The presence of LTR retrotransposons has been associated with genome expansion because it can potentially give rise to numerous daughter copies [50].

Phylogenetic analysis revealed that all fish species were well clustered according to their taxonomic level. Matrilineal inheritance of mtDNA was observed in C. altivelis × E. lanceolatus and other hybrids (E. fuscoguttatus × E. lanceolatus, E. moara × E. lanceolatus and E. coioides × E. lanceolatus). All hybrids were closely clustered with their female parent. The genetic distance between the hybrid grouper and female parent (C. altivelis) was lower than those between the hybrid grouper and male parent (E. lanceolatus) in 5S rDNA and mtDNA. These results indicate that the hybrid grouper had a closer genetic relationship with female parent. The genetic distance of 5S rDNA between type I of the hybrids and type II of the hybrids was 0.043. Type I and E. lanceolatus had the minimum genetic distance (0.014), type II and C. altivelis had the minimum genetic distance (0.007). These observations have suggested the presence of genetic material from both parents in the hybrid grouper.

The synthesis of rRNA genes is closely related to the complex process of ribosome biogenesis, and can be regulated by changes of the transcription rate and the number of gene copies that are transcribed [51]. DNA methylation of rRNA gene promoter is an epigenetic switch, which can regulate rRNA transcription by controlling the number of rRNA genes in the on or off state [52]. It is widely known that hypomethylation of rRNA genes promoter can facilitates the binding of transcription factor on chromatin templates and drive the synthesis of rRNA in order for ribosomes biogenesis [51, 53]. In our results, methylation level of sites 107 and 178 presented the significant difference between the hybrid grouper and its parents. Among them, the methylation level of the hybrid grouper was in the middle compared with that of parents. Hybridization often results in dramatic genome reconfigurations including epigenetic changes [54]. Scholars have proposed that genome epigenetic regulation is an important factor for the formation of heterosis [55, 56]. In this study, results of geometric morphometrics indicated that the growth rate of 2-year-old the hybrid grouper was higher than that of C. altivelis, exhibiting the growth advantage. As it is known, grouper generally reaches sexual maturity at 3—5 years. Before reaching sexual maturity, the demand for energy is mainly used for growth and development. The protein synthesis is indispensable in the process of growth and development. Therefore, we infer that DNA methylation of 5S rDNA IGS sequences may indirectly affect the growth and development of grouper by regulating the rate of protein synthesis, but it still requires extensive evidences to support this viewpoint.

Conclusions

The work in this study firstly confirmed the existence of diploid the hybrid grouper with 48 chromosomes. The molecular and phylogenetic analysis of rDNA and mtDNA recovered that the hybrid grouper has inherited the genetic materials from both of parents. The morphological trait data showed the hybrid grouper exhibits improved growth compared with C. altivelis. The DNA methylation of 5S rDNA IGS sequences in the hybrid grouper was also lower than that of C. altivelis, which provided a clue to understanding the growth advantage in the hybrid grouper.

Availability of data and materials

All raw data of bisulfite sequencing in Cromileptes altivelis × Epinephelus lanceolatus and parents have been uploaded to GenBank database (BioProject ID: PRJNA1009716, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1009716/). All 5S rDNA sequences generated during the current study are available in the NCBI (accession number: OM289957-OM289959, OP244353).

References

  1. Pierre S, Gaillard S, Prévot-D’Alvise N, Aubert J, Rostaing-Capaillon O, Leung-Tack D, Grillasca JP. Grouper aquaculture: Asian success and Mediterranean trials. Aquat Conserv Mar Freshwat Ecosyst. 2008;18:297–308.

    Article  Google Scholar 

  2. Guo CY, Huang YH, Wei SN, Ouyang ZL, Yan Y, Huang XH, Qin QW. Establishment of a new cell line from the heart of giant grouper, Epinephelus lanceolatus (Bloch), and its application in toxicology and virus susceptibility. J Fish Dis. 2013;38:175–86.

    Article  PubMed  Google Scholar 

  3. Ding SX, Wang YH, Wang J, Zhuang X, Su YQ, You YZ, Li QF. Molecular phylogenetic relationships of 30 grouper species in China Seas based on 16S rDNA fragment sequences. Acta Zoologica Sinica. 2006;52:504–13 (in Chinese).

    CAS  Google Scholar 

  4. Gao Y, Lu S, Wu M, Wei Y, Jin Z, Wu X. Effects of dietary protein levels on growth, feed utilization and expression of growth related genes of juvenile giant grouper (Epinephelus lanceolatus). Aquaculture. 2019;504:369–74.

    Article  CAS  Google Scholar 

  5. Myoung JG, Kang CB, Yoo JM, Lee EK, Kim S, Jeong CH, Kim BI. First Record of the Giant Grouper Epinephelus lanceolatus (Perciformes: Serranidae: Epinephelinae) from Jeju Island, South Korea. Fish Aquat Sci. 2013;16:49–52.

    Google Scholar 

  6. Ou YJ, Li JE, Chen FH. Introduction, Acclimation and Induction of Gonad Development and Reproduction of High Finned Grouper, Cromileptes altivelis. J Zhanjiang Ocean Univ. 1999;3:20–32 (in Chinese).

    Google Scholar 

  7. Urbanelli S, Porretta D, Mastrantonio V, Bellini R, Pieraccini G, Romoli R, Crasta G, Nascetti G. Hybridization, natural selection, and evolution of reproductive isolation: a 25-years survey of an artificial sympatric area between two mosquito sibling species of the Aedes mariae complex. Evolution. 2014;68:3030–8.

    Article  PubMed  Google Scholar 

  8. Liu Y, Tian YS, Wang LN, Li ZT, Zhang JJ, Wu YP, Chen S, Li LL, Li WS, Ma WH, Lin HW, Wang QB. Construction of high-density linkage maps and QTL mapping for growth-related traits in F1 hybrid Yunlong grouper (Epinephelus moara♀ × E. lanceolatus♂). Aquaculture. 2021; 548:737698.

  9. Li ZT, Tian YS, Cheng Ml, Wang L, Zhang JJ, Wu YP, et al. The complete mitochondrial genome of the hybrid grouper Epinephelus moara (♀)×Epinephelus tukula (♂), and phylogenetic analysis in subfamily Epinephelinae. Acta Oceanol Sin. 2020;39, 65–75.

  10. Wu S, Li J, Ou Y, Lv G, Liu J. Allometric growth of hybrid grouper (Epinephelus coioides♀×E. lanceolatus♂) larvae and juveniles. J Fish Sci China. 2014;21: 503–510. (in Chinese).

  11. Strauss RE. Evolutionary Allometry and Variation in Body Form in the South American Catfish Genus Corydoras (Callichthyidae). Syst Biol. 1985;34:381–96.

    Article  Google Scholar 

  12. Cao L, Qin Q, Xiao Q, Yin H, Wen J, Liu Q, et al. Nucleolar dominance in a tetraploidy hybrid lineage derived from Carassius auratus red var. (♀) × Megalobrama amblycephala (♂). Front Genet. 2018;9:386.

  13. Cao L, Zhao C, Wang C, Qin H, Liu S. Evolutionary dynamics of 18S and 5S rDNA in autotriploid Carassius auratus. Gene. 2020;737: 144433.

    Article  CAS  PubMed  Google Scholar 

  14. Gong L, Shi W, Yang M, Luo H. Variations in the conserved 18S and 5.8S reveal the putative pseudogenes in 18S-ITS1–5.8S rDNA of Cynoglossus melampetalus (Pleuronectiformes: Cynoglossidae). Biochem Biophys Res Commun. 2021;534:233–239.

  15. Ganley A, Kobayashi T. Highly efficient concerted evolution in the ribosomal DNA repeats: total rDNA repeat variation revealed by whole-genome shotgun. Genome Res. 2007;17:184–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Mayer C, Schmitz K, Grummt I. Santoro, R. Intergenic transcripts regulate the epigenetic state of rRNA genes. Mol. Cell. 2006;22:351–361.

  17. Bird AP. CpG-rich islands and the function of DNA methylation. Nature. 1986;321:209–13.

    Article  CAS  PubMed  Google Scholar 

  18. Ou M, Mao HH, Luo Q, Zhao J, Liu HY, Zhu XP, et al. The DNA methylation level is associated with the superior growth of the hybrid fry in snakehead fish (Channa argus × Channa maculata). Gene. 2019;703:125–33.

    Article  CAS  PubMed  Google Scholar 

  19. Jiang CD, Deng CY, Xiong YZ. Differences of cytosine methylation in parental lines and F1 hybrids of Large White × Meishan crosses and their effects on F1 performance. J Agric Biotechnol. 2007;2:195–200.

    Google Scholar 

  20. Yang H, Li Q. The DNA methylation level is associated with the superior growth of the hybrid crosses in the Pacific oyster Crassostrea gigas. Aquaculture. 2021;547:737421.

  21. Harpending HC, Batzer MA, Gurven M, Jorde LB, Rogers AR, Sherry ST. Genetic traces of ancient demography. Proc Natl Acad Sci USA. 1998;95:1961–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Manchado M, Catanese G, Ponce M, Funes V, Infante C. The complete mitochondrial genome of the Senegal sole, Solea senegalensis Kaup. comparative analysis of tandem repeats in the control region among soles. DNA Seq. 2007;18:169–175.

  23. Yang B, Chen XY, Yang JX. Structure of the Mitochondrial DNA Control Region and Population Genetic Diversity Analysis of Anabarilius grahami check for this species in other resources (Regan). Zool Res. 2008;29:379–85.

    Article  CAS  Google Scholar 

  24. Homola JJ, Kanefsky J, Scribner KT, Kalish TG, Tonello MA. Genetic identification of two putative world record Michigan salmonids resolves stakeholder and manager questions. J Great Lakes Res. 2012;38:176–9.

    Article  CAS  Google Scholar 

  25. Mu XD, Wang XJ, Liu Y, Song HM, Liu C, Gu DG, Wei H, Luo JR, Hu YC. An unusual mitochondrial genome structure of the tonguefish, Cynoglossus trigrammus: Control region translocation and a long additional non-coding region inversion. Gene. 2015;573:216–24.

    Article  CAS  PubMed  Google Scholar 

  26. Alvarez MC, Otis J, Amores A, Guise K. Short-term cell culture technique for obtaining chromosomes in marine and freshwater fish. J Fish Biol. 2010;39:817–24.

    Article  Google Scholar 

  27. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.

  28. Hall TA. BioEdit: A User-Friendly Biological Sequence Alignment Editor and Analysis Program for Windows 95/98/NT. Nuclc Acids Symposium Series. 1999;41:95–8.

    CAS  Google Scholar 

  29. Gao F, Liang H, Lu H, Wang J, Xia M, Yuan Z, Yao Y, Wang T, Tan X, Laurence A, et al. Global analysis of DNA methylation in hepatocellular carcinoma by a liquid hybridization capture-based bisulfite sequencing approach. Clin Epigenetics. 2015;7:86.

  30. Gao F, Zhang J, Jiang P, Gong D, Wang JW, Xia YD, et al. Marked methylation changes in intestinal genes during the perinatal period of preterm neonates. BMC Genomics. 2014;15:716.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Bookstein FL. Combining the tools of geometric morphometrics. In: Marcus LF, Corti M, Loy A, Naylor GJP, Slice DE. (eds) Advances in morphometrics. NATO ASI Series. 1996.

  32. Bookstein FL. Shape and the Information in Medical Images: A Decade of the Morphometric Synthesis. Comput Vis Image Underst. 1997;66:97–118.

    Article  Google Scholar 

  33. Chen J, Ye Z, Yu Z, Wang J, Li P, Chen X, et al. The complete mitochondrial genome of the hybrid grouper (Cromileptes altivelis♀ × Epinephelus lanceolatus♂) with phylogenetic consideration. Mitochondrial DNA Part B. 2017;2:171–2.

  34. Ou YJ, Xie J. Karyotype analysis of Cromileptes altivelis. South China Fish Sci. 2007;3:49–53 (in Chinese).

    Google Scholar 

  35. Liu L, Zhang Y, Chen C, Li Y, Kong X, Yu H, et al. Preparation of Chromosomes and Analysis of Karyotypes of Hybrid F_1(Epinephelus fuscoguttatus♀×Epinephelus lanceolatus♂)Groupers. Prog Fish Sci. 2016;37:33–40 (in Chinese).

    Google Scholar 

  36. Cheng M, Tian Y, Wu Y, Li Z, Zhang J, Wang L, et al. Chromosome Karyotype Analysis of Hybrids of Epinephelus moara♀×Epinephelus lanceolatus♂. Prog Fish Sci. 2018;39:67–75 (in Chinese).

    Google Scholar 

  37. Nei M, Rooney AP. Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005;39:121.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Qin Q, He W, Liu S, Wang J, Liu Y. Analysis of 5S rDNA organization and variation in polyploid hybrids from crosses of different fish subfamilies. J Exp Zool B Mol Dev Evol. 2010;314B:403–11.

    Article  CAS  Google Scholar 

  39. Martins C, Galetti PM Jr. Organization of 5S rDNA in species of the fish Leporinus: two different genomic locations are characterized by distinct nontranscribed spacers. Genome. 2001;44:903–10.

    Article  CAS  PubMed  Google Scholar 

  40. Wasko AP, Martins C, Wright JM, Galetti PM Jr. Molecular organization of 5S rDNA in fishes of the genus <i>Brycon</i>. Genome. 2001;44:893–902.

    Article  CAS  PubMed  Google Scholar 

  41. Negi MS, Rajagopal JY, Chauhan N, Cronn R, Lakshmikumaran M. Length and sequence heterogeneity in 5S rDNA of Populus deltoides. Genome. 2002;45:1181–8.

    Article  CAS  PubMed  Google Scholar 

  42. Volkov RA, Panchuk II, Borisjuk NV, Hosiawa-Baranska M, Maluszynska J, Hemleben V. Evolutional dynamics of 45S and 5S ribosomal DNA in ancient allohexaploid Atropa belladonna. BMC Plant Biol. 2017;17:21.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Martins C, Galeti PM Jr. Two 5S rDNA arrays in neotropical fish species: is it a general rule for fishes? Genetica. 2001;111:439–46.

    Article  CAS  PubMed  Google Scholar 

  44. Cloix C, Tutois S, Mathieu O, Cuvillier C, Espagnol MC, Picard G, et al. Analysis of 5S rDNA Arrays in Arabidopsis thaliana: Physical Mapping and Chromosome-Specific Polymorphisms. Genome Res. 2000;10:679.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Rebordinos L, Cross I, Merlo A. High Evolutionary Dynamism in 5S rDNA of Fish: State of the Art. Cytogenet Genome Res. 2013;141:103–13.

    Article  CAS  PubMed  Google Scholar 

  46. Bardella VB, Cabral-De-Mello DC. Uncovering the molecular organization of unusual highly scattered 5S rDNA: The case of Chariesterus armatus (Heteroptera). Gene. 2018;646:153–8.

    Article  CAS  PubMed  Google Scholar 

  47. Yang M, Kong X, Wei S, Gong L, Luo H, Wu B. Remarkable sequence polymorphisms in 18S rDNA of Pleuronichthys cornutus (Pleuronectiformes Pleuronectidae). Gene. 2018;677:251–8.

    Article  CAS  PubMed  Google Scholar 

  48. Merlo MA, Cross I, Manchado M, Cárdenas S, Rebordinos L. The 5S rDNA High Dynamism in Diplodus sargus is a Transposon-Mediated Mechanism. Comparison with Other Multigene Families and Sparidae Species. J Mol Evol. 2013;76:83–97.

  49. Merlo MA, Pacchiarini T, Portela-Bens S, Cross I, Manchado M, Rebordinos L. Genetic characterization of Plectorhinchus mediterraneus yields important clues about genome organization and evolution of multigene families. BMC Genet. 2012;13:33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Ungerer MC, Strakosh SC, Stimpson KM. Proliferation of Ty3/gypsy-like retrotransposons in hybrid sunflower taxa inferred from phylogenetic data. BMC Biol. 2009;7:1–13.

    Article  Google Scholar 

  51. Grummt I. Different epigenetic layers engage in complex crosstalk to define the epigenetic state of mammalian rRNA genes. Hum Mol Genet. 2007;16(Spec No 1):R21–27.

  52. Lawrence RJ, Pikaard CS. Chromatin turn ons and turn offs of ribosomal RNA genes. Cell Cycle. 2004;3:880–3.

    Article  CAS  PubMed  Google Scholar 

  53. Deaton AM, Bird A. CpG islands and the regulation of transcription. Genes Dev. 2011;25:1010–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Koroma AP, Jones R, Michalak P. Snapshot of DNA methylation changes associated with hybridization in Xenopus. Physiol Genomics. 2011;43:1276–80.

    Article  CAS  PubMed  Google Scholar 

  55. Romagnoli S, Maddaloni M, Livini C, Motto M. Relationship between gene expression and hybrid vigor in primary root tips of young maize (Zea mays L.) plantlets. Theor Appl Genet. 1990;80(6):769–775.

  56. Tsaftaris AS, Kafka M. Mechanisms of Heterosis in Crop Plants. J Crop Prod. 1997;1(1):95–111.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to sincerely thank many researchers who help to complete this manuscript.

Funding

The research was supported by the specific research fund of National Natural Science Foundation of China (32160861), Hainan Provincial Natural Science Foundation of China (321QN263), The Major Science and Technology plan of Hainan Province (ZDKJ2021017), the Youth Project of Yazhou Bay Innovation Institute of Hainan Tropical Ocean University (2022CXYQNXM06), State Key Laboratory of Developmental Biology of Freshwater Fish (2020KF001), National Natural Science Foundation of China (32002389) and Scientific Research Foundation of Hainan Tropical Ocean University (RHDRC202010).

Author information

Authors and Affiliations

Authors

Contributions

H.H contributed to the conception and designed the study, L.C. carried out the experimental work, participated in rafted the manuscript. Y.L. analyzed sequences and P.C., X.H., J.M. and N.Y. participated in interpretation and discussion of the results. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hai Huang.

Ethics declarations

Ethics approval and consent to participate

This study was approved by Animal Ethical and Welfare Committee of Hainan Tropical Ocean University (HTOU) and followed the guidelines statement of the Administration of Affairs Concerning Animal Experimentation of China. This study does not involve the use of any human data or tissue. The animals used in the study were obtained from Delin Chengxin Aquaculture Co., Ltd, and Hainan Chenhai Aquatic Co., Ltd from which we have obtained consent to use these animals in our research.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. The primers used for PCR and BSP-PCR. Table S2. GenBank accession numbers for mitochondrial genome sequences.

Additional file 2:

Supplementary Figure 1. Representative sequences of 5S rDNA intergenic spacers (IGS) sequence.

Additional file 3:

Supplementary Figure 2. Representative sequences of COI gene from Cromileptes altivelis, Epinephelus lanceolatus, the hybrid grouper.

Additional file 4:

Supplementary Figure 3. Representative partial sequences of D-loop region from Cromileptes altivelis, Epinephelus lanceolatus, the hybrid grouper.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cao, L., Chen, P., Hou, X. et al. rDNA and mtDNA analysis for the identification of genetic characters in the hybrid grouper derived from hybridization of Cromileptes altivelis (female) × Epinephelus lanceolatus (male). BMC Genom Data 25, 5 (2024). https://doi.org/10.1186/s12863-023-01188-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12863-023-01188-5

Keywords