A second-generation genetic linkage map for bighead carp (Aristichthys nobilis) based on microsatellite markers

Construction of a high-density genetic map and QTL localization of body weight and wool production related traits in Alpine Merino sheep based on WGR

BACKGROUND

The Alpine Merino is a new breed of fine-wool sheep adapted to the cold and arid climate of the plateau in the world. It has been popularized in Northwest China due to its superior adaptability as well as excellent production performance. Those traits related to body weight, wool yield, and wool fiber characteristics, which are economically essential traits in Alpine Merino sheep, are controlled by QTL (Quantitative Trait Loci). Therefore, the identification of QTL and genetic markers for these key economic traits is a critical step in establishing a MAS (Marker-Assisted Selection) breeding program.

RESULTS

In this study, we constructed the high-density genetic linkage map of Alpine Merino sheep by sequencing 110 F1 generation individuals using WGR (Whole Genome Resequencing) technology. 14,942 SNPs (Single Nucleotide Polymorphism) were identified and genotyped. The map spanned 2,697.86 cM, with an average genetic marker interval of 1.44 cM. A total of 1,871 high-quality SNP markers were distributed across 27 linkage groups, with an average of 69 markers per LG (Linkage Group). Among them, the smallest genetic distance is 19.62 cM for LG2, while the largest is 237.19 cM for LG19. The average genetic distance between markers in LGs ranged from 0.24 cM (LG2) to 3.57 cM (LG17). The marker density in the LGs ranged from LG14 (39 markers) to LG1 (150 markers).

CONCLUSIONS

The first genetic map of Alpine Merino sheep we constructed included 14,942 SNPs, while 46 QTLs associated with body weight, wool yield and wool fiber traits were identified, laying the foundation for genetic studies and molecular marker-assisted breeding. Notably, there were QTL intervals for overlapping traits on LG4 and LG8, providing potential opportunities for multi-trait co-breeding and further theoretical support for selection and breeding of ultra-fine and meaty Alpine Merino sheep.

Updated Genome Assembly of Bighead Carp (Hypophthalmichthys nobilis) and Its Differences Between Male and Female on Genomic, Transcriptomic, and Methylation Level

Cyprinidae is one of the largest family in freshwater fishes, and it is most intensively cultured fish taxon of the world. However, studies about sex determination in this large family is still rear, and one of the reasons is lack of high quality and complete genome. Here, we used nanopore to sequence the genome of a male bighead carp, obtaining contig N50 = 24.25 Mb, which is one of the best assemblies in Cyprinidae. Five males and five females were re-sequenced, and a male-specific region on LG19 was confirmed. We find this region holds many male-specific markers in other Cyprinidae fishes, such as grass carp and silver carp. Transcriptome analyses of hypothalamus and pituitary tissues showed that several sex-specific differentially expressed genes were associated with steroid biosynthesis. The UCH64E gene, located in the male-specific region on LG19, showed higher expression levels in male than female tissues of bighead carp. The methyl-RAD of hypothalamus tissues between males and females indicated that the sexual methylation differences are significant in bighead carp. We also compared the methylation sites recognized using methyl-RAD and nanopore raw reads and found that approximately 73% of the methylation sites identified using methyl-RAD were within nanopore CpG sites.

Molecular cloning and functional characterization of the hypoxia-inducible factor-1α in bighead carp (Aristichthys nobilis)

HIF-l is the earliest documented and most widely studied hypoxia-inducible factor (HIF) and plays a key role in the cell hypoxia signal transduction pathway. Particularly, the HIF-1α protein is sensitive to oxygen and plays a critical role in hypoxia regulation. This study is the first to report on the molecular cloning and characterization of HIF-1α in bighead carp (Aristichthys nobilis; anHIF-1α). The full-length cDNA of anHIF-1α was 2361 bp, and encodes an estimated 674 amino acids with a predicted molecular mass of 76.10 kDa and a theoretical isoelectric point of 7.72. Moreover, the conserved basic Helix-Loop-Helix domain along with two Per-ARNT-Sim domains (A/B), and C-TAD were identified in this protein. Interestingly, the tertiary structure of the anHIF-1α protein was found to be extremely similar to that of mice. Multiple comparison and phylogenetic tree results demonstrated that anHIF-1α was highly conserved. Under normoxic conditions, anHIF-1α mRNA transcripts could be detected in all tissues examined with the highest expression level in the heart. With gradually decreasing oxygen concentrations, anHIF-1α mRNA level was upregulated significantly in the gill, liver, kidney, spleen, intestine, brain, and muscle tissues (P < 0.05). Similarly, anHIF-1α was expressed in all examined bighead carp tissues, and the results suggested that the upregulation of anHIF-1α at the transcriptional level may be an important stress response adaptation to hypoxia in bighead carp. Finally, based on the tertiary structure comparative analyses between anHIF-1α with mouse HIF-1α, we think the physiological function, and protein structure of HIF-1α could be compared between fish and mammal in the future.

