Genetic Mapping of Head Size Related Traits in Common Carp (Cyprinus carpio)

Genome-wide association study analysis to resolve the key regulatory mechanism of biomineralization in Pinctada fucata martensii

Biomineralization-controlled exo-/endoskeleton growth contributes to body growth and body size diversity. Molluscan shells undergo ectopic biomineralization to form the exoskeleton and biocalcified "pearl" involved in invading defence. Notably, exo-/endoskeletons have a common ancestral origin, but their regulation and body growth are largely unknown. This study employed the pearl oyster, Pinctada fucata marntensii, a widely used experimental model for biomineralization in invertebrates, to perform whole-genome resequencing of 878 individuals from wild and breeding populations. This study characterized the genetic architecture of biomineralization-controlled growth and ectopic biomineralization. The insulin-like growth factor (IGF) endocrine signal interacted with ancient single-copy transcription factors to form the regulatory network. Moreover, the "cross-phylum" regulation of key long noncoding RNA (lncRNA) in bivalves and mammals indicated the conserved genetic and epigenetic regulation in exo-/endoskeleton growth. Thyroid hormone signal and apoptosis regulation in pearl oysters affected ectopic biomineralization in pearl oyster. These findings provide insights into the mechanism underlying the evolution and regulation of biomineralization in exo-/endoskeleton animals and ectopic biomineralization.

Genome-Wide Association Study Reveals SNPs and Candidate Genes Related to Growth and Body Shape in Bighead Carp (Hypophthalmichthys nobilis)

Growth is an economically important trait in bighead carp and other aquaculture species that affects production efficiency. Interestingly, the head of the bighead carp has a high market value in China; therefore, it is important to study the genetic bases of both growth and body shape traits. A genome-wide association study was performed based on 2b-RAD sequencing of 776 individuals to identify SNPs associated with growth and body shape traits, including body weight, body length, body height, and deheaded body length. In total, 26 significant and 19 suggestive SNPs were identified, and more than half of these significant SNPs were clustered in LG16. Two LGs (LG16 and LG21) contained QTLs associated with body weight. Fourteen SNPs of LG16 and two LG21 SNPs were found to be associated with body length. For body height, 12 significantly associated SNPs were identified in LG16. Additionally, 12 SNPs of LG16 and 3 SNPs of LG21 were found to be associated with deheaded body length. Forty-three genes were significantly or suggestively associated with body shape/growth traits based on GWAS results, 18 of which were candidate genes for all BW, BL, BH, and DBL traits. One of these genes, fndc5b, was selected for further analyses. Association analysis revealed that one SNP (g.245 C > T) in the introns of fndc5b was significantly associated with growth-related traits in growth-extreme samples. The mRNA levels of fndc5b in the brains of the lightweight group were significantly higher than those of the heavy-weight group. This study helps to reveal the genetic structure of growth and body development in fish and provides candidate genes for future molecular marker-assisted selection for fast growth and better body conformation in bighead carp.

Comparative Transcriptome Analysis of Organ-Specific Adaptive Responses to Hypoxia Provides Insights to Human Diseases

The common carp is a hypoxia-tolerant fish, and the understanding of its ability to live in low-oxygen environments has been applied to human health issues such as cancer and neuron degeneration. Here, we investigated differential gene expression changes during hypoxia in five common carp organs including the brain, the gill, the head kidney, the liver, and the intestine. Based on RNA sequencing, gene expression changes under hypoxic conditions were detected in over 1800 genes in common carp. The analysis of these genes further revealed that all five organs had high expression-specific properties. According to the results of the GO and KEGG, the pathways involved in the adaptation to hypoxia provided information on responses specific to each organ in low oxygen, such as glucose metabolism and energy usage, cholesterol synthesis, cell cycle, circadian rhythm, and dopamine activation. DisGeNET analysis showed that some human diseases such as cancer, diabetes, epilepsy, metabolism diseases, and social ability disorders were related to hypoxia-regulated genes. Our results suggested that common carp undergo various gene regulations in different organs under hypoxic conditions, and integrative bioinformatics may provide some potential targets for advancing disease research.

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.

