Genome-wide association study reveals genes associated with the absence of intermuscular bones in tambaqui (Colossoma macropomum)

Screening for IBs-relative genes by transcriptome analysis and generation IBs-less mutants in Culter alburnus

Intermuscular bones (IBs), distributed specifically in the myosepta on both sides of lower teleosts, negatively affect palatability and processing. Recent research in zebrafish and several economically important farmed fishes has led to the breakthrough discovery of the mechanism of IBs formation and generation of IBs-loss mutants. This study explored the ossification patterns of IBs in juvenile Culter alburnus. Besides, some key genes and bone-related signaling pathways were identified by transcriptomic data. Furthermore, PCR microarray validation revealed that claudin1 could potentially regulate IBs formation. Additionally, we created several IBs-reduced mutants of C. alburnus by loss of the function of bone morphogenetic proteins 6 (bmp6) gene using CRISPR/Cas9 editing. These results suggested that CRISPR/Cas9-mediated bmp6 knockout was promising approach for breeding IBs-free strain in other cyprinids.

Transcriptomic analysis of intermuscular bone development in barbel steed (Hemibarbus labeo)

Intermuscular bones (IBs), which are little, bony spicules in muscle, are embedded in lower teleosts' myosepta. Despite the importance of studying IB development in freshwater aquaculture species, the genes associated with IB development need to be further explored. In the present study, we identified four stages of IB development in barbel steed (Hemibarbus labeo), namely stage 1: IBs have not emerged, stage 2: a few small IBs have emerged in the tail, stage 3: longer IBs gradually emerged in the tail and stage 4: all of the IBs in the tail are mature and long, via Alizarin red staining. Subsequently, we used the HiseqXTen platform to sequence and de novo assemble the transcriptome of epaxial muscle (between 35th and 40th myomere) of barbel steed at 29 days (stage 1) and 42 days (stage 3) after hatching. A total of 190,814 unigenes were obtained with an average length and N50 of 648 bp and 1027 bp, respectively. We found 2174 differentially expressed genes (DEGs) between stages 1 and 3, of which 378 and 1796 were up- and down-regulated, respectively. Functional enrichment analysis showed that several DEGs functioned in ossification, positive regulation of osteoblast differentiation, osteoblast differentiation, and BMP signaling pathway, and were further enriched in signal pathway, including osteoclast differentiation, TGF-β signaling pathway, cytokine-cytokine receptor interaction, Jak-STAT signaling pathway, and other KEEG pathways. In conclusion, we identified genes that may be related to IB development, such as kazal type serine peptidase inhibitor domain 1 (KAZALD1), extracellular matrix protein 1 (ECM1), tetranectin, bone morphogenetic protein 1 (bmp1), acid phosphatase 5 (ACP5), collagen type XI alpha 1 chain (COL11A1), matrix metallopeptidase 9 (MMP9), pannexin-3 (PANX3), sp7 transcription factor (Sp7), and c-x-c motif chemokine ligand 8 (CXCL8), by comparing the transcriptomes of epaxial muscle before and after IB ossification. This study provided a theoretical basis for identifying the molecular mechanisms underlying IB development in fish.

Single-cell transcriptomes and runx2b-/- mutants reveal the genetic signatures of intermuscular bone formation in zebrafish

