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

FGFs function in regulating myoblasts differentiation in spotted sea bass (Lateolabrax maculatus)

Fibroblast growth factors (FGFs) are a family of structurally related peptides that regulate processes such as cell proliferation, differentiation, and damage repair. In our previous study, fibroblast growth factor receptor 4 (fgfr4) was detected in the most significant quantitative trait loci (QTL), when identified of QTLs and genetic markers for growth-related traits in spotted sea bass. However, knowledge of the function of fgfr4 is lacking, even the legends to activate the receptor is unknown in fish. To remedy this problem, in the present study, a total of 33 fgfs were identified from the genomic and transcriptomic databases of spotted sea bass, of which 10 were expressed in the myoblasts. According to the expression pattern during myoblasts proliferation and differentiation, fgf6a, fgf6b and fgf18 were selected for further prokaryotic expression and purification. The recombinant proteins FGF6a, FGF6b and FGF18 were found to inhibit myoblast differentiation. Overall, our results provide a theoretical basis for the molecular mechanisms of growth regulation in economic fish such as spotted sea bass.

Investigating the role of genetic variation in vgll3 and six6 in the domestication of gilthead seabream (Sparus aurata Linnaeus) and European seabass (Dicentrarchus labrax Linnaeus)

Gene function conservation is crucial in molecular ecology, especially for key traits like growth and maturation in teleost fish. The vgll3 and six6 genes are known to influence age-at-maturity in Atlantic salmon, but their impact on other fish species is poorly understood. Here, we investigated the association of vgll3 and six6 in the domestication of gilthead seabream and European seabass, both undergoing selective breeding for growth-related traits in the Mediterranean. We analysed two different sets of samples using two different genotyping approaches. The first dataset comprised farmed and wild populations from Greece, genotyped for SNPs within the two genes ('gene-level genotyping'). The second dataset examined 300-600 k SNPs located in the chromosomes of the two genes, derived from a meta-analysis of a Pool-Seq experiment involving farmed and wild populations distributed widely across the Mediterranean ('chromosome-level genotyping'). The gene-level analysis revealed a statistically significant allele frequency differences between farmed and wild populations on both genes in each species. This finding was partially supported by the chromosome-level analysis, identifying highly differentiated regions may be involved in the domestication process at varying distances from the candidate genes. Noteworthy genomic features were found, such as a CpG island in gilthead seabream and novel candidate genes in European seabass, warranting further investigation. These findings support a putative role of vgll3 and six6 in the maturation and growth of gilthead seabream and European seabass, emphasizing the need for further research on their conserved function.

Genome-Wide RADseq Reveals Genetic Differentiation of Wild and Cultured Populations of Large Yellow Croaker

Larimichthys crocea (also known as the large yellow croaker) is one of the most economically important marine fishes in China, and research on the ecology and genetics of this species is of immense significance. In this study, we performed restriction site-associated DNA sequencing (RAD-seq) of 54 individuals collected from four sites in China to analyze the genetic structure and diversity of large yellow croaker at the genome level. It revealed that the large yellow croaker populations in the Ningde and Zhoushan coastal waters can be clearly distinguished. Different genetic diversity indices were used to analyze the genetic diversity of the large yellow croaker, which showed that there was a differentiation trend between the wild and farmed populations in Ningde. Moreover, we identified genetically differentiated genomic regions between the populations. GO gene enrichment analysis identified genes that are related to fatty acid metabolism and growth. These findings enhance our understanding of genetic differentiation and adaptation to different living environments, providing a theoretical basis for the preservation and restoration of the genetic resources of the large yellow croaker.

Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers

This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.

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.

Genome-Wide Association Study of Growth Traits in a Four-Way Crossbred Pig Population

