Chromosomal mapping of microsatellite repeats in the rock bream fish Oplegnathus fasciatus, with emphasis of their distribution in the neo-Y chromosome

Genome-wide sequence divergence of satellite DNA could underlie meiotic failure in male hybrids of bighead catfish and North African catfish (Clarias, Clariidae)

Hybrid sterility, a hallmark of postzygotic isolation, arises from parental genome divergence disrupting meiosis. While chromosomal incompatibility is often implicated, the underlying mechanisms remain unclear. This study investigated meiotic behavior and genome-wide divergence in bighead catfish (C. macrocephalus), North African catfish (C. gariepinus), and their sterile male hybrids (important in aquaculture). Repetitive DNA analysis using bioinformatics and cytogenetics revealed significant divergence in satellite DNA (satDNA) families between parental species. Notably, one hybrid exhibited successful meiosis and spermatozoa production, suggesting potential variation in sterility expression. Our findings suggest that genome-wide satDNA divergence, rather than chromosome number differences, likely contributes to meiotic failure and male sterility in these catfish hybrids.

Innovative approach for high-throughput exploiting sex-specific markers in Japanese parrotfish Oplegnathus fasciatus

BACKGROUND

The use of sex-specific molecular markers has become a prominent method in enhancing fish production and economic value, as well as providing a foundation for understanding the complex molecular mechanisms involved in fish sex determination. Over the past decades, research on male and female sex identification has predominantly employed molecular biology methodologies such as restriction fragment length polymorphism, random amplification of polymorphic DNA, simple sequence repeat, and amplified fragment length polymorphism. The emergence of high-throughput sequencing technologies, particularly Illumina, has led to the utilization of single nucleotide polymorphism and insertion/deletion variants as significant molecular markers for investigating sex identification in fish. The advancement of sex-controlled breeding encounters numerous challenges, including the inefficiency of current methods, intricate experimental protocols, high costs of development, elevated rates of false positives, marker instability, and cumbersome field-testing procedures. Nevertheless, the emergence and swift progress of PacBio high-throughput sequencing technology, characterized by its long-read output capabilities, offers novel opportunities to overcome these obstacles.

FINDINGS

Utilizing male/female assembled genome information in conjunction with short-read sequencing data survey and long-read PacBio sequencing data, a catalog of large-segment (>100 bp) insertion/deletion genetic variants was generated through a genome-wide variant site-scanning approach with bidirectional comparisons. The sequence tagging sites were ranked based on the long-read depth of the insertion/deletion site, with markers exhibiting lower long-read depth being considered more effective for large-segment deletion variants. Subsequently, a catalog of bulk primers and simulated PCR for the male/female variant loci was developed, incorporating primer design for the target region and electronic PCR (e-PCR) technology. The Japanese parrotfish (Oplegnathus fasciatus), belonging to the Oplegnathidae family within the Centrarchiformes order, holds significant economic value as a rocky reef fish indigenous to East Asia. The criteria for rapid identification of male and female differences in Japanese parrotfish were established through agarose gel electrophoresis, which revealed 2 amplified bands for males and 1 amplified band for females. A high-throughput identification catalog of sex-specific markers was then constructed using this method, resulting in the identification of 3,639 (2,786 INS/853 DEL, ♀ as reference) and 3,672 (2,876 INS/833 DEL, ♂ as reference) markers in conjunction with 1,021 and 894 high-quality genetic sex identification markers, respectively. Sixteen differential loci were randomly chosen from the catalog for validation, with 11 of them meeting the criteria for male/female distinctions. The implementation of cost-effective and efficient technological processes would facilitate the rapid advancement of genetic breeding through expediting the high-throughput development of sex genetic markers for various species.

CONCLUSIONS

Our study utilized assembled genome information from male and female individuals obtained from PacBio, in addition to data from short-read sequencing data survey and long-read PacBio sequencing data. We extensively employed genome-wide variant site scanning and identification, high-throughput primer design of target regions, and e-PCR batch amplification, along with statistical analysis and ranking of the long-read depth of the variant sites. Through this integrated approach, we successfully compiled a catalog of large insertion/deletion sites (>100 bp) in both male and female Japanese parrotfish.

