Construction of a High-Density Genetic Map and Quantitative Trait Locus Mapping in the Manila clam Ruditapes philippinarum

Characterization of Eighty-Eight Single-Nucleotide Polymorphism Markers in the Manila Clam Ruditapes philippinarum Based on High-Resolution Melting (HRM) Analysis

Single-nucleotide polymorphisms (SNPs) are the most commonly used DNA markers in population genetic studies. We used the Illumina HiSeq4000 platform to develop single-nucleotide polymorphism (SNP) markers for Manila clam Ruditapes philippinarum using restriction site-associated DNA sequencing (RAD-seq) genotyping. Eighty-eight SNP markers were successfully developed by using high-resolution melting (HRM) analysis, with a success rate of 44%. SNP markers were analyzed for genetic diversity in two clam populations. The observed heterozygosity per locus ranged from 0 to 0.9515, while the expected heterozygosity per locus ranged from 0.0629 to 0.4997. The value of FIS was estimated to be from -0.9643 to 1.0000. The global Fst value was 0.1248 (p < 0.001). After Bonferroni correction, 15 loci deviated significantly from the Hardy-Weinberg equilibrium (p < 0.0006). These SNP markers provide a valuable resource for population and conservation genetics studies in this commercially important species.

Age-Dependent Growth-Related QTL Variations in Pacific Abalone, Haliotis discus hannai

Pacific abalone is a high-value, commercially important marine invertebrate. It shows low growth as well as individual and yearly growth variation in aquaculture. Marker-assisted selection breeding could potentially resolve the problem of low and variable growth and increase genetic gain. Expression of quantitative trait loci (QTLs) for growth-related traits, viz., body weight, shell length, and shell width were analyzed at the first, second, and third year of age using an F1 cross population. A total of 37 chromosome-wide QTLs were identified in linkage groups 01, 02, 03, 04, 06, 07, 08, 10, 11, 12, and 13 at different ages. None of the QTLs detected at any one age were expressed in all three age groups. This result suggests that growth-related traits at different ages are influenced by different QTLs in each year. However, multiple-trait QTLs (where one QTL affects all three traits) were detected each year that are also age-specific. Eleven multiple-trait QTLs were detected at different ages: two QTLs in the first year; two QTLs in the second year; and seven QTLs in the third year. As abalone hatcheries use three-year-old abalone for breeding, QTL-linked markers that were detected at the third year of age could potentially be used in marker-assisted selection breeding programs.

Examination of the potential roles of insulin-like peptide receptor in regulating the growth of Manila clam Ruditapes philippinarum

The role of insulin/insulin-like growth factor (IGF) signaling pathway in the growth regulation of marine invertebrates has not been fully studied. In this study, the economically important species Ruditapes philippinarum was sacrificed to clarify the role of IGF system in the growth regulation of R. philippinarum by real-time quantitative PCR. Systematic bioinformatics analysis can identify the major genes of IGF signaling pathway and insulin-like peptide receptor (ILPR) - mediated signaling pathway in R. philippinarum. The expression levels of IGF and its downstream signaling pathway genes in larger clams were significantly higher than those in small clams, indicating that they were involved in the growth regulation of R. philippinarum. These results suggest that IGF signaling pathway and ILPR mediated signaling pathway to regulate the growth of R. philippinarum.

The first high-density genetic map of common cockle (Cerastoderma edule) reveals a major QTL controlling shell color variation

Shell color shows broad variation within mollusc species and despite information on the genetic pathways involved in shell construction and color has recently increased, more studies are needed to understand its genetic architecture. The common cockle (Cerastoderma edule) is a valuable species from ecological and commercial perspectives which shows important variation in shell color across Northeast Atlantic. In this study, we constructed a high-density genetic map, as a tool for screening common cockle genome, which was applied to ascertain the genetic basis of color variation in the species. The consensus genetic map comprised 19 linkage groups (LGs) in accordance with the cockle karyotype (2n = 38) and spanned 1073 cM, including 730 markers per LG and an inter-marker distance of 0.13 cM. Five full-sib families showing segregation for several color-associated traits were used for a genome-wide association study and a major QTL on chromosome 13 associated to different color-traits was detected. Mining on this genomic region revealed several candidate genes related to shell construction and color. A genomic region previously reported associated with divergent selection in cockle distribution overlapped with this QTL suggesting its putative role on adaptation.

Animal board invited review: Widespread adoption of genetic technologies is key to sustainable expansion of global aquaculture

