Construction of a high-density linkage map and fine mapping of QTLs for growth and gonad related traits in blunt snout bream

Linkage mapping and QTL analysis of growth traits in Rhopilema esculentum

R. esculentum is a popular seafood in Asian countries and an economic marine fishery resource in China. However, the genetic linkage map and growth-related molecular markers are still lacking, hindering marker assisted selection (MAS) for genetic improvement of R. esculentum. Therefore, we firstly used 2b-restriction site-associated DNA (2b-RAD) method to sequence 152 R. esculentum specimens and obtained 9100 single nucleotide polymorphism (SNP) markers. A 1456.34 cM linkage map was constructed using 2508 SNP markers with an average interval of 0.58 cM. Then, six quantitative trait loci (QTLs) for umbrella diameter and body weight were detected by QTL analysis based on the new linkage map. The six QTLs are located on four linkage groups (LGs), LG4, LG13, LG14 and LG15, explaining 9.4% to 13.4% of the phenotypic variation. Finally, 27 candidate genes in QTLs regions of LG 14 and 15 were found associated with growth and one gene named RE13670 (sushi, von Willebrand factor type A, EGF and pentraxin domain-containing protein 1-like) may play an important role in controlling the growth of R. esculentum. This study provides valuable information for investigating the growth mechanism and MAS breeding in R. esculentum.

Quantitative trait loci for growth-related traits in Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing

This study aimed to identify quantitative trait loci (QTLs) for growth-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers in Japanese quail. A QTL mapping population of 277 F₂ birds was obtained from an intercross between a male of a large-sized strain and three females of a normal-sized strain. Body weight (BW) was measured weekly from hatching to 16 weeks of age. Non-linear regression growth models of Weibull, Logistic, Gompertz, Richards, and Brody were analyzed, and growth curve parameters of Richards was selected as the best model to describe the quail growth curve of the F₂ birds. Restriction-site associated DNA sequencing developed 125 SNP markers that were informative between their parental strains. The SNP markers were distributed on 16 linkage groups that spanned 795.9 centiMorgan (cM) with an average marker interval of 7.3 cM. QTL analysis of phenotypic traits revealed four main-effect QTLs. Detected QTLs were located on chromosomes 1 and 3 and were associated with BW from 4 to 16 weeks of age and asymptotic weight of Richards model at genome-wide significant at 1% or 5% level. No QTL was detected for BW from 0 to 3 weeks of age. This is the first report identified QTLs for asymptotic weight of the Richards parameter in Japanese quail. These results highlight that the combination of QTL studies and the RAD-seq method will aid future breeding programs identify genes underlying the QTL and the application of marker-assisted selection in the poultry industry, particularly the Japanese quail.

Mapping of Quantitative Trait Loci Controlling Egg-Quality and -Production Traits in Japanese Quail (Coturnix japonica) Using Restriction-Site Associated DNA Sequencing

This research was conducted to identify quantitative trait loci (QTL) associated with egg-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers using restriction-site associated DNA sequencing (RAD-seq) in Japanese quail. A total of 138 F₂ females were produced by full-sib mating of F₁ birds derived from an intercross between a male of the large-sized strain with three females of the normal-sized strain. Eggs were investigated at two different stages: the beginning stage of egg-laying and at 12 weeks of age (second stage). Five eggs were analyzed for egg weight, lengths of the long and short axes, egg shell strength and weight, yolk weight and diameter, albumen weight, egg equator thickness, and yolk color (L*, a*, and b* values) at each stage. Moreover, the age at first egg, the cumulative number of eggs laid, and egg production rate were recorded. RAD-seq developed 118 SNP markers and mapped them to 13 linkage groups using the Map Manager QTX b20 software. Markers were spanned on 776.1 cM with an average spacing of 7.4 cM. Nine QTL were identified on chromosomes 2, 4, 6, 10, 12, and Z using the simple interval mapping method in the R/qtl package. The QTL detected affected 10 egg traits of egg weight, lengths of the long and short axes of egg, egg shell strength, yolk diameter and weight, albumen weight, and egg shell weight at the beginning stage, yellowness of the yolk color at the second stage, and age at first egg. This is the first report to perform a quail QTL analysis of egg-related traits using RAD-seq. These results highlight the effectiveness of RAD-seq associated with targeted QTL and the application of marker-assisted selection in the poultry industry, particularly in the Japanese quail.

