Three EST-SSR markers associated with QTL for the growth of the clam Meretrix meretrix revealed by selective genotyping

Mapping Genetic Loci for Quantitative Traits of Golden Shell Color, Mineral Element Contents, and Growth-Related Traits in Pacific Oyster (Crassostrea gigas)

Golden shell color and mineral content are important economic traits of Pacific oyster (Crassostrea gigas). In this study, we mapped a series of quantitative trait loci (QTLs) that control zinc (Zn) and magnesium (Mg) content, shell color and growth performance to two sex-averaged linkage maps from the FAM-A and FAM-B families. In total, ten QTLs were identified in seven linkage groups (LGs) in the FAM-B family, and seven QTLs were identified in four linkage groups in the FAM-A family. Two QTLs affecting the trait of golden shell color were identified in LG8 of the FAM-A and LG10 of the FAM-B families, which could explain 20.2 and 10.5% of the phenotypic variations, respectively. Two QTLs for Zn content were identified that could contribute to 17.9 and 34.44% of the phenotypic variations in FAM-A. Six QTLs for Zn and Mg contents were identified in four LGs (LG1, LG2, LG5, and LG9) in FAM-B, which explained 13.5-26.7% of the phenotypic variations. In addition, seven QTLs related to oyster growth were recognized in both FAM-A and FAM-B families accounting for 14.6-36.7% of the phenotypic variations. All of the DNA markers in QTL regions were blasted and 14 genes associated with above traits were identified. The mRNA expression of these genes was determined by quantitative RT-PCR. These QTLs and candidate genes could be used as potential targets for marker-assisted selection in C. gigas breeding.

Identification of SNP markers associated with tolerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in Litopenaeus vannamei

The high concentration of ammonia from deteriorated aquaculture environments and the intensive culture system could increase the susceptibility to pathogens and even cause high mortality in Litopenaeus vannamei. In addition, we have revealed that the ammonia-tolerant shrimp also have high disease resistance in L. vannamei. In the present study, in order to identify SNP markers associated with tolerance to ammonia toxicity, we developed and characterized SNPs from our previous transcriptome sequencing data of ammonia-stressed and control groups, and a marker-trait association analysis was performed for marker-assisted selection (MAS) to increase production in L. vannamei. A total of 318,919 SNPs were identified from the transcriptome sequences, and 25,772 SNPs were found from the 1826 ammonia-responsive genes with functional annotation. We selected 49 SNPs from 26 ammonia-responsive genes that had strong homologies to known genes in the shrimp and probably involved in immune function as candidate markers for genotyping, among which 39 SNPs were polymorphic for further marker-trait association analysis with the ammonia-tolerant (AT) and ammonia-sensitive (AS) groups. Finally, 12 out of the 49 SNP markers were identified to be associated with ammonia tolerance, containing 10 loci with significantly different allele frequencies and 10 loci with significantly different genotyping frequencies between the AT and AS groups. Among the associated markers, the G allele of TSP-1 (the first locus from the thrombospondin gene), the A allele of TSP-3, and the C allele of XBP1-5 (the fifth locus from X-box binding protein 1) only presented in the AT groups, but they were absent from the AS groups, which would be the preference of the MAS for the ammonia-tolerant shrimp. In addition, when the 12 associated SNP markers were used for analysis, the genetic diversity of the AT groups were significantly higher than that of the AS groups, but when the 39 loci were used there was no difference. This is the first report for the markers associated with ammonia tolerance in this species, indirectly with disease resistance, which provided important potential for genetic selection to increase survival rate and production in shrimp farming.

Isolation of novel microsatellite markers and their application for genetic diversity and parentage analyses in sika deer

Every part of the sika deer (Cervus nippon) body is valuable traditional Chinese medicine. And sika deer is the most important semi-domestic medicinal animal that is widely bred in Jilin province northeast of China. But few studies had been conducted to characterize the microsatellite markers derived from sika deer. We firstly used IlluminaHiSeq™2500 sequencing technology obtained 125Mbp genomic data of sika deer. Using microsatellite identification tool (MISA), 22,479 microsatellites were identified. From these data, 100 potential primers were selected for further polymorphic validation, finally, 76 primer pairs were successfully amplified and 29 primer pairs were found to be obvious polymorphic in 8 different individuals. Using those polymorphic microsatellite markers, we analyzed the genetic diversity of Jilin sika deer population. The mean number of alleles of the 29 loci is 9.31 based on genotyping blood DNA from 96 Jilin sika deer; The mean expected heterozygosity and polymorphic information content (PIC) value of the 29 loci is 0.72 and 0.68 respectively, and among which 26 loci are highly polymorphic (PIC>0.50). According to the electrophoretic results and PIC value of these 29 loci, 10 loci with combined paternity exclusion probabilities>99.99% were selected to use in parentage verification for 16 sika deer. All the offspring of a family could be successfully assigned to their biological father. These microsatellite markers generated in this study could greatly facilitate future studies of molecular breeding in sika deer.