•

The extent of application of genetic technologies to aquaculture production varies widely by species and geography.

•

Achieving a more universal application of seed derived from scientifically based breeding programmes is an important goal in order to meet increasing global demands for seafood production.

•

This article reviews the status of genetic technologies across the world’s top 10 highly produced species.

•

Opportunities and barriers to achieving broad-scale uptake of genetic technologies in global aquaculture are discussed.

•

A future outlook for potential disruptive genetic technologies and how they might affect global aquaculture production is given.

Aquaculture production comprises a diverse range of species, geographies, and farming systems. The application of genetics and breeding technologies towards improved production is highly variable, ranging from the use of wild-sourced seed through to advanced family breeding programmes augmented by genomic techniques. This technical variation exists across some of the most highly produced species globally, with several of the top ten global species by volume generally lacking well-managed breeding programmes. Given the well-documented incremental and cumulative benefits of genetic improvement on production, this is a major missed opportunity. This short review focusses on (i) the status of application of selective breeding in the world’s most produced aquaculture species, (ii) the range of genetic technologies available and the opportunities they present, and (iii) a future outlook towards realising the potential contribution of genetic technologies to aquaculture sustainability and global food security.

Comparative transcriptomic analysis revealed dynamic changes of distinct classes of genes during development of the Manila clam (Ruditapes philippinarum)

Background

The Manila clam Ruditapesphilippinarum is one of the most economically important marine shellfish. However, the molecular mechanisms of early development in Manila clams are largely unknown. In this study, we collected samples from 13 stages of early development in Manila clam and compared the mRNA expression pattern between samples by RNA-seq techniques.

Results

We applied RNA-seq technology to 13 embryonic and larval stages of the Manila clam to identify critical genes and pathways involved in their development and biological characteristics. Important genes associated with different morphologies during the early fertilized egg, cell division, cell differentiation, hatching, and metamorphosis stages were identified. We detected the highest number of differentially expressed genes in the comparison of the pediveliger and single pipe juvenile stages, which is a time when biological characteristics greatly change during metamorphosis. Gene Ontology (GO) enrichment analysis showed that expression levels of microtubule protein-related molecules and Rho genes were upregulated and that GO terms such as ribosome, translation, and organelle were enriched in the early development stages of the Manila clam. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the foxo, wnt, and transforming growth factor-beta pathways were significantly enriched during early development. These results provide insights into the molecular mechanisms at work during different periods of early development of Manila clams.

Conclusion

These transcriptomic data provide clues to the molecular mechanisms underlying the development of Manila clam larvae. These results will help to improve Manila clam reproduction and development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08813-0.

Genotyping-by-Sequencing of the regional Pacific abalone (Haliotis discus) genomes reveals population structures and patterns of gene flow

Continuous monitoring of the present genetic status is essential to preserve the genetic resource of wild populations. In this study, we sequenced regional Pacific abalone Haliotis discus samples from three different locations around the Korean peninsula to assess population structure, utilizing Genotyping-by-Sequencing (GBS) method. Using PstI enzyme for genome reduction, we demonstrated the resultant library represented the whole genome region with even spacing, and as a result 16,603 single nucleotide variants (SNVs) were produced. Genetic diversity and population structure were investigated using several methods, and a strong genetic heterogeneity was observed in the Korean abalone populations. Additionally, by comparison of the variant sets among population groups, we were able to discover 26 Korean abalone population-specific SNVs, potentially associated with phenotype differences. This is the first study demonstrating the feasibility of GBS for population genetic study on H. discus. Our results will provide valuable data for the genetic conservation and management of wild abalone populations in Korea and help future GBS studies on the marine mollusks.

Transcriptomic analysis provides insights into candidate genes and molecular pathways involved in growth of Manila clam Ruditapes philippinarum