A comparative cytogenetic study of Hypsibarbusmalcolmi and H.wetmorei (Cyprinidae, Poropuntiini)

Cyprininae are a highly diversified but demonstrably monophyletic lineage of cypriniform fishes. Here, the karyotype and chromosomal characteristics of Hypsibarbusmalcolmi (Smith, 1945) and H.wetmorei (Smith, 1931) were examined using conventional, nucleolus organizing regions (NORs) and molecular cytogenetic protocols. The diploid chromosome number (2n) of H.malcolmi was 50, the fundamental number (FN) was equal to 62, and the karyotype displayed 8m + 4sm + 38a with NORs located at the centromeric and telomeric positions of the short arms of chromosome pairs 1 and 2, respectively. 2n of H.wetmorei was 50, FN 78, karyotype 14m + 14sm + 22a with the NORs at the telomeric position of the short arm of chromosome pair 2. 2n and FN in males and females were identical. Fluorescence in situ hybridization using different microsatellite motifs as probes also showed substantial genomic divergence between both studied species. In H.wetmorei, (CAG)n and (CAC)n microsatellites accumulated in the telomeric regions of all chromosomes, while in H.malcolmi, they had scattered signals on all chromosomes. Besides, the (GAA)n microsatellites were distributed along all chromosomes of H.malcolmi, but there was a strong hybridization pattern in the centromeric region of a single pair in H.wetmorei. These cytogenomic difference across the genomes of these Hypsibarbus Rainboth, 1996 species are markers for specific evolutionary differentiation within these two species.

Characterization of Early Gonadal Differentiation and Induction of Sex Reversal in the Rock Bream Oplegnathus fasciatus

Rock bream (Oplegnathus fasciatus) is a typical fish that has a unique multiple sex chromosome system (♀X1X1X2X2/♂X1X2Y). We examined the early gonadal development in rock bream via continuous histological observations of the gonads at 40-120 days post hatching (dph). The fish was identified as a typical gonochorist, and female gonads were found to differentiate earlier than male gonads. The ovarian cavity of the female was initially observed at 80 dph, whereas the efferent duct of the male was not observed until 100 dph. Immunofluorescence with the vasa-antibody revealed that germ cells were predominantly distributed around the ovarian cavity in females and on the edge of the gonad in males during the early stages of sex differentiation. Sex reversal was induced via the oral administration of letrozole (LTZ), 17α-methyltestosterone (MT), and 17β-estradiol (E2), respectively, during the labile period of gonadal development. LTZ and MT induced 100% masculinization of genotype-females, whereas E2 induced only 50-60% feminization of genotype-males. Such findings suggest that the fish retained high sexual plasticity despite the existence of the neo-Y chromosome. MT and E2 had negative effect on fish growth, whereas LTZ did not exert such side effect. LTZ and MT could accelerate gonadal development in sex-reversed genotype-males, whereas E2 inhibited gonadal development in genotype-females of rock bream. These findings provide a basis for further research on the mechanisms of sex determination and differentiation in fishes with X1X2Y sex chromosome system and provide a sex reversal protocol for rock bream.

Sex-Specific Genomic Region Identification and Molecular Sex Marker Development of Rock Bream (Oplegnathus fasciatus)

Rock bream (Oplegnathus fasciatus) is a valuable commercial marine teleost species, which exhibits sexual dimorphism in growth performance. However, the absence of a rapid and cost-effective sex identification method based on sex-specific genetic marker has impeded study on sex determination mechanisms and breeding applications. In the present study, we firstly developed the PCR method for identifying potential sex-specific sequences in Oplegnathus fasciatus with the next-generation sequencing. Sex-specific genomic regions/loci for sex determination were discovered on Chr2 and Chr6 by genome-wide association analysis, sequencing depth, and heterozygosity comparison between females and males. Candidate sex-determining genes (CCDC63, ITR, WNT4) were furtherly detected in transcriptome data of testes and ovaries. Taken together, a male-specific 34-bp deletion on the Chr2 was identified and developed into molecular marker of sex for O. fasciatus. After validation in individuals with known phenotypic sexes, the accuracy was 100%. This study gives an insight into the mechanism of sex determination in O. fasciatus, and the gender marker is crucial both for future genomic research and for development of efficient and sustainable aquaculture practice.