A high-density genetic map construction and sex-related loci identification in Chinese Giant salamander

BACKGROUND

The Chinese giant salamander Andrias davidianus is an important amphibian species in China because of its increasing economic value, protection status and special evolutionary position from aquatic to terrestrial animal. Its large genome presents challenges to genetic research. Genetic linkage mapping is an important tool for genome assembly and determination of phenotype-related loci.

RESULTS

In this study, we constructed a high-density genetic linkage map using ddRAD sequencing technology to obtain SNP genotyping data of members from an full-sib family which sex had been determined. A total of 10,896 markers were grouped and oriented into 30 linkage groups, representing 30 chromosomes of A. davidianus. The genetic length of LGs ranged from 17.61 cM (LG30) to 280.81 cM (LG1), with a mean inter-locus distance ranging from 0.11(LG3) to 0.48 cM (LG26). The total genetic map length was 2643.10 cM with an average inter-locus distance of 0.24 cM. Three sex-related loci and four sex-related markers were found on LG6 and LG23, respectively.

CONCLUSION

We constructed the first High-density genetic linkage map and identified three sex-related loci in the Chinese giant salamander. Current results are expected to be a useful tool for future genomic studies aiming at the marker-assisted breeding of the species.

Whole genome sequencing of silver carp (Hypophthalmichthys molitrix) and bighead carp (Hypophthalmichthys nobilis) provide novel insights into their evolution and speciation

The edible silver carp (Hypophthalmichthys molitrix) and bighead carp (H. nobilis), which are two of the "Four Domesticated Fish" of China, are cultivated intensively worldwide. Here, we constructed 837- and 845-Mb draft genome assemblies for the silver carp and the bighead carp, respectively, including 24,571 and 24,229 annotated protein-coding genes. Genetic maps, anchoring 71.7% and 83.8% of all scaffolds, were obtained for the silver and bighead carp, respectively. Phylogenetic analysis showed that the bighead carp formed a clade with the silver carp, with an estimated divergence time of 3.6 million years ago; the time of divergence between the silver carp and zebrafish was 50.7 million years ago. An East Asian cyprinid genome-specific chromosome fusion took place ~9.2 million years after this clade diverged from the clade containing the common carp and Sinocyclocheilus. KEGG and GO analyses indicated that the expanded gene families in the silver and bighead carp were associated with diseases, the immune system and environmental adaptations. Genomic regions differentiating the silver and bighead carp populations were detected based on the whole-genome sequences of 42 individuals. Genes associated with the divergent regions were associated with reproductive system development and the development of primary female sexual characteristics. Thus, our results provided a novel systematic genomic analysis of the East Asian cyprinids, as well as the evolution and speciation of the silver carp and bighead carp.

A high-density SNP-based genetic map and several economic traits-related loci in Pelteobagrus vachelli

BACKGROUND

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection in aquaculture species. Pelteobagrus vachelli is a very popular commercial species in Asia. However, some specific characters hindered achievement of the traditional selective breeding based on phenotypes, such as lack of large-scale genomic resource and short of markers tightly associated with growth, sex determination and hypoxia tolerance related traits.

RESULTS