Growth traits are crucial economic traits in the commercial pig industry and have a substantial impact on pig production. However, the genetic mechanism of growth traits is not very clear. In this study, we performed a genome-wide association study (GWAS) based on the specific-locus amplified fragment sequencing (SLAF-seq) to analyze ten growth traits on 223 four-way intercross pigs. A total of 227,921 highly consistent single nucleotide polymorphisms (SNPs) uniformly dispersed throughout the entire genome were used to conduct GWAS. A total of 53 SNPs were identified for ten growth traits using the mixed linear model (MLM), of which 18 SNPs were located in previously reported quantitative trait loci (QTL) regions. Two novel QTLs on SSC4 and SSC7 were related to average daily gain from 30 to 60 kg (ADG30–60) and body length (BL), respectively. Furthermore, 13 candidate genes (ATP5O, GHRHR, TRIM55, EIF2AK1, PLEKHA1, BRAP, COL11A2, HMGA1, NHLRC1, SGSM1, NFATC2, MAML1, and PSD3) were found to be associated with growth traits in pigs. The GWAS findings will enhance our comprehension of the genetic architecture of growth traits. We suggested that these detected SNPs and corresponding candidate genes might provide a biological foundation for improving the growth and production performance of pigs in swine breeding.

Differential expression of MSTN, IGF2BP1, and FABP2 across different embryonic ages and sexes in white Muscovy ducks

To explore the effects of growth-related genes in both sexes and at different growth and development stages, male and female white Muscovy ducks at embryonic day E13, E17, E21, E25 and E29 were assessed in this study. RT-qPCR was used to determine the mRNA transcription levels of selected growth-related genes in the leg muscles of Muscovy ducks of both sexes and at different growth and developmental stages. MSTN, IGF2BP1 and FABP2 mRNAs were expressed in the leg muscles of male and female Muscovy ducks, but with different expression patterns. The MSTN and IGF2BP1 mRNA expression patterns were wavelike. MSTN mRNA expression was elevated at E13, increased at E17, decreased rapidly to the lowest level at E21, increased again at E25, and then decreased. IGF2BP1 mRNA expression was elevated at E13, increased at E17, decreased rapidly at E21, decreased rapidly to the lowest level at E25, and increased at E29. The expression trend of FABP2 mRNA was approximately "⊥" shape; the expression was the lowest at E13, increased slowly from E17 to E25, and increased extremely significantly at E29. In addition, the expression of MSTN in male Muscovy ducks was significantly higher than that in female ducks at E25 (P < 0.05). The expression of IGF2BP1 in male Muscovy ducks was extremely significantly higher than that in female ducks at E17 (P < 0.01). However, the expression of FABP2 in female Muscovy ducks was extremely significantly higher than that in male Muscovy ducks at E21 and E29 (P < 0.01). In conclusion, the mRNA expression of MSTN, IGF2BP1 and FABP2 in white Muscovy ducks is gestational age specific and sex specific. The differential gene expression patterns observed in this study provide a basis for understanding the physiological changes in white Muscovy ducks at different embryonic ages and in both sexes, supplementing the existing research on duck embryo muscle development. In addition, the findings provide a new framework for further discussion of poultry breeding.

Hybridization alters the shape of the genotypic fitness landscape, increasing access to novel fitness peaks during adaptive radiation

Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 single nucleotide polymorphisms (SNPs) that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation.

Identification of crucial factors involved in Cynoglossus semilaevis sexual size dimorphism by GWAS and demonstration of zbed1 regulatory network by DAP-seq

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species.

Functional Differentiation of BMP7 Genes in Zebrafish: bmp7a for Dorsal-Ventral Pattern and bmp7b for Melanin Synthesis and Eye Development