Identification and Characterization of Dimorphic Expression of Sex-Related Genes in Rock Bream, a Fish With Multiple Sex Chromosomes

The rock bream (Oplegnathus fasciatus) is a typical fish with a unique multiple sex chromosome system. In this study, we investigated the gene expression profiling in the gonads and brains of both males and females using RNA-Seq to identify sex-related genes and pathways. In accordance with the dimorphic expression profiles, combined with Gene ontology and KEGG enrichment analyses, a number of potential genes and pathways associated with sex determination were obtained from transcriptional analysis, especially some sex-biased genes and pathways. Next, we selected 18 candidate genes and analyzed their expression in different tissues and developmental stages. We found that the expression levels of Amh, Dmrt1, Sox9, Dmrtb1, and Nanos2 were significantly higher in the testis than those in the ovary or other tissues, whereas the expression levels of ZP4, Bouncer, RNF208, FoxH1, and TOB were significantly higher in the ovary than those in the testis. Furthermore, the expression levels of these genes in different developmental stages of gonads also showed sexually dimorphic patterns, suggesting that they might play important roles during gonadal development. These genes are useful markers for investigating sex determination and differentiation in rock bream. The findings of this study can provide insights into the molecular mechanisms of sex determination and differentiation in fish with multiple sex chromosome systems.

First genomic prediction and genome‐wide association for complex growth‐related traits in Rock Bream (Oplegnathus fasciatus)

Rock Bream (Oplegnathus fasciatus) is an important aquaculture species for offshore cage aquaculture and fish stocking of marine ranching in East Asia. Genomic selection has the potential to expedite genetic gain for the key target traits of a breeding program, but has not yet been evaluated in Oplegnathus. The purposes of the present study were to explore the performance of genomic selection to improve breeding value accuracy through real data analyses using six statistical models and to carry out genome‐wide association studies (GWAS) to dissect the genetic architecture of economically vital growth‐related traits (body weight, total length, and body depth) in the O. fasciatus population. After quality control, genotypes for 16,162 SNPs were acquired for 455 fish. Heritability was estimated to be moderate for the three traits (0.38 for BW, 0.33 for TL, and 0.24 for BD), and results of GWAS indicated that the underlying genetic architecture was polygenic. Six statistic models (GBLUP, BayesA, BayesB, BayesC, Bayesian Ridge‐Regression, and Bayesian LASSO) showed similar performance for the predictability of genomic estimated breeding value (GEBV). The low SNP density (around 1 K selected SNP based on GWAS) is sufficient for accurate prediction on the breeding value for the three growth‐related traits in the current studied population, which will provide a good compromise between genotyping costs and predictability in such standard breeding populations advanced. These consequences illustrate that the employment of genomic selection in O. fasciatus breeding could provide advantages for the selection of breeding candidates to facilitate complex economic growth traits.

Evolutionary Dynamics of the Repetitive DNA in the Karyotypes of Pipa carvalhoi and Xenopus tropicalis (Anura, Pipidae)