Growth is one of the most important traits of aquaculture breeding programs. Understanding the mechanisms underlying growth differences between individuals can contribute to improving growth rates through more efficient breeding schemes. Ruditapes philippinarum is an economically important marine bivalve. In order to gain insights into the molecular mechanisms to growth variability in marine shellfish, we conducted the transcriptome sequencing and examined the expression differences in growth-related gene and molecular pathways involved in growth trait of R. philippinarum. In this study, we investigated the molecular and gene expression differences in fast-growing and slow-growing Manila clam and focused on the analysis of the differential expression patterns of specific genes associated with growth by RNA-seq and qPCR analysis. A total of 61 differentially expressed genes (DEGs) were captured significantly differentially expressed, and were categorized into Ras signaling pathway, hedgehog signaling pathway, AMPK signaling pathway, p53 signaling pathway, regulation of actin cytoskeleton, focal adhesion, mTOR signaling pathway, VEGF signaling pathway, and TGF-beta signaling pathway. A total of 34 growth-related genes were validated significantly and up/downregulated at fast growing and slow growing of R. philippinarum. Functional enrichment analysis revealed the insulin signaling pathway, PI3K-Akt signaling pathway, and mTOR signaling pathway play pivotal roles in molecular function and regulation of growth trait in R. philippinarum. The growth-related genes and pathways obtained here provide important insights into the molecular basis of physiological acclimation, metabolic activity, and growth variability in marine bivalves.

Molecular Cloning and Polymorphism Analysis of PmFGF18 from Pinctada fucata martensii

Fibroblast growth factor 18 (FGF18) plays an important functional role in skeletal growth and development. The FGF18 gene was characterized in pearl oyster Pinctada fucata martensii (PmFGF18) with the full-length sequence containing an open reading frame of 714 bp encoding 237 amino acids. The domain analysis of PmFGF18 showed a distinctive FGF domain, with a high similarity to FGF18 protein sequences from Crassostrea gigas (43.35%) and C. virginica (37.43%). PmFGF18 expression was revealed in all analyzed tissues with a significantly higher expression level in the fast-growing group than the slow-growing group. The analysis of PmFGF18 polymorphism demonstrated 33 SNPs (single nucleotide polymorphisms) in the CDS and promoter region of PmFGF18 sequence. Association analysis revealed 19 SNPs (2 SNPs from CDS and 17 SNPs from the promoter region) associating significantly with growth traits. Among the associated SNPs, one SNP g.50918198 A > C was verified in the other breeding line. Therefore, PmFGF18 can be utilized as a candidate gene for growth, and its related SNPs could be used in selective breeding of P. f. martensii for the improvement of growth traits.

A proteomic study of resistance to Brown Ring disease in the Manila clam, Ruditapes philippinarum

Marine mollusk aquaculture has more than doubled over the past twenty years, accounting for over 15% of total aquaculture production in 2016. Infectious disease is one of the main limiting factors to the development of mollusk aquaculture, and the difficulties inherent to combating pathogens through antibiotic therapies or disinfection have led to extensive research on host defense mechanisms and host-pathogen relationships. It has become increasingly clear that characterizing the functional profiles of response to a disease is an essential step in understanding resistance mechanisms and moving towards more effective disease control. The Manila clam, Ruditapes philippinarum, is a main cultured bivalve species of economic importance which is affected by Brown Ring disease (BRD), an infection induced by the bacterium Vibrio tapetis. In this study, juvenile Manila clams were subjected to a 28-day controlled challenge with Vibrio tapetis, and visual and molecular diagnoses were carried out to distinguish two extreme phenotypes within the experimental clams: uninfected ("RES", resistant) and infected ("DIS", diseased) post-challenge. Total protein extractions were carried out for resistant and diseased clams, and proteins were identified using LC-MS/MS. Protein sequences were matched against a reference transcriptome of the Manila clam, and protein intensities based on label-free quantification were compared to reveal 49 significantly accumulated proteins in resistant and diseased clams. Proteins with known roles in pathogen recognition, lysosome trafficking, and various aspects of the energy metabolism were more abundant in diseased clams, whereas those with roles in redox homeostasis and protein recycling were more abundant in resistant clams. Overall, the comparison of the proteomic profiles of resistant and diseased clams after a month-long controlled challenge to induce the onset of Brown Ring disease suggests that redox homeostasis and maintenance of protein structure by chaperone proteins may play important and interrelated roles in resistance to infection by Vibrio tapetis in the Manila clam.

High-Density Genetic Map Construction and QTL Mapping of Leaf and Needling Traits in Ziziphus jujuba Mill