A Genome-Wide Association Study Reveals That Genes with Functions for Bone Development Are Associated with Body Conformation in Catfish

Body conformation is of great scientific and commercial interest for aquaculture fish species because it affects biological adaptation of the organism to environments, and is of economic importance to the aquaculture industry considering its direct effect on fillet yield. Catfish is the primary aquaculture species in the USA. Two major species used in the aquaculture industry, channel catfish and blue catfish, differ in body shape and therefore the backcross progenies serve as a good model for quantitative trait locus (QTL) analysis. Here, a genome-wide association study (GWAS) with hybrid catfish was conducted to identify the QTL for body conformation, including deheaded body length (DBL), body length (BL), body depth (BD), and body breadth (BB), which were all standardized by cubic root of body weight. Overall, the results indicate that the traits are polygenic. For DBL, linkage group (LG) 2 and LG 24 contain significant QTL, and LG 13 and LG 26 contain suggestively associated QTL (-log₁₀(P value) > 4.5). Compared with DBL, additional SNPs were identified to be associated with body length on LG 2, LG 7, and LG 18. Although no significant QTL for body depth was found, three suggestively associated QTLs were identified on LG 5, LG 13, and LG 14. No SNP for body breadth reached the threshold for suggestive association. Genes close to the associated SNPs were determined, many of which are known to be involved in bone development. This work therefore provides the basis for future identification of causal genes for the control of body conformation.

Mapping QTL for Resistance Against Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), leads to mass mortality in mariculture. However, phenotypic selection for resistance against VNN is very difficult. To facilitate marker-assisted selection (MAS) for resistance against VNN and understanding of the genetic architecture underlying the resistance against this disease, we mapped quantitative trait loci (QTL) for resistance against VNN in Asian seabass. We challenged fingerlings at 37 days post-hatching (dph), from a single back-cross family, with NNV at a concentration of 9 × 10(6) TCID50/ml for 2 h. Daily mortalities were recorded and collected. A panel of 330 mortalities and 190 surviving fingerlings was genotyped using 149 microsatellites with 145 successfully mapped markers covering 24 linkage groups (LGs). Analysis of QTL for both resistance against VNN and survival time was conducted using interval mapping. Five significant QTL located in four LGs and eight suggestive QTL in seven LGs were identified for resistance. Another five significant QTL in three LGs and five suggestive QTL in three LGs were detected for survival time. One significant QTL, spanning 3 cM in LG20, was identified for both resistance and survival time. These QTL explained 2.2-4.1% of the phenotypic variance for resistance and 2.2-3.3% of the phenotypic variance for survival time, respectively. Our results suggest that VNN resistance in Asian seabass is controlled by many loci with small effects. Our data provide information for fine mapping of QTL and identification of candidate genes for a better understanding of the mechanism of disease resistance.

Characterization of a long-chain fatty acid-CoA ligase 1 gene and association between its SNPs and growth traits in the clam Meretrix meretrix

Long-chain fatty acid-CoA ligases (ACSLs) play crucial roles in fatty acid (FA) metabolism. They convert free long-chain FA into acyl-CoAs, which are key intermediates in both anabolic and catabolic pathways. A long-chain fatty acid-CoA ligase gene was cloned in the clam Meretrix meretrix (MmeACSL1), with a full-length cDNA of 1865 bp encoding 475 amino acids. Its expression was only detected in hepatopancreas by semi-quantitative reverse transcription PCR. Expression level of MmeACSL1 exhibited a significant increase in a starvation experiment (P<0.05). This indicates that MmeACSL1 plays an important role in normal metabolism of M. meretrix and may be involved in energy supply and storage. Two exon SNPs and six intron SNPs were developed in this gene by direct sequencing. A marker-trait association analysis showed that five of these SNPs were significantly associated with growth traits (P<0.05). Also, haplotypes comprised of the five SNPs were revealed to be significantly growth-related (P<0.05), which further corroborates the relationship of this gene to growth traits of M. meretrix.

Identification of a Serum amyloid A gene and the association of SNPs with Vibrio-resistance and growth traits in the clam Meretrix meretrix