Intermuscular bones (IBs) are mineralized spicules, present in the myosepta of many, but not all, teleost species. IBs are often small and sharp, and they consequently limit how the fish can be processed; the IBs may cause injury or trauma if lodged in consumers' throats or mouths, and therefore affect the appeal of the fish to many consumers. The development of IBs in teleosts is still not fully understood and the molecular basis of IB development remains to be established. Here, the characteristics of IB tissue are evaluated based on single-cell transcriptomics in wild-type zebrafish. The analysis defined 18 distinct cell types. Differentiation trajectories showed that IBs are derived from tendons and that a core tendon-osteoblast cell lineage is related to IB formation. In particular, the functions of 10 candidate genes were evaluated via CRISPR-Cas9 mutants. Among those, runx2b^-/- mutants completely lost IBs, while swimming performance, growth and bone mineral density were not significantly different from runx2b^+/+ zebrafish. Comparative single-cell RNA sequencing (scRNA-seq) analysis in runx2b^-/- and runx2b^+/+ zebrafish revealed the role of osteoblasts in IB formation. In addition, differentially expressed genes were enriched in the transforming growth factor β/bone morphogenetic protein (TGF-β/BMP) pathway after runx2b deletion. This study provides evidence for the crucial role of runx2b regulation in IB formation. Genetic breeding can target runx2b regulation and generate strains of commercial fish species without IBs, which can improve the safe consumption and economic value of many farmed fish species.

Bulked Segregant Analysis and Association Analysis Identified the Polymorphisms Related to the Intermuscular Bones in Common Carp (Cyprinus carpio)

Simple Summary

Many widely cultured freshwater fish species, such as common carp, belong to the Cyprinidae family. However, most cyprinids have numerous and complex intermuscular bones (IBs), resulting in an adverse effect on cyprinid fish meat processing and consumption. Numerous studies have been trying to understand the development mechanism of IBs and to identify the SNPs associated with the total IB number. However, the SNPs associated with different forms of IBs have been studied less thoroughly. The joint effects of the SNPs on IB development also remain poorly understood. The common carp has numerous geographical populations and domesticated strains, diversifying its phenotypes. The question of whether consensus IB-related SNPs or genes exist among multiple strains of common carp has also not yet been answered. Selective breeding of IB-reduced common carp has been hindered due to a lack of effective molecular markers. To answer these questions, we performed bulked segregant analysis (BSA) to detect the consensus SNPs in three strains. The consensus BSA-SNPs and the other SNPs in their flanking regions were validated in additional individuals. The SNPs associated with the frequency of different IB types were identified. We examined the joint effects of significant SNPs on the numbers of different types of IBs. The identified genetic markers may benefit future selective breeding and reduce the IB number in common carp.

Abstract

The allotetraploid common carp is one of the most important freshwater food fish. However, the IBs found in allotetraploid common carp increase the difficulty in fish meat processing and consumption. Although candidate genes associated with the total IB number have been identified, the SNPs associated with the numbers of the total IBs and different forms of IBs have not yet been identified, hindering the breeding of IB-reduced common carp. Herein, the numbers of different types of IBs in three common carp strains were measured. Using whole-genome resequencing and bulked segregant analysis in three pairs of IB-more and IB-less groups, we identified the consensus nonsynonymous SNPs in three strains of common carp. Screening the flanking regions of these SNPs led to the detection of other SNPs. Association study detected 21 SNPs significantly associated with the number of total IBs, epineural-IBs, and ten detailed types of IBs. We observed the joint effects of multiple SNPs on each associated IB number with an improved explained percentage of phenotypic variation. The resulting dataset provides a resource to understand the molecular mechanisms of IB development in different common carp strains. These SNPs are potential markers for future selection to generate IB-reduced common carp.

GBS-MeDIP: A protocol for parallel identification of genetic and epigenetic variation in the same reduced fraction of genomes across individuals

The GBS-MeDIP protocol combines two previously described techniques, Genotype-by-Sequencing (GBS) and Methylated-DNA-Immunoprecipitation (MeDIP). Our method allows for parallel and cost-efficient interrogation of genetic and methylomic variants in the DNA of many reduced genomes, taking advantage of the barcoding of DNA samples performed in the GBS and the subsequent creation of DNA pools, then used as an input for the MeDIP. The GBS-MeDIP is particularly suitable to identify genetic and methylomic biomarkers when resources for whole genome interrogation are lacking.

•

Our protocol combines two well-known techniques, GBS and MeDIP.

•