Dynamic mRNA and miRNA expression of the head during early development in bighead carp (Hypophthalmichthys nobilis)

Background

Head of fish species, an exquisitely complex anatomical system, is important not only for studying fish evolution and development, but also for economic values. Currently, although some studies have been made on fish growth and body shapes, very limited information is available on the molecular mechanism of head development.

Results

In this study, RNA sequencing (RNA–Seq) and small RNA sequencing (sRNA–Seq) technologies were used to conduct integrated analysis for the head of bighead carp at different development stages, including 1, 3, 5, 15 and 30 Dph (days post hatch). By RNA-Seq data, 26 pathways related to growth and bone formation were identified as the main physiological processes during early development. Coupling this to sRNA–Seq data, we picked out six key pathways that may be responsible for head development, namely ECM receptor interaction, TNF signaling pathway, osteoclast differentiation, PI3K–Akt signaling pathway, Neuroactive ligand–receptor interaction and Jak–STAT signaling pathway. Totally, 114 important candidate genes from the six pathways were obtained. Then we found the top 20 key genes according to the degree value by cytohubba, which regulated cell growth, skeletal formation and blood homeostasis, such as pik3ca, pik3r1, egfr, vegfa, igf1 and itga2b. Finally, we also acquired 19 key miRNAs playing multiple roles in the perfection of various tissues in the head (such as brain, eye and mouth) and mineralization of head bone system, such as let–7e, miR–142a–5p, miR–144–3p, miR–23a–3p and miR–223.

Conclusions

Results of this study will be informative for genetic mechanisms of head development and also provide potential candidate targets for the interaction regulation during early growth in bighead carp.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08387-x.

Construction of Three High-Density Genetic Linkage Maps and Dynamic QTL Mapping of Growth Traits in Yellow River Carp (Cyprinus carpio haematopterus)

To provide the theoretical basis for researching growth, development, and molecular marker-assisted breeding of the economically important Yellow River carp (Cyprinus carpio haematopterus) using dynamic quantitative trait locus (QTL) mapping, we constructed three genetic linkage maps from 207 progeny using a new modified genotyping-by-sequencing method. The three maps contained 16,886, 16,548, and 7482 single nucleotide polymorphism markers, respectively, with an average interval of 0.36 cM, 0.45 cM, and 1.00 cM. We identified 148 QTLs related to four growth traits that were located on 25 chromosomes from three growth stages of Yellow River carp. A total of 32, 36, 43, and 37 QTLs were associated with body length, height, width, and weight, respectively. Among them, 47 QTLs were detected for only one growth trait in one stage, but all of the other QTLs were co-localized. Of the 14 main QTLs, 13 were located on chromosome 12, which suggests the presence of growth-related genes on this chromosome. We then detected 17 candidate genes within 50 K upstream and downstream of the 14 main QTLs. This is the first report of the dynamic QTL mapping of growth traits of Yellow River carp, and the results can be used in future studies of growth, development, and molecular-assisted breeding of this species.

Genome-Wide Association Study of Body Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker (Larimichthys crocea) is one of the most important cultured marine fish on the southeast coast of China. Its body shape is important for the aquaculture industry since it affects the behavior such as swimming, ingesting, and evading, as well as customer preference. Due to the greater consumer demand of small head, slender body large yellow croaker, selecting and breeding of slender individuals with the assistance of genetic markers will benefit the industry quickly. In this study, several traits were employed to represent body shape, including body depth/body length (BD/BL), body thickness/body length (BT/BL), caudal peduncle depth/caudal peduncle length (CPDLR), tail length/body length (TL/BL), and body area/head area (BA/HA). Genome-wide association study was conducted with a panmictic population of 280 individuals to identify SNP and genes potentially associated with body shape. A set of 20 SNPs on 12 chromosomes were identified to be significantly associated with body shape-related traits. Besides, 5 SNPs were identified to be suggestive associated with CPDLR and BT/BL. Surrounding these SNPs, we found some body shape-related candidate genes, including fabp1, acrv1, bcor, mstn, bambi, and neo1, which involved in lipid metabolism, TGF-β signaling, and BMP pathway and other important regulatory pathways. These results will be useful for the understanding of the genetic basis of body shape formation and helpful for body shape controlling of large yellow croaker by using marker-assisted selection.