By making use of 5059 ddRAD markers in P. vachelli, a high-resolution genetic linkage map was successfully constructed. The map' length was 4047.01 cM by using an interval of 0.11 cm, which is an average marker standard. Comparative genome mapping revealed that a high proportion (83.2%) of markers with a one-to-one correspondence were observed between P. vachelli and P. fulvidraco. Based on the genetic map, 8 significant genome-wide QTLs for 4 weight, 1 body proportion, 2 sex determination, and 1 hypoxia tolerance related traits were detected on 4 LGs. Some SNPs from these significant genome-wide QTLs were observably associated with these phenotypic traits in other individuals by Kompetitive Allele Specific PCR. In addition, two candidate genes for weight, Sipa1 and HSD11B2, were differentially expressed between fast-, medium- and slow-growing P. vachelli. Sema7a, associated with hypoxia tolerance, was induced after hypoxia exposure and reoxygenation.

CONCLUSIONS

We mapped a set of suggestive and significant QTLs as well as candidate genes for 12 growth, 1 sex determination and 1 hypoxia tolerance related traits based on a high-density genetic linkage map by making use of SNP markers for P. fulvidraco. Our results have offered a valuable method about the much more efficient production of all-male, fast growth and hypoxia tolerance P. vachelli for the aquaculture industry.

A first genetic linage map construction and QTL mapping for growth traits in Larimichthys polyactis

Larimichthys polyactis is a commercially important marine fish species in Eastern Asia, yet very few genetic resources exist. In particular, genetic linkage maps are critical tools for genetic breeding. In this study, we generated a high resolution linkage map from a family of 110 individuals and their parents by resequencing the individuals. 3,802 effective SNPs were mapped to 24 linkage groups (LGs). The map spanned 2,567.39 cm, with an average marker interval of 0.66 cm. We used the map to conduct QTL analysis for growth traits, and found 31 markers were significantly associated with growth-related traits. Specifically, three SNPs were identified for total length, nineteen SNPs for body length, and nine SNPs for body weight. The identified SNPs could explain 15.2-22.6% of the phenotypic variation. SNPs associated with growth traits were distributed on LG6 and LG11, and candidate genes included, kif26b, bat1, gna1, gbgt1, and amfr, which may regulate growth. The linkage map and mapped QTLs would be useful for improving the quality of L. polyactis via marker-assisted selection.

ddRADseq-assisted construction of a high-density SNP genetic map and QTL fine mapping for growth-related traits in the spotted scat (Scatophagus argus)

BACKGROUND

Scatophagus argus is a popular farmed fish in several countries of Southeast Asia, including China. Although S. argus has a highly promising economic value, a significant lag of breeding research severely obstructs the sustainable development of aquaculture industry. As one of the most important economic traits, growth traits are controlled by multiple gene loci called quantitative trait loci (QTLs). It is urgently needed to launch a marker assisted selection (MAS) breeding program to improve growth and other pivotal traits. Thus a high-density genetic linkage map is necessary for the fine mapping of QTLs associated with target traits.

RESULTS

Using restriction site-associated DNA sequencing, 6196 single nucleotide polymorphism (SNP) markers were developed from a full-sib mapping population for genetic map construction. A total of 6193 SNPs were grouped into 24 linkage groups (LGs), and the total length reached 2191.65 cM with an average marker interval of 0.35 cM. Comparative genome mapping revealed 23 one-to-one and 1 one-to-two syntenic relationships between S. argus LGs and Larimichthys crocea chromosomes. Based on the high-quality linkage map, a total of 44 QTLs associated with growth-related traits were identified on 11 LGs. Of which, 19 significant QTLs for body weight were detected on 9 LGs, explaining 8.8-19.6% of phenotypic variances. Within genomic regions flanking the SNP markers in QTL intervals, we predicted 15 candidate genes showing potential relationships with growth, such as Hbp1, Vgll4 and Pim3, which merit further functional exploration.

CONCLUSIONS

The first SNP genetic map with a fine resolution of 0.35 cM for S. argus has been developed, which shows a high level of syntenic relationship with L. crocea genomes. This map can provide valuable information for future genetic, genomic and evolutionary studies. The QTLs and SNP markers significantly associated with growth-related traits will act as useful tools in gene mapping, map-based cloning and MAS breeding to speed up the genetic improvement in important traits of S. argus. The interesting candidate genes are promising for further investigations and have the potential to provide deeper insights into growth regulation in the future.