Serum amyloid A (SAA), an acute response protein as well as an apolipoprotein, is considered to play crucial roles in both innate immunity and lipid metabolism. In this study, a SAA gene (MmSAA) was identified in the clam Meretrix meretrix. The full length DNA of MmSAA was 1407bp, consisting of three exons and two introns. The distribution of MmSAA in clam tissues was examined with the highest expression in hepatopancreas. In response to the Vibrio parahaemolyticus challenge, MmSAA mRNA showed significantly higher expression at 24 h post-challenge in experimental clams (P < 0.05). Forty-eight single nucleotide polymorphisms (SNPs) in the DNA partial sequence of MmSAA were discovered and examined for their association with Vibrio-resistance and growth traits, respectively. The single SNP association analysis indicated that five single SNPs (g.42, g.72, g.82, g.147 and g.165) were significantly associated with Vibrio-resistance (P < 0.05). Haplotype analysis produced additional support for association with the Chi-square values 6.393 (P = 0.012). Among the five selected SNPs, the effect of a missense mutation (g.82, A → G) was detected by site-directed mutagenesis with fusion expression of protein assay, and the result showed that the recombinant plasmids containing wild-type pET30a-MmSAA had more inhibition effect than the mutant ones on the growth rate of the host bacteria. In addition, four growth traits of the clams in 09G3SPSB population were recorded and the SNP g.176 was found to be significantly associated with the growth traits with the Global score value 0.790 (P = 0.015). Our findings suggested that common genetic variation in MmSAA might contribute to the risk of susceptibility to Vibrio infection and might be associated with the growth traits in the clams M. meretrix, and more works are still needed to validate these SNPs as potential markers for actual selective breeding.

Development of Vibrio spp. infection resistance related SNP markers using multiplex SNaPshot genotyping method in the clam Meretrix meretrix

The clam Meretrix meretrix is a commercially important mollusc species in the coastal areas of South and Southeast Asia. In the present study, large-scale SNPs were genotyped by the Multiplex SNaPshot genotyping method among the stocks of M. meretrix with different Vibrio spp. infection resistance profile. Firstly, the AUTOSNP software was applied to mine SNPs from M. meretrix transcriptome, and 323 SNP loci (including 120 indels) located on 64 contigs were selected based on Uniprot-GO associations. Then, 38 polymorphic SNP loci located on 15 contigs were genotyped successfully in the clam stocks with different resistance to Vibrio parahaemolyticus infection (11-R and 11-S groups). Pearson's Chi-square test was applied to compare the allele and genotype frequency distributions of the SNPs between the different stocks, and seven SNP markers located on three contigs were found to be associated with V. parahaemolyticus infection resistance trait. Haplotype-association analysis showed that six haplotypes had significantly different frequency distributions in 11-S and 11-R (P < 0.05). With selective genotyping between 09-R and 09-C populations, which had different resistance to Vibrio harveyi infection, four out of the seven selected SNPs had significantly different distributions (P < 0.05) and therefore they were considered to be associated with Vibrio spp. infection resistance. Sequence alignments and annotations indicated that the contigs containing the associated SNPs had high similarity to the immune related genes. All these results would be useful for the future marker-assisted selection of M. meretrix strains with high Vibrio spp. infection resistance.

Genetic mapping and QTL analysis of growth-related traits in Pinctada fucata using restriction-site associated DNA sequencing

The pearl oyster, Pinctada fucata (P. fucata), is one of the marine bivalves that is predominantly cultured for pearl production. To obtain more genetic information for breeding purposes, we constructed a high-density linkage map of P. fucata and identified quantitative trait loci (QTL) for growth-related traits. One F1 family, which included the two parents, 48 largest progeny and 50 smallest progeny, was sampled to construct a linkage map using restriction site-associated DNA sequencing (RAD-Seq). With low coverage data, 1956.53 million clean reads and 86,342 candidate RAD loci were generated. A total of 1373 segregating SNPs were used to construct a sex-average linkage map. This spanned 1091.81 centimorgans (cM), with 14 linkage groups and an average marker interval of 1.41 cM. The genetic linkage map coverage, Coa, was 97.24%. Thirty-nine QTL-peak loci, for seven growth-related traits, were identified using the single-marker analysis, nonparametric mapping Kruskal-Wallis (KW) test. Parameters included three for shell height, six for shell length, five for shell width, four for hinge length, 11 for total weight, eight for soft tissue weight and two for shell weight. The QTL peak loci for shell height, shell length and shell weight were all located in linkage group 6. The genotype frequencies of most QTL peak loci showed significant differences between the large subpopulation and the small subpopulation (P<0.05). These results highlight the effectiveness of RAD-Seq as a tool for generation of QTL-targeted and genome-wide marker data in the non-model animal, P. fucata, and its possible utility in marker-assisted selection (MAS).

A new simple method for improving QTL mapping under selective genotyping