The large amphibian genomes contain numerous repetitive DNA components that have played an important role in the karyotypic diversification of this vertebrate group. Hypotheses based on the presumable primitive karyotype (2n = 20) of the anurans of the family Pipidae suggest that they have evolved principally through intrachromosomal rearrangements. Pipa is the only South American pipid, while all the other genera are found in Africa. The divergence of the South American lineages from the African ones occurred at least 136 million years ago and is thought to have had a strong biogeographic component. Here, we tested the potential of the repetitive DNA to enable a better understanding of the differentiation of the karyotype among the family Pipidae and to expand our capacity to interpret the chromosomal evolution in this frog family. Our results indicate a long history of conservation in the chromosome bearing the H3 histone locus, corroborating inferences on the chromosomal homologies between the species in pairs 6, 8, and 9. The chromosomal distribution of the microsatellite motifs also provides useful markers for comparative genomics at the chromosome level between Pipa carvalhoi and Xenopus tropicalis, contributing new insights into the evolution of the karyotypes of these species. We detected similar patterns in the distribution and abundance of the microsatellite arrangements, which reflect the shared organization in the terminal/subterminal region of the chromosomes between these two species. By contrast, the microsatellite probes detected a differential arrangement of the repetitive DNA among the chromosomes of the two species, allowing longitudinal differentiation of pairs that are identical in size and morphology, such as pairs 1, 2, 4, and 5. We also found evidence of the distinctive composition of the repetitive motifs of the centromeric region between the species analyzed in the present study, with a clear enrichment of the (CA) and (GA) microsatellite motifs in P. carvalhoi. Finally, microsatellite enrichment in the pericentromeric region of chromosome pairs 6, 8, and 9 in the P. carvalhoi karyotype, together with interstitial telomeric sequences (ITS), validate the hypothesis that pericentromeric inversions occurred during the chromosomal evolution of P. carvalhoi and reinforce the role of the repetitive DNA in the remodeling of the karyotype architecture of the Pipidae.

Deciphering the Origin and Evolution of the X1X2Y System in Two Closely-Related Oplegnathus Species (Oplegnathidae and Centrarchiformes)

Oplegnathus fasciatus and O. punctatus (Teleostei: Centrarchiformes: Oplegnathidae), are commercially important rocky reef fishes, endemic to East Asia. Both species present an X1X2Y sex chromosome system. Here, we investigated the evolutionary forces behind the origin and differentiation of these sex chromosomes, with the aim to elucidate whether they had a single or convergent origin. To achieve this, conventional and molecular cytogenetic protocols, involving the mapping of repetitive DNA markers, comparative genomic hybridization (CGH), and whole chromosome painting (WCP) were applied. Both species presented similar 2n, karyotype structure and hybridization patterns of repetitive DNA classes. 5S rDNA loci, besides being placed on the autosomal pair 22, resided in the terminal region of the long arms of both X1 chromosomes in females, and on the X1 and Y chromosomes in males. Furthermore, WCP experiments with a probe derived from the Y chromosome of O. fasciatus (OFAS-Y) entirely painted the X1 and X2 chromosomes in females and the X1, X2, and Y chromosomes in males of both species. CGH failed to reveal any sign of sequence differentiation on the Y chromosome in both species, thereby suggesting the shared early stage of neo-Y chromosome differentiation. Altogether, the present findings confirmed the origin of the X1X2Y sex chromosomes via Y-autosome centric fusion and strongly suggested their common origin.

Genome sequence of the barred knifejaw Oplegnathus fasciatus (Temminck & Schlegel, 1844): the first chromosome-level draft genome in the family Oplegnathidae

BACKGROUND

The barred knifejaw (Oplegnathus fasciatus), a member of the Oplegnathidae family of the Centrarchiformes, is a commercially important rocky reef fish native to East Asia. Oplegnathus fasciatus has become an important fishery resource for offshore cage aquaculture and fish stocking of marine ranching in China, Japan, and Korea. Recently, sexual dimorphism in growth with neo-sex chromosome and widespread biotic diseases in O. fasciatus have been increasing concern in the industry. However, adequate genome resources for gaining insight into sex-determining mechanisms and establishing genetically resistant breeding systems for O. fasciatus are lacking. Here, we analyzed the entire genome of a female O. fasciatus fish using long-read sequencing and Hi-C data to generate chromosome-length scaffolds and a highly contiguous genome assembly.

FINDINGS