The Chinese jujube (Ziziphus jujuba Mill., 2n = 2x = 24), one of the most popular fruit trees in Asia, is widely cultivated and utilized in China, where it is traditionally consumed as both a fresh and dried food resource. A high-density genetic map can provide the necessary framework for quantitative trait loci (QTL) analyses and map-based gene cloning and molecular breeding. In this study, we constructed a new high-density genetic linkage map via a genotyping-by-sequencing approach. For the consensus linkage map, a total of 3,792 markers spanning 2,167.5 cM were mapped onto 12 linkage groups, with an average marker interval distance of 0.358 cM. The genetic map anchored 301 Mb (85.7%) of scaffolds from the sequenced Z. jujuba "Junzao" genome. Based on this genetic map, 30 potential QTLs were detected, including 27 QTLs for leaf traits and 3 QTLs for needling length. This high-density genetic map and the identified QTLs for relevant agronomic traits lay the groundwork for functional genetic mapping, map-based cloning, and marker-assisted selection in jujube.

Construction of a high-resolution genetic map and identification of quantitative trait loci for salt tolerance in jute (Corchous spp.)

BACKGROUND

Jute (Corchorus spp.) is the most important natural fiber crop after cotton in terms of cultivation area and production. Salt stress greatly restricts plant development and growth. A high-density genetic linkage map is the basis of quantitative trait locus (QTLs) mapping. Several high-density genetic maps and QTLs mapping related to salt tolerance have been developed through next-generation sequencing in many crop species. However, such studies are rare for jute. Only several low-density genetic maps have been constructed and no salt tolerance-related QTL has been mapped in jute to date.

RESULTS

We developed a high-density genetic map with 4839 single nucleotide polymorphism markers spanning 1375.41 cM and an average distance of 0.28 cM between adjacent markers on seven linkage groups (LGs) using an F2 jute population, LGs ranged from LG2 with 299 markers spanning 113.66 cM to LG7 with 1542 markers spanning 350.18 cM. In addition, 99.57% of gaps between adjacent markers were less than 5 cM. Three obvious and 13 minor QTLs involved in salt tolerance were identified on four LGs explaining 0.58-19.61% of the phenotypic variance. The interval length of QTL mapping varied from 1.3 to 20.2 cM. The major QTL, qJST-1, was detected under two salt stress conditions that explained 11.81 and 19.61% of the phenotypic variation, respectively, and peaked at 19.3 cM on LG4.

CONCLUSIONS

We developed the first high-density and the most complete genetic map of jute to date using a genotyping-by-sequencing approach. The first QTL mapping related to salt tolerance was also carried out in jute. These results should provide useful resources for marker-assisted selection and transgenic breeding for salt tolerance at the germination stage in jute.

Quantitative Genetic Mapping and Genome Assembly in the Lesser Wax Moth Achroia grisella

Specific characteristics of the male Achroia grisella acoustic mating signal determine a male’s attractiveness toward females. These features are genetically variable in populations, and mapping experiments have been used to identify loci contributing to song variation, and understand the evolutionary forces acting on this important sexual trait. Here we built on this foundation and carried out QTL (Quantitative Trait Locus) mapping using >1,000 recombinant individuals, genotyping this large cohort at thousands of sequence-based markers covering the entire collection of 30 A. grisella chromosomes. This dense marker set, coupled with our development of an annotated, draft genome of A. grisella, allowed us to link >3,000 genome scaffolds, >10,000 predicted genes, and close to 275Mb of genome sequence to chromosomes. Our QTL mapping confirmed a fraction of the QTL identified in a previous study, and additionally revealed novel loci. Collectively, QTL explained only small fractions of the phenotypic variance, suggesting many more causative factors remain below the detection threshold of our study. A surprising, and ultimately challenging feature of our study was the low level of intrachromosomal recombination present in our mapping population. This led to difficulty ordering markers along linkage groups, necessitating a chromosome-by-chromosome mapping approach, rather than true interval mapping, and precluded confident ordering/orienting of scaffolds along each chromosome. Nonetheless, our study increased the genomic resources available for the A. grisella system. Enabled by ever more powerful technologies, future investigators will be able to leverage our data to provide more detailed genetic dissection of male song variation in A. grisella.

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 genetic map: genotyping by sequencing (GBS) to map purple seed coat color (Psc) in hulless barley

Background

Colored hulless barley are more suitable in food processing compared to normal (yellow) varieties because it is rich in bioactive compounds and produces higher extraction pearling fractions. Therefore, seed coat color is an important agronomic trait for the breeding and study of hulless barley.

Results