A high-density genetic linkage map for Chinese perch (Siniperca chuatsi) using 2.3K genotyping-by-sequencing SNPs

Chinese perch, Siniperca chuatsi (Basilewsky), is one of the most commercially important cultured fishes in China. In the present study, a high-density genetic linkage map of Chinese perch was constructed by genotyping-by-sequencing technique with an F1 mapping panel containing 190 progenies. A total of 2328 SNPs were assigned to 24 linkage groups (LGs), agreeing with the chromosome haploid number in this species (n = 24). The sex-averaged map covered 97.9% of the Chinese perch genome, with the length of 1694.3 cM and a marker density of 0.7 cM/locus. The number of markers per LG ranged from 57 to 222, with a mean of 97. The length of LGs varied from 43.2 to 108.2 cM, with a mean size of 70.6 cM. The recombination rate of females was 1.5:1, which was higher than that of males. To better understand the distribution pattern of segregation distortion between the two sexes of Chinese perch, the skewed markers were retained and used to reconstruct the sex-specific maps. The 16 segregation distortion regions were identified on 10 LGs of the female map, while 12 segregation distortion regions on eight LGs of the male map. Among these LGs, six LGs matched between the sex-specific maps. This high-density linkage map could provide a solid basis for identifying QTL associated with economically important traits, and for implementing marker-assisted selection breeding of Chinese perch.

Construction of High-Resolution RAD-Seq Based Linkage Map, Anchoring Reference Genome, and QTL Mapping of the Sex Chromosome in the Marine Medaka Oryzias melastigma

Medaka (Oryzias sp.) is an important fish species in ecotoxicology and considered as a model species due to its biological features including small body size and short generation time. Since Japanese medaka Oryzias latipes is a freshwater species with access to an excellent genome resource, the marine medaka Oryzias melastigma is also applicable for the marine ecotoxicology. In genome era, a high-density genetic linkage map is a very useful resource in genomic research, providing a means for comparative genomic analysis and verification of de novo genome assembly. In this study, we developed a high-density genetic linkage map for O. melastigma using restriction-site associated DNA sequencing (RAD-seq). The genetic map consisted of 24 linkage groups with 2,481 single nucleotide polymorphism (SNP) markers. The total map length was 1,784 cM with an average marker space of 0.72 cM. The genetic map was integrated with the reference-assisted chromosome assembly (RACA) of O. melastigma, which anchored 90.7% of the assembled sequence onto the linkage map. The values of complete Benchmarking Universal Single-Copy Orthologs were similar to RACA assembly but N50 (23.74 Mb; total genome length 779.4 Mb; gap 5.29%) increased to 29.99 Mb (total genome length 778.7 Mb; gap 5.2%). Using MapQTL analysis with SNP markers, we identified a major quantitative trait locus for sex traits on the Om10. The integration of the genetic map with the reference genome of marine medaka will serve as a good resource for studies in molecular toxicology, genomics, CRISPR/Cas9, and epigenetics.

Comparative analysis of spleen transcriptome detects differences in evolutionary adaptation of immune defense functions in bighead carp and silver carp

The evolutionary divergence of the immune defense functions in bighead carp (A. nobilis) and silver carp (H. molitrix) is still not understood at the molecular level. Here, we obtained 48,821,754 and 55,054,480 clean reads from spleen tissue libraries prepared for bighead carp and silver carp using Illumina paired-end sequencing technology, respectively, and identified 4976 orthologous genes from the transcriptome data sets by comparative analysis. Adaptive evolutionary analysis showed that 212 orthologous genes and 255 Gene Ontology (GO) terms were subjected to positive selection(Ka/Ks values > 1) only in bighead carp, and 195 orthologous genes and 309 GO terms only in silver carp. Among immune defense functions with significant evolutionary divergence, the positively selected biological processes in bighead carp mainly included B cell-mediated immunity, chemokine-mediated signaling pathway, and immunoglobulin mediated immune response, whereas those in silver carp mainly included the antigen processing and presentation, defense response to fungus, and detection of bacteria. Moreover, we found 2974 genes expressed only in spleen of bighead carp and 3494 genes expressed only in spleen of silver carp, where these genes were mostly enriched in the same biological processes or pathways. These results provide a better understanding of the differences in resistance to some diseases by bighead carp and silver carp, and also facilitate the identification of candidate genes related to disease resistance.

Restriction site-associated DNA sequencing for SNP discovery and high-density genetic map construction in southern catfish (Silurus meridionalis)