Bone morphogenetic protein 7 (BMP7) belongs to the transforming growth factor β (TGF-β) family, which not only induces cartilage and bone formation, but also regulates eye development and melanoma tumorigenesis in mammals. In teleosts, BMP7 differentiates into two subtypes, bmp7a and bmp7b, which have clearly differentiated structures. To fully understand the functional differentiation of bmp7a and bmp7b in fish species, we successfully constructed bmp7a and bmp7b gene deletion mutants in zebrafish using CRISPR/Cas9-mediated gene editing technology. Our results showed that bmp7a mutation caused abnormal development of the embryo’s dorsal-ventral pattern that led to death; bmp7b mutation induced growth inhibition and increased melanin production in the skin and eye of mutants. Histological analysis revealed that melanin in the retina of the eyes in bmp7b mutants increased, and behavioral observation showed that the vision and sensitivity to food of the mutants were reduced. Transcriptome analysis of the skin and eye tissues showed that the expression changes of wnt7ba and gna14 in bmp7b mutants might promote the increase of melanin. Additionally, the eye transcriptome analysis indicated that changes in the structure of the eyes in bmp7b mutants led to defects in phototransduction, and seven DEGs (rgs9a, rgs9b, rcvrn2, guca1d, grk1b, opn1mw4, and gc2) were identified as key candidate genes that affected the photonic response of the eyes. The study revealed the functional differentiation of bmp7a and bmp7b in teleosts and the first report about the inhibitory effect of bmp7b on melanogenesis may provide useful information for the future research on human melanoma-related diseases.

Genomic prediction of growth in a commercially, recreationally, and culturally important marine resource, the Australian snapper (Chrysophrys auratus)

Growth is one of the most important traits of an organism. For exploited species, this trait has ecological and evolutionary consequences as well as economical and conservation significance. Rapid changes in growth rate associated with anthropogenic stressors have been reported for several marine fishes, but little is known about the genetic basis of growth traits in teleosts. We used reduced genome representation data and genome-wide association approaches to identify growth-related genetic variation in the commercially, recreationally, and culturally important Australian snapper (Chrysophrys auratus, Sparidae). Based on 17,490 high-quality single-nucleotide polymorphisms and 363 individuals representing extreme growth phenotypes from 15,000 fish of the same age and reared under identical conditions in a sea pen, we identified 100 unique candidates that were annotated to 51 proteins. We documented a complex polygenic nature of growth in the species that included several loci with small effects and a few loci with larger effects. Overall heritability was high (75.7%), reflected in the high accuracy of the genomic prediction for the phenotype (small vs large). Although the single-nucleotide polymorphisms were distributed across the genome, most candidates (60%) clustered on chromosome 16, which also explains the largest proportion of heritability (16.4%). This study demonstrates that reduced genome representation single-nucleotide polymorphisms and the right bioinformatic tools provide a cost-efficient approach to identify growth-related loci and to describe genomic architectures of complex quantitative traits. Our results help to inform captive aquaculture breeding programs and are of relevance to monitor growth-related evolutionary shifts in wild populations in response to anthropogenic pressures.

Construction of Genetic Linkage Maps From a Hybrid Family of Large Yellow Croaker (Larimichthys crocea)

Consensus and sex-specific genetic linkage maps for large yellow croaker (Larimichthys crocea) were constructed using samples from an F₁ family produced by crossing a Daiqu female and a Mindong male. A total of 20,147 single nucleotide polymorphisms (SNPs) by restriction site associated DNA sequencing were assigned to 24 linkage groups (LGs). The total length of the consensus map was 1757.4 centimorgan (cM) with an average marker interval of 0.09 cM. The total length of female and male linkage map was 1533.1 cM and 1279.2 cM, respectively. The average female-to-male map length ratio was 1.2 ± 0.23. Collapsed markers in the genetic maps were re-ordered according to their relative positions in the ASM435267v1 genome assembly to produce integrated genetic linkage maps with 9885 SNPs distributed across the 24 LGs. The recombination pattern of most LGs showed sigmoidal patterns of recombination, with higher recombination in the middle and suppressed recombination at both ends, which corresponds with the presence of sub-telocentric and acrocentric chromosomes in the species. The average recombination rate in the integrated female and male maps was respectively 3.55 cM/Mb and 3.05 cM/Mb. In most LGs, higher recombination rates were found in the integrated female map, compared to the male map, except in LG12, LG16, LG21, LG22, and LG24. Recombination rate profiles within each LG differed between the male and the female, with distinct regions indicating potential recombination hotspots. Separate quantitative trait loci (QTL) and association analyses for growth related traits in 6 months fish were performed, however, no significant QTL was detected. The study indicates that there may be genetic differences between the two strains, which may have implications for the application of DNA-information in the further breeding schemes.