Genotyping-by-sequencing single-nucleotide polymorphism (GBS-SNP) analysis of a doubled haploid (DH) mapping population (Nierumuzha × Kunlun10) was conducted to map the purple seed coat color genes (Psc). A high-density genetic map of hulless barley was constructed, which contains 3662 efficient SNP markers with 1129 bin markers. Seven linkage groups were resolved, which had a total length of 645.56 cM. Chromosome length ranged from 60.21 cM to 127.21 cM, with average marker density of 0.57 cM. A total of five loci accounting for 3.79% to 23.86% of the observed phenotypic variation for Psc were detected using this high-density map. Five structural candidate genes (F3’M, HID, UF3GT, UFGT and 5MAT) and one regulatory factor (Ant1) related to flavonoid or anthocyanin biosynthesis were identified..

Conclusions

Five structural candidate genes and one regulatory factor related to flavonoid or anthocyanin biosynthesis have been identified using a high-density genetic map of hulless barley. This study lays the foundation for map-based cloning of Psc but provides a valuable tool for studying marker-trait associations and its application to marker-assisted breeding of hulless barley.

Electronic supplementary material

The online version of this article (10.1186/s41065-018-0072-6) contains supplementary material, which is available to authorized users.

Highly dense linkage maps from 31 full-sibling families of turbot (Scophthalmus maximus) provide insights into recombination patterns and chromosome rearrangements throughout a newly refined genome assembly

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

Genomic Tools and Selective Breeding in Molluscs

The production of most farmed molluscs, including mussels, oysters, scallops, abalone, and clams, is heavily dependent on natural seed from the plankton. Closing the lifecycle of species in hatcheries can secure independence from wild stocks and enables long-term genetic improvement of broodstock through selective breeding. Genomic techniques have the potential to revolutionize hatchery-based selective breeding by improving our understanding of the characteristics of mollusc genetics that can pose a challenge for intensive aquaculture and by providing a new suite of tools for genetic improvement. Here we review characteristics of the life history and genetics of molluscs including high fecundity, self-fertilization, high genetic diversity, genetic load, high incidence of deleterious mutations and segregation distortion, and critically assess their impact on the design and effectiveness of selective breeding strategies. A survey of the results of current breeding programs in the literature show that selective breeding with inbreeding control is likely the best strategy for genetic improvement of most molluscs, and on average growth rate can be improved by 10% per generation and disease resistance by 15% per generation across the major farmed species by implementing individual or family-based selection. Rapid advances in sequencing technology have resulted in a wealth of genomic resources for key species with the potential to greatly improve hatchery-based selective breeding of molluscs. In this review, we catalog the range of genomic resources currently available for molluscs of aquaculture interest and discuss the bottlenecks, including lack of high-quality reference genomes and the relatively high cost of genotyping, as well as opportunities for applying genomics-based selection.

A Microarray Study of Carpet-Shell Clam (Ruditapes decussatus) Shows Common and Organ-Specific Growth-Related Gene Expression Differences in Gills and Digestive Gland

Growth rate is one of the most important traits from the point of view of individual fitness and commercial production in mollusks, but its molecular and physiological basis is poorly known. We have studied differential gene expression related to differences in growth rate in adult individuals of the commercial marine clam Ruditapes decussatus. Gene expression in the gills and the digestive gland was analyzed in 5 fast-growing and five slow-growing animals by means of an oligonucleotide microarray containing 14,003 probes. A total of 356 differentially expressed genes (DEG) were found. We tested the hypothesis that differential expression might be concentrated at the growth control gene core (GCGC), i.e., the set of genes that underlie the molecular mechanisms of genetic control of tissue and organ growth and body size, as demonstrated in model organisms. The GCGC includes the genes coding for enzymes of the insulin/insulin-like growth factor signaling pathway (IIS), enzymes of four additional signaling pathways (Raf/Ras/Mapk, Jnk, TOR, and Hippo), and transcription factors acting at the end of those pathways. Only two out of 97 GCGC genes present in the microarray showed differential expression, indicating a very little contribution of GCGC genes to growth-related differential gene expression. Forty eight DEGs were shared by both organs, with gene ontology (GO) annotations corresponding to transcription regulation, RNA splicing, sugar metabolism, protein catabolism, immunity, defense against pathogens, and fatty acid biosynthesis. GO term enrichment tests indicated that genes related to growth regulation, development and morphogenesis, extracellular matrix proteins, and proteolysis were overrepresented in the gills. In the digestive gland overrepresented GO terms referred to gene expression control through chromatin rearrangement, RAS-related small GTPases, glucolysis, and energy metabolism. These analyses suggest a relevant role of, among others, some genes related to the IIS, such as the ParaHox gene Xlox, CCAR and the CCN family of secreted proteins, in the regulation of growth in bivalves.