Genome Wide Analysis for Growth at Two Growth Stages in A New Fast-Growing Common Carp Strain (Cyprinus carpio L.)

In order to identify candidate genes or loci associated with growth performance of the newly established common carp strain, Xinlong, we conducted a genome-wide association analysis using 2b-RAD technology on 123 individuals. We constructed two sets of libraries associated with growth-related parameters (weight, length, width and depth) measured at two different grow-out stages. Among the 413,059 SNPs identified using SOAP SNP calling, 147,131 were tested for GWAS after quality filtering. Finally, 39 overlapping SNPs, assigned to four genomic locations, were associated with growth traits in two stages. These loci were assigned to functional classes related to immune response, response to stress, neurogenesis, cholesterol metabolism and development, and proliferation and differentiation of cells. By overlapping results of Plink and EMMAX analyses, we identified three genes: TOX, PLK2 and CD163 (both methods P < 0.05). Our study results could be used for marker-assisted selection to further improve the growth of the Xinlong strain, and illustrate that largely different sets of genes drive the growth of carp in the early and late grow-out stages.

Construction of a High-Density Genetic Linkage Map and QTL Mapping for Growth-Related Traits in Takifugu bimaculatus

Takifugu bimaculatus is a euryhaline species, distributed ranging from the southern Yellow Sea to the South China Sea. Their tolerance to a wide range of salinity and temperature, coupled with a desirable firm texture, makes T. bimaculatus a strong candidate for Takifugu aquaculture in subtropics areas. Due to the increasing demand in markets and emerging of the Takifugu aquaculture industry, close attention has been paid to improvement on the T. bimaculatus production. In aquaculture, the great effort has been put into marker-assisted selective breeding, and efficient improvement was realized. However, few genetic resources on T. bimaculatus are provided so far. Aiming at understanding the genetic basis underlying important economic growth traits, facilitating genetic improvement and enriching the genetic resource in T. bimaculatus, we constructed the first genetic linkage map for T. bimaculatus via double digestion restriction-site association DNA sequencing and conducted quantitative traits locus (QTL) mapping for growth-related traits. The map comprised 1976 single nucleotide polymorphism markers distributed on 22 linkage groups (LG), with a total genetic distance of 2039.74 cM. Based on the linkage map, a chromosome-level assembly was constructed whereby we carried out comparative genomics analysis, verifying the high accuracy on contigs ordering of the linkage map. On the other hand, 18 QTLs associated with growth traits were detected on LG6, LG7, LG8, LG10, LG20, and LG21 with phenotypical variance ranging from 15.1 to 56.4%. Candidate genes participating in cartilage development, fat accumulation, and other growth-related regulation activities were identified from these QTLs, including col11a1, foxa2, and thrap3. The linkage map provided a solid foundation for chromosomes assembly and refinement. QTLs reported here unraveled the genomic architecture of some growth traits, which will advance the investigation of aquaculture breeding efforts in T. bimaculatus.

High-Density Linkage Map and Mapping for Sex and Growth-Related Traits of Largemouth Bass (Micropterus salmoides)