Fine Mapping of the High-pH Tolerance and Growth Trait-Related Quantitative Trait Loci (QTLs) and Identification of the Candidate Genes in Pacific White Shrimp (Litopenaeus vannamei)

High-pH tolerance and growth are important traits for the shrimp culture industry in areas with saline-alkali water. In the present study, an F1 full-sib family of Pacific white shrimp (Litopenaeus vannamei) was generated with a new "semidirectional cross" method, and double-digest restriction site-associated DNA sequencing (ddRAD-Seq) technology was applied to genotype the 2 parents and 148 progenies. A total of 3567 high-quality markers were constructed for the genetic linkage map, and the total map length was 4161.555 centimorgans (cM), showing 48 linkage groups (LGs) with an average interlocus length of 1.167 cM. With a constrained logarithm of odds (LOD) score ≥ 2.50, 12 high-pH tolerance and 2 growth (body weight) QTLs were located. L. vannamei genomic scaffolds were used to assist with the detection of 21 stress- and 5 growth-related scaffold genes. According to the high-pH transcriptome data of our previous study, 6 candidate high-pH response genes were discovered, and 5 of these 6 genes were consistently expressed with the high-pH transcriptome data, validating the locations of the high-pH tolerance trait-related QTLs in this study. This paper is the first report of fine-mapping high-pH tolerance and growth (body weight) trait QTLs in one L. vannamei genetic map. Our results will further benefit marker-assisted selection work and might be useful for promoting genomic research on the shrimp L. vannamei.

Quantitative Trait Loci (QTL) Mapping

Quantitative trait loci (QTL) are genetic regions that influence phenotypic variation of a complex trait, often through genetic interactions with each other and the environment. These are commonly identified through a statistical genetic analysis known as QTL mapping. Here, I present a step-by-step, practical approach to QTL mapping along with a sample data file. I focus on methods commonly used and discoveries that have been made in fishes, and utilize a multiple QTL mapping (MQM) approach in the free software package R/qtl.

An SNP-Based Genetic Map and QTL Mapping for Growth Traits in the Red-Spotted Grouper (Epinephelus akaara)

The red-spotted grouper (Epinephelus akaara) is one of the most commercially important aquatic species in China. However, its seedstock has low larval survival rates, and its stability is confronted with the danger of overexploitation. In this study, a high-density genetic map was constructed using 3435 single nucleotide polymorphisms (SNPs) from 142 first generation (F₁) full-sib offspring and two parents of a red-spotted grouper population. The total genetic length of the map was 2300.12 cM with an average intermarker distance of 0.67 cM. Seventeen genome-wide significant quantitative trait loci (QTLs) for growth-related traits were detected on 24 linkage groups, including 5 QTLs for full length, 7 QTLs for body length, and 5 QTLs for body weight. The contribution values of explained phenotypic variance ranged from 10.7% to 12.9%. Moreover, 13 potential candidate genes for growth-related traits were identified. Collectively, these findings will be useful for conducting marker-assisted selection of the red-spotted grouper in future studies.

A High-Density Genetic Linkage Map and QTL Mapping for Sex in Black Tiger Shrimp (Penaeus monodon)

The black tiger shrimp, Penaeus monodon, is important in both fishery and aquaculture and is the second-most widely cultured shrimp species in the world. However, the current strains cannot meet the market needs in various cultural environments, and the genome resources for P. monodon are still lacking. Restriction-site associated DNA sequencing (RADseq) has been widely used in genetic linkage map construction and in quantitative trait loci (QTL) mapping. We constructed a high-density genetic linkage map with RADseq in a full-sib family. This map contained 6524 single nucleotide polymorphisms (SNPs) and 2208 unique loci. The total length was 3275.4 cM, and the genetic distance was estimated to be 1.1 Mb/cM. The sex trait is a dichotomous phenotype, and the same interval was detected as a QTL using QTL mapping and genome-wide association analysis. The most significant locus explained 77.4% of the phenotype variance. The sex locus was speculated to be the same in this species based on the sequence alignments in Mozambique, India, and Hawaii populations. The constructed genetic linkage map provided a valuable resource for QTL mapping, genome assembly, and genome comparison for shrimp. The demonstrated common sex locus is a step closer to locating the underlying gene.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.