We assembled the O. fasciatus genome with a total of 245.0 Gb of raw reads that were generated using both Pacific Bioscience (PacBio) Sequel and Illumina HiSeq 2000 platforms. The final draft genome assembly was approximately 778.7 Mb, which reached a high level of continuity with a contig N50 of 2.1 Mb. The genome size was consistent with the estimated genome size (777.5 Mb) based on k-mer analysis. We combined Hi-C data with a draft genome assembly to generate chromosome-length scaffolds. Twenty-four scaffolds corresponding to the 24 chromosomes were assembled to a final size of 768.8 Mb with a contig N50 of 2.1 Mb and a scaffold N50 of 33.5 Mb using 1,372 contigs. The identified repeat sequences accounted for 33.9% of the entire genome, and 24 003 protein-coding genes with an average of 10.1 exons per gene were annotated using de novo methods, with RNA sequencing data and homologies to other teleosts. According to phylogenetic analysis using protein-coding genes, O. fasciatus is closely related to Larimichthys crocea, with O. fasciatus diverging from their common ancestor approximately 70.5-88.5 million years ago.

CONCLUSIONS

We generated a high-quality draft genome for O. fasciatus using long-read PacBio sequencing technology, which represents the first chromosome-level reference genome for Oplegnathidae species. Assembly of this genome assists research into fish sex-determining mechanisms and can serve as a resource for accelerating genome-assisted improvements in resistant breeding systems.

Chromosomes of Asian cyprinid fishes: cytogenetic analysis of two representatives of small paleotetraploid tribe Probarbini

Background

Polyploidy, although still poorly explored, represents an important evolutionary event in several cyprinid clades. Herein, Catlocarpio siamensis and Probarbus jullieni - representatives of the paleotetraploid tribe Probarbini, were characterized both by conventional and molecular cytogenetic methods.

Results

Alike most other paleotetraploid cyprinids (with 2n = 100), both species studied here shared 2n = 98 but differed in karyotypes: C. siamensis displayed 18m + 34sm + 46st/a; NF = 150, while P. jullieni exhibited 26m + 14sm + 58st/a; NF = 138. Fluorescence in situ hybridization (FISH) with rDNA probes revealed two (5S) and eight (18S) signals in C. siamensis, respectively, and six signals for both probes in P. jullieni. FISH with microsatellite motifs evidenced substantial genomic divergence between both species. The almost doubled size of the chromosome pairs #1 in C. siamensis and #14 in P. jullieni compared to the rest of corresponding karyotypes indicated chromosomal fusions.

Conclusion

Based on our findings, together with likely the same reduced 2n = 98 karyotypes in the remainder Probarbini species, we hypothesize that the karyotype 2n = 98 might represent a derived character, shared by all members of the Probarbini clade. Besides, we also witnessed considerable changes in the amount and distribution of certain repetitive DNA classes, suggesting complex post-polyploidization processes in this small paleotetraploid tribe.

The testis and ovary transcriptomes of the rock bream (Oplegnathus fasciatus): A bony fish with a unique neo Y chromosome

The rock bream (Oplegnathus fasciatus) is considerably one of the most economically important marine fish in East Asia and has a unique neo-Y chromosome system that is a good model to study the sex determination and differentiation in fish. In the present study, we used Illumina sequencing technology (HiSeq2000) to sequence, assemble and annotate the transcriptome of the testis and ovary tissues of rock bream. A total of 40,004,378 (NCBI SRA database SRX1406649) and 53,108,992 (NCBI SRA database SRX1406648) high quality reads were obtained from testis and ovary RNA sequencing, respectively, and 60,421 contigs (with average length of 1301 bp) were obtained after de novo assembling with Trinity software. Digital gene expression analysis reveals 14,036 contigs that show gender-enriched expressional profile with either testis-enriched (237 contigs) or ovary-enriched (581 contigs) with RPKM > 100. There are 237 male- and 582 female-abundant expressed genes that show sex dimorphic expression. We hope that the gonad transcriptome and those gender-enriched transcripts of rock bream can provide some insight into the understanding of genome-wide transcriptome profile of teleost gonad tissue and give useful information in fish gonad development.