Whole-genome resequencing of large yellow croaker (Larimichthys crocea) reveals the population structure and signatures of environmental adaptation

Large yellow croaker is an economically important fish in China and East Asia. Despite its economic importance, genome-wide adaptions of domesticated large yellow croaker are largely unknown. Here, we performed whole-genome resequencing of 198 individuals of large yellow croaker obtained in the sea or from farmers in Zhoushan or Ningde. Population genomics analyses revealed the genetic population structure of our samples, reflecting the living environment. Each effective population size is estimated to be declining over generations. Moreover, we identified genetically differentiated genomic regions between the sea-captured population in the Zhoushan Sea area and that of the Ningde Sea area or between the sea-captured population and the farmed population in either area. Gene ontology analyses revealed the gene groups under selective sweep for the adaptation to the domesticated environment. All these results suggest that individuals of the large yellow croaker populations show genomic signatures of adaptation to different living environments.

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.

Genome-Scale Phylogenetic and Population Genetic Studies Provide Insight Into Introgression and Adaptive Evolution of Takifugu Species in East Asia

As a typical marine adaptive radiation species, most Takifugu species are widely distributed in East Asian offshore, which have diversified morphological characteristics and different ecological habits. The phylogenetic relationship and population structure of the Takifugu species was complicated because of incomplete lineage sorting, widespread hybridization and introgression. Therefore, to systematically clarify the phylogenetic relationships of Takifugu genus, explore the introgression and natural hybridization between different Takifugu species, and detect the selective signatures in the adaptive evolution of diversified traits, whole-genome resequencing was used in 122 Takifugu samples from 10 species. Phylogenetic analysis showed solid sister-group relationships between Takifugu bimaculatus and Takifugu flavidus, Takifugu oblongus, and Takifugu niphobles, Takifugu rubripes, and Takifugu obscurus, Takifugu xanthoptreus, and Takifugu ocellatus. Further admixture analysis indicated the divergence of T. obscurus population and the bidirectional gene flow between T. bimaculatus and T. flavidus. Using species-specific homozygous genetic variance sites, we detected the asymmetric introgression between T. bimaculatus and T. flavidus at China East sea and southern Taiwan Strait. By genome-scale genetic diversity scanning, we detected two copies of syt1, zar1 and tgfbr1 related to the semilunar reproduction rhythm in T. niphobles, involved in memory formation, embryo maturation and female reproduction. Furthermore, we also found lots of T. niphobles specific mutations in CDS region of circadian rhythm related genes and endocrine hormone genes. For Takifugu species, our research provides reliable genetic resources and results for the phylogeny, introgression, hybridization and adaptive evolution, and could be used as a guide for the formulation of the protection and proliferation release policies.

Genome-Wide Association Study Identifies Genomic Loci of Sex Determination and Gonadosomatic Index Traits in Large Yellow Croaker (Larimichthys crocea)

Larimichthys crocea is one of the traditional marine culture fishes in China, widely distributed in South China Sea, East Sea, and southern Yellow Sea. Sex dimorphism is evident in this species that females present a substantial growth strength than males, suggesting breeding females could obtain more economic benefits in L. crocea aquaculture industry. With the continuous expansion of aquaculture industry, both identifying sex-associated genome region and understanding the genetic basis underlying gonad differentiation and development matter to not only sex control aquaculture but also breeding industry. Thus, genome-wide association analysis (GWAS) of sex determination was conducted with a random breeding population of 905 individuals (including 463 females and 442 males) by ddRAD sequencing. For sex determination, 21 significant single nucleotide polymorphisms (SNPs) in chromosome (Chr) 22 were identified. Surrounding these SNPs, we founded 14 candidate genes, including dmrt1, dmrt3, and piwil2, fam102a, and odf2. The sex-associated region was narrowed down further to 2.4 Mb on Chr22 through F_st scanning and insertion-deletion (InDel) analysis. Besides, 3 SNPs in the supposed sex-determining region on Chr22 were identified as highly associated with gonad differentiation through GWAS on gonadosomatic index (GSI) in 350 males and 231 females. Because of the significant difference of GSI between females and males of L. crocea, GWAS on GSI of different genders was also conducted independently. Finally, we identified a SNP in Chr18 showing genome-wide significant association with male GSI (MGSI) and three genes axl, cyp2a10, and cyp2g1 involved in the gonadal development regulation process of aromatase. Overall, this study explored the genetic basis of sex determination mechanism and provided novel insights into gonad differentiation and development, offering solid genetic support for sex control breeding, marker-assisted selection, and marine resources conservation.