Construction of a High-Density Linkage Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Channel Catfish (Ictalurus punctatus)

A high-density genetic linkage map is of particular importance in the fine mapping for important economic traits and whole genome assembly in aquaculture species. The channel catfish (Ictalurus punctatus), a species native to North America, is one of the most important commercial freshwater fish in the world. Outside of the United States, China has become the major producer and consumer of channel catfish after experiencing rapid development in the past three decades. In this study, based on restriction site associated DNA sequencing (RAD-seq), a high-density genetic linkage map of channel catfish was constructed by using single nucleotide polymorphisms (SNPs) in a F₁ family composed of 156 offspring and their two parental individuals. A total of 4,768 SNPs were assigned to 29 linkage groups (LGs), and the length of the linkage map reached 2,480.25 centiMorgans (cM) with an average distance of 0.55 cM between loci. Based on this genetic linkage map, 223 genomic scaffolds were anchored to the 29 LGs of channel catfish, and a total length of 704.66 Mb was assembled. Quantitative trait locus (QTL) mapping and genome-wide association analysis identified 10 QTLs of sex-related and six QTLs of growth-related traits at LG17 and LG28, respectively. Candidate genes associated with sex dimorphism, including spata2, spata5, sf3, zbtb38, and fox, were identified within QTL intervals on the LG17. A sex-linked marker with simple sequence repeats (SSR) in zbtb38 gene of the LG17 was validated for practical verification of sex in the channel catfish. Thus, the LG17 was considered as a sex-related LG. Potential growth-related genes were also identified, including important regulators such as megf9, npffr1, and gas1. In a word, we constructed the high-density genetic linkage map and developed the sex-linked marker in channel catfish, which are important genetic resources for future marker-assisted selection (MAS) of this economically important teleost.

Dissecting Key Adaptation Traits in the Polyploid Perennial Medicago sativa Using GBS-SNP Mapping

Understanding key adaptation traits is crucial to developing new cultivars with broad adaptations. The main objective of this research is to understand the genetic basis of winter hardiness (WH) and fall dormancy (FD) in alfalfa and the association between the two traits. QTL analysis was conducted in a pseudo-testcross F1 population developed from two cultivars contrasting in FD (3010 with FD = 2 and CW 1010 with FD = 10). The mapping population was evaluated in three replications at two locations (Watkinsville and Blairsville, GA). FD levels showed low to moderate correlations with WH (0.22-0.57). Assessing dormancy in winter is more reliable than in the fall in southern regions with warm winters. The mapping population was genotyped using Genotyping-by-sequencing (GBS). Single dose allele SNPs (SDA) were used for constructing linkage maps. The parental map (CW 1010) consisted of 32 linkage groups spanning 2127.5 cM with 1377 markers and an average marker density of 1.5 cM/SNP. The maternal map (3010) had 32 linkage groups spanning 2788.4 cM with 1837 SDA SNPs with an average marker density of 1.5 cM/SNP. Forty-five significant (P < 0.05) QTLs for FD and 35 QTLs for WH were detected on both male and female linkage maps. More than 75% (22/28) of the dormancy QTL detected from the 3010 parent did not share genomic regions with WH QTLs and more than 70% (12/17) dormancy QTLs detected from CW 1010 parent were localized in different genomic regions than WH QTLs. These results suggest that the two traits have independent inheritance and therefore can be improved separately in breeding programs.