Genomic Organization of Repetitive DNA Elements and Its Implications for the Chromosomal Evolution of Channid Fishes (Actinopterygii, Perciformes)

Channid fishes, commonly referred to as "snakeheads", are currently very important in Asian fishery and aquaculture due to the substantial decline in natural populations because of overexploitation. A large degree of chromosomal variation has been found in this family, mainly through the use of conventional cytogenetic investigations. In this study, we analyzed the karyotype structure and the distribution of 7 repetitive DNA sequences in several Channa species from different Thailand river basins. The aim of this study was to investigate the chromosomal differentiation among species and populations to improve upon the knowledge of its biodiversity and evolutionary history. Rearrangements, such as pericentric inversions, fusions and polyploidization, appear to be important events during the karyotypic evolution of this genus, resulting in the chromosomal diversity observed among the distinct species and even among populations of the same species. In addition, such variability is also increased by the genomic dynamism of repetitive elements, particularly by the differential distribution and accumulation of rDNA sequences on chromosomes. This marked diversity is likely linked to the lifestyle of the snakehead fishes and their population fragmentation, as already identified for other fish species. The karyotypic features highlight the biodiversity of the channid fishes and justify a taxonomic revision of the genus Channa, as well as of the Channidae family as a whole, as some nominal species may actually constitute species complexes.

Structurally Complex Organization of Repetitive DNAs in the Genome of Cobia (Rachycentron canadum)

Repetitive DNAs comprise the largest fraction of the eukaryotic genome. They include microsatellites or simple sequence repeats (SSRs), which play an important role in the chromosome differentiation among fishes. Rachycentron canadum is the only representative of the family Rachycentridae. This species has been focused on several multidisciplinary studies in view of its important potential for marine fish farming. In the present study, distinct classes of repetitive DNAs, with emphasis on SSRs, were mapped in the chromosomes of this species to improve the knowledge of its genome organization. Microsatellites exhibited a diversified distribution, both dispersed in euchromatin and clustered in the heterochromatin. The multilocus location of SSRs strengthened the heterochromatin heterogeneity in this species, as suggested by some previous studies. The colocalization of SSRs with retrotransposons and transposons pointed to a close evolutionary relationship between these repetitive sequences. A number of heterochromatic regions highlighted a greater complex organization than previously supposed, harboring a diversity of repetitive elements. In this sense, there was also evidence of colocalization of active genetic regions and different classes of repetitive DNAs in a common heterochromatic region, which offers a potential opportunity for further researches regarding the interaction of these distinct fractions in fish genomes.

Genomic organization of repetitive DNAs and its implications for male karyotype and the neo-Y chromosome differentiation in Erythrinus erythrinus (Characiformes, Erythrinidae)

Studies have demonstrated the effective participation of repetitive DNA sequences in the origin and differentiation of the sex chromosomes in some biological groups. In this study several microsatellites and retrotranposable sequences were cytogenetically mapped in the Erythrinus erythrinus (Bloch & Schneider, 1801) male genome (karyomorph C), focusing on the distribution of these sequences in the sex chromosomes and in the evolutionary processes related to their differentiation. Males of E. erythrinus - karyomorph C - present 2n = 51 chromosomes (7m + 2sm + 6st + 36a), including the X1X2Y sex chromosomes. The C-positive heterochromatin has a predominant localization on the centromeric region of most chromosome pairs, but also in some telomeric regions. The 5S rDNA sites are located in the centromeric region of 27 chromosomes, including 26 acrocentric ones and the metacentric Y chromosome. The retrotransposons Rex 1 and Rex 6 show a dispersed pattern in the karyotype, contrasting with the Rex 3 distribution which is clearly co-localized with all the 27 5S rDNA sites. The microsatellite sequences show a differential distribution, some of them restricted to telomeric and/or interstitial regions and others with a scattered distribution on the chromosomes. However, no preferential accumulation of these elements were observed in the neo-Y chromosome, in contrast to what usually occurs in simple sex chromosome systems.