The largemouth bass is an important species, and its culture has risen sharply with the surge in fish aquaculture in China. Due to the lack of selective breeding technology for the largemouth bass, the growth rate and disease resistance are low, its sexual maturation is slow, and other serious problems are contributing to a sharp decline in the safety and quality of largemouth bass products in recent decades. Therefore, comprehensive breeding programs to improve the economic performance and promote the modern industrial development of largemouth bass must be considered a priority. Here, a total of 152 adult largemouth bass, including two parents and 150 progenies, were selected to produce the genetic mapping family. Then, a high-density linkage map was constructed based on restriction site–associated DNA sequencing using 6,917 single-nucleotide polymorphisms (SNPs) located in 24 linkage groups (LGs). The total genetic length of the linkage map was 1,261.96 cM, and the length of each LG varied from 24.72 cM for LG02 to 117.53 cM for LG16, with an average length of 52.58 cM and an average SNP number of 286. Thirteen significant quantitative trait loci (QTLs) for sex determination were located on LG04, LG05, LG08, LG12, LG15, LG21, and LG23. An informative QTL cluster that included six QTLs was detected on LG12. However, one notable QTL, which accounted for 71.48% of the total phenotypic variation, was located in the region of 1.85 cM on LG05. In addition, 32 identified QTLs were related to growth, including body weight, body length, body height, and head length. The QTLs for these growth-related traits are located in 13 LG regions and have little effect on phenotypic variation. This high-density genetic linkage map will enable the fine-mapping of economic traits and support the future genome assembly of the largemouth bass. Additionally, our study will be useful for future selective culture of largemouth bass and could potentially be used in molecular-assisted breeding of largemouth bass for aquaculture.

Genome-Wide Association Study of Growth and Body-Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea) Using ddRAD Sequencing

Large yellow croaker (Larimichthys crocea) is an economically important marine fish species of China. Due to overfishing and marine pollution, the wild stocks of this croaker have collapsed in the past decades. Meanwhile, the cultured croaker is facing the difficulties of reduced genetic diversity and low growth rate. To explore the molecular markers related to the growth traits of croaker and providing the related SNPs for the marker-assisted selection, we used double-digest restriction-site associated DNA (ddRAD) sequencing to dissect the genetic bases of growth traits in a cultured population and identify the SNPs that associated with important growth traits by GWAS. A total of 220 individuals were genotyped by ddRAD sequencing. After quality control, 27,227 SNPs were identified in 220 samples and used for GWAS analysis. We identified 13 genome-wide significant associated SNPs of growth traits on 8 chromosomes, and the beta P of these SNPs ranged from 0.01 to 0.86. Through the definition of candidate regions and gene annotation, candidate genes related to growth were identified, including important regulators such as fgf18, fgf1, nr3c1, cyp8b1, fabp2, cyp2r1, ppara, and ccm2l. We also identified SNPs and candidate genes that significantly associated with body shape, including bmp7, col1a1, col11a2, and col18a1, which are also economically important traits for large yellow croaker aquaculture. The results provided insights into the genetic basis of growth and body shape in large yellow croaker population and would provide reliable genetic markers for molecular marker-assisted selection in the future. Meanwhile, the result established a basis for our subsequent fine mapping and related gene study.

Constructing a High-Density Genetic Linkage Map for Large Yellow Croaker (Larimichthys crocea) and Mapping Resistance Trait Against Ciliate Parasite Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is the most economically important marine cage-farming fish in China in the past decade. However, the sustainable development of large yellow croaker aquaculture has been severely hampered by several diseases, of which, the white spot disease caused by ciliate protozoan parasite Cryptocaryon irritans ranks the most damaging disease in large yellow croaker cage farms. To better understand the genetic basis of parasite infection and disease resistance to C. irritans, it is vital to map the traits and localize the underlying candidate genes in L. crocea genome. Here, we constructed a high-density genetic linkage map using double-digest restriction-site associated DNA (ddRAD)-based high-throughput SNP genotyping data of a F1 mapping family, which had been challenged with C. irritans for resistant trait measure. A total of 5261 SNPs was grouped and oriented into 24 linkage groups (LGs), representing 24 chromosomes of L. crocea. The total genetic map length was 1885.67 cM with an average inter-locus distance of 0.36 cM. Quantitative trait loci (QTL) mapping identified seven significant QTLs in four LGs linked to C. irritans disease resistance. Candidate genes underlying disease resistance were identified from the reference genome, including ifnar1, ifngr2, ikbke, and CD112. Comparative genomic analysis between large yellow croaker and the four closely related species revealed high evolutionary conservation of chromosomes, though inter-chromosomal rearrangements do exist. Especially, the croaker genome structure was closer to the medaka genome than stickleback, indicating that the croaker genome might retain the teleost ancestral genome structure. The high-density genetic linkage map provides an important tool and resource for fine mapping, comparative genome analysis, and molecular selective breeding of large yellow croaker.