Development and Evaluation of a High-Throughput Single-Nucleotide Polymorphism Array for Large Yellow Croaker (Larimichthys crocea)

High-density single-nucleotide polymorphism (SNP) genotyping array is an essential tool for genetic analyses of animals and plants. Large yellow croaker (Larimichthys crocea) is one of the most commercially important marine fish species in China. Although plenty of SNPs have been identified in large yellow croaker, no high-throughput genotyping array is available. In this study, a high-throughput SNP array named NingXin-I with 600K SNPs was developed and evaluated. A set of 82 large yellow croakers were collected from different locations of China and re-sequenced. A total of 9.34M SNPs were identified by mapping sequence reads to the large yellow croaker reference genome. About 1.98M candidate SNPs were selected for further analyses by using criteria such as SNP quality score and conversion performance in the final array. Finally, 579.5K SNPs evenly distributed across the large yellow croaker genome with an average spacing of 1.19 kb were proceeded to array production. The performance of NingXin-I array was evaluated in 96 large yellow croaker individuals from five populations, and 83.38% SNPs on the array were polymorphic sites. A further test of the NingXin-I array in five closely related species in Sciaenidae identified 26.68–56.23% polymorphic SNP rate across species. A phylogenetic tree inferred by using the genotype data generated by NingXin-I confirmed the phylogenetic distance of the species in Sciaenidae. The performance of NingXin-I in large yellow croaker and the other species in Sciaenidae suggested high accuracy and broad application. The NingXin-I array should be valuable for quantitative genetic studies, such as genome-wide association studies (GWASs), high-density linkage map construction, haplotype analysis, and genome-based selection.

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.

Genetic analysis of whole mitochondrial genome of Lateolabrax maculatus (Perciformes: Moronidae) indicates the presence of two populations along the Chinese coast

The whole mitochondrial genome ofLateolabrax maculatus(Cuvier, 1828) was used to investigate the reasons for the observed patterns of genetic differentiation among 12 populations in northern and southern China. The haplotype diversity and nucleotide diversity ofL. maculatuswere 0.998 and 0.00169, respectively. Pairwise FSTvalues between populations ranged from 0.001 to 0.429, correlating positively with geographic distance. Genetic structure analysis and haplotype network analysis indicated that these populations were split into two groups, in agreement with geographic segregation and environment. Tajima’s D values, Fu’s Fs tests and Bayesian skyline plot (BSP) indicated that a demographic expansion event may have occurred in the history ofL. maculatus. Through selection pressure analysis, we found evidence of significant negative selection at the ATP6, ND3, Cytb, COX3, COX2 and COX1 genes. In our hypotheses, this study implied that demographic events and selection of local environmental conditions, including temperature, are responsible for population divergence. These findings are a step forward toward the understanding of the genetic basis of differentiation and adaptation, as well as conservation ofL. maculatus.

The Polymorphism of LvMMD2 and Its Association with Growth Traits in Litopenaeus vannamei