The selective genotyping approach, where only individuals from the high and low extremes of the trait distribution are selected for genotyping and the remaining individuals are not genotyped, has been known as a cost-saving strategy to reduce genotyping work and can still maintain nearly equivalent efficiency to complete genotyping in QTL mapping. We propose a novel and simple statistical method based on the normal mixture model for selective genotyping when both genotyped and ungenotyped individuals are fitted in the model for QTL analysis. Compared to the existing methods, the main feature of our model is that we first provide a simple way for obtaining the distribution of QTL genotypes for the ungenotyped individuals and then use it, rather than the population distribution of QTL genotypes as in the existing methods, to fit the ungenotyped individuals in model construction. Another feature is that the proposed method is developed on the basis of a multiple-QTL model and has a simple estimation procedure similar to that for complete genotyping. As a result, the proposed method has the ability to provide better QTL resolution, analyze QTL epistasis, and tackle multiple QTL problem under selective genotyping. In addition, a truncated normal mixture model based on a multiple-QTL model is developed when only the genotyped individuals are considered in the analysis, so that the two different types of models can be compared and investigated in selective genotyping. The issue in determining threshold values for selective genotyping in QTL mapping is also discussed. Simulation studies are performed to evaluate the proposed methods, compare the different models, and study the QTL mapping properties in selective genotyping. The results show that the proposed method can provide greater QTL detection power and facilitate QTL mapping for selective genotyping. Also, selective genotyping using larger genotyping proportions may provide roughly equivalent power to complete genotyping and that using smaller genotyping proportions has difficulties doing so. The R code of our proposed method is available on http://www.stat.sinica.edu.tw/chkao/.

Three vibrio-resistance related EST-SSR markers revealed by selective genotyping in the clam Meretrix meretrix

The clam Meretrix meretrix is an important commercial bivalve distributed in the coastal areas of South and Southeast Asia. In this study, marker-trait association analyses were performed based on the stock materials of M. meretrix with different vibrio-resistance profile obtained by selective breeding. Forty-eight EST-SSR markers were screened and 27 polymorphic SSRs of them were genotyped in the clam stocks with different resistance to Vibrio parahaemolyticus (11-R and 11-S) and to Vibrio harveyi (09-R and 09-C). Allele frequency distributions of the SSRs among different stocks were compared using Pearson's Chi-square test, and three functional EST-SSR markers (MM959, MM4765 and MM8364) were found to be associated with vibrio-resistance trait. The 140-bp allele of MM959 and 128-bp allele of MM4765 had significantly higher frequencies in resistant groups (11-R and 09-R) than in susceptive/control groups (11-S and 09-C) (P < 0.01), which suggested that the clams carrying these two alleles have stronger resistance against vibrio. Clam individuals of 11-S were divided into three subgroups based on the survival time post-challenge, and the multi-dimensional scaling (MDS) analysis showed that clusters generated by genetic similarity revealed by the three SSR markers were consistent with the three subgroups distinctions. The putative functions of contig959, contig4765 and contig8364 also suggested that the three SSR-involved genes might play important roles in immunity of M. meretrix. All these results supported that EST-SSR markers MM959, MM4765 and MM8364 were associated with vibrio-resistance and would be useful for marker-assisted selection (MAS) in M. meretrix genetic breeding.

Development of molecular resources for an intertidal clam, Sinonovacula constricta, using 454 transcriptome sequencing

Background

The razor clam Sinonovacula constricta is a benthic intertidal bivalve species with important commercial value. Despite its economic importance, knowledge of its transcriptome is scarce. Next generation sequencing technologies offer rapid and efficient tools for generating large numbers of sequences, which can be used to characterize the transcriptome, to develop effective molecular markers and to identify genes associated with growth, a key breeding trait.

Results

Total RNA was isolated from the mantle, gill, liver, siphon, gonad and muscular foot tissues. High-throughput deep sequencing of S. constricta using 454 pyrosequencing technology yielded 859,313 high-quality reads with an average read length of 489 bp. Clustering and assembly of these reads produced 16,323 contigs and 131,346 singletons with average lengths of 1,376 bp and 458 bp, respectively. Based on transcriptome sequencing, 14,615 sequences had significant matches with known genes encoding 147,669 predicted proteins. Subsequently, previously unknown growth-related genes were identified. A total of 13,563 microsatellites (SSRs) and 13,634 high-confidence single nucleotide polymorphism loci (SNPs) were discovered, of which almost half were validated.

Conclusion

De novo sequencing of the razor clam S. constricta transcriptome on the 454 GS FLX platform generated a large number of ESTs. Candidate growth factors and a large number of SSRs and SNPs were identified. These results will impact genetic studies of S. constricta.