Single-nucleotide polymorphism (SNP) markers and high-density genetic maps are important resources for marker-assisted selection, mapping of quantitative trait loci (QTLs) and genome structure analysis. Although linkage maps in certain catfish species have been obtained, high-density maps remain unavailable in the economically important southern catfish (Silurus meridionalis). Recently developed restriction site-associated DNA (RAD) markers have proved to be a promising tool for SNP detection and genetic map construction. The objective of the present study was to construct a high-density linkage map using SNPs generated by next-generation RAD sequencing in S. meridionalis for future genetic and genomic studies. An F1 population of 100 individuals was obtained by intraspecific crossing of two wild heterozygous individuals. In total, 77 634 putative high-quality bi-allelic SNPs between the parents were discovered by mapping the parents' paired-end RAD reads onto the reference contigs from both parents, of which 54.7% were transitions and 45.3% were transversions (transition/transversion ratio of 1.2). Finally, 26 714 high-quality RAD markers were grouped into 29 linkage groups by using de novo clustering methods (Stacks). Among these markers, 4514 were linked to the female genetic map, 23 718 to the male map and 6715 effective loci were linked to the integrated map spanning 5918.31 centimorgans (cM), with an average marker interval of 0.89 cM. High-resolution genetic maps are a useful tool for both marker-assisted breeding and various genome investigations in catfish, such as sequence assembly, gene localization, QTL detection and genome structure comparison. Hence, such a high-density linkage map will serve as a valuable resource for comparative genomics and fine-scale QTL mapping in catfish species.

A High-Density Genetic Linkage Map and QTL Fine Mapping for Body Weight in Crucian Carp (Carassius auratus) Using 2b-RAD Sequencing

A high-resolution genetic linkage map is essential for a wide range of genetics and genomics studies such as comparative genomics analysis and QTL fine mapping. Crucian carp (Carassius auratus) is widely distributed in Eurasia, and is an important aquaculture fish worldwide. In this study, a high-density genetic linkage map was constructed for crucian carp using 2b-RAD technology. The consensus map contains 8487 SNP markers, assigning to 50 linkage groups (LGs) and spanning 3762.88 cM, with an average marker interval of 0.44 cM and genome coverage of 98.8%. The female map had 4410 SNPs, and spanned 3500.42 cM (0.79 cM/marker), while the male map had 4625 SNPs and spanned 3346.33 cM (0.72 cM/marker). The average recombination ratio of female to male was 2.13:1, and significant male-biased recombination suppressions were observed in LG47 and LG49. Comparative genomics analysis revealed a clear 2:1 syntenic relationship between crucian carp LGs and chromosomes of zebrafish and grass carp, and a 1:1 correspondence, but extensive chromosomal rearrangement, between crucian carp and common carp, providing evidence that crucian carp has experienced a fourth round of whole genome duplication (4R-WGD). Eight chromosome-wide QTL for body weight at 2 months after hatch were detected on five LGs, explaining 10.1-13.2% of the phenotypic variations. Potential candidate growth-related genes, such as an EGF-like domain and TGF-β, were identified within the QTL intervals. This high-density genetic map and QTL analysis supplies a basis for genome evolutionary studies in cyprinid fishes, genome assembly, and QTL fine mapping for complex traits in crucian carp.

The genetic map of goldfish (Carassius auratus) provided insights to the divergent genome evolutions in the Cyprinidae family

A high-density linkage map of goldfish (Carassius auratus) was constructed using RNA-sequencing. This map consists of 50 linkage groups with 8,521 SNP markers and an average resolution of 0.62 cM. Approximately 84% of markers are in protein-coding genes orthologous to zebrafish proteins. We performed comparative genome analysis between zebrafish and medaka, common carp, grass carp, and goldfish to study the genome evolution events in the Cyprinidae family. The comparison revealed large synteny blocks among Cyprinidae fish and we hypothesized that the Cyprinidae ancestor undergone many inter-chromosome rearrangements after speciation from teleost ancestor. The study also showed that goldfish genome had one more round of whole genome duplication (WGD) than zebrafish. Our results illustrated that most goldfish markers were orthologous to genes in common carp, which had four rounds of WGD. Growth-related regions and genes were identified by QTL analysis and association study. Function annotations of the associated genes suggested that they might regulate development and growth in goldfish. This first genetic map enables us to study the goldfish genome evolution and provides an important resource for selective breeding of goldfish.

A high-density genetic map and growth related QTL mapping in bighead carp (Hypophthalmichthys nobilis)