The Pacific white shrimp Litopenaeus vannamei is one of the major economic aquaculture species. The growth trait is considered as the most important trait in L. vannamei aquaculture. Identification of the genetic components underlying growth-related traits in L. vannamei could be useful for the selective breeding of growth trait. Our previous work identified several growth-related SNPs by genome-wide association study (GWAS). Based on the assembled genome, we identified a new candidate gene (LvMMD2) beside the associated marker. This gene encodes the progestin and AdipoQ receptor 10 (PAQR10) protein. We further investigate the polymorphisms of LvMMD2 and their association with body weight of L. vannamei. By resequencing the coding region of LvMMD2, a total of 8 SNPs were identified, including 6 synonymous mutations and 2 nonsynonymous mutations. Association analyses based on a population of 322 individuals revealed that several SNPs located in the coding region of LvMMD2 were significantly associated with the body weight, especially the nonsynonymous mutation named as MMD_5 contributed the most association to the trait and it could explain 10.5% of phenotypic variance. In addition, several genes involved in growth and development have been identified as LvMMD2-interacting genes. These findings strongly suggested that LvMMD2 might be an important gene regulating the shrimp growth. More importantly, the MMD_5 could be a promising candidate locus for marker-assisted selection (MAS) of the body weight in L. vannamei.

First High-Density Linkage Map and QTL Fine Mapping for Growth-Related Traits of Spotted Sea bass (Lateolabrax maculatus)

Possessing powerful adaptive capacity and a pleasant taste, spotted sea bass (Lateolabrax maculatus) has a broad natural distribution and is one of the most popular mariculture fish in China. However, the genetic improvement program for this fish is still in its infancy. Growth is the most economically important trait and is controlled by quantitative trait loci (QTL); thus, the identification of QTLs and genetic markers for growth-related traits is an essential step for the establishment of marker-assisted selection (MAS) breeding programs. In this study, we report the first high-density linkage map of spotted sea bass constructed by sequencing 333 F1 generation individuals in a full-sib family using 2b-RAD technology. A total of 6883 SNP markers were anchored onto 24 linkage groups, spanning 2189.96 cM with an average marker interval of 0.33 cM. Twenty-four growth-related QTLs, including 13 QTLs for body weight and 11 QTLs for body length, were successfully detected, with phenotypic variance explained (PVE) ranging from 5.1 to 8.6%. Thirty potential candidate growth-related genes surrounding the associated SNPs were involved in cell adhesion, cell proliferation, cytoskeleton reorganization, calcium channels, and neuromodulation. Notably, the fgfr4 gene was detected in the most significant QTL; this gene plays a pivotal role in myogenesis and bone growth. The results of this study may facilitate marker-assisted selection for breeding populations and establish the foundation for further genomic and genetic studies investigating spotted sea bass.

Identification of Candidate Growth-Related SNPs and Genes Using GWAS in Brown-Marbled Grouper (Epinephelus fuscoguttatus)

Brown-marbled grouper, Epinephelus fuscoguttatus, is not only an important commercial fish species, but also an important crossbreeding parent in grouper industry. Improvement of growth traits of this species contributes to the development of grouper breeding. Currently, the development of molecular marker associated with growth of brown-marbled grouper is rare. Thus, we performed the first genome-wide association study (GWAS) for five growth traits in 172 brown-marbled groupers with 43,688 SNPs detected by ddRAD-seq. We identified a total of 5 significant and 18 suggestive QTLs located in multiple chromosomes associated with growth traits. In the 20 kb window of the significant SNPs and suggestive SNPs, 5 and 14 potential candidate genes affecting growth were detected, respectively. Five potential candidate genes near the significantly associated SNPs were selected for expression analysis. Among of which, bmp2k, wasf1, and acyp2 involved in bone development, maintenance of mitochondrion structure, and metabolism were differentially expressed. Interestingly, the SNP 23:29601315 located in the intron of bmp2k was significantly associated with body weight, body length, body height, and body thickness and suggestively associated with total length. We verified the locus using another new group including 123 individuals. The results showed that individuals with CC genotype have better growth traits comparing other individuals. Our findings not only contribute to understanding the molecular mechanism of growth regulation, but also promote the advance of marker-assisted selection in brown-marbled grouper.