The GBS-MeDIP allows for parallel interrogation of genetic/methylomic variants.

•

Allows for cost-efficient analysis of reduced genomes from multiple individuals.

•

Ideal for the identification of genetic and/or methylomic markers at reduced cost.

Genomic selection signatures in farmed Colossoma macropomum from tropical and subtropical regions in South America

Tambaqui or cachama (Colossoma macropomum) is one of the most important neotropical freshwater fish used for aquaculture in South America, and its production is concentrated at low latitudes (close to the Equator, 0°), where the water temperature is warm. Therefore, understanding how selection shapes genetic variations and structure in farmed populations is of paramount importance in evolutionary biology. High‐throughput sequencing to generate genome‐wide data for fish species allows for elucidating the genomic basis of adaptation to local or farmed conditions and uncovering genes that control the phenotypes of interest. The present study aimed to detect genomic selection signatures and analyze the genetic variability in farmed populations of tambaqui in South America using single‐nucleotide polymorphism (SNP) markers obtained with double‐digest restriction site‐associated DNA sequencing. Initially, 199 samples of tambaqui farmed populations from different locations (located in Brazil, Colombia, and Peru), a wild population (Amazon River, Brazil), and the base population of a breeding program (Aquaculture Center, CAUNESP, Jaboticabal, SP, Brazil) were genotyped. Observed and expected heterozygosity was 0.231–0.350 and 0.288–0.360, respectively. Significant genetic differentiation was observed using global FST analyses of SNP loci (FST = 0.064, p < 0.050). Farmed populations from Colombia and Peru that differentiated from the Brazilian populations formed distinct groups. Several regions, particularly those harboring the genes of significance to aquaculture, were identified to be under positive selection, suggesting local adaptation to stress under different farming conditions and management practices. Studies aimed at improving the knowledge of genomics of tambaqui farmed populations are essential for aquaculture to gain deeper insights into the evolutionary history of these fish and provide resources for the establishment of breeding programs.

Development of a multi-species SNP array for serrasalmid fish Colossoma macropomum and Piaractus mesopotamicus

Scarce genomic resources have limited the development of breeding programs for serrasalmid fish Colossoma macropomum (tambaqui) and Piaractus mesopotamicus (pacu), the key native freshwater fish species produced in South America. The main objectives of this study were to design a dense SNP array for this fish group and to validate its performance on farmed populations from several locations in South America. Using multiple approaches based on different populations of tambaqui and pacu, a final list of 29,575 and 29,612 putative SNPs was selected, respectively, to print an Axiom AFFYMETRIX (THERMOFISHER) SerraSNP array. After validation, 74.17% (n = 21,963) and 71.25% (n = 21,072) of SNPs were classified as polymorphic variants in pacu and tambaqui, respectively. Most of the SNPs segregated within each population ranging from 14,199 to 19,856 in pacu; and from 15,075 to 20,380 in tambaqui. Our results indicate high levels of genetic diversity and clustered samples according to their hatchery origin. The developed SerraSNP array represents a valuable genomic tool approaching in-depth genetic studies for these species.

Molecular mechanisms of intermuscular bone development in fish: a review

Intermuscular bones (IBs) are slender linear bones embedded in muscle, which ossify from tendons through a process of intramembranous ossification, and only exist in basal teleosts. IBs are essential for fish swimming, but they present a choking risk during human consumption, especially in children, which can lead to commercial risks that have a negative impact on the aquaculture of these fish. In this review, we discuss the morphogenesis and functions of IBs, including their underlying molecular mechanisms, as well as the advantages and disadvantages of different methods for IB studies and techniques for breeding and generating IB-free fish lines. This review reveals that the many key genes involved in tendon development, osteoblast differentiation, and bone formation, e.g., scxa, msxC, sost, twist, bmps, and osterix, also play roles in IB development. Thus, this paper provides useful information for the breeding of new fish strains without IBs via genome editing and artificial selection.