Growth related traits in fish are controlled by quantitative trait loci (QTL), but no QTL for growth have been detected in bighead carp (Hypophthalmichthys nobilis) due to the lack of high-density genetic map. In this study, an ultra-high density genetic map was constructed with 3,121 SNP markers by sequencing 117 individuals in a F₁ family using 2b-RAD technology. The total length of the map was 2341.27 cM, with an average marker interval of 0.75 cM. A high level of genomic synteny between our map and zebrafish was detected. Based on this genetic map, one genome-wide significant and 37 suggestive QTL for five growth-related traits were identified in 6 linkage groups (i.e. LG3, LG11, LG15, LG18, LG19, LG22). The phenotypic variance explained (PVE) by these QTL varied from 15.4% to 38.2%. Marker within the significant QTL region was surrounded by CRP1 and CRP2, which played an important role in muscle cell division. These high-density map and QTL information provided a solid base for QTL fine mapping and comparative genomics in bighead carp.

Comparative mapping for bighead carp (Aristichthys nobilis) against model and non-model fishes provides insights into the genomic evolution of cyprinids

Comparative mapping provides an efficient method to connect genomes of non-model and model fishes. In this study, we used flanking sequences of the 659 microsatellites on a genetic map of bighead carp (Aristichthys nobilis) to comprehensively study syntenic relationships between bighead carp and nine model and non-model fishes. Of the five model and two food fishes with whole genome data, Cyprinus carpio showed the highest rate of positive BLAST hits (95.3 %) with bighead carp map, followed by Danio rerio (70.9 %), Oreochromis niloticus (21.7 %), Tetraodon nigroviridis (6.4 %), Gasterosteus aculeatus (5.2 %), Oryzias latipes (4.7 %) and Fugu rubripes (3.5 %). Chromosomal syntenic analyses showed that inversion was the basic chromosomal rearrangement during genomic evolution of cyprinids, and the extent of inversions and translocations was found to be positively correlated with evolutionary relationships among fishes studied. Among the five investigated cyprinids, linkage groups (LGs) of bighead carp, Hypophthalmichthys molitrix and Ctenopharyngodon idella exhibited a one-to-one relationship. Besides, LG 9 of bighead carp and homologous LGs of silver carp and grass carp all corresponded to the chromosomes 10 and 22 of zebrafish, suggesting that chromosomal fission may have occurred in the ancestor of zebrafish. On the other hand, LGs of bighead carp and common carp showed an approximate one-to-two relationship with extensive translocations, confirming the occurrence of a 4th whole genome duplication in common carp. This study provides insights into the understanding of genome evolution among cyprinids and would aid in transferring positional and functional information of genes from model fish like zebrafish to non-model fish like bighead carp.

Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish

The traits of cultured fish must continually be genetically improved to supply high-quality animal protein for human consumption. Economically important fish traits are controlled by multiple gene quantitative trait loci (QTL), most of which have minor effects, but a few genes may have major effects useful for molecular breeding. In this review, we chose relevant studies on some of the most intensively cultured fish and concisely summarize progress on identifying and verifying QTLs for such traits as growth, disease and stress resistance and sex in recent decades. The potential applications of these major-effect genes and their associated markers in marker-assisted selection and molecular breeding, as well as future research directions are also discussed. These genetic and genomic analyses will be valuable for elucidating the mechanisms modulating economically important traits and to establish more effective molecular breeding techniques in fish.

Centromere localization for Bighead Carp (Aristichthys nobilis) through half-tetrad analysis in diploid gynogenetic families

Gene-centromere (G-C) mapping provides insights into structural and behavioural properties of chromosomes. In this study, G-C mapping using microsatellite markers and meiogynogenetic (meiotic gynogenetic) families were performed in bighead carp (Aristichthys nobilis, 2N = 48), which belongs to Cyprinidae. A total of 218 microsatellites were selected across 24 linkage groups (LGs) of a recently well-defined genetic linkage map for bighead carp, with 151 being heterozygous in at least one of six dams in diploid meiogynogenetic families. After tests for Mendelian segregation in two diploid control families, 103 microsatellites were used for G-C distance calculation in 383 gynogens. The second division segregation frequency (y) was computed through half-tetrad analyses, and the values ranged from 0 to 0.97 (mean 0.40). High G-C recombination frequencies (over 0.667) were observed in 18 (17.5%) of the loci examined, which revealed a low level of chiasma interferences compared with other fishes studied previously. Distribution of G-C distances across LGs ranged from 0 cM to 48.5 cM (mean 20 cM) under the assumption of complete interference. All 24 centromeres were localized according to their closest-related microsatellites at 95% confident intervals. The average distance between centromeres and their closest-linked markers was 6.1 cM with 15 out of 24 LGs having a distance below 5 cM. Based on the centromere positions in this study, we proposed a formula of 24 m/sm+24 t/st chromosomes with 92 arms for bighead carp, which was mostly in accordance with a previously reported karyotype for bighead carp (24 m/sm+24 st). These results of centromere localization provide a basic framework and important resources for genetics and comparative genomics studies in bighead carp and its closely-related cyprinid species.